Genetic polymorphism of CYP2U1, a cytochrome P450 involved in fatty acids hydroxylation

Transcriptomic profiling of hepatic tissues for drug metabolism genes in nonalcoholic fatty liver disease: A study of human and animals

Background

Drug metabolism genes are involved in the in vivo metabolic processing of drugs. In previous research, we found that a high-fat diet affected the transcript levels of mouse hepatic genes responsible for drug metabolism.

Aims

Our research intends to discover the drug metabolism genes that are dysregulated at the transcriptome level in nonalcoholic fatty liver disease (NAFLD).

Methods

We analyzed the transcriptome for drug metabolism genes of 35 human liver tissues obtained during laparoscopic cholecystectomy. Additionally, we imported transcriptome data from mice fed a high-fat diet in previous research and two open-access Gene Expression Omnibus (GEO) datasets (GSE63067 and GSE89632). Then, using quantitative real-time polymerase chain reaction (qRT-PCR), we cross-linked the differentially expressed genes (DEGs) in clinical and animal samples and validated the common genes.

Results

In this study, we identified 35 DEGs, of which 33 were up-regulated and two were down-regulated. Moreover, we found 71 DEGs (39 up- and 32 down-regulated), 276 DEGs (157 up- and 119 down-regulated), and 158 DEGs (117 up- and 41 down-regulated) in the GSE63067, GSE89632, and high-fat diet mice, respectively. Of the 35 DEGs, nine co-regulated DEGs were found in the Venn diagram (CYP20A1, CYP2U1, SLC9A6, SLC26A6, SLC31A1, SLC46A1, SLC46A3, SULT1B1, and UGT2A3).

Conclusion

Nine significant drug metabolism genes were identified in NAFLD. Future research should investigate the impacts of these genes on drug dose adjustment in patients with NAFLD.

Clinical Trial Registration

http://www.chictr.org.cn, identifier ChiCTR2100041714.

Detoxification Cytochrome P450s (CYPs) in Families 1-3 Produce Functional Oxylipins from Polyunsaturated Fatty Acids

This manuscript reviews the CYP-mediated production of oxylipins and the current known function of these diverse set of oxylipins with emphasis on the detoxification CYPs in families 1-3. Our knowledge of oxylipin function has greatly increased over the past 3-7 years with new theories on stability and function. This includes a significant amount of new information on oxylipins produced from linoleic acid (LA) and the omega-3 PUFA-derived oxylipins such as α-linolenic acid (ALA), docosahexaenoic acid (DHA), and eicosapentaenoic acid (EPA). However, there is still a lack of knowledge regarding the primary CYP responsible for producing specific oxylipins, and a lack of mechanistic insight for some clinical associations between outcomes and oxylipin levels. In addition, the role of CYPs in the production of oxylipins as signaling molecules for obesity, energy utilization, and development have increased greatly with potential interactions between diet, endocrinology, and pharmacology/toxicology due to nuclear receptor mediated CYP induction, CYP inhibition, and receptor interactions/crosstalk. The potential for diet-diet and diet-drug/chemical interactions is high given that these promiscuous CYPs metabolize a plethora of different endogenous and exogenous chemicals.

Human Orphan Cytochrome P450 2U1 Catalyzes the ω-Hydroxylation of Leukotriene B4

Cytochrome P450 2U1 (CYP2U1) identified from the human genome remains poorly known since few data are presently available on its physiological function(s) and substrate(s) specificity. CYP2U1 mutations are associated with complicated forms of hereditary spastic paraplegia, alterations of mitochondrial architecture and bioenergetics. In order to better know the biological roles of CYP2U1, we used a bioinformatics approach. The analysis of the data invited us to focus on leukotriene B₄ (LTB₄), an important inflammatory mediator. Here, we show that CYP2U1 efficiently catalyzes the hydroxylation of LTB₄ predominantly on its ω-position. We also report docking experiments of LTB₄ in a 3D model of truncated CYP2U1 that are in agreement with this hydroxylation regioselectivity. The involvement of CYP2U1 in the metabolism of LTB₄ could have strong physiological consequences in cerebral pathologies including ischemic stroke because CYP2U1 is predominantly expressed in the brain.

Genetic diversity and population structure of muscovy duck (Cairina moschata) from Nigeria

