A correlative study of the genomic underpinning of virulence traits and drug tolerance of Candida auris.

Candida auris is an opportunistic fungal pathogen with high mortality rates which presents a clear threat to public health. The risk of C. auris infection is high because it can colonize the body, resist antifungal treatment, and evade the immune system. The genetic mechanisms for these traits are not well known. Identifying them could lead to new targets for new treatments. To this end, we present an analysis of the genetics and gene expression patterns of C. auris carbon metabolism, drug resistance, and macrophage interaction. We chose to study two C. auris isolates simultaneously, one drug sensitive (B11220 from Clade II) and one drug resistant (B11221 from Clade III). Comparing the genomes, we confirm the previously reported finding that B11220 was missing a 12.8 kb region on chromosome VI. This region contains a gene cluster encoding proteins related to alternative sugar utilization. We show that B11221, which has the gene cluster, readily assimilates and utilizes D-galactose and L-rhamnose as compared to B11220, which harbors the deletion. B11221 exhibits increased adherence and drug resistance compared to B11220 when grown in these sugars. Transcriptomic analysis of both isolates grown on glucose or galactose showed that the gene cluster was upregulated when grown on D-galactose. These findings reinforce growing evidence of a link between metabolism and drug tolerance. B11221 resists phagocytosis by macrophages and exhibits decreased ß-1,3-glucan exposure, a key determinant that allows Candida to evade the host immune system, as compared to B11220. In a transcriptomic analysis of both isolates co-cultured with macrophages, we find upregulation of genes associated with transport and transcription factors in B11221. Our studies show a positive correlation between membrane composition and immune evasion, alternate sugar utilization, and drug tolerance in C. auris.

Factors associated with medical consumable availability in level 1 facilities in Malawi: a secondary analysis of a facility census.

BACKGROUND: Medical consumable stock-outs negatively affect health outcomes not only by impeding or delaying the effective delivery of services but also by discouraging patients from seeking care. Consequently, supply chain strengthening is being adopted as a key component of national health strategies. However, evidence on the factors associated with increased consumable availability is limited. METHODS: In this study, we used the 2018-19 Harmonised Health Facility Assessment data from Malawi to identify the factors associated with the availability of consumables in level 1 facilities, ie, rural hospitals or health centres with a small number of beds and a sparsely equipped operating room for minor procedures. We estimate a multilevel logistic regression model with a binary outcome variable representing consumable availability (of 130 consumables across 940 facilities) and explanatory variables chosen based on current evidence. Further subgroup analyses are carried out to assess the presence of effect modification by level of care, facility ownership, and a categorisation of consumables by public health or disease programme, Malawi's Essential Medicine List classification, whether the consumable is a drug or not, and level of average national availability. FINDINGS: Our results suggest that the following characteristics had a positive association with consumable availability-level 1b facilities or community hospitals had 64% (odds ratio [OR] 1·64, 95% CI 1·37-1·97) higher odds of consumable availability than level 1a facilities or health centres, Christian Health Association of Malawi and private-for-profit ownership had 63% (1·63, 1·40-1·89) and 49% (1·49, 1·24-1·80) higher odds respectively than government-owned facilities, the availability of a computer had 46% (1·46, 1·32-1·62) higher odds than in its absence, pharmacists managing drug orders had 85% (1·85, 1·40-2·44) higher odds than a drug store clerk, proximity to the corresponding regional administrative office (facilities greater than 75 km away had 21% lower odds [0·79, 0·63-0·98] than facilities within 10 km of the district health office), and having three drug order fulfilments in the 3 months before the survey had 14% (1·14, 1·02-1·27) higher odds than one fulfilment in 3 months. Further, consumables categorised as vital in Malawi's Essential Medicine List performed considerably better with 235% (OR 3·35, 95% CI 1·60-7·05) higher odds than other essential or non-essential consumables and drugs performed worse with 79% (0·21, 0·08-0·51) lower odds than other medical consumables in terms of availability across facilities. INTERPRETATION: Our results provide evidence on the areas of intervention with potential to improve consumable availability. Further exploration of the health and resource consequences of the strategies discussed will be useful in guiding investments into supply chain strengthening. FUNDING: UK Research and Innovation as part of the Global Challenges Research Fund (Thanzi La Onse; reference MR/P028004/1), the Wellcome Trust (Thanzi La Mawa; reference 223120/Z/21/Z), the UK Medical Research Council, the UK Department for International Development, and the EU (reference MR/R015600/1).

