The origin of new genes: glimpses from the young and old

Predicting Protein Function in the AI and Big Data Era

It is an exciting time for researchers working to link proteins to their functions. Most techniques for extracting functional information from genomic sequences were developed several years ago, with major progress driven by the availability of big data. Now, groundbreaking advances in deep-learning and AI-based methods have enriched protein databases with three-dimensional information and offer the potential to predict biochemical properties and biomolecular interactions, providing key functional insights. This progress is expected to increase the proportion of functionally bright proteins in databases and deepen our understanding of life at the molecular level.

New perspectives on Drosophila melanogaster de novo gene origination revealed by investigation of ancient African genetic variation

De novo genes can be defined as sequences producing evolutionarily derived transcripts that are not homologous to transcripts produced in an ancestor. While they appear to be taxonomically widespread, there is little agreement regarding their abundance, their persistence times in genomes, the population genetic processes responsible for their spread or loss, or their possible functions. In Drosophila melanogaster, 2 approaches have been used to discover these genes and investigate their properties. One uses traditional comparative approaches and existing genomic resources and annotations. A second approach uses raw transcriptome data to discover unannotated genes for which there is no evidence of presence in related species. Investigations using the second approach have focused on D. melanogaster genotypes from recently established cosmopolitan populations. However, most of the genetic variation in the species is found in African populations, suggesting the possibility that fuller understanding of genetic novelties in the species may follow from studies of these populations. Here, we investigate de novo gene candidates expressed in testis and accessory glands in a sample of flies from Zambia and compare them with candidate de novo genes expressed in North American populations. We report a large number of previously undiscovered de novo gene candidates, most of which are expressed polymorphically. Many are predicted to code for secreted proteins. In spite of much different levels of genomic variation in Zambian and North American populations, they express similar numbers of candidate de novo genes. We find evidence from genetic analysis of Raleigh inbred lines that a fraction of rarely expressed gene candidates in this population represent deleterious transcription promoted by inbreeding depression. Many de novo gene candidates are expressed in multiple tissues and both sexes, raising questions about how they may interact with natural selection. The relative importance of positive and negative selection, however, remains unclear.

Genomic Analysis Reveals the Role of New Genes in Venom Regulatory Network of Parasitoid Wasps

New genes play a critical role in phenotypic diversity and evolutionary innovation. Parasitoid wasps, a highly abundant and diverse group of insects, parasitize other arthropods and exhibit remarkable evolutionary adaptations, such as evading host immune responses and exploiting host resources. However, the specific contributions of new genes to their unique traits remain poorly understood. Here, we identified 480 new genes that emerged after the Nasonia-Pteromalus divergence. Among these, 272 (56.7%) originated through DNA-mediated duplication, representing the largest proportion, followed by 77 (16.0%) derived from RNA-mediated duplication and 131 (27.3%) that arose de novo. Comparative analysis revealed that these new genes generally have shorter coding sequences and fewer exons compared to single-copy older genes conserved in the seven parasitoid wasps. These new genes are predominantly expressed in the reproductive glands and exhibit venom gland-biased expression. Notably, gene co-expression network analysis further identified that a new gene may act as a hub by interacting with older genes to regulate venom-related networks rather than directly encoding venom proteins. Together, our findings provide novel insights into the role of new genes in driving venom innovation in parasitoid wasps.

Genome-wide identification and expression analysis of PP2C gene families in two Chimonanthus species

The protein phosphatase 2 C (PP2C) plays a crucial role in the growth and development of plants. However, limited information on the PP2C genes in Chimonanthus spp. is available. Therefore, the comprehensive genome-wide identification and analysis of the PP2Cs gene family is necessary in Chimonanthus spp. to provide basic information for further study. In this work, 93 members of the CsPP2C gene family and 85 members of the CpPP2C gene family were identified. Based on the phylogenetic analysis, 93 CsPP2Cs and 85 CpPP2Cs genes were classified into 13 subgroups. Based on RNA-Seq data, specific expression patterns of CsPP2Cs in different tissues were identified. CsPP2C28, CsPP2C55, and CsPP2C17 showed high expression during leaf senescence, and combined with cis-element analysis, it is speculated that they may participate in regulating plant senescence. CsPP2C47, CsPP2C27, CsPP2C42, and CsPP2C41 may play an important role in responding to abscisic acid during seed dormancy and germination. These results contribute to a deeper understanding of the functions of the CsPP2C gene family and CpPP2C gene family, and providing candidate genes for genetic engineering and breeding to enhance important traits such as stress resistance and growth development in C. salicifolius and C. praecox.

DeNoFo: a file format and toolkit for standardised, comparable de novo gene annotation

Motivation

De novo genes emerge from previously non-coding regions of the genome, challenging the traditional view that new genes primarily arise through duplication and adaptation of existing ones. Characterised by their rapid evolution and their novel structural properties or functional roles, de novo genes represent a young area of research. Therefore, the field currently lacks established standards and methodologies, leading to inconsistent terminology and challenges in comparing and reproducing results.

Results

This work presents a standardised annotation format to document the methodology of de novo gene datasets in a reproducible way. We developed DeNoFo, a toolkit to provide easy access to this format that simplifies annotation of datasets and facilitates comparison across studies. Unifying the different protocols and methods in one standardised format, while providing integration into established file formats, such as fasta or gff, ensures comparability of studies and advances new insights in this rapidly evolving field.

Availability and Implementation

DeNoFo is available through the official Python Package Index (PyPI) and at https://github.com/EDohmen/denofo . All tools have a graphical user interface and a command line interface. The toolkit is implemented in Python3, available for all major platforms and installable with pip and uv.

Unveiling the functional nature of retrogenes in dinoflagellates

Retroposition is a gene duplication mechanism that uses RNA molecules as intermediaries to generate new gene copies. Dinoflagellates are proposed as an ideal model for exploring this process due to the tagging of retrogenes with DNA-encoded remnants of the dinoflagellate-specific splice-leader motif at their 5' end. We conducted a comprehensive search for retrogenes in dinoflagellate transcriptomes to uncover their functional nature and the processes underlying their redundancy. We obtained a high-confidence set of hypothetical functional retrogenes widespread through the dinoflagellate lineage. Through annotations and gene ontology enrichment analysis, we found that the functional diversity of retrogenes reflects the most prevalent and active processes during stress periods, particularly those involving post-translational modifications and cell signalling pathways. Additionally, the significant presence of retrogenes linked to specific biological processes involved in symbiosis and toxin production underscores the role of retrogenes in adaptation. The expression profile and codon composition similar to protein-coding genes confirm the operational status of retrogenes and strengthen the idea that retrogenes recapitulate parental gene expression and function. This study provides new evidence supporting widespread gene retroposition across dinoflagellates and highlights the functional link of retrogenes with the core activity of the cell.

Evolutionary landscape of plant chalcone isomerase-fold gene families

Flavonoids are crucial for plant survival and adaptive evolution, and chalcone isomerase (CHI) genes serve as key rate-limiting gene in the flavonoid biosynthesis pathway. It is important for plant adaptive evolution to comprehensively study the evolution and diversity of the CHI gene families. However, the CHI gene families in many plant lineages remain elusive. This study systematically identified CHI genes from 259 species including algae, bryophytes, ferns, gymnosperms, and angiosperms. A total of 1,738 CHI gene family members were discovered. We analyzed the diversity, distribution trajectory, and the driving forces of gene duplication during the evolution of the plant lineages. The present study is the first to identify potential type II and type IV CHI genes in the extant liverwort model species Marchantia polymorpha. The distribution pattern of CHI genes across the plant kingdom reveals that the origin of type II CHI can be traced back to the last common ancestor of bryophytes and vascular plants, and type III CHI may represent the ancestral form of the CHI gene family. The identification of conserved motifs showed significant differences in motif distribution among different CHI gene types. It was found that the drivers of gene duplication varied across plant lineages: dispersed duplications (DSD) were predominant in algae and bryophytes, whole-genome duplication (WGD) was the main driver in basal angiosperms and monocots, while tandem duplications (TD) predominating in eudicots. Structural clustering analysis demonstrated the 3-layer sandwich structure in the CHI-fold proteins remained conserved in the central region, while repeated loss of N-terminal sequences contributed to structural diversity. This study provides a deeper understanding of the evolution and diversity of the CHI-fold proteins and lays a theoretical foundation for further studies of their function and the identification of new functional CHI genes.

