Search

C1q deficiency is a rare inborn error of immunity characterized by increased susceptibility to infections and autoimmune manifestations mimicking SLE, with an associated morbidity and mortality. Because C1q is synthesized by monocytes, to date, four patients treated with allogeneic HSCT have been reported, with a positive outcome in three. We conducted an international retrospective study to assess the outcome of HSCT in C1q deficiency. Eighteen patients, fourteen previously unreported, from eleven referral centres, were included. Two patients had two HSCTs, thus 20 HSCTs were performed in total, at a median age of 10 years (range 0.9-19). Indications for HSCT were autoimmune manifestations not controlled by ongoing treatment in seventeen, and early development of MALT lymphoma in one patient. Overall survival (OS) was 71% and event-free survival was 59% at two years (considering an event as acute GvHD ≥ grade III, disease recurrence and death). In eleven patients HSCT led to resolution of autoimmune features and discontinuation of immunosuppressive treatments (follow-up time range 3-84 months). Five patients died due to transplant-related complications. Patients with a severe autoimmune phenotype, defined as neurological and/or renal involvement, had the worst OS (40% vs 84%; p = 0.034). Reviewing data of 69 genetically confirmed C1q deficient patients, we found that anti-Ro antibodies are associated with neurologic involvement, and anti-RNP and anti-DNA antibodies with renal involvement. In conclusion, HSCT may be a valid curative option for C1q deficiency, but careful selection of patients, with an accurate assessment of risk and benefit, is mandatory.

BackgroundActivated phosphoinositide 3-kinase delta syndrome (APDS) is a combined immunodeficiency with a heterogeneous phenotype considered reversible by allogeneic hematopoietic cell transplantation (HCT).ObjectivesThis study sought to characterize HCT outcomes in APDS.MethodsRetrospective data were collected on 57 patients with APDS1/2 (median age, 13 years; range, 2-66 years) who underwent HCT.ResultsPre-HCT comorbidities such as lung, gastrointestinal, and liver pathology were common, with hematologic malignancy in 26%. With median follow-up of 2.3 years, 2-year overall and graft failure-free survival probabilities were 86% and 68%, respectively, and did not differ significantly by APDS1 versus APDS2, donor type, or conditioning intensity. The 2-year cumulative incidence of graft failure following first HCT was 17% overall but 42% if mammalian target of rapamycin inhibitor(s) (mTORi) were used in the first year post-HCT, compared with 9% without mTORi. Similarly, 2-year cumulative incidence of unplanned donor cell infusion was overall 28%, but 65% in the context of mTORi receipt and 23% without. Phenotype reversal occurred in 96% of evaluable patients, of whom 17% had mixed chimerism. Vulnerability to renal complications continued post-HCT, adding new insights into potential nonimmunologic roles of phosphoinositide 3-kinase not correctable through HCT.ConclusionsGraft failure, graft instability, and poor graft function requiring unplanned donor cell infusion were major barriers to successful HCT. Post-HCT mTORi use may confer an advantage to residual host cells, promoting graft instability. Longer-term post-HCT follow-up of more patients is needed to elucidate the kinetics of immune reconstitution and donor chimerism, establish approaches that reduce graft instability, and assess the completeness of phenotype reversal over time.