The domestic Muscovy duck (Cairina moschata) provide unique genetic resources patterned by both tropical environmental conditions and human activities, the evaluation of their genetic diversity and population structure will shade light on the mechanism of their remarkable adaptive capacities. We therefore analyzed the variation in mtDNA cytochrome b and nuclear DNA CYP2U1 sequences of 378 Nigerian Muscovy ducks (comprising of 287 de novo and 91 downloaded) plus 80 published sequences of Muscovy ducks from India. The results showed high haplotype diversity (0.800 ± 0.023) among Nigerian Muscovy duck populations with 91 distinct haplotypes for the nuclear DNA CYP2U1 gene but low (0.266 ± 0.033) for cytochrome b with 31 haplotypes. The median-joining networks of both markers grouped Nigerian Muscovy ducks into two; the first group consisting of only Nigerian Muscovy duck populations, and the second group Nigerian with Indian populations. Neutrality test results indicated that Nigerian populations experienced recent population expansion and/or genetic hitchhiking. A geographic signal was absent in line with previously studied poultry species in Nigeria. The most prominent haplotype dominated across all regions in Nigeria, which may be due to extensive genetic intermixing except for the Indian population (F _ST = 0.02550, P = 0.01075). This indicated low genetic differentiation between and within Nigerian Muscovy duck as revealed by the suitability of the nuclear DNA CYP2U1 gene.

Chromosome-Scale Genome of Masked Palm Civet (Paguma larvata) Shows Genomic Signatures of Its Biological Characteristics and Evolution

The masked palm civet (Paguma larvata) is a small carnivore with distinct biological characteristics, that likes an omnivorous diet and also serves as a vector of pathogens. Although this species is not an endangered animal, its population is reportedly declining. Since the severe acute respiratory syndrome (SARS) epidemic in 2003, the public has been particularly concerned about this species. Here, we present the first genome of the P. larvata, comprising 22 chromosomes assembled using single-tube long fragment read (stLFR) and Hi-C technologies. The genome length is 2.41 Gb with a scaffold N50 of 105.6 Mb. We identified the 107.13 Mb X chromosome and one 1.34 Mb Y-linked scaffold and validated them by resequencing 45 P. larvata individuals. We predicted 18,340 protein-coding genes, among which 18,333 genes were functionally annotated. Interestingly, several biological pathways related to immune defenses were found to be significantly expanded. Also, more than 40% of the enriched pathways on the positively selected genes (PSGs) were identified to be closely related to immunity and survival. These enriched gene families were inferred to be essential for the P. larvata for defense against the pathogens. However, we did not find a direct genomic basis for its adaptation to omnivorous diet despite multiple attempts of comparative genomic analysis. In addition, we evaluated the susceptibility of the P. larvata to the SARS-CoV-2 by screening the RNA expression of the ACE2 and TMPRSS2/TMPRSS4 genes in 16 organs. Finally, we explored the genome-wide heterozygosity and compared it with other animals to evaluate the population status of this species. Taken together, this chromosome-scale genome of the P. larvata provides a necessary resource and insights for understanding the genetic basis of its biological characteristics, evolution, and disease transmission control.

The Functions of Cytochrome P450 ω-hydroxylases and the Associated Eicosanoids in Inflammation-Related Diseases

The cytochrome P450 (CYP) ω-hydroxylases are a subfamily of CYP enzymes. While CYPs are the main metabolic enzymes that mediate the oxidation reactions of many endogenous and exogenous compounds in the human body, CYP ω-hydroxylases mediate the metabolism of multiple fatty acids and their metabolites via the addition of a hydroxyl group to the ω- or (ω-1)-C atom of the substrates. The substrates of CYP ω-hydroxylases include but not limited to arachidonic acid, docosahexaenoic acid, eicosapentaenoic acid, epoxyeicosatrienoic acids, leukotrienes, and prostaglandins. The CYP ω-hydroxylases-mediated metabolites, such as 20-hyroxyleicosatrienoic acid (20-HETE), 19-HETE, 20-hydroxyl leukotriene B4 (20-OH-LTB₄), and many ω-hydroxylated prostaglandins, have pleiotropic effects in inflammation and many inflammation-associated diseases. Here we reviewed the classification, tissue distribution of CYP ω-hydroxylases and the role of their hydroxylated metabolites in inflammation-associated diseases. We described up-regulation of CYP ω-hydroxylases may be a pathogenic mechanism of many inflammation-associated diseases and thus CYP ω-hydroxylases may be a therapeutic target for these diseases. CYP ω-hydroxylases-mediated eicosanods play important roles in inflammation as pro-inflammatory or anti-inflammatory mediators, participating in the process stimulated by cytokines and/or the process stimulating the production of multiple cytokines. However, most previous studies focused on 20-HETE,and further studies are needed for the function and mechanisms of other CYP ω-hydroxylases-mediated eicosanoids. We believe that our studies of CYP ω-hydroxylases and their associated eicosanoids will advance the translational and clinal use of CYP ω-hydroxylases inhibitors and activators in many diseases.

