Retention of duplicated genes in evolution

Pot-pollen DNA barcoding as a tool to determine the diversity of plant species visited by Ecuadorian stingless bees

Identifying the main species of plants from where Ecuadorian stingless bees collect pollen is one of the key objectives of management and conservation improvement for these insects. This study aims to determine the botanical origin of pot-pollen using two barcodes, comparing two methodologies (DNA barcoding versus electron microscopy and morphometric tools) and determine the genus and species of pollen source plants of the main honey-producing stingless bees in Ecuador. As main results, Prockia crucis, Coffea canephora, Miconia nervosa, Miconia notabilis, Laurus nobilis, Cecropia ficifolia, Theobroma sp., Artocarpus sp., Croton sp., Euphorbia sp., Mikania sp., and Ophryosporus sp., were the genera and species with the highest presence in the nests (n = 35) of three genera of stingless bees of two provinces located in different climatic regions inside the continental Ecuador. Plant species richness in both areas was statistically similar (p-value = 0.21). We concluded that floral sources' molecular identification with the ITS2 region had a higher number of genera and species detected, than the rbcL gene and microscopy tools, for the Ecuadorian landscapes. We confirmed that the foraging behavior of Melipona sp., Scaptotrigona sp., and Tetragonisca sp., could include non-native flora (27%, 12/44 identifications) that provide a rich source of pollen. Stingless beekeepers could use this information to create flower calendars and establish a schedule for better management of stingless bees in secondary and modified environments.

Transcriptome analysis reveals the mechanism of tolerance to copper toxicity in the white rot fungus Trametes hirsuta AH28-2

Heavy metals, such as copper (Cu), are prevalent in the environment and pose a substantial threat to human health. White rot fungi, especially Trametes spp., display prominent Cu tolerance and removal capacity. However, how Trametes responds to environmental Cu stress remains poorly understood. Here, we found that Trametes hirsuta AH28-2 exhibits Cu removal efficiencies varying from 80.8 % at 1.25 mg/L to 57.6 % at 37 mg/L. Comparative transcriptome analysis identified 812, 1898, and 2110 differentially expressed genes (DEGs) at the Cu concentrations of 1.25, 12.5, and 25 mg/L, respectively. Some DEGs were associated with antioxidant defense systems, secondary metabolite biosynthesis (terpenoids and polyketides), transmembrane transport, and glutathione metabolism, potentially enhancing Cu tolerance. The activities of antioxidant enzymes such as superoxide dismutase, catalase, and laccase were increased under Cu stress. qRT-PCR confirmed the alterations in gene expression and demonstrated that glutathione S-transferases, catalases, cytochrome P450s, and laccases were involved in counteracting Cu-induced stress. Gene silencing experiments further confirmed the crucial roles of laccases in this process. Many transcription factors were enriched under Cu stress, including the Zn2Cys6 family transcription factor GME8421_g (TH8421), which was significantly upregulated at the Cu concentration of 12.5 mg/L. ChIP-seq identified five antioxidant enzyme-encoding genes as direct targets of TH8421, forming a regulatory network that protects against Cu stress. These findings offer insights into the molecular mechanisms driving Cu toxicity tolerance in Trametes fungi.

Evolutional, expressional and functional analysis of WRKY gene family reveals that PbeWRKY16 and PbeWRKY31 contribute to the Valsa canker resistance in Pyrus betulifolia

The WRKY transcription factor family plays a crucial role in regulating plant growth and stress responses. However, there are few studies on the regulation of resistance to Valsa canker. In this study, a comprehensive analysis of WRKY genes across 19 plant species was conducted. The potential members of Valsa canker resistance regulation were identified via functional validation. A total of 1641 WRKY genes could be categorized into seven groups. WRKY family members show subfamily- and species-specific expansions. In Rosaceae, Group II-d and II-e were rapidly expanded, which mainly originated based from whole genome duplication (WGD). Cis-element analysis and protein interaction network prediction underscored that most WRKYs respond to stress signals. Based on expressional investigation and Weighted Gene Co-expression Network Analysis (WGCNA), 9 WRKY genes in Pyrus betulaefolia were screened as candidates for Valsa canker resistance regulation. Functional analysis further demonstrated that PbeWRKY16 and PbeWRKY31 regulate the expression of genes involved in salicylic acid (SA) biosynthesis and transport, thereby enhancing resistance to Valsa canker and activating immune responses. Our results provide a foundation for understanding the evolutionary mechanisms of the WRKY gene family and screened potential family members on Valsa canker resistance regulation.

Unraveling the adaptive evolution and functional diversification of Rab GTPases in Salix matsudana under salt condition

Soil salinity, exacerbated by human activities and climate change, is a growing threat to agricultural and forestry productivity. Willow known for their extensive root systems and rapid growth plays a crucial role in soil stabilization and water regulation. This study investigates the Rab GTPase family in S. matsudana, which is vital for intracellular vesicular trafficking and plant responses to salt stress. We identified 120 SmRabs in S. matsudana, revealing a wide range of physicochemical properties. Phylogenetic analysis categorized these into eight clades (RabA-H), with RabA being the most populous. Duplication events were discovered, with 182 syntenic SmRab pairs, suggesting whole-genome or segmental duplications. Ka/Ks ratios indicated purifying selection, with a few exceptions suggesting positive selection. Quantitative real-time PCR analysis revealed differential expression of SmRabs under salt stresses, suggesting their roles in stress response. Heterologous expression of SmRabD2g in yeast enhanced salt tolerance, indicating a protective role. Arabidopsis transformation with SmRabD2g construct showed improved salt tolerance, with transgenic plants exhibiting reduced damage under salt stress. This comprehensive study provides insights into the evolutionary dynamics and functional roles of SmRabs in S. matsudana, offering potential for genetic engineering to enhance salt tolerance in willows, contributing to sustainable production in saline soils.

Studying sleep orthologs in Epsilonproteobacteria through an evolutionary lens: investigating sleep mysteries through phylogenomics

The current study employed phylogenomic methods to examine sleep-related genes' evolutionary role and significance in Sulfurimonas paralvinellae of the Epsilonproteobacteria class. This has facilitated the identification of conserved sleep orthologs, including DnaK (Hsp70), serine hydroxymethyltransferase (SHMT), and potassium channel family proteins, exhibiting sequence similarities ranging from 39.13% to 61.45%. These findings align with prior research indicating that chaperones and ion channels are conserved during sleep. This was demonstrated by the observation that proteins with fewer domains exhibited more significant conservation than others, such as adenylate kinase (AK). Distinct adaptations in bifunctional protein-serine/threonine kinases and phosphatases were linked to S. paralvinellae, an extremophilic organism adapted to high-pressure and/or high-temperature conditions, indicating functional divergence influenced by the organism's environment. The Gene Ontology study results indicated catalytic activity, potassium channel function, and cellular processes, underscoring the significance of ion channels in regulating the sleep-wake cycle. Furthermore, the categories not recognized as particularly significant for the over-represented genes encompassed metabolic and signal transduction categories, suggesting enhanced functional flexibility within this protein subfamily. The findings emphasize that orthologous interactions are complex and influenced by subfunctionalization and neofunctionalization of ecology and evolution. These findings enhance the existing understanding of the evolution of sleep-related genes and their association with metabolic and environmental changes, providing a foundation for subsequent experimental investigations and cross-taxonomic comparisons.

Comprehensive analysis of the NAC transcription factor gene family in Sophora tonkinensis Gagnep

BACKGROUND

Sophora tonkinensis Gagnep. has long been utilized in the treatment of anti-inflammatory and pain-relieving, with its principal active compounds being alkaloids and flavonoids. NAC transcription factors, a large family of plant-specific regulators, play pivotal roles in growth, development, stress responses, and secondary metabolism. However, comprehensive genome-wide characterization of S. tonkinensis NAC gene family (StNAC) remains unexplored.

RESULTS

