Cytochrome P450 Allele CYP3A7*1C Associates with Adverse Outcomes in Chronic Lymphocytic Leukemia, Breast, and Lung Cancer

The long-term effectiveness and mechanism of oncolytic virotherapy combined with anti-PD-L1 antibody in colorectal cancer patient

Colorectal cancer (CRC) is known to be resistant to immunotherapy. In our phase-I clinical trial, one patient achieved a 313-day prolonged response during the combined treatment of oncolytic virotherapy and immunotherapy. To gain a deeper understanding of the potential molecular mechanisms, we performed a comprehensive multi-omics analysis on this patient and three non-responders. Our investigation unveiled that, initially, the tumor microenvironment (TME) of this responder presented minimal infiltration of T cells and natural killer cells, along with a relatively higher presence of macrophages compared to non-responders. Remarkably, during treatment, there was a progressive increase in CD4+ T cells, CD8+ T cells, and B cells in the responder's tumor tissue. This was accompanied by a significant upregulation of transcription factors associated with T-cell activation and cytotoxicity, including GATA3, EOMES, and RUNX3. Furthermore, dynamic monitoring of peripheral blood samples from the responder revealed a rapid decrease in circulating tumor DNA (ctDNA), suggesting its potential as an early blood biomarker of treatment efficacy. Collectively, our findings demonstrate the effectiveness of combined oncolytic virotherapy and immunotherapy in certain CRC patients and provide molecular evidence that virotherapy can potentially transform a "cold" TME into a "hot" one, thereby improving sensitivity to immunotherapy.

Genetically predicted serum testosterone and risk of gynecological disorders: a Mendelian randomization study

Background

Testosterone plays a key role in women, but the associations of serum testosterone level with gynecological disorders risk are inconclusive in observational studies.

Methods

We leveraged public genome-wide association studies to analyze the effects of four testosterone related exposure factors on nine gynecological diseases. Causal estimates were calculated by inverse variance-weighted (IVW), MR-Egger and weighted median methods. The heterogeneity test was performed on the obtained data through Cochrane's Q value, and the horizontal pleiotropy test was performed on the data through MR-Egger intercept and MR-PRESSO methods. "mRnd" online analysis tool was used to evaluate the statistical power of MR estimates.

Results

The results showed that total testosterone and bioavailable testosterone were protective factors for ovarian cancer (odds ratio (OR) = 0.885, P = 0.012; OR = 0.871, P = 0.005) and endometriosis (OR = 0.805, P = 0.020; OR = 0.842, P = 0.028) but were risk factors for endometrial cancer (OR = 1.549, P < 0.001; OR = 1.499, P < 0.001) and polycystic ovary syndrome (PCOS) (OR = 1.606, P = 0.019; OR = 1.637, P = 0.017). dehydroepiandrosterone sulfate (DHEAS) is a protective factor against endometriosis (OR = 0.840, P = 0.016) and premature ovarian failure (POF) (OR = 0.461, P = 0.046) and a risk factor for endometrial cancer (OR= 1.788, P < 0.001) and PCOS (OR= 1.970, P = 0.014). sex hormone-binding globulin (SHBG) is a protective factor against endometrial cancer (OR = 0.823, P < 0.001) and PCOS (OR = 0.715, P = 0.031).

Conclusion

Our analysis suggested causal associations between serum testosterone level and ovarian cancer, endometrial cancer, endometriosis, PCOS, POF.

Identification of Selective CYP3A7 and CYP3A4 Substrates and Inhibitors Using a High-Throughput Screening Platform

Cytochrome P450 (CYP) 3A7 is one of the major xenobiotic metabolizing enzymes in human embryonic, fetal, and newborn liver. CYP3A7 expression has also been observed in a subset of the adult population, including pregnant women, as well as in various cancer patients. The characterization of CYP3A7 is not as extensive as other CYPs, and health authorities have yet to provide guidance towards DDI assessment. To identify potential CYP3A7-specific molecules, we used a P450-Glo CYP3A7 enzyme assay to screen a library of ∼5,000 compounds, including FDA-approved drugs and drug-like molecules, and compared these screening data with that from a P450-Glo CYP3A4 assay. Additionally, a subset of 1,000 randomly selected compounds were tested in a metabolic stability assay. By combining the data from the qHTS P450-Glo and metabolic stability assays, we identified several chemical features important for CYP3A7 selectivity. Halometasone was chosen for further evaluation as a potential CYP3A7-selective inhibitor using molecular docking. From the metabolic stability assay, we identified twenty-two CYP3A7-selective substrates over CYP3A4 in supersome setting. Our data shows that CYP3A7 has ligand promiscuity, much like CYP3A4. Furthermore, we have established a large, high-quality dataset that can be used in predictive modeling for future drug metabolism and interaction studies.

Cross-ancestry Genome-wide Association Studies of Sex Hormone Concentrations in Pre- and Postmenopausal Women

OBJECTIVE

Concentrations of circulating sex hormones have been associated with a variety of diseases in women and are strongly influenced by menopausal status. We investigated the genetic architectures of circulating concentrations of estradiol, testosterone, and SHBG by menopausal status in women of European and African ancestry.

METHODS

