Therapy-induced secretion of spliceosomal components mediates pro-survival crosstalk between ovarian cancer cells.

Ovarian cancer often develops resistance to conventional therapies, hampering their effectiveness. Here, using ex vivo paired ovarian cancer ascites obtained before and after chemotherapy and in vitro therapy-induced secretomes, we show that molecules secreted by ovarian cancer cells upon therapy promote cisplatin resistance and enhance DNA damage repair in recipient cancer cells. Even a short-term incubation of chemonaive ovarian cancer cells with therapy-induced secretomes induces changes resembling those that are observed in chemoresistant patient-derived tumor cells after long-term therapy. Using integrative omics techniques, we find that both ex vivo and in vitro therapy-induced secretomes are enriched with spliceosomal components, which relocalize from the nucleus to the cytoplasm and subsequently into the extracellular vesicles upon treatment. We demonstrate that these molecules substantially contribute to the phenotypic effects of therapy-induced secretomes. Thus, SNU13 and SYNCRIP spliceosomal proteins promote therapy resistance, while the exogenous U12 and U6atac snRNAs stimulate tumor growth. These findings demonstrate the significance of spliceosomal network perturbation during therapy and further highlight that extracellular signaling might be a key factor contributing to the emergence of ovarian cancer therapy resistance.

Exploring the diversity of cancer-associated fibroblasts: insights into mechanisms of drug resistance.

Introduction: Among the various stromal cell types within the tumor microenvironment, cancer-associated fibroblasts (CAFs) emerge as the predominant constituent, exhibiting a diverse array of oncogenic functions not intrinsic to normal fibroblasts. Their involvement spans across all stages of tumorigenesis, encompassing initiation, progression, and metastasis. Current understanding posits the coexistence of distinct subpopulations of CAFs within the tumor microenvironment across a spectrum of solid tumors, showcasing both pro- and antitumor activities. Recent advancements in single-cell transcriptomics have revolutionized our ability to meticulously dissect the heterogeneity inherent to CAF populations. Furthermore, accumulating evidence underscores the pivotal role of CAFs in conferring therapeutic resistance to tumors against various drug modalities. Consequently, efforts are underway to develop pharmacological agents specifically targeting CAFs. Methods: This review embarks on a comprehensive analysis, consolidating data from 36 independent single-cell RNA sequencing investigations spanning 17 distinct human malignant tumor types. Results: Our exploration centers on elucidating CAF population markers, discerning their prognostic relevance, delineating their functional contributions, and elucidating the underlying mechanisms orchestrating chemoresistance. Discussion: Finally, we deliberate on the therapeutic potential of harnessing CAFs as promising targets for intervention strategies in clinical oncology.

Splicing Factor PQBP1 Curtails BAX Expression to Promote Ovarian Cancer Progression.

Splicing factor polyglutamine binding protein-1 (PQBP1) is abundantly expressed in the central nervous system during development, and mutations in the gene cause intellectual disability. However, the roles of PQBP1 in cancer progression remain largely unknown. Here, it is shown that PQBP1 overexpression promotes tumor progression and indicates worse prognosis in ovarian cancer. Integrative analysis of spyCLIP-seq and RNA-seq data reveals that PQBP1 preferentially binds to exon regions and modulates exon skipping. Mechanistically, it is shown that PQBP1 regulates the splicing of genes related to the apoptotic signaling pathway, including BAX. PQBP1 promotes BAX exon 2 skipping to generate a truncated isoform that undergoes degradation by nonsense-mediated mRNA decay, thus making cancer cells resistant to apoptosis. In contrast, PQBP1 depletion or splice-switching antisense oligonucleotides promote exon 2 inclusion and thus increase BAX expression, leading to inhibition of tumor growth. Together, the results demonstrate an oncogenic role of PQBP1 in ovarian cancer and suggest that targeting the aberrant splicing mediated by PQBP1 has therapeutic potential in cancer treatment.

The Nitro Group Reshapes the Effects of Pyrido[3,4-g]quinazoline Derivatives on DYRK/CLK Activity and RNA Splicing in Glioblastoma Cells.