Genome assembly of Stewartia sinensis reveals origin and evolution of orphan genes in Theaceae

Orphan genes play crucial roles in diverse biological processes, but the evolutionary trajectories and functional divergence remain largely unexplored. The Theaceae family, including the economically and culturally important tea plant, offers a distinctive model to examine these aspects. Here, we integrated Nanopore long-read sequencing, Illumina short-read sequencing, and Hi-C methods to decode a pseudo-chromosomal genome assembly of Stewartia sinensis, from the earliest-diverging tribe of Theaceae, spanning 2.95 Gb. Comparative genomic analysis revealed the absence of recent whole-genome duplication events in the Theaceae ancestor, highlighting tandem duplications as the predominant mechanism of gene expansion. We identified 31,331 orphan genes, some of which appear to have ancient origins, suggesting early emergence with frequent gains and losses, while others seem more specific and recent. Notably, orphan genes are distinguished by shorter lengths, fewer exons and functional domains compared to genes that originate much earlier, like transcription factors. Moreover, tandem duplication contributes significantly to the adaptive evolution and characteristic diversity of Theaceae, and it is also a major mechanism driving the origination of orphan genes. This study illuminates the evolutionary dynamics of orphan genes, providing a valuable resource for understanding the origin and evolution of tea plant flavor and enhancing genetic breeding efforts.

Domain-Shuffling in the Evolution of Cyclostomes and Gnathostomes

Vertebrates acquired various novel traits that were pivotal in their morphological evolution. Domain shuffling, rearrangements of functional domains between genes, is a key molecular mechanism in deuterostome evolution. However, comprehensive studies focusing on early vertebrates are lacking. With advancements in genomic studies, the genomes of early vertebrate groups and cyclostomes are now accessible, enabling detailed comparative analysis while considering the timing of gene acquisition during evolution. Here, we compared 22 metazoans, including four cyclostomes, to identify genes containing novel domain architectures acquired via domain-shuffling (DSO-Gs), in the common ancestor of vertebrates, gnathostomes, and cyclostomes. We found that DSO-Gs in the common ancestor of vertebrates were associated with novel vertebrate characteristics and those in the common ancestor of gnathostomes correlated with gnathostome-specific traits. Notably, several DSO-Gs acquired in common ancestors of vertebrates have been linked to myelination, a distinct characteristic of gnathostomes. Additionally, in situ hybridization revealed specific expression patterns for the three vertebrate DSO-Gs in cyclostomes, supporting their potential functions. Our findings highlight the significance of DSO-Gs in the emergence of novel traits in the common ancestors of vertebrates, gnathostomes, and cyclostomes.

Gene novelty and gene family expansion in the early evolution of Lepidoptera

BACKGROUND

Almost 10% of all known animal species belong to Lepidoptera: moths and butterflies. To understand how this incredible diversity evolved we assess the role of gene gain in driving early lepidopteran evolution. Here, we compared the complete genomes of 115 insect species, including 99 Lepidoptera, to search for novel genes coincident with the emergence of Lepidoptera.

RESULTS

We find 217 orthogroups or gene families which emerged on the branch leading to Lepidoptera; of these 177 likely arose by gene duplication followed by extensive sequence divergence, 2 are candidates for origin by horizontal gene transfer, and 38 have no known homology outside of Lepidoptera and possibly arose via de novo gene genesis. We focus on two new gene families that are conserved across all lepidopteran species and underwent extensive duplication, suggesting important roles in lepidopteran biology. One encodes a family of sugar and ion transporter molecules, potentially involved in the evolution of diverse feeding behaviours in early Lepidoptera. The second encodes a family of unusual propeller-shaped proteins that likely originated by horizontal gene transfer from Spiroplasma bacteria; we name these the Lepidoptera propellin genes.

CONCLUSION

We provide the first insights into the role of genetic novelty in the early evolution of Lepidoptera. This gives new insight into the rate of gene gain during the evolution of the order as well as providing context on the likely mechanisms of origin. We describe examples of new genes which were retained and duplicated further in all lepidopteran species, suggesting their importance in Lepidoptera evolution.

Genome-wide identification and expression analysis of orphan genes in twelve Musa (sub)species

Orphan genes (OGs), also known as lineage-specific genes, are species-specific genes that play a crucial role in species-specific adaptations to various stresses. Although OGs have been identified in several plant species, there is no information on OGs in banana genomes. This study aimed to systematically identify OGs in twelve banana (sub)species using comparative genomics. The results showed that OG content varied widely among these (sub)species, from 0.4% in Musa itinerans to 7.3% in Ensete glaucum. Genetic structure analysis showed that banana OGs have significantly shorter protein lengths, smaller molecular weight, fewer exons, and shorter exon lengths than non-orphan genes (NOGs). Subcellular localization predictions showed that banana OGs are mainly found in the chloroplast, nucleus, and cytosol, and are evenly distributed across chromosomes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses suggested that OGs may be involved in cellular processes, metabolic processes, and molecular transport. The transcriptome analysis of 9 AAA cultivars against 4 M. acuminata subspecies genomes showed the OGs content. Analysis of gene expression in M. acuminata subsp. malaccensis showed 75 differentially expressed (DE) OGs in response to abiotic stresses and 46 DE OGs related to biotic stresses, indicating that these OGs might play important roles in response to abiotic and biotic stresses. This study provides a foundation for further in-depth research into the functions of OGs in bananas.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-025-04213-9.

Functional innovation through new genes as a general evolutionary process

In the past decade, our understanding of how new genes originate in diverse organisms has advanced substantially, and more than a dozen molecular mechanisms for generating initial gene structures were identified, in addition to gene duplication. These new genes have been found to integrate into and modify pre-existing gene networks primarily through mutation and selection, revealing new patterns and rules with stable origination rates across various organisms. This progress has challenged the prevailing belief that new proteins evolve from pre-existing genes, as new genes may arise de novo from noncoding DNA sequences in many organisms, with high rates observed in flowering plants. New genes have important roles in phenotypic and functional evolution across diverse biological processes and structures, with detectable fitness effects of sexual conflict genes that can shape species divergence. Such knowledge of new genes can be of translational value in agriculture and medicine.

Genome sequences of four Ixodes species expands understanding of tick evolution

BACKGROUND

Ticks, hematophagous Acari, pose a significant threat by transmitting various pathogens to their vertebrate hosts during feeding. Despite advances in tick genomics, high-quality genomes were lacking until recently, particularly in the genus Ixodes, which includes the main vectors of Lyme disease.

RESULTS

Here, we present the genome sequences of four tick species, derived from a single female individual, with a particular focus on the European species Ixodes ricinus, achieving a chromosome-level assembly. Additionally, draft assemblies were generated for the three other Ixodes species, I. persulcatus, I. pacificus, and I. hexagonus. The quality of the four genomes and extensive annotation of several important gene families have allowed us to study the evolution of gene repertoires at the level of the genus Ixodes and of the tick group. We have determined gene families that have undergone major amplifications during the evolution of ticks, while an expression atlas obtained for I. ricinus reveals striking patterns of specialization both between and within gene families. Notably, several gene family amplifications are associated with a proliferation of single-exon genes-most strikingly for fatty acid elongases and sulfotransferases.

CONCLUSIONS

The integration of our data with existing genomes establishes a solid framework for the study of gene evolution, improving our understanding of tick biology. In addition, our work lays the foundations for applied research and innovative control targeting these organisms.

