Dynamics of karyotype evolution.

In the evolution of species, the karyotype changes with a timescale of tens to hundreds of thousand years. In the development of cancer, the karyotype often is modified in cancerous cells over the lifetime of an individual. Characterizing these changes and understanding the mechanisms leading to them has been of interest in a broad range of disciplines including evolution, cytogenetics, and cancer genetics. A central issue relates to the relative roles of random vs deterministic mechanisms in shaping the changes. Although it is possible that all changes result from random events followed by selection, many results point to other non-random factors that play a role in karyotype evolution. In cancer, chromosomal instability leads to characteristic changes in the karyotype, in which different individuals with a specific type of cancer display similar changes in karyotype structure over time. Statistical analyses of chromosome lengths in different species indicate that the length distribution of chromosomes is not consistent with models in which the lengths of chromosomes are random or evolve solely by simple random processes. A better understanding of the mechanisms underlying karyotype evolution should enable the development of quantitative theoretical models that combine the random and deterministic processes that can be compared to experimental determinations of the karyotype in diverse settings.

Genomic analysis and phylogenetic characterization of Himalayan snow trout, Schizothorax esocinus based on mitochondrial protein-coding genes.

BACKGROUND: Mitochondrial DNA (mtDNA) has become a significant tool for exploring genetic diversity and delineating evolutionary links across diverse taxa. Within the group of cold-water fish species that are native to the Indian Himalayan region, Schizothorax esocinus holds particular importance due to its ecological significance and is potentially vulnerable to environmental changes. This research aims to clarify the phylogenetic relationships within the Schizothorax genus by utilizing mitochondrial protein-coding genes. METHODS: Standard protocols were followed for the isolation of DNA from S. esocinus. For the amplification of mtDNA, overlapping primers were used, and then subsequent sequencing was performed. The genetic features were investigated by the application of bioinformatic approaches. These approaches covered the evaluation of nucleotide composition, codon usage, selective pressure using nonsynonymous substitution /synonymous substitution (Ka/Ks) ratios, and phylogenetic analysis. RESULTS: The study specifically examined the 13 protein-coding genes of Schizothorax species which belongs to the Schizothoracinae subfamily. Nucleotide composition analysis showed a bias towards A + T content, consistent with other cyprinid fish species, suggesting evolutionary conservation. Relative Synonymous Codon Usage highlighted leucine as the most frequent (5.18%) and cysteine as the least frequent (0.78%) codon. The positive AT-skew and the predominantly negative GC-skew indicated the abundance of A and C. Comparative analysis revealed significant conservation of amino acids in multiple genes. The majority of amino acids were hydrophobic rather than polar. The purifying selection was revealed by the genetic distance and Ka/Ks ratios. Phylogenetic study revealed a significant genetic divergence between S. esocinus and other Schizothorax species with interspecific K2P distances ranging from 0.00 to 8.87%, with an average of 5.76%. CONCLUSION: The present study provides significant contributions to the understanding of mitochondrial genome diversity and genetic evolution mechanisms in Schizothoracinae, hence offering vital insights for the development of conservation initiatives aimed at protecting freshwater fish species.

Pervasive relaxed selection in termite genomes.

Genetic changes that enabled the evolution of eusociality have long captivated biologists. More recently, attention has focussed on the consequences of eusociality on genome evolution. Studies have reported higher molecular evolutionary rates in eusocial hymenopteran insects compared with their solitary relatives. To investigate the genomic consequences of eusociality in termites, we analysed nine genomes, including newly sequenced genomes from three non-eusocial cockroaches. Using a phylogenomic approach, we found that termite genomes have experienced lower rates of synonymous substitutions than those of cockroaches, possibly as a result of longer generation times. We identified higher rates of non-synonymous substitutions in termite genomes than in cockroach genomes, and identified pervasive relaxed selection in the former (24-31% of the genes analysed) compared with the latter (2-4%). We infer that this is due to reductions in effective population size, rather than gene-specific effects (e.g. indirect selection of caste-biased genes). We found no obvious signature of increased genetic load in termites, and postulate efficient purging of deleterious alleles at the colony level. Additionally, we identified genomic adaptations that may underpin caste differentiation, such as genes involved in post-translational modifications. Our results provide insights into the evolution of termites and the genomic consequences of eusociality more broadly.

The molecular evolution of cancer associated genes in mammals.