This study identified 85 NAC proteins from the S. tonkinensis genome database. Phylogenetic analysis revealed that StNAC proteins were categorized into 15 subgroups based on their homology with Arabidopsis thaliana NAC proteins. Gene structure analysis demonstrated a variation in intron numbers ranging from 1 to 7, with a majority of StNAC genes containing 2-3 introns. Chromosomal distribution analysis indicated an uneven spread of StNAC genes across 9 chromosomes, with the highest number of StNAC genes on Chr3. Detection of 4 tandem duplicates and 32 segmental duplicates revealed that segmental duplication primarily drive StNAC genes amplification. Prediction of cis-regulatory elements suggested the involvement of StNAC genes in growth, stress responses, and hormone regulation. Gene expression analysis showed substantial variability expression of StNAC genes across different tissues. Notably, eight StNAC genes were identified as significantly associated alkaloid and flavonoid levels. qRT-PCR validation indicated that five genes were highly expressed in tissues, corroborating transcriptome data.

CONCLUSION

These findings offer valuable insights for further functional characterization of NAC genes and their potential roles in alkaloid and flavonoid biosynthesis in S. tonkinensis.

How enzyme functions evolve: genetic, structural, and kinetic perspectives

Understanding the emergence or loss of enzyme functions comprises several approaches, such as genetic, structural, and kinetic studies. Promiscuous enzyme activities have been proposed as starting points for the emergence of novel enzyme functions, for example, through genetic models such as neofunctionalization and subfunctionalization. In both cases, neutral evolution would fix gene redundancy, critical in relaxing functional constraints and allowing specific mutations to drive innovation. The evolution of enzyme activities has a structural basis, with genetic mutations modifying the active site architecture, conformational dynamics, or interaction networks, which leads to the creation, enhancement, or restriction of enzyme functions where epistatic interactions are crucial. These structural changes impact the described kinetic mechanisms like ground-state stabilization (affinity), transition-state stabilization (catalysis), or a combination of both. Case studies across diverse enzyme families illustrate these principles, emphasizing the interplay between genetic, structural, and kinetic approaches. Finally, we discuss the importance of understanding evolutionary mechanisms and their impact on protein engineering and drug design for biomedical and industrial applications. However, these studies highlight that further experimental evolutionary data collection is necessary to enable the training of advanced machine learning models for use in biotechnological applications.

Selection Pressure Regulates the Evolution-Structure-Function Paradigm of Monocyte Chemoattractant Protein Family

Monocyte chemoattractant proteins (MCPs) are involved in monocyte trafficking during severe inflammation by modulating the chemokine-glycosaminoglycan-receptor signaling axis. MCPs comprise a family of four chemokines (CCL2, CCL7, CCL8, and CCL13/12) that exhibit differential expression patterns in mammals, functional diversity, and receptor/glycosaminoglycan (GAG) binding promiscuity. In this context, the evolution-structure-function paradigm of MCP chemokines in mammals was established by assessing phylogeny, functional divergence, selection pressure, and coevolution in correlation with structural and surface characteristics. Comprehensive analyses were performed using an array of evolutionary and structural bioinformatic methods including molecular dynamics simulations. Our findings demonstrate that substitutions in receptor/GAG-interacting residues mediate episodic diversification and functional diversity in MCP chemokines. Additionally, a balanced interplay of selection pressures has driven the functional changes observed among MCP paralogs, with positive selection at various receptor/GAG-binding sites contributing to their promiscuous receptor/GAG interactions. Meanwhile, processes like purifying selection and coevolution maintain the classical chemokine structure and preserve the ancestral functions of MCP chemokines. Overall, this study suggests that selection pressure on sites within the N-terminal region [N-loop and 310-helix] and 40S loop of MCP chemokines alters surface properties to fine-tune the molecular interactions and functional characteristics without altering the overall chemokine structure.

Towards a Conceptual Framework to Better Understand the Advantages and Limitations of Model Organisms

Model organisms (MO) are widely used in neuroscience to study brain processes, behavior, and the biological foundation of human diseases. However, the use of MO has also been criticized for low reliability and insufficient success rate in the development of therapeutic approaches, because the success of MO use also led to overoptimistic and simplistic applications, which sometimes resulted in wrong conclusions. Here, we develop a conceptual framework of MO to support scientists in their practical work and to foster discussions about their power and limitations. For this purpose, we take advantage of concepts developed in the philosophy of science and adjust them for practical application by neuroscientists. We suggest that MO can be best understood as tools that are used to gain information about a group of species or a phenomenon in a species of interest. These learning processes are made possible by some properties of MO, which facilitate the process of acquisition of understanding or provide practical advantages, and the possibility to transfer information between species. However, residual uncertainty in the reliability of information transfer remains, and incorrect generalizations can be side-effects of epistemic benefits, which we consider as representational and epistemic risks. This suggests that to use MO most effectively, scientists should analyze the similarity relation between the involved species, weigh advantages and risks of certain epistemic benefits, and invest in carefully designed validation experiments. Altogether, our analysis illustrates how scientists can benefit from philosophical concepts for their research practice.

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.

Transcriptional atlas of Daphnia magna

Transcriptomics studies are more likely to achieve predictive results when they rely on tissue- and cell-specific transcriptional data. Identification of cell types in novel model organisms by their transcriptional profiles is difficult without data on transcriptional differences among major tissues and anatomical features. Here we report the first dataset on tissue- and organ-specific transcriptomics in freshwater plankton crustacean Daphnia magna, reporting markers of embryos, hemocytes, gut, carapace, antennae-2, and head, as well as the remaining carcass. Embryos are the most transcriptionally different from adults' features, with antennae and carapace being the most differentiated among them. We demonstrate that transcriptional markers of embryos vs. adults and of various adult anatomical features can be used to provide validation and functional explanation to published differential expression in response to environmental factors like infection, hypoxia, toxicants, or kairomones; to annotate Daphnia single cell data; and to ask questions about transcriptional diversification within extended gene families.

The chromosome-scale genomes of two Tinospora species reveal differential regulation of the MEP pathway in terpenoid biosynthesis

BACKGROUND

The relationship between gene family expansion and the resulting changes in plant phenotypes has shown remarkable complexity during the evolution. The gene family expansion has contributed to the diversity in plant phenotypes, specifically metabolites through neo-functionalization and sub-functionalization. However, the negative regulatory effects associated with the gene family expansion remain poorly understood.

RESULTS

Here, we present the chromosome-scale genomes of Tinospora crispa and Tinospora sinensis. Comparative genomic analyses demonstrated conserved chromosomal evolution within the Menispermaceae family. KEGG analysis revealed a significant enrichment of genes related to terpenoid biosynthesis in T. sinensis. However, T. crispa exhibited a higher abundance of terpenoids compared to T. sinensis. Detailed analysis revealed the expansion of genes encoding 1-hydroxy-2-methyl 2-(E)-butenyl 4-diphosphate synthase (HDS), a key enzyme in the 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway of terpenoid biosynthesis in T. sinensis. TsiHDS4 retained the ancestral function of converting methylerythritol cyclic diphosphate (MEcPP) to (E)-4-hydroxy-3-methylbut-2-enyl diphosphate (HMBPP). However, the noncanonical CDS-derived small peptide TsiHDS5 was shown to interact with TsiHDS4, inhibiting its catalytic activity. This interaction reduced the levels of HMBPP and isopentenyl pyrophosphate (IPP), which represent key substrates for downstream terpenoid biosynthesis.

CONCLUSIONS

These findings offer clues to decipher the variations in the MEP pathway of terpenoid biosynthesis between T. crispa and T. sinensis and form a basis for further detailed research on the negative regulation of expanded genes.

Proteostasis modulates gene dosage evolution in antibiotic-resistant bacteria