A Zea genus-specific micropeptide controls kernel dehydration in maize

Kernel dehydration rate (KDR) is a crucial production trait that affects mechanized harvesting and kernel quality in maize; however, the underlying mechanisms remain unclear. Here, we identified a quantitative trait locus (QTL), qKDR1, as a non-coding sequence that regulates the expression of qKDR1 REGULATED PEPTIDE GENE (RPG). RPG encodes a 31 amino acid micropeptide, microRPG1, which controls KDR by precisely modulating the expression of two genes, ZmETHYLENE-INSENSITIVE3-like 1 and 3, in the ethylene signaling pathway in the kernels after filling. microRPG1 is a Zea genus-specific micropeptide and originated de novo from a non-coding sequence. Knockouts of microRPG1 result in faster KDR in maize. By contrast, overexpression or exogenous application of the micropeptide shows the opposite effect both in maize and Arabidopsis. Our findings reveal the molecular mechanism of microRPG1 in kernel dehydration and provide an important tool for future crop breeding.

A chromosome-scale genome assembly of the pioneer plant Stylosanthes angustifolia: insights into genome evolution and drought adaptation

BACKGROUND

Drought is a major limiting factor for plant survival and crop productivity. Stylosanthes angustifolia, a pioneer plant, exhibits remarkable drought tolerance, yet the molecular mechanisms driving its drought resistance remain largely unexplored.

RESULTS

We present a chromosome-scale reference genome of S. angustifolia, which provides insights into its genome evolution and drought tolerance mechanisms. The assembled genome is 645.88 Mb in size, containing 319.98 Mb of repetitive sequences and 36,857 protein-coding genes. The high quality of this genome assembly is demonstrated by the presence of 99.26% BUSCO and a 19.49 long terminal repeat assembly index. Evolutionary analyses revealed that S. angustifolia shares a whole-genome duplication (WGD) event with other legumes but lacks recent WGD. Additionally, S. angustifolia has undergone gene expansion through tandem duplication approximately 12.31 million years ago. Through integrative multiomics analyses, we identified 4 gene families-namely, xanthoxin dehydrogenase, 2-hydroxyisoflavanone dehydratase, patatin-related phospholipase A, and stachyose synthetase-that underwent tandem duplication and were significantly upregulated under drought stress. These gene families contribute to the biosynthesis of abscisic acid, genistein, daidzein, jasmonic acid, and stachyose, thereby enhancing drought tolerance.

CONCLUSIONS

The genome assembly of S. angustifolia represents a significant advancement in understanding the genetic mechanisms underlying drought tolerance in this pioneer plant species. This genomic resource provides critical insights into the evolution of drought resistance and offers valuable genetic information for breeding programs aimed at improving drought resistance in crops.

Identification of Retrocopies in Lepidoptera and Impact on Domestication of Silkworm

BACKGROUND

During the domestication of silkworm, an economic insect, its physiological characteristics have changed greatly. RNA-based gene duplication, known as retrocopy, plays an important role in the formation of new genes and genome evolution, but the retrocopies of lepidopteran insects have not been fully identified and analyzed, which not only severely limits researchers from exploring the effects of retrocopies on lepidopteran insects but also affects the studies on the domestication of silkworm.

METHODS

We compared the genomes and proteomes of eight lepidopteran insects and used a series of screening criteria for auxiliary screening to obtain the retrocopies in lepidopteran insects and explored their characteristics. In addition, based on the silkworm transcriptome data from the SilkDB3.0 website, we explored the functions of the retrocopies on the domestication of the silkworm.

RESULTS

A total of 1993 retrocopies and 1208 parental genes in lepidopteran insects were obtained. We revealed that the retrocopies in Lepidoptera do not conform to the "out of X" hypothesis but fit the "out of testis" hypothesis. These retrocopies were subject to strong functional constraints and performed important functions in growth and development. Transcriptome analysis revealed that the expression pattern of the retrocopies and their parental genes were irrelevant. Through the analysis of the retrocopies in silkworm generated after domestication and located in the candidate domestication regions, the possible universal connection between the retrocopies and the domestication of silkworm were found.

CONCLUSIONS

Our study pioneered the exploration of retrocopies in multiple Lepidoptera species and found the potential association between the retrocopies and the domestication of silkworm.

Tissue-specific, temporal, and core gene-dependent expression patterns of Hsp70s reveal functional allocation in Chlamys farreri under heat stress

Heat shock proteins 70 KDa (Hsp70s) engage in a broad spectrum of cellular functions in response to various stressors. Marine bivalves face substantial threats from the rising seawater temperature attributed to global warming. In the present study, expression patterns of Hsp70s in Zhikong scallop Chlamys farreri (CfHsp70s) were determined in embryos and larvae at all developmental stages, in healthy adult tissues, and across four various tissues exposed to high temperature for acute and chronic periods through in silico analysis. Spatiotemporal expressions suggested CfHsp70s performed specific functional differentiations in scallop's development and growth. Regulatory expression patterns of CfHsp70s, characterized by predominant down-regulation in the mantle, gill and hemocytes, as well as contrasting up-regulation in the heart, suggest differential functional allocation of CfHsp70s among tissues in response to heat stress. Particularly, a core set of 14 CfHsp70s, especially the nine members of the Hsp70B2s, characterized by gene expansion, intron-less structure, shorter gene length, preference for hydrophilic amino acids, and coordinated expression profiles, was predominantly responsible for the inducible up-regulations observed across all four tissue types. Collectively, the tissue-specific, temporal and core gene-dependent expression patterns of CfHsp70s illustrate the functional allocation and molecular evolution of Hsp70 family members in Zhikong scallops.

Origins of lineage-specific elements via gene duplication, relocation, and regional rearrangement in Neurospora crassa

The origin of new genes has long been a central interest of evolutionary biologists. However, their novelty means that they evade reconstruction by the classical tools of evolutionary modelling. This evasion of deep ancestral investigation necessitates intensive study of model species within well-sampled, recently diversified, clades. One such clade is the model genus Neurospora, members of which lack recent gene duplications. Several Neurospora species are comprehensively characterized organisms apt for studying the evolution of lineage-specific genes (LSGs). Using gene synteny, we documented that 78% of Neurospora LSG clusters are located adjacent to the telomeres featuring extensive tracts of non-coding DNA and duplicated genes. Here, we report several instances of LSGs that are likely from regional rearrangements and potentially from gene rebirth. To broadly investigate the functions of LSGs, we assembled transcriptomics data from 68 experimental data points and identified co-regulatory modules using Weighted Gene Correlation Network Analysis, revealing that LSGs are widely but peripherally involved in known regulatory machinery for diverse functions. The ancestral status of the LSG mas-1, a gene with roles in cell-wall integrity and cellular sensitivity to antifungal toxins, was investigated in detail alongside its genomic neighbours, indicating that it arose from an ancient lysophospholipase precursor that is ubiquitous in lineages of the Sordariomycetes. Our discoveries illuminate a "rummage region" in the N. crassa genome that enables the formation of new genes and functions to arise via gene duplication and relocation, followed by fast mutation and recombination facilitated by sequence repeats and unconstrained non-coding sequences.

De novo gene integration into regulatory networks via interaction with conserved genes in peach

De novo genes can evolve "from scratch" from noncoding sequences, acquiring novel functions in organisms and integrating into regulatory networks during evolution to drive innovations in important phenotypes and traits. However, identifying de novo genes is challenging, as it requires high-quality genomes from closely related species. According to the comparison with nine closely related Prunus genomes, we determined at least 178 de novo genes in P. persica "baifeng". The distinct differences were observed between de novo and conserved genes in gene characteristics and expression patterns. Gene ontology enrichment analysis suggested that Type I de novo genes originated from sequences related to plastid modification functions, while Type II genes were inferred to have derived from sequences related to reproductive functions. Finally, transcriptome sequencing across different tissues and developmental stages suggested that de novo genes have been evolutionarily recruited into existing regulatory networks, playing important roles in plant growth and development, which was also supported by WGCNA analysis and quantitative trait loci data. This study lays the groundwork for future research on the origins and functions of genes in Prunus and related taxa.