Using data on 229 966 women from the UK Biobank, we conducted genome-wide association studies (GWASs) of circulating concentrations of estradiol, testosterone, and SHBG in premenopausal and postmenopausal women. We tested for evidence of heterogeneity of genetic effects by menopausal status and genetic ancestry. We conducted gene-based enrichment analyses to identify tissues in which genes with GWAS-enriched signals were expressed.

RESULTS

We identified 4 loci (5q35.2, 12q14.3, 19q13.42, 20p12.3) that were associated with detectable concentrations of estradiol in both pre- and postmenopausal women of European ancestry. Heterogeneity analysis identified 1 locus for testosterone (7q22.1) in the CYP3A7 gene and 1 locus that was strongly associated with concentrations of SHBG in premenopausal women only (10q15.1) near the AKR1C4 gene. Gene-based analysis of testosterone revealed evidence of enrichment of GWAS signals in genes expressed in adipose tissue for postmenopausal women. We did not find any evidence of ancestry-specific genetic effects for concentrations of estradiol, testosterone, or SHBG.

CONCLUSIONS

We identified specific loci that showed genome-wide significant evidence of heterogeneity by menopausal status for testosterone and SHBG. We also observed support for a more prominent role of genetic variants located near genes expressed in adipose tissue in determining testosterone concentrations among postmenopausal women.

Allele and genotype frequencies of CYP3A4, CYP3A5, CYP3A7, and GSTP1 gene polymorphisms among mainland Tibetan, Mongolian, Uyghur, and Han Chinese populations

Over 50% prescribed drugs are metabolised by cytochrome P450 3A (CYP3A) and glutathione S-transferase pi (GSTP1) adds a glutathione to the oxidative products by CYP3A, which increases the hydrophilic property of metabolites and facilitates the excretion. Single nucleotide polymorphisms (SNPs) of CYP3A and GSTP1 show a diverse allele and genotype frequencies distribution among the world populations. The present study aimed to investigate the genotype and allele frequency distribution patterns of CYP3A4, CYP3A5, CYP3A7 and GSTP1 polymorphisms among healthy participants in mainland Tibetan, Mongolian, Uyghur, and Han Chinese populations. Blood samples were collected from 842 unrelated healthy subjects (323 Tibetan, 134 Mongolian, 162 Uyghur, and 223 Han) for genotyping analysis. Variant allele frequencies of CYP3A4 rs2242480, CYP3A5 rs776746, CYP3A7 rs2257401, and GSTP1 Ile105Val were observed in Han (0.253, 0.686, 0.312 and 0.188), Tibetan (0.186, 0.819, 0.192 and 0.173), Mongolian (0.198, 0.784, 0.228 and 0.235) and Uyghur (0.179, 0.858, 0.182 and 0.250) respectively. The allele frequency of CYP3A7*1C in Uyghur (0.019) was higher than that in Tibetan (0.002, p < 0.01). There was a strong linkage disequilibrium between CYP3A4 rs2242480, CYP3A5 rs776746, and CYP3A7 rs2257401 among the four ethnic groups. The results might be useful for the precise medication in the Chinese populations.

OpenCYP: An open source database exploring human variability in activities and frequencies of polymophisms for major cytochrome P-450 isoforms across world populations

The open source database "OpenCYP database" has been developed based on the results of extensive literature searches from the peer-reviewed literature. OpenCYP provides data on human variability on baseline of activities and polymophism frequencies for selected cytochrome P-450 isoforms (CYP1A2, CYP2A6, CYP2D6, CYP3A4/3A5 and CYP3A7) in healthy adult populations from world populations. CYP enzymatic activities were generally expressed as the metabolic ratio (MR) between an unchanged probe drug and its metabolite(s) in urine or plasma measured in healthy adults. Data on other age groups were very limited and fragmented, constituting an important data gap. Quantitative comparisons were often hampered by the different experimental conditions used. However, variability was quite limited for CYP1A2, using caffeine as a probe substrate, with a symmetrical distribution of metabolic activity values. For CYP3A4, human variability was dependent on the probe substrate itself with low variability when data considering the dextromethorphan/demethilathed metabolite MR were used and large variability when the urinary 6β-hydroxycortisol/cortisol ratio was used. The largest variability in CYP activity was shown for CYP2D6 activity, after oral dosing of dextromethorphan, for which genetic polymorphisms are well characterised and constitute a significant source of variability. It is foreseen that the OpenCYP database can contribute to promising tools to support the further development of QIVIVE and PBK models for human risk assessment of chemicals particularly when combined with information on isoform-specific content in cells using proteomic approaches.

Promotion of rs3746804 (p. L267P) polymorphism to intracellular SLC52A3a trafficking and riboflavin transportation in esophageal cancer cells