Alterations of Cytochrome P450s and UDP-Glucuronosyltransferases in Brain Under Diseases and Their Clinical Significances

Cytochrome P450s (CYPs) and UDP-glucuronosyltransferases (UGTs) are both greatly important metabolic enzymes in various tissues, including brain. Although expressions of brain CYPs and UGTs and their contributions to drug disposition are much less than liver, both CYPs and UGTs also mediate metabolism of endogenous substances including dopamine and serotonin as well as some drugs such as morphine in brain, demonstrating their important roles in maintenance of brain homeostasis or pharmacological activity of drugs. Some diseases such as epilepsy, Parkinson's disease and Alzheimer's disease are often associated with the alterations of CYPs and UGTs in brain, which may be involved in processes of these diseases via disturbing metabolism of endogenous substances or resisting drugs. This article reviewed the alterations of CYPs and UGTs in brain, the effects on endogenous substances and drugs and their clinical significances. Understanding the roles of CYPs and UGTs in brain provides some new strategies for the treatment of central nervous system diseases.

Differential Expression of P450 Genes and nAChR Subunits Associated With Imidacloprid Resistance in Laodelphax striatellus (Hemiptera: Delphacidae)

Imidacloprid is a key insecticide used for controlling sucking insect pests, including the small brown planthopper (Laodelphax striatellus, Fallén) (Hemiptera: Delphacidae), an important agricultural pest of rice. A strain of L. striatellus (YN-ILR) developed 21-fold resistance when selected with imidacloprid on a susceptible YN strain. An in vitro study on piperonyl butoxide synergism indicated that enhanced detoxification mediated by cytochrome P450s contributed to imidacloprid resistance to some extent, and multiple P450 genes showed altered expression in the imidacloprid-resistant YN-ILR strain compared with the susceptible YN strain (CYP425B1-CYP6BD10 had 1.51- to 11.45-fold higher expression, CYP4CE2-CYP4DD1V2 had 0.12- to 0.57-fold lower expression). While there were no mutations in target nicotinic acetylcholine receptor (nAChR) genes, subunits of Lsα1, Lsβ1, and Lsβ3 in the YN-ILR strain showed 3.86-, 4.39-, and 2.59-fold higher expression and Lsa8 displayed 0.38-fold lower expression than the YN strain. Moreover, 21-fold moderate imidacloprid resistance in individuals of L. striatellus did not produce a fitness cost. The findings suggest that L. striatellus has the capacity to develop resistance to imidacloprid through P450 detoxification and potential target nAChR expression changes, and moderate imidacloprid resistance was not associated with a fitness cost.

20-HETE promotes glucose-stimulated insulin secretion in an autocrine manner through FFAR1

The long-chain fatty acid receptor FFAR1 is highly expressed in pancreatic β-cells. Synthetic FFAR1 agonists can be used as antidiabetic drugs to promote glucose-stimulated insulin secretion (GSIS). However, the physiological role of FFAR1 in β-cells remains poorly understood. Here we show that 20-HETE activates FFAR1 and promotes GSIS via FFAR1 with higher potency and efficacy than dietary fatty acids such as palmitic, linoleic, and α-linolenic acid. Murine and human β-cells produce 20-HETE, and the ω-hydroxylase-mediated formation and release of 20-HETE is strongly stimulated by glucose. Pharmacological inhibition of 20-HETE formation and blockade of FFAR1 in islets inhibits GSIS. In islets from type-2 diabetic humans and mice, glucose-stimulated 20-HETE formation and 20-HETE-dependent stimulation of GSIS are strongly reduced. We show that 20-HETE is an FFAR1 agonist, which functions as an autocrine positive feed-forward regulator of GSIS, and that a reduced glucose-induced 20-HETE formation contributes to inefficient GSIS in type-2 diabetes.

FFAR1 receptor is highly expressed in beta cells and its activation has been suggested as therapy against type-2 diabetes. Here, Tunaru et al. show that 20-hydroxyeicosatetraenoic acid, produced within the islets upon glucose stimulation, acts in an autocrine manner to stimulate insulin secretion via FFAR1 activation.