Specialization Restricts the Evolutionary Paths Available to Yeast Sugar Transporters

Functional innovation at the protein level is a key source of evolutionary novelties. The constraints on functional innovations are likely to be highly specific in different proteins, which are shaped by their unique histories and the extent of global epistasis that arises from their structures and biochemistries. These contextual nuances in the sequence-function relationship have implications both for a basic understanding of the evolutionary process and for engineering proteins with desirable properties. Here, we have investigated the molecular basis of novel function in a model member of an ancient, conserved, and biotechnologically relevant protein family. These Major Facilitator Superfamily sugar porters are a functionally diverse group of proteins that are thought to be highly plastic and evolvable. By dissecting a recent evolutionary innovation in an α-glucoside transporter from the yeast Saccharomyces eubayanus, we show that the ability to transport a novel substrate requires high-order interactions between many protein regions and numerous specific residues proximal to the transport channel. To reconcile the functional diversity of this family with the constrained evolution of this model protein, we generated new, state-of-the-art genome annotations for 332 Saccharomycotina yeast species spanning ∼400 My of evolution. By integrating phylogenetic and phenotypic analyses across these species, we show that the model yeast α-glucoside transporters likely evolved from a multifunctional ancestor and became subfunctionalized. The accumulation of additive and epistatic substitutions likely entrenched this subfunction, which made the simultaneous acquisition of multiple interacting substitutions the only reasonably accessible path to novelty.

De Novo Genes

Although the majority of annotated new genes in a given genome appear to have arisen from duplication-related mechanisms, recent studies have shown that genes can also originate de novo from ancestrally nongenic sequences. Investigating de novo-originated genes offers rich opportunities to understand the origin and functions of new genes, their regulatory mechanisms, and the associated evolutionary processes. Such studies have uncovered unexpected and intriguing facets of gene origination, offering novel perspectives on the complexity of the genome and gene evolution. In this review, we provide an overview of the research progress in this field, highlight recent advancements, identify key technical and conceptual challenges, and underscore critical questions that remain to be addressed.

Testis-specific H2B.W1 disrupts nucleosome integrity by reducing DNA-histone interactions

Multiple testis-specific histone variants are involved in the dynamic chromatin transitions during spermatogenesis. H2B.W1 (previously called H2BFWT) is an H2B variant specific to primate testis with hitherto unclear functions, although its single-nucleotide polymorphisms (SNPs) are closely associated with male non-obstructive infertility. Here, we found that H2B.W1 is only expressed in the mid-late spermatogonia stages, and H2B.W1 nucleosomes are defined by a more flexible structure originating from weakened interactions between histones and DNA. Furthermore, one of its SNPs, H2B.W1-H100R, which is associated with infertility, further destabilizes the nucleosomes and increases the nucleosome unwrapping rate by interfering with the R100 and H4 K91/R92 interaction. Our results suggest that destabilizing H2B.W1 containing nucleosomes might change the chromatin structure of spermatogonia, and that H2B.W1-H100R enhances the nucleosome-destabilizing effects, leading to infertility.

Single-cell analysis of the amphioxus hepatic caecum and vertebrate liver reveals genetic mechanisms of vertebrate liver evolution

The evolution of the vertebrate liver is a prime example of the evolution of complex organs, yet the driving genetic factors behind it remain unknown. Here we study the evolutionary genetics of liver by comparing the amphioxus hepatic caecum and the vertebrate liver, as well as examining the functional transition within vertebrates. Using in vivo and in vitro experiments, single-cell/nucleus RNA-seq data and gene knockout experiments, we confirm that the amphioxus hepatic caecum and vertebrate liver are homologous organs and show that the emergence of ohnologues from two rounds of whole-genome duplications greatly contributed to the functional complexity of the vertebrate liver. Two ohnologues, kdr and flt4, play an important role in the development of liver sinusoidal endothelial cells. In addition, we found that liver-related functions such as coagulation and bile production evolved in a step-by-step manner, with gene duplicates playing a crucial role. We reconstructed the genetic footprint of the transfer of haem detoxification from the liver to the spleen during vertebrate evolution. Together, these findings challenge the previous hypothesis that organ evolution is primarily driven by regulatory elements, underscoring the importance of gene duplicates in the emergence and diversification of a complex organ.

Chloroplast genomes of Simarouba Aubl., molecular evolution and comparative analyses within Sapindales

Simarouba, a neotropical genus in the family Simaroubaceae, currently lacks comprehensive genomic data in existing databases. This study aims to fill this gap by providing genomic resources for three Simarouba species, S. amara, S. versicolor, and S. glauca. It also aims to perform comparative molecular evolutionary analyses in relation to other species within the order Sapindales. The analysis of these three Simarouba species revealed the presence of the typical quadripartite structure expected in plastomes. However, some pseudogenization events were identified in the psbC, infA, rpl22, and ycf1 genes. In particular, the CDS of the psbC gene in S. amara was reduced from 1422 bp to 584 bp due to a premature stop codon. Nucleotide diversity data pointed to gene and intergenic regions as promising candidates for species and family discrimination within the group, specifically matK, ycf1, ndhF, rpl32, petA-psbJ, and trnS-trnG. Selection signal analyses showed strong evidence for positive selection on the rpl23 gene. Phylogenetic analyses indicated that S. versicolor and S. glauca have a closer phylogenetic relationship than S. amara. We provide chloroplast genomes of three Simaruba species and use them to elucidate plastome evolution, highlight the presence of pseudogenization, and identify potential DNA barcode regions.

New genes driven by segmental duplications share a testis-specific expression pattern in the chromosome-level genome assembly of tree sparrow

Based on a chromosome-level genome assembly, a burst of new genes with different structures but a similar testis-specific expression pattern was detected in tree sparrow.

The role of heterochronic gene expression and regulatory architecture in early developmental divergence

New developmental programs can evolve through adaptive changes to gene expression. The annelid Streblospio benedicti has a developmental dimorphism, which provides a unique intraspecific framework for understanding the earliest genetic changes that take place during developmental divergence. Using comparative RNAseq through ontogeny, we find that only a small proportion of genes are differentially expressed at any time, despite major differences in larval development and life history. These genes shift expression profiles across morphs by either turning off any expression in one morph or changing the timing or amount of gene expression. We directly connect the contributions of these mechanisms to differences in developmental processes. We examine F1 offspring - using reciprocal crosses - to determine maternal mRNA inheritance and the regulatory architecture of gene expression. These results highlight the importance of both novel gene expression and heterochronic shifts in developmental evolution, as well as the trans-acting regulatory factors in initiating divergence.

Genomic and Transcriptomic Insights into the Evolution of C4 Photosynthesis in Grasses

C4 photosynthesis has independently evolved over 62 times within 19 angiosperm families. The recurrent evolution of C4 photosynthesis appears to contradict the complex anatomical and biochemical modifications required for the transition from C3 to C4 photosynthesis. In this study, we conducted an integrated analysis of genomics and transcriptomics to elucidate the molecular underpinnings of convergent C4 evolution in the grass family. Our genome-wide exploration of C4-related gene families suggests that the expansion of these gene families may have played an important role in facilitating C4 evolution in the grass family. A phylogenomic synteny network analysis uncovered the emergence of C4 genes in various C4 grass lineages from a common ancestral gene pool. Moreover, through a comparison between non-C4 and C4  PEPCs, we pinpointed 14 amino acid sites exhibiting parallel adaptations. These adaptations, occurring post the BEP-PACMAD divergence, shed light on why all C4 origins in grasses are confined to the PACMAD clade. Furthermore, our study revealed that the ancestor of Chloridoideae grasses possessed a more favorable molecular preadaptation for C4 functions compared to the ancestor of Panicoideae grasses. This molecular preadaptation potentially explains why C4 photosynthesis evolved earlier in Chloridoideae than in Panicoideae and why the C3-to-C4 transition occurred once in Chloridoideae but multiple times in Panicoideae. Additionally, we found that C4 genes share similar cis-elements across independent C4 lineages. Notably, NAD-ME subtype grasses may have retained the ancestral regulatory machinery of the C4  NADP-ME gene, while NADP-ME subtype grasses might have undergone unique cis-element modifications.