Riboflavin is an essential micronutrient for normal cellular growth and function. Lack of dietary riboflavin is associated with an increased risk for esophageal squamous cell carcinoma (ESCC). Previous studies have identified that the human riboflavin transporter SLC52A3a isoform (encoded by SLC52A3) plays a prominent role in esophageal cancer cell riboflavin transportation. Furthermore, SLC52A3 gene single nucleotide polymorphisms rs3746804 (T>C, L267P) and rs3746803 (C >T, T278M) are associated with ESCC risk. However, whether SLC52A3a (p.L267P) and (p.T278M) act in riboflavin transportation in esophageal cancer cell remains inconclusive. Here, we constructed the full-length SLC52A3a protein fused to green fluorescent protein (GFP-SLC52A3a-WT and mutants L267P, T278M, and L267P/T278M). It was confirmed by immunofluorescence-based confocal microscopy that SLC52A3a-WT, L267P, T278M, and L267P/T278M expressed in cell membrane, as well as in a variety of intracellular punctate structures. The live cell confocal imaging showed that SLC52A3a-L267P and L267P/T278M increased the intracellular trafficking of SLC52A3a in ESCC cells. Fluorescence recovery after photobleaching of GFP-tagged SLC52A3a meant that intracellular trafficking of SLC52A3a-L267P and L267P/T278M was rapid dynamics process, leading to its stronger ability to transport riboflavin. Taken together, the above results indicated that the rs3746804 (p.L267P) polymorphism promoted intracellular trafficking of SLC52A3a and riboflavin transportation in ESCC cells.

Evaluation of the Association of Polymorphisms With Palbociclib-Induced Neutropenia: Pharmacogenetic Analysis of PALOMA-2/-3

BACKGROUND

The most frequently reported treatment-related adverse event in clinical trials with the cyclin-dependent kinase 4/6 (CDK4/6) inhibitor palbociclib is neutropenia. Allelic variants in ABCB1 and ERCC1 might be associated with early occurrence (i.e., end of week 2 treatment) of grade 3/4 neutropenia. Pharmacogenetic analyses were performed to uncover associations between single nucleotide polymorphisms (SNPs) in these genes, patient baseline characteristics, and early occurrence of grade 3/4 neutropenia.

MATERIALS AND METHODS

ABCB1 (rs1045642, rs1128503) and ERCC1 (rs3212986, rs11615) were analyzed in germline DNA from palbociclib-treated patients from PALOMA-2 (n = 584) and PALOMA-3 (n = 442). SNP, race, and cycle 1 day 15 (C1D15) absolute neutrophil count (ANC) data were available for 652 patients. Univariate and multivariable analyses evaluated associations between SNPs, patient baseline characteristics, and early occurrence of grade 3/4 neutropenia. Analyses were stratified by Asian (n = 122) and non-Asian (n = 530) ethnicity. Median progression-free survival (mPFS) was estimated using the Kaplan-Meier method. The effect of genetic variants on palbociclib pharmacokinetics was analyzed.

RESULTS

ABCB1 and ERCC1_rs11615 SNP frequencies differed between Asian and non-Asian patients. Multivariable analysis showed that low baseline ANC was a strong independent risk factor for C1D15 grade 3/4 neutropenia regardless of race (Asians: odds ratio [OR], 6.033, 95% confidence interval [CI], 2.615-13.922, p < .0001; Non-Asians: OR, 6.884, 95% CI, 4.138-11.451, p < .0001). ABCB1_rs1128503 (C/C vs. T/T: OR, 0.57, 95% CI, 0.311-1.047, p = .070) and ERCC1_rs11615 (A/A vs. G/G: OR, 1.75, 95% CI, 0.901-3.397, p = .098) were potential independent risk factors for C1D15 grade 3/4 neutropenia in non-Asian patients. Palbociclib mPFS was consistent across genetic variants; exposure was not associated with ABCB1 genotype.

CONCLUSION

This is the first comprehensive assessment of pharmacogenetic data in relationship to exposure to a CDK4/6 inhibitor. Pharmacogenetic testing may inform about potentially increased likelihood of patients developing severe neutropenia (NCT01740427, NCT01942135).

IMPLICATIONS FOR PRACTICE

Palbociclib plus endocrine therapy improves hormone receptor-positive/human epidermal growth factor receptor 2-negative advanced breast cancer outcomes, but is commonly associated with neutropenia. Genetic variants in ABCB1 may influence palbociclib exposure, and in ERCC1 are associated with chemotherapy-induced severe neutropenia. Here, the associations of single nucleotide polymorphisms in these genes and baseline characteristics with neutropenia were assessed. Low baseline absolute neutrophil count was a strong risk factor (p < .0001) for grade 3/4 neutropenia. There was a trend indicating that ABCB1_rs1128503 and ERCC1_rs11615 were potential risk factors (p < .10) for grade 3/4 neutropenia in non-Asian patients. Pharmacogenetic testing could inform clinicians about the likelihood of severe neutropenia with palbociclib.

Structural Basis for the Diminished Ligand Binding and Catalytic Ability of Human Fetal-Specific CYP3A7

Cytochrome P450 3A7 (CYP3A7) is a fetal/neonatal liver enzyme that participates in estriol synthesis, clearance of all-trans retinoic acid, and xenobiotic metabolism. Compared to the closely related major drug-metabolizing enzyme in adult liver, CYP3A4, the ligand binding and catalytic capacity of CYP3A7 are substantially reduced. To better understand the structural basis for these functional differences, the 2.15 Å crystal structure of CYP3A7 has been solved. Comparative analysis of CYP3A enzymes shows that decreased structural plasticity rather than the active site microenvironment defines the ligand binding ability of CYP3A7. In particular, a rotameric switch in the gatekeeping amino acid F304 triggers local and long-range rearrangements that transmit to the F-G fragment and alter its interactions with the I-E-D-helical core, resulting in a more rigid structure. Elongation of the β₃-β₄ strands, H-bond linkage in the substrate channel, and steric constraints in the C-terminal loop further increase the active site rigidity and limit conformational ensemble. Collectively, these structural distinctions lower protein plasticity and change the heme environment, which, in turn, could impede the spin-state transition essential for optimal reactivity and oxidation of substrates.

