Genetic variation in aldo-keto reductase 1D1 (AKR1D1) affects the expression and activity of multiple cytochrome P450s

Cecum microbiome and metabolism characteristics of Silky Fowl and White Leghorn chicken in late laying stages

Cecal microflora plays a key role in the production performance and immune function of chickens. White Leghorn (WL) is a well-known commercial layer line chicken with high egg production rate. In contrast, Silky Fowl (SF), a Chinese native chicken variety, has a low egg production rate, but good immune performance. This study analyzed the composition of cecal microbiota, metabolism, and gene expression in intestinal tissue of these varieties and the correlations among them. Significant differences were observed in the cecal microbes: Bacteroides was significantly enriched in WL, whereas Veillonellaceae and Parabacteroides were significantly enriched in SF. Carbohydrate biosynthesis and metabolism pathways were significantly upregulated in WL cecum, which might provide more energy to the host, leading to persistently high levels of egg production. The higher Parabacteroides abundance in SF increased volicitin content, enhanced α-linolenic acid metabolism, and significantly negatively correlated with metabolites of propanoate metabolism and carbohydrate metabolism. Genes related to lipid metabolism, immunity, and melanogenesis were significantly upregulated in the SF cecum, regulating lipid metabolism, and participating in the immune response, while genes related to glucose metabolism and bile acid metabolism were expressed at higher levels in WL, benefiting energy support. This study provided a mechanism for intestinal microorganisms and metabolic pathways to regulate chicken egg-laying performance and immunity.

Genome-wide association study reveals candidate genes relevant to body weight in female turkeys (Meleagris gallopavo)

The underlying genetic mechanisms affecting turkey growth traits have not been widely investigated. Genome-wide association studies (GWAS) is a powerful approach to identify candidate regions associated with complex phenotypes and diseases in livestock. In the present study, we performed GWAS to identify regions associated with 18-week body weight in a turkey population. The data included body weight observations for 24,989 female turkeys genotyped based on a 65K SNP panel. The analysis was carried out using a univariate mixed linear model with hatch-week-year and the 2 top principal components fitted as fixed effects and the accumulated polygenic effect of all markers captured by the genomic relationship matrix as random. Thirty-three significant markers were observed on 1, 2, 3, 4, 7 and 12 chromosomes, while 26 showed strong linkage disequilibrium extending up to 410 kb. These significant markers were mapped to 37 genes, of which 13 were novel. Interestingly, many of the investigated genes are known to be involved in growth and body weight. For instance, genes AKR1D1, PARP12, BOC, NCOA1, ADCY3 and CHCHD7 regulate growth, body weight, metabolism, digestion, bile acid biosynthetic and development of muscle cells. In summary, the results of our study revealed novel candidate genomic regions and candidate genes that could be managed within a turkey breeding program and adapted in fine mapping of quantitative trait loci to enhance genetic improvement in this species.

Mechanistic Insights into Side Effects of Troglitazone and Rosiglitazone Using a Novel Inverse Molecular Docking Protocol

Thiazolidinediones form drugs that treat insulin resistance in type 2 diabetes mellitus. Troglitazone represents the first drug from this family, which was removed from use by the FDA due to its hepatotoxicity. As an alternative, rosiglitazone was developed, but it was under the careful watch of FDA for a long time due to suspicion, that it causes cardiovascular diseases, such as heart failure and stroke. We applied a novel inverse molecular docking protocol to discern the potential protein targets of both drugs. Troglitazone and rosiglitazone were docked into predicted binding sites of >67,000 protein structures from the Protein Data Bank and examined. Several new potential protein targets with successfully docked troglitazone and rosiglitazone were identified. The focus was devoted to human proteins so that existing or new potential side effects could be explained or proposed. Certain targets of troglitazone such as 3-oxo-5-beta-steroid 4-dehydrogenase, neutrophil collagenase, stromelysin-1, and VLCAD were pinpointed, which could explain its hepatoxicity, with additional ones indicating that its application could lead to the treatment/development of cancer. Results for rosiglitazone discerned its interaction with members of the matrix metalloproteinase family, which could lead to cancer and neurodegenerative disorders. The concerning cardiovascular side effects of rosiglitazone could also be explained. We firmly believe that our results deepen the mechanistic understanding of the side effects of both drugs, and potentially with further development and research maybe even help to minimize them. On the other hand, the novel inverse molecular docking protocol on the other hand carries the potential to develop into a standard tool to predict possible cross-interactions of drug candidates potentially leading to adverse side effects.