Adapting the rhizome concept to an extended definition of viral quasispecies and the implications for molecular evolution

The rhizome concept proposed by Gilles Deleuze and Félix Guattari offers a novel perspective on the organization and interdependence of complex constellations of heterogeneous entities, their mapping and their ruptures. The emphasis of the present study is placed on the dynamics of contacts and communication among such entities that arise from experimentation, without any favored hierarchy or origin. When applied to biological evolution, the rhizome concept integrates all types of heterogeneity resulting from "symbiotic" relationships among living beings (or their genomic material), horizontal genetic transfer, recombination and mutation, and breaks away from the approach that gives rise to the phylogenetic tree of life. It has already been applied to describe the dynamics and evolution of RNA viruses. Thus, here we introduce a novel framework for the interpretation the viral quasispecies concept, which explains the evolution of RNA virus populations as the result of dynamic interconnections and multifaceted interdependence between highly heterogeneous viral sequences and its inherently heterogeneous host cells. The rhizome network perspective underlines even further the medical implications of the broad mutant spectra of viruses that are in constant flow, given the multiple pathways they have available for fitness loss and gain.

Characterization of GPX Gene Family in Pepper (Capsicum annuum L.) under Abiotic Stress and ABA Treatment

Plant glutathione peroxidases (GPXs) are important enzymes for removing reactive oxygen species in plant cells and are closely related to the stress resistance of plants. This study identified the GPX gene family members of pepper (Capsicum annuum L.), "CM333", at the whole-genome level to clarify their expression patterns and enzyme activity changes under abiotic stress and ABA treatment. The results showed that eight CaGPX genes were unevenly distributed across four chromosomes and one scaffold of the pepper genome, and their protein sequences had Cys residues typical of the plant GPX domains. The analysis of collinearity, phylogenetic tree, gene structure, and conserved motifs indicated that the CaGPX gene sequence is conserved, structurally similar, and more closely related to the sequence structure of Arabidopsis. Meanwhile, many cis elements involved in stress, hormones, development, and light response were found in the promoter region of the CaGPX gene. In addition, CaGPX1/4 and CaGPX6 were basically expressed in all tissues, and their expression levels were significantly upregulated under abiotic stress and ABA treatment. Subcellular localization showed that CaGPX1 and CaGPX4 are localized in chloroplasts. Additionally, the variations in glutathione peroxidase activity (GSH-Px) mostly agreed with the variations in gene expression. In summary, the CaGPXs gene may play an important role in the development of peppers and their response to abiotic stress and hormones.

How antisense transcripts can evolve to encode novel proteins

Protein coding features can emerge de novo in non coding transcripts, resulting in emergence of new protein coding genes. Studies across many species show that a large fraction of evolutionarily novel non-coding RNAs have an antisense overlap with protein coding genes. The open reading frames (ORFs) in these antisense RNAs could also overlap with existing ORFs. In this study, we investigate how the evolution an ORF could be constrained by its overlap with an existing ORF in three different reading frames. Using a combination of mathematical modeling and genome/transcriptome data analysis in two different model organisms, we show that antisense overlap can increase the likelihood of ORF emergence and reduce the likelihood of ORF loss, especially in one of the three reading frames. In addition to rationalising the repeatedly reported prevalence of de novo emerged genes in antisense transcripts, our work also provides a generic modeling and an analytical framework that can be used to understand evolution of antisense genes.

Diverse origins of near-identical antifreeze proteins in unrelated fish lineages provide insights into evolutionary mechanisms of new gene birth and protein sequence convergence

Determining the origins of novel genes and the genetic mechanisms underlying the emergence of new functions is challenging yet crucial for understanding evolutionary innovations. The convergently evolved fish antifreeze proteins provide excellent opportunities to investigate evolutionary origins and pathways of new genes. Particularly notable is the near-identical type I antifreeze proteins (AFPI) in four phylogenetically divergent fish taxa. This study tested the hypothesis of protein sequence convergence beyond functional convergence in three unrelated AFPI-bearing fish lineages, revealing different paths by which a similar protein arose from diverse genomic resources. Comprehensive comparative analyses of de novo sequenced genome of the winter flounder and grubby sculpin, available high-quality genome of the cunner and 14 other relevant species found that the near-identical AFPI originated from a distinct genetic precursor in each lineage. Each independently evolved a coding region for the novel ice-binding protein while retaining sequence identity in the regulatory regions with their respective ancestor. The deduced evolutionary processes and molecular mechanisms are consistent with the Innovation-Amplification-Divergence (IAD) model applicable to AFPI formation in all three lineages, a new Duplication-Degeneration-Divergence (DDD) model we propose for the sculpin lineage, and a DDD model with gene fission for the cunner lineage. This investigation illustrates the multiple ways by which a novel functional gene with sequence convergence at the protein level could evolve across divergent species, advancing our understanding of the mechanistic intricacies in new gene formation.

Cytochrome P450 Enzyme Design by Constraining the Catalytic Pocket in a Diffusion Model

Although cytochrome P450 enzymes are the most versatile biocatalysts in nature, there is insufficient comprehension of the molecular mechanism underlying their functional innovation process. Here, by combining ancestral sequence reconstruction, reverse mutation assay, and progressive forward accumulation, we identified 5 founder residues in the catalytic pocket of flavone 6-hydroxylase (F6H) and proposed a "3-point fixation" model to elucidate the functional innovation mechanisms of P450s in nature. According to this design principle of catalytic pocket, we further developed a de novo diffusion model (P450Diffusion) to generate artificial P450s. Ultimately, among the 17 non-natural P450s we generated, 10 designs exhibited significant F6H activity and 6 exhibited a 1.3- to 3.5-fold increase in catalytic capacity compared to the natural CYP706X1. This work not only explores the design principle of catalytic pockets of P450s, but also provides an insight into the artificial design of P450 enzymes with desired functions.

PHACTboost: A Phylogeny-Aware Pathogenicity Predictor for Missense Mutations via Boosting

Most algorithms that are used to predict the effects of variants rely on evolutionary conservation. However, a majority of such techniques compute evolutionary conservation by solely using the alignment of multiple sequences while overlooking the evolutionary context of substitution events. We had introduced PHACT, a scoring-based pathogenicity predictor for missense mutations that can leverage phylogenetic trees, in our previous study. By building on this foundation, we now propose PHACTboost, a gradient boosting tree-based classifier that combines PHACT scores with information from multiple sequence alignments, phylogenetic trees, and ancestral reconstruction. By learning from data, PHACTboost outperforms PHACT. Furthermore, the results of comprehensive experiments on carefully constructed sets of variants demonstrated that PHACTboost can outperform 40 prevalent pathogenicity predictors reported in the dbNSFP, including conventional tools, metapredictors, and deep learning-based approaches as well as more recent tools such as AlphaMissense, EVE, and CPT-1. The superiority of PHACTboost over these methods was particularly evident in case of hard variants for which different pathogenicity predictors offered conflicting results. We provide predictions of 215 million amino acid alterations over 20,191 proteins. PHACTboost is available at https://github.com/CompGenomeLab/PHACTboost. PHACTboost can improve our understanding of genetic diseases and facilitate more accurate diagnoses.

An Orphan Gene Enhances Male Reproductive Success in Plutella xylostella