Crohn's disease-associated ATG16L1 T300A genotype is associated with improved survival in gastric cancer

BACKGROUND

A non-synonymous single nucleotide polymorphism of the ATG16L1 gene, T300A, is a major Crohn's disease (CD) susceptibility allele, and is known to be associated with increased apoptosis induction in the small intestinal crypt base in CD subjects and mouse models. We hypothesized that ATG16L1 T300A genotype also correlates with increased tumor apoptosis and therefore could lead to superior clinical outcome in cancer subjects.

METHODS

T300A genotyping by Taqman assay was performed for gastric carcinoma subjects who underwent resection from two academic medical centers. Transcriptomic analysis was performed by RNA-seq on formalin-fixed paraffin-embedded cancerous tissue. Tumor apoptosis and autophagy were determined by cleaved caspase-3 and p62 immunohistochemistry, respectively. The subjects' genotypes were correlated with demographics, various histopathologic features, transcriptome, and clinical outcome.

FINDINGS

Of the 220 genotyped subjects, 163 (74%) subjects carried the T300A allele(s), including 55 (25%) homozygous and 108 (49%) heterozygous subjects. The T300A/T300A subjects had superior overall survival than the other groups. Their tumors were associated with increased CD-like lymphoid aggregates and increased tumor apoptosis without concurrent increase in tumor mitosis or defective autophagy. Transcriptomic analysis showed upregulation of WNT/β-catenin signaling and downregulation of PPAR, EGFR, and inflammatory chemokine pathways in tumors of T300A/T300A subjects.

INTERPRETATION

Gastric carcinoma of subjects with the T300A/T300A genotype is associated with repressed EGFR and PPAR pathways, increased tumor apoptosis, and improved overall survival. Genotyping gastric cancer subjects may provide additional insight for clinical stratification.

Associations Between Cytochrome P450 (CYP) Gene Single-Nucleotide Polymorphisms and Second-to-Fourth Digit Ratio in Chinese University Students

BACKGROUND Cytochrome P450 (CYP) genes are necessary for the production or metabolism of fetal sex hormones during pregnancy. The second-to-fourth digit ratio (2D: 4D) is formed in the early stage of human fetal development and considered an indicator reflecting prenatal sex steroids levels. We explored the association between 2D: 4D and single-nucleotide polymorphisms (SNPs) of CYP. MATERIAL AND METHODS Correlation analysis between 2D: 4D and 8 SNPs, rs2687133 (CPY3A7), rs7173655 (CYP11A1), rs1004467, rs17115149, and rs2486758 (CYP17A1), and rs4646, rs2255192, rs4275794 (CYP19A1), was performed using data from 426 female and 412 male Chinese university students. SNP genotyping was conducted using PCR. Digit lengths were photographed and measured by image processing software. RESULTS rs2486758 (CYP17A1) correlated with left hand 2D: 4D in men (P=0.026), and rs1004467 (CYP17A1) correlated with right hand 2D: 4D in men (P=0.008) and the whole population (P=0.032). In men, allele G rs1004467 decreased right hand 2D: 4D, while allele C of rs2486758 increased left hand 2D: 4D. In women, left hand 2D: 4D was higher in genotypes with allele A of SNP rs4646 (CYP19A1) under the dominant genetic model; female DR-L was higher in genotypes with allele T of rs17115149 (CYP11A1). SNPs rs2687133 (CYP3A7) and rs1004467 (CYP17A1) were significantly correlated with right hand 2D: 4D (P=0.0107). CONCLUSIONS SNPs rs1004467 and rs2486758 of CYP17A1 are significant in the relationship between 2D: 4D and CYP gene polymorphisms under different conditions. SNP interactions between CYP genes probably impact 2D: 4D. The correlation between 2D: 4D and some sex hormone-related diseases may be due to the effect of CYP variants on the 2 phenotypes.

CYP3A7*1C allele: linking premenopausal oestrone and progesterone levels with risk of hormone receptor-positive breast cancers

BACKGROUND

Epidemiological studies provide strong evidence for a role of endogenous sex hormones in the aetiology of breast cancer. The aim of this analysis was to identify genetic variants that are associated with urinary sex-hormone levels and breast cancer risk.

METHODS

We carried out a genome-wide association study of urinary oestrone-3-glucuronide and pregnanediol-3-glucuronide levels in 560 premenopausal women, with additional analysis of progesterone levels in 298 premenopausal women. To test for the association with breast cancer risk, we carried out follow-up genotyping in 90,916 cases and 89,893 controls from the Breast Cancer Association Consortium. All women were of European ancestry.

RESULTS

