Global transcriptome modulation by xenobiotics: the role of alternative splicing in adaptive responses to chemical exposures

BACKGROUND

Xenobiotic exposures can extensively influence the expression and alternative splicing of drug-metabolizing enzymes, including cytochromes P450 (CYPs), though their transcriptome-wide impact on splicing remains underexplored. This study used a well-characterized splicing event in the Cyp2b2 gene to validate a sandwich-cultured primary rat hepatocyte model for studying global splicing in vitro. Using endpoint PCR, RNA sequencing, and bioinformatics tools (rSeqDiff, rMATs, IGV), we analyzed differential gene expression and splicing in CYP and nuclear receptor genes, as well as the entire transcriptome, to understand how xenobiotic exposures shape alternative splicing and activate xenosensors.

METHODS

Primary rat hepatocytes in sandwich culture were exposed to two methylenedioxybenzene (MDB) congeners and carbamazepine, with gene expression and splicing assessed. A 3D-clustergram integrating KEGG pathway analysis with differential gene expression provided distinct splicing landscapes for each xenobiotic, showing that splicing diversity does not always align with gene expression changes.

RESULTS

Endpoint PCR revealed a Cyp2b2v to wild-type Cyp2b2 splicing ratio near 1:1 (100%) under most conditions, while RNA-seq showed a stable baseline closer to 40%. C6-MDB reduced this ratio to ~ 50% by PCR and ~ 39% by RNA-seq, showing slight method-dependent variations yet consistent trends. In contrast, exon 6 skipping in Cyp1a1 occurred only with MDB exposure, implicating AHR activation. Xenobiotic treatments also induced alternative splicing in defensome and stress-responsive genes, including the phase II enzyme Gstm3, Albumin, Orm1, and Fxyd1, highlighting their roles in xenobiotic response modulation. Significant splicing changes in factors such as SRSF1, SRSF7, and METTL3 suggest a coordinated feedback loop involving epitranscriptomic modulation and cross-talk within SR protein networks, refining splice site selection, transcript stability, and cellular fate.

CONCLUSIONS

This study demonstrates how xenobiotic structural features influence gene expression and splicing, revealing splicing patterns that expand our understanding of transcriptome diversity and function. By identifying regulatory mechanisms, including AHR activation, epitranscriptomic modulation, and crosstalk within SR protein networks, that shape adaptive responses to xenobiotic stress, this work offers insights into the splicing and transcriptional networks that maintain cellular homeostasis. These findings provide predictive biomarkers for toxic exposures and highlight the potential of splicing profiles as diagnostic tools for assessing the health impacts of chemical exposure.

Aliens in the CYPome of the black fungus gnat, Bradysia coprophila

The diverse cytochrome P450 enzymes of insects play essential physiological roles and also play important roles in the metabolism of environmental chemicals such as insecticides. We manually curated the complement of P450 (CYP) genes, or CYPome, of the black fungus gnat, Bradysia (Sciara) coprophila (Diptera, Sciaroidea), a species with a variable number of chromosomes. This CYPome carries two types of "alien" P450 genes. The first type of alien P450s was found among the 163 CYP genes of the core genome (autosomes and X). They consist of 28 sequences resulting from horizontal gene transfer, with closest sequences not found in insects, but in other arthropods, often Collembola. These genes are not contaminants, because they are expressed genes with introns, found in synteny with regular dipteran genes, also found in B. odoriphaga and B. hygida. Two such "alien" genes are representatives of CYP clans not otherwise found in insects, a CYP53 sequence related to fungal CYP53 genes, and a CYP19-like sequence similar to some collembolan sequences but of unclear origin. The second type of alien P450s are represented by 99 sequences from germline-restricted chromosomes (GRC). While most are P450 pseudogenes, 33 are apparently intact, with half being more closely related to P450s from Cecidomyiidae than from Sciaridae, thus supporting the hypothesis of a cross-family hybridization origin of the GRC.

An Updated Overview of the Role of CYP450 during Xenobiotic Metabolization in Regulating the Acute Myeloid Leukemia Microenvironment

Oxidative stress is associated with several acute and chronic disorders, including hematological malignancies such as acute myeloid leukemia, the most prevalent acute leukemia in adults. Xenobiotics are usually harmless compounds that may be detrimental, such as pharmaceuticals, environmental pollutants, cosmetics, and even food additives. The storage of xenobiotics can serve as a defense mechanism or a means of bioaccumulation, leading to adverse effects. During the absorption, metabolism, and cellular excretion of xenobiotics, three steps may be distinguished: (i) inflow by transporter enzymes, (ii) phases I and II, and (iii) phase III. Phase I enzymes, such as those in the cytochrome P450 superfamily, catalyze the conversion of xenobiotics into more polar compounds, contributing to an elevated acute myeloid leukemia risk. Furthermore, genetic polymorphism influences the variability and susceptibility of related myeloid neoplasms, infant leukemias associated with mixed-lineage leukemia (MLL) gene rearrangements, and a subset of de novo acute myeloid leukemia. Recent research has shown a sustained interest in determining the regulators of cytochrome P450, family 2, subfamily E, member 1 (CYP2E1) expression and activity as an emerging field that requires further investigation in acute myeloid leukemia evolution. Therefore, this review suggests that CYP2E1 and its mutations can be a therapeutic or diagnostic target in acute myeloid leukemia.

A CYPome-wide study reveals new potential players in the pathogenesis of Parkinson's disease