Plutella xylostella exhibits exceptional reproduction ability, yet the genetic basis underlying the high reproductive capacity remains unknown. Here, we demonstrate that an orphan gene, lushu, which encodes a sperm protein, plays a crucial role in male reproductive success. Lushu is located on the Z chromosome and is prevalent across different P. xylostella populations worldwide. We subsequently generated lushu mutants using transgenic CRISPR/Cas9 system. Knockout of Lushu results in reduced male mating efficiency and accelerated death in adult males. Furthermore, our findings highlight that the deficiency of lushu reduced the transfer of sperms from males to females, potentially resulting in hindered sperm competition. Additionally, the knockout of Lushu results in disrupted gene expression in energy-related pathways and elevated insulin levels in adult males. Our findings reveal that male reproductive performance has evolved through the birth of a newly evolved, lineage-specific gene with enormous potentiality in fecundity success. These insights hold valuable implications for identifying the target for genetic control, particularly in relation to species-specific traits that are pivotal in determining high levels of fecundity.

The Rapid Evolution of De Novo Proteins in Structure and Complex

Recent studies in the rice genome-wide have established that de novo genes, evolving from noncoding sequences, enhance protein diversity through a stepwise process. However, the pattern and rate of their evolution in protein structure over time remain unclear. Here, we addressed these issues within a surprisingly short evolutionary timescale (<1 million years for 97% of Oryza de novo genes) with comparative approaches to gene duplicates. We found that de novo genes evolve faster than gene duplicates in the intrinsically disordered regions (such as random coils), secondary structure elements (such as α helix and β strand), hydrophobicity, and molecular recognition features. In de novo proteins, specifically, we observed an 8% to 14% decay in random coils and intrinsically disordered region lengths and a 2.3% to 6.5% increase in structured elements, hydrophobicity, and molecular recognition features, per million years on average. These patterns of structural evolution align with changes in amino acid composition over time as well. We also revealed higher positive charges but smaller molecular weights for de novo proteins than duplicates. Tertiary structure predictions showed that most de novo proteins, though not typically well folded on their own, readily form low-energy and compact complexes with other proteins facilitated by extensive residue contacts and conformational flexibility, suggesting a faster-binding scenario in de novo proteins to promote interaction. These analyses illuminate a rapid evolution of protein structure in de novo genes in rice genomes, originating from noncoding sequences, highlighting their quick transformation into active, protein complex-forming components within a remarkably short evolutionary timeframe.

AcMYB266, a key regulator of the red coloration in pineapple peel: a case of subfunctionalization in tandem duplicated genes

Red fruit peel is an attractive target for pineapple breeding. Various pineapple accessions with distinct red coloration patterns exist; however, the precise molecular mechanism accounting for these differences remains unknown, which hinders the pineapple breeding process from combining high fruit quality with red peel. In this study, we characterized a transcription factor, AcMYB266, which is preferentially expressed in pineapple peel and positively regulates anthocyanin accumulation. Transgenic pineapple, Arabidopsis, and tobacco plants overexpressing AcMYB266 exhibited significant anthocyanin accumulation. Conversely, transient silencing of this gene led to decreased anthocyanin accumulation in pineapple red bracts. In-depth analysis indicated that variations of AcMYB266 sequences in the promoter instead of the protein-coding region seem to contribute to different red coloration patterns in peels of three representative pineapple varieties. In addition, we found that AcMYB266 was located in a cluster of four MYB genes exclusive to and conserved in Ananas species. Of this cluster, each was proved to regulate anthocyanin synthesis in different pineapple tissues, illustrating an interesting case of gene subfunctionalization after tandem duplication. In summary, we have characterized AcMYB266 as a key regulator of pineapple red fruit peel and identified an MYB cluster whose members were subfunctionalized to specifically regulate the red coloration of different pineapple tissues. The present study will assist in establishing a theoretical mechanism for pineapple breeding for red fruit peel and provide an interesting case for the investigation of gene subfunctionalization in plants.

Genome-Wide Analysis of BBX Gene Family in Three Medicago Species Provides Insights into Expression Patterns under Hormonal and Salt Stresses

BBX protein is a class of zinc finger transcription factors that have B-box domains at the N-terminus, and some of these proteins contain a CCT domain at the C-terminus. It plays an important role in plant growth, development, and metabolism. However, the expression pattern of BBX genes in alfalfa under hormonal and salt stresses is still unclear. In this study, we identified a total of 125 BBX gene family members by the available Medicago reference genome in diploid alfalfa (Medicago sativa spp. Caerulea), a model plant (M. truncatula), and tetraploid alfalfa (M. sativa), and divided these members into five subfamilies. We found that the conserved motifs of BBXs of the same subfamily reveal similarities. We analyzed the collinearity relationship and duplication mode of these BBX genes and found that the expression pattern of BBX genes is specific in different tissues. Analysis of the available transcriptome data suggests that some members of the BBX gene family are involved in multiple abiotic stress responses, and the highly expressed genes are often clustered together. Furthermore, we identified different expression patterns of some BBX genes under salt, ethylene, salt and ethylene, salicylic acid, and salt and salicylic acid treatments, verified by qRT-PCR, and analyzed the subcellular localization of MsBBX2, MsBBX17, and MsBBX32 using transient expression in tobacco. The results showed that BBX genes were localized in the nucleus. This study systematically analyzed the BBX gene family in Medicago plants, which provides a basis for the study of BBX gene family tolerance to abiotic stresses.

Plastid DNA is a major source of nuclear genome complexity and of RNA genes in the orphan crop moringa

BACKGROUND

Unlike Transposable Elements (TEs) and gene/genome duplication, the role of the so-called nuclear plastid DNA sequences (NUPTs) in shaping the evolution of genome architecture and function remains poorly studied. We investigate here the functional and evolutionary fate of NUPTs in the orphan crop Moringa oleifera (moringa), featured by the highest fraction of plastid DNA found so far in any plant genome, focusing on (i) any potential biases in their distribution in relation to specific nuclear genomic features, (ii) their contribution to the emergence of new genes and gene regions, and (iii) their impact on the expression of target nuclear genes.

RESULTS

In agreement with their potential mutagenic effect, NUPTs are underrepresented among structural genes, although their overall transcription levels and broadness were only lower when involved exonic regions; the occurrence of plastid DNA generally did not result in a broader expression, except among those affected in introns by older NUPTs. In contrast, we found a strong enrichment of NUPTs among specific superfamilies of retrotransposons and several classes of RNA genes, including those participating in the protein biosynthetic machinery (i.e., rRNA and tRNA genes) and a specific class of regulatory RNAs. A significant fraction of NUPT RNA genes was found to be functionally expressed, thus potentially contributing to the nuclear pool.

CONCLUSIONS

Our results complete our view of the molecular factors driving the evolution of nuclear genome architecture and function, and support plastid DNA in moringa as a major source of (i) genome complexity and (ii) the nuclear pool of RNA genes.

Unveiling the Genetic Blueprint of a Desert Scorpion: A Chromosome-level Genome of Hadrurus arizonensis Provides the First Reference for Parvorder Iurida

Over 400 million years old, scorpions represent an ancient group of arachnids and one of the first animals to adapt to life on land. Presently, the lack of available genomes within scorpions hinders research on their evolution. This study leverages ultralong nanopore sequencing and Pore-C to generate the first chromosome-level assembly and annotation for the desert hairy scorpion, Hadrurus arizonensis. The assembled genome is 2.23 Gb in size with an N50 of 280 Mb. Pore-C scaffolding reoriented 99.6% of bases into nine chromosomes and BUSCO identified 998 (98.6%) complete arthropod single copy orthologs. Repetitive elements represent 54.69% of the assembled bases, including 872,874 (29.39%) LINE elements. A total of 18,996 protein-coding genes and 75,256 transcripts were predicted, and extracted protein sequences yielded a BUSCO score of 97.2%. This is the first genome assembled and annotated within the family Hadruridae, representing a crucial resource for closing gaps in genomic knowledge of scorpions, resolving arachnid phylogeny, and advancing studies in comparative and functional genomics.