BackgroundCommon chromosome 9p21 single nucleotide polymorphisms (SNPs) increase coronary heart disease risk, independent of traditional lipid risk factors. However, lipids comprise large numbers of structurally related molecules not measured in traditional risk measurements, and many have inflammatory bioactivities. Here, we applied lipidomic and genomic approaches to 3 model systems to characterize lipid metabolic changes in common Chr9p21 SNPs, which confer ≈30% elevated coronary heart disease risk associated with altered expression of ANRIL, a long ncRNA.MethodsUntargeted and targeted lipidomics was applied to plasma from NPHSII (Northwick Park Heart Study II) homozygotes for AA or GG in rs10757274, followed by correlation and network analysis. To identify candidate genes, transcriptomic data from shRNA downregulation of ANRIL in HEK-293 cells was mined. Transcriptional data from vascular smooth muscle cells differentiated from induced pluripotent stem cells of individuals with/without Chr9p21 risk, nonrisk alleles, and corresponding knockout isogenic lines were next examined. Last, an in-silico analysis of miRNAs was conducted to identify how ANRIL might control lysoPL (lysophosphospholipid)/lysoPA (lysophosphatidic acid) genes.ResultsElevated risk GG correlated with reduced lysoPLs, lysoPA, and ATX (autotaxin). Five other risk SNPs did not show this phenotype. LysoPL-lysoPA interconversion was uncoupled from ATX in GG plasma, suggesting metabolic dysregulation. Significantly altered expression of several lysoPL/lysoPA metabolizing enzymes was found in HEK cells lacking ANRIL. In the vascular smooth muscle cells data set, the presence of risk alleles associated with altered expression of several lysoPL/lysoPA enzymes. Deletion of the risk locus reversed the expression of several lysoPL/lysoPA genes to nonrisk haplotype levels. Genes that were altered across both cell data sets were DGKA, MBOAT2, PLPP1, and LPL. The in-silico analysis identified 4 ANRIL-regulated miRNAs that control lysoPL genes as miR-186-3p, miR-34a-3p, miR-122-5p, and miR-34a-5p.ConclusionsA Chr9p21 risk SNP associates with complex alterations in immune-bioactive phospholipids and their metabolism. Lipid metabolites and genomic pathways associated with coronary heart disease pathogenesis in Chr9p21 and ANRIL-associated disease are demonstrated.

BackgroundCommon chromosome 9p21 single nucleotide polymorphisms (SNPs) increase coronary heart disease risk, independent of traditional lipid risk factors. However, lipids comprise large numbers of structurally related molecules not measured in traditional risk measurements, and many have inflammatory bioactivities. Here, we applied lipidomic and genomic approaches to 3 model systems to characterize lipid metabolic changes in common Chr9p21 SNPs, which confer ≈30% elevated coronary heart disease risk associated with altered expression of ANRIL, a long ncRNA.MethodsUntargeted and targeted lipidomics was applied to plasma from NPHSII (Northwick Park Heart Study II) homozygotes for AA or GG in rs10757274, followed by correlation and network analysis. To identify candidate genes, transcriptomic data from shRNA downregulation of ANRIL in HEK-293 cells was mined. Transcriptional data from vascular smooth muscle cells differentiated from induced pluripotent stem cells of individuals with/without Chr9p21 risk, nonrisk alleles, and corresponding knockout isogenic lines were next examined. Last, an in-silico analysis of miRNAs was conducted to identify how ANRIL might control lysoPL (lysophosphospholipid)/lysoPA (lysophosphatidic acid) genes.ResultsElevated risk GG correlated with reduced lysoPLs, lysoPA, and ATX (autotaxin). Five other risk SNPs did not show this phenotype. LysoPL-lysoPA interconversion was uncoupled from ATX in GG plasma, suggesting metabolic dysregulation. Significantly altered expression of several lysoPL/lysoPA metabolizing enzymes was found in HEK cells lacking ANRIL. In the vascular smooth muscle cells data set, the presence of risk alleles associated with altered expression of several lysoPL/lysoPA enzymes. Deletion of the risk locus reversed the expression of several lysoPL/lysoPA genes to nonrisk haplotype levels. Genes that were altered across both cell data sets were DGKA, MBOAT2, PLPP1, and LPL. The in-silico analysis identified 4 ANRIL-regulated miRNAs that control lysoPL genes as miR-186-3p, miR-34a-3p, miR-122-5p, and miR-34a-5p.ConclusionsA Chr9p21 risk SNP associates with complex alterations in immune-bioactive phospholipids and their metabolism. Lipid metabolites and genomic pathways associated with coronary heart disease pathogenesis in Chr9p21 and ANRIL-associated disease are demonstrated.