De novo transcriptome sequencing and comparative analysis to discover genes involved in ovarian maturity in Strongylocentrotus nudus

Strongylocentrotus nudus is an edible sea urchin, mainly harvested in China. Correlation studies indicated that S. nudus with larger diameter have a prolonged marketing time and better palatability owing to their precocious gonads and extended maturation process. However, the molecular mechanism underlying this phenomenon is still unknown. Here, transcriptome sequencing was applied to study the ovaries of adult S. nudus with different shell diameters to explore the possible mechanism. In this study, four independent cDNA libraries were constructed, including two from the big size urchins and two from the small ones using a HiSeq™2500 platform. A total of 88,581 unigenes were acquired with a mean length of 1354bp, of which 66,331 (74.88%) unigenes could be annotated using six major publicly available databases. Comparative analysis revealed that 353 unigenes were differentially expressed (with log2(ratio)≥1, FDR≤0.001) between the two groups. Of these, 20 differentially expressed genes (DEGs) were selected to confirm the accuracy of RNA-seq data by quantitative real-time RT-PCR. Furthermore, gene ontology and KEGG pathway enrichment analyses were performed to find the putative genes and pathways related to ovarian maturity. Eight unigenes were identified as significant DEGs involved in reproduction related pathways; these included Mos, Cdc20, Rec8, YP30, cytochrome P450 2U1, ovoperoxidase, proteoliaisin, and rendezvin. Our research fills the gap in the studies on the S. nudus ovaries using transcriptome analysis.

Cytochrome P450 2U1, a very peculiar member of the human P450s family

Cytochrome P450 2U1 (CYP2U1) exhibits several distinctive characteristics among the 57 human CYPs, such as its presence in almost all living organisms with a highly conserved sequence, its particular gene organization with only five exons, its major location in thymus and brain, and its protein sequence involving an unusually long N-terminal region containing 8 proline residues and an insert of about 20 amino acids containing 5 arginine residues after the transmembrane helix. Few substrates, including fatty acids, N-arachidonoylserotonin (AS), and some drugs, have been reported so far. However, its biological roles remain largely unknown, even though CYP2U1 mutations have been involved in some pathological situations, such as complicated forms of hereditary spastic paraplegia. These data together with its ability to hydroxylate some fatty acids and AS suggest its possible role in lipid metabolism.

Genetic diversity and relationship of Indian Muscovy duck populations

Assessment of genetic diversity within and between populations is a prerequisite for sustainable utilization of domestic species. The domestic Muscovy duck (Cairina moschata) is an economically important species around the world for its unique meat taste and low-caloric content. It is one of the important domestic species in India as it ensures food security to the rural sectors. In this study, we have analyzed the genetic diversity and relationship of four Muscovy duck populations collected from different states (Assam, Mizoram, Odisha and Kerala) of India using mtDNA cytochrome b and nuclear DNA CYP2U1 genes. The results showed low genetic diversity among populations for both the genes. Kerala population showed significant genetic differences from the other three populations. The median joining network of cytochrome b gene suggested that the domestic Muscovy ducks present in India are the product of a single domestication event and probably introduced to India several years ago, as reported elsewhere. This study has also showed the suitability of nuclear DNA CYP2U1 gene in genetic diversity analysis.

Epidermal growth factor gene is a newly identified candidate gene for gout

Chromosome 4q25 has been identified as a genomic region associated with gout. However, the associations of gout with the genes in this region have not yet been confirmed. Here, we performed two-stage analysis to determine whether variations in candidate genes in the 4q25 region are associated with gout in a male Chinese Han population. We first evaluated 96 tag single nucleotide polymorphisms (SNPs) in eight inflammatory/immune pathway- or glucose/lipid metabolism-related genes in the 4q25 region in 480 male gout patients and 480 controls. The SNP rs12504538, located in the elongation of very-long-chain-fatty-acid-like family member 6 gene (Elovl6), was found to be associated with gout susceptibility (P_adjusted = 0.00595). In the second stage of analysis, we performed fine mapping analysis of 93 tag SNPs in Elovl6 and in the epidermal growth factor gene (EGF) and its flanking regions in 1017 male patients gout and 1897 healthy male controls. We observed a significant association between the T allele of EGF rs2298999 and gout (odds ratio = 0.77, 95% confidence interval = 0.67–0.88, P_adjusted = 6.42 × 10⁻³). These results provide the first evidence for an association between the EGF rs2298999 C/T polymorphism and gout. Our findings should be validated in additional populations.