Ensemble Detection of DNA Engineering Signatures.

Synthetic biology is creating genetically engineered organisms at an increasing rate for many potentially valuable applications, but this potential comes with the risk of misuse or accidental release. To begin to address this issue, we have developed a system called GUARDIAN that can automatically detect signatures of engineering in DNA sequencing data, and we have conducted a blinded test of this system using a curated Test and Evaluation (T&E) data set. GUARDIAN uses an ensemble approach based on the guiding principle that no single approach is likely to be able to detect engineering with perfect accuracy. Critically, ensembling enables GUARDIAN to detect sequence inserts in 13 target organisms with a high degree of specificity that requires no subject matter expert (SME) review.

Endothelial SMAD1/5 signaling couples angiogenesis to osteogenesis in juvenile bone.

Skeletal development depends on coordinated angiogenesis and osteogenesis. Bone morphogenetic proteins direct bone formation in part by activating SMAD1/5 signaling in osteoblasts. However, the role of SMAD1/5 in skeletal endothelium is unknown. Here, we found that endothelial cell-conditional SMAD1/5 depletion in juvenile mice caused metaphyseal and diaphyseal hypervascularity, resulting in altered trabecular and cortical bone formation. SMAD1/5 depletion induced excessive sprouting and disrupting the morphology of the metaphyseal vessels, with impaired anastomotic loop formation at the chondro-osseous junction. Endothelial SMAD1/5 depletion impaired growth plate resorption and, upon long-term depletion, abrogated osteoprogenitor recruitment to the primary spongiosa. Finally, in the diaphysis, endothelial SMAD1/5 activity was necessary to maintain the sinusoidal phenotype, with SMAD1/5 depletion inducing formation of large vascular loops and elevated vascular permeability. Together, endothelial SMAD1/5 activity sustains skeletal vascular morphogenesis and function and coordinates growth plate remodeling and osteoprogenitor recruitment dynamics in juvenile mouse bone.

The changes in health service utilisation in Malawi during the COVID-19 pandemic.

INTRODUCTION: The COVID-19 pandemic and the restriction policies implemented by the Government of Malawi may have disrupted routine health service utilisation. We aimed to find evidence for such disruptions and quantify any changes by service type and level of health care. METHODS: We extracted nationwide routine health service usage data for 2015-2021 from the electronic health information management systems in Malawi. Two datasets were prepared: unadjusted and adjusted; for the latter, unreported monthly data entries for a facility were filled in through systematic rules based on reported mean values of that facility or facility type and considering both reporting rates and comparability with published data. Using statistical descriptive methods, we first described the patterns of service utilisation in pre-pandemic years (2015-2019). We then tested for evidence of departures from this routine pattern, i.e., service volume delivered being below recent average by more than two standard deviations was viewed as a substantial reduction, and calculated the cumulative net differences of service volume during the pandemic period (2020-2021), in aggregate and within each specific facility. RESULTS: Evidence of disruptions were found: from April 2020 to December 2021, services delivered of several types were reduced across primary and secondary levels of care-including inpatient care (-20.03% less total interactions in that period compared to the recent average), immunisation (-17.61%), malnutrition treatment (-34.5%), accidents and emergency services (-16.03%), HIV (human immunodeficiency viruses) tests (-27.34%), antiretroviral therapy (ART) initiations for adults (-33.52%), and ART treatment for paediatrics (-41.32%). Reductions of service volume were greatest in the first wave of the pandemic during April-August 2020, and whereas some service types rebounded quickly (e.g., outpatient visits from -17.7% to +3.23%), many others persisted at lower level through 2021 (e.g., under-five malnutrition treatment from -15.24% to -42.23%). The total reduced service volume between April 2020 and December 2021 was 8 066 956 (-10.23%), equating to 444 units per 1000 persons. CONCLUSION: We have found substantial evidence for reductions in health service delivered in Malawi during the COVID-19 pandemic which may have potential health consequences, the effect of which should inform how decisions are taken in the future to maximise the resilience of healthcare system during similar events.

SVR Chemometrics to Quantify ß-Lactoglobulin and α-Lactalbumin in Milk Using MIR.