Cancer is a disease that many multicellular organisms have faced for millions of years, and species have evolved various tumour suppression mechanisms to control oncogenesis. Although cancer occurs across the tree of life, cancer related mortality risks vary across mammalian orders, with Carnivorans particularly affected. Evolutionary theory predicts different selection pressures on genes associated with cancer progression and suppression, including oncogenes, tumour suppressor genes and immune genes. Therefore, we investigated the evolutionary history of cancer associated gene sequences across 384 mammalian taxa, to detect signatures of selection across categories of oncogenes (GRB2, FGL2 and CDC42), tumour suppressors (LITAF, Casp8 and BRCA2) and immune genes (IL2, CD274 and B2M). This approach allowed us to conduct a fine scale analysis of gene wide and site-specific signatures of selection across mammalian lineages under the lens of cancer susceptibility. Phylogenetic analyses revealed that for most species the evolution of cancer associated genes follows the species' evolution. The gene wide selection analyses revealed oncogenes being the most conserved, tumour suppressor and immune genes having similar amounts of episodic diversifying selection. Despite BRCA2's status as a key caretaker gene, episodic diversifying selection was detected across mammals. The site-specific selection analyses revealed that the two apoptosis associated domains of the Casp8 gene of bats (Chiroptera) are under opposing forces of selection (positive and negative respectively), highlighting the importance of site-specific selection analyses to understand the evolution of highly complex gene families. Our results highlighted the need to critically assess different types of selection pressure on cancer associated genes when investigating evolutionary adaptations to cancer across the tree of life. This study provides an extensive assessment of cancer associated genes in mammals with highly representative, and substantially large sample size for a comparative genomic analysis in the field and identifies various avenues for future research into the mechanisms of cancer resistance and susceptibility in mammals.

Polyphyletic origin of saxitoxin biosynthesis genes in the marine dinoflagellate Alexandrium revealed by comparative transcriptomics.

The marine dinoflagellate Alexandrium is known to form harmful algal blooms, and at least 14 species within the genus can produce saxitoxins (STXs). STX biosynthesis genes (sxt) are individually revealed in toxic dinoflagellates; however, the evolutionary history remains controversial. Herein, we determined the transcriptome sequences of toxic Alexandrium (A. catenella and A. pacificum) and non-toxic Alexandrium (A. fraterculus and A. fragae) and characterized their sxt by focusing on evolutionary events and STX production. Comparative transcriptome analysis revealed higher homology of the sxt in toxic Alexandrium than in non-toxic species. Notably, non-toxic Alexandrium spp. were found to have lost two sxt core genes, namely sxtA4 and sxtG. Expression levels of 28 transcripts related to eight sxt core genes showed that sxtA, sxtG, and sxtI were relatively high (>1.5) in the toxic group compared to the non-toxic group. In contrast, the non-toxic group showed high expression levels in sxtU (1.9) and sxtD (1.7). Phylogenetic tree comparisons revealed distinct evolutionary patterns between 28S rDNA and sxtA, sxtB, sxtI, sxtD, and sxtU. However, similar topology was observed between 28S rDNA, sxtS, and sxtH/T. In the sxtB and sxtI phylogeny trees, toxic Alexandrium and cyanobacteria were clustered together, separating from non-toxic species. These suggest that Alexandrium may acquire sxt genes independently via horizontal gene transfer from toxic cyanobacteria and other multiple sources, demonstrating monocistronic transcripts of sxt in dinoflagellates.

Lung adenocarcinomas without driver genes converge to common adaptive strategies through diverse genetic, epigenetic, and niche construction evolutionary pathways.

Somatic evolution selects cancer cell phenotypes that maximize survival and proliferation in dynamic environments. Although cancer cells are molecularly heterogeneous, we hypothesized convergent adaptive strategies to common host selection forces can be inferred from patterns of epigenetic and genetic evolutionary selection in similar tumors. We systematically investigated gene mutations and expression changes in lung adenocarcinomas with no common driver genes (n = 313). Although 13,461 genes were mutated in at least one sample, only 376 non-synonymous mutations evidenced positive evolutionary selection with conservation of 224 genes, while 1736 and 2430 genes exhibited ≥ two-fold increased and ≥ 50% decreased expression, respectively. Mutations under positive selection are more frequent in genes with significantly altered expression suggesting they often "hardwire" pre-existing epigenetically driven adaptations. Conserved genes averaged 16-fold higher expression in normal lung tissue compared to those with selected mutations demonstrating pathways necessary for both normal cell function and optimal cancer cell fitness. The convergent LUAD phenotype exhibits loss of differentiated functions and cell-cell interactions governing tissue organization. Conservation with increased expression is found in genes associated with cell cycle, DNA repair, p53 pathway, epigenetic modifiers, and glucose metabolism. No canonical driver gene pathways exhibit strong positive selection, but extensive down-regulation of membrane ion channels suggests decreased transmembrane potential may generate persistent proliferative signals. NCD LUADs perform niche construction generating a stiff, immunosuppressive microenvironment through selection of specific collagens and proteases. NCD LUADs evolve to a convergent phenotype through a network of interconnected genetic, epigenetic, and ecological pathways.