For pregnanediol-3-glucuronide, there were no genome-wide significant associations; for oestrone-3-glucuronide, we identified a single peak mapping to the CYP3A locus, annotated by rs45446698. The minor rs45446698-C allele was associated with lower oestrone-3-glucuronide (-49.2%, 95% CI -56.1% to -41.1%, P = 3.1 × 10-18); in follow-up analyses, rs45446698-C was also associated with lower progesterone (-26.7%, 95% CI -39.4% to -11.6%, P = 0.001) and reduced risk of oestrogen and progesterone receptor-positive breast cancer (OR = 0.86, 95% CI 0.82-0.91, P = 6.9 × 10-8).

CONCLUSIONS

The CYP3A7*1C allele is associated with reduced risk of hormone receptor-positive breast cancer possibly mediated via an effect on the metabolism of endogenous sex hormones in premenopausal women.

Pharmacogenomics of anticancer drugs: Personalising the choice and dose to manage drug response

The field of pharmacogenomics has made great strides in oncology over the last 20 years and indeed a significant number of pre-emptive genetic tests are now routinely undertaken prior to anticancer drug administration. Many of these gene-drug interactions are the fruits of candidate gene and genome-wide association studies, which have largely focused on common genetic variants (allele frequency>1%). Examples where there is clinical utility include genotyping or phenotyping for G6PD to prevent rasburicase-induced RBC haemolysis, and TPMT to prevent thiopurine-induced bone marrow suppression. Other associations such as CYP2D6 status in determining the efficacy of tamoxifen are more controversial because of contradictory evidence from different sources, which has led to variability in the implementation of testing. As genomic technology becomes ever cheaper and more accessible, we must look to the additional data our genome can provide to explain interindividual variability in anticancer drug response. Clearly genes do not act on their own and it is therefore important to investigate genetic factors in conjunction with clinical factors, interacting concomitant drug therapies and other factors such as the microbiome, which can all affect drug disposition. Taking account of all of these factors, in conjunction with the somatic genome, is more likely to provide better predictive accuracy in determining anticancer drug response, both efficacy and safety. This review summarises the existing knowledge related to the pharmacogenomics of anticancer drugs and discusses areas of opportunity for further advances in personalisation of therapy in order to improve both drug safety and efficacy.

A Genome-wide Association Study of Circulating Levels of Atorvastatin and Its Major Metabolites

Atorvastatin (ATV) is frequently prescribed and generally well  tolerated, but can lead to myotoxicity, especially at higher doses. A genome-wide association study of circulating levels of ATV, 2-hydroxy (2-OH) ATV, ATV lactone (ATV L), and 2-OH ATV L was performed in 590 patients who had been hospitalized with a non-ST elevation acute coronary syndrome 1 month earlier and were on high-dose ATV (80 mg or 40 mg daily). The UGT1A locus (lead single nucleotide polymorphism, rs887829) was strongly associated with both increased 2-OH ATV/ATV (P = 7.25 × 10^-16 ) and 2-OH ATV L/ATV L (P = 3.95 × 10^-15 ) metabolic ratios. Moreover, rs45446698, which tags CYP3A7*1C, was nominally associated with increased 2-OH ATV/ATV (P = 6.18 × 10^-7 ), and SLCO1B1 rs4149056 with increased ATV (P = 2.21 × 10^-6 ) and 2-OH ATV (P = 1.09 × 10^-6 ) levels. In a subset of these patients whose levels of ATV and metabolites had also been measured at 12 months after hospitalization (n = 149), all of these associations remained, except for 2-OH ATV and rs4149056 (P = 0.057). Clinically, rs4149056 was associated with increased muscular symptoms (odds ratio (OR) 3.97; 95% confidence interval (CI) 1.29-12.27; P = 0.016) and ATV intolerance (OR 1.55; 95% CI 1.09-2.19; P = 0.014) in patients (n = 870) primarily discharged on high-dose ATV. In summary, both novel and recognized genetic associations have been identified with circulating levels of ATV and its major metabolites. Further study is warranted to determine the clinical utility of genotyping rs4149056 in patients on high-dose ATV.

Quality Control Measures and Validation in Gene Association Studies: Lessons for Acute Illness

Acute illness is a complex constellation of responses involving dysregulated inflammatory and immune responses, which are ultimately associated with multiple organ dysfunction. Gene association studies have associated single-nucleotide polymorphisms (SNPs) with clinical and pharmacological outcomes in a variety of disease states, including acute illness. With approximately 4 to 5 million SNPs in the human genome and recent studies suggesting that a large portion of SNP studies are not reproducible, we suggest that the ultimate clinical utility of SNPs in acute illness depends on validation and quality control measures. To investigate this issue, in December 2018 and January 2019 we searched the literature for peer-reviewed studies reporting data on associations between SNPs and clinical outcomes and between SNPs and pharmaceuticals (i.e., pharmacogenomics) published between January 2011 to February 2019. We review key methodologies and results from a variety of clinical and pharmacological gene association studies, including trauma and sepsis studies, as illustrative examples on current SNP association studies. In this review article, we have found three key points which strengthen the potential accuracy of SNP association studies in acute illness and other diseases: providing evidence of following a protocol quality control method such as the one in Nature Protocols or the OncoArray QC Guidelines; enrolling enough patients to have large cohort groups; and validating the SNPs using an independent technique such as a second study using the same SNPs with new patient cohorts. Our survey suggests the need to standardize validation methods and SNP quality control measures in medicine in general, and specifically in the context of complex disease states such as acute illness.

