Evolution by gene duplication: an update

Two homoeoallelic gene expression of TaCHLIs ensures normal chlorophyll biosynthesis in Hexaploid wheat

Being polyploid has a fitness advantage but is physically complex. During polyploid plant evolution, some duplicate genes retain their ancestral function, which affected the plant phenotype in allelic dosage or functional redundancy. However, how duplicated genes whose products needed to form functional complexes coped with deleterious mutations remained unclear. Here, we report a yellow green leaf-2 (ygl2) mutant with yellow-green leaves derived from a cross of Shaan 3025 (S3025) and Shi 4185 (S4185) that was controlled by a combination of Tachli-7A null and Tachli-7B truncation, whereas Tachli-7A null or Tachli-7B truncation individually resulted in normal leaf colour. Our results indicated genetic complementarity between TaCHLI-7A and TaCHLI-7B is responsible for normal chloroplast development. Furthermore, TaCHLI-7D was conserved in ygl2, S3025 and S4185 at both sequence and expression levels. Furthermore, two-thirds of the total mRNA abundance in S4185 with Tachli-7A null was sufficient for chlorophyll synthesis, indicating that redundant mRNA dosage was the reason for genetic complementarity. Particularly, Tachli-7A null can be retained in several modern cultivars with no disadvantage under field conditions, probably because the redundant mRNA dosage is expected to buffer the gene imbalance caused by the imperfect relationship between different copies of TaCHLIs and their molecular interactors. Furthermore, the loss of TaCHLI-7A seems to preserve the minimum dosage and maximise simplification. Our findings provide evidence of homoeologs loss and functional mechanism during polyploid evolution.

A comprehensive analysis of the O-methyltransferase gene family in the chili pepper (Capsicum annuum) identifies COMT36 involved in capsaicinoids biosynthesis

O-methyltransferases (OMTs) play vital roles in secondary metabolites methylation. While chili peppers (Capsicum annuum) are rich in specialized metabolites, OMT functional characterization remains limited in this species. In this study, we identified 58 OMT genes through genome-wide analysis, revealing uneven chromosomal distribution with frequent tandem duplications contributing to gene family expansion. Phylogenetic and conserved motif analysis classified these OMTs into caffeoyl-CoA O-methyltransferases (CCoAOMT) and caffeic acid O-methyltransferases (COMT) subfamilies. Differential expression patterns across tissues and promoter cis-element variations suggested functional diversification. Notably, CaCOMT36 was demonstrated to participate in capsaicinoid biosynthesis. These findings provide valuable insights into the molecular characteristics and biological functions of the OMT gene family in chili pepper, particularly the potential mechanisms associated with O-methylation during capsaicinoid biosynthesis.

P-Type Pentatricopeptide Repeat Proteins YS1 and YS2 Function in Splicing of petB Intron to Maintain Chloroplast Homeostasis During Rice Seedling Development

Regulating chloroplast gene expression is crucial for maintaining chloroplast function and plant development. Pentatricopeptide repeat (PPR) proteins form a vast protein family that regulates organelle genes and has multiple functions during plant development. Here, we found that two P-type PPR proteins, YS1 (yellow-green seedling 1) and YS2, jointly regulated seedling development in rice. The loss of YS1 and YS2 exhibited the collapsed chloroplast thylakoids and decreased photosynthetic activity, leading to the yellowing and death of rice seedlings. YS1 and YS2 could directly bind to the transcript of the psbH-petB intergenic region to facilitate the splicing of petB intron, thereby affecting the splicing efficiency of petD, which is located downstream of petB in the five-cistronic transcription unit psbB-psbT-psbH-petB-petD. The mutations in YS1 and YS2 led to decreased mature transcripts of petB and petD after splicing, significantly reducing the protein levels of PetB and PetD. This further led to deficiencies in the cytochrome b6/f and photosystem I complexes of the electron transport chain (ETC), ultimately resulting in decreased ETC-produced NADPH and reduced contents of carbohydrates in ys mutants. Moreover, transcriptome sequencing analysis revealed that YS1 and YS2 were vital for chloroplast organization and carbohydrate metabolism, as well as chloroplast RNA processing. In previous studies, the mechanism of petB intron splicing in the five-cistronic transcription unit psbB-psbT-psbH-petB-petD of rice is unclear. Our study revealed that the two highly conserved proteins YS1 and YS2 were functionally redundant and played critical roles in photosynthesis and seedling development through their involvement in petB intron splicing to maintain chloroplast homeostasis in rice. This work broadened the perspective on PPR-mediated chloroplast development and laid a foundation for exploring the biofunctions of duplicated genes in higher plants.

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.

Exploring species taxonomic kingdom using information entropy and nucleotide compositional features of coding sequences based on machine learning methods

The flow of genetic information from DNA to protein is governed by the central dogma of molecular biology. Genetic drift and mutations usually lead to changes in DNA composition, thereby affecting the coding sequences (CDS) that encode functional proteins. Analyzing the nucleotide distribution in the coding regions of species is crucial for understanding their evolution. In this study, we applied Markov processes to analyze codon formation in 37,031,061 CDSs across 3,735 species genomes, spanning viruses, archaea, bacteria, and eukaryotes, to explore compositional changes. Our results revealed species preferences for different nucleotides. Information entropies and Markov information densities show that eukaryotes exhibit higher redundancy, followed by viruses, suggesting more gene duplication in eukaryotes and high mutation rates in viruses. Evolutionary trends showed an increase in information entropy and a decrease in Markov entropy, with negative correlations between first- and second-order Markov information densities. Furthermore, uniform manifold approximation and projection (UMAP) was used to reduce information redundancy for revealing unique evolutionary patterns in species classification. The machine learning methods demonstrated excellent performance in species classification accuracy, providing profound insights into CDS evolution and protein synthesis.

Phylogenetic Insights into the Evolutionary History of the RSPO Gene Family in Metazoa

Background: The RSPO gene family encodes secreted glycoproteins that are rich in cysteine, which generally serve as activators of the Wnt signaling pathway in animals. Four types of this family have been identified in a few model species. However, the evolution of the family remains unclear. Methods: In this study, we identified a total of 1496 RSPO homologs through an extensive survey of the RSPO genes in 430 animals. Gene family clustering and phylogenetic analysis identified four major subtypes of the family (RSPO1-RSPO4) and clarified their distribution of copy number in different species. Results and Conclusions: Members of the RSPO4 subfamily that were closest to ancestral forms existed in both Deuterostomes and Protostomates, and we speculate that representatives of this subfamily already existed in Urbilatera, the last common ancestor of Deuterostomes. Particularly, in some RSPO3 subtypes of Actinopterygii (ray-finned fishes), an FU repeated motif with three conserved cysteines was identified. Further conservative analysis of amino acids and alignment of tertiary protein structure revealed the potential functional sites for each subgroup. The results provide insight into the phylogenetic relationships and evolutionary patterns of conserved motifs of RSPO family genes in animal kingdoms, which will guide further studies on the biological functions of RSPO in other non-model species.

Interred mechanisms of resistance and host immune evasion revealed through network-connectivity analysis of M. tuberculosis complex graph pangenome

Mycobacterium tuberculosis complex successfully adapts to environmental pressures through mechanisms of rapid adaptation which remain poorly understood despite knowledge gained through decades of research. In this study, we used 110 reference-quality, complete de novo assembled, long-read sequenced clinical genomes to study patterns of structural adaptation through a graph-based pangenome analysis, elucidating rarely studied mechanisms that enable enhanced clinical phenotypes offering a novel perspective to the species' adaptation. Across isolates, we identified a pangenome of 4,325 genes (3,767 core and 558 accessory), revealing 290 novel genes, and a substantially more complete account of difficult-to-sequence esx/pe/pgrs/ppe genes. Seventy-four percent of core genes were deemed non-essential in vitro, 38% of which support the pathogen's survival in vivo, suggesting a need to broaden current perspectives on essentiality. Through information-theoretic analysis, we reveal the ppe genes that contribute most to the species' diversity-several with known consequences for antigenic variation and immune evasion. Construction of a graph pangenome revealed topological variations that implicate genes known to modulate host immunity (Rv0071-73, Rv2817c, cas2), defense against phages/viruses (cas2, csm6, and Rv2817c-2821c), and others associated with host tissue colonization. Here, the prominent trehalose transport pathway stands out for its involvement in caseous granuloma catabolism and the development of post-primary disease. We show paralogous duplications of genes implicated in bedaquiline (mmpL5 in all L1 isolates) and ethambutol (embC-A) resistance, with a paralogous duplication of its regulator (embR) in 96 isolates. We provide hypotheses for novel mechanisms of immune evasion and antibiotic resistance through gene dosing that can escape detection by molecular diagnostics.IMPORTANCEM. tuberculosis complex (MTBC) has killed over a billion people in the past 200 years alone and continues to kill nearly 1.5 million annually. The pathogen has a versatile ability to diversify under immune and drug pressure and survive, even becoming antibiotic persistent or resistant in the face of harsh chemotherapy. For proper diagnosis and design of an appropriate treatment regimen, a full understanding of this diversification and its clinical consequences is desperately needed. A mechanism of diversification that is rarely studied systematically is MTBC's ability to structurally change its genome. In this article, we have de novo assembled 110 clinical genomes (the largest de novo assembled set to date) and performed a pangenomic analysis. Our pangenome provides structural variation-based hypotheses for novel mechanisms of immune evasion and antibiotic resistance through gene dosing that can compromise molecular diagnostics and lead to further emergence of antibiotic resistance.

