Toxicological and pharmacokinetic properties of sucralose-6-acetate and its parent sucralose: in vitro screening assays.

The purpose of this study was to determine the toxicological and pharmacokinetic properties of sucralose-6-acetate, a structural analog of the artificial sweetener sucralose. Sucralose-6-acetate is an intermediate and impurity in the manufacture of sucralose, and recent commercial sucralose samples were found to contain up to 0.67% sucralose-6-acetate. Studies in a rodent model found that sucralose-6-acetate is also present in fecal samples with levels up to 10% relative to sucralose which suggest that sucralose is also acetylated in the intestines. A MultiFlow® assay, a high-throughput genotoxicity screening tool, and a micronucleus (MN) test that detects cytogenetic damage both indicated that sucralose-6-acetate is genotoxic. The mechanism of action was classified as clastogenic (produces DNA strand breaks) using the MultiFlow® assay. The amount of sucralose-6-acetate in a single daily sucralose-sweetened drink might far exceed the threshold of toxicological concern for genotoxicity (TTCgenotox) of 0.15 µg/person/day. The RepliGut® System was employed to expose human intestinal epithelium to sucralose-6-acetate and sucralose, and an RNA-seq analysis was performed to determine gene expression induced by these exposures. Sucralose-6-acetate significantly increased the expression of genes associated with inflammation, oxidative stress, and cancer with greatest expression for the metallothionein 1 G gene (MT1G). Measurements of transepithelial electrical resistance (TEER) and permeability in human transverse colon epithelium indicated that sucralose-6-acetate and sucralose both impaired intestinal barrier integrity. Sucralose-6-acetate also inhibited two members of the cytochrome P450 family (CYP1A2 and CYP2C19). Overall, the toxicological and pharmacokinetic findings for sucralose-6-acetate raise significant health concerns regarding the safety and regulatory status of sucralose itself.

Genome and transcriptome sequencing of the green bottle fly, Lucilia sericata, reveals underlying factors of sheep flystrike and maggot debridement therapy.

The common green bottle blow fly Lucilia sericata (family, Calliphoridae) is widely used for maggot debridement therapy, which involves the application of sterile maggots to wounds. The larval excretions and secretions are important for consuming necrotic tissue and inhibiting bacterial growth in wounds of patients. Lucilia sericata is also of importance as a pest of sheep and in forensic studies to estimate a postmortem interval. Here we report the assembly of a 565.3 Mb genome from long read PacBio DNA sequencing of genomic DNA. The genome contains 14,704 predicted protein coding genes and 1709 non-coding genes. Targeted annotation and transcriptional analyses identified genes that are highly expressed in the larval salivary glands (secretions) and Malpighian tubules (excretions) under normal growth conditions and following heat stress. The genomic resources will underpin future genetic studies and in development of engineered strains for genetic control of L. sericata and for biotechnology-enhanced maggot therapy.

High-resolution characterization of the structural features and genetic variation of six feline leukocyte antigen class I loci via single molecule, real-time (SMRT) sequencing.

Of the 12 full-length feline leukocyte antigen class I (FLAI) loci, 3 are presumed to be classical: FLAI-E, FLAI-H, and FLAI-K. As diversity is a class Ia hallmark, multi-allelism is an important surrogate supporting a classical designation, in the absence of direct demonstration of T-cell restriction. Conversely, limited polymorphism at an expressed locus suggests regulation of immune effectors with invariant receptors, and non-classical status. FLAI-A, FLAI-J, FLAI-L, and FLAI-O are putative class Ib genes in cats. For both classes, identifying prevalent variants across outbred populations can illuminate specific genotypes to be prioritized for immune studies, as shared alleles direct shared responses. Since variation is concentrated in exons 2 and 3, which encode the antigen-binding domains, partial-length cloning/sequencing can be used for allele discovery, but is laborious and occasionally ambiguous. Here we develop a targeted approach to FLAI genotyping, using the single-molecule real-time (SMRT) platform, which allows full-length (3.4-kb) reads without assembly. Consensus sequences matched full-length Sanger references. Thirty-one new class Ia genes were found in 17 cats. Alleles segregated strongly by loci, and the origins of formerly difficult-to-assign sequences were resolved. Although not targeted, FLAI-L and FLAI-J, and the pseudogene FLAI-F, were also returned. Eighteen class Ib alleles were identified. Diversity was restricted and outside hypervariable regions. Both class Ib genes were transcriptionally active. Novel alternative splicing of FLAI-L was observed. SMRT sequencing of FLAI amplicons is useful for full-length genotyping at feline class Ia loci. High-throughput sequencing could allow highly accurate allele surveys in large cat cohorts.