Evolution of gene expression frequently drives antibiotic resistance in bacteria. We had previously (Patel and Matange, eLife, 2021) shown that, in Escherichia coli, mutations at the mgrB locus were beneficial under trimethoprim exposure and led to overexpression of dihydrofolate reductase (DHFR), encoded by the folA gene. Here, we show that DHFR levels are further enhanced by spontaneous duplication of a genomic segment encompassing folA and spanning hundreds of kilobases. This duplication was rare in wild-type E. coli. However, its frequency was elevated in a lon-knockout strain, altering the mutational landscape early during trimethoprim adaptation. We then exploit this system to investigate the relationship between trimethoprim pressure and folA copy number. During long-term evolution, folA duplications were frequently reversed. Reversal was slower under antibiotic pressure, first requiring the acquisition of point mutations in DHFR or its promoter. Unexpectedly, despite resistance-conferring point mutations, some populations under high trimethoprim pressure maintained folA duplication to compensate for low abundance DHFR mutants. We find that evolution of gene dosage depends on expression demand, which is generated by antibiotic and exacerbated by proteolysis of drug-resistant mutants of DHFR. We propose a novel role for proteostasis as a determinant of copy number evolution in antibiotic-resistant bacteria.

The role of unbalanced segmental duplication in sex chromosome evolution in Australian ridge-tailed goannas

Varanids are known for conserved sex chromosomes, but there are differences in the size of the W chromosome but not in morphology among species representing varying stages of sex chromosome evolution. We tested for homology of the ZW sex chromosome system with size differences in varanids among four species from two lineages in Australia, the Odatria and the Gouldii. We found that while DNA sequences of the sex chromosomes are conserved in the species we tested, we also identified a homologous region on an enlarged autosomal microchromosome that shares sequences with the W chromosome in some isolated populations of V. acanthurus and V. citrinus from the Odatria lineage. The enlarged microchromosome was unpaired in all individuals tested and is likely an unbalanced segmental duplication translocated between chromosome 1, the W, and another microchromosome. This suggests an ancient balanced duplication homologous to the W and the terminal region of the long arm of chromosome 1. The most parsimonious explanation is that the duplicated region likely originated on chromosome 1. We hypothesised in our reconstruction that genes and related DNA sequences associated with the sex-linkage group have likely originated on an autosome. Subsequently, the sequences may have undergone duplication and translocation to the W chromosome, followed by the accumulation of lineage specific repeat elements and amplifications on the W at different rates in various lineages. Lastly, these sequences are likely to have undergone duplication and translocation to another autosomal microchromosome. Given the role of segmental duplications and translocations as important evolutionary drivers of speciation in other taxa, together with the rapid speciation that has occurred in Australian varanids, our findings provide broader insight into the evolutionary pathway leading to rapid chromosomal and genic divergence of species.

Chromosome-level genome assembly of Pontederia cordata L. provides insights into its rapid adaptation and variation of flower colours

Pontederia cordata L. is an aquatic ornamental plant native to the Americas but has been widely distributed in South Asia, Australia, and Europe. The genetic mechanisms behind its rapid adaptation and spread have not yet been well understood. To understand the mechanisms for its rapid adaptation, this study assembled the first chromosome-level genome of P. cordata. The genome assembly, which spans 527.5 Mb, is anchored on 8 pseudochromosomes with a scaffold N50 of 48 Mb and encompasses 29,389 protein-coding genes. Further analyses revealed that P. cordata had experienced 3 whole-genome duplications (WGDs) events. These WGDs are associated with gene family expansion and increased numbers of resistance gene analogs and transcription factors. Positive selection analysis indicated that genes derived from tandem duplication (TD) and proximal duplication were more likely to undergo positive selection, and were enriched in plant defense and disease resistance. These results implied that WGDs, TD, and positive selection enhanced the environmental adaptability of P. cordata. In addition, we found that down-regulation of F3'5'H, DFR, ANS, and UFGT likely caused the flower colour variation for P. cordata from violet to white. The first chromosome-level genome of P. cordata here provides a valuable genomic resource for investigating the rapid adaptation and flower colour variation of the species.

Multilevel Gene Expression Changes in Lineages Containing Adaptive Copy Number Variants

Copy number variants (CNVs) are an important class of genetic variation that can mediate rapid adaptive evolution. Whereas, CNVs can increase the relative fitness of the organism, they can also incur a cost due to the associated increased gene expression and repetitive DNA. We previously evolved populations of Saccharomyces cerevisiae over hundreds of generations in glutamine-limited (Gln-) chemostats and observed the recurrent evolution of CNVs at the GAP1 locus. To understand the role that gene expression plays in adaptation, both in relation to the adaptation of the organism to the selective condition and as a consequence of the CNV, we measured the transcriptome, translatome, and proteome of 4 strains of evolved yeast, each with a unique CNV, and their ancestor in Gln- chemostats. We find CNV-amplified genes correlate with higher mRNA abundance; however, this effect is reduced at the level of the proteome, consistent with post-transcriptional dosage compensation. By normalizing each level of gene expression by the abundance of the preceding step we were able to identify widespread differences in the efficiency of each level of gene expression. Genes with significantly different translational efficiency were enriched for potential regulatory mechanisms including either upstream open reading frames, RNA-binding sites for Ssd1, or both. Genes with lower protein expression efficiency were enriched for genes encoding proteins in protein complexes. Taken together, our study reveals widespread changes in gene expression at multiple regulatory levels in lineages containing adaptive CNVs highlighting the diverse ways in which genome evolution shapes gene expression.

Lineage-specific duplication and functional diversification of DOPA-decarboxylase genes in the Gryllidae family, as revealed in Gryllus bimaculatus

The DOPA-decarboxylase (DDC) gene is crucial for dopamine synthesis and influences various biological functions in insects, including body coloration, behavior, learning, and sleep. However, its evolutionary impact remains largely unexplored. This study reports on the tandem duplication of two bona fide ddc genes (ddc1 and ddc2) in the Gryllidae cricket family. We herein investigated the function of ddc1 and ddc2 using Gryllus bimaculatus (Gb) as a model. Our results revealed that Gb'ddc1 was expressed systemically, with its expression being higher immediately after molting compared to the stage following melanin pigmentation. In homozygous knockout mutants of Gb'ddc1, generated via CRISPR/Cas9, reduced body color pigmentation and had translucent cuticles, decreased dopamine levels, and over-accumulated DOPA. These mutants died shortly after hatching, likely due to cuticle defects, underscoring the essential role of dopamine, produced by Gb'ddc1, in melanin synthesis. Conversely, Gb'ddc2 expression was confined to the ovary and was not up-regulated after molting. Homozygous knockout mutants of Gb'ddc2 exhibited no body color defects, whereas hatchability and embryonic development rates were significantly reduced. Interestingly, dopamine levels in the ovaries were significantly elevated in Gb'ddc2 mutants. This suggests that normal ovarian dopamine levels, modulated by Gb'ddc2, are vital for fertility maintenance. The function of Gb'ddc2 differs from that of typical ddc, indicating neofunctionalization through evolutionary duplication. Overall, Gb'ddc1 and Gb'ddc2 have distinct functions, and precise regulation of ovarian dopamine levels using these two ddc genes may have enhanced cricket fertility.

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 de novo Gene Promotes Seed Germination Under Drought Stress in Arabidopsis

The origin of genes from noncoding sequences is a long-term and fundamental biological question. However, how de novo genes originate and integrate into the existing pathways to regulate phenotypic variations is largely unknown. Here, we selected 7 genes from 782 de novo genes for functional exploration based on transcriptional and translational evidence. Subsequently, we revealed that Sun Wu-Kong (SWK), a de novo gene that originated from a noncoding sequence in Arabidopsis thaliana, plays a role in seed germination under osmotic stress. SWK is primarily expressed in dry seed, imbibing seed and silique. SWK can be fully translated into an 8 kDa protein, which is mainly located in the nucleus. Intriguingly, SWK was integrated into an extant pathway of hydrogen peroxide content (folate synthesis pathway) via the upstream gene cytHPPK/DHPS, an Arabidopsis-specific gene that originated from the duplication of mitHPPK/DHPS, and downstream gene GSTF9, to improve seed germination in osmotic stress. In addition, we demonstrated that the presence of SWK may be associated with drought tolerance in natural populations of Arabidopsis. Overall, our study highlights how a de novo gene originated and integrated into the existing pathways to regulate stress adaptation.