Serine-threonine protein kinases of the DYRK and CLK families regulate a variety of vital cellular functions. In particular, these enzymes phosphorylate proteins involved in pre-mRNA splicing. Targeting splicing with pharmacological DYRK/CLK inhibitors emerged as a promising anticancer strategy. Investigation of the pyrido[3,4-g]quinazoline scaffold led to the discovery of DYRK/CLK binders with differential potency against individual enzyme isoforms. Exploring the structure-activity relationship within this chemotype, we demonstrated that two structurally close compounds, pyrido[3,4-g]quinazoline-2,10-diamine 1 and 10-nitro pyrido[3,4-g]quinazoline-2-amine 2, differentially inhibited DYRK1-4 and CLK1-3 protein kinases in vitro. Unlike compound 1, compound 2 efficiently inhibited DYRK3 and CLK4 isoenzymes at nanomolar concentrations. Quantum chemical calculations, docking and molecular dynamic simulations of complexes of 1 and 2 with DYRK3 and CLK4 identified a dramatic difference in electron donor-acceptor properties critical for preferential interaction of 2 with these targets. Subsequent transcriptome and proteome analyses of patient-derived glioblastoma (GBM) neurospheres treated with 2 revealed that this compound impaired CLK4 interactions with spliceosomal proteins, thereby altering RNA splicing. Importantly, 2 affected the genes that perform critical functions for cancer cells including DNA damage response, p53 signaling and transcription. Altogether, these results provide a mechanistic basis for the therapeutic efficacy of 2 previously demonstrated in in vivo GBM models.

Four sides to the story: A proteomic comparison of liquid-phase and matrix-bound extracellular vesicles in 2D and 3D cell cultures.

Multipotent mesenchymal stromal cells (MSCs)-derived extracellular vesicles (EVs) play important roles in cellular communication and are extensively studied as promising therapeutic agents. While there is a substantial pool of studies on liquid-phase EVs, data on EVs bound to the extracellular matrix (ECM) is lacking. There is also an emerging trend of accumulating and comparing data on characteristics of EVs obtained in different culturing conditions. Aiming to reveal proteomic signatures of EVs obtained from conditioned media and ECM of MSCs cultured in 2D and 3D conditions, we performed liquid chromatography with tandem mass spectrometry. Bioinformatic analysis revealed common patterns in proteomic composition of liquid-phase EVs and matrix-bound vesicles (MBVs), namely extracellular environment organization, immune, and transport pathways enrichment. However, extracellular environmental organization pathways are more enriched in liquid-phase EVs than in MBVs, while MBVs proteins noticeably enrich enzymatic pathways. Furthermore, each type of EVs from 2D and 3D cultures has a unique differential abundance profile. We have also performed comparative functional assays, namely scratch assay to assess EVs effect on cell migration and tubulogenesis assay to evaluate EVs angiogenic potential. We found that both liquid-phase EVs and MBVs enhance cell migration, while angiogenic potential is higher in MBVs. Results of the present study suggest that while both liquid-phase EVs and MBVs have therapeutic potential, some unique features of each subgroup may determine optimal areas of their application.

Expression level of CD117 (KIT) on ovarian cancer extracellular vesicles correlates with tumor aggressiveness.

Ovarian cancer is known to be the most lethal malignancy among all gynecological cancers affecting a large number of women worldwide. The treatment of ovarian cancer is challenging due to the high recurrence rate of the disease and is further complicated by acquired chemoresistance. Most ovarian cancer deaths are the result of the metastatic spread of drug-resistant cells. The theory of cancer stem cells (CSC) suggests that both tumor initiation and progression are driven by a population of undifferentiated capable of self-renewal, tumor initiation and development of chemoresistance. The CD117 mast/stem cell growth factor receptor (KIT) is the most commonly used marker for ovarian CSCs. Here, we analyze the correlation between CD117 expression and histological tumor type in ovarian cancer cell lines (SK-OV-3 and MES-OV) and in small/medium extracellular vesicles (EVs) isolated from the urine of ovarian cancer patients. We have demonstrated that the abundance of CD117 on cells and EVs is correlated with tumor grade and therapy resistance status. Moreover, using small EVs isolated from ovarian cancer ascites, it was shown that recurrent disease is characterized by a much higher abundance of CD117 on EVs than primary tumor.

Non-canonical functions of spliceosome components in cancer progression.