Mutation characteristics and molecular evolution of ovarian metastasis from gastric cancer and potential biomarkers for paclitaxel treatment.

Ovarian metastasis is one of the major causes of treatment failure in patients with gastric cancer (GC). However, the genomic characteristics of ovarian metastasis in GC remain poorly understood. In this study, we enroll 74 GC patients with ovarian metastasis, with 64 having matched primary and metastatic samples. Here, we show a characterization of the mutation landscape of this disease, alongside an investigation into the molecular heterogeneity and pathway mutation enrichments between synchronous and metachronous metastasis. We classify patients into distinct clonal evolution patterns based on the distribution of mutations in paired samples. Notably, the parallel evolution group exhibits the most favorable prognosis. Additionally, by analyzing the differential response to chemotherapy, we identify potential biomarkers, including SALL4, CCDC105, and CLDN18, for predicting the efficacy of paclitaxel treatment. Furthermore, we validate that CLDN18 fusion mutations improve tumor response to paclitaxel treatment in GC with ovarian metastasis in vitro and vivo.

Interspecific transfer of genetic information through polyploid bridges.

Hybridization blurs species boundaries and leads to intertwined lineages resulting in reticulate evolution. Polyploidy, the outcome of whole genome duplication (WGD), has more recently been implicated in promoting and facilitating hybridization between polyploid species, potentially leading to adaptive introgression. However, because polyploid lineages are usually ephemeral states in the evolutionary history of life it is unclear whether WGD-potentiated hybridization has any appreciable effect on their diploid counterparts. Here, we develop a model of cytotype dynamics within mixed-ploidy populations to demonstrate that polyploidy can in fact serve as a bridge for gene flow between diploid lineages, where introgression is fully or partially hampered by the species barrier. Polyploid bridges emerge in the presence of triploid organisms, which despite critically low levels of fitness, can still allow the transfer of alleles between diploid states of independently evolving mixed-ploidy species. Notably, while marked genetic divergence prevents polyploid-mediated interspecific gene flow, we show that increased recombination rates can offset these evolutionary constraints, allowing a more efficient sorting of alleles at higher-ploidy levels before introgression into diploid gene pools. Additionally, we derive an analytical approximation for the rate of gene flow at the tetraploid level necessary to supersede introgression between diploids with nonzero introgression rates, which is especially relevant for plant species complexes, where interspecific gene flow is ubiquitous. Altogether, our results illustrate the potential impact of polyploid bridges on the (re)distribution of genetic material across ecological communities during evolution, representing a potential force behind reticulation.

Lineage-specific gene duplication and expansion of DUF1216 gene family in Brassicaceae.

Proteins containing domain of unknown function (DUF) are prevalent in eukaryotic genome. The DUF1216 proteins possess a conserved DUF1216 domain resembling to the mediator protein of Arabidopsis RNA polymerase II transcriptional subunit-like protein. The DUF1216 family are specifically existed in Brassicaceae, however, no comprehensive evolutionary analysis of DUF1216 genes have been performed. We performed a first comprehensive genome-wide analysis of DUF1216 proteins in Brassicaceae. Totally 284 DUF1216 genes were identified in 27 Brassicaceae species and classified into four subfamilies on the basis of phylogenetic analysis. The analysis of gene structure and conserved motifs revealed that DUF1216 genes within the same subfamily exhibited similar intron/exon patterns and motif composition. The majority members of DUF1216 genes contain a signal peptide in the N-terminal, and the ninth position of the signal peptide in most DUF1216 is cysteine. Synteny analysis revealed that segmental duplication is a major mechanism for expanding of DUF1216 genes in Brassica oleracea, Brassica juncea, Brassica napus, Lepidium meyneii, and Brassica carinata, while in Arabidopsis thaliana and Capsella rubella, tandem duplication plays a major role in the expansion of the DUF1216 gene family. The analysis of Ka/Ks (non-synonymous substitution rate/synonymous substitution rate) ratios for DUF1216 paralogous indicated that most of gene pairs underwent purifying selection. DUF1216 genes displayed a specifically high expression in reproductive tissues in most Brassicaceae species, while its expression in Brassica juncea was specifically high in root. Our studies offered new insights into the phylogenetic relationships, gene structures and expressional patterns of DUF1216 members in Brassicaceae, which provides a foundation for future functional analysis.

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

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