Genome-wide identification, phylogenetic investigation and abiotic stress responses analysis of the PP2C gene family in litchi (Litchi chinensis Sonn.)

As an important regulatory protein phosphatase in the abscisic acid (ABA) signal transduction pathway and mitogen-activated protein kinases (MAPK) cascade, type-2C protein phosphatase (PP2C) plays crucial roles in plant responses to abiotic stresses. However, the PP2C gene family's responses to abiotic stress in litchi (Litchi chinensis Sonn.) have not been systematically studied. In this study, we predicted the 68 PP2C (designated LcPP2C) genes randomly distributed across fourteen chromosomes in the litchi genome. Phylogenetic tree analysis among litchi, Arabidopsis (Arabidopsis thaliana), and rice (Oryza sativa) revealed that the phylogenetic tree was divided into thirteen groups (A, B, C, D, E, F1, F2, G, H, I, J, K, and L). Closely linked LcPP2C genes within the same group exhibited various similarities in gene structures and motif compositions. Collinearity analysis demonstrated that segmental duplication (SD) events were the main dramatically increasing numbers in the LcPP2C gene family members. Cis-acting element analysis revealed that the 68 LcPP2C genes contained hormone and stress response elements with varying quantities, implying their potential in litchi stress resistance. Expression analysis showed that all the LcPP2C genes exhibited varying expression levels across nine different litchi tissues, more than 50% of genes within each group displayed similar tissue-specific expression patterns. The expression intensity, duration and regulation direction (up- or down-regulation) of the LcPP2C genes were varied under different abiotic stresses (cold, heat, and drought). The physiological and biochemical tests indicated that eight activation indexes (peroxidase (POD), catalase (CAT), superoxide dismutase (SOD), malondialdehyde (MDA), proline (PRO), soluble protein (SP), hydrogen peroxide (H2O2), and soluble sugar (SS)) increase at different level. Additionally, we analyzed physicochemical properties, subcellular locations, and secondary structures of the LcPP2C family members. Notably, the extensive connectivity of LcPP2C32/60/9/37 underscored their vital roles in orchestrating and regulating biomolecular networks. These results provide valuable information for the identification of the LcPP2C genes and ideas for the cultivation of its transgenic induction lines in litchi.

Strategic targeting of Cas9 nickase expands tandem gene arrays

Expanding tandem gene arrays facilitates adaptation through dosage effects and gene family formation via sequence diversification. However, experimental induction of such expansions remains challenging. Here, we introduce a method termed break-induced replication (BIR)-mediated tandem repeat expansion (BITREx) to address this challenge. BITREx places Cas9 nickase adjacent to a tandem gene array to break the replication fork that has just replicated the array, forming a single-ended double-strand break. This break is subsequently end-resected to become single stranded. Since there is no repeat unit downstream of the break, the single-stranded DNA often invades an upstream unit to initiate ectopic BIR, resulting in array expansion. BITREx has successfully expanded gene arrays in budding yeast, with the CUP1 array reaching ∼1 Mb. Furthermore, appropriate splint DNAs allow BITREx to generate tandem arrays de novo from single-copy genes. We have also demonstrated BITREx in mammalian cells. Therefore, BITREx will find various unique applications in genome engineering.

Genome-wide identification and expression profiling of the hemoglobin gene family under hypoxia stress in Sillago sihama

Hemoglobin (Hb) plays a pivotal role in oxygen transport and is essential for the adaptive response to hypoxic stress. Studying the molecular evolution characteristics and expression pattern of S. sihama Hb gene family can provide theoretical basis for further revealing the roles of the Hb genes against hypoxia stress. This study aimed to identify the molecular evolution characteristics and expression patterns of Hb genes in S. sihama. Utilizing bioinformatics methods, we identified ten members of the Hb gene family within the whole genome of S. sihama, located on chromosomes LG16 and LG21 in a tandem repeat pattern. Phylogenetic analysis classified the Hb gene family into two subfamilies (alpha and beta), showing close relationships with the Hb gene family members of the large yellow croaker (Larimichthys crocea). Both phylogenetic and collinearity analysis indicated the conservation of Hb genes throughout evolution. The nonsynonymous substitution rate / synonymous substitution rate (Ka/Ks) results showed that its evolution was purifying selection. RNA-seq and qRT-PCR analyses revealed that five Hb genes were highlighted high expression levels in the heart and liver of S. sihama, with sshba1.4, sshbad, and sshba1 showing significant up-regulation in response to hypoxic stress. In summary, this study highlights the vital role of the Hb gene family in the hypoxic stress response of S. sihama, suggesting their potential as crucial candidate genes for regulating environmental stress in this species.

Genome-wide identification and characterization of ALKB homolog gene family in wheat (Triticum aestivum L.)

Introduction

N6-methyladenosine (m6A) is the most prevalent posttranscriptional modification in eukaryotic mRNAs. AlkB homologs (ALKBHs) are involved in plant responses to stress by modulating m6A methylation. However, homologous genes in wheat remain largely uncharacterized.

Methods and results

In this study, 30 ALKBH genes were identified in wheat, and analyzed their physicochemical properties. The phylogenetic analysis allowed the classification of these genes into seven distinct subfamilies. Additionally, their conserved domains, motif compositions, gene structures, chromosomal localization, and synteny, and the predicted cis-acting elements within their promoters were examined. Expression analysis revealed that TaALKBH9B-5 exhibited the highest expression and its demethylase activity was investigated. Furthermore, TaALKBH9B-5 was significantly upregulated in response to abscisic acid treatment and cold stress, indicating a positive regulatory trend.

Discussion

In conclusion, this study provides a comprehensive genomic assessment of the TaALKBH gene family and offers a theoretical framework for understanding the role of TaALKBH9B in the response to abiotic stress in wheat.

Characterization of avian β-defensin genes in Galliformes reveals widespread evolutionary diversification and distinct evolutionary relationships with infection risk

BACKGROUND

Avian β-defensins (AvBDs) represent a key family of antimicrobial host defense peptides in birds. Accumulating evidence suggests that the evolutionary trajectory of β-defensin genes is specific to the gene, timescale, and species involved, implying that species-specific ecological and life-history differences drive divergent selective pressures on these genes. However, their evolutionary dynamics, particularly the interactions with ecological factors and life-history traits, remain insufficiently explored.

RESULTS

Through a comprehensive survey of 25 species spanning all major clades of Galliformes, 354 AvBD genes were identified. Comparative sequence analysis, genomic organization, and phylogenetic studies collectively reveal significant evolutionary diversification characterized by gene duplication, pseudogenization, and gene loss across these species. Notably, chicken AvBD3 exhibits significant differences in its coding regions, while AvBD6 and AvBD7 appear to have copy number variations, with species-specific paralogs of AvBD6 being especially prominent. Moreover, positive selection was more frequently observed in recently diverged gene lineages compared to ancestral ones. Using 70 samples from eight galliform species, the study further identified the prevalence of species-specific amino acid alleles. Phylogenetic comparative analysis demonstrated that the evolution of nine AvBD genes (AvBD2, -4, -5, -8, -9, -10, -11, -12, and -14) is significantly associated with specific ecological factors and life-history characteristics. Additionally, the evolutionary rates of these genes showed distinct relationship with inferred infection risk, likely reflecting the multifunctionality of β-defensins and potential trade-offs between immune defense and other biological functions.

CONCLUSIONS

This cross-species identification and systematic evolutionary analysis of AvBDs in Galliformes deepen our understanding of the co-evolution of host defense peptides, offering valuable insights into their natural biology and evolution, and paving the way for future applications as alternatives to traditional antibiotics.