Population genetic diversity in zebrafish lines.

Toxicological and pharmacological researchers have seized upon the many benefits of zebrafish, including the short generation time, well-characterized development, and early maturation as clear embryos. A major difference from many model organisms is that standard husbandry practices in zebrafish are designed to maintain population diversity. While this diversity is attractive for translational applications in human and ecological health, it raises critical questions on how interindividual genetic variation might contribute to chemical exposure or disease susceptibility differences. Findings from pooled samples of zebrafish support this supposition of diversity yet cannot directly measure allele frequencies for reference versus alternate alleles. Using the Tanguay lab Tropical 5D zebrafish line (T5D), we performed whole genome sequencing on a large group (n = 276) of individual zebrafish embryos. Paired-end reads were collected on an Illumina 3000HT, then aligned to the most recent zebrafish reference genome (GRCz10). These data were used to compare observed population genetic variation across species (humans, mice, zebrafish), then across lines within zebrafish. We found more single nucleotide polymorphisms (SNPs) in T5D than have been reported in SNP databases for any of the WIK, TU, TL, or AB lines. We theorize that some subset of the novel SNPs may be shared with other zebrafish lines but have not been identified in other studies due to the limitations of capturing population diversity in pooled sequencing strategies. We establish T5D as a model that is representative of diversity levels within laboratory zebrafish lines and demonstrate that experimental design and analysis can exert major effects when characterizing genetic diversity in heterogeneous populations.

A Cytoplasmic RNA Virus Alters the Function of the Cell Splicing Protein SRSF2.

To replicate efficiently, viruses must create favorable cell conditions and overcome cell antiviral responses. We previously reported that the reovirus protein µ2 from strain T1L, but not strain T3D, represses one antiviral response: alpha/beta interferon signaling. We report here that T1L, but not T3D, µ2 localizes to nuclear speckles, where it forms a complex with the mRNA splicing factor SRSF2 and alters its subnuclear localization. Reovirus replicates in cytoplasmic viral factories, and there is no evidence that reovirus genomic or messenger RNAs are spliced, suggesting that T1L µ2 might target splicing of cell RNAs. Indeed, RNA sequencing revealed that reovirus T1L, but not T3D, infection alters the splicing of transcripts for host genes involved in mRNA posttranscriptional modifications. Moreover, depletion of SRSF2 enhanced reovirus replication and cytopathic effect, suggesting that T1L µ2 modulation of splicing benefits the virus. This provides the first report of viral antagonism of the splicing factor SRSF2 and identifies the viral protein that determines strain-specific differences in cell RNA splicing.IMPORTANCE Efficient viral replication requires that the virus create favorable cell conditions. Many viruses accomplish this by repressing specific antiviral responses. We demonstrate here that some mammalian reoviruses, RNA viruses that replicate strictly in the cytoplasm, express a protein variant that localizes to nuclear speckles, where it targets a cell mRNA splicing factor. Infection with a reovirus strain that targets this splicing factor alters splicing of cell mRNAs involved in the maturation of many other cell mRNAs. Depletion of this cell splicing factor enhances reovirus replication and cytopathic effect. Our results provide the first evidence of viral antagonism of this splicing factor and suggest that downstream consequences to the cell are global and benefit the virus.

Identification of Cytauxzoon felis antigens via protein microarray and assessment of expression library immunization against cytauxzoonosis.

BACKGROUND: Cytauxzoonosis is a disease of felids in North America caused by the tick-transmitted apicomplexan parasite Cytauxzoon felis. Cytauxzoonosis is particularly virulent for domestic cats, but no vaccine currently exists. The parasite cannot be cultivated in vitro, presenting a significant limitation for vaccine development. METHODS: Recent sequencing of the C. felis genome has identified over 4300 putative protein-encoding genes. From this pool we constructed a protein microarray containing 673 putative C. felis proteins. This microarray was probed with sera from C. felis-infected and naïve cats to identify differentially reactive antigens which were incorporated into two expression library vaccines, one polyvalent and one monovalent. We assessed the efficacy of these vaccines to prevent of infection and/or disease in a tick-challenge model. RESULTS: Probing of the protein microarray resulted in identification of 30 differentially reactive C. felis antigens that were incorporated into the two expression library vaccines. However, expression library immunization failed to prevent infection or disease in cats challenged with C. felis. CONCLUSIONS: Protein microarray facilitated high-throughput identification of novel antigens, substantially increasing the pool of characterized C. felis antigens. These antigens should be considered for development of C. felis vaccines, diagnostics, and therapeutics.