Identification of the FBN gene family in tomato and functional analysis of SlFBN11 in the electron transport under low night temperature

FBNs are lipid-associated proteins that play a critical role in plant growth and stress response. However, the mechanisms of how FBNs proteins participate in the low night temperature response in tomato still unclear. Here we conducted a comprehensive genome-wide analysis of the FBN gene family in Solanum lycopersicum. In total, 14 SlFBN genes were identified, and information on their gene structures, protein motifs, phylogenetic relationships, and stress-related cis-regulatory elements (CREs) was provided. Among these, SlFBN11 was selected as a promising candidate for further functional characterization. The silencing of SlFBN11 destroys the redox balance of the PSI reaction center under low night temperature (LNT) stress, which led to increased ROS accumulation. Surprisingly, we found that the silencing of SlFNR2 also displayed an imbalance in electron transport of the PSI under LNT stress. Further experiments showed SlFBN11 can interact with SlFNR2 to positively response electron transport low night temperature. Collectively, the study provides a comprehensive analysis of the FBN genes family in Solanum lycopersicum and provides a theoretical basis for our understanding of the function of FBN genes in adaptation to LNT stresses.

Zebrafish: unraveling genetic complexity through duplicated genes

The zebrafish is an invaluable model organism for genetic, developmental, and disease research. Although its high conservation with humans is often cited as justification for its use, the zebrafish harbors oft-ignored genetic characteristics that may provide unique insights into gene structure and function. Zebrafish, along with other teleost fish, underwent an additional round of whole genome duplication after their split from tetrapods-resulting in an abundance of duplicated genes when compared to other vertebrates. These duplicated genes have evolved in distinct ways over the ensuing 350 million years. Thus, each gene within a duplicated gene pair has nuanced differences that create a unique identity. By investigating both members of the gene pair together, we can elucidate the mechanisms that underly protein structure and function and drive the complex interplay within biological systems, such as signal transduction cascades, genetic regulatory networks, and evolution of tissue and organ function. It is crucial to leverage such studies to explore these molecular dynamics, which could have far-reaching implications for both basic science and therapeutic development. Here, we will review the role of gene duplications and the existing models for gene divergence and retention following these events. We will also highlight examples within each of these models where studies comparing duplicated genes in the zebrafish have yielded key insights into protein structure, function, and regulation.

The cells are all-right: Regulation of the Lefty genes by separate enhancers in mouse embryonic stem cells

Enhancers play a critical role in regulating precise gene expression patterns essential for development and cellular identity; however, how gene-enhancer specificity is encoded within the genome is not clearly defined. To investigate how this specificity arises within topologically associated domains (TAD), we performed allele-specific genome editing of sequences surrounding the Lefty1 and Lefty2 paralogs in mouse embryonic stem cells. The Lefty genes arose from a tandem duplication event and these genes interact with each other in chromosome conformation capture assays which place these genes within the same TAD. Despite their physical proximity, we demonstrate that these genes are primarily regulated by separate enhancer elements. Through CRISPR-Cas9 mediated deletions to remove the intervening chromatin between the Lefty genes, we reveal a distance-dependent dosage effect of the Lefty2 enhancer on Lefty1 expression. These findings indicate a role for chromatin distance in insulating gene expression domains in the Lefty locus in the absence of architectural insulation.

Pleiotropy increases with gene age in six model multicellular eukaryotes

Fundamental traits of genes, including function, length and GC content, all vary with gene age. Pleiotropy, where a single gene affects multiple traits, arises through selection for novel traits and is expected to be removed from the genome through subfunctionalization following duplication events. It is unclear, however, how these opposing forces shape the prevalence of pleiotropy through time. We hypothesized that the prevalence of pleiotropy would be lowest in young genes, peak in middle aged genes, and then either decrease to a middling level in ancient genes or stay near the middle-aged peak, depending on the balance between exaptation and subfunctionalization. To address this question, we have calculated gene age and pleiotropic status for several model multicellular eukaryotes, including Homo sapiens, Mus musculus, Danio rerio, Drosophila melanogaster, Caenorhabditis elegans, and Arabidopsis thaliana. Gene age was determined by finding the most distantly related species that shared an ortholog using the Open Tree of Life and the Orthologous Matrix Database (OMAdb). Pleiotropic status was determined using both protein-protein interactions (STRINGdb) and associated biological processes (Gene Ontology). We found that middle-aged and ancient genes tend to be more pleiotropic than young genes, and that this relationship holds across all species evaluated and across both modalities of measuring pleiotropy. We also found absolute differences in the degree of pleiotropy based on gene functional class, but only when looking at biological process count. From these results we propose that there is a fundamental relationship between pleiotropy and gene age and further study of this relationship may shed light on the mechanism behind the functional changes genes undergo as they age.

Independent neofunctionalization of Dxo1 in Saccharomyces and Candida led to 25S rRNA processing function

Eukaryotic genomes typically encode one member of the DXO/Dxo1/Rai1 family of enzymes, which can hydrolyze the 5' ends of RNAs with a variety of structures that deviate from the canonical 7mGpppN. In contrast, the Saccharomyces genome encodes two family members and the second copy, Dxo1, is a distributive 5' exoribonuclease that is required for the final maturation of the 5' end of 25S rRNA from a 25S' precursor. Here we show that this 25S rRNA maturation function is not conserved across kingdoms, but arose in the budding yeasts. Interestingly, the origin of 25S processing capacity coincides with the duplication of this gene, and this capacity is absent in the nonduplicated genes. Strikingly, two different clades of budding yeasts have undergone parallel evolution: Both duplicated their DXO/Dxo1/Rai1 gene, and in both cases, one copy gained the 25S processing function. This was accompanied by many parallel sequence changes, a remarkable case of reproducible neofunctionalization.

Dosage constraint of the ribosome-associated molecular chaperone drives the evolution and fates of its duplicates in bacteria

Gene duplication events happen prevalently during evolution, and the mechanisms governing the loss or retention of duplicated genes are mostly elusive. Our genome scanning analysis revealed that trigger factor (TF), the one and only bacterial ribosome-associated molecular chaperone, is singly copied in virtually every bacterium except for a very few that possess two or more copies. However, even in these exceptions, only one complete TF copy exists, while other homologs lack the N-terminal domain that contains the conserved ribosome binding site (RBS) motif. Consistently, we demonstrated that the overproduction of the N-terminal complete TF proteins is detrimental to the cell, which can be rescued by removing the N-terminal domain. Our findings also indicated that TF overproduction leads to a decrease in protein productivity and profile changes in proteome due to its characteristic ribosome binding and holdase activities. Additionally, these N-terminal deficient TF homologs in bacteria with multiple TF homologs partition the function of TF via subfunctionalization. Our results revealed that TF is subjected to a dosage constraint that originates from its own intrinsic functions, which may drive the evolution and fates of duplicated TFs in bacteria.

IMPORTANCE

Gene duplication events presumably occur in tig, which encodes the ribosome-associated molecular chaperone trigger factor (TF). However, TF is singly copied in virtually every bacterium, and these exceptions with multiple TF homologs always retain only one complete copy while other homologs lack the N-terminal domain. Here, we reveal the manner and mechanism underlying the evolution and fates of TF duplicates in bacteria. We discovered that the mutation-to-loss or retention-to-sub/neofunctionalization of TF duplicates is associated with the dosage constraint of N-terminal complete TF. The dosage constraint of TF is attributed to its characteristic ribosome binding and substrate-holding activities, causing a decrease in protein productivity and profile changes in cellular proteome.

Substrate specificity and protein stability drive the divergence of plant-specific DNA methyltransferases