Genetic and environmental factors lead to the manifestation of Parkinson's disease (PD) but related mechanisms are only rudimentarily understood. Cytochromes P450 (P450s) are involved in the biotransformation of toxic compounds and in many physiological processes and thus predestinated to be involved in PD. However, so far only SNPs (single nucleotide polymorphisms) in CYP2D6 and CYP2E1 have been associated with the susceptibility of PD. Our aim was to evaluate the role of all 57 human P450s and their redox partners for the etiology and pathophysiology of PD and to identify novel potential players which may lead to the identification of new biomarkers and to a causative treatment of PD. The PPMI (Parkinson's Progression Markers Initiative) database was used to extract the gene sequences of all 57 P450s and their three redox partners to analyze the association of SNPs with the occurrence of PD. Applying statistical analyses of the data, corresponding odds ratios (OR) and confidence intervals (CI) were calculated. We identified SNPs significantly over-represented in patients with a genetic predisposition for PD (GPD patients) or in idiopathic PD (IPD patients) compared to HC (healthy controls). Xenobiotic-metabolizing P450s show a significant accumulation of SNPs in PD patients compared with HC supporting the role of toxic compounds in the pathogenesis of PD. Moreover, SNPs with high OR values (>5) in P450s catalyzing the degradation of cholesterol (CYP46A1, CY7B1, CYP39A1) indicate a prominent role of cholesterol metabolism in the brain for PD risk. Finally, P450s participating in the metabolism of eicosanoids show a strong over-representation of SNPs in PD patients underlining the effect of inflammation on the pathogenesis of PD. Also, the redox partners of P450 show SNPs with OR > 5 in PD patients. Taken together, we demonstrate that SNPs in 26 out of 57 P450s are at least 5-fold over-represented in PD patients suggesting these P450s as new potential players in the pathogenesis of PD. For the first time exceptionally high OR values (up to 12.9) were found. This will lead to deeper insight into the origin and development of PD and may be applied to develop novel strategies for a causative treatment of this disease.

The Role of CYP3A in Health and Disease

CYP3A is an enzyme subfamily in the cytochrome P450 (CYP) superfamily and includes isoforms CYP3A4, CYP3A5, CYP3A7, and CYP3A43. CYP3A enzymes are indiscriminate toward substrates and are unique in that these enzymes metabolize both endogenous compounds and diverse xenobiotics (including drugs); almost the only common characteristic of these compounds is lipophilicity and a relatively large molecular weight. CYP3A enzymes are widely expressed in human organs and tissues, and consequences of these enzymes' activities play a major role both in normal regulation of physiological levels of endogenous compounds and in various pathological conditions. This review addresses these aspects of regulation of CYP3A enzymes under physiological conditions and their involvement in the initiation and progression of diseases.

CYP35 family in Caenorhabditis elegans biological processes: fatty acid synthesis, xenobiotic metabolism, and stress responses

With more than 80 cytochrome P450 (CYP) encoding genes found in the nematode Caenorhabditis elegans (C. elegans), the cyp35 genes are one of the important genes involved in many biological processes such as fatty acid synthesis and storage, xenobiotic stress response, dauer and eggshell formation, and xenobiotic metabolism. The C. elegans CYP35 subfamily consisted of A, B, C, and D, which have the closest homolog to human CYP2 family. C. elegans homologs could answer part of the hunt for human disease genes. This review aims to provide an overview of CYP35 in C. elegans and their human homologs, to explore the roles of CYP35 in various C. elegans biological processes, and how the genes of cyp35 upregulation or downregulation are influenced by biological processes, upon exposure to xenobiotics or changes in diet and environment. The C. elegans CYP35 gene expression could be upregulated by heavy metals, pesticides, anti-parasitic and anti-chemotherapeutic agents, polycyclic aromatic hydrocarbons (PAHs), nanoparticles, drugs, and organic chemical compounds. Among the cyp35 genes, cyp-35A2 is involved in most of the C. elegans biological processes regulation. Further venture of cyp35 genes, the closest homolog of CYP2 which is the largest family of human CYPs, may have the power to locate cyps gene targets, discovery of novel therapeutic strategies, and possibly a successful medical regime to combat obesity, cancers, and cyps gene-related diseases.

A genome-wide association study of plasma concentrations of warfarin enantiomers and metabolites in sub-Saharan black-African patients

Diversity in pharmacogenomic studies is poor, especially in relation to the inclusion of black African patients. Lack of funding and difficulties in recruitment, together with the requirement for large sample sizes because of the extensive genetic diversity in Africa, are amongst the factors which have hampered pharmacogenomic studies in Africa. Warfarin is widely used in sub-Saharan Africa, but as in other populations, dosing is highly variable due to genetic and non-genetic factors. In order to identify genetic factors determining warfarin response variability, we have conducted a genome-wide association study (GWAS) of plasma concentrations of warfarin enantiomers/metabolites in sub-Saharan black-Africans. This overcomes the issue of non-adherence and may have greater sensitivity at genome-wide level, to identify pharmacokinetic gene variants than focusing on mean weekly dose, the usual end-point used in previous studies. Participants recruited at 12 outpatient sites in Uganda and South Africa on stable warfarin dose were genotyped using the Illumina Infinium H3Africa Consortium Array v2. Imputation was conducted using the 1,000 Genomes Project phase III reference panel. Warfarin/metabolite plasma concentrations were determined by high-performance liquid chromatography with tandem mass spectrometry. Multivariable linear regression was undertaken, with adjustment made for five non-genetic covariates and ten principal components of genetic ancestry. After quality control procedures, 548 participants and 17,268,054 SNPs were retained. CYP2C9*8, CYP2C9*9, CYP2C9*11, and the CYP2C cluster SNP rs12777823 passed the Bonferroni-adjusted replication significance threshold (p < 3.21E-04) for warfarin/metabolite ratios. In an exploratory GWAS analysis, 373 unique SNPs in 13 genes, including CYP2C9*8, passed the Bonferroni-adjusted genome-wide significance threshold (p < 3.846E-9), with 325 (87%, all located on chromosome 10) SNPs being associated with the S-warfarin/R-warfarin outcome (top SNP rs11188082, CYP2C19 intron variant, p = 1.55E-17). Approximately 69% of these SNPs were in linkage disequilibrium (r² > 0.8) with CYP2C9*8 (n = 216) and rs12777823 (n = 8). Using a pharmacokinetic approach, we have shown that variants other than CYP2C9*2 and CYP2C9*3 are more important in sub-Saharan black-Africans, mainly due to the allele frequencies. In exploratory work, we conducted the first warfarin pharmacokinetics-related GWAS in sub-Saharan Africans and identified novel SNPs that will require external replication and functional characterization before they can be considered for inclusion in warfarin dosing algorithms.