BackgroundImmunodysregulation polyendocrinopathy enteropathy x-linked (IPEX) syndrome is a monogenic autoimmune disease caused by FOXP3 mutations. Because it is a rare disease, the natural history and response to treatments, including allogeneic hematopoietic stem cell transplantation (HSCT) and immunosuppression (IS), have not been thoroughly examined.ObjectiveThis analysis sought to evaluate disease onset, progression, and long-term outcome of the 2 main treatments in long-term IPEX survivors.MethodsClinical histories of 96 patients with a genetically proven IPEX syndrome were collected from 38 institutions worldwide and retrospectively analyzed. To investigate possible factors suitable to predict the outcome, an organ involvement (OI) scoring system was developed.ResultsWe confirm neonatal onset with enteropathy, type 1 diabetes, and eczema. In addition, we found less common manifestations in delayed onset patients or during disease evolution. There is no correlation between the site of mutation and the disease course or outcome, and the same genotype can present with variable phenotypes. HSCT patients (n = 58) had a median follow-up of 2.7 years (range, 1 week-15 years). Patients receiving chronic IS (n = 34) had a median follow-up of 4 years (range, 2 months-25 years). The overall survival after HSCT was 73.2% (95% CI, 59.4-83.0) and after IS was 65.1% (95% CI, 62.8-95.8). The pretreatment OI score was the only significant predictor of overall survival after transplant (P = .035) but not under IS.ConclusionsPatients receiving chronic IS were hampered by disease recurrence or complications, impacting long-term disease-free survival. When performed in patients with a low OI score, HSCT resulted in disease resolution with better quality of life, independent of age, donor source, or conditioning regimen.

BACKGROUND:Immunodysregulation polyendocrinopathy enteropathy x-linked (IPEX) syndrome is a monogenic autoimmune disease caused by FOXP3 mutations. Because it is a rare disease, the natural history and response to treatments, including allogeneic hematopoietic stem cell transplantation (HSCT) and immunosuppression (IS), have not been thoroughly examined. OBJECTIVE:This analysis sought to evaluate disease onset, progression, and long-term outcome of the 2 main treatments in long-term IPEX survivors. METHODS:Clinical histories of 96 patients with a genetically proven IPEX syndrome were collected from 38 institutions worldwide and retrospectively analyzed. To investigate possible factors suitable to predict the outcome, an organ involvement (OI) scoring system was developed. RESULTS:We confirm neonatal onset with enteropathy, type 1 diabetes, and eczema. In addition, we found less common manifestations in delayed onset patients or during disease evolution. There is no correlation between the site of mutation and the disease course or outcome, and the same genotype can present with variable phenotypes. HSCT patients (n = 58) had a median follow-up of 2.7 years (range, 1 week-15 years). Patients receiving chronic IS (n = 34) had a median follow-up of 4 years (range, 2 months-25 years). The overall survival after HSCT was 73.2% (95% CI, 59.4-83.0) and after IS was 65.1% (95% CI, 62.8-95.8). The pretreatment OI score was the only significant predictor of overall survival after transplant (P = .035) but not under IS. CONCLUSIONS:Patients receiving chronic IS were hampered by disease recurrence or complications, impacting long-term disease-free survival. When performed in patients with a low OI score, HSCT resulted in disease resolution with better quality of life, independent of age, donor source, or conditioning regimen.