Dysregulation of pre-mRNA splicing is a common hallmark of cancer cells and it is associated with altered expression, localization, and mutations of the components of the splicing machinery. In the last few years, it has been elucidated that spliceosome components can also influence cellular processes in a splicing-independent manner. Here, we analyze open source data to understand the effect of the knockdown of splicing factors in human cells on the expression and splicing of genes relevant to cell proliferation, migration, cell cycle regulation, DNA repair, and cell death. We supplement this information with a comprehensive literature review of non-canonical functions of splicing factors linked to cancer progression. We also specifically discuss the involvement of splicing factors in intercellular communication and known autoregulatory mechanisms in restoring their levels in cells. Finally, we discuss strategies to target components of the spliceosome machinery that are promising for anticancer therapy. Altogether, this review greatly expands understanding of the role of spliceosome proteins in cancer progression.

Ratiometric i-Motif-Based Sensor for Precise Long-Term Monitoring of pH Micro Alterations in the Nucleoplasm and Interchromatin Granules.

DNA-intercalated motifs (iMs) are facile scaffolds for the design of various pH-responsive nanomachines, including biocompatible pH sensors. First, DNA pH sensors relied on complex intermolecular scaffolds. Here, we used a simple unimolecular dual-labeled iM scaffold and minimized it by replacing the redundant loop nucleosides with abasic or alkyl linkers. These modifications improved the thermal stability of the iM and increased the rates of its pH-induced conformational transitions. The best effects were obtained upon the replacement of all three native loops with short and flexible linkers, such as the propyl one. The resulting sensor showed a pH transition value equal to 6.9 ± 0.1 and responded rapidly to minor acidification (tau1/2 <1 s for 7.2 → 6.6 pH jump). We demonstrated the applicability of this sensor for pH measurements in the nuclei of human lung adenocarcinoma cells (pH = 7.4 ± 0.2) and immortalized embryonic kidney cells (pH = 7.0 ± 0.2). The sensor stained diffusely the nucleoplasm and piled up in interchromatin granules. These findings highlight the prospects of iMs in the studies of normal and pathological pH-dependent processes in the nucleus, including the formation of biomolecular condensates.

α-Diimine Cisplatin Derivatives: Synthesis, Structure, Cyclic Voltammetry and Cytotoxicity.

Three new Pt(II) complexes [(dpp-DAD)PtCl2] (I), [(Mes-DAD(Me)2)PtCl2] (II) and [(dpp-DAD(Me)2)PtCl2] (III) were synthesized by the direct reaction of [(CH3CN)2PtCl2] and corresponding redox-active 1,4-diaza-1,3-butadienes (DAD). The compounds were isolated in a single crystal form and their molecular structures were determined by X-ray diffraction. The purity of the complexes and their stability in solution was confirmed by NMR analysis. The Pt(II) ions in all compounds are in a square planar environment. The electrochemical reduction of complexes I-III proceeds in two successive cathodic stages. The first quasi-reversible reduction leads to the relatively stable monoanionic complexes; the second cathodic stage is irreversible. The coordination of 1,4-diaza-1,3-butadienes ligands with PtCl2 increases the reduction potential and the electron acceptor ability of the DAD ligands. The synthesized compounds were tested in relation to an adenocarcinoma of the ovary (SKOV3).

Splicing factor BUD31 promotes ovarian cancer progression through sustaining the expression of anti-apoptotic BCL2L12.

Dysregulated expression of splicing factors has important roles in cancer development and progression. However, it remains a challenge to identify the cancer-specific splicing variants. Here we demonstrate that spliceosome component BUD31 is increased in ovarian cancer, and its higher expression predicts worse prognosis. We characterize the BUD31-binding motif and find that BUD31 preferentially binds exon-intron regions near splicing sites. Further analysis reveals that BUD31 inhibition results in extensive exon skipping and a reduced production of long isoforms containing full coding sequence. In particular, we identify BCL2L12, an anti-apoptotic BCL2 family member, as one of the functional splicing targets of BUD31. BUD31 stimulates the inclusion of exon 3 to generate full-length BCL2L12 and promotes ovarian cancer progression. Knockdown of BUD31 or splice-switching antisense oligonucleotide treatment promotes exon 3 skipping and results in a truncated isoform of BCL2L12 that undergoes nonsense-mediated mRNA decay, and the cells subsequently undergo apoptosis. Our findings reveal BUD31-regulated exon inclusion as a critical factor for ovarian cancer cell survival and cancer progression.

Alternative RNA splicing modulates ribosomal composition and determines the spatial phenotype of glioblastoma cells.