Cytochrome P450 1B1: A Key Regulator of Ocular Iron Homeostasis and Oxidative Stress

Cytochrome P450 (CYP) 1B1 belongs to the superfamily of heme-containing monooxygenases. Unlike other CYP enzymes, which are highly expressed in the liver, CYP1B1 is predominantly found in extrahepatic tissues, such as the brain, and ocular tissues including retina and trabecular meshwork. CYP1B1 metabolizes exogenous chemicals such as polycyclic aromatic hydrocarbons. CYP1B1 also metabolizes endogenous bioactive compounds including estradiol and arachidonic acid. These metabolites impact various cellular and physiological processes during development and pathological processes. We previously showed that CYP1B1 deficiency mitigates ischemia-mediated retinal neovascularization and drives the trabecular meshwork dysgenesis through increased levels of oxidative stress. However, the underlying mechanisms responsible for CYP1B1-deficiency-mediated increased oxidative stress remain largely unresolved. Iron is an essential element and utilized as a cofactor in a variety of enzymes. However, excess iron promotes the production of hydroxyl radicals, lipid peroxidation, increased oxidative stress, and cell damage. The retinal endothelium is recognized as a major component of the blood-retinal barrier, which controls ocular iron levels through the modulation of proteins involved in iron regulation present in retinal endothelial cells, as well as other ocular cell types including trabecular meshwork cells. We previously showed increased levels of reactive oxygen species and lipid peroxidation in the absence of CYP1B1, and in the retinal vasculature and trabecular meshwork, which was reversed by administration of antioxidant N-acetylcysteine. Here, we review the important role CYP1B1 expression and activity play in maintaining retinal redox homeostasis through the modulation of iron levels by retinal endothelial cells. The relationship between CYP1B1 expression and activity and iron levels has not been previously delineated. We review the potential significance of CYP1B1 expression, estrogen metabolism, and hepcidin-ferroportin regulatory axis in the local regulation of ocular iron levels.

Genome-wide identification, evolution analysis of cytochrome P450 monooxygenase multigene family and their expression patterns during the early somatic embryogenesis in Dimocarpus longan Lour

Cytochrome P450 (CYP), a multi-gene superfamily, is involved in a broad range of physiological processes, including hormone responses and secondary metabolism throughout the plant life cycle. Longan (Dimocarpus longan), a subtropical and tropical evergreen fruit tree, its embryonic development is closely related to the yield and quality of fruits. And a large number of secondary metabolites, such as flavonoids and carotenoids, are also produced during the longan somatic embryogenesis (SE). It is important, therefore, to study potential functions of CYPs in longan. However, the knowledge of longan CYPs is still very limited. Here, a total of 327 DlCYPs were identified using the genome-search method, which could be classified into nine clans. The expansion of the DlCYP family was mainly caused by tandem duplication (TD) events. Promoter cis-acting elements analysis elucidated that DlCYPs played important roles in hormonal responses. A total of 246 DlCYPs exhibited six different expression patterns during the early SE based on longan transcriptomic data. Eight DlCYPs underwent alternative splicing (AS) events, and they might produce one to six isoforms. And the AS transcript of DlCYP97C1 might act as an alternative to the full-length transcript in ICpEC and GE stages. Finally, protein-protein interaction (PPI) networks and miRNA target prediction elucidated that DlCYPs might be involved in the phenylpropanoid metabolic pathway and primarily regulated and targeted by miR413. In summary, our results provided valuable inventory for understanding the classification and biological functions of DlCYPs and provided insight into further functional verification of DlCYPs during the longan early SE.

[Splice variants of mRNA of cytochrome P450 genes: analysis by the nanopore sequencing method in human liver tissue and HepG2 cell line]

The analysis of cytochrome P450 transcripts was carried out by the nanopore sequencing in liver tissue samples of three donors and HepG2 line cells. It has been demonstrated that direct mRNA sequencing with a MinION nanopore sequencer (Oxford Nanopore Technologies) allows one to obtained quantitative profiles for transcripts (and their splice variants) of cytochrome P450 superfamily genes encoding isoforms involved in metabolism of the large (~80%) part of drugs. The splice variant profiles substantially differ for donors. The cytochrome P450 gene expression at the transcript level is significantly weaker in cells of the HepG2 line compared with that in the normal liver tissue. This limits the capability of the direct mRNA nanopore sequencing for studying alternative splicing of cytochrome P450 transcripts in HepG2 cells. Both quantitative and qualitative profiles of the cytochrome P450 gene expression at the transcript level are notably differ in human liver tissue and HepG2 cells.

How paediatric drug development and use could benefit from OMICs: a c4c expert group white paper