Immunopathological investigation and genetic evolution of Avian leukosis virus Subgroup-J associated with myelocytomatosis in broiler flocks in Egypt.

BACKGROUND: Avian leukosis virus Subgroup-J (ALV-J) is a rapidly oncogenic evolving retrovirus infecting a variety of avian species; causing severe economic losses to the local poultry industry. METHODS: To investigate ALV-J, a total of 117 blood samples and 57 tissue specimens of different organs were collected for virological, and pathological identification, serological examinations, molecular characterization, and sequencing analysis. To the best of our knowledge, this is the first detailed report recorded in broiler flocks in Egypt. The present study targets the prevalence of a viral tumor disease circulating in broiler flocks in the El-Sharqia, El-Dakahliya, and Al-Qalyubiyya Egyptian governorates from 2021 to 2023 using different diagnostic techniques besides ALV-J gp85 genetic diversity determination. RESULT: We first isolated ALV-J on chicken embryo rough cell culture; showing aggregation, rounding, and degeneration. Concerning egg inoculation, embryonic death, stunting, and curling were observed. Only 79 serum samples were positive for ALV-J (67.52%) based on the ELISA test. Histopathological investigation showed tumors consist of uniform masses, usually well-differentiated myelocytes, lymphoid cells, or both in the liver, spleen, and kidneys. Immunohistochemical examination showed that the myelocytomatosis-positive signals were in the spleen, liver, and kidney. The PCR assay of ALV-J gp85 confirmed 545 base pairs with only 43 positive samples (75.4%). Two positive samples were sequenced and submitted to the Genbank with accession numbers (OR509852-OR509853). Phylogenetic analysis based on the gp85 gene showed that the ALV-J Dakahlia-2 isolate is genetically related to ALV-EGY/YA 2021.3, ALV-EGY/YA 2021.4, ALV-EGY/YA 2021.14, and ALV-EGY/YA 2021.9 with amino acid identity percentage 96%, 97%; 96%, 96%; respectively. Furthermore, ALV-J Sharqia-1 isolate is highly genetically correlated to ALV-EGY/YA 2021.14, and ALV-EGY/YA 2021.9, ALV-J isolate QL1, ALV-J isolate QL4, ALV-J isolate QL3, ALV-EGY/YA 2021.4 with amino acid identity percentage 97%, 97%; 98%, 97%, 97%, 95%; respectively. CONCLUSIONS: This study confirmed that ALV-J infection had still been prevalent in broilers in Egypt, and the genetic characteristics of the isolates are diverse.

A phylogenetic method linking nucleotide substitution rates to rates of continuous trait evolution.

Genomes contain conserved non-coding sequences that perform important biological functions, such as gene regulation. We present a phylogenetic method, PhyloAcc-C, that associates nucleotide substitution rates with changes in a continuous trait of interest. The method takes as input a multiple sequence alignment of conserved elements, continuous trait data observed in extant species, and a background phylogeny and substitution process. Gibbs sampling is used to assign rate categories (background, conserved, accelerated) to lineages and explore whether the assigned rate categories are associated with increases or decreases in the rate of trait evolution. We test our method using simulations and then illustrate its application using mammalian body size and lifespan data previously analyzed with respect to protein coding genes. Like other studies, we find processes such as tumor suppression, telomere maintenance, and p53 regulation to be related to changes in longevity and body size. In addition, we also find that skeletal genes, and developmental processes, such as sprouting angiogenesis, are relevant.

wgd v2: a suite of tools to uncover and date ancient polyploidy and whole-genome duplication.