Evolutionary and functional divergence of Sfx, a plasmid-encoded H-NS homolog, underlies the regulation of IncX plasmid conjugation

Conjugative plasmids are widespread among prokaryotes, highlighting their evolutionary success. Conjugation systems on most natural plasmids are repressed by default. The negative regulation of F-plasmid conjugation is partially mediated by the chromosomal nucleoid-structuring protein (H-NS). Recent bioinformatic analyses have revealed that plasmid-encoded H-NS homologs are widespread and exhibit high sequence diversity. However, the functional roles of most of these homologs and the selective forces driving their phylogenetic diversification remain unclear. In this study, we characterized the functionality and evolution of Sfx, a H-NS homolog encoded by the model IncX2 plasmid R6K. We demonstrate that Sfx, but not chromosomal H-NS, can repress R6K conjugation. Notably, we find evidence of positive selection acting on the ancestral Sfx lineage. Positively selected sites are located in the dimerization, oligomerization, and DNA-binding interfaces, many of which contribute to R6K repression activity-indicating that adaptive evolution drove the functional divergence of Sfx. We additionally show that Sfx can physically interact with various chromosomally encoded proteins, including H-NS, StpA, and Hha. Hha enhances the ability of Sfx to regulate R6K conjugation, suggesting that Sfx retained functionally important interactions with chromosomal silencing proteins. Surprisingly, the loss of Sfx does not negatively affect the stability or dissemination of R6K in laboratory conditions, reflecting the complexity of selective pressures favoring conjugation repression. Overall, our study sheds light on the functional and evolutionary divergence of a plasmid-borne H-NS-like protein, highlighting how these loosely specific DNA-binding proteins evolved to specifically regulate different plasmid functions.IMPORTANCEConjugative plasmids play a crucial role in spreading antimicrobial resistance and virulence genes. Most natural conjugative plasmids conjugate only under specific conditions. Therefore, studying the molecular mechanisms underlying conjugation regulation is essential for understanding antimicrobial resistance and pathogen evolution. In this study, we characterized the conjugation regulation of the model IncX plasmid R6K. We discovered that Sfx, a H-NS homolog carried by the plasmid, represses conjugation. Molecular evolutionary analyses combined with gain-of-function experiments indicate that positive selection underlies the conjugation repression activity of Sfx. Additionally, we demonstrate that the loss of Sfx does not adversely affect R6K maintenance under laboratory conditions, suggesting additional selective forces favoring Sfx carriage. Overall, this work underscores the impact of protein diversification on plasmid biology, enhancing our understanding of how molecular evolution affects broader plasmid ecology.

Subcellular Enrichment Patterns of New Genes in Drosophila Evolution

The evolutionary patterns of proteins within subcellular compartments underlie the innovation and diversification foundation of the living eukaryotic organism. The location of proteins in subcellular compartments promotes the formation of network interaction modules, which in turn reshape the architecture of higher-level protein-protein interaction networks. Here, we conducted the most up-to-date gene age dating of Drosophila melanogaster by employing recently available long-read sequencing genomes as references. We found that an elevated gene fixation in the most recent common ancestor of Drosophila genus predated the divergence between two Drosophila subgenera, and a significant tendency of these genes in D. melanogaster encode proteins that localize to the extracellular matrix, accompanying the adaptive radiation of Drosophila species. Proteins encoded by genes located in the extracellular space exhibit higher sequence divergence, suggesting a rapid evolutionary process. We also observed that proteins encoded by genes originating from the same evolutionary branches tend to co-localize in the same subcellular compartments, and proteins in the same subcellular compartment tend to interact with each other. The proteins encoded by genes that have persisted through deeper branches exhibit broader localization across multiple subcellular compartments, enhancing the likelihood of their integration into various protein or gene regulatory networks, thereby increasing functional diversity. These evolutionary patterns not only contribute to understanding the evolution of subcellular localization in proteins encoded by genes originating from different branches, but also provide insights into the evolution of protein-protein networks driven by the emergence of new genes.

Decoding resilience: ecology, regulation, and evolution of biosynthetic gene clusters

Secondary metabolism is crucial for plant survival and can generate chemistry with nutritional, therapeutic, and industrial value. Biosynthetic genes of selected secondary metabolites cluster within localised chromosomal regions. The arrangement of these biosynthetic gene clusters (BGCs) challenges the long-held model of random gene order in eukaryotes, raising questions about their regulation, ecological significance, and evolution. In this review, we address these questions by exploring the contribution of BGCs to ecologically relevant plant-biotic interactions, while also evaluating the molecular-(epi)genetic mechanisms controlling their coordinated stress- and tissue-specific expression. Based on evidence that BGCs have distinct chromatin signatures and are enriched with transposable elements (TEs), we integrate emerging hypotheses into an updated evolutionary model emphasising how stress-induced epigenetic processes have shaped BGC formation.

DNA methylation dynamics in gymnosperm duplicate genes: implications for genome evolution and stress adaptation

Duplicate genes are pivotal in driving evolutionary innovation, often exhibiting expression divergence that offers a system to investigate the role of DNA methylation in transcriptional regulation. However, previous studies have predominantly focused on angiosperms, leaving the methylation patterns in major lineages of land plants still unclear. This study explores DNA methylation evolution in duplicate genes across representative gymnosperm species with large genomes, spanning over 300 million years, using genomic, transcriptomic, and high-depth DNA methylomic data. We observed variations in DNA methylation levels along gene bodies, flanking regions, and methylation statuses of coding regions across different duplication types. Biased divergences in DNA methylation and gene expression frequently occurred between duplicate copies. Specifically, methylation divergences in the 2-kb downstream regions negatively correlated with gene expression. Both CG and CHG DNA methylation in gene bodies were positively correlated with gene length, suggesting these methylation types may function as an epigenomic buffer to mitigate the adverse impact of gene length on expression. Duplicate genes exhibiting both methylation and expression divergences were notably enriched in adaptation-related biological processes, suggesting that DNA methylation may aid adaptive evolution in gymnosperms by regulating stress response genes. Changes in expression levels correlated with switches in methylation status within coding regions of transposed duplicates. Specifically, depletion for CG methylation or enrichment for non-CG methylation significantly reduced the expression of translocated copies. This correlation suggests that DNA methylation may reduce genetic redundancy by silencing translocated copies. Our study highlights the significance of DNA methylation in plant genome evolution and stress adaptation.

Functionally diversified BiP orthologs control body growth, reproduction, stress resistance, aging, and ER-Phagy in Caenorhabditis elegans

Cellular systems that govern protein folding rely on a delicate balance of functional redundancy and diversification to maintain protein homeostasis (proteostasis). Here, we use Caenorhabditis elegans to demonstrate how both overlapping and divergent activities of two homologous endoplasmic reticulum (ER)-resident HSP70 family chaperones, HSP-3 and HSP-4, orchestrate ER proteostasis and contribute to organismal physiology. We identify tissue-, age-, and stress-specific protein expression patterns and find both redundant and distinct functions for HSP-3 and HSP-4 in ER stress resistance, reproduction, and body size regulation. We show that only HSP-3 overexpression is sufficient to improve longevity and that loss of HSP-3 or HSP-4 during distinct stages of the worm cycle or specific tissues have opposing effects on worm lifespan. Furthermore, we find that loss of HSP-4, but not HSP-3, improves tolerance to protein aggregation induced-stress by activating ER-Phagy through the engagement of IRE-1 and the putative ER-Phagy receptor, C18E9.2. Mechanistically, we show that de-repression of IRE-1 via HSP-4 dissociation allows for direct inhibition of C18E9.2- mediated ER-Phagy and demonstrate that a conserved orthologous mechanism involving the respective human orthologs, BiP, Sec-62, and IRE-1, contributes to ER proteostasis regulation in human cells. Taken as a whole, our study demonstrates that functional diversification of orthologous proteins within a single organelle is an efficient mechanism to maximize stress resilience while also defining a novel link between ER- phagy and proteostasis regulation.

The Rich Evolutionary History of the Reactive Oxygen Species Metabolic Arsenal Shapes Its Mechanistic Plasticity at the Onset of Metazoan Regeneration