The safety and efficacy of pharmacotherapy in children, particularly preterms, neonates, and infants, is limited by a paucity of good quality data from prospective clinical drug trials. A specific challenge is the establishment of valid biomarkers. OMICs technologies may support these efforts, by complementary information about targeted and non-targeted molecules through systematic characterization and quantitation of biological samples. OMICs technologies comprise at least genomics, epigenomics, transcriptomics, proteomics, metabolomics, and microbiomics in addition to the patient's phenotype. OMICs technologies are in part hypothesis-generating allowing an in depth understanding of disease pathophysiology and pharmacological mechanisms. Application of OMICs technologies in paediatrics faces major challenges before routine adoption. First, developmental processes need to be considered, including a sub-division into specific age groups as developmental changes clearly impact OMICs data. Second, compared to the adult population, the number of patients is limited as well as type and amount of necessary biomaterial, especially in neonates and preterms. Thus, advanced trial designs and biostatistical methods, non-invasive biomarkers, innovative biobanking concepts including data and samples from healthy children, as well as analytical approaches (e.g. liquid biopsies) should be addressed to overcome these obstacles. The ultimate goal is to link OMICs technologies with innovative analysis tools, like artificial intelligence at an early stage. The use of OMICs data based on a feasible approach will contribute to identify complex phenotypes and subpopulations of patients to improve development of medicines for children with potential economic advantages.

Genome-wide Analysis of Alternative Gene Splicing Associated with Virulence in the Brown Planthopper Nilaparvata lugens (Hemiptera: Delphacidae)

Alternative splicing of protein coding genes plays a profound role in phenotypic variation for many eukaryotic organisms. The development of high-throughput sequencing and bioinformatics algorithms provides the possibility of genome-wide identification of alternative splicing events in eukaryotes. However, for the brown planthopper Nilaparvata lugens, a destructive pest of rice crops, whole-genome distribution of alternative splicing events and the role of alternative splicing in the phenotypic plasticity of virulence have not previously been estimated. Here, we developed an analysis pipeline to identify alternative splicing events in the genome of N. lugens. Differential expression analysis and functional annotation were performed on datasets related to different virulence phenotypes. In total, 27,880 alternative splicing events corresponding to 9,787 multi-exon genes were detected in N. lugens. Among them, specifically expressed alternative splicing transcripts in the virulent Mudgo population were enriched in metabolic process categories, while transcripts in the avirulent TN1 population were enriched in regulator activity categories. In addition, genes encoding odorant receptor, secreted saliva protein and xenobiotic metabolic P450 monooxygenase showed different splicing patterns between Mudgo population and TN1 population. Host change experiment also revealed that an isoform of a P450 gene could be specially induced by the stimulation of resistant rice variety Mudgo. This research pioneered a genome-wide study of alternative gene splicing in the rice brown planthopper. Differences in alternative splicing between virulent and avirulent populations indicated that alternative splicing might play an important role in the formation of virulence phenotypes in N. lugens.

The Central Role of Cytochrome P450 in Xenobiotic Metabolism-A Brief Review on a Fascinating Enzyme Family

Human Cytochrome P450 (CYP) enzymes constitute a superfamily of membrane-bound hemoproteins that are responsible for the metabolism of a wide variety of clinically, physiologically, and toxicologically important compounds. These heme-thiolate monooxygenases play a pivotal role in the detoxification of xenobiotics, participating in the metabolism of many structurally diverge compounds. This short-review is intended to provide a summary on the major roles of CYPs in Phase I xenobiotic metabolism. The manuscript is focused on eight main topics that include the most relevant aspects of past and current CYP research. Initially, (I) a general overview of the main aspects of absorption, distribution, metabolism, and excretion (ADME) of xenobiotics are presented. This is followed by (II) a background overview on major achievements in the past of the CYP research field. (III) Classification and nomenclature of CYPs is briefly reviewed, followed by (IV) a summary description on CYP's location and function in mammals. Subsequently, (V) the physiological relevance of CYP as the cornerstone of Phase I xenobiotic metabolism is highlighted, followed by (VI) reviewing both genetic determinants and (VI) nongenetic factors in CYP function and activity. The last topic of the review (VIII) is focused on the current challenges of the CYP research field.

The significance of CYP11A1 expression in skin physiology and pathology

CYP11A1, a member of the cytochrome P450 family, plays several key roles in the human body. It catalyzes the first and rate-limiting step in steroidogenesis, converting cholesterol to pregnenolone. Aside from the classical steroidogenic tissues such as the adrenals, gonads and placenta, CYP11A1 has also been found in the brain, gastrointestinal tract, immune systems, and finally the skin. CYP11A1 activity in the skin is regulated predominately by StAR protein and hence cholesterol levels in the mitochondria. However, UVB, UVC, CRH, ACTH, cAMP, and cytokines IL-1, IL-6 and TNFα can also regulate its expression and activity. Indeed, CYP11A1 plays several critical roles in the skin through its initiation of local steroidogenesis and specific metabolism of vitamin D, lumisterol, and 7-dehydrocholesterol. Products of these pathways regulate the protective barrier and skin immune functions in a context-dependent fashion through interactions with a number of receptors. Disturbances in CYP11A1 activity can lead to skin pathology.

Antisense oligonucleotide development for the selective modulation of CYP3A5 in renal disease

