Liang Li

Cytochrome P450 3A5 (CYP3A5) is a critical drug-metabolizing enzyme in human hepatocytes, playing a key role in the metabolism of various clinically relevant drugs, with considerable variability in gene expression across individuals. Long non-coding RNAs (lncRNAs) are likely to be important regulators within the CYP3A5 regulatory network in the liver. In this study, we employed a computational approach to estimate Sobol's sensitivity indices (SSI) under generalized linear models, applied to liver RNA expression microarray data (GTEx v8). The SSI-based analysis revealed that the long non-coding RNA LINC02499 exhibits the highest SSI value in relation to CYP3A5 expression in the liver. Furthermore, we conducted a comprehensive evaluation of LINC02499's biological characteristics and confirmed its regulatory role in CYP3A5 expression through real-time quantitative PCR and Western blotting experiments at the cellular level. We also demonstrated the direct interaction between miR-329-5p and both LINC02499 and CYP3A5 mRNA using dual-luciferase reporter assays and electrophoretic mobility shift assays (EMSA). In conclusion, our findings outline a regulatory network in which LINC02499 functions as a molecular sponge, mitigating the post-transcriptional inhibitory effect of miR-329-5p on CYP3A5 expression.

Purpose: The aim of this study was to analyze the effects of various genetic polymorphisms and clinical factors on tacrolimus (TAC) concentration in the convalescence period (CP) and stabilization ...

CYP3A4, CYP3A5, and CYP3A7, which are located in a multigene locus (CYP3A), play crucial roles in drug metabolism. To understand the highly variable hepatic expression of CYP3As, regulatory network analyses have focused on transcription factors (TFs). Since long non-coding RNAs (lncRNAs) likely contribute to such networks, we assessed the regulatory effects of both TFs and lncRNAs on CYP3A expression in the human liver and small intestine, main organs of CYP3A expression. Using weighted gene co-expression network analysis (WGCNA) of GTEx v8 RNA expression data and multiple stepwise regression analysis, we constructed TF-lncRNA-CYP3A co-expression networks. Multiple lncRNAs and TFs displayed robust associations with CYP3A expression that differed between liver and small intestines (LINC02499, HNF4A-AS1, AC027682.6, LOC102724153, and RP11-503C24.6), indicating that lncRNAs contribute to variance in CYP3A expression in both organs. Of these, HNF4A-AS1 had been experimentally demonstrated to affect CYP3A expression. Incorporating ncRNAs into CYP3A expression regulatory network revealed additional candidate TFs associated with CYP3A expression. These results serve as a guide for experimental studies on lncRNA-TF regulation of CYP3A expression in the liver and small intestines.

The cytochrome P450 3As (CYP3As) are abundantly expressed in the liver and metabolize many commonly prescribed medications. Their expression is highly variable between individuals with little known genetic cause. Despite extensive investigation, cis‐acting genetic elements that control the expression of the CYP3As remain uncharacterized. Using chromatin conformation capture (4C assays), we detected reciprocal interaction between a distal regulatory region (DRR) and the CYP3A4 promoter. The DRR colocalizes with a variety of enhancer marks and was found to promote transcription in reporter assays. CRISPR‐mediated deletion of the DRR decreased expression of CYP3A4, CYP3A5, and CYP3A7, supporting its role as a shared enhancer regulating the expression of three CYP3A genes. Using reporter gene assays, we identified two single‐nucleotide polymorphisms (rs115025140 and rs776744/rs776742) that increased DRR‐driven luciferase reporter expression. In a liver cohort (n = 246), rs115025140 was associated with increased expression of CYP3A4 mRNA (1.8‐fold) and protein (1.6‐fold) and rs776744/rs776742 was associated with 1.39‐fold increased expression of CYP3A5 mRNA. The rs115025140 is unique to the African population and in a clinical cohort of African Americans taking statins for lipid control rs115025140 carriers showed a trend toward reduced statin‐mediated lipid reduction. In addition, using a published cohort of Chinese patients who underwent renal transplantation taking tacrolimus, rs776744/rs776742 carriers were associated with reduced tacrolimus concentration after adjusting for CYP3A5*3. Our results elucidate a complex regulatory network controlling expression of three CYP3A genes and identify two novel regulatory variants with potential clinical relevance for predicting CYP3A4 and CYP3A5 expression.