Regeneration, the ability to restore body parts after injury, is widespread in metazoans; however, the underlying molecular and cellular mechanisms involved in this process remain largely unknown, and its evolutionary history is consequently unresolved. Recently, reactive oxygen species (ROS) have been shown in several metazoan models to be triggers of apoptosis and cell proliferation that drive regenerative success. However, it is not known whether the contribution of ROS to regeneration relies on conserved mechanisms. Here we performed a comparative genomic analysis of ROS metabolism actors across metazoans, and carried out a comparative study of the deployment and roles of ROS during regeneration in two different metazoan models: the annelid Platynereis dumerilii and the cnidarian Nematostella vectensis. We established that the vast majority of metazoans encode a core redox kit allowing for the production and detoxification of ROS, and overall regulation of ROS levels. However, the precise composition of the redox arsenal can vary significantly from species to species, suggesting that evolutionary constraints apply to ROS metabolism functions rather than precise actors. We found that while ROS are necessary for regeneration in both Platynereis and Nematostella, the two species deploy different enzymatic activities controlling ROS dynamics, and display distinct effects of ROS signaling on injury-induced apoptosis and cell proliferation. We conclude that, while ROS are a common feature of metazoan regeneration, their production and contribution to this phenomenon may depend on different molecular mechanisms highlighting the overall plasticity of the machinery.

The molecular toll pathway repertoire in anopheline mosquitoes

Innate immunity in mosquitoes has received much attention due to its potential impact on vector competence for vector-borne disease pathogens, including malaria parasites. The nuclear factor (NF)-κB-dependent Toll pathway is a major regulator of innate immunity in insects. In mosquitoes, this pathway controls transcription of the majority of the known canonical humoral immune effectors, mediates anti-bacterial, anti-fungal and anti-viral immune responses, and contributes to malaria parasite killing. However, besides initial gene annotation of putative Toll pathway members and genetic analysis of the contribution of few key components to immunity, the molecular make-up and function of the Toll pathway in mosquitoes is largely unexplored. To facilitate functional analyses of the Toll pathway in mosquitoes, we report here manually annotated and refined gene models of Toll-like receptors and all putative components of the intracellular signal transduction cascade across 19 anopheline genomes, and in two culicine genomes. In addition, based on phylogenetic analyses, we identified differing levels of evolutionary constraint across the intracellular Toll pathway members, and identified a recent radiation of TOLL1/5 within the Anopheles gambiae complex. Together, this study provides insight into the evolution of TLRs and the putative members of the intracellular signal transduction cascade within the genus Anopheles.

Plant Nuclear Factor Y (NF-Y) Transcription Factors: Evolving Insights into Biological Functions and Gene Expansion

Gene expansion is a common phenomenon in plant transcription factor families; however, the underlying molecular mechanisms remain elusive. Examples of gene expansion in transcription factors are found in all eukaryotes. One example is plant nuclear factor Y (NF-Y) transcription factors. NF-Y is ubiquitous to eukaryotes and comprises three independent protein families: NF-YA, NF-YB, and NF-YC. While animals and fungi mostly have one of each NF-Y subunit, NF-Y is greatly expanded in plants. For example, humans have one each of NF-YA, NF-YB, and NF-YC, while the model plant Arabidopsis has ten each of NF-YA, NF-YB, and NF-YC. Our understanding of the plant NF-Y, including its biological roles, molecular mechanisms, and gene expansion, has improved over the past few years. Here we will review its biological roles and focus on studies demonstrating that NF-Y can serve as a model for plant gene expansion. These studies show that NF-Y can be classified into ancestrally related subclasses. Further, the primary structure of each NF-Y contains a conserved core domain flanked by non-conserved N- and C-termini. The non-conserved N- and C-termini, under pressure for diversifying selection, may provide clues to this gene family's retention and functional diversification following gene duplication. In summary, this review demonstrates that NF-Y expansion has the potential to be used as a model to study the gene expansion and retention of transcription factor families.

Identifying candidate genes for sugar accumulation in sugarcane: an integrative approach

BACKGROUND

Elucidating the intricacies of the sugarcane genome is essential for breeding superior cultivars. This economically important crop originates from hybridizations of highly polyploid Saccharum species. However, the large size (10 Gb), high degree of polyploidy, and aneuploidy of the sugarcane genome pose significant challenges to complete genome sequencing, assembly, and annotation. One successful strategy for identifying candidate genes linked to agronomic traits, particularly those associated with sugar accumulation, leverages synteny and potential collinearity with related species.

RESULTS

In this study, we explored synteny between sorghum and sugarcane. Genes from a sorghum Brix QTL were used to screen bacterial artificial chromosome (BAC) libraries from two Brazilian sugarcane varieties (IACSP93-3046 and SP80-3280). The entire region was successfully recovered, confirming synteny and collinearity between the species. Manual annotation identified 51 genes in the hybrid varieties that were subsequently confirmed to be present in Saccharum spontaneum. This study employed a multifaceted approach to identify candidate genes for sugar accumulation, including retrieving the genomic region of interest, performing a gene-by-gene analysis, analyzing RNA-seq data for internodes from Saccharum officinarum and S. spontaneum accessions, constructing a coexpression network to examine the expression patterns of genes within the studied region and their neighbors, and finally identifying differentially expressed genes (DEGs).

CONCLUSIONS

This comprehensive approach led to the discovery of three candidate genes potentially involved in sugar accumulation: an ethylene-responsive transcription factor (ERF), an ABA 8'-hydroxylase, and a prolyl oligopeptidase (POP). These findings could be valuable for identifying additional candidate genes for other important agricultural traits and directly targeting candidate genes for further work in molecular breeding.

Genome-Wide Identification and Evolutionary Analysis of Functional BBM-like Genes in Plant Species

Background/Objectives: BABY BOOM (BBM), a transcription factor from the APETALA2 (AP2) protein family, plays a critical role in somatic embryo induction and apomixis. BBM has now been widely applied to induce apomixis or enhance plant transformation and regeneration efficiency through overexpression or ectopic expression. However, the structural and functional evolutionary history of BBM genes in plants is still not well understood. Methods: The protein sequences of 10 selected plant species were used to locate the branch of BBM-Like by key domain identification and phylogenetic tree construction. The identified BBML genes were used for further conserved motif identification, gene structural analysis, miRNA binding site prediction, cis-acting element prediction, collinear analysis, protein-protein interaction network construction, three-dimensional structure modeling, molecular docking, and expression pattern analysis. Results: A total of 24 BBML proteins were identified from 10 representative plant species. Phylogenetic relationship analysis displayed that BBML proteins from eudicots and monocots were divided into two clusters, with monocots exhibiting a higher number of BBMLs. Gene duplication events indicated that whole genome/segmental duplication were the primary drivers of BBML genes' evolution in the tested species, with purifying selection playing a key role during evolution processes. Comparative analysis of motif, domains, and gene structures revealed that most BBMLs were highly evolutionarily conserved. The expression patterns of BBML genes revealed significant tissue specificity, particularly in the root and embryo. We also constructed protein-protein interaction networks and molecular docking models to identify functional pathways and key amino acid residues of BBML proteins. The functions of BBMLs may differ between monocots and eudicots, as suggested by the functional enrichment of interacting proteins. Conclusions: Our research delved into the molecular mechanism, evolutionary relationships, functional differentiation, and expression patterns of BBML genes across plants, laying the groundwork for further investigations into the molecular properties and biological roles of BBMLs.

Asymmetrical Evolution of Promoter Methylation of Mammalian Genes after Duplication

Even though gene duplication is a key source of new genes and evolutionary innovation, it is unclear how duplicates survive the period immediately following gene duplication, in which both copies are functionally redundant. In the absence of epigenetic silencing, the abundance of the gene product would double after gene duplication, which would often have deleterious effects. However, recent duplicates exhibit low expression levels, which could be at least partially explained by high levels of promoter methylation. What evolutionary paths lead to duplicate hypermethylation, and does it affect both duplicates or only one? Here, we compare levels of promoter methylation in 10 human and 16 mouse tissues, between singletons and duplicates and among human-mouse orthologs of different kinds (one-to-one, one-to-many, many-to-one, and many-to-many). Our results indicate that: (i) on average, duplicates are more methylated than singletons in mouse, but less methylated than singletons in human, (ii) recently duplicated genes tend to exhibit high levels of promoter methylation, (iii) genes that undergo duplication tend to be highly methylated before duplication, (iv) after gene duplication, one of the copies (the daughter copy, i.e. the one that relocates to a new genomic context) tends to undergo an additional increase in promoter methylation, whereas the other (the parental copy, which remains in the original genomic location) tends to retain preduplication methylation levels, and (v) daughter copies tend to be lowly expressed. These observations support a model in which daughter copies are repressed via promoter hypermethylation and can thus survive the filter of purifying selection until both copies diverge functionally.

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.

Forward-in-time simulation of chromosomal rearrangements: The invisible backbone that sustains long-term adaptation