Relative Quantification and Higher-Order Modeling of the Plasma Glycan Cancer Burden Ratio in Ovarian Cancer Case-Control Samples.

An early-stage, population-wide biomarker for ovarian cancer (OVC) is essential to reverse its high mortality rate. Aberrant glycosylation by OVC has been reported, but studies have yet to identify an N-glycan with sufficiently high specificity. We curated a human biorepository of 82 case-control plasma samples, with 27%, 12%, 46%, and 15% falling across stages I-IV, respectively. For relative quantitation, glycans were analyzed by the individuality normalization when labeling with glycan hydrazide tags (INLIGHT) strategy for enhanced electrospray ionization, MS/MS analysis. Sixty-three glycan cancer burden ratios (GBRs), defined as the log10 ratio of the case-control extracted ion chromatogram abundances, were calculated above the limit of detection. The final GBR models, built using stepwise forward regression, included three significant terms: OVC stage, normalized mean GBR, and tag chemical purity; glycan class, fucosylation, or sialylation were not significant variables. After Bonferroni correction, seven N-glycans were identified as significant (p < 0.05), and after false discovery rate correction, an additional four glycans were determined to be significant (p < 0.05), with one borderline (p = 0.05). For all N-glycans, the vectors of the effects from stages II-IV were sequentially reversed, suggesting potential biological changes in OVC morphology or in host response.

Questionnaire-based exposome-wide association studies for common diseases in the Personalized Environment and Genes Study.

The exposome collectively refers to all exposures, beginning in utero and continuing throughout life, and comprises not only standard environmental exposures such as point source pollution and ozone levels but also exposures from diet, medication, lifestyle factors, stress, and occupation. The exposome interacts with individual genetic and epigenetic characteristics to affect human health and disease, but large-scale studies that characterize the exposome and its relationships with human disease are limited. To address this gap, we used extensive questionnaire data from the diverse North Carolina-based Personalized Environment and Genes Study (PEGS, n = 9, 429) to evaluate exposure associations in relation to common diseases. We performed an exposome-wide association study (ExWAS) to examine single exposure models and their associations with 11 common complex diseases, namely allergic rhinitis, asthma, bone loss, fibroids, high cholesterol, hypertension, iron-deficient anemia, ovarian cysts, lower GI polyps, migraines, and type 2 diabetes. Across diseases, we found associations with lifestyle factors and socioeconomic status as well as asbestos, various dust types, biohazardous material, and textile-related exposures. We also found disease-specific associations such as fishing with lead weights and migraines. To differentiate between a replicated result and a novel finding, we used an AI-based literature search and database tool that allowed us to examine the current literature. We found both replicated findings, especially for lifestyle factors such as sleep and smoking across diseases, and novel findings, especially for occupational exposures and multiple diseases.

Exposure to PFAS chemicals induces sex-dependent alterations in key rate-limiting steps of lipid metabolism in liver steatosis.