A newly evolved rice-specific gene JAUP1 regulates jasmonate biosynthesis and signalling to promote root development and multi-stress tolerance

Root architecture and function are critical for plants to secure water and nutrient supply from the soil, but environmental stresses alter root development. The phytohormone jasmonic acid (JA) regulates plant growth and responses to wounding and other stresses, but its role in root development for adaptation to environmental challenges had not been well investigated. We discovered a novel JA Upregulated Protein 1 gene (JAUP1) that has recently evolved in rice and is specific to modern rice accessions. JAUP1 regulates a self-perpetuating feed-forward loop to activate the expression of genes involved in JA biosynthesis and signalling that confers tolerance to abiotic stresses and regulates auxin-dependent root development. Ectopic expression of JAUP1 alleviates abscisic acid- and salt-mediated suppression of lateral root (LR) growth. JAUP1 is primarily expressed in the root cap and epidermal cells (EPCs) that protect the meristematic stem cells and emerging LRs. Wound-activated JA/JAUP1 signalling promotes crosstalk between the root cap of LR and parental root EPCs, as well as induces cell wall remodelling in EPCs overlaying the emerging LR, thereby facilitating LR emergence even under ABA-suppressive conditions. Elevated expression of JAUP1 in transgenic rice or natural rice accessions enhances abiotic stress tolerance and reduces grain yield loss under a limited water supply. We reveal a hitherto unappreciated role for wound-induced JA in LR development under abiotic stress and suggest that JAUP1 can be used in biotechnology and as a molecular marker for breeding rice adapted to extreme environmental challenges and for the conservation of water resources.

Evolution of a Human-Specific De Novo Open Reading Frame and Its Linked Transcriptional Silencer

In the human genome, two short open reading frames (ORFs) separated by a transcriptional silencer and a small intervening sequence stem from the gene SMIM45. The two ORFs show different translational characteristics, and they also show divergent patterns of evolutionary development. The studies presented here describe the evolution of the components of SMIM45. One ORF consists of an ultra-conserved 68 amino acid (aa) sequence, whose origins can be traced beyond the evolutionary age of divergence of the elephant shark, ~462 MYA. The silencer also has ancient origins, but it has a complex and divergent pattern of evolutionary formation, as it overlaps both at the 68 aa ORF and the intervening sequence. The other ORF consists of 107 aa. It develops during primate evolution but is found to originate de novo from an ancestral non-coding genomic region with root origins within the Afrothere clade of placental mammals, whose evolutionary age of divergence is ~99 MYA. The formation of the complete 107 aa ORF during primate evolution is outlined, whereby sequence development is found to occur through biased mutations, with disruptive random mutations that also occur but lead to a dead-end. The 107 aa ORF is of particular significance, as there is evidence to suggest it is a protein that may function in human brain development. Its evolutionary formation presents a view of a human-specific ORF and its linked silencer that were predetermined in non-primate ancestral species. The genomic position of the silencer offers interesting possibilities for the regulation of transcription of the 107 aa ORF. A hypothesis is presented with respect to possible spatiotemporal expression of the 107 aa ORF in embryonic tissues.

Genome-wide identification and characterization of protein phosphatase 2C (PP2C) gene family in sunflower (Helianthus annuus L.) and their expression profiles in response to multiple abiotic stresses

Plant protein phosphatase 2C (PP2C) plays vital roles in responding to various stresses, stimulating growth factors, phytohormones, and metabolic activities in many important plant species. However, the PP2C gene family has not been investigated in the economically valuable plant species sunflower (Helianthus annuus L.). This study used comprehensive bioinformatics tools to identify and characterize the PP2C gene family members in the sunflower genome (H. annuus r1.2). Additionally, we analyzed the expression profiles of these genes using RNA-seq data under four different stress conditions in both leaf and root tissues. A total of 121 PP2C genes were identified in the sunflower genome distributed unevenly across the 17 chromosomes, all containing the Type-2C phosphatase domain. HanPP2C genes are divided into 15 subgroups (A-L) based on phylogenetic tree analysis. Analyses of conserved domains, gene structures, and motifs revealed higher structural and functional similarities within various subgroups. Gene duplication and collinearity analysis showed that among the 53 HanPP2C gene pairs, 48 demonstrated segmental duplications under strong purifying selection pressure, with only five gene pairs showing tandem duplications. The abundant segmental duplication was observed compared to tandem duplication, which was the major factor underlying the dispersion of the PP2C gene family in sunflowers. Most HanPP2C proteins were localized in the nucleus, cytoplasm, and chloroplast. Among the 121 HanPP2C genes, we identified 71 miRNAs targeting 86 HanPP2C genes involved in plant developmental processes and response to abiotic stresses. By analyzing cis-elements, we identified 63 cis-regulatory elements in the promoter regions of HanPP2C genes associated with light responsiveness, tissue-specificity, phytohormone, and stress responses. Based on RNA-seq data from two sunflower tissues (leaf and root), 47 HanPP2C genes exhibited varying expression levels in leaf tissue, while 49 HanPP2C genes showed differential expression patterns in root tissue across all stress conditions. Transcriptome profiling revealed that nine HanPP2C genes (HanPP2C12, HanPP2C36, HanPP2C38, HanPP2C47, HanPP2C48, HanPP2C53, HanPP2C54, HanPP2C59, and HanPP2C73) exhibited higher expression in leaf tissue, and five HanPP2C genes (HanPP2C13, HanPP2C47, HanPP2C48, HanPP2C54, and HanPP2C95) showed enhanced expression in root tissue in response to the four stress treatments, compared to the control conditions. These results suggest that these HanPP2C genes may be potential candidates for conferring tolerance to multiple stresses and further detailed characterization to elucidate their functions. From these candidates, 3D structures were predicted for six HanPP2C proteins (HanPP2C47, HanPP2C48, HanPP2C53, HanPP2C54, HanPP2C59, and HanPP2C73), which provided satisfactory models. Our findings provide valuable insights into the PP2C gene family in the sunflower genome, which could play a crucial role in responding to various stresses. This information can be exploited in sunflower breeding programs to develop improved cultivars with increased abiotic stress tolerance.

The genome of Haberlea rhodopensis provides insights into the mechanisms for tolerance to multiple extreme environments

Haberlea rhodopensis, a resurrection species, is the only plant known to be able to survive multiple extreme environments, including desiccation, freezing temperatures, and long-term darkness. However, the molecular mechanisms underlying tolerance to these stresses are poorly studied. Here, we present a high-quality genome of Haberlea and found that ~ 23.55% of the 44,306 genes are orphan. Comparative genomics analysis identified 89 significantly expanded gene families, of which 25 were specific to Haberlea. Moreover, we demonstrated that Haberlea preserves its resurrection potential even in prolonged complete darkness. Transcriptome profiling of plants subjected to desiccation, darkness, and low temperatures revealed both common and specific footprints of these stresses, and their combinations. For example, PROTEIN PHOSPHATASE 2C (PP2C) genes were substantially induced in all stress combinations, while PHYTOCHROME INTERACTING FACTOR 1 (PIF1) and GROWTH RESPONSE FACTOR 4 (GRF4) were induced only in darkness. Additionally, 733 genes with unknown functions and three genes encoding transcription factors specific to Haberlea were specifically induced/repressed upon combination of stresses, rendering them attractive targets for future functional studies. The study provides a comprehensive understanding of the genomic architecture and reports details of the mechanisms of multi-stress tolerance of this resurrection species that will aid in developing strategies that allow crops to survive extreme and multiple abiotic stresses.

Origin of functional de novo genes in humans from "hopeful monsters"