While chromosomal rearrangements are ubiquitous in all domains of life, very little is known about their evolutionary significance, mostly because, apart from a few specifically studied and well-documented mechanisms (interaction with recombination, gene duplication, etc.), very few models take them into account. As a consequence, we lack a general theory to account for their direct and indirect contributions to evolution. Here, we propose Aevol, a forward-in-time simulation platform specifically dedicated to unravelling the evolutionary significance of chromosomal rearrangements (CR) compared to local mutations (LM). Using the platform, we evolve populations of organisms in four conditions characterized by an increasing diversity of mutational operators-from substitutions alone to a mix of substitutions, InDels and CR-but with a constant global mutational rate. Despite being almost invisible in the phylogeny owing to the scarcity of their fixation in the lineages, we show that CR make a decisive contribution to the evolutionary dynamics by comparing the outcome in these four conditions. As expected, chromosomal rearrangements allow fast expansion of the gene repertoire through gene duplication, but they also reduce the effect of diminishing-returns epistasis, hence sustaining adaptation on the long-run. At last, we show that chromosomal rearrangements tightly regulate the size of the genome through indirect selection for reproductive robustness. Overall, these results confirm the need to improve our theoretical understanding of the contribution of chromosomal rearrangements to evolution and show that dedicated platforms like Aevol can efficiently contribute to this agenda.

Evolutionary analysis of Quinone Reductases 1 and 2 suggests that NQO2 evolved to function as a pseudoenzyme

Quinone reductases 1 and 2 (NQO1 and NQO2) are paralogous FAD-linked enzymes found in all amniotes. NQO1 and NQO2 have similar structures, and both catalyze the reduction of quinones and other electrophiles; however, the two enzymes differ in their cosubstrate preference. While NQO1 can use both redox couples NADH and NADPH, NQO2 is almost inactive with these cosubstrates and instead must use dihydronicotinamide riboside (NRH) and small synthetic cosubstrates such as N-benzyl-dihydronicotinamide (BNAH) for efficient catalysis. We used ancestral sequence reconstruction to investigate the catalytic properties of a predicted common ancestor and two additional ancestors from each of the evolutionary pathways to extant NQO1 and NQO2. In all cases, the small nicotinamide cosubstrates NRH and BNAH were good cosubstrates for the common ancestor and the enzymes along both the NQO1 and NQO2 lineages. In contrast, with NADH as cosubstrate, extant NQO1 evolved to a catalytic efficiency 100 times higher than the common ancestor, while NQO2 has evolved to a catalytic efficiency 3000 times lower than the common ancestor. The evolutionary analysis combined with site-directed mutagenesis revealed a potential site of interaction for the ADP portion of NAD(P)H in NQO1 that is altered in charge and structure in NQO2. The results indicate that while NQO1 evolved to have greater efficiency with NAD(P)H, befitting an enzymatic function in cells, NQO2 was under selective pressure to acquire extremely low catalytic efficiency with NAD(P)H. These divergent trajectories have implications for the functions of both enzymes.

Benchmarking recent computational tools for DNA-binding protein identification

Identification of DNA-binding proteins (DBPs) is a crucial task in genome annotation, as it aids in understanding gene regulation, DNA replication, transcriptional control, and various cellular processes. In this paper, we conduct an unbiased benchmarking of 11 state-of-the-art computational tools as well as traditional tools such as ScanProsite, BLAST, and HMMER for identifying DBPs. We highlight the data leakage issue in conventional datasets leading to inflated performance. We introduce new evaluation datasets to support further development. Through a comprehensive evaluation pipeline, we identify potential limitations in models, feature extraction techniques, and training methods, and recommend solutions regarding these issues. We show that combining the predictions of the two best computational tools with BLAST-based prediction significantly enhances DBP identification capability. We provide this consensus method as user-friendly software. The datasets and software are available at https://github.com/Rafeed-bot/DNA_BP_Benchmarking.

Evolution of aromatic amino acid metabolism in plants: a key driving force behind plant chemical diversity in aromatic natural products

A diverse array of plant aromatic compounds contributes to the tremendous chemical diversity in the plant kingdom that cannot be seen in microbes or animals. Such chemodiversity of aromatic natural products has emerged, occasionally in a lineage-specific manner, to adopt to challenging environmental niches, as various aromatic specialized metabolites play indispensable roles in plant development and stress responses (e.g. lignin, phytohormones, pigments and defence compounds). These aromatic natural products are synthesized from aromatic amino acids (AAAs), l-tyrosine, l-phenylalanine and l-tryptophan. While amino acid metabolism is generally assumed to be conserved between animals, microbes and plants, recent phylogenomic, biochemical and metabolomic studies have revealed the diversity of the AAA metabolism that supports efficient carbon allocation to downstream biosynthetic pathways of AAA-derived metabolites in plants. This review showcases the intra- and inter-kingdom diversification and origin of committed enzymes involved in plant AAA biosynthesis and catabolism and their potential application as genetic tools for plant metabolic engineering. I also discuss evolutionary trends in the diversification of plant AAA metabolism that expands the chemical diversity of AAA-derived aromatic natural products in plants. This article is part of the theme issue 'The evolution of plant metabolism'.

Evolution of plant metabolism: the state-of-the-art

Immense chemical diversity is one of the hallmark features of plants. This chemo-diversity is mainly underpinned by a highly complex and biodiverse biochemical machinery. Plant metabolic enzymes originated and were inherited from their eukaryotic and prokaryotic ancestors and further diversified by the unprecedentedly high rates of gene duplication and functionalization experienced in land plants. Unlike prokaryotic microbes, which display frequent horizontal gene transfer events and multiple inputs of energy and organic carbon, land plants predominantly rely on organic carbon generated from CO2 and have experienced relatively few gene transfers during their recent evolutionary history. As such, plant metabolic networks have evolved in a stepwise manner using existing networks as a starting point and under various evolutionary constraints. That said, until recently, the evolution of only a handful of metabolic traits had been extensively investigated and as such, the evolution of metabolism has received a fraction of the attention of, the evolution of development, for example. Advances in metabolomics and next-generation sequencing have, however, recently led to a deeper understanding of how a wide range of plant primary and specialized (secondary) metabolic pathways have evolved both as a consequence of natural selection and of domestication and crop improvement processes. This article is part of the theme issue 'The evolution of plant metabolism'.

The Two Sides of Indoleamine 2,3-Dioxygenase 2 (IDO2)

Indoleamine 2,3-dioxygenase 1 (IDO1) and IDO2 originated from gene duplication before vertebrate divergence. While IDO1 has a well-defined role in immune regulation, the biological role of IDO2 remains unclear. Discovered in 2007, IDO2 is located near the IDO1 gene. Because of their high sequence similarity, IDO2 was initially thought to be a tryptophan (Trp)-degrading enzyme like IDO1. Differently from what expected, IDO2 displays extremely low catalytic activity toward Trp. Nevertheless, many studies, often contradictory, have tried to demonstrate that IDO2 modulates immune responses by catabolizing Trp into kynurenine, an unconvincing hypothesis linked to an incomplete understanding of IDO2's activity. In this study, we review IDO2's functional role beyond Trp metabolism. IDO2's evolutionary persistence across species, despite being almost inactive as an enzyme, suggests it has some relevant biological importance. IDO2 expression in human normal cells is poor, but significant in various cancers, with two prevalent SNPs. Overall, the comparison of IDO2 to IDO1 as a Trp-degrading enzyme may have led to misunderstandings about IDO2's true physiological and pathological roles. New insights suggest that IDO2 might function more as a signaling molecule, particularly in cancer contexts, and further studies could reveal its potential as a target for cancer therapy.

Genome-Wide Characterization of Solanum tuberosum CCoAOMT Gene Family and Identification of StCCoAOMT Genes Involved in Anthocyanin Biosynthesis

BACKGROUND

The caffeoyl-CoA-O methyltransferase (CCoAOMT) family plays essential roles in the methylation of various secondary metabolites, including anthocyanins. Despite the wide identification of the CCoAOMT family in plants, the characterization and function of CCoAOMT protein members in Solanum tuberosum remain poorly understood.

METHODS AND RESULTS

In this study, a total of 12 StCCoAOMT members were identified in the genome of S. tuberosum using the Blastp and HMM search and were unevenly located on eight chromosomes. Collinearity analysis revealed that four tandem duplicated gene pairs and two segmental duplicated gene pairs existed in the S. tuberosum genome, demonstrating that duplication events play a key role in the expansion of the CCoAOMT family. All StCCoAOMTs were clustered into group I and group II based on phylogenetic analysis, which was further verified by the conserved motifs and gene structures analysis. The cis-acting elements analysis illustrated that StCCoAOMTs might be important for photosynthesis, hormone responses, and abiotic stress. Expression analysis demonstrated that StCCoAOMT genes have diverse transcript levels in various tissues and that StCCoAOMT10 was significantly expressed in purple potatoes with abundant anthocyanin content according to RNA-seq data and qRT-PCR assays. In addition, the subcellular localization assay validated that the StCCoAOMT10 protein was mainly localized in the cytoplasm and nucleus.