Protein content variation in milk can impact the quality and consistency of dairy products, necessitating access to in-line real time monitoring. Here, we present a chemometric approach for the qualitative and quantitative monitoring of ß-lactoglobulin and α-lactalbumin, using mid-infrared spectroscopy (MIR). In this study, we employed Hotelling T2 and Q-residual for outlier detection, automated preprocessing using nippy, conducted wavenumber selection with genetic algorithms, and evaluated four chemometric models, including partial least squares, support vector regression (SVR), ridge, and logistic regression to accurately predict the concentrations of ß-lactoglobulin and α-lactalbumin in milk. For the quantitative analysis of these two whey proteins, SVR performed the best to interpret protein concentration from 197 MIR spectra originating from 42 Cornell University samples of preserved pasteurized modified milk. The R2 values obtained for ß-lactoglobulin and α-lactalbumin using leave one out cross-validation (LOOCV) are 92.8% and 92.7%, respectively, which is the highest correlation reported to date. Our approach introduced a combination of preprocessing automation, genetic algorithm-based wavenumber selection, and used Optuna to optimize the framework for tuning hyperparameters of the chemometric models, resulting in the best chemometric analysis of MIR data to quantitate ß-lactoglobulin and α-lactalbumin to date.

YAP and TAZ couple osteoblast precursor mobilization to angiogenesis and mechanoregulation in murine bone development.

Fetal bone development occurs through the conversion of avascular cartilage to vascularized bone at the growth plate. This requires coordinated mobilization of osteoblast precursors with blood vessels. In adult bone, vessel-adjacent osteoblast precursors are maintained by mechanical stimuli; however, the mechanisms by which these cells mobilize and respond to mechanical cues during embryonic development are unknown. Here, we show that the mechanoresponsive transcriptional regulators Yes-associated protein (YAP) and transcriptional co-activator with PDZ-binding motif (TAZ) spatially couple osteoblast precursor mobilization to angiogenesis, regulate vascular morphogenesis to control cartilage remodeling, and mediate mechanoregulation of embryonic murine osteogenesis. Mechanistically, YAP and TAZ regulate a subset of osteoblast-lineage cells, identified by single-cell RNA sequencing as vessel-associated osteoblast precursors, which regulate transcriptional programs that direct blood vessel invasion through collagen-integrin interactions and Cxcl12. Functionally, in 3D human cell co-culture, CXCL12 treatment rescues angiogenesis impaired by stromal cell YAP/TAZ depletion. Together, these data establish functions of the vessel-associated osteoblast precursors in bone development.

Modeling Contraception and Pregnancy in Malawi: A Thanzi La Onse Mathematical Modeling Study.

Malawi has high unmet need for contraception with a costed national plan to increase contraception use. Estimating how such investments might impact future population size in Malawi can help policymakers understand effects and value of policies to increase contraception uptake. We developed a new model of contraception and pregnancy using individual-level data capturing complexities of contraception initiation, switching, discontinuation, and failure by contraception method, accounting for differences by individual characteristics. We modeled contraception scale-up via a population campaign to increase initiation of contraception (Pop) and a postpartum family planning intervention (PPFP). We calibrated the model without new interventions to the UN World Population Prospects 2019 medium variant projection of births for Malawi. Without interventions Malawi's population passes 60 million in 2084; with Pop and PPFP interventions. it peaks below 35 million by 2100. We compare contraception coverage and costs, by method, with and without interventions, from 2023 to 2050. We estimate investments in contraception scale-up correspond to only 0.9 percent of total health expenditure per capita though could result in dramatic reductions of current pressures of very rapid population growth on health services, schools, land, and society, helping Malawi achieve national and global health and development goals.

Quantitative analysis of the UDP-glucuronosyltransferase transcriptome in human tissues.

UDP-glucuronosyltransferases (UGTs) are phase II drug metabolizing enzymes that play important roles in the detoxification of endogenous and exogenous substrates. The 22 human UGTs belong to four families (UGT1, UGT2, UGT3, and UGT8) and differ in their expression, substrate specificity, UDP-sugar preference, and physiological functions. Differential expression/activity of the UGTs contributes to interperson variability in drug responses and toxicity, hormone homeostasis, and disease/cancer risks. However, in normal tissues, the tissue-specific expression profiles and transcriptional regulation of the UGTs are still not fully understood. In this study, we comprehensively analyzed the transcriptome of 22 UGTs in 54 human tissues/regions using RNAseq data from GTEx. We then validated the findings in the liver and small intestine samples using real-time PCR. Our results showed large interindividual variability across tissues in the expression of each UGT and the overall composition of UGT pools, consisting of different UGTs and their splice isoforms. Our results also revealed coexpression of the UGTs, Cytochrome P450s, and many transcription factors in the liver, suggesting potential coregulation or functional coordination. Our results provide the groundwork for future studies to detail further the regulation of the expression and activity of the UGTs.