Glioblastoma (GBM) is characterized by exceptionally high intratumoral heterogeneity. However, the molecular mechanisms underlying the origin of different GBM cell populations remain unclear. Here, we found that the compositions of ribosomes of GBM cells in the tumour core and edge differ due to alternative RNA splicing. The acidic pH in the core switches before messenger RNA splicing of the ribosomal gene RPL22L1 towards the RPL22L1b isoform. This allows cells to survive acidosis, increases stemness and correlates with worse patient outcome. Mechanistically, RPL22L1b promotes RNA splicing by interacting with lncMALAT1 in the nucleus and inducing its degradation. Contrarily, in the tumour edge region, RPL22L1a interacts with ribosomes in the cytoplasm and upregulates the translation of multiple messenger RNAs including TP53. We found that the RPL22L1 isoform switch is regulated by SRSF4 and identified a compound that inhibits this process and decreases tumour growth. These findings demonstrate how distinct GBM cell populations arise during tumour growth. Targeting this mechanism may decrease GBM heterogeneity and facilitate therapy.

Deeper insights into transcriptional features of cancer-associated fibroblasts: An integrated meta-analysis of single-cell and bulk RNA-sequencing data.

Cancer-associated fibroblasts (CAFs) have long been known as one of the most important players in tumor initiation and progression. Even so, there is an incomplete understanding of the identification of CAFs among tumor microenvironment cells as the list of CAF marker genes varies greatly in the literature, therefore it is imperative to find a better way to identify reliable markers of CAFs. To this end, we summarized a large number of single-cell RNA-sequencing data of multiple tumor types and corresponding normal tissues. As a result, for 9 different types of cancer, we identified CAF-specific gene expression signatures and found 10 protein markers that showed strongly positive staining of tumor stroma according to the analysis of IHC images from the Human Protein Atlas database. Our results give an insight into selecting the most appropriate combination of cancer-associated fibroblast markers. Furthermore, comparison of different approaches for studying differences between cancer-associated and normal fibroblasts (NFs) illustrates the superiority of transcriptome analysis of fibroblasts obtained from fresh tissue samples. Using single-cell RNA sequencing data, we identified common differences in gene expression patterns between normal and cancer-associated fibroblasts, which do not depend on the type of tumor.

Interplay between A-to-I Editing and Splicing of RNA: A Potential Point of Application for Cancer Therapy.

Adenosine-to-inosine RNA editing is a system of post-transcriptional modification widely distributed in metazoans which is catalyzed by ADAR enzymes and occurs mostly in double-stranded RNA (dsRNA) before splicing. This type of RNA editing changes the genetic code, as inosine generally pairs with cytosine in contrast to adenosine, and this expectably modulates RNA splicing. We review the interconnections between RNA editing and splicing in the context of human cancer. The editing of transcripts may have various effects on splicing, and resultant alternatively spliced isoforms may be either tumor-suppressive or oncogenic. Dysregulated RNA splicing in cancer often causes the release of excess amounts of dsRNA into cytosol, where specific dsRNA sensors provoke antiviral-like responses, including type I interferon signaling. These responses may arrest cell division, causing apoptosis and, externally, stimulate antitumor immunity. Thus, small-molecule spliceosome inhibitors have been shown to facilitate the antiviral-like signaling and are considered to be potential cancer therapies. In turn, a cytoplasmic isoform of ADAR can deaminate dsRNA in cytosol, thereby decreasing its levels and diminishing antitumor innate immunity. We propose that complete or partial inhibition of ADAR may enhance the proapoptotic and cytotoxic effects of splicing inhibitors and that it may be considered a promising addition to cancer therapies targeting RNA splicing.

Highly Sensitive Nanomagnetic Quantification of Extracellular Vesicles by Immunochromatographic Strips: A Tool for Liquid Biopsy.

Extracellular vesicles (EVs) are promising agents for liquid biopsy-a non-invasive approach for the diagnosis of cancer and evaluation of therapy response. However, EV potential is limited by the lack of sufficiently sensitive, time-, and cost-efficient methods for their registration. This research aimed at developing a highly sensitive and easy-to-use immunochromatographic tool based on magnetic nanoparticles for EV quantification. The tool is demonstrated by detection of EVs isolated from cell culture supernatants and various body fluids using characteristic biomarkers, CD9 and CD81, and a tumor-associated marker-epithelial cell adhesion molecules. The detection limit of 3.7 × 105 EV/µL is one to two orders better than the most sensitive traditional lateral flow system and commercial ELISA kits. The detection specificity is ensured by an isotype control line on the test strip. The tool's advantages are due to the spatial quantification of EV-bound magnetic nanolabels within the strip volume by an original electronic technique. The inexpensive tool, promising for liquid biopsy in daily clinical routines, can be extended to other relevant biomarkers.