CYP3A5 is the primary CYP3A subfamily enzyme expressed in the human kidney and its aberrant expression may contribute to a broad spectrum of renal disorders. Pharmacogenetic studies have reported inconsistent linkages between CYP3A5 expression and hypertension, however, most investigators have considered CYP3A5*1 as active and CYP3A5*3 as an inactive allele. Observations of gender specific differences in CYP3A5*3/*3 protein expression suggest additional complexity in gene regulation that may underpin an environmentally responsive role for CYP3A5 in renal function. Reconciliation of the molecular mechanism driving conditional restoration of functional CYP3A5*3 expression from alternatively spliced transcripts, and validation of a morpholino-based approach for selectively suppressing renal CYP3A5 expression, is the focus of this work. Morpholinos targeting a cryptic splice acceptor created by the CYP3A5*3 mutation in intron 3 rescued functional CYP3A5 expression in vitro, and salt-sensitive cellular mechanisms regulating splicing and conditional expression of CYP3A5*3 transcripts are reported. The potential for a G-quadruplex (G4) in intron 3 to mediate restored splicing to exon 4 in CYP3A5*3 transcripts was also investigated. Finally, a proximal tubule microphysiological system (PT-MPS) was used to evaluate the safety profile of morpholinos in proximal tubule epithelial cells, highlighting their potential as a therapeutic platform for the treatment of renal disease.

Characteristics of In2G Variant in Congenital Adrenal Hyperplasia Due to 21-Hydroxylase Deficiency

Substantial research has been performed during the last decades on the clinical and genetic variability of congenital adrenal hyperplasia (CAH) and its most common form, 21-hydroxylase deficiency (21OHD). CAH is one of the most prevalent autosomal recessive diseases in humans, and it can be divided into classic-further subdivided into salt wasting (SW) and simple virilizing (SV)-and non-classic (NC) forms. Pathogenic variants of CYP21A2 gene, encoding the 21-hydroxylase enzyme, have been reported with variable prevalence in different populations. NM_000500.9:c.293-13C/A>G (In2G) variant represents the most common CYP21A2 gene changes related to the classic 21OHD form. However, the phenotype of In2G carriers is variable depending on the variant homozygous/heterozygous status and combination with other CYP21A2 pathogenic variants. In addition, identical genotypes, harboring the homozygous In2G variant, can present with variable phenotypes including the SW and SV or rarely NC form of the disease. Here, we analyze and present the clinical aspects, genotype/phenotype correlations, and other characteristics related to the CYP21A2 In2G variant.

Cytochrome P450 expression and regulation in the brain

The regulation of brain cytochrome P450 enzymes (CYPs) is different compared with respective hepatic enzymes. This may result from anatomical bases and physiological functions of the two organs. The brain is composed of a variety of functional structures built of different interconnected cell types endowed with specific receptors that receive various neuronal signals from other brain regions. Those signals activate transcription factors or alter functioning of enzyme proteins. Moreover, the blood-brain barrier (BBB) does not allow free penetration of all substances from the periphery into the brain. Differences in neurotransmitter signaling, availability to endogenous and exogenous active substances, and levels of transcription factors between neuronal and hepatic cells lead to differentiated expression and susceptibility to the regulation of CYP genes in the brain and liver. Herein, we briefly describe the CYP enzymes of CYP1-3 families, their distribution in the brain, and discuss brain-specific regulation of CYP genes. In parallel, a comparison to liver CYP regulation is presented. CYP enzymes play an essential role in maintaining the levels of bioactive molecules within normal ranges. These enzymes modulate the metabolism of endogenous neurochemicals, such as neurosteroids, dopamine, serotonin, melatonin, anandamide, and exogenous substances, including psychotropics, drugs of abuse, neurotoxins, and carcinogens. The role of these enzymes is not restricted to xenobiotic-induced neurotoxicity, but they are also involved in brain physiology. Therefore, it is crucial to recognize the function and regulation of CYP enzymes in the brain to build a foundation for future medicine and neuroprotection and for personalized treatment of brain diseases.

Diversity and evolution of the P450 family in arthropods

The P450 family (CYP genes) of arthropods encodes diverse enzymes involved in the metabolism of foreign compounds and in essential endocrine or ecophysiological functions. The P450 sequences (CYPome) from 40 arthropod species were manually curated, including 31 complete CYPomes, and a maximum likelihood phylogeny of nearly 3000 sequences is presented. Arthropod CYPomes are assembled from members of six CYP clans of variable size, the CYP2, CYP3, CYP4 and mitochondrial clans, as well as the CYP20 and CYP16 clans that are not found in Neoptera. CYPome sizes vary from two dozen genes in some parasitic species to over 200 in species as diverse as collembolans or ticks. CYPomes are comprised of few CYP families with many genes and many CYP families with few genes, and this distribution is the result of dynamic birth and death processes. Lineage-specific expansions or blooms are found throughout the phylogeny and often result in genomic clusters that appear to form a reservoir of catalytic diversity maintained as heritable units. Among the many P450s with physiological functions, six CYP families are involved in ecdysteroid metabolism. However, five so-called Halloween genes are not universally represented and do not constitute the unique pathway of ecdysteroid biosynthesis. The diversity of arthropod CYPomes has only partially been uncovered to date and many P450s with physiological functions regulating the synthesis and degradation of endogenous signal molecules (including ecdysteroids) and semiochemicals (including pheromones and defense chemicals) remain to be discovered. Sequence diversity of arthropod P450s is extreme, and P450 sequences lacking the universally conserved Cys ligand to the heme have evolved several times. A better understanding of P450 evolution is needed to discern the relative contributions of stochastic processes and adaptive processes in shaping the size and diversity of CYPomes.

Altered expression of cytochrome P450 enzymes involved in metabolism of androgens and vitamin D in the prostate as a risk factor for prostate cancer

Prostate cancer is the most common malignant disease among men. The signaling pathways, regulated by the androgen and vitamin D receptors, play a key role in prostate cancer. The intracellular level of androgens and vitamin D determines not only receptor functionality, but also the efficacy of cellular processes regulated by them (cell proliferation, apoptosis, differentiation etc.). It is known that several androgen-metabolizing P450s (CYP3A4/5/43 and CYP2B6) and P450 enzymes (CYP2R1, CYP27A1, CYP27B1, CYP24A1, CYP3A4, CYP2J2), which are necessary for vitamin D metabolism, are expressed in the prostate. It was shown that alterations in an expression pattern of the certain cytochrome P450s might lead to the development of castration-resistant cancer (CYP3A4, CYP2J2, CYP24A1), and to chemo-resistance (CYP3A4, CYP3A5, CYP2B6) and early mortality (CYP2B6, CYP27A1, CYP24A1). Moreover, steroidogenic CYPs (CYP17A1, CYP11A1) are not expressed in normal prostate tissue. Alterations in their expression levels in steroidogenic tissues are closely associated with carcinogenesis, and, most importantly, with the development of aggressive forms of prostate cancer. Hence, it is important, to study how expression of CYPs in the prostate might be regulated, to understand the mechanisms of disease development and to improve the effectiveness of therapy. Several CYPs (CYP3A43, CYP2B6, CYP27A1, CYP24A1) can be considered as prognostic and diagnostic markers of prostate cancer. To propose personalized treatment, individual differences in CYP expression should be taken into account.

Liquid Chromatography-Tandem Mass Spectrometry (LC-MS/MS)-Based Proteomics of Drug-Metabolizing Enzymes and Transporters

Liquid chromatography-tandem mass spectrometry (LC-MS/MS)-based proteomics is a powerful tool for identifying and quantifying proteins in biological samples, outperforming conventional antibody-based methods in many aspects. LC-MS/MS-based proteomics studies have revealed the protein abundances of many drug-metabolizing enzymes and transporters (DMETs) in tissues relevant to drug metabolism and disposition. Previous studies have consistently demonstrated marked interindividual variability in DMET protein expression, suggesting that varied DMET function is an important contributing factor for interindividual variability in pharmacokinetics (PK) and pharmacodynamics (PD) of medications. Moreover, differential DMET expression profiles were observed across different species and in vitro models. Therefore, caution must be exercised when extrapolating animal and in vitro DMET proteomics findings to humans. In recent years, DMET proteomics has been increasingly utilized for the development of physiologically based pharmacokinetic models, and DMET proteins have also been proposed as biomarkers for prediction of the PK and PD of the corresponding substrate drugs. In sum, despite the existence of many challenges in the analytical technology and data analysis methods of LC-MS/MS-based proteomics, DMET proteomics holds great potential to advance our understanding of PK behavior at the individual level and to optimize treatment regimens via the DMET protein biomarker-guided precision pharmacotherapy.

The cytochrome P450 enzyme CYP24A1 increases proliferation of mutant KRAS-dependent lung adenocarcinoma independent of its catalytic activity

We previously reported that overexpression of cytochrome P450 family 24 subfamily A member 1 (CYP24A1) increases lung cancer cell proliferation by activating RAS signaling and that CYP24A1 knockdown inhibits tumor growth. However, the mechanism of CYP24A1-mediated cancer cell proliferation remains unclear. Here, we conducted cell synchronization and biochemical experiments in lung adenocarcinoma cells, revealing a link between CYP24A1 and anaphase-promoting complex (APC), a key cell cycle regulator. We demonstrate that CYP24A1 expression is cell cycle-dependent; it was higher in the G₂-M phase and diminished upon G₁ entry. CYP24A1 has a functional destruction box (D-box) motif that allows binding with two APC adaptors, CDC20-homologue 1 (CDH1) and cell division cycle 20 (CDC20). Unlike other APC substrates, however, CYP24A1 acted as a pseudo-substrate, inhibiting CDH1 activity and promoting mitotic progression. Conversely, overexpression of a CYP24A1 D-box mutant compromised CDH1 binding, allowing CDH1 hyperactivation, thereby hastening degradation of its substrates cyclin B1 and CDC20, and accumulation of the CDC20 substrate p21, prolonging mitotic exit. These activities also occurred with a CYP24A1 isoform 2 lacking the catalytic cysteine (Cys-462), suggesting that CYP24A1's oncogenic potential is independent of its catalytic activity. CYP24A1 degradation reduced clonogenic survival of mutant KRAS-driven lung cancer cells, and calcitriol treatment increased CYP24A1 levels and tumor burden in Lsl-KRAS^G12D mice. These results disclose a catalytic activity-independent growth-promoting role of CYP24A1 in mutant KRAS-driven lung cancer. This suggests that CYP24A1 could be therapeutically targeted in lung cancers in which its expression is high.

Molecular cloning, characterization and gene expression analysis of aminolevulinic acid synthase in Litopenaeus vannamei

5-Aminolevulinic acid synthase (ALAS) is the rate-limiting enzyme in the biosynthesis of heme, a prosthetic group that is found in hemoproteins, including those involved in molting. To better understand the roles of ALAS in L. vannamei (LvALAS), we analyzed its sequence and tissue distribution, the effects of age and bacterial infection on its gene expression, and the effects of LvALAS gene silencing. We also examined the expressions of three hemoproteins, the cytochrome oxidase subunit I (COX I) and subunit IV (COX IV) and catalase. Three LvALAS splicing variants were found in the hepatopancreas, with the main splicing variant having an open reading frame that encodes 532 aa. LvALAS transcripts were found in each of the eleven tissues tested in this study, with the highest gene expression in the intestine. The transcript abundances of LvALAS, COX I and COX IV in the hepatopancreas and stomach tended to decrease with age. LvALAS and catalase gene expressions significantly increased in the stomach after V. parahaemolyticus infection. LvALAS gene expression in the hepatopancreas, stomach and intestine (12- and 24-hours post-injection) was relatively lower in dsALAS-injected shrimp than in PBS-injected shrimp. All the PBS-injected shrimp molted after 8-10 days while no molting activity was observed in the dsALAS-injected shrimp group within the 14 days post-injection period. Our results provide evidence that (1) only the housekeeping form of ALAS exists in L. vannamei; LvALAS gene expression (2) decreases with age and (3) increases after bacterial infection; and (4) an ALAS-dependent pathway is necessary for proper molting in L. vannamei.

Subcellular dynamics of variants of SG2NA in NIH3T3 fibroblasts

SG2NA, a WD40 repeat protein of the Striatin subfamily, has four splicing and one messenger RNA edit variants. It is fast emerging as a scaffold for multimeric signaling complexes with roles in tissue development and disease. The green fluorescent protein (GFP)-tagged variants of SG2NA were ectopically expressed in NIH3T3 cells and their modulation by serum and GSK3β-ERK signaling were monitored. The 87, 78, and 35 kDa variants showed a biphasic modulation by serum till 24 h but the 52 kDa variant remained largely unresponsive. Inhibition of phosphatases by okadaic acid increased the levels of the endogenous 78 kDa and the ectopically expressed GFP-tagged 87 and 78 kDa SG2NAs. Contrastingly, okadaic acid treatment reduced the level of GFP-tagged 35 kDa SG2NA, suggesting differential modes of their stability through phosphorylation-dephosphorylation. The inhibition of GSK3β by LiCl showed a gradual decrease in the levels of 78 kDa. In the case of the other variants viz, GFP-tagged 35, 52, and 87 kDa, inhibition of GSK3β caused an initial increase followed by a decrease with a subtle difference in kinetics and intensities. Similar results were also seen upon inhibition of GSK3β by small interfering RNA. All the variants showed an increase followed by a decrease upon inhibition of extracellular-signal-regulated-kinase (ERK). These variants are localized in the plasma membrane, endoplasmic reticulum, mitochondria, and the nucleus with different propensities and no discernable subcellular distribution was seen upon stimulation by serum and the inhibition of phosphatases, GSK3β, and ERK. Taken together, the variants of SG2NA are modulated by the kinase-phosphatase network in a similar but characteristic manner.

The UGTome: The expanding diversity of UDP glycosyltransferases and its impact on small molecule metabolism

The UDP glycosyltransferase (UGT) superfamily of enzymes is responsible for the metabolism and clearance of thousands of lipophilic chemicals including drugs, toxins and endogenous signaling molecules. They provide a protective interface between the organism and its chemical-rich environment, as well as controlling critical signaling pathways to maintain healthy tissue function. UGTs are associated with drug responses and interactions, as well as a wide range of diseases including cancer. The human genome contains 22 UGT genes; however as befitting their exceptionally diverse substrate ranges and biological activities, the output of these UGT genes is functionally diversified by multiple processes including alternative splicing, post-translational modification, homo- and hetero-oligomerization, and interactions with other proteins. All UGT genes are subject to extensive alternative splicing generating variant/truncated UGT proteins with altered functions including the capacity to dominantly modulate/inhibit cognate full-length forms. Heterotypic oligomerization of different UGTs can alter kinetic properties relative to monotypic complexes, and potentially produce novel substrate specificities. Moreover, the recently profiled interactions of UGTs with non-UGT proteins may facilitate coordination between different metabolic processes, as well as providing opportunities for UGTs to engage in novel 'moonlighting' functions. Herein we provide a detailed and comprehensive review of all known modes of UGT functional diversification and propose a UGTome model to describe the resulting expansion of metabolic capacity and its potential to modulate drug/xenobiotic responses and cell behaviours in normal and disease contexts.

Cytochrome P450: genotype to phenotype

The cytochromes P450 comprise a family of enzymes that are responsible for around three-quarters of all drug metabolism reactions that occur in human populations. Many isoforms of cytochrome P450 exist but most reactions are undertaken by CYP2C9, CYP2C19, CYP2D6 and CYP3A4. This brief review focusses on the first three isozymes which exhibit polymorphism of phenotype.If there is a wide variation in drug metabolising capacity within the population, this may precipitate clinical consequences and influence the drug treatment of patients. Such problems range from a lack of efficacy to unanticipated toxicity. In order to minimise untoward events and "personalise" a patient's treatment, efforts have been made to discover an individual's drug metabolism status. This requires knowledge of the subject's phenotype at the time of clinical treatment. Since such testing is difficult, time-consuming and costly, the simpler approach of genotyping has been advocated.However, the correlation between genotype and phenotype is not good, with values of up to 50% misprediction being reported. Genotype-assisted forecasts cannot therefore be used with confidence to replace actual phenotype measurements. Obfuscating factors discussed include gene splicing, single nucleotide polymorphisms, epigenetics and microRNA, transcription regulation and multiple gene copies.

Alternative splicing of the vitamin D receptor modulates target gene expression and promotes ligand-independent functions

Alternative splicing modulates gene function by creating splice variants with alternate functions or non-coding RNA activity. Naturally occurring variants of nuclear receptor (NR) genes with dominant negative or gain-of-function phenotypes have been documented, but their cellular roles, regulation, and responsiveness to environmental stress or disease remain unevaluated. Informed by observations that class I androgen and estrogen receptor variants display ligand-independent signaling in human cancer tissues, we questioned whether the function of class II NRs, like the vitamin D receptor (VDR), would also respond to alternative splicing regulation. Artificial VDR constructs lacking exon 3 (Dex3-VDR), encoding part of the DNA binding domain (DBD), and exon 8 (Dex8-VDR), encoding part of the ligand binding domain (LBD), were transiently transfected into DU-145 cells and stably-integrated into Caco-2 cells to study their effect on gene expression and cell viability. Changes in VDR promoter signaling were monitored by the expression of target genes (e.g. CYP24A1, CYP3A4 and CYP3A5). Ligand-independent VDR signaling was observed in variants lacking exon 8, and a significant loss of gene suppressor function was documented for variants lacking exon 3. The gain-of-function behavior of the Dex8-VDR variant was recapitulated in vitro using antisense oligonucleotides (ASO) that induce the skipping of exon 8 in wild-type VDR. ASO targeting the splice acceptor site of exon 8 significantly stimulated ligand-independent VDR reporter activity and the induction of CYP24A1 above controls. These results demonstrate how alternative splicing can re-program NR gene function, highlighting novel mechanisms of toxicity and new opportunities for the use of splice-switching oligonucleotides (SSO) in precision medicine.

Effects of 5-Aminolevulinic Acid on Gene Expression, Immunity, and ATP Levels in Pacific White Shrimp, Litopenaeus vannamei

With the emergence of several infectious diseases in shrimp aquaculture, there is a growing interest in the use of feed additives to enhance shrimp immunity. Recently, the use of 5-aminolevulinic acid (5-ALA), a non-protein amino acid that plays a rate-limiting role in heme biosynthesis, has received attention for its positive effect on immunity in livestock animals. To evaluate the effect of 5-ALA in the Pacific white shrimp, Litopenaeus vannamei, we conducted microarray analysis, a Vibrio parahaemolyticus immersion challenge test, an ATP level assay, and gene expression analysis of some hemoproteins and genes associated with heme synthesis and degradation. Out of 15,745 L. vannamei putative genes on the microarray, 101 genes were differentially expressed by more than fourfold (p < 0.05) between 5-ALA-supplemented and control shrimp hepatopancreas. 5-ALA upregulated 99 of the 101 genes, 41 of which were immune- and defense-related genes based on sequence homology. Compared to the control, the 5-ALA-supplemented group had a higher survival rate in the challenge test, higher transcript levels of porphobilinogen synthase, ferrochelatase, catalase, nuclear receptor E75, and heme oxygenase-1 and higher levels of ATP. These findings suggest that dietary 5-ALA enhanced the immune response of L. vannamei to V. parahaemolyticus, upregulated immune- and defense-related genes, and enhanced aerobic energy metabolism, respectively. Further studies are needed to elucidate the extent of 5-ALA use in shrimp culture.

N-Acetyl-2-Aminofluorene (AAF) Processing in Adult Rat Hepatocytes in Primary Culture Occurs by High-Affinity Low-Velocity and Low-Affinity High-Velocity AAF Metabolite-Forming Systems

N-acetyl-2-aminofluorene (AAF) is a procarcinogen used widely in physiological investigations of chemical hepatocarcinogenesis. Its metabolic pathways have been described extensively, yet little is known about its biochemical processing, growth cycle expression, and pharmacological properties inside living hepatocytes-the principal cellular targets of this hepatocarcinogen. In this report, primary monolayer adult rat hepatocyte cultures and high specific-activity [ring G-3 H]-N-acetyl-2-aminofluorene were used to extend previous observations of metabolic activation of AAF by highly differentiated, proliferation-competent hepatocytes in long-term cultures. AAF metabolism proceeded by zero-order kinetics. Hepatocytes processed significant amounts of procarcinogen (≈12 μg AAF/106 cells/day). Five ring-hydroxylated and one deacetylated species of AAF were secreted into the culture media. Extracellular metabolite levels varied during the growth cycle (days 0-13), but their rank quantitative order was time invariant: 5-OH-AAF > 7-OH-AAF > 3-OH-AAF > N-OH-AAF > aminofluorene (AF) > 1-OH-AAF. Lineweaver-Burk analyses revealed two principal classes of metabolism: System I (high-affinity and low-velocity), Km[APPARENT] = 1.64 × 10-7  M and VMAX[APPARENT] = 0.1 nmol/106 cells/day and System II (low-affinity and high-velocity), Km[APPARENT] = 3.25 × 10-5  M and VMAX[APPARENT] = 1000 nmol/106 cells/day. A third system of metabolism of AAF to AF, with Km[APPARENT] and VMAX[APPARENT] constants of 9.6 × 10-5  M and 4.7 nmol/106 cells/day, was also observed. Evidence provided in this report and its companion paper suggests selective roles and intracellular locations for System I- and System II-mediated AAF metabolite formation during hepatocarcinogenesis, although some of the molecules and mechanisms responsible for multi-system processing remain to be fully defined.

Prologue

All living things are associated with a boundary defined ecological niche. Steady state conditions are rarely constant but evolutionary adaptation is too slow to adapt to daily threats so a surrogate variation mechanism is necessary. The genome defines the most basic instructions for life so that a molecular biology perspective provides the foundation for understanding resilience. Variations in the expression of RNA offers rapid variation and this book proposes this is the basis of resilience. This book attempts to illuminate mechanisms of resilience beginning with elaborating threats leading to disruption in steady state conditions. Recognition of threats and defense systems are described followed by adaptive changes in gene expression that refine responses. Finally, environmental conditions are discussed that serve to dampen the adaptive response oscillator to disruptive threats at the level of RNA expression.

This prologue is intended to acquaint the reader with my background and the genesis of optimism for an idea that the benefit of transcriptome plasticity is resilience. I grew up in several National Parks, remote regions of the United States that are set aside to preserve natural environments. I attended 12 schools by the time I graduated from high school, a fact that forced me to develop personal resilience. My career path as a scientist followed a path from ecologist to pharmacologist to molecular biologist. I was a professor that transitioned to biotechnology ensuring research subjects involving very diverse in subject matter so I appreciate the value of plasticity.