BackgroundRare DNA breakage repair disorders predispose to infection and lymphoreticular malignancies. Hematopoietic cell transplantation (HCT) is curative, but coadministered chemotherapy or radiotherapy is damaging because of systemic radiosensitivity. We collected HCT outcome data for Nijmegen breakage syndrome, DNA ligase IV deficiency, Cernunnos-XRCC4-like factor (Cernunnos-XLF) deficiency, and ataxia-telangiectasia (AT).MethodsData from 38 centers worldwide, including indication, donor, conditioning regimen, graft-versus-host disease, and outcome, were analyzed. Conditioning was classified as myeloablative conditioning (MAC) if it contained radiotherapy or alkylators and reduced-intensity conditioning (RIC) if no alkylators and/or 150 mg/m2 fludarabine or less and 40 mg/kg cyclophosphamide or less were used.ResultsFifty-five new, 14 updated, and 18 previously published patients were analyzed. Median age at HCT was 48 months (range, 1.5-552 months). Twenty-nine patients underwent transplantation for infection, 21 had malignancy, 13 had bone marrow failure, 13 received pre-emptive transplantation, 5 had multiple indications, and 6 had no information. Twenty-two received MAC, 59 received RIC, and 4 were infused; information was unavailable for 2 patients. Seventy-three of 77 patients with DNA ligase IV deficiency, Cernunnos-XLF deficiency, or Nijmegen breakage syndrome received conditioning. Survival was 53 (69%) of 77 and was worse for those receiving MAC than for those receiving RIC (P = .006). Most deaths occurred early after transplantation, suggesting poor tolerance of conditioning. Survival in patients with AT was 25%. Forty-one (49%) of 83 patients experienced acute GvHD, which was less frequent in those receiving RIC compared with those receiving MAC (26/56 [46%] vs 12/21 [57%], P = .45). Median follow-up was 35 months (range, 2-168 months). No secondary malignancies were reported during 15 years of follow-up. Growth and developmental delay remained after HCT; immune-mediated complications resolved.ConclusionRIC HCT resolves DNA repair disorder-associated immunodeficiency. Long-term follow-up is required for secondary malignancy surveillance. Routine HCT for AT is not recommended.

Severe combined immunodeficiency (SCID) is effectively treated with hematopoietic cell transplantation (HCT), with overall survival approaching 90% in contemporary reports. However, survivors are at risk for developing late complications because of the variable durability of high-quality immune function, underlying genotype of SCID, comorbidities due to infections in the pretransplantation and post-transplantation periods, and use of conditioning before transplantation. An international group of transplantation experts was convened in 2016 to review the current knowledge of late effects seen in SCID patients after HCT and to develop recommendations for screening and monitoring for late effects. This report provides recommendations for screening and management of pediatric and adult SCID patients treated with HCT.