CONCLUSIONS

These results will be of great importance for a better understanding of the features of CCoAOMT family members, especially of the candidate genes involved in the methylation of anthocyanins in S. tuberosum, and also for improving the nutritional quality of S. tuberosum.

Discovery of Five Classes of Bacterial Defensins: Ancestral Precursors of Defensins from Eukarya?

Defensins are present in many organisms and are divided into two evolutionary groups, termed cis- and trans-defensins. Cis-defensins have only recently been reported in bacteria, and knowledge of these defensins is limited, with no family classification. Here, we describe the identification of 74 cis-defensins from bacteria and propose five classes for their classification. We also report the first NMR structure determination of a Myxoccocus xanthus defensin, as well as its in silico expression analysis. Xanthusin-1 has a unique structure among the published defensins, which could indicate that the proposed class II peptides constitute a separate group of defensins. Xanthusin-1 gene expression was observed in casitone-based and Streptomyces coelicolor coculture-grown media. Our results demonstrate a wider distribution of defensins outside the Eukarya domain, shedding light on the origin and distribution of defensins. The sharing of three disulfide defensins between bacteria and eukaryotes points to a possible prokaryotic origin of the CSαβ motif. Moreover, the identification of defensins in Gram-positive and Gram-negative bacteria indicates an early origin but with many gene losses during the evolutionary process, similar to findings for eukaryotic defensins.

Genome-Wide Analysis of the Common Fig (Ficus carica L.) R2R3-MYB Genes Reveals Their Structure, Evolution, and Roles in Fruit Color Variation

The R2R3-MYB transcription factor (TF) family is crucial for regulating plant growth, stress response, and fruit ripening. Although this TF family has been examined in a multitude of plants, the R2R3-MYB TFs in Ficus carica, a Mediterranean fruit species, have yet to be characterized. This study identified and classified 63 R2R3-MYB genes (FcMYB1 to FcMYB63) in the F. carica genome. We analyzed these genes for physicochemical properties, conserved motifs, phylogenetic relationships, gene architecture, selection pressure, and gene expression profiles and networks. The genes were classified into 29 clades, with members of the same clade showing similar exon-intron structures and motif compositions. Of the 54 orthologous gene pairs shared with mulberry (Morus notabilis), 52 evolved under negative selection, while two pairs (FcMYB55/MnMYB20 and FcMYB59/MnMYB31) experienced diversifying selection. RNA-Seq analysis showed that FcMYB26, FcMYB33, and FcMYB34 were significantly overexpressed in fig fruit peel during maturation phase III. Weighted gene co-expression network analysis (WGCNA) indicated that these genes are part of an expression module associated with the anthocyanin pathway. RT-qPCR validation confirmed these findings and revealed that the Tunisian cultivars 'Zidi' and 'Soltani' have cultivar-specific R2R3-FcMYB genes highly overexpressed during the final stage of fruit maturation and color acquisition. These genes likely influence cultivar-specific pigment synthesis. This study provides a comprehensive overview of the R2R3-MYB TF family in fig, offering a framework for selecting genes related to fruit peel color in breeding programs.

Molecular identification and expression patterns of sweet cherry HIPPs and functional analysis of PavHIPP16 in cold stress

MAIN CONCLUSION

The HIPP proteins are involved in low-temperature stress, the growth of sweet cherry, and may be potential targets for genetic improvement. PavHIPP16 improved cold resistance in Arabidopsis. In response to abiotic stressors, the heavy metal-associated isoprenylated plant protein (HIPP) proteins play a crucial regulatory role. Although the function of HIPP has been identified in some plants, there have been fewer systematic studies conducted on sweet cherry (Prunus avium L.). Therefore, we performed a comprehensive analysis and expression profiling of PavHIPPs using bioinformatics, RT-PCR, and qRT-PCR techniques. Our findings revealed that 28 PavHIPP genes were unevenly distributed across eight chromosomes. We predicted nine motifs in PavHIPP proteins and observed similar gene structures among highly homologous proteins. The promoter sequences of PavHIPPs contained numerous regulatory elements associated with an adversity of stress. The expression levels of some members showed varying degrees of change under low-temperature treatment. These genes were differentially expressed during flower and fruit development. Arabidopsis overexpressing the PavHIPP16 (OE) gene showed significantly lower relative conductivity and malondialdehyde (MDA) content compared with the wild-type (WT) plants under cold environment. Conversely, peroxidase (POD) activity, superoxide dismutase (SOD) activity, and proline content were significantly higher in OE Arabidopsis than in WT plants. Overall, our results suggest that PavHIPP16 OE Arabidopsis thaliana exhibited enhanced adaptability compared to WT plants under cold conditions. This study provides a foundation for future investigations of the pathways regulating sweet cherry growth and development mediated by the HIPP genes.

Global analysis of neuropeptide receptor conservation across phylum Nematoda

BACKGROUND

The phylum Nematoda is incredibly diverse and includes many parasites of humans, livestock, and plants. Peptide-activated G protein-coupled receptors (GPCRs) are central to the regulation of physiology and numerous behaviors, and they represent appealing pharmacological targets for parasite control. Efforts are ongoing to characterize the functions and define the ligands of nematode GPCRs, with already most peptide GPCRs known or predicted in Caenorhabditis elegans. However, comparative analyses of peptide GPCR conservation between C. elegans and other nematode species are limited, and many nematode GPCRs remain orphan. A phylum-wide perspective on peptide GPCR profiles will benefit functional and applied studies of nematode peptide GPCRs.

RESULTS

We constructed a pan-phylum resource of C. elegans peptide GPCR orthologs in 125 nematode species using a semi-automated pipeline for analysis of predicted proteome datasets. The peptide GPCR profile varies between nematode species of different phylogenetic clades and multiple C. elegans peptide GPCRs have orthologs across the phylum Nematoda. We identified peptide ligands for two highly conserved orphan receptors, NPR-9 and NPR-16, that belong to the bilaterian galanin/allatostatin A (Gal/AstA) and somatostatin/allatostatin C (SST/AstC) receptor families. The AstA-like NLP-1 peptides activate NPR-9 in cultured cells and are cognate ligands of this receptor in vivo. In addition, we discovered an AstC-type peptide, NLP-99, that activates the AstC-type receptor NPR-16. In our pan-phylum resource, the phylum-wide representation of NPR-9 and NPR-16 resembles that of their cognate ligands more than those of allatostatin-like peptides that do not activate these receptors.

CONCLUSIONS

The repertoire of C. elegans peptide GPCR orthologs varies across phylogenetic clades and several peptide GPCRs show broad conservation in the phylum Nematoda. Our work functionally characterizes the conserved receptors NPR-9 and NPR-16 as the respective GPCRs for the AstA-like NLP-1 peptides and the AstC-related peptide NLP-99. NLP-1 and NLP-99 are widely conserved in nematodes and their representation matches that of their receptor in most species. These findings demonstrate the conservation of a functional Gal/AstA and SST/AstC signaling system in nematodes. Our dataset of C. elegans peptide GPCR orthologs also lays a foundation for further functional studies of peptide GPCRs in the widely diverse nematode phylum.

Dosage sensitivity shapes balanced expression and gene longevity of homoeologs after whole-genome duplications in angiosperms

Following whole-genome duplication (WGD), duplicate gene pairs (homoeologs) can evolve varying degrees of expression divergence. However, the determinants influencing these relative expression level differences (RFPKM) between homoeologs remain elusive. In this study, we analyzed the RFPKM between homoeologs in 3 angiosperms, Nymphaea colorata, Nelumbo nucifera, and Acorus tatarinowii, all having undergone a single WGD since the origin of angiosperms. Our results show significant positive correlations in RFPKM of homoeologs among tissues within the same species, and among orthologs across these 3 species, indicating convergent expression balance/bias between homoeologous gene copies following independent WGDs. We linked RFPKM between homoeologs to gene attributes associated with dosage-balance constraints, such as protein-protein interactions, lethal-phenotype scores in Arabidopsis (Arabidopsis thaliana) orthologs, domain numbers, and expression breadth. Notably, homoeologs with lower RFPKM often had more interactions and higher lethal-phenotype scores, indicating selective pressures favoring balanced expression. Also, homoeologs with lower RFPKM were more likely to be retained after WGDs in angiosperms. Within Nelumbo, greater RFPKM between homoeologs correlated with increased cis- and trans-regulatory differentiation between species, highlighting the ongoing escalation of gene expression divergence. We further found that expression degeneration in 1 copy of homoeologs is inclined toward nonfunctionalization. Our research highlights the importance of balanced expression, shaped by dosage-balance constraints, in the evolutionary retention of homoeologs in plants.