Quantitative analysis of mRNA expression levels of aldo-keto reductase and short-chain dehydrogenase/reductase isoforms in human livers

The aldo-keto reductase (AKR) and short-chain dehydrogenase/reductase (SDR) superfamilies are responsible for the reduction in compounds containing the aldehyde, ketone, and quinone groups. In humans, 12 AKR isoforms (AKR1A1, AKR1B1, AKR1B10, AKR1B15, AKR1C1, AKR1C2, AKR1C3, AKR1C4, AKR1D1, AKR1E2, AKR7A2, and AKR7A3) and 6 SDR isoforms (CBR1, CBR3, CBR4, HSD11B1, DHRS4, and DCXR) have been found to catalyze the reduction in xenobiotics, but their hepatic expression levels are unclear. The purpose of this study is to determine the absolute mRNA expression levels of these 18 isoforms in the human liver. In 22 human livers, all isoforms, except for AKR1B15, are expressed, and AKR1C2 (on average 1.6 × 10⁶ copy/μg total RNA), AKR1C3 (1.3 × 10⁶), AKR1C1 (1.3 × 10⁶), CBR1 (9.7 × 10⁵), and HSD11B1 (1.1 × 10⁶) are abundant, representing 67% of the total expression of reductases in the liver. The expression levels of AKR1C2, AKR1C3, AKR1C1, CBR1, and HSD11B1 are significantly correlated with each other, except between AKR1C2 and CBR1, suggesting that they might be regulated by common factor(s). In conclusion, this study comprehensively determined the absolute expression of mRNA expression of each AKR and SDR isoform in the human liver.

Glucocorticoids regulate AKR1D1 activity in human liver in vitro and in vivo

Steroid 5β-reductase (AKR1D1) is highly expressed in human liver where it inactivates endogenous glucocorticoids and catalyses an important step in bile acid synthesis. Endogenous and synthetic glucocorticoids are potent regulators of metabolic phenotype and play a crucial role in hepatic glucose metabolism. However, the potential of synthetic glucocorticoids to be metabolised by AKR1D1 as well as to regulate its expression and activity has not been investigated. The impact of glucocorticoids on AKR1D1 activity was assessed in human liver HepG2 and Huh7 cells; AKR1D1 expression was assessed by qPCR and Western blotting. Genetic manipulation of AKR1D1 expression was conducted in HepG2 and Huh7 cells and metabolic assessments were made using qPCR. Urinary steroid metabolite profiling in healthy volunteers was performed pre- and post-dexamethasone treatment, using gas chromatography-mass spectrometry. AKR1D1 metabolised endogenous cortisol, but cleared prednisolone and dexamethasone less efficiently. In vitro and in vivo, dexamethasone decreased AKR1D1 expression and activity, further limiting glucocorticoid clearance and augmenting action. Dexamethasone enhanced gluconeogenic and glycogen synthesis gene expression in liver cell models and these changes were mirrored by genetic knockdown of AKR1D1 expression. The effects of AKR1D1 knockdown were mediated through multiple nuclear hormone receptors, including the glucocorticoid, pregnane X and farnesoid X receptors. Glucocorticoids down-regulate AKR1D1 expression and activity and thereby reduce glucocorticoid clearance. In addition, AKR1D1 down-regulation alters the activation of multiple nuclear hormone receptors to drive changes in gluconeogenic and glycogen synthesis gene expression profiles, which may exacerbate the adverse impact of exogenous glucocorticoids.