Isolation and sequencing of the spoilage yeast Kregervanrija delftensis.

New sequencing workflows can enable the genomics of microbes isolated from craft beverages. Here, we use hybrid, short, and long-read sequencing to assemble the genome of a yeast isolated from cider, Kregervanrija delftensis NCC-J. The Kregervanrija genus has few genomes available; thus, this contributes to yeast genomics.

A Comprehensive Evaluation of the Effects of RNA-Editing Proteins ADAR and ADARB1 on the Expression of the Drug-Metabolizing Enzymes in HepaRG Cells.

Two RNA-editing proteins, the adenosine deaminase acting on RNA, ADAR, and ADARB1, broadly regulate gene expression in editing-dependent and editing-independent manners. Previous studies showed that the expression of the drug-metabolizing cytochrome P450s (P450s) and UDP glucuronosyltransferases (UGTs) changes upon knockdown (KD) of ADAR or ADARB1 in different hepatic cell lines. To systematically survey the effects of these two ADARs on the expression of P450s and UGTs, we used small interfering RNA in HepaRG cells and tested the association between the expression of the P450s and ADARs in a liver sample cohort (n = 246). KD of ADAR increased the expression of the CYP3As and CYP2C9 and reduced the expression of the others, whereas KD of ADARB1 reduced the expression of nearly all genes tested. ADAR KD also suppressed the induction of most P450s, whereas ADARB1 KD had mixed effects depending on the inducer/gene combination. P450 expression was positively associated with both ADARs in liver samples, consistent with the KD results. However, after adjusting for the expression of transcription factors (TFs) known to regulate P450 expression, the associations disappeared, indicating that the effects of ADAR or ADARB1 primarily occur through TFs. Moreover, we found that the expression of normally spliced CYP3A5 transcripts is increased in both KDs, indicating a direct effect of the ADARs on promoting the usage of the cryptic splice site generated by CYP3A5*3. Taken together, our results revealed the nonoverlapping regulatory effects of ADAR and ADARB1 and supported their broad roles in controlling the expression of drug-metabolizing enzymes in the liver. SIGNIFICANCE STATEMENT: Here, this study systematically surveyed the roles of ADAR and ADARB1 in both basal and induced expression of drug-metabolizing enzymes and assessed their coexpression in liver samples. This study's results support that ADAR and ADARB1 regulate the expression of the drug-metabolizing enzymes in the liver, suggesting that factors affecting ADAR expression also have the potential to impact drug metabolism.

Draft genome of the Ogataea polymorpha type strain CBS 4732.

The methylotrophic yeast Ogataea polymorpha is of significant biotechnological interest, particularly in high-temperature fermentations and for recombinant protein production. Here, we present a high-quality genome assembly for the O. polymorpha type strain CBS 4732.

A frequent CYP2D6 variant promotes skipping of exon 3 and reduces CYP2D6 protein expression in human liver samples.

CYP2D6 is one of the most polymorphic drug-metabolizing enzymes in the liver. While genetic CYP2D6 variants serve as clinical biomarkers to predict CYP2D6 activity, large inter-person variability in CYP2D6 expression remains unaccounted for. Previous results suggest that there is variable expression of a CYP2D6 splice isoform with an in-frame deletion of exon 3 (CYP2D6ΔE3) encoding a protein lacking numerous active site residues. Here, using fragment analysis and RT-qPCR, we revealed that rs1058164 G (MAF = 27%-43%) is associated with increased formation of CYP2D6∆E3 in human liver samples (1.4-2.5-fold) and transfected cells. Furthermore, western blots showed that rs1058164 G was associated with a 50% decrease in full-length hepatic CYP2D6 protein expression. In addition, by studying a larger liver cohort, we confirmed our previous results that rs16947 (CYP2D6*2) reduces full-length CYP2D6 mRNA by increasing the production of an unstable splice isoform lacking exon 6 (CYP2D6ΔE6) and that the impact of CYP2D6ΔE6 is offset in carriers of the downstream enhancer variant rs5758550. The three frequent SNPs (rs1058164, rs16947, and rs5758550) form various 3-SNP-haplotypes, each with distinct CYP2D6 expression characteristics. Using an expression score (ES) system, we tested the impact of the 3-SNP-haplotype on improving the standard model to predict hepatic CYP2D6 protein expression based on genotype. A model that incorporates the 3-SNP-haplotype provided the best fit for CYP2D6 expression and also accounted for more variability in CYP2D6 protein levels (59%) than a model based on the accepted standard (36%) or one that only adds rs16947 and rs5758550 (42%). Clinical studies are needed to determine whether including the 3-SNP-haplotype alongside current standard CYP2D6 models improves the predictive value of CYP2D6 panels.