MOTIVATION: Major improvements in sequencing technologies and genome sequence assembly have led to a huge increase in the number of available genome sequences. In turn, these genome sequences form an invaluable source for evolutionary, ecological, and comparative studies. One kind of analysis that has become routine is the search for traces of ancient polyploidy, particularly for plant genomes, where whole-genome duplication (WGD) is rampant. RESULTS: Here, we present a major update of a previously developed tool wgd, namely wgd v2, to look for remnants of ancient polyploidy, or WGD. We implemented novel and improved previously developed tools to (a) construct KS age distributions for the whole-paranome (collection of all duplicated genes in a genome), (b) unravel intragenomic and intergenomic collinearity resulting from WGDs, (c) fit mixture models to age distributions of gene duplicates, (d) correct substitution rate variation for phylogenetic placement of WGDs, and (e) date ancient WGDs via phylogenetic dating of WGD-retained gene duplicates. The applicability and feasibility of wgd v2 for the identification and the relative and absolute dating of ancient WGDs is demonstrated using different plant genomes. AVAILABILITY AND IMPLEMENTATION: wgd v2 is open source and available at https://github.com/heche-psb/wgd.

Molecular evolution and characterization of type III polyketide synthase gene family in Aquilaria sinensis.

2-(2-Phenylethyl) chromone (PEC) and its derivatives are markers of agarwood formation and are also related to agarwood quality. However, the biosynthetic and regulatory mechanisms of PECs still remain mysterious. Several studies suggested that type III polyketide synthases (PKSs) contribute to PEC biosynthesis in Aquilaria sinensis. Furthermore, systematic studies on the evolution of PKSs in A. sinensis have rarely been reported. Herein, we comprehensively analyzed PKS genes from 12 plant genomes and characterized the AsPKSs in detail. A unique branch contained only AsPKS members was identified through evolutionary analysis, including AsPKS01 that was previously indicated to participate in PEC biosynthesis. AsPKS07 and AsPKS08, two tandem-duplicated genes of AsPKS01 and lacking orthologous genes in evolutionary models, were selected for their transient expression in the leaves of Nicotiana benthamiana. Subsequently, PECs were detected in the extracts of N. benthamiana leaves, suggesting that AsPKS07 and AsPKS08 promote PEC biosynthesis. The interaction between the promoters of AsPKS07, AsPKS08 and five basic leucine zippers (bZIPs) from the S subfamily indicated that their transcripts could be regulated by these transcription factors (TFs) and might further contribute to PECs biosynthesis in A. sinensis. Our findings provide valuable insights into the molecular evolution of the PKS gene family in A. sinensis and serve as a foundation for advancing PEC production through the bioengineering of gene clusters. Ultimately, this contribution is expected to shed light on the mechanism underlying agarwood formation.

Shared single copy genes are generally reliable for inferring phylogenetic relationships among polyploid taxa.

Polyploidy, or whole-genome duplication, is expected to confound the inference of species trees with phylogenetic methods for two reasons. First, the presence of retained duplicated genes requires the reconciliation of the inferred gene trees to a proposed species tree. Second, even if the analyses are restricted to shared single copy genes, the occurrence of reciprocal gene loss, where the surviving genes in different species are paralogs from the polyploidy rather than orthologs, will mean that such genes will not have evolved under the corresponding species tree and may not produce gene trees that allow inference of that species tree. Here we analyze three different ancient polyploidy events, using synteny-based inferences of orthology and paralogy to infer gene trees from nearly 17,000 sets of homologous genes. We find that the simple use of single copy genes from polyploid organisms provides reasonably robust phylogenetic signals, despite the presence of reciprocal gene losses. Such gene trees are also most often in accord with the inferred species relationships inferred from maximum likelihood models of gene loss after polyploidy: a completely distinct phylogenetic signal present in these genomes. As seen in other studies, however, we find that methods for inferring phylogenetic confidence yield high support values even in cases where the underlying data suggest meaningful conflict in the phylogenetic signals.

Essentiality, protein-protein interactions and evolutionary properties are key predictors for identifying cancer-associated genes using machine learning.

The distinctive nature of cancer as a disease prompts an exploration of the special characteristics the genes implicated in cancer exhibit. The identification of cancer-associated genes and their characteristics is crucial to further our understanding of this disease and enhanced likelihood of therapeutic drug targets success. However, the rate at which cancer genes are being identified experimentally is slow. Applying predictive analysis techniques, through the building of accurate machine learning models, is potentially a useful approach in enhancing the identification rate of these genes and their characteristics. Here, we investigated gene essentiality scores and found that they tend to be higher for cancer-associated genes compared to other protein-coding human genes. We built a dataset of extended gene properties linked to essentiality and used it to train a machine-learning model; this model reached 89% accuracy and > 0.85 for the Area Under Curve (AUC). The model showed that essentiality, evolutionary-related properties, and properties arising from protein-protein interaction networks are particularly effective in predicting cancer-associated genes. We were able to use the model to identify potential candidate genes that have not been previously linked to cancer. Prioritising genes that score highly by our methods could aid scientists in their cancer genes research.