Graphene Oxide Causes Disordered Zonation Due to Differential Intralobular Localization in the Liver

The liver is the primary organ to sequester nanodrugs, representing a substantial hurdle for drug delivery and raising toxicity concerns. However, the mechanistic details underlying the liver sequestration and effects on the liver are still elusive. The difficulty in studying the liver lies in its complexity, which is structured with stringently organized anatomical units called lobules. Graphene oxide (GO) has attracted attention for its applications in biomedicine, especially as a nanocarrier; however, its sequestration and effects in the liver, the major enrichment and metabolic organ, are less understood. Herein, we unveiled the differential distribution of GO in lobules in the liver, with a higher amount surrounding portal triad zones than the central vein zones. Strikingly, liver zonation patterns also changed, as reflected by changes in vital zonated genes involved in hepatocyte integrity and metabolism, leading to compromised hepatic functions. RNA-Seq and DNA methylation sequencing analyses unraveled that GO-induced changes in liver functional zonation could be ascribed to dysregulation of key signaling pathways governing liver zonation at not only mRNA transcriptions but also DNA methylation imprinting patterns, partially through TET-dependent signaling. Together, this study reveals the differential GO distribution pattern in liver lobules and pinpoints the genetic and epigenetic mechanisms in GO-induced liver zonation alterations.

Influence of Genetic Variants on Steady-State Etonogestrel Concentrations Among Contraceptive Implant Users

OBJECTIVE

To identify genetic variants that influence steady-state etonogestrel concentrations among contraceptive implant users.

METHODS

We enrolled healthy, reproductive-age women in our pharmacogenomic study using etonogestrel implants for 12-36 months without concomitant use of hepatic enzyme inducers or inhibitors. We collected participant characteristics, measured serum etonogestrel concentrations, and genotyped each participant for 120 single nucleotide variants in 14 genes encoding proteins involved in steroid hormone (ie, estrogens, progestins) metabolism, regulation, or function. We performed generalized linear modeling to identify genetic variants associated with steady-state etonogestrel concentrations.

RESULTS

We enrolled 350 women, who had a median serum etonogestrel concentration of 137.4 pg/mL (range 55.8-695.1). Our final generalized linear model contained three genetic variants associated with serum etonogestrel concentrations: NR1I2(PXR) rs2461817 (β=13.36, P=.005), PGR rs537681 (β=-29.77, P=.007), and CYP3A7*1C (β=-35.06, P=.025). Variant allele frequencies were 69.4%, 84.9%, and 5.1%, respectively. Our linear model also contained two nongenetic factors associated with etonogestrel concentrations: body mass index (BMI) (β=-3.08, P=7.0×10) and duration of implant use (β=-1.60, P=5.8×10); R for the model =0.17.

CONCLUSION

Only BMI and duration of implant use remained significantly associated with steady-state etonogestrel concentrations. Of the three novel genetic associations found, one variant associated with increased etonogestrel metabolism (CYP3A7*1C) causes adult expression of fetal CYP3A7 proteins and can consequently alter steroid hormone metabolism. Women with this variant may potentially have increased metabolism of all steroid hormones, as 27.8% (5/18) of CYP3A7*1C carriers had serum etonogestrel concentrations that fell below the threshold for consistent ovulatory suppression (less than 90 pg/mL). More pharmacogenomic investigations are needed to advance our understanding of how genetic variation can influence the effectiveness and safety of hormonal contraception, and lay the groundwork for personalized medicine approaches in women's health.

CLINICAL TRIAL REGISTRATION

ClinicalTrials.gov, NCT03092037.

Neonatal cytochrome P450 CYP3A7: A comprehensive review of its role in development, disease, and xenobiotic metabolism

The human cytochrome P450 CYP3A7, once thought to be an enzyme exclusive to fetal livers, has more recently been identified in neonates and developing infants as old as 24 months post-gestational age. CYP3A7 has been demonstrated to metabolize two endogenous compounds that are known to be important in the growth and development of the fetus and neonate, namely dehydroepiandrosterone sulfate (DHEA-S) and all-trans retinoic acid (atRA). In addition, it is also known to metabolize a variety of drugs and xenobiotics, albeit generally to a lesser extent relative to CYP3A4/5. CYP3A7 is an important component in the development and protection of the fetal liver and additionally plays a role in certain disease states, such as cancer and adrenal hyperplasia. Ultimately, a full understanding of the expression, regulation, and metabolic properties of CYP3A7 is needed to provide neonates with appropriate individualized pharmacotherapy. This article summarizes the current state of knowledge of CYP3A7, including its discovery, distribution, alleles, RNA splicing, expression and regulation, metabolic properties, substrates, and inhibitors.

Stereoselective Oxidation Kinetics of Deoxycholate in Recombinant and Microsomal CYP3A Enzymes: Deoxycholate 19-Hydroxylation Is an In Vitro Marker of CYP3A7 Activity