For a long time, it was believed that new genes arise only from modifications of preexisting genes, but the discovery of de novo protein-coding genes that originated from noncoding DNA regions demonstrates the existence of a "motherless" origination process for new genes. However, the features, distributions, expression profiles, and origin modes of these genes in humans seem to support the notion that their origin is not a purely "motherless" process; rather, these genes arise preferentially from genomic regions encoding preexisting precursors with gene-like features. In such a case, the gene loci are typically not brand new. In this short review, we will summarize the definition and features of human de novo genes and clarify their process of origination from ancestral non-coding genomic regions. In addition, we define the favored precursors, or "hopeful monsters," for the origin of de novo genes and present a discussion of the functional significance of these young genes in brain development and tumorigenesis in humans. This article is categorized under: RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution.

Evolution of chemosensory tissues and cells across ecologically diverse Drosophilids

Chemosensory tissues exhibit significant between-species variability, yet the evolution of gene expression and cell types underlying this diversity remain poorly understood. To address these questions, we conducted transcriptomic analyses of five chemosensory tissues from six Drosophila species and integrated the findings with single-cell datasets. While stabilizing selection predominantly shapes chemosensory transcriptomes, thousands of genes in each tissue have evolved expression differences. Genes that have changed expression in one tissue have often changed in multiple other tissues but at different past epochs and are more likely to be cell type-specific than unchanged genes. Notably, chemosensory-related genes have undergone widespread expression changes, with numerous species-specific gains/losses including novel chemoreceptors expression patterns. Sex differences are also pervasive, including a D. melanogaster-specific excess of male-biased expression in sensory and muscle cells in its forelegs. Together, our analyses provide new insights for understanding evolutionary changes in chemosensory tissues at both global and individual gene levels.

Divergent Evolutionary Rates of Primate Brain Regions as Revealed by Genomics and Transcriptomics

Although the primate brain contains numerous functionally distinct structures that have experienced diverse genetic changes during the course of evolution and development, these changes remain to be explored in detail. Here we utilize two classic metrics from evolutionary biology, the evolutionary rate index (ERI) and the transcriptome age index (TAI), to investigate the evolutionary alterations that have occurred in each area and developmental stage of the primate brain. We observed a higher evolutionary rate for those genes expressed in the non-cortical areas during primate evolution, particularly in human, with the highest rate of evolution being exhibited at brain developmental stages between late infancy and early childhood. Further, the transcriptome age of the non-cortical areas was lower than that of the cerebral cortex, with the youngest age apparent at brain developmental stages between late infancy and early childhood. Our exploration of the evolutionary patterns manifest in each brain area and developmental stage provides important reference points for further research into primate brain evolution.

Germline HPF1 retrogene insertion in RB1 gene involved in cancer predisposition

About half of the human genome is composed of repeated sequences derived from mobile elements, mainly retrotransposons, generally without pathogenic effect. Familial forms of retinoblastoma are caused by germline pathogenic variants in RB1 gene. Here, we describe a family with retinoblastoma affecting a father and his son. No pathogenic variant was identified after DNA analysis of RB1 gene coding sequence and exon-intron junctions. However, RB1 mRNA analysis showed a chimeric transcript with insertion of 114 nucleotides from HPF1 gene inside RB1 gene. This chimeric transcript led to an insertion of 38 amino acids in functional domain of retinoblastoma protein. Subsequent DNA analysis in RB1 intron 17 revealed the presence of a full-length HPF1 retrogene insertion in opposite orientation. Functional assay shows that this insertion has a deleterious impact on retinoblastoma protein function. This is the first report of a full-length retrogene insertion involved in human Mendelian disease leading to a chimeric transcript and a non-functional chimeric protein. Some retrogene insertions may be missed by standard diagnostic genetic testing, so contribution of retrogene insertions to human disease may be underestimated. The increasing use of whole genome sequencing in diagnostic settings will help to get a more comprehensive view of retrogenes.

Divergent Retention of Sucrose Metabolism Genes after Whole Genome Triplication in the Tomato (Solanum lycopersicum)

Sucrose, the primary carbon transport mode and vital carbohydrate for higher plants, significantly impacts plant growth, development, yield, and quality formation. Its metabolism involves three key steps: synthesis, transport, and degradation. Two genome triplication events have occurred in Solanaceae, which have resulted in massive gene loss. In this study, a total of 48 and 65 genes from seven sucrose metabolism gene families in Vitis vinifera and Solanum lycopersicum were identified, respectively. The number of members comprising the different gene families varied widely. And there were significant variations in the pattern of gene duplication and loss in the tomato following two WGD events. Tandem duplication is a major factor in the expansion of the SWEET and Acid INV gene families. All the genes are irregularly distributed on the chromosomes, with the majority of the genes showing collinearity with the grape, particularly the CIN family. And the seven gene families were subjected to a purifying selection. The expression patterns of the different gene families exhibited notable variations. This study presents basic information about the sucrose metabolism genes in the tomato and grape, and paves the way for further investigations into the impact of SCT events on the phylogeny, gene retention duplication, and function of sucrose metabolism gene families in the tomato or Solanaceae, and the adaptive evolution of the tomato.

Analysis of Plant-Specific ANTH Domain-Containing Protein in Marchantia polymorpha

Membrane trafficking is a fundamental mechanism for protein and lipid transport in eukaryotic cells and exhibits marked diversity among eukaryotic lineages with distinctive body plans and lifestyles. Diversification of the membrane trafficking system is associated with the expansion and secondary loss of key machinery components, including RAB GTPases, soluble N-ethylmaleimide-sensitive factor attachment protein receptors (SNAREs) and adaptor proteins, during plant evolution. The number of AP180 N-terminal homology (ANTH) proteins, an adaptor family that regulates vesicle formation and cargo sorting during clathrin-mediated endocytosis, increases during plant evolution. In the genome of Arabidopsis thaliana, 18 genes for ANTH proteins have been identified, a higher number than that in yeast and animals, suggesting a distinctive diversification of ANTH proteins. Conversely, the liverwort Marchantia polymorpha possesses a simpler repertoire; only two genes encoding canonical ANTH proteins have been identified in its genome. Intriguingly, a non-canonical ANTH protein is encoded in the genome of M. polymorpha, which also harbors a putative kinase domain. Similar proteins have been detected in sporadic lineages of plants, suggesting their ancient origin and multiple secondary losses during evolution. We named this unique ANTH group phosphatidylinositol-binding clathrin assembly protein-K (PICALM-K) and characterized it in M. polymorpha using genetic, cell biology-based and artificial intelligence (AI)-based approaches. Our results indicate a flagella-related function of MpPICALM-K in spermatozoids, which is distinct from that of canonical ANTH proteins. Therefore, ANTH proteins have undergone significant functional diversification during evolution, and PICALM-K represents a plant-unique ANTH protein that is delivered by neofunctionalization through exon shuffling.

Generation of de novo miRNAs from template switching during DNA replication

The mechanisms generating novel genes and genetic information are poorly known, even for microRNA (miRNA) genes with an extremely constrained design. All miRNA primary transcripts need to fold into a stem-loop structure to yield short gene products ([Formula: see text]22 nt) that bind and repress their mRNA targets. While a substantial number of miRNA genes are ancient and highly conserved, short secondary structures coding for entirely novel miRNA genes have been shown to emerge in a lineage-specific manner. Template switching is a DNA-replication-related mutation mechanism that can introduce complex changes and generate perfect base pairing for entire hairpin structures in a single event. Here, we show that the template-switching mutations (TSMs) have participated in the emergence of over 6,000 suitable hairpin structures in the primate lineage to yield at least 18 new human miRNA genes, that is 26% of the miRNAs inferred to have arisen since the origin of primates. While the mechanism appears random, the TSM-generated miRNAs are enriched in introns where they can be expressed with their host genes. The high frequency of TSM events provides raw material for evolution. Being orders of magnitude faster than other mechanisms proposed for de novo creation of genes, TSM-generated miRNAs enable near-instant rewiring of genetic information and rapid adaptation to changing environments.