Oncotype DX results increase concordance in adjuvant chemotherapy recommendations for early-stage breast cancer.

Adjuvant chemotherapy recommendations for ER+/HER2- early-stage breast cancers (eBC) involve integrating prognostic and predictive information which rely on physician judgment; this can lead to discordant recommendations. In this study we aim to evaluate whether Oncotype DX improves confidence and agreement among oncologists in adjuvant chemotherapy recommendations. We randomly select 30 patients with ER+/HER2- eBC and recurrence score (RS) available from an institutional database. We ask 16 breast oncologists with varying years of clinical practice in Italy and the US to provide recommendation for the addition of chemotherapy to endocrine therapy and their degree of confidence in the recommendation twice; first, based on clinicopathologic features only (pre-RS), and then with RS result (post-RS). Pre-RS, the average rate of chemotherapy recommendation is 50.8% and is higher among junior (62% vs 44%; p < 0.001), but similar by country. Oncologists are uncertain in 39% of cases and recommendations are discordant in 27% of cases (interobserver agreement K 0.47). Post-RS, 30% of physicians change recommendation, uncertainty in recommendation decreases to 5.6%, and discordance decreases to 7% (interobserver agreement K 0.85). Interpretation of clinicopathologic features alone to recommend adjuvant chemotherapy results in 1 out of 4 discordant recommendations and relatively high physician uncertainty. Oncotype DX results decrease discordancy to 1 out of 15, and reduce physician uncertainty. Genomic assay results reduce subjectivity in adjuvant chemotherapy recommendations for ER +/HER2- eBC.

Genetic basis for probiotic yeast phenotypes revealed by nanopore sequencing.

Probiotic yeasts are emerging as preventative and therapeutic solutions for disease. Often ingested via cultured foods and beverages, they can survive the harsh conditions of the gastrointestinal tract and adhere to it, where they provide nutrients and inhibit pathogens like Candida albicans. Yet, little is known of the genomic determinants of these beneficial traits. To this end, we have sequenced 2 food-derived probiotic yeast isolates that mitigate fungal infections. We find that the first strain, KTP, is a strain of Saccharomyces cerevisiae within a small clade that lacks any apparent ancestry from common European/wine S. cerevisiae strains. Significantly, we show that S. cerevisiae KTP genes involved in general stress, pH tolerance, and adherence are markedly different from S. cerevisiae S288C but are similar to the commercial probiotic yeast species S. boulardii. This suggests that even though S. cerevisiae KTP and S. boulardii are from different clades, they may achieve probiotic effect through similar genetic mechanisms. We find that the second strain, ApC, is a strain of Issatchenkia occidentalis, one of the few of this family of yeasts to be sequenced. Because of the dissimilarity of its genome structure and gene organization, we infer that I. occidentalis ApC likely achieves a probiotic effect through a different mechanism than the Saccharomyces strains. Therefore, this work establishes a strong genetic link among probiotic Saccharomycetes, advances the genomics of Issatchenkia yeasts, and indicates that probiotic activity is not monophyletic and complimentary mixtures of probiotics could enhance health benefits beyond a single species.

YAP and TAZ couple osteoblast precursor mobilization to angiogenesis and mechanoregulated bone development.