Toxicants with the potential to bioaccumulate in humans and animals have long been a cause for concern, particularly due to their association with multiple diseases and organ injuries. Per- and polyfluoro alkyl substances (PFAS) and polycyclic aromatic hydrocarbons (PAH) are two such classes of chemicals that bioaccumulate and have been associated with steatosis in the liver. Although PFAS and PAH are classified as chemicals of concern, their molecular mechanisms of toxicity remain to be explored in detail. In this study, we aimed to identify potential mechanisms by which an acute exposure to PFAS and PAH chemicals can induce lipid accumulation and whether the responses depend on chemical class, dose, and sex. To this end, we analyzed mechanisms beginning with the binding of the chemical to a molecular initiating event (MIE) and the consequent transcriptomic alterations. We collated potential MIEs using predictions from our previously developed ToxProfiler tool and from published steatosis adverse outcome pathways. Most of the MIEs are transcription factors, and we collected their target genes by mining the TRRUST database. To analyze the effects of PFAS and PAH on the steatosis mechanisms, we performed a computational MIE-target gene analysis on high-throughput transcriptomic measurements of liver tissue from male and female rats exposed to either a PFAS or PAH. The results showed peroxisome proliferator-activated receptor (PPAR)-α targets to be the most dysregulated, with most of the genes being upregulated. Furthermore, PFAS exposure disrupted several lipid metabolism genes, including upregulation of fatty acid oxidation genes (Acadm, Acox1, Cpt2, Cyp4a1-3) and downregulation of lipid transport genes (Apoa1, Apoa5, Pltp). We also identified multiple genes with sex-specific behavior. Notably, the rate-limiting genes of gluconeogenesis (Pck1) and bile acid synthesis (Cyp7a1) were specifically downregulated in male rats compared to female rats, while the rate-limiting gene of lipid synthesis (Scd) showed a PFAS-specific upregulation. The results suggest that the PPAR signaling pathway plays a major role in PFAS-induced lipid accumulation in rats. Together, these results show that PFAS exposure induces a sex-specific multi-factorial mechanism involving rate-limiting genes of gluconeogenesis and bile acid synthesis that could lead to activation of an adverse outcome pathway for steatosis.

Interactive data sharing for multiple questionnaire-based exposome-wide association studies and exposome correlations in the Personalized Environment and Genes Study.

The correlations among individual exposures in the exposome, which refers to all exposures an individual encounters throughout life, are important for understanding the landscape of how exposures co-occur, and how this impacts health and disease. Exposome-wide association studies (ExWAS), which are analogous to genome-wide association studies (GWAS), are increasingly being used to elucidate links between the exposome and disease. Despite increased interest in the exposome, tools and publications that characterize exposure correlations and their relationships with human disease are limited, and there is a lack of data and results sharing in resources like the GWAS catalog. To address these gaps, we developed the PEGS Explorer web application to explore exposure correlations in data from the diverse North Carolina-based Personalized Environment and Genes Study (PEGS) that were rigorously calculated to account for differing data types and previously published results from ExWAS. Through globe visualizations, PEGS Explorer allows users to explore correlations between exposures found to be associated with complex diseases. The exposome data used for analysis includes not only standard environmental exposures such as point source pollution and ozone levels but also exposures from diet, medication, lifestyle factors, stress, and occupation. The web application addresses the lack of accessible data and results sharing, a major challenge in the field, and enables users to put results in context, generate hypotheses, and, importantly, replicate findings in other cohorts. PEGS Explorer will be updated with additional results as they become available, ensuring it is an up-to-date resource in exposome science.

Computational and phylogenetic validation of nematode horizontal gene transfer.

Sequencing of expressed genes has shown that nematodes, particularly the plant-parasitic nematodes, have genes purportedly acquired from other kingdoms by horizontal gene transfer. The prevailing orthodoxy is that such transfer has been a driving force in the evolution of niche specificity, and a recent paper in BMC Evolutionary Biology that presents a detailed phylogenetic analysis of cellulase genes in the free-living nematode Pristionchus pacificus at the species, genus and family levels substantiates this hypothesis.

The Genome of the Poecilogonous Annelid Streblospio benedicti.

Streblospio benedicti is a common marine annelid that has become an important model for developmental evolution. It is the only known example of poecilogony (where two distinct developmental modes occur within a single species) that is due to a heritable difference in egg size. The dimorphic developmental programs and life-histories exhibited in this species depend on differences within the genome, making it an optimal model for understanding the genomic basis of developmental divergence. Studies using S. benedicti have begun to uncover the genetic and genomic principles that underlie developmental uncoupling, but until now they have been limited by the lack of availability of genomic tools. Here, we present an annotated chromosomal-level genome assembly of S. benedicti generated from a combination of Illumina reads, Nanopore long reads, Chicago and Hi-C chromatin interaction sequencing, and a genetic map from experimental crosses. At 701.4 Mb, the S. benedicti genome is the largest annelid genome to date that has been assembled to chromosomal scaffolds. The complete genome of S. benedicti is valuable for functional genomic analyses of development and evolution, as well as phylogenetic comparison within the annelida and the Lophotrochozoa. Despite having two developmental modes, there is no evidence of genome duplication or substantial gene number expansions. Instead, lineage-specific repeats account for much of the expansion of this genome compared with other annelids.

Proteomic and bioinformatic analysis of the root-knot nematode Meloidogyne hapla: the basis for plant parasitism.

