Marwa Tantawy
Danxin Wang- Center for Pharmacogenomics, Department of Pharmacotherapy and Translational Research, College of Pharmacy, University of Florida, Gainesville, FL, United States
Cytochrome P450 (CYP) drug metabolizing enzymes are responsible for the metabolism of over 70% of currently used medications with the CYP3A family being the most important CYP enzymes in the liver. Large inter-person variability in expression/activity of the CYP3As greatly affects drug exposure and treatment outcomes, yet the cause of such variability remains elusive. Micro-RNAs (miRNAs) are small noncoding RNAs that negatively regulate gene expression and are involved in diverse cellular processes including metabolism of xenobiotics and therapeutic outcomes. Target prediction and in vitro functional assays have linked several miRNAs to the control of CYP3A4 expression. Yet, their co-expression with CYP3As in the liver remain unclear. In this study, we used genome-wide miRNA profiling in liver samples to identify miRNAs associated with the expression of the CYP3As. We identified and validated both miR-107 and miR-1260 as strongly associated with the expression of CYP3A4, CYP3A5, and CYP3A43. Moreover, we found associations between miR-107 and nine transcription factors (TFs) that regulate CYP3A expression, with estrogen receptor alpha (ESR1) having the largest effect size. Including ESR1 and the other TFs in the regression model either diminished or abolished the associations between miR-107 and the CYP3As, indicating that the role of miR-107 in CYP3A expression may be indirect and occur through these key TFs. Indeed, testing the other nine CYPs previously shown to be regulated by ESR1 identified similar miR-107 associations that were dependent on the exclusion of ESR1 and other key TFs in the regression model. In addition, we found significant differences in miRNA expression profiles in liver samples between race and sex. Together, our results identify miR-107 as a potential epigenetic regulator that is strongly associated with the expression of many CYPs, likely via impacting the CYP regulatory network controlled by ESR1 and other key TFs. Therefore, both genetic and epigenetic factors that alter the expression of miR-107 may have a broad influence on drug metabolism.
Introduction
The cytochrome P450s (CYPs) are comprised of a superfamily of enzymes involved in the biotransformation of xenobiotics and endogenous substances. Members of the CYP1, CYP2, and CYP3 families are the main drug metabolizing enzymes in the liver, responsible for ∼70% of all clinically used medications (Wienkers and Heath, 2005). Among them, the CYP3A subfamily is particularly important in adults because of their broad substrate specificity and prominent expression in liver (Danielson, 2002), where many medications are metabolized. The four CYP3A genes (CYP3A4, 3A5, 3A7, and 3A43) are encoded in a cluster located on chromosome 7q22.1. CYP3A4 is mainly expressed in the adult liver and is responsible for metabolizing 30%–40% of currently used drugs (Danielson, 2002). CYP3A5 shares many substrates with CYP3A4 and is also expressed outside the liver, but due to the frequent loss of function genetic variant CYP3A5*3, many individuals do not express CYP3A5 (Kuehl et al., 2001). CYP3A7 is primarily expressed in fetal livers with low expression levels observed in adults, while CYP3A43 is expressed in liver and extra-hepatic tissues like testis and prostate, but its substrates remain unclear. Overall, CYP3A-mediated drug metabolism is highly variable between individuals, greatly affecting drug exposure and treatment outcomes.
Many factors are known to influence the expression and function of the CYP3A enzymes. The genetic variant CYP3A4*22 (frequency 4%–8% in Whites and <1% in Blacks) reduces CYP3A4 expression 2–6-fold via affecting alternative splicing and transcription (Wang et al., 2011; Wang and Sadee, 2016; Collins and Wang, 2020) and is considered the most clinically relevant genetic variant in CYP3A4, associating with many phenotypes related to CYP3A4 metabolism (Mulder et al., 2021). Besides cis-acting genetic polymorphisms, the expression of the CYP3As is also regulated by trans-acting transcription factors (TFs), epigenetic modifications, and non-genetic factors (Zanger and Schwab, 2013). At the transcriptional level, the expression of many liver enriched TFs is correlated with the expression of CYP3As and other CYPs (Yang et al., 2010). The most significantly correlated TFs include: NR1I3 (constitutive androstane receptor or CAR), RXR (retinoid X receptor), NR1I2 (pregnane X receptor or PXR), HNF4A (hepatocyte nuclear factor 4 alpha), PPARA (peroxisome proliferator activated receptor alpha), FOXA2 (forkhead box A2), AHR (aryl hydrocarbon receptor) and ARNT (aryl hydrocarbon receptor nuclear translocator), as well as our recently identified master regulator ESR1 (estrogen receptor alpha) (Wang et al., 2019; Collins and Wang, 2021a). Moreover, ESR1 appears to be the top regulator for the expression of the CYP3As and other CYP enzymes, explaining up to 63% of the variability in expression of CYP3A4 and the other CYPs (Collins and Wang, 2021a). The activity of these TFs is known to be influenced by endogenous ligands and xenobiotic inducers, and it is likely that underlying genetic and/or epigenetic factors that alter their expression will also contribute to variation in CYP3A expression (Zanger and Schwab, 2013).
Micro-RNAs (miRNAs) are small noncoding RNAs of 21–25 nucleotides (nts) that negatively regulate gene expression. This effect is typically mediated by their binding to the 3′-untranslated regions (3′-UTRs) of target mRNAs, thereby interfering with mRNA translation or leading to degradation of the mRNA. There is also evidence for transcriptional control of gene expression via miRNA-targeting (Catalanotto et al., 2016). Through these mechanisms, miRNAs play an essential role in gene regulation and have been well-documented for a variety of biological processes, including: differentiation, development, proliferation, apoptosis, and necrosis (Londin et al., 2015; Yang et al., 2017). miRNAs also regulate the expression of genes related to pharmacodynamics and pharmacokinetics (Rieger et al., 2013), thereby altering drug toxicity (Wang et al., 2009; Szkolnicka et al., 2016) and therapeutic outcomes (Mishra et al., 2007; Sharma et al., 2020). Moreover, the expression of miRNAs is subject to change under drug treatment (Gufford et al., 2018) or environmental chemical exposure (Hou et al., 2011), which in turn regulates the expression of the CYP enzymes (Li et al., 2019). The interplay between dysregulation of miRNAs, changes in CYP expression, and CYP-dependent bioactivation is considered to be the underlying cause for environmental toxicology and carcinogenesis (Kalscheuer et al., 2008; Li et al., 2019). However, the regulatory roles of miRNAs in regulating constitutive expression of the CYP enzymes remain unclear.
A few studies have reported miRNA-mediated regulation of the CYP enzymes, including CYP3A4, by either direct 3′UTR binding or through indirect targeting of key TFs controlling CYP expression (Pan et al., 2009; Li et al., 2019). These previous studies focused on several candidate miRNAs that were predicted to bind CYP3A4 or its regulators and used in vitro cell models to evaluate the function of the miRNAs on their targets (Li et al., 2019). However, the relationship between the expression levels of these miRNAs and CYP3A4 has not been tested in the liver, with the exception of a single study illustrating the relationship between four miRNAs in a small (27 sample) liver cohort from the Chinese Han population (Wei et al., 2014). Therefore, the role of miRNA-mediated regulation of the CYP3A genes in human livers warrants further investigation. In this study, we applied an untargeted genome-wide approach to profile miRNA expression in a large cohort of human liver samples from both European and African American donors to search for miRNAs that may regulate the expression of the CYP3As.
Materials and methods
Human liver samples
Human liver samples from African American (AA) (n = 104) and European American (EA) (n = 127) donors were obtained from The Cooperative Human Tissue Network (CHTN). Demographics of liver donors are in Supplementary Table S1. The median ages are AA—56 years (range 0–97) and EA—61 years (range 14–83) and the percentage of males for each are AA—44% and EA—46%. The University of Florida Institutional Review Board approved the human tissue study.
RNA preparation and quantitative gene expression analysis
Total RNA was extracted from the liver tissue samples using the Direct-zol RNA miniprep plus kit (ZYMO Research, United States) followed by cDNA synthesis as described previously (Collins and Wang, 2021a). mRNA levels of CYP3As, CYP2Cs, CYP1As, ESR1, NR1I3, NR1I2, RXR, HNF4A, PPARA, FOXA2, AHR, and ARNT, and the internal control β-actin were measured with real-time PCR using TaqMan assays (Thermo Fisher Scientific) or SYBR Green with gene-specific primers as described (Collins and Wang, 2021a). The relative expression of each gene was calculated using the following formula: expression level of tested gene = antilog2 (mean Ct value of internal control—Ct value of tested gene) *106. After Log10 transformation, the expression data of all genes tested followed normal distribution as described previously (Collins and Wang, 2021a). The expression profiles of 12 CYPs and 9 TFs, the relationship between each CYP, the association between CYP and TFs, and the racial differences in expression between AA and EA have been described previously (Collins and Wang, 2021a).
miRNA profiling in the discovery phase
We used the TaqMan OpenArray Human MicroRNA panels and QuantStudio 12K Flex (Applied Biosystem, Thermo Fisher Scientific, CA, United States) for initial genome wide miRNA profiling in 96 liver samples (48 AA and 48 EA). The TaqMan OpenArray Human MicroRNA Panel contains 754 miRNAs with validated TaqMan miRNA assays derived from Sanger miRbase release v.14. Total RNA (10 ng) was reverse-transcribed using MegaplexTM pool RT-primers. Pre-amplification of cDNA with MegaplexTM PreAmp primers was performed according to the manufacturer’s protocol and recommendations (Thermo Fisher Scientific, United States). The pre-amplified products were diluted and mixed with TaqMan OpenArray Real-time PCR master mix and added to 384-well OpenArray Sample Loading plates. TaqMan OpenArray Human MicroRNA panels were automatically loaded by the AccuFill System and then placed in QuantStudio TM 12K Flex Real-Time PCR system for cycling.
Open array data quality control and statistical analysis
Cycle threshold (Ct) values with amplification scores (Amp Scores) below 1.24 and Cq confidence (Cq Conf) < 0.8 were filtered out. miRNAs with missing Ct values or Ct > 35 in >50% samples and samples with low miRNA detection were excluded. Data was normalized using global normalization by subtracting the individual Ct from the mean Ct of all miRNAs (Luo, 2012). Differential expression of miRNAs between race and sex was detected using the 2–ΔΔCT method (Livak and Schmittgen, 2001) and differences that exceeded 2-fold (FC > 2 or FC < 0.5) and p value ≤ 0.05 were considered statistically significant. Linear regression was used to test the association between the miRNAs and CYP3A expression, adjusting for race, sex, and age. All statistical analyses were carried out in the statistical programming environment R, using SAS v9.4 (Cary, NC) and Statistical Package for the Social Sciences SPSS v.26 (IBM Corp., United States).
miRNA targeted TaqMan assays for validation and data analysis
The TaqMan Advanced miRNA assays (Thermo Fisher Scientific, United States) were used to determine the expression of the two most significant miRNAs (miR-107 and miR-1260) in 231 liver samples for validation. Total RNA (10 ng) was reverse transcribed and quantitative real-time PCR was performed according to the manufacturer’s protocol and recommendations (Thermo Fisher Scientific, United States). From the discovery data set, we found that the combination of miR-132 and miR-484 was the most appropriate internal control for data normalization, as they had the most stable value of 0.003 according to Normfinder software (Andersen et al., 2004). Thus, the expression levels of miR-132 and mir-484 were also measured for normalization purposes. A multiple linear regression model was used to test the association between miRNAs and the expression of the CYP3As. We used forward and backward stepwise regression to select the best set of predictors in the multiple linear regression models with a cut-off p-value ≤ 0.05. Covariates considered for inclusion were: sex, race, age, and genotypes (CYP3A4*22 for CYP3A4 and CYP3A5*3 for CYP3A5, respectively), and nine TFs known to regulate CYP3A gene expression (Collins and Wang, 2021a). Modulated modularity clustering (MMC) (Stone and Ayroles, 2009) was used to explore the inter-miRNA relationships. MMC was designed to detect latent structure of the variance-covariance matrix using weighted graphs. The method searches for optimal community structure and detects the magnitude of pairwise relationships. The optimal number of clusters and the optimal cluster size were selected by using Spearman correlation.
Genotyping assays (rs62471956 and rs776746)
The single nucleotide polymorphism rs776746 (CYP3A5*3) and rs62471956, which is in complete linkage disequilibrium (LD) with rs35599367 (CYP3A4*22) (Collins and Wang, 2020), were genotyped using the OpenArray genotyping platform according to the manufacturer’s protocol (Life Technology, CA, United States).
Quantification of CYP3A4 protein in liver samples
CYP3A4 protein levels in 154 liver samples were measured using the capillary western blotting Jess system as described (Collins and Wang, 2021b). Briefly, total tissue lysates were prepared (Collins and Wang, 2021b) and a mouse anti-CYP3A4 antibody (R&D MAB 9079, 1:20 dilution) and NIR-conjugated anti-mouse secondary antibody (Biotechne, San Jose, CA, United States, 1:20 dilution) were used for detection. The total protein loaded in each lane was measured using the total protein channel in the Jess system. CYP3A4 protein (pmol) per mg total protein was calculated from a standard curve generated from a purified GST-fusion CYP3A4 protein (FisherScientific) measured on the Jess system. After log10 transformation, CYP3A4 protein levels followed a normal distribution as previously described (Collins and Wang, 2021b).
Results
Genome-wide expression profile of miRNAs in human liver samples
A total of 252 miRNAs in 91 samples (AA = 46, EA = 45) passed quality control. Three internal controls (U6rRNA, RNU44, and RNU48) and two miRNAs predicted to be tRNA fragments (HSA-MIR-1274B and HSA-MIR-720) were excluded, resulting in 247 miRNAs being included for further analysis. The expression levels between the 247 miRNAs varied drastically, with median Ct values ranging from 14.3 to 30.1 (>32,000-fold difference) (Supplementary Figure S1). Supplementary Table S2 lists the top 20 highest and lowest expressed miRNAs. Of these, miR-122, miR-24, and miR-19 showed the highest expression levels in liver, consistent with previous results (Rieger et al., 2013).
There was also significant variation in the expression of the miRNAs between the 91 samples. The inter-quantile range (IQR) for the Ct values ranged from 0.28 to 14.24 (Supplementary Figure S2). Nearly 60% (n = 149) of the miRNAs showed an IQR < 2 (<4-fold difference), while 10% (n = 26) of the miRNAs showed an IQR of 5–14 (32–16,384-fold). The miRNAs miR-1260, miR-29c, and miR-215 had the largest IQR of 11–14 (2048–16,384-fold).
To explore the relationship between the expression of the 247 miRNAs we used MMC (Stone and Ayroles, 2009). MMC identified seven clusters (Supplementary Table S3) with the number of miRNA members in each cluster ranging from 11 to 146. Cluster 1 had the strongest correlation between miRNAs (n = 15, average correlation = 0.72), while clusters 2–6 showed moderate correlation (n = 86, average correlation = 0.49–0.58). The majority of the miRNAs (n = 146) were clustered in Cluster 7 and were only weakly correlated (average correlation = 0.17).
Expression of the miRNAs differed significantly between samples originating from either AA or EA donors. Thirteen miRNAs were lower in AA compared to EA (fold-change < 0.5, p < 0.05) and 79 were higher (fold change > 2, p < 0.05) (Figure 1A). Similarly, the expression of the miRNAs differed between males and females, with nine miRNAs being lower (fold-change < 0.5, p < 0.05) and five being higher (fold-change > 2, p < 0.05) in females (Figure 1B). In addition, two miRNAs were positively correlated with age, while eight were negatively correlated (Supplementary Table S4). Thus, age, sex, and race were included as covariates for further association analyses.
FIGURE 1. Volcano plots showing differential expression of miRNAs between race (A) and sex (B) in liver samples. The plots indicate expression of significantly higher (blue dots) and lower (red dots) miRNAs in (A) AA vs. EA or (B) female vs. male. Negative log10 p-values are plotted on the y-axis, and log2 normalized fold change in expression is on the x-axis. For a miRNA to be considered as both significantly and differentially expressed a 2-fold difference and p ≤ 0.05 was required.
The association between the expression of miRNAs and the CYP3As
In a previous study, we showed mRNA expression for 12 CYP genes in the same liver cohort (Collins and Wang, 2021a). Using that dataset, we tested for association between the expression of the 247 miRNAs and four CYP3A genes using linear regression. As shown in Figure 2A, 21 miRNAs were associated with the expression of CYP3A4 (p < 0.05), however, after Bonferroni correction, only miR-107 and miR-1260 remained significant (p < 0.0002). Similar results were found with CYP3A5 and CYP3A43 (Figures 2B,C), but the associations between CYP3A5/miR-1260 and CYP3A43/miR-107 became insignificant after Bonferroni correction. In contrast, none of the miRNAs showed significant association with CYP3A7 after correction (Figure 2D).
FIGURE 2. Association between miRNAs and mRNA expression of four CYP3A genes in human liver samples. (A) CYP3A4; (B) CYP3A5; (C) CYP3A43 and (D) CYP3A7. miRNAs were measured using The TaqMan OpenArray Human MicroRNA panels. Linear regression was used to test the association between the 247 measured miRNAs and mRNA expression of the four CYP3As. Only miRNAs with an association p ≤ 0.05 are shown.
Since miRNAs can regulate gene expression via transcriptional or post-transcriptional pathways we also tested for association between the miRNAs and CYP3A4 protein. Previously, we reported that in the 179 human liver samples the average protein concentration is 101 pmol/mg, with a range of 6.9–257 pmol/mg, and that CYP3A4 mRNA and protein levels are moderately correlated (correlation coefficient = 0.778) (Collins and Wang, 2021b). We associated the expression of the miRNAs to the previous CYP3A4 protein dataset and found 18 miRNAs that were significantly associated (p < 0.05) with CYP3A4 protein levels (Figure 3), nine of which (miR-107, miR-1260, miR-766, miR-622, miR-30d, miR-92a, miR-193b, miR-1180 and miR-497) were also associated with CYP3A4 mRNA as shown in Figure 2. Notably, miR-107 was also found to be the most significant miRNA associated with CYP3A4 protein levels. Conversely, miR-1260 was not found to be significantly associated with CYP3A4 protein expression after Bonferroni correction.
FIGURE 3. Association between miRNAs and CYP3A4 protein expression in human liver samples. miRNAs were measured using The TaqMan OpenArray Human MicroRNA panels. Linear regression was used to test the association between the 247 measured miRNAs and CYP3A4 protein expression. Only miRNAs with an association p ≤ 0.05 are shown.
Validation of the association between miR-107/miR-1260 and the CYP3As
The genome-wide approach identified miR-107 and miR-1260 as the top two miRNAs associated with the expression of CYP3A4, CYP3A5, and CYP3A43 in human liver samples. To validate these results, we measured the expression of both miR-107 and miR-1260 using targeted TaqMan assays in our entire liver sample cohort (n = 231). After normalization by an internal control (combination of miR-132 and miR-484, see Methods), the expression levels of miR-107 and miR-1260 followed a normal distribution (Supplementary Figure S3). We then tested the association between these miRNAs and CYP3A expression. As expression of several miRNAs was affected by race, sex, and age, we included these covariates in our association models. Also, two known variants that reduce expression of the CYP3As, CYP3A4*22 (Wang et al., 2011) and CYP3A5*3 (Kuehl et al., 2001) were included to account for their impact on CYP3A expression. In agreement with our results from the discovery phase, miR-107 and miR-1260 were negatively associated with the expression levels of CYP3A4, CYP3A5, and CYP3A43, but not CYP3A7 (Table 1).
TABLE 1. Association between the expression levels of miR-107/miR-1260 and the CYP3A genes in human liver samples (n = 231) using linear regression analysis.
As expression of the CYP3As is known to be controlled by multiple TFs, we also tested for association between miR-107 and miR-1260 and the expression of nine key CYP3A TFs (ESR1, NR1I3, NR1I2, RXR, HNF4A, PPARA, FOXA2, AHR, and ARNT) that correlate with expression of the CYP3As (Collins and Wang, 2021a). Forward and backward stepwise regression (cut-off p ≤ 0.05) was used to select which TFs should be included in the multiple linear regression models, which varied for the different CYP3A members (Table 1). Interestingly, when the selected TFs were included in the models, the association between either miR-107 or miR-1260 with the CYP3As either became insignificant (CYP3A4 and CYP3A5) or much less significant (CYP3A43) (Table 1).
Consistent with our results showing miR-107 association with CYP3A4 mRNA, miR-107 is associated with CYP3A4 protein levels after adjusting for age, sex, race, and CYP3A4*22 genotype, and the strength of association is reduced after further adjusting for TFs (Table 2). Also, there was no association between expression levels of miR-1260 and CYP3A4 protein, in agreement with our genome-wide association results.
TABLE 2. Association between the expression levels of miR-107/miR-1260 and CYP3A4 protein levels in human liver samples (n = 154) using linear regression analysis.
These results indicate that the association between miR-107/miR-1260 and the CYP3A genes may be indirectly mediated through their action on TFs. Indeed, miR-107 is negatively associated with all nine TFs tested in this study, and miR-1260 is also negatively associated with the expression of NR1I3, NR1I2, FOXA2, and HNF4A (Table 3). As many of these TFs also regulate the other CYP enzymes (Wang et al., 2019; Collins and Wang, 2021a), we also tested for association between miR-107/miR-1260 and the other nine drug-metabolizing CYPs. As expected, miR-107 is negatively associated with all nine CYP enzymes in the CYP1 and CYP2 families (Table 4) after adjusting by age, sex, and race. Similar to what we observed for the CYP3As, when TFs were included in the regression model, the association between miR-107 and the CYPs became either less significant or insignificant. miR-1260 was only associated with the expression of the CYP3As, CYP1A1, and CYP1A2, but not the other seven CYP enzymes tested (Supplementary Table S5). After including TFs in the model, the associations between miR-1260 and two CYP1A enzymes also became less significant.
TABLE 3. Association between the expression levels of miR-107/miR-1260 and nine transcription factors tested in this study in human liver samples (n = 231) using linear regression analysis adjusting for covariates age, race, and sex.
TABLE 4. Association between expression levels of miR-107 and nine CYP genes in human liver samples (n = 231) using linear regression analysis with or without adjusting for transcription factors as indicated. Age, sex, and race were included as covariates in all analysis.
Discussion
Using an untargeted genome-wide miRNA expression profiling approach we identified miR-107 and miR-1260 as two miRNAs associating with the expression of CYP3A4, CYP3A5, and CYP3A43 in 91 liver samples (Figure 2). The results were then validated in a larger liver cohort of 231 samples (Tables 1, 2). Furthermore, we found an association between miR-107 and nine TFs known to regulate the expression of the CYP genes (Tables 3, 4). These results indicate a potential broad regulatory role of miR-107 on the expression of CYP enzymes via regulating the expression of key CYP controlling TFs. Moreover, we found profound differences in liver miRNA expression between race and sex. To our knowledge, this is the first study to explore genome-wide liver miRNA expression in a large number of samples originating from both European and African Americans and to also compare differences between the two racial backgrounds, gender, and to analyze associations between miRNAs and the CYP3A genes.
Potential regulatory roles of miR-107 and miR-1260 on the expression of CYP3A4 and the other drug metabolizing cytochrome P450 enzymes
Previous miRNA binding site prediction analysis followed by functional studies in cell models has linked 21 miRNAs to the expression of CYP3A4 via direct binding to its 3′UTR (Pan et al., 2009; Wei et al., 2014; Shi et al., 2015; Liu et al., 2016; Gill et al., 2017; Yan et al., 2017; Yu D et al, 2018) (Supplementary Table S6). However, it is unknown whether these miRNAs are expressed in the liver and whether their expression levels are correlated with CYP3A4. Probes to detect all 21 previously identified miRNAs were present in the TaqMan OpenArray Human MicroRNA panels used here. However, our study found that 13/21 (62%) were not readily detectable in liver samples and were therefore excluded from our further analyses. Among the eight miRNAs also measured in our samples, only one (miR-27b) was significantly associated with the expression of CYP3A4 mRNA (p = 0.035), while two others were significantly associated with CYP3A4 protein levels (miR-548-3p, p = 0.0086 and miR-629, p = 0.018) (Supplementary Table S6). These results indicate that most of the previously identified CYP3A4-relevant miRNAs play an uncertain role in the liver or that our approach did not readily detect their effects. Based on our results, miR-107 is the most significant miRNA associated with expression of CYP3A4, as well as CYP3A5, CYP3A43, and all nine CYP genes tested. These findings suggest a potential broad regulatory role of miR-107 on the expression of the CYP enzymes, and therefore, drug metabolism.
We also identified and validated miR-1260 as associated with the expression of CYP3A4, CYP3A5, CYP3A43, CYP1A1, and CYP1A2, but not the other seven CYP genes tested. Therefore, compared to miR-107, which appears to play a broad role in CYP regulation, the role of miR-1260 appears to be limited to the CYP3A and CYP1A families. Also, compared to miR-107, the expression level of miR-1260 in the liver is 97-fold lower, the strength of its association with the CYPs is weaker, and there is no association between miR-1260 and CYP3A4 protein levels. Thus, the in vivo impact of miR-1260 on CYP expression and drug metabolism may be less significant than miR-107.
The potential role of transcription factors and estrogen receptor alpha in mediating the association between miR-107 and the cytochrome P450s
Previously published computational analyses using multiple miRNA target prediction algorithms did not identify miR-107 binding sites in the 3′UTR of CYP3A4 (Ramamoorthy and Skaar, 2011; Swart and Dandara, 2014). Therefore, miR-107 regulation of CYP3A4 is unlikely to be mediated by direct binding to its 3′UTR. Instead, based on our results, it appears that miR-107 may regulate CYP3A4 expression indirectly through modulating the expression of key CYP-controlling TFs, as it is negatively associated with all nine TFs tested. Furthermore, inclusion of these TFs in our regression models caused the association between miR-107 and both CYP3A4 mRNA and protein to become insignificant (Tables 1, 2).
Of the key CYP-controlling TFs, the negative effect of miR-107 appears to be strongest for ESR1 (Table 3). Previously, we identified ligand-free ESR1 (i.e., ESR1 function in the absence of estrogen) as a master regulator mediating the expression of CYP3A4 and many other CYPs in hepatocytes and liver tissue (Wang et al., 2019; Collins and Wang, 2021a; Collins et al., 2021). Interestingly, miR-107 is predicted to directly bind to the 3′UTR of ESR1 and caused reduced ESR1 mRNA and protein expression in a cell model (Bao et al., 2019). It seems possible then, that the widespread negative association of miR-107 with the CYPs may be mediated by its direct, negative effect on ESR1. In support of this model, miR-107 is not associated with CYP3A7, which is a fetal enzyme that is not well associated with ESR1 or the other TFs (Collins and Wang, 2021a). Therefore, the effects of miR-107 appear to be limited to CYP enzymes regulated by ESR1 and the other TFs tested here. It is worth mentioning that ESR1 also regulates several of the other CYP-controlling TFs (e.g., NR1I2, FOXA2, HNF4A, and PPARA) (Wang et al., 2019) and may be the underlying reason why these TFs are also associated with miR-107 levels. However, there are likely non-ESR1 mediated miR-107 effects as well. For instance, NR1I3 was strongly associated with miR-107 expression, while ESR1 does not appear to regulate its expression (Wang et al., 2019).
In contrast to its apparent indirect effect on CYP3A4, miR-107 is reported to directly regulate CYP2C8, and potentially all three CYP2Cs, as miR-107 binding motifs are also in CYP2C9 and CYP2C19 (Zhang et al., 2012). Our results showed that all three CYP2C transcripts are strongly and negatively associated with miR-107 and that these associations disappear when including the TFs in the regression models (Table 4). These results indicate that miR-107 may regulate the expression of the CYP2C enzymes through both direct and indirect mechanisms and at both the transcriptional and post-transcriptional levels. As miRNAs have many potential targets, and regulation of the CYPs involves numerous regulatory proteins, the role of miR-107 in CYP regulation warrants further investigation.
Like miR-107, the regulatory role of miR-1260 also appears to be indirectly mediated by TFs, since including TFs in the regression models also abolished or reduced the association between miR-1260 and the CYPs. miR-1260 is associated with the expression of NR1I2, NR1I3, FOXA2, and HNF4A, but not ESR1. Consistently, miR-1260 is most strongly associated with the expression of CYP1A1 and CYP1A2, two CYPs that are only minimally associated with expression of ESR1 (Collins and Wang, 2021a). Therefore, it appears that miR-1260 may have a role in a different CYP regulatory network than miR-107.
Clinical implications and differences in miRNA profiles between race, sex, and age
miRNA expression is subject to change under many physiological conditions, disease states, and environmental effects (Gulyaeva and Kushlinskiy, 2016), potentially altering downstream pathways in response. miR-107 belongs to a highly conserved miR-15/miR-107 gene family that shares a common “AGCAGC” sequence near the end of the mature miRNAs’ 5′UTR (Finnerty et al., 2010). miR-107 is transcribed from the host gene PANK1, which is broadly expressed in many tissues and has been implicated in cancers (Turco et al., 2020), chemosensitivity (Chen et al., 2021), ischemic stroke (Yang et al., 2014), insulin resistance (Yang et al., 2020), and neuronal differentiation (Croizier et al., 2018). Dysregulation of miR-107 is associated with diverse pathways, including: obesity and diabetes (Foley and O'Neill, 2012), Alzheimer’s disease (Chang et al., 2017), cancers (Turco et al., 2020; Chen et al., 2021), and drug treatment and environmental exposure (Chou et al., 2017; Chen et al., 2020). In parallel, CYP substrate drug clearance and expression of the CYP enzymes also change in many of these events, such as obesity (Brill et al., 2012), diabetes (Darakjian et al., 2021), and cancers (Fahy et al., 2012). Whether miR-107 is contributing to altered expression of the CYPs under these adverse conditions is an interesting future direction.
In comparison, relatively little has been reported concerning miR-1260. Recent studies have shown that miR-1260 is upregulated in numerous cancers, including hepatic carcinoma (Yan et al., 2013) and implicated in chemosensitivity (Zhao et al., 2018), and hypoxia-induced vascular smooth muscle cell proliferation (Seong and Kang, 2020). Interestingly, both the CYP3As (Yu et al, 2018b) and CYP1As (McKinnon et al., 1991) are down-regulated in hepatic cancer cells, possibly mediated by increased miR-1260 expression in these cells. Also, in HepaRG cells (hepatic cell line) treatment with the well-characterized CYP3A inducer rifampicin also caused significant reduction in miR-1260 expression (Swart and Dandara, 2019). Rifampicin-induced CYP3A expression is mediated by NR1I2 (PXR) (Chai et al., 2013), which we found to be significantly and negatively associated with miR-1260 expression (Table 4). Therefore, while miR-1260 appears to function in a non-ESR1 regulatory network, its perturbation in response to cancer progression or drug treatment may play a significant role in regulation of the CYPs.
Our results showed large differences in the expression levels of miRNAs between liver samples from AA and EA donors (Figure 1). A previous study also identified 33 miRNAs that are differentially expressed in lymphoblastoid cell lines (LCLs) between residents of Utah (CEU) and Nigerians (YRI), and these differentially expressed miRNAs were significantly and inversely associated with the expression of at least one mRNA, indicating racial differences in miRNA expression may link to phenotypic differences between populations (Huang et al., 2011). Interestingly, of the 33 miRNAs showing racial expression differences in LCLs, 17 are expressed in the liver, 13 of which (76%) also showed racial expression differences in our data, indicating there may be common mechanisms underlying racial differences in miRNA expression across different tissues. Our previous results show that many CYPs and CYP-regulating TFs are differentially expressed between AA and EA (Collins and Wang, 2021a). However, there was no differential expression of miR-107 (p = 0.652) or miR-1260 (p = 0.192) between AA and EA in this study, indicating that miR-107/miR-1260 are unlikely to contribute to differential expression of the CYPs and TFs between AA and EA.
Sex-biased gene expression is pronounced in the liver (Zhang et al., 2011) and miRNAs may be contributing to these biases. We also found several miRNAs that differ in expression between the sexes (n = 13), but miR-107 (p = 0.644) and miR-1260 (p = 0.323) were not among these. Therefore, it appears that these two miRNAs do not contribute to previously reported sex-biased expression of the CYPs (Wolbold et al., 2003; Lamba et al., 2010).
Differential expression of several miRNAs in the liver has been associated with age (Rieger et al., 2013), and our results identified eight miRNAs to be weakly associated with age (Supplementary Table S4). Compared with the previous study (Rieger et al., 2013), only one miRNA (miR-31) was also identified. The reason for this discrepancy is unknown and requires further investigation. However, some of liver samples in previous study were derived from donors with liver diseases (Rieger et al., 2013) and it is unclear whether disease state may confound miRNA-age correlation.
In conclusion, we have identified both miR-107 and miR-1260 as two top miRNAs associated with the expression of the CYP enzymes. The effect of miR-107 appears to be through alteration of key CYP-regulating TFs, possibly through an ESR1-directed regulatory network. Thus, if validated, expression of miR-107 may serve as an indicator for the expression and activity of many relevant drug-metabolizing enzymes in the liver. Therefore, any physiological, pathophysiological, or environmental conditions that alter miR-107 have the potential to broadly influence drug metabolism and treatment outcomes.
Data availability statement
The miRNA expression datasets presented in this study can be found in Supplementary Table S7.
Ethics statement
Deidentified human liver samples were obtained from The Cooperative Human Tissue Network (CHTN) with written informed consent. University of Florida Institutional review board reviewed and approved the use of human tissue in this study.
Author contributions
MT designed and preformed experiments, analyzed data, and wrote manuscript; JC performed experiments and wrote manuscript; DW designed experiments, analyzed data, and wrote manuscript.
Funding
This work was supported by National Institute of Health (R35 GM140845 to DW). The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2022.943538/full#supplementary-material
References
Andersen, C. L., Jensen, J. L., and Ørntoft, T. F. (2004). Normalization of real-time quantitative reverse transcription-PCR data: A model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res. 64 (15), 5245–5250. doi:10.1158/0008-5472.CAN-04-0496
Bao, W., Zhang, Y., Li, S., Fan, Q., Qiu, M., Wang, Y., et al. (2019). miR-107-5p promotes tumor proliferation and invasion by targeting estrogen receptor-α in endometrial carcinoma. Oncol. Rep. 41 (3), 1575–1585. doi:10.3892/or.2018.6936
Brill, M. J., Diepstraten, J., van Rongen, A., van Kralingen, S., van den Anker, J. N., Knibbe, C. A., et al. (2012). Impact of obesity on drug metabolism and elimination in adults and children. Clin. Pharmacokinet. 51 (5), 277–304. doi:10.2165/11599410-000000000-00000
Catalanotto, C., Cogoni, C., and Zardo, G. (2016). MicroRNA in control of gene expression: An overview of nuclear functions. Int. J. Mol. Sci. 17 (10), E1712. doi:10.3390/ijms17101712
Chai, X., Zeng, S., and Xie, W. (2013). Nuclear receptors PXR and CAR: Implications for drug metabolism regulation, pharmacogenomics and beyond. Expert Opin. Drug Metab. Toxicol. 9 (3), 253–266. doi:10.1517/17425255.2013.754010
Chang, W. S., Wang, Y. H., Zhu, X. T., and Wu, C. J. (2017). Genome-wide profiling of miRNA and mRNA expression in alzheimer's disease. Med. Sci. Monit. 23, 2721–2731. doi:10.12659/msm.905064
Chen, H. A., Li, C. C., Lin, Y. J., Wang, T. F., Chen, M. C., Su, Y. H., et al. (2021). Hsa-miR-107 regulates chemosensitivity and inhibits tumor growth in hepatocellular carcinoma cells. Aging (Albany NY) 13 (8), 12046–12057. doi:10.18632/aging.202908
Chen, Y., Gong, W., Zhou, Y., Fan, R., Wu, Y., Pei, W., et al. (2020). Metformin up-regulated miR-107 expression and enhanced the inhibitory effect of miR-107 on gastric cancer growth. Transl. Cancer Res. 9 (4), 2941–2950. doi:10.21037/tcr.2020.03.25
Chou, W. C., Lee, P. H., Tan, Y. Y., Lin, H. C., Yang, C. W., Chen, K. H., et al. (2017). An integrative transcriptomic analysis reveals bisphenol A exposure-induced dysregulation of microRNA expression in human endometrial cells. Toxicol. Vitro 41, 133–142. doi:10.1016/j.tiv.2017.02.012
Collins, J. M., Huo, Z., and Wang, D. (2021). ESR1 ChIP-seq identifies distinct ligand-free ESR1 genomic binding sites in human hepatocytes and liver tissue. Int. J. Mol. Sci. 22 (3), 1461. doi:10.3390/ijms22031461
Collins, J. M., and Wang, D. (2020). Cis-acting regulatory elements regulating CYP3A4 transcription in human liver. Pharmacogenet. Genomics 30, 107–116. doi:10.1097/FPC.0000000000000402
Collins, J. M., and Wang, D. (2021a). Co-Expression of drug metabolizing cytochrome P450 enzymes and estrogen receptor alpha (ESR1) in human liver: Racial differences and the regulatory role of ESR1. Drug Metab. Pers. Ther. 36, 205–214. doi:10.1515/dmpt-2020-0160
Collins, J. M., and Wang, D. (2021b). Cytochrome P450 3A4 (CYP3A4) protein quantification using capillary Western blot technology and total protein normalization. J. Pharmacol. Toxicol. Methods 112, 107117. doi:10.1016/j.vascn.2021.107117
Croizier, S., Park, S., Maillard, J., and Bouret, S. G. (2018). Central Dicer-miR-103/107 controls developmental switch of POMC progenitors into NPY neurons and impacts glucose homeostasis. Elife 7, e40429. doi:10.7554/eLife.40429
Danielson, P. B. (2002). The cytochrome P450 superfamily: Biochemistry, evolution and drug metabolism in humans. Curr. Drug Metab. 3 (6), 561–597. doi:10.2174/1389200023337054
Darakjian, L., Deodhar, M., Turgeon, J., and Michaud, V. (2021). Chronic inflammatory status observed in patients with type 2 diabetes induces modulation of cytochrome P450 expression and activity. Int. J. Mol. Sci. 22 (9), 4967. doi:10.3390/ijms22094967
Fahy, B. N., Guo, T., and Ghose, R. (2012). Impact of hepatic malignancy on CYP3A4 gene expression. J. Surg. Res. 178 (2), 768–772. doi:10.1016/j.jss.2012.06.008
Finnerty, J. R., Wang, W. X., Hébert, S. S., Wilfred, B. R., Mao, G., Nelson, P. T., et al. (2010). The miR-15/107 group of microRNA genes: Evolutionary biology, cellular functions, and roles in human diseases. J. Mol. Biol. 402 (3), 491–509. doi:10.1016/j.jmb.2010.07.051
Foley, N. H., and O'Neill, L. A. (2012). miR-107: a toll-like receptor-regulated miRNA dysregulated in obesity and type II diabetes. J. Leukoc. Biol. 92 (3), 521–527. doi:10.1189/jlb.0312160
Gill, P., Bhattacharyya, S., McCullough, S., Letzig, L., Mishra, P. J., Luo, C., et al. (2017). MicroRNA regulation of CYP 1A2, CYP3A4 and CYP2E1 expression in acetaminophen toxicity. Sci. Rep. 7 (1), 12331. doi:10.1038/s41598-017-11811-y
Gufford, B. T., Robarge, J. D., Eadon, M. T., Gao, H., Lin, H., Liu, Y., et al. (2018). Rifampin modulation of xeno- and endobiotic conjugating enzyme mRNA expression and associated microRNAs in human hepatocytes. Pharmacol. Res. Perspect. 6 (2), e00386. doi:10.1002/prp2.386
Gulyaeva, L. F., and Kushlinskiy, N. E. (2016). Regulatory mechanisms of microRNA expression. J. Transl. Med. 14 (1), 143. doi:10.1186/s12967-016-0893-x
Hou, L., Wang, D., and Baccarelli, A. (2011). Environmental chemicals and microRNAs. Mutat. Res. 714 (1-2), 105–112. doi:10.1016/j.mrfmmm.2011.05.004
Huang, R. S., Gamazon, E. R., Ziliak, D., Wen, Y., Im, H. K., Zhang, W., et al. (2011). Population differences in microRNA expression and biological implications. RNA Biol. 8 (4), 692–701. doi:10.4161/rna.8.4.16029
Kalscheuer, S., Zhang, X., Zeng, Y., and Upadhyaya, P. (2008). Differential expression of microRNAs in early-stage neoplastic transformation in the lungs of F344 rats chronically treated with the tobacco carcinogen 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone. Carcinogenesis 29 (12), 2394–2399. doi:10.1093/carcin/bgn209
Kuehl, P., Zhang, J., Lin, Y., Lamba, J., Assem, M., Schuetz, J., et al. (2001). Sequence diversity in CYP3A promoters and characterization of the genetic basis of polymorphic CYP3A5 expression. Nat. Genet. 27 (4), 383–391. doi:10.1038/86882
Lamba, V., Panetta, J. C., Strom, S., and Schuetz, E. G. (2010). Genetic predictors of interindividual variability in hepatic CYP3A4 expression. J. Pharmacol. Exp. Ther. 332 (3), 1088–1099. doi:10.1124/jpet.109.160804
Li, D., Tolleson, W. H., Yu, D., Chen, S., Guo, L., Xiao, W., et al. (2019). Regulation of cytochrome P450 expression by microRNAs and long noncoding RNAs: Epigenetic mechanisms in environmental toxicology and carcinogenesis. J. Environ. Sci. Health. C Environ. Carcinog. Ecotoxicol. Rev. 37 (3), 180–214. doi:10.1080/10590501.2019.1639481
Liu, J. E., Ren, B., Tang, L., Tang, Q. J., Liu, X. Y., Li, X., et al. (2016). The independent contribution of miRNAs to the missing heritability in CYP3A4/5 functionality and the metabolism of atorvastatin. Sci. Rep. 6, 26544. doi:10.1038/srep26544
Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 25 (4), 402–408. doi:10.1006/meth.2001.1262
Londin, E., Loher, P., Telonis, A. G., Quann, K., Clark, P., Jing, Y., et al. (2015). Analysis of 13 cell types reveals evidence for the expression of numerous novel primate- and tissue-specific microRNAs. Proc. Natl. Acad. Sci. U. S. A. 112 (10), E1106–E1115. doi:10.1073/pnas.1420955112
Luo, S. (2012). MicroRNA expression analysis using the Illumina microRNA-Seq Platform. Methods Mol. Biol. 822, 183–188. doi:10.1007/978-1-61779-427-8_12
McKinnon, R. A., Hall, P. D., Quattrochi, L. C., Tukey, R. H., and McManus, M. E. (1991). Localization of CYP1A1 and CYP1A2 messenger RNA in normal human liver and in hepatocellular carcinoma by in situ hybridization. Hepatology 14 (5), 848–856. doi:10.1002/hep.1840140517
Mishra, P. J., Humeniuk, R., Longo-Sorbello, G. S., Banerjee, D., and Bertino, J. R. (2007). A miR-24 microRNA binding-site polymorphism in dihydrofolate reductase gene leads to methotrexate resistance. Proc. Natl. Acad. Sci. U. S. A. 104 (33), 13513–13518. doi:10.1073/pnas.0706217104
Mulder, T. A. M., van Eerden, R. A. G., de With, M., Elens, L., Hesselink, D. A., Matic, M., et al. (2021). Genotyping in clinical practice: Ready for implementation? Front. Genet. 12, 711943. doi:10.3389/fgene.2021.711943
Pan, Y. Z., Gao, W., and Yu, A. M. (2009). MicroRNAs regulate CYP3A4 expression via direct and indirect targeting. Drug Metab. Dispos. 37 (10), 2112–2117. doi:10.1124/dmd.109.027680
Ramamoorthy, A., and Skaar, T. C. (2011). In silico identification of microRNAs predicted to regulate the drug metabolizing cytochrome P450 genes. Drug Metab. Lett. 5 (2), 126–131. doi:10.2174/187231211795305258
Rieger, J. K., Klein, K., Winter, S., and Zanger, U. M. (2013). Expression variability of absorption, distribution, metabolism, excretion-related microRNAs in human liver: Influence of nongenetic factors and association with gene expression. Drug Metab. Dispos. 41 (10), 1752–1762. doi:10.1124/dmd.113.052126
Seong, M., and Kang, H. (2020). Hypoxia-induced miR-1260b regulates vascular smooth muscle cell proliferation by targeting GDF11. BMB Rep. 53 (4), 206–211. doi:10.5483/bmbrep.2020.53.4.136
Sharma, A. R., Vohra, M., Shukla, V., Guddattu, V., Razak Uk, A., Shetty, R., et al. (2020). Coding SNPs in hsa-miR-1343-3p and hsa-miR-6783-3p target sites of CYP2C19 modulates clopidogrel response in individuals with cardiovascular diseases. Life Sci. 245, 117364. doi:10.1016/j.lfs.2020.117364
Shi, Y., Liu, Y., Wei, Z., Zhang, Y., Zhang, L., Jiang, S., et al. (2015). Hsa-miR-27a is involved in the regulation of CYP3A4 expression in human livers from Chinese Han population. Pharmacogenomics 16 (12), 1379–1386. doi:10.2217/pgs.15.82
Stone, E. A., and Ayroles, J. F. (2009). Modulated modularity clustering as an exploratory tool for functional genomic inference. PLoS Genet. 5 (5), e1000479. doi:10.1371/journal.pgen.1000479
Swart, M., and Dandara, C. (2014). Genetic variation in the 3'-UTR of CYP1A2, CYP2B6, CYP2D6, CYP3A4, NR1I2, and UGT2B7: Potential effects on regulation by microRNA and pharmacogenomics relevance. Front. Genet. 5, 167. doi:10.3389/fgene.2014.00167
Swart, M., and Dandara, C. (2019). MicroRNA mediated changes in drug metabolism and target gene expression by efavirenz and rifampicin in vitro: Clinical implications. OMICS 23 (10), 496–507. doi:10.1089/omi.2019.0122
Szkolnicka, D., Lucendo-Villarin, B., Moore, J. K., Simpson, K. J., Forbes, S. J., Hay, D. C., et al. (2016). Reducing hepatocyte injury and necrosis in response to paracetamol using noncoding RNAs. Stem Cells Transl. Med. 5 (6), 764–772. doi:10.5966/sctm.2015-0117
Turco, C., Donzelli, S., and Fontemaggi, G. (2020). miR-15/107 microRNA gene group: Characteristics and functional implications in cancer. Front. Cell. Dev. Biol. 8, 427. doi:10.3389/fcell.2020.00427
Wang, D., Guo, Y., Wrighton, S. A., Cooke, G. E., and Sadee, W. (2011). Intronic polymorphism in CYP3A4 affects hepatic expression and response to statin drugs. Pharmacogenomics J. 11 (4), 274–286. doi:10.1038/tpj.2010.28
Wang, D., Lu, R., Rempala, G., and Sadee, W. (2019). Ligand-Free Estrogen Receptor alpha (ESR1) as master regulator for the expression of CYP3A4 and other cytochrome P450 enzymes in the human liver. Mol. Pharmacol. 96 (4), 430–440. doi:10.1124/mol.119.116897
Wang, D., and Sadee, W. (2016). CYP3A4 intronic SNP rs35599367 (CYP3A4*22) alters RNA splicing. Pharmacogenet. Genomics 26 (1), 40–43. doi:10.1097/FPC.0000000000000183
Wang, K., Zhang, S., Marzolf, B., Troisch, P., Brightman, A., Hu, Z., et al. (2009). Circulating microRNAs, potential biomarkers for drug-induced liver injury. Proc. Natl. Acad. Sci. U. S. A. 106 (11), 4402–4407. doi:10.1073/pnas.0813371106
Wei, Z., Jiang, S., Zhang, Y., Wang, X., Peng, X., Meng, C., et al. (2014). The effect of microRNAs in the regulation of human CYP3A4: A systematic study using a mathematical model. Sci. Rep. 4, 4283. doi:10.1038/srep04283
Wienkers, L. C., and Heath, T. G. (2005). Predicting in vivo drug interactions from in vitro drug discovery data. Nat. Rev. Drug Discov. 4 (10), 825–833. doi:10.1038/nrd1851
Wolbold, R., Klein, K., Burk, O., Nüssler, A. K., Neuhaus, P., Eichelbaum, M., et al. (2003). Sex is a major determinant of CYP3A4 expression in human liver. Hepatology 38 (4), 978–988. doi:10.1053/jhep.2003.50393
Yan, H., Wang, S., Yu, H., Zhu, J., and Chen, C. (2013). Molecular pathways and functional analysis of miRNA expression associated with paclitaxel-induced apoptosis in hepatocellular carcinoma cells. Pharmacology 92 (3-4), 167–174. doi:10.1159/000354585
Yan, L., Liu, J., Zhao, Y., Nie, Y., Ma, X., Kan, Q., et al. (2017). Suppression of miR-628-3p and miR-641 is involved in rifampin-mediated CYP3A4 induction in HepaRG cells. Pharmacogenomics 18 (1), 57–64. doi:10.2217/pgs-2016-0088
Yang, L., Zhang, B., Wang, X., Liu, Z., Li, J., Zhang, S., et al. (2020). P53/PANK1/miR-107 signalling pathway spans the gap between metabolic reprogramming and insulin resistance induced by high-fat diet. J. Cell. Mol. Med. 24 (6), 3611–3624. doi:10.1111/jcmm.15053
Yang, S. J., Yang, S. Y., Wang, D. D., Chen, X., Shen, H. Y., Zhang, X. H., et al. (2017). The miR-30 family: Versatile players in breast cancer. Tumour Biol. 39 (3), 1010428317692204. doi:10.1177/1010428317692204
Yang, X., Zhang, B., Molony, C., Chudin, E., Hao, K., Zhu, J., et al. (2010). Systematic genetic and genomic analysis of cytochrome P450 enzyme activities in human liver. Genome Res. 20 (8), 1020–1036. doi:10.1101/gr.103341.109
Yang, Z. B., Zhang, Z., Li, T. B., Lou, Z., Li, S. Y., Yang, H., et al. (2014). Up-regulation of brain-enriched miR-107 promotes excitatory neurotoxicity through down-regulation of glutamate transporter-1 expression following ischaemic stroke. Clin. Sci. 127 (12), 679–689. doi:10.1042/CS20140084
Yu, D., Wu, L., Gill, P., Tolleson, W. H., Chen, S., Sun, J., et al. (2018a). Multiple microRNAs function as self-protective modules in acetaminophen-induced hepatotoxicity in humans. Arch. Toxicol. 92 (2), 845–858. doi:10.1007/s00204-017-2090-y
Yu, T., Wang, X., Zhu, G., Han, C., Su, H., Liao, X., et al. (2018b). The prognostic value of differentially expressed CYP3A subfamily members for hepatocellular carcinoma. Cancer Manag. Res. 10, 1713–1726. doi:10.2147/CMAR.S159425
Zanger, U. M., and Schwab, M. (2013). Cytochrome P450 enzymes in drug metabolism: Regulation of gene expression, enzyme activities, and impact of genetic variation. Pharmacol. Ther. 138 (1), 103–141. doi:10.1016/j.pharmthera.2012.12.007
Zhang, S. Y., Surapureddi, S., Coulter, S., Ferguson, S. S., and Goldstein, J. A. (2012). Human CYP2C8 is post-transcriptionally regulated by microRNAs 103 and 107 in human liver. Mol. Pharmacol. 82 (3), 529–540. doi:10.1124/mol.112.078386
Zhang, Y., Klein, K., Sugathan, A., Nassery, N., Dombkowski, A., Zanger, U. M., et al. (2011). Transcriptional profiling of human liver identifies sex-biased genes associated with polygenic dyslipidemia and coronary artery disease. PLoS One 6 (8), e23506. doi:10.1371/journal.pone.0023506
Keywords: cytochrome P450, drug metabolism, microRNA, estrogen receptor alpha (ERα), gene expression, CYP3A4, miR-107
Citation: Tantawy M, Collins JM and Wang D (2022) Genome-wide microRNA profiles identify miR-107 as a top miRNA associating with expression of the CYP3As and other drug metabolizing cytochrome P450 enzymes in the liver. Front. Pharmacol. 13:943538. doi: 10.3389/fphar.2022.943538
Received: 13 May 2022; Accepted: 06 July 2022;
Published: 17 August 2022.
Edited by:
Brahim Achour, University of Rhode Island, United StatesReviewed by:
Meijuan Tu, UC Davis Medical Center, United StatesXiaobo Zhong, University of Connecticut, United States
Copyright © 2022 Tantawy, Collins and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Danxin Wang, RGFueGluLndhbmdAY29wLnVmbC5lZHU=