The primary bile acids (BAs) synthesized from cholesterol in the liver are converted to secondary BAs by gut microbiota. It was recently disclosed that the major secondary BA, deoxycholate (DCA) species, is stereoselectively oxidized to tertiary BAs exclusively by CYP3A enzymes. This work subsequently investigated the in vitro oxidation kinetics of DCA at C-1β, C-3β, C-4β, C-5β, C-6α, C-6β, and C-19 in recombinant CYP3A enzymes and naive enzymes in human liver microsomes (HLMs). The stereoselective oxidation of DCA fit well with Hill kinetics at 1-300 μM in both recombinant CYP3A enzymes and pooled HLMs. With no contributions or trace contributions from CYP3A5, CYP3A7 favors oxidation at C-19, C-4β, C-6α, C-3β, and C-1β, whereas CYP3A4 favors the oxidation at C-5β and C-6β compared with each other. Correlation between DCA oxidation and testosterone 6β-hydroxylation in 14 adult single-donor HLMs provided proof-of-concept evidence that DCA 19-hydroxylation is an in vitro marker reaction for CYP3A7 activity, whereas oxidation at other sites represents mixed indicators for CYP3A4 and CYP3A7 activities. Deactivation caused by DCA-induced cytochrome P450-cytochrome P420 conversion, as shown by the spectral titrations of isolated CYP3A proteins, was observed when DCA levels were near or higher than the critical micelle concentration (about 1500 μM). Unlike CYP3A4, CYP3A7 showed abnormally elevated activities at 500 and 750 μM, which might be associated with an altered affinity for DCA multimers. The disclosed kinetic and functional roles of CYP3A isoforms in disposing of the gut bacteria-derived DCA may help in understanding the structural and functional mechanisms of CYP3A.

Continuum of Host-Gut Microbial Co-metabolism: Host CYP3A4/3A7 are Responsible for Tertiary Oxidations of Deoxycholate Species

The gut microbiota modifies endogenous primary bile acids (BAs) to produce exogenous secondary BAs, which may be further metabolized by cytochrome P450 enzymes (P450s). Our primary aim was to examine how the host adapts to the stress of microbe-derived secondary BAs by P450-mediated oxidative modifications on the steroid nucleus. Five unconjugated tri-hydroxyl BAs that were structurally and/or biologically associated with deoxycholate (DCA) were determined in human biologic samples by liquid chromatography-tandem mass spectrometry in combination with enzyme-digestion techniques. They were identified as DCA-19-ol, DCA-6β-ol, DCA-5β-ol, DCA-6α-ol, DCA-1β-ol, and DCA-4β-ol based on matching in-laboratory synthesized standards. Metabolic inhibition assays in human liver microsomes and recombinant P450 assays revealed that CYP3A4 and CYP3A7 were responsible for the regioselective oxidations of both DCA and its conjugated forms, glycodeoxycholate (GDCA) and taurodeoxycholate (TDCA). The modification of secondary BAs to tertiary BAs defines a host liver (primary BAs)-gut microbiota (secondary BAs)-host liver (tertiary BAs) axis. The regioselective oxidations of DCA, GDCA, and TDCA by CYP3A4 and CYP3A7 may help eliminate host-toxic DCA species. The 19- and 4β-hydroxylation of DCA species demonstrated outstanding CYP3A7 selectivity and may be useful as indicators of CYP3A7 activity.

Analytical Validation of Variants to Aid in Genotype-Guided Therapy for Oncology

The Clinical Laboratory Improvement Amendments of 1988 require that pharmacogenetic genotyping methods need to be established according to technical standards and laboratory practice guidelines before testing can be offered to patients. Testing methods for variants in ABCB1, CBR3, COMT, CYP3A7, C8ORF34, FCGR2A, FCGR3A, HAS3, NT5C2, NUDT15, SBF2, SEMA3C, SLC16A5, SLC28A3, SOD2, TLR4, and TPMT were validated in a Clinical Laboratory Improvement Amendments-accredited laboratory. Because no known reference materials were available, existing DNA samples were used for the analytical validation studies. Pharmacogenetic testing methods developed here were shown to be accurate and 100% analytically sensitive and specific. Other Clinical Laboratory Improvement Amendments-accredited laboratories interested in offering pharmacogenetic testing for these genetic variants, related to genotype-guided therapy for oncology, could use these publicly available samples as reference materials when developing and validating new genetic tests or refining current assays.

Genome-wide association study of offspring birth weight in 86 577 women identifies five novel loci and highlights maternal genetic effects that are independent of fetal genetics

Genome-wide association studies of birth weight have focused on fetal genetics, whereas relatively little is known about the role of maternal genetic variation. We aimed to identify maternal genetic variants associated with birth weight that could highlight potentially relevant maternal determinants of fetal growth. We meta-analysed data on up to 8.7 million SNPs in up to 86 577 women of European descent from the Early Growth Genetics (EGG) Consortium and the UK Biobank. We used structural equation modelling (SEM) and analyses of mother–child pairs to quantify the separate maternal and fetal genetic effects. Maternal SNPs at 10 loci (MTNR1B, HMGA2, SH2B3, KCNAB1, L3MBTL3, GCK, EBF1, TCF7L2, ACTL9, CYP3A7) were associated with offspring birth weight at P < 5 × 10⁻⁸. In SEM analyses, at least 7 of the 10 associations were consistent with effects of the maternal genotype acting via the intrauterine environment, rather than via effects of shared alleles with the fetus. Variants, or correlated proxies, at many of the loci had been previously associated with adult traits, including fasting glucose (MTNR1B, GCK and TCF7L2) and sex hormone levels (CYP3A7), and one (EBF1) with gestational duration. The identified associations indicate that genetic effects on maternal glucose, cytochrome P450 activity and gestational duration, and potentially on maternal blood pressure and immune function, are relevant for fetal growth. Further characterization of these associations in mechanistic and causal analyses will enhance understanding of the potentially modifiable maternal determinants of fetal growth, with the goal of reducing the morbidity and mortality associated with low and high birth weights.