BackgroundX-linked hyper-IgM syndrome (XHIGM) is a primary immunodeficiency with high morbidity and mortality compared with those seen in healthy subjects. Hematopoietic cell transplantation (HCT) has been considered a curative therapy, but the procedure has inherent complications and might not be available for all patients.ObjectivesWe sought to collect data on the clinical presentation, treatment, and follow-up of a large sample of patients with XHIGM to (1) compare long-term overall survival and general well-being of patients treated with or without HCT along with clinical factors associated with mortality and (2) summarize clinical practice and risk factors in the subgroup of patients treated with HCT.MethodsPhysicians caring for patients with primary immunodeficiency diseases were identified through the Jeffrey Modell Foundation, United States Immunodeficiency Network, Latin American Society for Immunodeficiency, and Primary Immune Deficiency Treatment Consortium. Data were collected with a Research Electronic Data Capture Web application. Survival from time of diagnosis or transplantation was estimated by using the Kaplan-Meier method compared with log-rank tests and modeled by using proportional hazards regression.ResultsTwenty-eight clinical sites provided data on 189 patients given a diagnosis of XHIGM between 1964 and 2013; 176 had valid follow-up and vital status information. Sixty-seven (38%) patients received HCT. The average follow-up time was 8.5 ± 7.2 years (range, 0.1-36.2 years). No difference in overall survival was observed between patients treated with or without HCT (P = .671). However, risk associated with HCT decreased for diagnosis years 1987-1995; the hazard ratio was significantly less than 1 for diagnosis years 1995-1999. Liver disease was a significant predictor of overall survival (hazard ratio, 4.9; 95% confidence limits, 2.2-10.8; P

BackgroundIntravenous busulfan combined with therapeutic drug monitoring to guide dosing improves outcomes after allogeneic haemopoietic cell transplantation (HCT). The best method to estimate busulfan exposure and optimum exposure in children or young adults remains unclear. We therefore assessed three approaches to estimate intravenous busulfan exposure (expressed as cumulative area under the curve [AUC]) and associated busulfan AUC with clinical outcomes in children or young adults undergoing allogeneic HCT.MethodsIn this retrospective analysis, patients from 15 centres in the Netherlands, USA, Canada, Switzerland, UK, Italy, Germany, and Australia who received a busulfan-based conditioning regimen between March 18, 2001, and Feb 12, 2015, were included. Cumulative AUC was calculated by numerical integration using non-linear mixed effect modelling (AUCNONMEM), non-compartmental analysis (AUC from 0 to infinity [AUC0-∞] and to the next dose [AUC0-τ]), and by individual centres using various approaches (AUCcentre). The main outcome of interest was event-free survival. Other outcomes of interest were graft failure or relapse, or both; transplantation-related mortality; acute toxicity (veno-occlusive disease or acute graft versus-host disease [GvHD]); chronic GvHD; overall survival; and chronic-GvHD-free event-free survival. We used propensity-score-adjusted Cox proportional hazard models, Weibull models, and Fine-Gray competing risk regressions for statistical analyses.Findings790 patients were enrolled, 674 of whom were included: 274 (41%) with malignant and 400 (59%) with non-malignant disease. Median age was 4·5 years (IQR 1·4-10·7). The median busulfan AUCNONMEM was 74·4 mg × h/L (95% CI 31·1-104·6), which correlated with the standardised method AUC0-∞ (r2=0·74), but the latter correlated poorly with AUCcentre (r2=0·35). Estimated 2-year event-free survival was 69·7% (95% CI 66·2-73·0). Event-free survival at 2 years was 77·0% (95% CI 72·1-82·9) in the 257 patients with an optimum intravenous busulfan AUC of 78-101 mg × h/L compared with 66·1% (60·9-71·4) in the 235 patients at the low historical target of 58-86 mg × h/L and 49·5% (29·2-66·0) in the 44 patients with a high (>101 mg × h/L) busulfan AUC (p=0·011). Compared with the low AUC group, graft failure or relapse occurred less frequently in the optimum AUC group (hazard ratio [HR] 0·57, 95% CI 0·39-0·84; p=0·0041). Acute toxicity (HR 1·69, 1·12-2·57; p=0·013) and transplantation-related mortality (2·99, 1·82-4·92; p101 mg × h/L) than in the low AUC group (

BackgroundPatients with severe combined immunodeficiency disease who have matched sibling donors (MSDs) can proceed to hematopoietic cell transplantation (HCT) without conditioning chemotherapy.ObjectiveWe sought to determine whether the results of HCT without chemotherapy-based conditioning from matched unrelated donors (URDs), either from volunteer adults or umbilical cord blood, are comparable with those from MSDs.MethodsWe performed a multicenter survey of severe combined immunodeficiency transplantation centers in North America, Europe, and Australia to compile retrospective data on patients who have undergone unconditioned HCT from either URDs (n = 37) or MSDs (n = 66).ResultsMost patients undergoing URD HCT (92%) achieved donor T-cell engraftment compared with 97% for those with MSDs; however, estimated 5-year overall and event-free survival were worse for URD recipients (71% and 60%, respectively) compared with MSD recipients (92% and 89%, respectively; P < .01 for both). URD recipients who received pre-HCT serotherapy had similar 5-year overall survival (100%) to MSD recipients. The incidences of grade II to IV acute and chronic graft-versus-host disease were higher in URD (50% and 39%, respectively) compared with MSD (22% and 5%, respectively) recipients (P < .01 for both). In the surviving patients there was no difference in T-cell reconstitution at the last follow-up between the URD and MSD recipients; however, MSD recipients were more likely to achieve B-cell reconstitution (72% vs 17%, P < .001).ConclusionUnconditioned URD HCT achieves excellent rates of donor T-cell engraftment similar to that seen in MSD recipients, and reconstitution rates are adequate. However, only a minority will have myeloid and B-cell reconstitution, and attention must be paid to graft-versus-host disease prophylaxis. This approach might be safer in children ineligible for intense regimens to spare the potential complications of chemotherapy.

BackgroundPatients with severe combined immunodeficiency disease who have matched sibling donors (MSDs) can proceed to hematopoietic cell transplantation (HCT) without conditioning chemotherapy.ObjectiveWe sought to determine whether the results of HCT without chemotherapy-based conditioning from matched unrelated donors (URDs), either from volunteer adults or umbilical cord blood, are comparable with those from MSDs.MethodsWe performed a multicenter survey of severe combined immunodeficiency transplantation centers in North America, Europe, and Australia to compile retrospective data on patients who have undergone unconditioned HCT from either URDs (n = 37) or MSDs (n = 66).ResultsMost patients undergoing URD HCT (92%) achieved donor T-cell engraftment compared with 97% for those with MSDs; however, estimated 5-year overall and event-free survival were worse for URD recipients (71% and 60%, respectively) compared with MSD recipients (92% and 89%, respectively; P < .01 for both). URD recipients who received pre-HCT serotherapy had similar 5-year overall survival (100%) to MSD recipients. The incidences of grade II to IV acute and chronic graft-versus-host disease were higher in URD (50% and 39%, respectively) compared with MSD (22% and 5%, respectively) recipients (P < .01 for both). In the surviving patients there was no difference in T-cell reconstitution at the last follow-up between the URD and MSD recipients; however, MSD recipients were more likely to achieve B-cell reconstitution (72% vs 17%, P < .001).ConclusionUnconditioned URD HCT achieves excellent rates of donor T-cell engraftment similar to that seen in MSD recipients, and reconstitution rates are adequate. However, only a minority will have myeloid and B-cell reconstitution, and attention must be paid to graft-versus-host disease prophylaxis. This approach might be safer in children ineligible for intense regimens to spare the potential complications of chemotherapy.

BACKGROUND: Severe combined immunodeficiency (SCID) is a syndrome characterized by profound T-cell deficiency. BCG vaccine is contraindicated in patients with SCID. Because most countries encourage BCG vaccination at birth, a high percentage of patients with SCID are vaccinated before their immune defect is detected. OBJECTIVES: We sought to describe the complications and risks associated with BCG vaccination in patients with SCID. METHODS: An extensive standardized questionnaire evaluating complications, therapeutics, and outcomes regarding BCG vaccination in patients given a diagnosis of SCID was widely distributed. Summary statistics and association analysis was performed. RESULTS: Data on 349 BCG-vaccinated patients with SCID from 28 centers in 17 countries were analyzed. Fifty-one percent of the patients had BCG-associated complications, 34% disseminated and 17% localized (a 33,000- and 400-fold increase, respectively, over the general population). Patients receiving early vaccination (≤1 month) showed an increased prevalence of complications (P = .006) and death caused by BCG-associated complications (P < .0001). The odds of experiencing complications among patients with T-cell numbers of 250/μL or less at diagnosis was 2.1 times higher (95% CI, 1.4-3.4 times higher; P = .001) than among those with T-cell numbers of greater than 250/μL. BCG-associated complications were reported in 2 of 78 patients who received antimycobacterial therapy while asymptomatic, and no deaths caused by BCG-associated complications occurred in this group. In contrast, 46 BCG-associated deaths were reported among 160 patients treated with antimycobacterial therapy for a symptomatic BCG infection (P < .0001). CONCLUSIONS: BCG vaccine has a very high rate of complications in patients with SCID, which increase morbidity and mortality rates. Until safer and more efficient antituberculosis vaccines become available, delay in BCG vaccination should be considered to protect highly vulnerable populations from preventable complications.

Little information is currently available regarding the pharmacokinetics (PK) of busulfan in infants and small children to help guide decisions for safe and efficacious drug therapy. The objective of this study was to develop an algorithm for individualized dosing of i.v. busulfan in infants and children weighing ≤12 kg, that would achieve targeted exposure with the first dose of busulfan. Population PK modeling was conducted using intensive time-concentration data collected through the routine therapeutic drug monitoring of busulfan in 149 patients from 8 centers. Busulfan PK was well described by a 1-compartment base model with linear elimination. The important clinical covariates affecting busulfan PK were actual body weight and age. Based on our model, the predicted clearance of busulfan increases approximately 1.7-fold between 6 weeks to 2 years of life. For infants age

Little information is currently available regarding the pharmacokinetics (PK) of busulfan in infants and small children to help guide decisions for safe and efficacious drug therapy. The objective of this study was to develop an algorithm for individualized dosing of i.v. busulfan in infants and children weighing ≤12 kg, that would achieve targeted exposure with the first dose of busulfan. Population PK modeling was conducted using intensive time-concentration data collected through the routine therapeutic drug monitoring of busulfan in 149 patients from 8 centers. Busulfan PK was well described by a 1-compartment base model with linear elimination. The important clinical covariates affecting busulfan PK were actual body weight and age. Based on our model, the predicted clearance of busulfan increases approximately 1.7-fold between 6 weeks to 2 years of life. For infants age

Little information is currently available regarding the pharmacokinetics (PK) of busulfan in infants and small children to help guide decisions for safe and efficacious drug therapy. The objective of this study was to develop an algorithm for individualized dosing of i.v. busulfan in infants and children weighing ≤12 kg, that would achieve targeted exposure with the first dose of busulfan. Population PK modeling was conducted using intensive time-concentration data collected through the routine therapeutic drug monitoring of busulfan in 149 patients from 8 centers. Busulfan PK was well described by a 1-compartment base model with linear elimination. The important clinical covariates affecting busulfan PK were actual body weight and age. Based on our model, the predicted clearance of busulfan increases approximately 1.7-fold between 6 weeks to 2 years of life. For infants age

Liver cytochrome P450s (P450s) play critical roles in drug metabolism, toxicology, and metabolic processes. Despite rapid progress in the understanding of these enzymes, a systematic investigation of the full spectrum of functionality of individual P450s, the interrelationship or networks connecting them, and the genetic control of each gene/enzyme is lacking. To this end, we genotyped, expression-profiled, and measured P450 activities of 466 human liver samples and applied a systems biology approach via the integration of genetics, gene expression, and enzyme activity measurements. We found that most P450s were positively correlated among themselves and were highly correlated with known regulators as well as thousands of other genes enriched for pathways relevant to the metabolism of drugs, fatty acids, amino acids, and steroids. Genome-wide association analyses between genetic polymorphisms and P450 expression or enzyme activities revealed sets of SNPs associated with P450 traits, and suggested the existence of both cis-regulation of P450 expression (especially for CYP2D6) and more complex trans-regulation of P450 activity. Several novel SNPs associated with CYP2D6 expression and enzyme activity were validated in an independent human cohort. By constructing a weighted coexpression network and a Bayesian regulatory network, we defined the human liver transcriptional network structure, uncovered subnetworks representative of the P450 regulatory system, and identified novel candidate regulatory genes, namely, EHHADH, SLC10A1, and AKR1D1. The P450 subnetworks were then validated using gene signatures responsive to ligands of known P450 regulators in mouse and rat. This systematic survey provides a comprehensive view of the functionality, genetic control, and interactions of P450s.

Liver cytochrome P450s (P450s) play critical roles in drug metabolism, toxicology, and metabolic processes. Despite rapid progress in the understanding of these enzymes, a systematic investigation of the full spectrum of functionality of individual P450s, the interrelationship or networks connecting them, and the genetic control of each gene/enzyme is lacking. To this end, we genotyped, expression-profiled, and measured P450 activities of 466 human liver samples and applied a systems biology approach via the integration of genetics, gene expression, and enzyme activity measurements. We found that most P450s were positively correlated among themselves and were highly correlated with known regulators as well as thousands of other genes enriched for pathways relevant to the metabolism of drugs, fatty acids, amino acids, and steroids. Genome-wide association analyses between genetic polymorphisms and P450 expression or enzyme activities revealed sets of SNPs associated with P450 traits, and suggested the existence of both cis-regulation of P450 expression (especially for CYP2D6) and more complex trans-regulation of P450 activity. Several novel SNPs associated with CYP2D6 expression and enzyme activity were validated in an independent human cohort. By constructing a weighted coexpression network and a Bayesian regulatory network, we defined the human liver transcriptional network structure, uncovered subnetworks representative of the P450 regulatory system, and identified novel candidate regulatory genes, namely, EHHADH, SLC10A1, and AKR1D1. The P450 subnetworks were then validated using gene signatures responsive to ligands of known P450 regulators in mouse and rat. This systematic survey provides a comprehensive view of the functionality, genetic control, and interactions of P450s.