DNA methylation is an important epigenetic mechanism essential for transposon silencing and genome integrity. Across evolution, the substrates of DNA methylation have diversified between kingdoms. In plants, chromomethylase3 (CMT3) and CMT2 mediate CHG and CHH methylation, respectively. However, how these two methyltransferases diverge on substrate specificities during evolution remains unknown. Here, we reveal that CMT2 originates from a duplication of an evolutionarily ancient CMT3 in flowering plants. Lacking a key arginine residue recognizing CHG in CMT2 impairs its CHG methylation activity in most flowering plants. An engineered V1200R mutation empowers CMT2 to restore CHG and CHH methylations in Arabidopsis cmt2cmt3 mutant, testifying a loss-of-function effect for CMT2 during evolution. CMT2 has evolved a long and unstructured amino terminus critical for protein stability, especially under heat stress, and is plastic to tolerate various natural mutations. Together, this study reveals the mechanism of chromomethylase divergence for context-specific DNA methylation in plants and sheds important lights on DNA methylation evolution and function.

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.

Expression and functional divergence of a type-A response regulator paralog pair formed by dispersed duplication during Populus deltoides evolution

Gene duplication and divergence are essential to plant evolution. The Arabidopsis type-A response regulator (ARR) family, negative regulators in cytokinin signaling, exemplifies gene expansion and differential retention. Despite extensive research, the understanding of type-A RR homologs in woody plants remains limited. In this study, the evolution history of type-A RR gene families across four rosids and one monocot has been comprehensively investigated. Focusing on Populus deltoides, a unique pair of dispersed duplicates, PdRR8 and PdFERR, is identified, and their duplication is estimated to have occurred in the common ancestor of the four rosids. The duplication remnants corresponding to PdRR8 have been retained in all rosids but the counterpart of PdFERR has been lost. In poplar, PdRR8 shows the highest expression levels in leaves, while PdFERR is specifically expressed in female floral buds. Among various external stimuli, cold strongly represses PdRR8 promoter activity, whereas 6-BA markedly inhibits that of PdFERR. Overexpression of PdRR8 in the Arabidopsis arr16arr17 double-mutant fully complements the reduced hydrotropic response. In contrast, PdFERR fails to rescue the hydrotropic defects of the mutant. Results of evolutionary, expression and functional analyses indicate that PdRR8, rather than PdFERR, is the true ortholog of the ARR16-ARR17 paralogs. Though PdRR8 and PdFERR originate from a common ancestral gene and evolve under strong negative selection, these two dispersed duplicates have exhibited differential expression and some degree of functional divergence.

Quest for Orthologs in the Era of Biodiversity Genomics

The era of biodiversity genomics is characterized by large-scale genome sequencing efforts that aim to represent each living taxon with an assembled genome. Generating knowledge from this wealth of data has not kept up with this pace. We here discuss major challenges to integrating these novel genomes into a comprehensive functional and evolutionary network spanning the tree of life. In summary, the expanding datasets create a need for scalable gene annotation methods. To trace gene function across species, new methods must seek to increase the resolution of ortholog analyses, e.g. by extending analyses to the protein domain level and by accounting for alternative splicing. Additionally, the scope of orthology prediction should be pushed beyond well-investigated proteomes. This demands the development of specialized methods for the identification of orthologs to short proteins and noncoding RNAs and for the functional characterization of novel gene families. Furthermore, protein structures predicted by machine learning are now readily available, but this new information is yet to be integrated with orthology-based analyses. Finally, an increasing focus should be placed on making orthology assignments adhere to the findable, accessible, interoperable, and reusable (FAIR) principles. This fosters green bioinformatics by avoiding redundant computations and helps integrating diverse scientific communities sharing the need for comparative genetics and genomics information. It should also help with communicating orthology-related concepts in a format that is accessible to the public, to counteract existing misinformation about evolution.

The emerging role of fat-inducing transcript 2 in endoplasmic reticulum proteostasis and lipoprotein biogenesis

PURPOSE OF REVIEW

This review examines the evolving role of the fat-inducing transcript 2 (FIT2) protein in lipid droplet (LD) biology and its broader implications in cellular physiology and disease. With recent advancements in understanding FIT2 function across various model systems, this review provides a timely synthesis of its mechanisms and physiological significance.

RECENT FINDINGS

FIT2, an endoplasmic reticulum (ER)-resident protein, has been established as a critical regulator of LD formation in diverse organisms, from yeast to mammals. It facilitates LD biogenesis by sequestering diacylglycerol (DAG) and potentially influencing ER membrane dynamics. Beyond its role in lipid metabolism, FIT2 intersects with the ER-associated degradation (ERAD), is critical for protein homeostasis, and is linked to the unfolded protein response (UPR). Dysregulation of FIT2 has also been linked to metabolic disorders such as insulin resistance and lipodystrophy, highlighting its clinical relevance.

SUMMARY

Insights into FIT2 function underscore its pivotal role in LD formation and lipid homeostasis. Understanding its involvement in ER proteostasis and very low density lipoprotein biogenesis has broad implications for metabolic diseases and cancer. Therapeutic strategies targeting FIT2 may offer novel approaches to modulate lipid metabolism and mitigate associated pathologies. Further research is needed to elucidate the full spectrum of FIT2's interactions within cellular lipid and protein networks, potentially uncovering new therapeutic avenues for metabolic and ER stress-related disorders.

Functional investigation of the two ClpPs and three ClpXs in Myxococcus xanthus DK1622

ClpXP is a protease complex that plays important roles in protein quality control and cell cycle regulation, but the functions of multiple ClpXs and multiple ClpPs in M. xanthus remain unknown. The genome of Myxococcus xanthus DK1622 contains two clpPs and three clpXs. The clpP1 and clpX1 genes are cotranscribed and are both essential, while the other copies are isolated in the genome and are deletable. The deletion of clpX2 caused the mutant to be deficient in fruiting body development, while the clpX3 gene is involved in resistance to thermal stress. Both ClpPs possess catalytic active sites, but only ClpP1 shows in vitro peptidase activity on the typical substrate Suc-LY-AMC. All of these clpP and clpX genes exhibit strong transcriptional upregulation in the stationary phase, and the transcription of the three clpX genes appears to be coordinated. Our results demonstrated that multiple ClpPs and multiple ClpXs are functionally divergent and may assist in the environmental adaptation and functional diversification of M. xanthus.IMPORTANCEClpXP is an important protease complex of bacteria and is involved in various physiological processes. Myxococcus xanthus DK1622 possesses two ClpPs and three ClpXs with unclear functions. We investigated the functions of these genes and demonstrated the essential roles of clpP1 and clpX1. Only ClpP1 has in vitro peptidase activity on Suc-LY-AMC, and the isolated clpX copies participate in distinct cellular processes. All of these genes exhibited significant transcriptional upregulation in the stationary phase. Divergent functions appear in multiple ClpPs and multiple ClpXs in M. xanthus DK1622.

Genome-Wide Identification and Characterization of MYB Transcription Factors in Sudan Grass under Drought Stress

Sudan grass (Sorghum sudanense S.) is a warm-season annual grass with high yield, rich nutritional value, good regeneration, and tolerance to biotic and abiotic stresses. However, prolonged drought affects the yield and quality of Sudan grass. As one of the largest families of multifunctional transcription factors in plants, MYB is widely involved in regulating plant growth and development, hormonal signaling, and stress responses at the gene transcription level. However, the regulatory role of MYB genes has not been well characterized in Sudan grass under abiotic stress. In this study, 113 MYB genes were identified in the Sudan grass genome and categorized into three groups by phylogenetic analysis. The promoter regions of SsMYB genes contain different cis-regulatory elements, which are involved in developmental, hormonal, and stress responses, and may be closely related to their diverse regulatory functions. In addition, collinearity analysis showed that the expansion of the SsMYB gene family occurred mainly through segmental duplications. Under drought conditions, SsMYB genes showed diverse expression patterns, which varied at different time points. Interaction networks of 74 SsMYB genes were predicted based on motif binding sites, expression correlations, and protein interactions. Heterologous expression showed that SsMYB8, SsMYB15, and SsMYB64 all significantly enhanced the drought tolerance of yeast cells. Meanwhile, the subcellular localization of all three genes is in the nucleus. Overall, this study provides new insights into the evolution and function of MYB genes and provides valuable candidate genes for breeding efforts in Sudan grass.

Problems with Paralogs: The Promise and Challenges of Gene Duplicates in Evo-Devo Research