Determination of major UDP-glucuronosyltransferase enzymes and their genotypes responsible for 20-HETE glucuronidation

The compound 20-HETE is involved in numerous physiological functions, including blood pressure and platelet aggregation. Glucuronidation of 20-HETE by UDP-glucuronosyltransferases (UGTs) is thought to be a primary pathway of 20-HETE elimination in humans. The present study identified major UGT enzymes responsible for 20-HETE glucuronidation and investigated their genetic influence on the glucuronidation reaction using human livers (n = 44). Twelve recombinant UGTs were screened to identify major contributors to 20-HETE glucuronidation. Based on these results, UGT2B7, UGT1A9, and UGT1A3 exhibited as major contributors to 20-HETE glucuronidation. The Km values of 20-HETE glucuronidation by UGT1A3, UGT1A9, and UGT2B7 were 78.4, 22.2, and 14.8 μM, respectively, while Vmax values were 1.33, 1.78, and 1.62 nmol/min/mg protein, respectively. Protein expression levels and genetic variants of UGT1A3, UGT1A9, and UGT2B7 were analyzed in human livers using Western blotting and genotyping, respectively. Glucuronidation of 20-HETE was significantly correlated with the protein levels of UGT2B7 (r(2) = 0.33, P < 0.001) and UGT1A9 (r(2) = 0.31, P < 0.001), but not UGT1A3 (r(2) = 0.02, P > 0.05). A correlation between genotype and 20-HETE glucuronidation revealed that UGT2B7 802C>T, UGT1A9 -118T9>T10, and UGT1A9 1399T>C significantly altered 20-HETE glucuronide formation (P < 0.05-0.001). Increased levels of 20-HETE comprise a risk factor for cardiovascular diseases, and the present data may increase our understanding of 20-HETE metabolism and cardiovascular complications.

Molecular aspects of hereditary spastic paraplegia

Hereditary spastic paraplegia (HSP) is a clinically and genetically heterogeneous group of neurodegenerative disorders characterized by progressive lower limbs spasticity and weakness. What was first thought to be a small group of rare Mendelian disorder has now become a large group that includes many complex syndromes. While large families with defined modes of inheritance were used for the initial HSP gene discovery, new sequencing technologies have recently allowed the study of small families, with the identification of many new disease causative genes. These discoveries are slowly leading to a better understanding of the molecular mechanisms underlying HSP with the identification of precise disease pathways. These insights may lead to new therapeutic strategies for what is a group of largely untreatable diseases. This review looks at the key players involved in HSP and where they act in their specific pathways.

Identification of cytochrome P450s involved in the metabolism of arachidonic acid in human platelets

Although cytochrome P450s (CYPs) have been identified in most human cells, identification of CYPs in human platelets remains poorly explored. CYP expressions in human platelets were screened by using reverse transcriptase-polymerase chain reaction and western blot analysis followed by functional assays using arachidonic acid (ARA). CYP1A1, 2U1, 2J2, 4A11, 4F2, and 5A1 were expressed as both proteins and mRNAs in platelets. Ethoxyresorufin-O-deethylase activity was observed in platelets and this activity was significantly decreased after treatment with the general P450 inhibitor SKF-525A and the CYP1A inhibitor, α-naphthoflavone (40-45%, P<0.001). Seventeen ARA metabolites were detected in ARA-treated platelets. Among these, the levels of 20-hydroxyeicosatetraenoic acid and epoxyeicosatrienoic acids were significantly decreased with the treatment of the P450 ω-hydroxylase inhibitor 17-octadecynoic acid (P<0.05-0.001). In summary, multiple ARA-metabolizing P450s were identified in human platelets. These findings may provide an important resource for understanding physiological function of platelet.

Expression of arachidonic acid-metabolizing cytochrome P450s in human megakaryocytic Dami cells