Emergence of enhancers at late DNA replicating regions.

Enhancers are fast-evolving genomic sequences that control spatiotemporal gene expression patterns. By examining enhancer turnover across mammalian species and in multiple tissue types, we uncover a relationship between the emergence of enhancers and genome organization as a function of germline DNA replication time. While enhancers are most abundant in euchromatic regions, enhancers emerge almost twice as often in late compared to early germline replicating regions, independent of transposable elements. Using a deep learning sequence model, we demonstrate that new enhancers are enriched for mutations that alter transcription factor (TF) binding. Recently evolved enhancers appear to be mostly neutrally evolving and enriched in eQTLs. They also show more tissue specificity than conserved enhancers, and the TFs that bind to these elements, as inferred by binding sequences, also show increased tissue-specific gene expression. We find a similar relationship with DNA replication time in cancer, suggesting that these observations may be time-invariant principles of genome evolution. Our work underscores that genome organization has a profound impact in shaping mammalian gene regulation.

Origin, evolution, and diversification of inositol 1,4,5-trisphosphate 3-kinases in plants and animals.

BACKGROUND: In Eukaryotes, inositol polyphosphates (InsPs) represent a large family of secondary messengers and play crucial roes in various cellular processes. InsPs are synthesized through a series of pohophorylation reactions catalyzed by various InsP kinases in a sequential manner. Inositol 1,4,5-trisphosphate 3-kinase (IP3 3-kinase/IP3K), one member of InsP kinase, plays important regulation roles in InsPs metabolism by specifically phosphorylating inositol 1,4,5-trisphosphate (IP3) to inositol 1,3,4,5-tetrakisphosphate (IP4) in animal cells. IP3Ks were widespread in fungi, plants and animals. However, its evolutionary history and patterns have not been examined systematically. RESULTS: A total of 104 and 31 IP3K orthologues were identified across 57 plant genomes and 13 animal genomes, respectively. Phylogenetic analyses indicate that IP3K originated in the common ancestor before the divergence of fungi, plants and animals. In most plants and animals, IP3K maintained low-copy numbers suggesting functional conservation during plant and animal evolution. In Brassicaceae and vertebrate, IP3K underwent one and two duplication events, respectively, resulting in multiple gene copies. Whole-genome duplication (WGD) was the main mechanism for IP3K duplications, and the IP3K duplicates have experienced functional divergence. Finally, a hypothetical evolutionary model for the IP3K proteins is proposed based on phylogenetic theory. CONCLUSION: Our study reveals the evolutionary history of IP3K proteins and guides the future functions of animal, plant, and fungal IP3K proteins.

Protein engineering a PhotoRNR chimera based on a unifying evolutionary apparatus among the natural classes of ribonucleotide reductases.

Ribonucleotide reductases (RNRs) are essential enzymes that catalyze the de novo transformation of nucleoside 5'-di(tri)phosphates [ND(T)Ps, where N is A, U, C, or G] to their corresponding deoxynucleotides. Despite the diversity of factors required for function and the low sequence conservation across RNRs, a unifying apparatus consolidating RNR activity is explored. We combine aspects of the protein subunit simplicity of class II RNR with a modified version of Escherichia coli class la photoRNRs that initiate radical chemistry with light to engineer a mimic of a class II enzyme. The design of this RNR involves fusing a truncated form of the active site containing α subunit with the functionally important C-terminal tail of the radical-generating ß subunit to render a chimeric RNR. Inspired by a recent cryo-EM structure, a [Re] photooxidant is located adjacent to Y356[ß], which is an essential component of the radical transport pathway in class I RNRs. Combination of this RNR photochimera with cytidine diphosphate (CDP), adenosine triphosphate (ATP), and light resulted in the generation of Y356• along with production of deoxycytidine diphosphate (dCDP) and cytosine. The photoproducts reflect an active site chemistry consistent with both the consensus mechanism of RNR and chemistry observed when RNR is inactivated by mechanism-based inhibitors in the active site. The enzymatic activity of the RNR photochimera in the absence of any ß metallocofactor highlights the adaptability of the 10-stranded αß barrel finger loop to support deoxynucleotide formation and accommodate the design of engineered RNRs.