Gene duplicates, or paralogs, serve as a major source of new genetic material and comprise seeds for evolutionary innovation. While originally thought to be quickly lost or nonfunctionalized following duplication, now a vast number of paralogs are known to be retained in a functional state. Daughter paralogs can provide robustness through redundancy, specialize via sub-functionalization, or neo-functionalize to play new roles. Indeed, the duplication and divergence of developmental genes have played a monumental role in the evolution of animal forms (e.g., Hox genes). Still, despite their prevalence and evolutionary importance, the precise detection of gene duplicates in newly sequenced genomes remains technically challenging and often overlooked. This presents an especially pertinent problem for evolutionary developmental biology, where hypothesis testing requires accurate detection of changes in gene expression and function, often in nontraditional model species. Frequently, these analyses rely on molecular reagents designed within coding sequences that may be highly similar in recently duplicated paralogs, leading to cross-reactivity and spurious results. Thus, care is needed to avoid erroneously assigning diverged functions of paralogs to a single gene, and potentially misinterpreting evolutionary history. This perspective aims to overview the prevalence and importance of paralogs and to shed light on the difficulty of their detection and analysis while offering potential solutions.

A Rigorous Framework to Classify the Postduplication Fate of Paralogous Genes

Gene duplication has a central role in evolution; still, little is known on the fates of the duplicated copies, their relative frequency, and on how environmental conditions affect them. Moreover, the lack of rigorous definitions concerning the fate of duplicated genes hinders the development of a global vision of this process. In this paper, we present a new framework aiming at characterizing and formally differentiating the fate of duplicated genes. Our framework has been tested via simulations, where the evolution of populations has been simulated using aevol, an in silico experimental evolution platform. Our results show several patterns that confirm some of the conclusions from previous studies, while also exhibiting new tendencies; this may open up new avenues to better understand the role of duplications as a driver of evolution.

Functional redundancy of the three insulin receptors of cockroaches

Gene duplication is a fundamental evolutionary process which provides opportunities to acquire new gene functions. In the case of the insulin receptors (InRs) in cockroaches and close-related insects, two successive duplications determined the occurrence of three InR genes: InR2, InR1 and InR3, the last two forming a sister cluster to InR2. The biological role of each of the gene duplicates and whether they resulted from neofunctionalization or subfunctionalization is still unclear. The analysis of the sequences from different lineages did not detect positive selection as driving the divergence of InR1 and InR3, discarding neofunctionalization, and suggesting that there is no functional divergence between both gene copies. Using the cockroach Blattella germanica as a model, we have determined that BgInR2 is the gene with the highest expression levels in all the tissues analyzed, both in adult females and males, as well as in nymphs and embryos. BgInR3 is second in expression levels while BgInR1 is expressed at lower levels and only in some tissues. The selective depletion by RNAi of each of the three InRs, analyzed in terms of phenotype and fat body transcriptomic profiles, resulted in essentially redundant effects, with a magnitude approximately proportional to the level of expression of the respective InR. Therefore, the results indicate that the InR duplicates likely experienced a subfunctionalization process, by which the three InRs maintained similar functions but contributing to those functions proportionally to their expression levels.

Exploring Evolutionary Pathways and Abiotic Stress Responses through Genome-Wide Identification and Analysis of the Alternative Oxidase (AOX) Gene Family in Common Oat (Avena sativa)

The alternative oxidase (AOX), a common terminal oxidase in the electron transfer chain (ETC) of plants, plays a crucial role in stress resilience and plant growth and development. Oat (Avena sativa), an important crop with high nutritional value, has not been comprehensively studied regarding the AsAOX gene family. Therefore, this study explored the responses and potential functions of the AsAOX gene family to various abiotic stresses and their potential evolutionary pathways. Additionally, we conducted a genome-wide analysis to explore the evolutionary conservation and divergence of AOX gene families among three Avena species (Avena sativa, Avena insularis, Avena longiglumis) and four Poaceae species (Avena sativa, Oryza sativa, Triticum aestivum, and Brachypodium distachyon). We identified 12 AsAOX, 9 AiAOX, and 4 AlAOX gene family members. Phylogenetic, motif, domain, gene structure, and selective pressure analyses revealed that most AsAOXs, AiAOXs, and AlAOXs are evolutionarily conserved. We also identified 16 AsAOX segmental duplication pairs, suggesting that segmental duplication may have contributed to the expansion of the AsAOX gene family, potentially preserving these genes through subfunctionalization. Chromosome polyploidization, gene structural variations, and gene fragment recombination likely contributed to the evolution and expansion of the AsAOX gene family as well. Additionally, we hypothesize that AsAOX2 may have potential function in resisting wounding and heat stresses, while AsAOX4 could be specifically involved in mitigating wounding stress. AsAOX11 might contribute to resistance against chromium and waterlogging stresses. AsAOX8 may have potential fuction in mitigating ABA-mediated stress. AsAOX12 and AsAOX5 are most likely to have potential function in mitigating salt and drought stresses, respectively. This study elucidates the potential evolutionary pathways of the AsAOXs gene family, explores their responses and potential functions to various abiotic stresses, identifies potential candidate genes for future functional studies, and facilitates molecular breeding applications in A. sativa.

Altered assembly paths mitigate interference among paralogous complexes

Protein complexes are fundamental to all cellular processes, so understanding their evolutionary history and assembly processes is important. Gene duplication followed by divergence is considered a primary mechanism for diversifying protein complexes. Nonetheless, to what extent assembly of present-day paralogous complexes has been constrained by their long evolutionary pathways and how cross-complex interference is avoided remain unanswered questions. Subunits of protein complexes are often stabilized upon complex formation, whereas unincorporated subunits are degraded. How such cooperative stability influences protein complex assembly also remains unclear. Here, we demonstrate that subcomplexes determined by cooperative stabilization interactions serve as building blocks for protein complex assembly. We further develop a protein stability-guided method to compare the assembly processes of paralogous complexes in cellulo. Our findings support that oligomeric state and the structural organization of paralogous complexes can be maintained even if their assembly processes are rearranged. Our results indicate that divergent assembly processes by paralogous complexes not only enable the complexes to evolve new functions, but also reinforce their segregation by establishing incompatibility against deleterious hybrid assemblies.

Evolutionary and Expression Analysis of the Pig MAGE Gene Family

The melanoma-associated antigen (MAGE) family found in eukaryotes plays a crucial role in cell proliferation and differentiation, spermatogenesis, neural development, etc. This study explored the validation and evolution of MAGE genes in eukaryotic genomes and their distribution and expression patterns in pigs. In total, 249 MAGE genes were found on 13 eukaryotic species. In total, 33, 25, and 18 genes were located on human, mouse, and pig genomes, respectively. We found eight, four, and three tandemly duplicated gene clusters on the human, mouse, and pig genomes, respectively. The majority of MAGE genes in mammals are located on the X chromosome. According to the phylogenetic analysis, the MAGE family genes were classified into 11 subfamilies. The NDN gene in zebrafish (DreNDN) was the root of this evolutionary tree. In total, 10 and 11 MAGE genes on human and mouse genomes, respectively, exhibited a collinearity relationship with the MAGE genes on pig genomes. Taking the MAGE family genes in pigs, the MAGE subfamilies had similar gene structures, protein motifs, and biochemical attributes. Using the RNA-seq data of Duroc pigs and Rongchang pigs, we detected that the expression of type I MAGE genes was higher in reproductive tissues, but type II MAGE genes were predominantly expressed in the brain tissue. These findings are a valuable resource for gaining insight into the evolution and expression of the MAGE family genes.

Substrate specificity and protein stability drive the divergence of plant-specific DNA methyltransferases