A developmental gradient reveals biosynthetic pathways to eukaryotic toxins in monocot geophytes

Numerous eukaryotic toxins that accumulate in geophytic plants are valuable in the clinic, yet their biosynthetic pathways have remained elusive. A notable example is the >150 Amaryllidaceae alkaloids (AmAs), including galantamine, an FDA-approved treatment for Alzheimer's disease. We show that while AmAs accumulate to high levels in many daffodil tissues, biosynthesis is localized to nascent, growing tissue at the leaf base. A similar trend is found in the production of steroidal alkaloids (e.g., cyclopamine) in corn lily. This model of active biosynthesis enabled the elucidation of a complete set of biosynthetic genes that can be used to produce AmAs. Taken together, our work sheds light on the developmental and enzymatic logic of diverse alkaloid biosynthesis in daffodils. More broadly, it suggests a paradigm for biosynthesis regulation in monocot geophytes, where plants are protected from herbivory through active charging of newly formed cells with eukaryotic toxins that persist as above-ground tissue develops.

Evolution of Sensory Receptors

Sensory receptors are at the interface between an organism and its environment and thus represent key sites for biological innovation. Here, we survey major sensory receptor families to uncover emerging evolutionary patterns. Receptors for touch, temperature, and light constitute part of the ancestral sensory toolkit of animals, often predating the evolution of multicellularity and the nervous system. In contrast, chemoreceptors exhibit a dynamic history of lineage-specific expansions and contractions correlated with the disparate complexity of chemical environments. A recurring theme includes independent transitions from neurotransmitter receptors to sensory receptors of diverse stimuli from the outside world. We then provide an overview of the evolutionary mechanisms underlying sensory receptor diversification and highlight examples where signatures of natural selection are used to identify novel sensory adaptations. Finally, we discuss sensory receptors as evolutionary hotspots driving reproductive isolation and speciation, thereby contributing to the stunning diversity of animals.

Zebrafish models of human-duplicated SRGAP2 reveal novel functions in microglia and visual system development

The expansion of the human SRGAP2 family, resulting in a human-specific paralog SRGAP2C, likely contributed to altered evolutionary brain features. The introduction of SRGAP2C in mouse models is associated with changes in cortical neuronal migration, axon guidance, synaptogenesis, and sensory-task performance. Truncated SRGAP2C heterodimerizes with the full-length ancestral gene product SRGAP2A and antagonizes its functions. However, the significance of SRGAP2 duplication beyond neocortex development has not been elucidated due to the embryonic lethality of complete Srgap2 knockout in mice. Using zebrafish, we show that srgap2 knockout results in viable offspring and that these larvae phenocopy "humanized" SRGAP2C larvae, including altered morphometric features (i.e., reduced body length and inter-eye distance) and differential expression of synapse-, axonogenesis-, and vision-related genes. Through single-cell transcriptome analysis, we demonstrate a skewed balance of excitatory and inhibitory neurons that likely contribute to increased susceptibility to seizures displayed by Srgap2 mutant larvae, a phenotype resembling SRGAP2 loss-of-function in a child with early infantile epileptic encephalopathy. Single-cell data also shows strong endogenous expression of srgap2 in microglia with mutants exhibiting altered membrane dynamics and likely delayed maturation of microglial cells. Microglia cells expressing srgap2 were also detected in the developing eye together with altered expression of genes related to axonogenesis in mutant retinal cells. Consistent with the perturbed gene expression in the retina, we found that SRGAP2 mutant larvae exhibited increased sensitivity to broad and fine visual cues. Finally, comparing the transcriptomes of relevant cell types between human (+SRGAP2C) and non-human primates (-SRGAP2C) revealed significant overlaps of gene alterations with mutant cells in our zebrafish models; this suggests that SRGAP2C plays a similar role altering microglia and the visual system in modern humans. Together, our functional characterization of conserved ortholog Srgap2 and human SRGAP2C in zebrafish uncovered novel gene functions and highlights the strength of cross-species analysis in understanding the development of human-specific features.

Innovation in viruses: fitness valley crossing, neutral landscapes, or just duplications?

Viruses evolve by periods of relative stasis interleaved with sudden, rapid series of mutation fixations, known as evolutionary bursts. These bursts can be triggered by external factors, such as environmental changes, antiviral therapies, or spill-overs from reservoirs into novel host species. However, it has also been suggested that bursts may result from the intrinsic evolutionary dynamics of viruses. Indeed, bursts could be caused by fitness valley crossing, or a neutral exploration of a fitness plateau until an escape mutant is found. In order to investigate the importance of these intrinsic causes of evolutionary bursts, we used a simulation software package to perform massive evolution experiments of viral-like genomes. We tested two conditions: (i) after an external change and (ii) in a constant environment, with the latter condition guaranteeing the absence of an external triggering factor. As expected, an external change was almost systematically followed by an evolutionary burst. However, we also observed bursts in the constant environment as well, albeit much less frequently. We analyzed how many of these bursts are triggered by deleterious, quasi-neutral, or beneficial mutations and show that, while bursts can occasionally be triggered by valley crossing or traveling along neutral ridges, many of them were triggered by chromosomal rearrangements and, in particular, segmental duplications. Our results suggest that combinatorial differences between the different mutation types lead to punctuated evolutionary dynamics, with long periods of stasis occasionally interrupted by short periods of rapid evolution, akin to what is observed in virus evolution.

The Molecular Toll Pathway Repertoire in Anopheline Mosquitoes

Innate immunity in mosquitoes has received much attention due to its potential impact on vector competence for vector-borne disease pathogens, including malaria parasites. The nuclear factor (NF)-κB-dependent Toll pathway is a major regulator of innate immunity in insects. In mosquitoes, this pathway controls transcription of the majority of the known canonical humoral immune effectors, mediates anti-bacterial, anti-fungal and anti-viral immune responses, and contributes to malaria parasite killing. However, besides initial gene annotation of putative Toll pathway members and genetic analysis of the contribution of few key components to immunity, the molecular make-up and function of the Toll pathway in mosquitoes is largely unexplored. To facilitate functional analyses of the Toll pathway in mosquitoes, we report here manually annotated and refined gene models of Toll-like receptors and all putative components of the intracellular signal transduction cascade across 19 anopheline genomes, and in two culicine genomes. In addition, based on phylogenetic analyses, we identified differing levels of evolutionary constraint across the intracellular Toll pathway members, and identified a recent radiation of TOLL1/5 within the An. gambiae complex. Together, this study provides insight into the evolution of TLRs and the putative members of the intracellular signal transduction cascade within the genus Anopheles.

A chromosome-level genome assembly of the disco clam, Ctenoides ales

The bivalve subclass Pteriomorphia, which includes the economically important scallops, oysters, mussels, and ark clams, exhibits extreme ecological, morphological, and behavioral diversity. Among this diversity are five morphologically distinct eye types, making Pteriomorphia an excellent setting to explore the molecular basis for the evolution of novel traits. Of pteriomorphian bivalves, Limida is the only order lacking genomic resources, greatly limiting the potential phylogenomic analyses related to eyes and phototransduction. Here, we present a limid genome assembly, the disco clam, Ctenoides ales (C. ales), which is characterized by invaginated eyes, exceptionally long tentacles, and a flashing light display. This genome assembly was constructed with PacBio long reads and Dovetail Omni-CTM proximity-ligation sequencing. The final assembly is ∼2.3Gb and over 99% of the total length is contained in 18 pseudomolecule scaffolds. We annotated 41,064 protein coding genes and reported a BUSCO completeness of 91.9% for metazoa_obd10. Additionally, we report a complete and annotated mitochondrial genome, which also had been lacking from Limida. The ∼20Kb mitogenome has 12 protein coding genes, 22 tRNAs, 2 rRNA genes, and a 1,589 bp duplicated sequence containing the origin of replication. The C. ales nuclear genome size is substantially larger than other pteriomorphian genomes, mainly accounted for by transposable element sequences. We inventoried the genome for opsins, the signaling proteins that initiate phototransduction, and found that, unlike its closest eyed-relatives, the scallops, C. ales lacks duplication of the rhabdomeric Gq-protein-coupled opsin that is typically used for invertebrate vision. In fact, C. ales has uncharacteristically few opsins relative to the other pteriomorphian families, all of which have unique expansions of xenopsins, a recently discovered opsin subfamily. This chromosome-level assembly, along with the mitogenome, is a valuable resource for comparative genomics and phylogenetics in bivalves and particularly for the understudied but charismatic limids.

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.