Cytochrome P450s (P450s) are involved in the metabolism of arachidonic acid (ARA), and ARA metabolites are associated with various cellular signaling pathways, such as blood hemostasis and inflammation. The present study demonstrates the expression of ARA-metabolizing P450s in the human megakaryocytic Dami cells using reverse transcriptase-polymerase chain reaction (RT-PCR) and immunublotting analysis followed by activity assays using ARA as a substrate. In addition to the previously identified CYP5A1, both protein and mRNAs of CYP1A1, 2U1, and 2J2 bands were detected. Ethoxyresorufin-O-deethylase (EROD) activity was observed in Dami cells, and its activity was significantly decreased after treatment with the P450 inhibitor SKF-525A when compared to the control groups (60% reduction, P < 0.001). CYP1A1 protein expression in Dami cells was induced by 3-methylenecholantheren. This increase in CYP1A1 protein level was correlated with enhanced EROD activity (fourfold increase vs. the control), as well as with increased metabolites, such as 20-hydroxyeicosatrienoic acid (20-HETE), 14, 15-EET (14-,15-epoxyeicosatrienoic acid), and 14, 15-dihydroxyeicosatrienoic acid (14, 15-DHET). The expression of soluble epoxide hydrolase, an enzyme responsible for the synthesis of DHETs from EETs, was confirmed by RT-PCR. Furthermore, 15 ARA metabolites, including 8,9-EET, 14,15-EET, and 20-HETE, were detected by LC-MS/MS in ARA-treated Dami cells, and their levels were decreased with the treatment of the SKF-525A. The present data suggest the possibility that the P450s play a role in the metabolism of ARA and other CYP-related substrates in human megakaryocytes and that P450 expression in megakaryocytic cell lines may predict their existences in platelets with functional activities.

Arachidonic acid ω-hydroxylase CYP4A11: inter-ethnic variations in the 8590T>C loss-of-function variant

The human Cytochrome P450 4A11 (CYP4A11) is a major ω-hydroxylase involved in the regulation of blood pressure in the kidney through the conversion of arachidonic acid into 20-hydroxyeicosatetraenoic acid (20-HETE). Previous studies have reported a significant association between the 8590T>C genetic variant of CYP4A11 and hypertension. Interestingly, several population-based studies have reported ethnic differences in the prevalence of hypertension, with the highest prevalence in African populations. The aim of this work was to determine the frequency and inter-ethnic comparison of the CYP4A11 (8590T>C) functional polymorphism, in five new ethnic groups: European (99 French Caucasians), African (36 Gabonese and 50 Senegalese), South American (60 Peruvians) and North African (53 Tunisians) populations, using polymerase chain reaction-single strand conformational polymorphism and sequencing strategies. We confirmed that the CYP4A11 (8590T>C) functional polymorphism exhibits inter-ethnic frequency differences. Noteworthy, the highest 8590C allele frequency was observed in the Tunisian (30.2%), followed by Senegalese (20%) populations. In addition, the CC genotype was only found in the Gabonese and Tunisian populations (5.6% and 8.4%, respectively). These populations may be of major interest to help to clarify the linkage between hypertension and CYP4A11 (8590T>C) genotype in African populations. These findings provide data for further studies that investigate the potential association of CYP4A11 (8590T>C) variant with an incidence of hypertension genesis in respect of ethnicity.

Genetic polymorphisms of glycine N-acyltransferase (GLYAT) in a French Caucasian population

In humans, the glycine N-acyltransferase enzyme (GLYAT) is thought to be important in the detoxification of endogenous and xenobiotic compounds which contain a carboxylic acid group, such as benzoic, isovaleric, or acetylsalicylic acids. The aim of this work was to report a comprehensive investigation of GLYAT genetic polymorphisms in DNA samples from 55 subjects of French Caucasian origin, using polymerase chain reaction-single-strand conformation polymorphism and sequencing strategies. Seven different polymorphisms of the GLYAT gene were identified, including two polymorphisms in the 5' flanking region of the gene (g.-8457C>T and g.-8010A>G), two polymorphisms in intron 5 (g.13931A>G and g.13944C>T) and three missense mutations in exon 2 (g.49T>A; p.Ser17Thr), exon 5 (g.13886A>G; p.Asn156Ser) and exon 6 (g.14435C>T; p.Arg199Cys). In addition to the wild-type allele GLYAT*1 (2.7%), four novel alleles were identified: GLYAT*2A (75.5%), *2B (4.5%), *3 (16.4%) and *4 (0.9%), and five different genotypes. Localisation of the p.Ser17Thr and p.Arg199Cys missense mutations in predicted secondary structures suggest that these variants might have a potential role on the GLYAT protein activity. These results could be helpful in investigating the potential association of GLYAT variants with an incidence of reduced efficiency in xenobiotic carboxylic acids detoxification in humans.