DNA methylation is an important epigenetic mechanism essential for transposon silencing and genome integrity. Across evolution, the substrates of DNA methylation have diversified between kingdoms to account for genome complexity. In plants, Chromomethylase3 (CMT3) and CMT2 are the major methyltransferases mediating CHG and CHH methylation, respectively. However, how these two enzymes diverge on substrate specificities during evolution remains unknown. Here, we reveal that CMT2 originates from a duplication of the evolutionarily more ancient CMT3 in flowering plants. Lacking a key arginine residue recognizing CHG in CMT2 impairs its CHG methylation activity in most flowering plants. An engineered V1200R mutation empowers CMT2 to restore both CHG and CHH methylation in Arabidopsis cmt2cmt3 mutant, testifying a loss-of-function effect for CMT2 after ∼200 million years of evolution. Interestingly, CMT2 has evolved a long and unstructured N-terminus critical for balancing protein stability, especially under heat stress. Furthermore, CMT2 N-terminus is plastic and can be tolerant to various natural mutations. Together, this study reveals the mechanism of chromomethylase divergence for context-specific DNA methylation in plants and sheds important lights on DNA methylation evolution and function.

Identification of homologs of the Chlamydia trachomatis effector CteG reveals a family of Chlamydiaceae type III secreted proteins that can be delivered into host cells

Chlamydiae are a large group of obligate endosymbionts of eukaryotes that includes the Chlamydiaceae family, comprising several animal pathogens. Among Chlamydiaceae, Chlamydia trachomatis causes widespread ocular and urogenital infections in humans. Like many bacterial pathogens, all Chlamydiae manipulate host cells by injecting them with type III secretion effector proteins. We previously characterized the C. trachomatis effector CteG, which localizes at the host cell Golgi and plasma membrane during distinct phases of the chlamydial infectious cycle. Here, we show that CteG is a Chlamydiaceae-specific effector with over 60 homologs phylogenetically categorized into two distinct clades (CteG I and CteG II) and exhibiting several inparalogs and outparalogs. Notably, cteG I homologs are syntenic to C. trachomatis cteG, whereas cteG II homologs are syntenic among themselves but not with C. trachomatis cteG. This indicates a complex evolution of cteG homologs, which is unique among C. trachomatis effectors, marked by numerous events of gene duplication and loss. Despite relatively modest sequence conservation, nearly all tested CteG I and CteG II proteins were identified as type III secretion substrates using Yersinia as a heterologous bacterial host. Moreover, most of the type III secreted CteG I and CteG II homologs were delivered by C. trachomatis into host cells, where they localized at the Golgi region and cell periphery. Overall, this provided insights into the evolution of bacterial effectors and revealed a Chlamydiaceae family of type III secreted proteins that underwent substantial divergence during evolution while conserving the capacity to localize at specific host cell compartments.

Multilevel gene expression changes in lineages containing adaptive copy number variants

Copy-number variants (CNVs) are an important class of recurrent variants that mediate adaptive evolution. While CNVs can increase the relative fitness of the organism, they can also incur a cost. We previously evolved populations of Saccharomyces cerevisiae over hundreds of generations in glutamine-limited (Gln-) chemostats and observed the recurrent evolution of CNVs at the GAP1 locus. To understand the role that expression plays in adaptation, both in relation to the adaptation of the organism to the selective condition, and as a consequence of the CNV, we measured the transcriptome, translatome, and proteome of 4 strains of evolved yeast, each with a unique CNV, and their ancestor in Gln- conditions. We find CNV-amplified genes correlate with higher RNA abundance; however, this effect is reduced at the level of the proteome, consistent with post-transcriptional dosage compensation. By normalizing each level of expression by the abundance of the preceding step we were able to identify widespread divergence in the efficiency of each step in the gene in the efficiency of each step in gene expression. Genes with significantly different translational efficiency were enriched for potential regulatory mechanisms including either upstream open reading frames, RNA binding sites for SSD1, or both. Genes with lower protein expression efficiency were enriched for genes encoding proteins in protein complexes. Taken together, our study reveals widespread changes in gene expression at multiple regulatory levels in lineages containing adaptive CNVs highlighting the diverse ways in which adaptive evolution shapes gene expression.

Gene networks governing the response of a calcareous sponge to future ocean conditions reveal lineage-specific XBP1 regulation of the unfolded protein response

Marine sponges are predicted to be winners in the future ocean due to their exemplary adaptive capacity. However, while many sponge groups exhibit tolerance to a wide range of environmental insults, calcifying sponges may be more susceptible to thermo-acidic stress. To describe the gene regulatory networks that govern the stress response of the calcareous sponge, Leucetta chagosensis (class Calcarea, order Clathrinida), individuals were subjected to warming and acidification conditions based on the climate models for 2100. Transcriptome analysis and gene co-expression network reconstruction revealed that the unfolded protein response (UPR) was activated under thermo-acidic stress. Among the upregulated genes were two lineage-specific homologs of X-box binding protein 1 (XBP1), a transcription factor that activates the UPR. Alternative dimerization between these XBP1 gene products suggests a clathrinid-specific mechanism to reversibly sequester the transcription factor into an inactive form, enabling the rapid regulation of pathways linked to the UPR in clathrinid calcareous sponges. Our findings support the idea that transcription factor duplication events may refine evolutionarily conserved molecular pathways and contribute to ecological success.

Characteristics of duplicated gene expression and DNA methylation regulation in different tissues of allopolyploid Brassica napus

Plant polyploidization increases the complexity of epigenomes and transcriptional regulation, resulting in genome evolution and enhanced adaptability. However, few studies have been conducted on the relationship between gene expression and epigenetic modification in different plant tissues after allopolyploidization. In this study, we studied gene expression and DNA methylation modification patterns in four tissues (stems, leaves, flowers and siliques) of Brassica napusand its diploid progenitors. On this basis, the alternative splicing patterns and cis-trans regulation patterns of four tissues in B. napus and its diploid progenitors were also analyzed. It can be seen that the number of alternative splicing occurs in the B. napus is higher than that in the diploid progenitors, and the IR type increases the most during allopolyploidy. In addition, we studied the fate changes of duplicated genes after allopolyploidization in B. napus. We found that the fate of most duplicated genes is conserved, but the number of neofunctionalization and specialization is also large. The genetic fate of B. napus was classified according to five replication types (WGD, PD, DSD, TD, TRD). This study also analyzed generational transmission analysis of expression and DNA methylation patterns. Our study provides a reference for the fate differentiation of duplicated genes during allopolyploidization.

Mutational robustness and the role of buffer genes in evolvability

Organisms rely on mutations to fuel adaptive evolution. However, many mutations impose a negative effect on fitness. Cells may have therefore evolved mechanisms that affect the phenotypic effects of mutations, thus conferring mutational robustness. Specifically, so-called buffer genes are hypothesized to interact directly or indirectly with genetic variation and reduce its effect on fitness. Environmental or genetic perturbations can change the interaction between buffer genes and genetic variation, thereby unmasking the genetic variation's phenotypic effects and thus providing a source of variation for natural selection to act on. This review provides an overview of our understanding of mutational robustness and buffer genes, with the chaperone gene HSP90 as a key example. It discusses whether buffer genes merely affect standing variation or also interact with de novo mutations, how mutational robustness could influence evolution, and whether mutational robustness might be an evolved trait or rather a mere side-effect of complex genetic interactions.

Decoding drought resilience: a comprehensive exploration of the cotton Eceriferum (CER) gene family and its role in stress adaptation

BACKGROUND

The cuticular wax serves as a primary barrier that protects plants from environmental stresses. The Eceriferum (CER) gene family is associated with wax production and stress resistance.

RESULTS