The prognostic value of differentially expressed CYP3A subfamily members for hepatocellular carcinoma

Objective

The activities of four cytochrome P3A (CYP3A) subfamily members (CYP3A4, CYP3A5, CYP3A7, and CYP3A43) are well documented in drug metabolism. However, the association between CYP3A subfamily members and hepatocellular carcinoma (HCC) remains unclear. This study investigated the prognostic value of CYP3A subfamily mRNA expression levels with HCC prognosis.

Materials and methods

Data from a total of 360 HCC patients were retrieved from The Cancer Genome Atlas database, and data from 231 HCC patients were retrieved from the Gene Expression Omnibus database. Kaplan–Meier analysis and Cox regression models were utilized to determine median survival, overall survival, and recurrence-free survival. Hazard ratios and 95% CI were calculated.

Results

Low expression of CYP3A4, CYP3A5, and CYP3A43 in the tumor tissue was associated with short median survival (crude p=0.004, 0.001, and 0.001; adjusted p=0.022, 0.005, and 0.013, respectively). Joint-effects combination analysis of CYP3A4, CYP3A5/CYP3A4, CYP3A43/CYP3A5, and CYP3A43 revealed that high expression groups of two genes (group C, group c, group 3) were associated with a reduced risk of death, as compared to low expression of two genes (group A, group a, group 1), and the adjusted p values were 0.001, 0.004, and 0.001, respectively. Joint-effects analysis of CYP3A4, CYP3A5, and CYP3A43 showed that groups III and IV had a reduced risk of death, as compared to group I (adjusted p=0.024 and 0.002, respectively).

Conclusion

CYP3A4, CYP3A5, and CYP3A43 mRNA expression levels are potential prognostic markers of HCC.

Machine Learning-Assisted Network Inference Approach to Identify a New Class of Genes that Coordinate the Functionality of Cancer Networks

Emerging evidence indicates the existence of a new class of cancer genes that act as "signal linkers" coordinating oncogenic signals between mutated and differentially expressed genes. While frequently mutated oncogenes and differentially expressed genes, which we term Class I cancer genes, are readily detected by most analytical tools, the new class of cancer-related genes, i.e., Class II, escape detection because they are neither mutated nor differentially expressed. Given this hypothesis, we developed a Machine Learning-Assisted Network Inference (MALANI) algorithm, which assesses all genes regardless of expression or mutational status in the context of cancer etiology. We used 8807 expression arrays, corresponding to 9 cancer types, to build more than 2 × 108 Support Vector Machine (SVM) models for reconstructing a cancer network. We found that ~3% of ~19,000 not differentially expressed genes are Class II cancer gene candidates. Some Class II genes that we found, such as SLC19A1 and ATAD3B, have been recently reported to associate with cancer outcomes. To our knowledge, this is the first study that utilizes both machine learning and network biology approaches to uncover Class II cancer genes in coordinating functionality in cancer networks and will illuminate our understanding of how genes are modulated in a tissue-specific network contribute to tumorigenesis and therapy development.

CYP3A7*1C allele is associated with reduced levels of 2-hydroxylation pathway oestrogen metabolites

Background:

Endogenous sex hormones are well-established risk factors for breast cancer; the contribution of specific oestrogen metabolites (EMs) and/or ratios of specific EMs is less clear. We have previously identified a CYP3A7*1C allele that is associated with lower urinary oestrone (E₁) levels in premenopausal women. The purpose of this analysis was to determine whether this allele was associated with specific pathway EMs.

Methods:

We measured successfully 12 EMs in mid-follicular phase urine samples from 30 CYP3A7*1C carriers and 30 non-carriers using HPLC-MS/MS.

Results:

In addition to having lower urinary E₁ levels, CYP3A7*1C carriers had significantly lower levels of four of the 2-hydroxylation pathway EMs that we measured (2-hydroxyestrone, P=1.1 × 10⁻¹²; 2-hydroxyestradiol, P=2.7 × 10⁻⁷; 2-methoxyestrone, P=1.9 × 10⁻¹²; and 2-methoxyestradiol, P=0.0009). By contrast, 16α-hydroxylation pathway EMs were slightly higher in carriers and significantly so for 17-epiestriol (P=0.002).

Conclusions:

The CYP3A7*1C allele is associated with a lower urinary E₁ levels, a more pronounced reduction in 2-hydroxylation pathway EMs and a lower ratio of 2-hydroxylation:16α-hydroxylation EMs in premenopausal women. To further characterise the association between parent oestrogens, EMs and subsequent risk of breast cancer, characterisation of additional genetic variants that influence oestrogen metabolism and large prospective studies of a broad spectrum of EMs will be required.