Horizontally transferred glycoside hydrolase 26 may aid hemipteran insects in plant tissue digestion

Glycoside hydrolases are enzymes that break down complex carbohydrates into simple sugars by catalyzing the hydrolysis of glycosidic bonds. There have been multiple instances of adaptive horizontal gene transfer of genes belonging to various glycoside hydrolase families from microbes to insects, as glycoside hydrolases can metabolize constituents of the carbohydrate-rich plant cell wall. In this study, we characterize the horizontal transfer of a gene from the glycoside hydrolase family 26 (GH26) from bacteria to insects of the order Hemiptera. Our phylogenies trace the horizontal gene transfer to the common ancestor of the superfamilies Pentatomoidea and Lygaeoidea, which include stink bugs and seed bugs. After horizontal transfer, the gene was assimilated into the insect genome as indicated by the gain of an intron, and a eukaryotic signal peptide. Subsequently, the gene has undergone independent losses and expansions in copy number in multiple lineages, suggesting an adaptive role of GH26s in some insects. Finally, we measured tissue-level gene expression of multiple stink bugs and the large milkweed bug using publicly available RNA-seq datasets. We found that the GH26 genes are highly expressed in tissues associated with plant digestion, especially in the principal salivary glands of the stink bugs. Our results are consistent with the hypothesis that this horizontally transferred GH26 was co-opted by the insect to aid in plant tissue digestion and that this HGT event was likely adaptive.

The link between gene duplication and divergent patterns of gene expression across a complex life cycle

The diversification of many lineages throughout natural history has frequently been associated with evolutionary changes in life cycle complexity. However, our understanding of the processes that facilitate differentiation in the morphologies and functions expressed by organisms throughout their life cycles is limited. Theory suggests that the expression of traits is decoupled across life stages, thus allowing for their evolutionary independence. Although trait decoupling between stages is well established, explanations of how said decoupling evolves have seldom been considered. Because the different phenotypes expressed by organisms throughout their life cycles are coded for by the same genome, trait decoupling must be mediated through divergence in gene expression between stages. Gene duplication has been identified as an important mechanism that enables divergence in gene function and expression between cells and tissues. Because stage transitions across life cycles require changes in tissue types and functions, we investigated the potential link between gene duplication and expression divergence between life stages. To explore this idea, we examined the temporal changes in gene expression across the monarch butterfly (Danaus plexippus) metamorphosis. We found that within homologous groups, more phylogenetically diverged genes exhibited more distinct temporal expression patterns. This relationship scaled such that more phylogenetically diverse homologous groups showed more diverse patterns of gene expression. Furthermore, we found that duplicate genes showed increased stage-specificity relative to singleton genes. Overall, our findings suggest an important link between gene duplication and the evolution of complex life cycles.

Facilitating stable gene integration expression and copy number amplification in Bacillus subtilis through a reversible homologous recombination switch

Strengthening the expression level of integrated genes on the genome is crucial for consistently expressing key enzymes in microbial cell factories for efficient bioproduction in synthetic biology. In comparison to plasmid-based multi-copy expression, the utilization of chromosomal multi-copy genes offers increased stability of expression level, diminishes the metabolic burden on host cells, and enhances overall genetic stability. In this study, we developed the "BacAmp", a stabilized gene integration expression and copy number amplification system for high-level expression in Bacillus subtilis, which was achieved by employing a combination of repressor and non-natural amino acids (ncAA)-dependent expression system to create a reversible switch to control the key gene recA for homologous recombination. When the reversible switch is turned on, genome editing and gene amplification can be achieved. Subsequently, the reversible switch was turned off therefore stabilizing the gene copy number. The stabilized gene amplification system marked by green fluorescent protein, achieved a 3-fold increase in gene expression by gene amplification and maintained the average gene copy number at 10 after 110 generations. When we implemented the gene amplification system for the regulation of N-acetylneuraminic acid (NeuAc) synthesis, the copy number of the critical gene increased to an average of 7.7, which yielded a 1.3-fold NeuAc titer. Our research provides a new avenue for gene expression in synthetic biology and can be applied in metabolic engineering in B. subtilis.

Addressing epistasis in the design of protein function

Mutations in protein active sites can dramatically improve function. The active site, however, is densely packed and extremely sensitive to mutations. Therefore, some mutations may only be tolerated in combination with others in a phenomenon known as epistasis. Epistasis reduces the likelihood of obtaining improved functional variants and dramatically slows natural and lab evolutionary processes. Research has shed light on the molecular origins of epistasis and its role in shaping evolutionary trajectories and outcomes. In addition, sequence- and AI-based strategies that infer epistatic relationships from mutational patterns in natural or experimental evolution data have been used to design functional protein variants. In recent years, combinations of such approaches and atomistic design calculations have successfully predicted highly functional combinatorial mutations in active sites. These were used to design thousands of functional active-site variants, demonstrating that, while our understanding of epistasis remains incomplete, some of the determinants that are critical for accurate design are now sufficiently understood. We conclude that the space of active-site variants that has been explored by evolution may be expanded dramatically to enhance natural activities or discover new ones. Furthermore, design opens the way to systematically exploring sequence and structure space and mutational impacts on function, deepening our understanding and control over protein activity.

R2R3-MYB Gene Family in Coptis teeta Wall.: Genome-Wide Identification, Phylogeny, Evolutionary Expansion, and Expression Analyses during Floral Development

The R2R3-MYB gene family represents a widely distributed class of plant transcription factors. This gene family plays an important role in many aspects of plant growth and development. However, the characterization of R2R3-MYB genes present in the genome of Coptis teeta has not been reported. Here, we describe the bioinformatic identification and characterization of 88 R2R3-MYB genes in this species, and the identification of members of the R2R3-MYB gene family in species within the order Ranales most closely related to Coptis teeta. The CteR2R3-MYB genes were shown to exhibit a higher degree of conservation compared to those of A. thaliana, as evidenced by phylogeny, conserved motifs, gene structure, and replication event analyses. Cis-acting element analysis confirmed the involvement of CteR2R3-MYB genes in a variety of developmental processes, including growth, cell differentiation, and reproduction mediated by hormone synthesis. In addition, through homology comparisons with the equivalent gene family in A. thaliana, protein regulatory network prediction and transcriptome data analysis of floral organs across three time periods of flower development, 17 candidate genes were shown to exhibit biased expression in two floral phenotypes of C. teeta. This suggests their potential involvement in floral development (anther development) in this species.

Strategic targeting of Cas9 nickase induces large segmental duplications

Gene/segmental duplications play crucial roles in genome evolution and variation. Here, we introduce paired nicking-induced amplification (PNAmp) for their experimental induction. PNAmp strategically places two Cas9 nickases upstream and downstream of a replication origin on opposite strands. This configuration directs the sister replication forks initiated from the origin to break at the nicks, generating a pair of one-ended double-strand breaks. If homologous sequences flank the two break sites, then end resection converts them to single-stranded DNAs that readily anneal to drive duplication of the region bounded by the homologous sequences. PNAmp induces duplication of segments as large as ∼1 Mb with efficiencies exceeding 10% in the budding yeast Saccharomyces cerevisiae. Furthermore, appropriate splint DNAs allow PNAmp to duplicate/multiplicate even segments not bounded by homologous sequences. We also provide evidence for PNAmp in mammalian cells. Therefore, PNAmp provides a prototype method to induce structural variations by manipulating replication fork progression.

Dual quorum-sensing control of purine biosynthesis drives pathogenic fitness of Enterococcus faecalis

Enterococcus faecalis is a resident of the human gut, though upon translocation to the blood or body tissues, it can be pathogenic. Here we discover and characterize two peptide-based quorum-sensing systems that transcriptionally modulate de novo purine biosynthesis in E. faecalis. Using a comparative genomic analysis, we find that most enterococcal species do not encode this system; E. moraviensis, E. haemoperoxidus and E. caccae, three species that are closely related to E. faecalis, encode one of the two systems, and only E. faecalis encodes both systems. We show that these systems are important for the intracellular survival of E. faecalis within macrophages and for the fitness of E. faecalis in a murine wound infection model. Taken together, we combine comparative genomics, microbiological, bacterial genetics, transcriptomics, targeted proteomics and animal model experiments to describe a paired quorum sensing mechanism that directly influences central metabolism and impacts the pathogenicity of E. faecalis.