On the basis of the complete genome sequence of the root-knot nematode Melodogyne hapla, we have deduced and annotated the entire proteome of this plant-parasite to create a database of 14,420 proteins. We have made this database, termed HapPep3, available from the Superfamily repository of model organism proteomes (http://supfam.mrc-lmb.cam.ac.uk/SUPERFAMILY). To experimentally confirm the HapPep3 assignments using proteomics, we applied a data-independent LC/MS(E) analysis to M. hapla protein extracts fractionated by SDS-PAGE. A total of 516 nonredundant proteins were identified with an average of 9 unique peptides detected per protein. Some proteins, including examples with complex gene organization, were defined by more than 20 unique peptide matches, thus, providing experimental confirmation of computational predictions of intron/exon structures. On the basis of comparisons of the broad physicochemical properties of the experimental and computational proteomes, we conclude that the identified proteins reflect a true and unbiased sampling of HapPep3. Conversely, HapPep3 appears to broadly cover the protein space able to be experimentally sampled. To estimate the false discovery rate, we queried human, plant, and bacterial databases for matches to the LC/MS(E)-derived peptides, revealing fewer than 1% of matches, most of which were to highly conserved proteins. To provide a functional comparison of the acquired and deduced proteomes, each was subjected to higher order annotation, including comparisons of Gene Ontology, protein domains, signaling, and localization predictions, further indicating concordance, although those proteins that did deviate seem to be highly significant. Approximately 20% of the experimentally sampled proteome was predicted to be secreted, and thus potentially play a role at the host-parasite interface. We examined reference pathways to determine the extent of proteome similarity of M. hapla to that of the free-living nematode, Caenorhabditis elegans, revealing significant similarities and differences. Collectively, the analyzed protein set provides an initial foundation to experimentally dissect the basis of plant parasitism by M. hapla.

Complete Genome Sequences of Six Lactobacilli Isolated from American Quarter Horses.

We report the complete circular genome sequences of six Lactobacillus strains and their plasmids, if any, from the fecal material of quarter horses at different ages.

Complete Genome Sequences of Lactobacillus Strains C25 and P38, Isolated from Chicken Cecum.

We report the complete circular genome sequences of Lactobacillus crispatus strain C25, its plasmid, and Lactobacillus animalis strain P38; both strains were isolated from the cecum of 4-week-old chickens. These isolates represent potential probiotic strains for poultry.

Genomic analyses of a livestock pest, the New World screwworm, find potential targets for genetic control programs.

The New World Screwworm fly, Cochliomyia hominivorax, is a major pest of livestock in South America and Caribbean. However, few genomic resources have been available for this species. A genome of 534 Mb was assembled from long read PacBio DNA sequencing of DNA from a highly inbred strain. Analysis of molecular evolution identified 40 genes that are likely under positive selection. Developmental RNA-seq analysis identified specific genes associated with each stage. We identify and analyze the expression of genes that are likely important for host-seeking behavior (chemosensory), development of larvae in open wounds in warm-blooded animals (heat shock protein, immune response) and for building transgenic strains for genetic control programs including gene drive (sex determination, germline). This study will underpin future experiments aimed at understanding the parasitic lifestyle of the screwworm fly and greatly facilitate future development of strains for efficient systems for genetic control of screwworm.

Elucidating Gene-by-Environment Interactions Associated with Differential Susceptibility to Chemical Exposure.

BACKGROUND: Modern societies are exposed to vast numbers of potentially hazardous chemicals. Despite demonstrated linkages between chemical exposure and severe health effects, there are limited, often conflicting, data on how adverse health effects of exposure differ across individuals. OBJECTIVES: We tested the hypothesis that population variability in response to certain chemicals could elucidate a role for gene-environment interactions (GxE) in differential susceptibility. METHODS: High-throughput screening (HTS) data on thousands of chemicals in genetically heterogeneous zebrafish were leveraged to identify a candidate chemical (Abamectin) with response patterns indicative of population susceptibility differences. We tested the prediction by generating genome-wide sequence data for 276 individual zebrafish displaying susceptible (Affected) vs. resistant (Unaffected) phenotypes following identical chemical exposure. RESULTS: We found GxE associated with differential susceptibility in the sox7 promoter region and then confirmed gene expression differences between phenotypic response classes. CONCLUSIONS: The results for Abamectin in zebrafish demonstrate that GxE associated with naturally occurring, population genetic variation play a significant role in mediating individual response to chemical exposure. https://doi.org/10.1289/EHP2662.

Polyphasic characterization of four soil-derived phenanthrene-degrading Acidovorax strains and proposal of Acidovorax carolinensis sp. nov.

Four bacterial strains identified as members of the Acidovorax genus were isolated from two geographically distinct but similarly contaminated soils in North Carolina, USA, characterized, and their genomes sequenced. Their 16S rRNA genes were highly similar to those previously recovered during stable-isotope probing (SIP) of one of the soils with the polycyclic aromatic hydrocarbon (PAH) phenanthrene. Heterotrophic growth of all strains occurred with a number of organic acids, as well as phenanthrene, but no other tested PAHs. Optimal growth occurred aerobically under mesophilic temperature, neutral pH, and low salinity conditions. Predominant fatty acids were C16:1ω7c/C16:1ω6c, C16:0, and C18:1ω7c, and were consistent with the genus. Genomic G+C contents ranged from 63.6 to 64.2%. A combination of whole genome comparisons and physiological analyses indicated that these four strains likely represent a single species within the Acidovorax genus. Chromosomal genes for phenanthrene degradation to phthalate were nearly identical to highly conserved regions in phenanthrene-degrading Delftia, Burkholderia, Alcaligenes, and Massilia species in regions flanked by transposable or extrachromosomal elements. The lower degradation pathway for phenanthrene metabolism was inferred by comparisons to described genes and proteins. The novel species Acidovorax carolinensis sp. nov. is proposed, comprising the four strains described in this study with strain NA3T as the type strain (=LMG 30136, =DSM 105008).

no blokes Is Essential for Male Viability and X Chromosome Gene Expression in the Australian Sheep Blowfly.

It has been hypothesized that the Drosophila 4th chromosome is derived from an ancient X chromosome [1]. In the Australian sheep blowfly, Lucilia cuprina, the heterochromatic X chromosome contains few active genes and orthologs of Drosophila X-linked genes are autosomal. Of 8 X-linked genes identified previously in L. cuprina, 6 were orthologs of Drosophila 4th-chromosome genes [2]. The X-linked genes were expressed equally in males and females. Here we identify an additional 51 X-linked genes and show that most are dosage compensated. Orthologs of 49 of the 59 X-linked genes are on the 4th chromosome in D. melanogaster. Because painting of fourth (Pof) is important for expression of Drosophila 4th-chromosome genes [3], we used Cas9 to make a loss-of-function knockin mutation in an L. cuprina Pof ortholog we call no blokes (nbl). Homozygous nbl males derived from homozygous nbl mothers die at the late pupal stage. Homozygous nbl females are viable, fertile, and live longer than heterozygous nbl females. RNA expression of most X-linked genes was reduced in homozygous nbl male pupae and to a lesser extent in nbl females compared to heterozygous siblings. The results suggest that NBL could be important for X chromosome dosage compensation in L. cuprina. NBL may also facilitate gene expression in the heterochromatic environment of the X chromosome in both sexes. This study supports the hypothesis on the origin of the Drosophila 4th chromosome and that a POF-like protein was required for normal gene expression on the ancient X chromosome.

Pharmacometabolomics Informs About Pharmacokinetic Profile of Methylphenidate.

Carboxylesterase 1 (CES1) metabolizes methylphenidate and other drugs. CES1 gene variation only partially explains pharmacokinetic (PK) variability. Biomarkers predicting the PKs of drugs metabolized by CES1 are needed. We identified lipids in plasma from 44 healthy subjects that correlated with CES1 activity as determined by PK parameters of methylphenidate including a ceramide (q value = 0.001) and a phosphatidylcholine (q value = 0.005). Carriers of the CES1 143E allele had decreased methylphenidate metabolism and altered concentration of this phosphatidylcholine (q value = 0.040) and several high polyunsaturated fatty acid lipids (PUFAs). The half-maximal inhibitory concentration (IC50 ) values of chenodeoxycholate and taurocholate were 13.55 and 19.51 µM, respectively, consistent with a physiological significance. In silico analysis suggested that bile acid inhibition of CES1 involved both binding to the active and superficial sites of the enzyme. We initiated identification of metabolites predicting PKs of drugs metabolized by CES1 and suggest lipids to regulate or be regulated by this enzyme.