New Insights into Therapy-Induced Progression of Cancer.

The malignant tumor is a complex heterogeneous set of cells functioning in a no less heterogeneous microenvironment. Like any dynamic system, cancerous tumors evolve and undergo changes in response to external influences, including therapy. Initially, most tumors are susceptible to treatment. However, remaining cancer cells may rapidly reestablish the tumor after a temporary remission. These new populations of malignant cells usually have increased resistance not only to the first-line agent, but also to the second- and third-line drugs, leading to a significant decrease in patient survival. Multiple studies describe the mechanism of acquired therapy resistance. In past decades, it became clear that, in addition to the simple selection of pre-existing resistant clones, therapy induces a highly complicated and tightly regulated molecular response that allows tumors to adapt to current and even subsequent therapeutic interventions. This review summarizes mechanisms of acquired resistance, such as secondary genetic alterations, impaired function of drug transporters, and autophagy. Moreover, we describe less obvious molecular aspects of therapy resistance in cancers, including epithelial-to-mesenchymal transition, cell cycle alterations, and the role of intercellular communication. Understanding these molecular mechanisms will be beneficial in finding novel therapeutic approaches for cancer therapy.

Chromatin Trapping of Factors Involved in DNA Replication and Repair Underlies Heat-Induced Radio- and Chemosensitization.

Hyperthermia has been used as an adjuvant treatment for radio- and chemotherapy for decades. In addition to its effects on perfusion and oxygenation of cancer tissues, hyperthermia can enhance the efficacy of DNA-damaging treatments such as radiotherapy and chemotherapy. Although it is believed that the adjuvant effects are based on hyperthermia-induced dysfunction of DNA repair systems, the mechanisms of these dysfunctions remain elusive. Here, we propose that elevated temperatures can induce chromatin trapping (c-trapping) of essential factors, particularly those involved in DNA repair, and thus enhance the sensitization of cancer cells to DNA-damaging therapeutics. Using mass spectrometry-based proteomics, we identified proteins that could potentially undergo c-trapping in response to hyperthermia. Functional analyses of several identified factors involved in DNA repair demonstrated that c-trapping could indeed be a mechanism of hyperthermia-induced transient deficiency of DNA repair systems. Based on our proteomics data, we showed for the first time that hyperthermia could inhibit maturation of Okazaki fragments and activate a corresponding poly(ADP-ribose) polymerase-dependent DNA damage response. Together, our data suggest that chromatin trapping of factors involved in DNA repair and replication contributes to heat-induced radio- and chemosensitization.

Investigation of Inter- and Intratumoral Heterogeneity of Glioblastoma Using TOF-SIMS.

Glioblastoma (GBM) is one of the most aggressive human cancers with a median survival of less than two years. A distinguishing pathological feature of GBM is a high degree of inter- and intratumoral heterogeneity. Intertumoral heterogeneity of GBM has been extensively investigated on genomic, methylomic, transcriptomic, proteomic and metabolomics levels, however only a few studies describe intratumoral heterogeneity because of the lack of methods allowing to analyze GBM samples with high spatial resolution. Here, we applied TOF-SIMS (Time-of-flight secondary ion mass spectrometry) for the analysis of single cells and clinical samples such as paraffin and frozen tumor sections obtained from 57 patients. We developed a technique that allows us to simultaneously detect the distribution of proteins and metabolites in glioma tissue with 800 nm spatial resolution. Our results demonstrate that according to TOF-SIMS data glioma samples can be subdivided into clinically relevant groups and distinguished from the normal brain tissue. In addition, TOF-SIMS was able to elucidate differences between morphologically distinct regions of GBM within the same tumor. By staining GBM sections with gold-conjugated antibodies against Caveolin-1 we could visualize border between zones of necrotic and cellular tumor and subdivide glioma samples into groups characterized by different survival of the patients. Finally, we demonstrated that GBM contains cells that are characterized by high levels of Caveolin-1 protein and cholesterol. This population may partly represent a glioma stem cells. Collectively, our results show that the technique described here allows to analyze glioma tissues with a spatial resolution beyond reach of most of other omics approaches and the obtained data may be used to predict clinical behavior of the tumor.