The influence of biological sex and sex hormones on bile acid synthesis and cholesterol homeostasis

Obesity and elevated serum lipids are associated with a threefold increase in the risk of developing atherosclerosis, a condition that underlies stroke, myocardial infarction, and sudden cardiac death. Strategies that aim to reduce serum cholesterol through modulation of liver enzymes have been successful in decreasing the risk of developing atherosclerosis and reducing mortality. Statins, which inhibit cholesterol biosynthesis in the liver, are considered among the most successful compounds developed for the treatment of cardiovascular disease. However, recent debate surrounding their effectiveness and safety prompts consideration of alternative cholesterol-lowering therapies, including increasing cholesterol catabolism through bile acid (BA) synthesis. Targeting the enzymes that convert cholesterol to BAs represents a promising alternative to other cholesterol-lowering approaches that treat atherosclerosis as well as fatty liver diseases and diabetes mellitus. Compounds that modify the activity of these pathways have been developed; however, there remains a lack of consideration of biological sex. This is necessary in light of strong evidence for sexual dimorphisms not only in the incidence and progression of the diseases they influence but also in the expression and activity of the proteins affected and in the manner in which men and women respond to drugs that modify lipid handling in the liver. A thorough understanding of the enzymes involved in cholesterol catabolism and modulation by biological sex is necessary to maximize their therapeutic potential.

The AKR1D1*36 (rs1872930) Allelic Variant Is Independently Associated With Clopidogrel Treatment Outcome

AIMS

The present observational cohort study evaluated the association between the AKR1D1*36 (rs1872930) allele and the risk of major adverse cardiovascular and cerebrovascular events (MACCE) in clopidogrel treated patients.

METHODS

We screened 198 consecutive cardiovascular patients on clopidogrel therapy admitted in October to November 2010 with cardiovascular or cerebrovascular symptoms; of these 118 met the study protocol entry criteria; the median age of the cohort was 62.5 years (IQR 57-66 years), and 55% were females.

RESULTS

The median follow up time was 38.5 (IQR 24-48) months; Kaplan-Meier/Log-rank analysis showed that patients carrying the AKR1D1*36 allelic variant have a shorter event-free-survival compared to wild type patients, hazard ratio = 2.193 (95% CI, 1.091 to 4.406); p = 0.0155. Multivariable Cox regression analysis confirmed the AKR1D1*36 allele as an independent risk factor (HR = 2.36; 95% CI, 1.34 to 4.18) and identified 3 other risk factors for MACCE; previous percutaneous interventions (PCI), HR = 2.78; (95% CI, 1.34 to 5.78), and a history of myocardial infarction, HR = 2.62; (95% CI, 1.48 to 4.64) at baseline and the previously reported CYP2C19*2 polymorphism (HR = 2.33; 95% CI, 1.33 to 4.06).

CONCLUSION

The AKR1D1*36 (rs1872930) variant is independently associated with a higher risk for MACCE and shorter event-free survival time.

AKR1D1*36 C>T (rs1872930) allelic variant is associated with variability of the CYP2C9 genotype predicted pharmacokinetics of ibuprofen enantiomers - a pilot study in healthy volunteers