Endochondral ossification requires coordinated mobilization of osteoblast precursors with blood vessels. During adult bone homeostasis, vessel adjacent osteoblast precursors respond to and are maintained by mechanical stimuli; however, the mechanisms by which these cells mobilize and respond to mechanical cues during embryonic development are unknown. Previously, we found that deletion of the mechanoresponsive transcriptional regulators, YAP and TAZ, from Osterix-expressing osteoblast precursors and their progeny caused perinatal lethality. Here, we show that embryonic YAP/TAZ signaling couples vessel-associated osteoblast precursor mobilization to angiogenesis in developing long bones. Osterix-conditional YAP/TAZ deletion impaired endochondral ossification in the primary ossification center but not intramembranous osteogenesis in the bone collar. Single-cell RNA sequencing revealed YAP/TAZ regulation of the angiogenic chemokine, Cxcl12, which was expressed uniquely in vessel-associated osteoblast precursors. YAP/TAZ signaling spatially coupled osteoblast precursors to blood vessels and regulated vascular morphogenesis and vessel barrier function. Further, YAP/TAZ signaling regulated vascular loop morphogenesis at the chondro-osseous junction to control hypertrophic growth plate remodeling. In human cells, mesenchymal stromal cell co-culture promoted 3D vascular network formation, which was impaired by stromal cell YAP/TAZ depletion, but rescued by recombinant CXCL12 treatment. Lastly, YAP and TAZ mediated mechanotransduction for load-induced osteogenesis in embryonic bone.

Endothelial SMAD1/5 signaling couples angiogenesis to osteogenesis during long bone growth.

Skeletal development depends on coordinated angiogenesis and osteogenesis. Bone morphogenetic proteins direct bone development by activating SMAD1/5 signaling in osteoblasts. However, the role of SMAD1/5 in skeletal endothelium is unknown. Here, we found that endothelial cell-conditional SMAD1/5 depletion in juvenile mice caused metaphyseal and diaphyseal hypervascularity, resulting in altered cancellous and cortical bone formation. SMAD1/5 depletion induced excessive sprouting, disrupting the columnar structure of the metaphyseal vessels and impaired anastomotic loop morphogenesis at the chondro-osseous junction. Endothelial SMAD1/5 depletion impaired growth plate resorption and, upon long term depletion, abrogated osteoprogenitor recruitment to the primary spongiosa. Finally, in the diaphysis, endothelial SMAD1/5 activity was necessary to maintain the sinusoidal phenotype, with SMAD1/5 depletion inducing formation of large vascular loops, featuring elevated endomucin expression, ectopic tip cell formation, and hyperpermeability. Together, endothelial SMAD1/5 activity sustains skeletal vascular morphogenesis and function and coordinates growth plate remodeling and osteoprogenitor recruitment dynamics during bone growth.

Health-related quality of life impacts upon 5-year survival after coronary artery bypass surgery.

BACKGROUND: Poor preoperative health-related quality of life (HRQoL) has been associated with reduced short-term survival after coronary artery bypass graft (CABG) surgery; however, its impact on long-term mortality is unknown. This study's objective was to determine if baseline HRQoL status predicts 5-year post-CABG mortality. METHODS: This prespecified, randomized on/off bypass follow-up study (ROOBY-FS) subanalysis compared baseline patient characteristics and HRQoL scores, obtained from the Seattle Angina Questionnaire (SAQ) and Veterans RAND Short Form-36 (VR-36), between 5-year post-CABG survivors and nonsurvivors. Standardized subscores were calculated for each questionnaire. Multivariable logistic regression assessed whether HRQoL survey subcomponents independently predicted 5-year mortality (p ≤ .05). RESULTS: Of the 2203 ROOBY-FS enrollees, 2104 (95.5%) completed baseline surveys. Significant differences between 5-year post-CABG deaths (n = 286) and survivors (n = 1818) included age, history of chronic obstructive pulmonary disease, stroke, peripheral vascular disease, renal dysfunction, diabetes, lower left ventricular ejection fraction, atrial fibrillation, depression, non-White race/ethnicity, lower education status, and off-pump CABG. Adjusting for these factors, baseline VR-36 physical component summary score (p = .01), VR-36 mental component summary score (p < .001), and SAQ physical limitation score (p = .003) were all associated with 5-year all-cause mortality. CONCLUSIONS: Pre-CABG HRQoL scores may provide clinically relevant prognostic information beyond traditional risk models and prove useful for patient-provider shared decision-making and enhancing pre-CABG informed consent.

Regulatory variants in a novel distal enhancer regulate the expression of CYP3A4 and CYP3A5.

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.

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.

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.