Peptidome profiling dataset of ovarian cancer and non-cancer proximal fluids: Ascites and blood sera.

Despite a large number of proteomic studies of biological fluids from ovarian cancer patients, there is a lack of sensitive screening methods in clinical practice (Kim et al., 2016) (DOI:https://doi.org/10.1111/cas.12987[1]). Low molecular weight endogenous peptides more easily diffuse across endothelial barriers than proteins and can be more relevant biomarker candidates (Meo et al., 2016) (DOI:https://doi.org/10.18632/oncotarget.8931[2], (Bery et al., 2014) DOI:https://doi.org/10.1186/1559-0275-11-13[3], (Huang et al., 2018) DOI:https://doi.org/10.1097/IGC.0000000000001166[4]). Detailed peptidomic analysis of 26 ovarian cancer and 15 non-cancer samples of biological fluids (ascites and sera) were performed using TripleTOF 5600+ mass-spectrometer. Prior to LC-MS/MS analysis, peptides were extracted from biological fluids using anion exchange sorbent with subsequent peptide desorption from the surface of highly abundant proteins. In total, we identified 4874 peptides; 3123 peptides were specific for the ovarian cancer samples. The mass-spectrometry peptidomics data presented in this data article have been deposited to the ProteomeXchange Consortium (Deutsch et al., 2017) (DOI:https://doi.org/10.1093/nar/gkw936[5]) via the PRIDE partner repository with the dataset identifier PXD009382 and https://doi.org/10.6019/PXD009382, http://www.ebi.ac.uk/pride/archive/projects/PXD009382.

Nucleotide Modifications Decrease Innate Immune Response Induced by Synthetic Analogs of snRNAs and snoRNAs.

Short nuclear regulatory RNAs play a key role in the main stages of maturation of the precursors of the major RNA species. Small nuclear RNAs (snRNAs) form the core of the spliceosome and are responsible for the splicing of pre-mRNA molecules. Small nucleolar RNAs (snoRNAs) direct post-transcriptional modification of pre-rRNAs. A promising strategy for the development of non-coding RNA (ncRNAs) mimicking molecules is the introduction of modified nucleotides, which are normally present in natural ncRNAs, into the structure of synthetic RNAs. We have created a set of snoRNAs and snRNA analogs and studied the effect of base modifications, specifically, pseudouridine (Ψ) and 5-methylcytidine (m5C), on the immune-stimulating and cytotoxic properties of these RNAs. Here, we performed a whole-transcriptome study of the influence of synthetic snoRNA analogs with various modifications on gene expression in human cells. Moreover, we confirmed the role of PKR in the recognition of snoRNA and snRNA analogs using the short hairpin RNA (shRNA) technique. We believe that the data obtained will contribute to the understanding of the role of nucleotide modification in ncRNA functions, and can be useful for creating the agents for gene regulation based on the structure of natural snoRNAs and snRNAs.

Genetically encodable bioluminescent system from fungi.

Bioluminescence is found across the entire tree of life, conferring a spectacular set of visually oriented functions from attracting mates to scaring off predators. Half a dozen different luciferins, molecules that emit light when enzymatically oxidized, are known. However, just one biochemical pathway for luciferin biosynthesis has been described in full, which is found only in bacteria. Here, we report identification of the fungal luciferase and three other key enzymes that together form the biosynthetic cycle of the fungal luciferin from caffeic acid, a simple and widespread metabolite. Introduction of the identified genes into the genome of the yeast Pichia pastoris along with caffeic acid biosynthesis genes resulted in a strain that is autoluminescent in standard media. We analyzed evolution of the enzymes of the luciferin biosynthesis cycle and found that fungal bioluminescence emerged through a series of events that included two independent gene duplications. The retention of the duplicated enzymes of the luciferin pathway in nonluminescent fungi shows that the gene duplication was followed by functional sequence divergence of enzymes of at least one gene in the biosynthetic pathway and suggests that the evolution of fungal bioluminescence proceeded through several closely related stepping stone nonluminescent biochemical reactions with adaptive roles. The availability of a complete eukaryotic luciferin biosynthesis pathway provides several applications in biomedicine and bioengineering.