The relative contribution of CYP2C9 allelic variants to the pharmacokinetics (PK) of ibuprofen (IBP) enantiomers has been studied extensively, but the potential clinical benefit of pharmacogenetically guided IBP treatment is not evident yet. The role of AKR1D1*36C>T (rs 1872930) allelic variant in interindividual variability of CYP450 mediated drug metabolism was recently elucidated. A total of 27 healthy male subjects, volunteers in IBP single-dose two-way cross-over bioequivalence studies were genotyped for CYP2C9*2, CYP2C9*3 and AKR1D1*36 polymorphisms. The correlation between CYP2C9 and AKR1D1 genetic profile and the PK parameters for S-(+) and R-(-)-IBP was evaluated. Remarkable changes in the PK values pointing to reduced CYP2C9 enzyme activity were detected only in the CYP2C9*2 allelic variant carriers. Statistically significant association between the AKR1D1*36 allele and the increased IBP metabolism (low AUC0-t and 0-∞, high Cltot and short tmax values for both enantiomers) was observed in subjects carrying the CYP2C9 *1/*3 or CYP2C9*1/*1 genotype. The clinical value of concomitant CYP2C9 and AKR1D1 genotyping has to be further verified.

Genome sequencing analysis of blood cells identifies germline haplotypes strongly associated with drug resistance in osteosarcoma patients

BACKGROUND

Osteosarcoma is the most common malignant bone tumor in children. Survival remains poor among histologically poor responders, and there is a need to identify them at diagnosis to avoid delivering ineffective therapy. Genetic variation contributes to a wide range of response and toxicity related to chemotherapy. The aim of this study is to use sequencing of blood cells to identify germline haplotypes strongly associated with drug resistance in osteosarcoma patients.

METHODS

We used sequencing data from two patient datasets, from Inova Hospital and the NCI TARGET. We explored the effect of mutation hotspots, in the form of haplotypes, associated with relapse outcome. We then mapped the single nucleotide polymorphisms (SNPs) in these haplotypes to genes and pathways. We also performed a targeted analysis of mutations in Drug Metabolizing Enzymes and Transporter (DMET) genes associated with tumor necrosis and survival.

RESULTS

We found intronic and intergenic hotspot regions from 26 genes common to both the TARGET and INOVA datasets significantly associated with relapse outcome. Among significant results were mutations in genes belonging to AKR enzyme family, cell-cell adhesion biological process and the PI3K pathways; as well as variants in SLC22 family associated with both tumor necrosis and overall survival. The SNPs from our results were confirmed using Sanger sequencing. Our results included known as well as novel SNPs and haplotypes in genes associated with drug resistance.

CONCLUSION

We show that combining next generation sequencing data from multiple datasets and defined clinical data can better identify relevant pathway associations and clinically actionable variants, as well as provide insights into drug response mechanisms.

Allele frequency and genotype distribution of aldo keto reductase 1D1 (AKR1D1) rs1872930 genetic variant in a Macedonian population

AKR1D1 has been recently recognized as novel candidate gene regulating CYP450 coexpression. The AKR1D1 3’-UTR SNP (rs1872930) functional genetic variation was considered as important source of phenotypic variation in the activity of CYP450. The minor allele frequency of rs1872930 varies among populations suggesting population specific pharmacogenetic implications. The aim of the present study was to determine the AKR1D1 (rs1872930) minor allele frequencies and genotypes distribution in a population from Republic of Macedonia. Four hundred and fifty unrelated Macedonian subjects were studied. AKR1D1 variant allele was detected by real time polymerase chain reaction. The distribution of AKR1D1*1/*1, *1/*36 and *36/*36 genotypes was 0.569, 0.360 and 0.071, respectively. The overall frequency of AKR1D1*36 variant allele was 0.251. Our study is the first to assess the frequency distribution of the AKR1D1*36 (rs1872930) genetic variant in population from Republic of Macedonia and the findings are in accordance with the frequency reported for the Caucasian population. Keywords: AKR1D1, rs1872930, CYP450, variability, gene expression, enzyme activity

Variation in CYP2A6 Activity and Personalized Medicine

The cytochrome P450 2A6 (CYP2A6) enzyme metabolizes several clinically relevant substrates, including nicotine-the primary psychoactive component in cigarette smoke. The gene that encodes the CYP2A6 enzyme is highly polymorphic, resulting in extensive interindividual variation in CYP2A6 enzyme activity and the rate of metabolism of nicotine and other CYP2A6 substrates including cotinine, tegafur, letrozole, efavirenz, valproic acid, pilocarpine, artemisinin, artesunate, SM-12502, caffeine, and tyrosol. CYP2A6 expression and activity are also impacted by non-genetic factors, including induction or inhibition by pharmacological, endogenous, and dietary substances, as well as age-related changes, or interactions with other hepatic enzymes, co-enzymes, and co-factors. As variation in CYP2A6 activity is associated with smoking behavior, smoking cessation, tobacco-related lung cancer risk, and with altered metabolism and resulting clinical responses for several therapeutics, CYP2A6 expression and enzyme activity is an important clinical consideration. This review will discuss sources of variation in CYP2A6 enzyme activity, with a focus on the impact of CYP2A6 genetic variation on metabolism of the CYP2A6 substrates.

Predictors of Variation in CYP2A6 mRNA, Protein, and Enzyme Activity in a Human Liver Bank: Influence of Genetic and Nongenetic Factors

Cytochrome P450 2A6 CYP2A6: metabolizes several clinically relevant substrates, including nicotine, the primary psychoactive component in cigarette smoke. Smokers vary widely in their rate of inactivation and clearance of nicotine, altering numerous smoking phenotypes. We aimed to characterize independent and shared impact of genetic and nongenetic sources of variation in CYP2A6 mRNA, protein, and enzyme activity in a human liver bank (n = 360). For the assessment of genetic factors, we quantified levels of CYP2A6, cytochrome P450 oxidoreductase (POR), and aldo-keto reductase 1D1 (AKR1D1) mRNA, and CYP2A6 and POR proteins. CYP2A6 enzyme activity was determined through measurement of cotinine formation from nicotine and 7-hydroxycoumarin formation from coumarin. Donor DNA was genotyped for CYP2A6, POR, and AKR1D1 genetic variants. Nongenetic factors assessed included gender, age, and liver disease. CYP2A6 phenotype measures were positively correlated to each other (r values ranging from 0.47-0.88, P < 0.001). Female donors exhibited higher CYP2A6 mRNA expression relative to males (P < 0.05). Donor age was weakly positively correlated with CYP2A6 protein (r = 0.12, P < 0.05) and activity (r = 0.20, P < 0.001). CYP2A6 reduced-function genotypes, but not POR or AKR1D1 genotypes, were associated with lower CYP2A6 protein (P < 0.001) and activity (P < 0.01). AKR1D1 mRNA was correlated with CYP2A6 mRNA (r = 0.57, P < 0.001), protein (r = 0.30, P < 0.001), and activity (r = 0.34, P < 0.001). POR protein was correlated with CYP2A6 activity (r = 0.45, P < 0.001). Through regression analyses, we accounted for 17% (P < 0.001), 37% (P < 0.001), and 77% (P < 0.001) of the variation in CYP2A6 mRNA, protein, and activity, respectively. Overall, several independent and shared sources of variation in CYP2A6 activity in vitro have been identified, which could translate to variable hepatic clearance of nicotine.

Prediction of Intravenous Busulfan Clearance by Endogenous Plasma Biomarkers Using Global Pharmacometabolomics

INTRODUCTION

High-dose busulfan (busulfan) is an integral part of the majority of hematopoietic cell transplantation conditioning regimens. Intravenous (IV) busulfan doses are personalized using pharmacokinetics (PK)-based dosing where the patient's IV busulfan clearance is calculated after the first dose and is used to personalize subsequent doses to a target plasma exposure. PK-guided dosing has improved patient outcomes and is clinically accepted but highly resource intensive.

OBJECTIVE

We sought to discover endogenous plasma biomarkers predictive of IV busulfan clearance using a global pharmacometabolomics-based approach.

METHODS

Using LC-QTOF, we analyzed 59 (discovery) and 88 (validation) plasma samples obtained before IV busulfan administration.

RESULTS

In the discovery dataset, we evaluated the association of the relative abundance of 1885 ions with IV busulfan clearance and found 21 ions that were associated with IV busulfan clearance tertiles (r² ≥ 0.3). Identified compounds were deoxycholic acid and/or chenodeoxycholic acid, and linoleic acid. We used these 21 ions to develop a parsimonious seven-ion linear predictive model that accurately predicted IV busulfan clearance in 93% (discovery) and 78% (validation) of samples.

CONCLUSION

IV busulfan clearance was significantly correlated with the relative abundance of 21 ions, seven of which were included in a predictive model that accurately predicted IV busulfan clearance in the majority of the validation samples. These results reinforce the potential of pharmacometabolomics as a critical tool in personalized medicine, with the potential to improve the personalized dosing of drugs with a narrow therapeutic index such as busulfan.

Genomewide comparison of the inducible transcriptomes of nuclear receptors CAR, PXR and PPARα in primary human hepatocytes

The ligand-activated nuclear receptor pregnane X receptor (PXR, NR1I2) and the constitutive androstane receptor (CAR, NR1I3) are two master transcriptional regulators of many important drug metabolizing enzymes and transporter genes (DMET) in response to xenobiotics including many drugs. The peroxisome proliferator-activated receptor alpha (PPARα, NR1C1), the target of lipid lowering fibrate drugs, primarily regulates fatty acid catabolism and energy-homeostasis. Recent research has shown that there are substantial overlaps in the regulated genes of these receptors. For example, both CAR and PXR also modulate the transcription of key enzymes involved in lipid and glucose metabolism and PPARα also functions as a direct transcriptional regulator of important DMET genes including cytochrome P450s CYP3A4 and CYP2C8. Despite their important and widespread influence on liver metabolism, comparative data are scarce, particularly at a global level and in humans. The major objective of this study was to directly compare the genome-wide transcriptional changes elucidated by the activation of these three nuclear receptors in primary human hepatocytes. Cultures from six individual donors were treated with the prototypical ligands for CAR (CITCO), PXR (rifampicin) and PPARα (WY14,643) or DMSO as vehicle control. Genomewide mRNA profiles determined with Affymetrix microarrays were analyzed for differentially expressed genes and metabolic functions. The results confirmed known prototype target genes and revealed strongly overlapping sets of coregulated but also distinctly regulated and novel responsive genes and pathways. The results further specify the role of PPARα as a regulator of drug metabolism and the role of the xenosensors PXR and CAR in lipid metabolism and energy homeostasis. This article is part of a Special Issue entitled: Xenobiotic nuclear receptors: New Tricks for An Old Dog, edited by Dr. Wen Xie.

Interindividual Variability in Cytochrome P450-Mediated Drug Metabolism

The cytochrome P450 (P450) enzymes are the predominant enzyme system involved in human drug metabolism. Alterations in the expression and/or activity of these enzymes result in changes in pharmacokinetics (and consequently the pharmacodynamics) of drugs that are metabolized by this set of enzymes. Apart from changes in activity as a result of drug-drug interactions (by P450 induction or inhibition), the P450 enzymes can exhibit substantial interindividual variation in basal expression and/or activity, leading to differences in the rates of drug elimination and response. This interindividual variation can result from a myriad of factors, including genetic variation in the promoter or coding regions, variation in transcriptional regulators, alterations in microRNA that affect P450 expression, and ontogenic changes due to exposure to xenobiotics during the developmental and early postnatal periods. Other than administering a probe drug or cocktail of drugs to obtain the phenotype or conducting a genetic analysis to determine genotype, methods to determine interindividual variation are limited. Phenotyping via a probe drug requires exposure to a xenobiotic, and genotyping is not always well correlated with phenotype, making both methodologies less than ideal. This article describes recent work evaluating the effect of some of these factors on interindividual variation in human P450-mediated metabolism and the potential utility of endogenous probe compounds to assess rates of drug metabolism among individuals.

Comparing gene expression during cadmium uptake and distribution: untreated versus oral Cd-treated wild-type and ZIP14 knockout mice

The nonessential metal cadmium (Cd) is toxic only after entering the cell. Proteins possibly relevant to intracellular Cd accumulation include the divalent metal transporter-1 (DMT1) and all 14 zinc-like iron-like protein (ZIP) importers, 10 zinc transporter (ZnT) exporters, and metallothionein chaperones MT1 and MT2. Comparing oral Cd-treated ZIP14 knockout (KO) with wild-type (WT) mice, we predicted Cd uptake and distribution would be diminished in the KO-because ZIP14 is very highly expressed in GI tract and liver; this was indeed observed for Cd content in liver. However, the reverse was found in kidney and lung from 6 or 12 h through 10 days of Cd exposure; at these times, Cd accumulation was unexpectedly greater in KO than WT mice; mRNA levels of the 27 above-mentioned genes were thus examined in proximal small intestine (PSI) versus kidney to see if these paradoxical effects could be explained by substantial alterations in any of the other 26 genes. PSI genes highly expressed in untreated WT animals included seven ZIP and five ZnT transporters, DMT1, MT1, and MT2; kidney genes included 11 ZIP and 7 ZnT transporters, DMT1, MT1, and MT2. Over 10 days of oral Cd, a bimodal response was seen for Cd content in PSI and for various mRNAs; initially, acute effects caused by the toxic metal; subsequently, the up- or down-regulation of important genes presumably to combat the sustained adversity. These data underscore the complex interplay between the gastrointestinal tract and renal proteins that might be relevant to Cd uptake and distribution in animals exposed to oral Cd.

Biologically active substances-enriched diet regulates gonadotrope cell activation pathway in liver of adult and old rats

According to the Hippocrates’ theorem “Let food be your medicine and medicine be your food”, dietary interventions may induce changes in the metabolic and inflammatory state by modulating the expression of important genes involved in the chronic disorders. The aim of the present study was to evaluate the influence of long-term (14 months) use of biologically active substances-enriched diet (BASE-diet) on transcriptomic profile of rats’ liver. The experiment was conducted on 36 Sprague–Dawley rats divided into two experimental groups (fed with control or BASE-diet, both n = 18). Control diet was a semi-synthetic diet formulated according to the nutritional requirements for laboratory animals. The BASE-diet was enriched with a mixture of polyphenolic compounds, β-carotene, probiotics, and n-3 and n-6 polyunsaturated fatty acids. In total, n = 3,017 differentially expressed (DE) genes were identified, including n = 218 DE genes between control and BASE groups after 3 months of feeding and n = 1,262 after 14 months. BASE-diet influenced the expression of genes involved particularly in the gonadotrope cell activation pathway and guanylate cyclase pathway, as well as in mast cell activation, gap junction regulation, melanogenesis and apoptosis. Especially genes involved in regulation of GnRH were strongly affected by BASE-diet. This effect was stronger with the age of animals and the length of diet use. It may suggest a link between the diet, reproductive system function and aging.

microRNA-34a is associated with expression of key hepatic transcription factors and cytochromes P450

microRNA (miRNA) mediated regulation of gene expression has emerged as a significant mechanism contributing to variation in gene expression. In this study, we evaluated the potential role of miRNAs in regulating expression of hepatic cytochromes P450 and their transcriptional regulatory genes. We screened the Targetscan database for high scoring miRNA binding site predictions in selected hepatic DMEs and transcription factors. Expression profiling for candidate miRNAs (n=22) and their target genes (n=20) was performed in 50 human liver samples (25 female, 25 male). Significant negative correlations were observed between expression levels of several CYPs/hepatic transcription factors and the hepatic miRNAs studied. Interestingly, hepatic miR-34a demonstrated significant negative correlation with expression levels of multiple hepatic transcription factors (including NR1I2 and HNF4α) and DMEs (CYP3A4, CYP2C19). miR-34a expression was also significantly higher in males than in females in congruence with previous observations of higher CYP3A4 expression in females versus males. A mediation analysis revealed that miR-34a was involved in significant mediation of the association observed between CYP2C19 and several hepatic transcription factors (HNF4α, NR1I2). miR-34a may thus play a key regulatory role and be a key contributory factor to the inter-individual variability observed in expression of key drug metabolizing genes in humans.