In a genome-wide identification study, a total of 52 members of the CER family were discovered in four Gossypium species: G. arboreum, G. barbadense, G. raimondii, and G. hirsutum. There were variations in the physicochemical characteristics of the Gossypium CER (GCER) proteins. Evolutionary analysis classified the identified GCERs into five groups, with purifying selection emerging as the primary evolutionary force. Gene structure analysis revealed that the number of conserved motifs ranged from 1 to 15, and the number of exons varied from 3 to 13. Closely related GCERs exhibited similar conserved motifs and gene structures. Analyses of chromosomal positions, selection pressure, and collinearity revealed numerous fragment duplications in the GCER genes. Additionally, nine putative ghr-miRNAs targeting seven G. hirsutum CER (GhCER) genes were identified. Among them, three miRNAs, including ghr-miR394, ghr-miR414d, and ghr-miR414f, targeted GhCER09A, representing the most targeted gene. The prediction of transcription factors (TFs) and the visualization of the regulatory TF network revealed interactions with GhCER genes involving ERF, MYB, Dof, bHLH, and bZIP. Analysis of cis-regulatory elements suggests potential associations between the CER gene family of cotton and responses to abiotic stress, light, and other biological processes. Enrichment analysis demonstrated a robust correlation between GhCER genes and pathways associated with cutin biosynthesis, fatty acid biosynthesis, wax production, and stress response. Localization analysis showed that most GCER proteins are localized in the plasma membrane. Transcriptome and quantitative reverse transcription-polymerase chain reaction (qRT-PCR) expression assessments demonstrated that several GhCER genes, including GhCER15D, GhCER04A, GhCER06A, and GhCER12D, exhibited elevated expression levels in response to water deficiency stress compared to control conditions. The functional identification through virus-induced gene silencing (VIGS) highlighted the pivotal role of the GhCER04A gene in enhancing drought resistance by promoting increased tissue water retention.

CONCLUSIONS

This investigation not only provides valuable evidence but also offers novel insights that contribute to a deeper understanding of the roles of GhCER genes in cotton, their role in adaptation to drought and other abiotic stress and their potential applications for cotton improvement.

ATAD3 Proteins: Unique Mitochondrial Proteins Essential for Life in Diverse Eukaryotic Lineages

ATPase family AAA domain-containing 3 (ATAD3) proteins are unique mitochondrial proteins that arose deep in the eukaryotic lineage but that are surprisingly absent in Fungi and Amoebozoa. These ∼600-amino acid proteins are anchored in the inner mitochondrial membrane and are essential in metazoans and Arabidopsis thaliana. ATAD3s comprise a C-terminal ATPases Associated with a variety of cellular Activities (AAA+) matrix domain and an ATAD3_N domain, which is located primarily in the inner membrane space but potentially extends to the cytosol to interact with the ER. Sequence and structural alignments indicate that ATAD3 proteins are most similar to classic chaperone unfoldases in the AAA+ family, suggesting that they operate in mitochondrial protein quality control. A. thaliana has four ATAD3 genes in two distinct clades that appear first in the seed plants, and both clades are essential for viability. The four genes are generally coordinately expressed, and transcripts are highest in growing apices and imbibed seeds. Plants with disrupted ATAD3 have reduced growth, aberrant mitochondrial morphology, diffuse nucleoids and reduced oxidative phosphorylation complex I. These and other pleiotropic phenotypes are also observed in ATAD3 mutants in metazoans. Here, we discuss the distribution of ATAD3 proteins as they have evolved in the plant kingdom, their unique structure, what we know about their function in plants and the challenges in determining their essential roles in mitochondria.

Mutational biases favor complexity increases in protein interaction networks after gene duplication

Biological systems can gain complexity over time. While some of these transitions are likely driven by natural selection, the extent to which they occur without providing an adaptive benefit is unknown. At the molecular level, one example is heteromeric complexes replacing homomeric ones following gene duplication. Here, we build a biophysical model and simulate the evolution of homodimers and heterodimers following gene duplication using distributions of mutational effects inferred from available protein structures. We keep the specific activity of each dimer identical, so their concentrations drift neutrally without new functions. We show that for more than 60% of tested dimer structures, the relative concentration of the heteromer increases over time due to mutational biases that favor the heterodimer. However, allowing mutational effects on synthesis rates and differences in the specific activity of homo- and heterodimers can limit or reverse the observed bias toward heterodimers. Our results show that the accumulation of more complex protein quaternary structures is likely under neutral evolution, and that natural selection would be needed to reverse this tendency.

Genome-wide identification of Aux/IAA gene family in white clover (Trifolium repens L.) and functional verification of TrIAA18 under different abiotic stress

BACKGROUND

White clover (Trifolium repens L.) is an excellent leguminous cool-season forage with a high protein content and strong nitrogen-fixing ability. Despite these advantages, its growth and development are markedly sensitive to environmental factors. Indole-3-acetic acid (IAA) is the major growth hormone in plants, regulating plant growth, development, and response to adversity. Nevertheless, the specific regulatory functions of Aux/IAA genes in response to abiotic stresses in white clover remain largely unexplored.

RESULTS

In this study, we identified 47 Aux/IAA genes in the white clover genome, which were categorized into five groups based on phylogenetic analysis. The TrIAAs promoter region co-existed with different cis-regulatory elements involved in developmental and hormonal regulation, and stress responses, which may be closely related to their diverse regulatory roles. Collinearity analysis showed that the amplification of the TrIAA gene family was mainly carried out by segmental duplication. White clover Aux/IAA genes showed different expression patterns in different tissues and under different stress treatments. In addition, we performed a yeast two-hybrid analysis to investigate the interaction between white clover Aux/IAA and ARF proteins. Heterologous expression indicated that TrIAA18 could enhance stress tolerance in both yeast and transgenic Arabidopsis thaliana.

CONCLUSION

These findings provide new scientific insights into the molecular mechanisms of growth hormone signaling in white clover and its functional characteristics in response to environmental stress.

Comparative genomics of sirenians reveals evolution of filaggrin and caspase-14 upon adaptation of the epidermis to aquatic life

The mammalian epidermis has evolved to protect the body in a dry environment. Genes of the epidermal differentiation complex (EDC), such as FLG (filaggrin), are implicated in the barrier function of the epidermis. Here, we investigated the molecular evolution of the EDC in sirenians (manatees and dugong), which have adapted to fully aquatic life, in comparison to the EDC of terrestrial mammals and aquatic mammals of the clade Cetacea (whales and dolphins). We show that the main subtypes of EDC genes are conserved or even duplicated, like late cornified envelope (LCE) genes of the dugong, whereas specific EDC genes have undergone inactivating mutations in sirenians. FLG contains premature stop codons in the dugong, and the ortholog of human CASP14 (caspase-14), which proteolytically processes filaggrin, is pseudogenized in the same species. As FLG and CASP14 have also been lost in whales, these mutations represent convergent evolution of skin barrier genes in different lineages of aquatic mammals. In contrast to the dugong, the manatee has retained functional FLG and CASP14 genes. FLG2 (filaggrin 2) is truncated in both species of sirenians investigated. We conclude that the land-to-water transition of sirenians was associated with modifications of the epidermal barrier at the molecular level.

A chromosome-level genome reveals genome evolution and molecular basis of anthraquinone biosynthesis in Rheum palmatum

BACKGROUND

Rhubarb is one of common traditional Chinese medicine with a diverse array of therapeutic efficacies. Despite its widespread use, molecular research into rhubarb remains limited, constraining our comprehension of the geoherbalism.

RESULTS

We assembled the genome of Rheum palmatum L., one of the source plants of rhubarb, to elucidate its genome evolution and unpack the biosynthetic pathways of its bioactive compounds using a combination of PacBio HiFi, Oxford Nanopore, Illumina, and Hi-C scaffolding approaches. Around 2.8 Gb genome was obtained after assembly with more than 99.9% sequences anchored to 11 pseudochromosomes (scaffold N50 = 259.19 Mb). Transposable elements (TE) with a continuous expansion of long terminal repeat retrotransposons (LTRs) is predominant in genome size, contributing to the genome expansion of R. palmatum. Totally 30,480 genes were predicted to be protein-coding genes with 473 significantly expanded gene families enriched in diverse pathways associated with high-altitude adaptation for this species. Two successive rounds of whole genome duplication event (WGD) shared by Fagopyrum tataricum and R. palmatum were confirmed. We also identified 54 genes involved in anthraquinone biosynthesis and other 97 genes entangled in flavonoid biosynthesis. Notably, RpALS emerged as a compelling candidate gene for the octaketide biosynthesis after the key residual screening.

CONCLUSION

Overall, our findings offer not only an enhanced understanding of this remarkable medicinal plant but also pave the way for future innovations in its genetic breeding, molecular design, and functional genomic studies.