Tacrolimus (TAC) is an immunosuppressant widely used in kidney transplantation. TAC displays considerable interindividual variability in pharmacokinetics (PKs). Genetic and clinical factors play important roles in TAC PKs. We enrolled a total of 251 Chinese renal transplant recipients and conducted a genomewide association study (GWAS), linkage disequilibrium (LD), and one-way analysis of variance (ANOVA) to find genetic variants affecting log-transformed TAC trough blood concentration/dose ratio (log[C-0/D]). In addition, we performed dual luciferase reporter gene assays and multivariate regression models to evaluate the effect of the genetic variants. The GWAS results showed that all 23 genomewide significant single-nucleotide polymorphisms (p < 5 x 10(-8)) were located on chromosome 7, including CYP3A5*3. LD, conditional association analysis, and one-way ANOVA showed that rs75125371 T > C independently influenced TAC log(C-0/D). Dual luciferase reporter gene assays indicated that rs75125371 minor allele (C) was significantly associated with increased normalized luciferase activity than the major allele (T) in the Huh7 cells (p = 1.2 x 10(-5)) and HepaRG cells (p = 0.0097). A model inclusive of age, sex, hematocrit, CYP3A5*3, and rs75125371 explained 37.34% variance in TAC C-0. These results suggest that rs75125371 T > C is a functional and population-specific variant affecting TAC C-0 in Chinese renal transplant recipients.

The aim of this study was to identify novel genetic variants affecting tacrolimus trough blood concentrations. We analyzed the association between 58 single nucleotide polymorphisms (SNPs) across the CYP3A gene cluster and the log-transformed tacrolimus concentration/dose ratio (log (C₀/D)) in 819 renal transplant recipients (Discovery cohort). Multivariate linear regression was used to test for associations between tacrolimus log (C₀/D) and clinical factors. Luciferase reporter gene assays were used to evaluate the functions of select SNPs. Associations of putative functional SNPs with log (C₀/D) were further tested in 631 renal transplant recipients (Replication cohort). Nine SNPs were significantly associated with tacrolimus log (C₀/D) after adjustment for CYP3A5*3 and clinical factors. Dual luciferase reporter assays indicated that the rs4646450 G allele and rs3823812 T allele were significantly associated with increased normalized luciferase activity ratios (p < 0.01). Moreover, CYP3A7*2 was associated with higher TAC log(C₀/D) in the group of CYP3A5 expressers. Age, serum creatinine and hematocrit were significantly associated with tacrolimus log (C₀/D). CYP3A7*2, rs4646450, and rs3823812 are proposed as functional SNPs affecting tacrolimus trough blood concentrations in Chinese renal transplant recipients. Clinical factors also significantly affect tacrolimus metabolism.

Aim: To evaluate genetic variants affecting mycophenolic acid (MPA) metabolism in Chinese renal transplant recipients. Methods: Total 11 SNPs of UGT1A9, UGT1A8, UGT2B7, ABCC2, ABCG2 and SLCO1B3 were genotyped in 408 Chinese renal transplant recipients. Associations between SNPs and MPA concentration/dose ratio (C0/D) were analyzed using different genetic models. Multivariate linear regression was used to analyze associations between log (C0/D) and clinical factors. Results: After adjustment by clinical factors, UGT2B7 rs7662029 was associated with log (C0/D) using a dominant (p = 0.041) and an additive (p = 0.038) model, ABCC2 rs717620 was associated with log (C0/D) using a recessive model (p = 0.019). Using additive model, SNP–SNP interactions were identified (p = 0.002) between ABCC2 rs717620 and UGT1A9 rs2741049, with interactions (p = 0.002) between ABCC2 rs717620 and UGT1A8 rs1042597. Age, albumin and serum creatinine were associated with log (C0/D). Conclusion: rs7662029 and rs717620 may affect MPA ...

Vascular endothelial cells play an important role in modulating anti-thrombus and maintaining the natural function of vascular by secreting many active substances. β-boswellic acid (β-BA) is an active triterpenoid compound from the extract of boswellia serrate. In this study, it is demonstrated that β-BA ameliorates plasma coagulation parameters, protects endothelium from blood stasis induced injury and prevents blood stasis induced impairment of endothelium-dependent vasodilatation. Moreover, it is found that β-BA significantly increases nitric oxide (NO) and cyclic guanosine 3’, 5’-monophosphate (cGMP) levels in carotid aortas of blood stasis rats. To stimulate blood stasis-like conditionsin vitro, human umbilical vein endothelial cells (HUVECs) were exposed to transient oxygen and glucose deprivation (OGD). Treatment of β-BA significantly increased intracellular NO level. Western blot and immunofluorescence as well as immunohistochemistry reveal that β-BA increases phosphorylation of enzyme nitric oxide synthase (eNOS) at Ser1177. In addition, β-BA mediated endothelium-dependent vasodilatation can be markedly blocked by eNOS inhibitor L-NAME in blood stasis rats. In OGD treated HUEVCs, the protective effect of β-BA is attenuated by knockdown of eNOS. In conclusion, the above findings provide convincing evidence for the protective effects of β-BA on blood stasis induced endothelial dysfunction by eNOS signaling pathway.