Regulatory elements in SEM1-DLX5-DLX6 (7q21.3) locus contribute to genetic control of coronal nonsyndromic craniosynostosis and bone density-related traits.

Purpose: The etiopathogenesis of coronal nonsyndromic craniosynostosis (cNCS), a congenital condition defined by premature fusion of 1 or both coronal sutures, remains largely unknown. Methods: We conducted the largest genome-wide association study of cNCS followed by replication, fine mapping, and functional validation of the most significant region using zebrafish animal model. Results: Genome-wide association study identified 6 independent genome-wide-significant risk alleles, 4 on chromosome 7q21.3 SEM1-DLX5-DLX6 locus, and their combination conferred over 7-fold increased risk of cNCS. The top variants were replicated in an independent cohort and showed pleiotropic effects on brain and facial morphology and bone mineral density. Fine mapping of 7q21.3 identified a craniofacial transcriptional enhancer (eDlx36) within the linkage region of the top variant (rs4727341; odds ratio [95% confidence interval], 0.48[0.39-0.59]; P = 1.2E-12) that was located in SEM1 intron and enriched in 4 rare risk variants. In zebrafish, the activity of the transfected human eDlx36 enhancer was observed in the frontonasal prominence and calvaria during skull development and was reduced when the 4 rare risk variants were introduced into the sequence. Conclusion: Our findings support a polygenic nature of cNCS risk and functional role of craniofacial enhancers in cNCS susceptibility with potential broader implications for bone health.

Pharmacogenetics in Italy: current landscape and future prospects.

Pharmacogenetics investigates sequence of genes that affect drug response, enabling personalized medication. This approach reduces drug-induced adverse reactions and improves clinical effectiveness, making it a crucial consideration for personalized medical care. Numerous guidelines, drawn by global consortia and scientific organizations, codify genotype-driven administration for over 120 active substances. As the scientific community acknowledges the benefits of genotype-tailored therapy over traditionally agnostic drug administration, the push for its implementation into Italian healthcare system is gaining momentum. This evolution is influenced by several factors, including the improved access to patient genotypes, the sequencing costs decrease, the growing of large-scale genetic studies, the rising popularity of direct-to-consumer pharmacogenetic tests, and the continuous improvement of pharmacogenetic guidelines. Since EMA (European Medicines Agency) and AIFA (Italian Medicines Agency) provide genotype information on drug leaflet without clear and explicit clinical indications for gene testing, the regulation of pharmacogenetic testing is a pressing matter in Italy. In this manuscript, we have reviewed how to overcome the obstacles in implementing pharmacogenetic testing in the clinical practice of the Italian healthcare system. Our particular emphasis has been on germline testing, given the absence of well-defined national directives in contrast to somatic pharmacogenetics.

Drug-Induced Liver Injury in Pregnancy: The U.S. Drug-Induced Liver Injury Network Experience.

There are limited data on the causative agents and characteristics of drug-induced liver injury in pregnant individuals. Data from patients with drug-induced liver injury enrolled in the ongoing multicenter Drug-Induced Liver Injury Network between 2004 and 2022 and occurring during pregnancy or 6 months postpartum were reviewed and compared with cases of drug-induced liver injury in nonpregnant women of childbearing age. Among 325 individuals of childbearing age in the Drug-Induced Liver Injury Network, 16 cases of drug-induced liver injury (5%) occurred during pregnancy or postpartum. Compared with drug-induced liver injury in nonpregnant women, pregnancy-related drug-induced liver injury was more severe ( P <.05). One elective termination and three miscarriages were documented; there were no maternal deaths. We recommend that isoniazid for latent tuberculosis be deferred to the postpartum period whenever feasible and that ß-blockers or calcium channel blockers rather than methyldopa be used for hypertension management during pregnancy.

Drug-induced liver injury (DILI) ascribed to non-steroidal anti-inflammatory drugs (NSAIDs) in the USA-Update with genetic correlations.

OBJECTIVE: To describe patients with NSAID-DILI, including genetic factors associated with idiosyncratic DILI. METHODS: In DILIN, subjects with presumed DILI are enrolled and followed for at least 6 months. Causality is adjudicated by a Delphic approach. HLA sequencing of multiethnic NSAID-DILI patients and HLA allele imputation of matching population controls were performed following overall, class and drug-based association analysis. Significant results were tested in a non-Hispanic White (NHW) case-control replication cohort. RESULTS: Between September 2004 and March 2022, causality was adjudicated in 2498, and 55 (41 [75%] women) were assessed as likely due to NSAIDs. Median age at onset was 55 y (range 22-83 y). Diclofenac was the causative drug in 29, celecoxib in 7, ibuprofen in 5, etodolac and meloxicam each in 4. Except for meloxicam and oxaprozin (n = 2), the liver injury was hepatocellular with median R 15-25. HLA-DRB1*04:03 and HLA-B*35:03 were significantly more frequent in NSAID-DILI patients than in non-NSAID DILI controls. Interestingly, 85% of the HLA-DRB1*04:03 carriers developed DILI due to the use of acetic acid derivative NSAIDs, supporting the hypothesis that HLA-DRB1*04:03 could be a drug and/or class risk factor. HLA-B*35:03 but not HLA-DRB1*04:03 association was confirmed in the independent NHW replication cohort, which was largely driven by diclofenac. CONCLUSIONS: Despite prevalent use, NSAID-DILI is infrequent in the United States. Diclofenac is the most commonly implicated, and adherence to warnings of risk and close observation are recommended. The increased frequency of HLA-B*35:03 and DRB1*04:03, driven by diclofenac, suggests the importance of immune-mediated responses.

HLA-B*53:01 Is a Significant Risk Factor of Liver Injury due to Phenytoin and Other Antiepileptic Drugs in African Americans.

INTRODUCTION: To investigate human leukocyte antigen alleles associated with liver injury due to antiepileptic drugs (AEDs) in African Americans (AA). METHODS: In this study, 21 AA with AED drug-induced liver injury (DILI), 176 AA with DILI due to non-AEDs, and 5816 AA population controls were included. RESULTS: HLA-B*53:01 was significantly associated with aromatic AED-DILI (odds ratio: 4.52, 95% confidence interval: 2.42-8.44, P = 1.46 × 10 -5 ). Phenytoin DILI showed the strongest association with HLA-B*53:01 (odds ratio: 9.17; 95% confidence interval: 3.61-23.28, P = 1.1 × 10 -5 ). The HLA-B*53:01 allele was carried by 8 of 9 AA phenytoin DILI cases. DISCUSSION: HLA-B*53:01 is a significant risk factor of liver injury due to antiepileptics, particularly phenytoin, in AA.

Clinical Characteristics and Outcomes of Drug-Induced Acute Kidney Injury Cases.

Introduction: Drug-induced acute kidney injury (DI-AKI) is a frequent adverse event. The identification of DI-AKI is challenged by competing etiologies, clinical heterogeneity among patients, and a lack of accurate diagnostic tools. Our research aims to describe the clinical characteristics and predictive variables of DI-AKI. Methods: We analyzed data from the Drug-Induced Renal Injury Consortium (DIRECT) study (NCT02159209), an international, multicenter, observational cohort study of enriched clinically adjudicated DI-AKI cases. Cases met the primary inclusion criteria if the patient was exposed to at least 1 nephrotoxic drug for a minimum of 24 hours prior to AKI onset. Cases were clinically adjudicated, and inter-rater reliability (IRR) was measured using Krippendorff's alpha. Variables associated with DI-AKI were identified using L1 regularized multivariable logistic regression. Model performance was assessed using the area under the receiver operating characteristic curve (ROC AUC). Results: A total of 314 AKI cases met the eligibility criteria for this analysis, and 271 (86%) cases were adjudicated as DI-AKI. The majority of the AKI cases were recruited from the United States (68%). The most frequent causal nephrotoxic drugs were vancomycin (48.7%), nonsteroidal antiinflammatory drugs (18.2%), and piperacillin/tazobactam (17.8%). The IRR for DI-AKI adjudication was 0.309. The multivariable model identified age, vascular capacity, hyperglycemia, infections, pyuria, serum creatinine (SCr) trends, and contrast media as significant predictors of DI-AKI with good performance (ROC AUC 0.86). Conclusion: The identification of DI-AKI is challenging even with comprehensive adjudication by experienced nephrologists. Our analysis identified key clinical characteristics and outcomes of DI-AKI compared to other AKI etiologies.

Identification of Reduced ERAP2 Expression and a Novel HLA Allele as Components of a Risk Score for Susceptibility to Liver Injury Due to Amoxicillin-Clavulanate.

BACKGROUND & AIMS: Drug-induced liver injury (DILI) due to amoxicillin-clavulanate (AC) has been associated with HLA-A∗02:01, HLA-DRB1∗15:01, and rs2476601, a missense variant in PTPN22. The aim of this study was to identify novel risk factors for AC-DILI and to construct a genetic risk score (GRS). METHODS: Transcriptome-wide association study and genome-wide association study analyses were performed on 444 AC-DILI cases and 10,397 population-based controls of European descent. Associations were confirmed in a validation cohort (n = 133 cases and 17,836 population-based controls). Discovery and validation AC-DILI cases were also compared with 1358 and 403 non-AC-DILI cases. RESULTS: Transcriptome-wide association study revealed a significant association of AC-DILI risk with reduced liver expression of ERAP2 (P = 3.7 × 10-7), coding for an aminopeptidase involved in antigen presentation. The lead eQTL single nucleotide polymorphism, rs1363907 (G), was associated with AC-DILI risk in the discovery (odds ratio [OR], 1.68; 95% CI, 1.23-1.66; P = 1.7 × 10-7) and validation cohorts (OR, 1.2; 95% CI, 1.04-2.05; P = .03), following a recessive model. We also identified HLA-B∗15:18 as a novel AC-DILI risk factor in both discovery (OR, 4.19; 95% CI, 2.09-8.36; P = 4.9 × 10-5) and validation (OR, 7.78; 95% CI, 2.75-21.99; P = .0001) cohorts. GRS, incorporating rs1363907, rs2476601, HLA-B∗15:18, HLA-A∗02:01, and HLA-DRB1∗15:01, was highly predictive of AC-DILI risk when cases were analyzed against both general population and non-AC-DILI control cohorts. GRS was the most significant predictor in a regression model containing known AC-DILI clinical risk characteristics and significantly improved the predictive model. CONCLUSIONS: We identified novel associations of AC-DILI risk with ERAP2 low expression and with HLA-B∗15:18. GRS based on the 5 risk variants may assist AC-DILI causality assessment and risk management.

Clinical features, outcomes, and HLA risk factors associated with nitrofurantoin-induced liver injury.

BACKGROUND & AIMS: Nitrofurantoin (NTF) is widely used for the treatment (short-term) and prevention (long-term) of urinary tract infections. We aimed to describe the clinical characteristics, outcomes, and HLA risk factors for NTF-induced liver injury (NTF-DILI) among individuals enrolled in the Drug Induced Liver Injury Network (DILIN). METHODS: Seventy-eight individuals with definite, highly likely, or probable NTF-DILI were enrolled into DILIN studies between 2004-2020. HLA alleles were compared between NTF-DILI and three control groups: population (n = 14,001), idiopathic autoimmune hepatitis (n = 231), and non-NTF DILI (n = 661). RESULTS: Liver injury was hepatocellular in 69% and icteric in 55%. AST > ALT was more common in the 44 long-exposure (≥1 year) NTF-DILI cases than in the 18 short (≤7 days) and 16 intermediate (>7 to <365 days) exposure cases (73% vs. 33% vs. 50%, respectively, p = 0.018), as was ANA or SMA positivity (91% vs. 44% vs. 50%, respectively, p <0.001), and corticosteroid use (61% vs. 27% vs. 44%, respectively, p = 0.06). In long-term NTF-DILI, bridging fibrosis, nodularity or cirrhosis, or clinical and imaging evidence for cirrhosis were present in 38%, with massive or sub-massive necrosis in 20%. No one in the short-term exposure group died or underwent transplantation, whereas 7 (12%) patients from the other groups died or underwent transplantation. After covariate adjustments, HLA-DRB1∗11:04 was significantly more frequent in NTF-DILI compared to population controls (odds ratio [OR] 4.29, p = 1.15 × 10-4), idiopathic autoimmune hepatitis (OR 11.77, p = 7.76 × 10-5), and non-NTF DILI (OR 3.34, p = 0.003). CONCLUSION: NTF-DILI can result in parenchymal necrosis, bridging fibrosis, cirrhosis, and death or liver transplantation, especially with long-term exposure, and is associated with HLA-DRB1∗11:04. To mitigate against serious liver injury associated with NTF, regulators should revise the prescribing information and consider other mitigation strategies. IMPACT AND IMPLICATIONS: Nitrofurantoin is a recognized cause of drug-induced liver injury (DILI). In this study consisting of a large cohort of well-phenotyped individuals with nitrofurantoin-induced liver injury, two distinct patterns of liver injury were identified: liver injury associated with short-term exposure, which is generally self-limiting, and liver injury associated with long-term exposure, which can lead to advanced fibrosis, cirrhosis and liver failure. HLA DRB1∗11:04 is a risk factor for liver injury due to long-term nitrofurantoin exposure. Our findings are important for regulators as well as physicians prescribing and pharmacists dispensing nitrofurantoin.

Long-read HiFi sequencing of NUDT15: Phased full-gene haplotyping and pharmacogenomic allele discovery.

To determine the phase of NUDT15 sequence variants for more comprehensive star (*) allele diplotyping, we developed a novel long-read single-molecule real-time HiFi amplicon sequencing method. A 10.5 kb NUDT15 amplicon assay was validated using reference material positive controls and additional samples for specimen type and blinded accuracy assessment. Triplicate NUDT15 HiFi sequencing of two reference material samples had nonreference genotype concordances of >99.9%, indicating that the assay is robust. Notably, short-read genome sequencing of a subset of samples was unable to determine the phase of star (*) allele-defining NUDT15 variants, resulting in ambiguous diplotype results. In contrast, long-read HiFi sequencing phased all variants across the NUDT15 amplicons, including a *2/*9 diplotype that previously was characterized as *1/*2 in the 1000 Genomes Project v3 data set. Assay throughput was also tested using 8.5 kb amplicons from 100 Ashkenazi Jewish individuals, which identified a novel NUDT15 *1 suballele (c.-121G>A) and a rare likely deleterious coding variant (p.Pro129Arg). Both novel alleles were Sanger confirmed and assigned as *1.007 and *20, respectively, by the PharmVar Consortium. Taken together, NUDT15 HiFi amplicon sequencing is an innovative method for phased full-gene characterization and novel allele discovery, which could improve NUDT15 pharmacogenomic testing and subsequent phenotype prediction.

Urinary Exosomes Identify Inflammatory Pathways in Vancomycin Associated Acute Kidney Injury.

BACKGROUND: Vancomycin is commonly used as a first line therapy for gram positive organisms such as methicillin resistant Staphylococcusaureus. Vancomycin-induced acute kidney injury (V-AKI) has been reported in up to 43% of patients, especially in those with higher targeted trough concentrations. The precise mechanism of injury in humans remains elusive, with recent evidence directed towards proximal tubule cell apoptosis. In this study, we investigated the protein contents of urinary exosomes in patients with V-AKI to further elucidate biomarkers of mechanisms of injury and potential responses. METHODS: Urine samples from patients with V-AKI who were enrolled in the DIRECT study and matched healthy controls from the UAB-UCSD O'Brien Center Biorepository were included in the analysis. Exosomes were extracted using solvent exclusion principle and polyethylene glycol induced precipitation. Protein identity and quantification was determined by label-free liquid chromatography mass spectrometry (LC/MS). The mean peak serum creatinine was 3.7 ± 1.4 mg/dL and time to kidney injury was 4.0 ± 3.0 days. At discharge, 90% of patients demonstrated partial recovery; 33% experienced full recovery by day 28. Proteomic analyses on five V-AKI and 7 control samples revealed 2009 proteins in all samples and 251 proteins significantly associated with V-AKI (Pi-score > 1). The top discriminatory proteins were complement C3, complement C4, galectin-3-binding protein, fibrinogen, alpha-2 macroglobulin, immunoglobulin heavy constant mu and serotransferrin. CONCLUSION: Urinary exosomes reveal up-regulation of inflammatory proteins after nephrotoxic injury in V-AKI. Further studies are necessary in a large patient sample to confirm these findings for elucidation of pathophysiologic mechanisms and validation of potential injury biomarkers.

Human Leukocyte Antigen B*14:01 and B*35:01 Are Associated With Trimethoprim-Sulfamethoxazole Induced Liver Injury.

BACKGROUND AND AIMS: Trimethoprim (TMP)-sulfamethoxazole (SMX) is an important cause of idiosyncratic drug-induced liver injury (DILI), but its genetic risk factors are not well understood. This study investigated the relationship between variants in the human leukocyte antigen (HLA) class 1 and 2 genes and well-characterized cases of TMP-SMX DILI. APPROACH AND RESULTS: European American and African American persons with TMP-SMX DILI were compared with respective population controls. HLA sequencing was performed by Illumina MiSeq (Illumina, San Diego, CA) for cases. The HLA genotype imputation with attribute bagging program was used to impute HLA alleles for controls. The allele frequency difference between case patients and controls was tested by Fisher's exact tests for each ethnic group. For European Americans, multivariable logistic regression with Firth penalization was used to test the HLA allelic effect after adjusting for age and the top two principal components. Molecular docking was performed to assess HLA binding with TMP and SMX. The European American subset had 51 case patients and 12,156 controls, whereas the African American subset had 10 case patients and 5,439 controls. Four HLA alleles were significantly associated in the European American subset, with HLA-B*14:01 ranking at the top (odds ratio, 9.20; 95% confidence interval, 3.16, 22.35; P = 0.0003) after covariate adjustment. All carriers of HLA-B*14:01 with TMP-SMX DILI possessed HLA-C*08:02, another significant allele (P = 0.0026). This pattern was supported by HLA-B*14:01-HLA-C*08:02 haplotype association (P = 1.33 × 10-5 ). For the African American patients, HLA-B*35:01 had 2.8-fold higher frequency in case patients than in controls, with 5 of 10 patients carrying this allele. Molecular docking showed cysteine at position 67 in HLA-B*14:01 and phenylalanine at position 67 in HLA-B*35:01 to be the predictive binding sites for SMX metabolites. CONCLUSIONS: HLA-B*14:01 is associated with TMP-SMX DILI in European Americans, and HLA-B*35:01 may be a potential genetic risk factor for African Americans.

Development and Analytical Validation of a 29 Gene Clinical Pharmacogenetic Genotyping Panel: Multi-Ethnic Allele and Copy Number Variant Detection.

To develop a novel pharmacogenetic genotyping panel, a multidisciplinary team evaluated available evidence and selected 29 genes implicated in interindividual drug response variability, including 130 sequence variants and additional copy number variants (CNVs). Of the 29 genes, 11 had guidelines published by the Clinical Pharmacogenetics Implementation Consortium. Targeted genotyping and CNV interrogation were accomplished by multiplex single-base extension using the MassARRAY platform (Agena Biosciences) and multiplex ligation-dependent probe amplification (MRC Holland), respectively. Analytical validation of the panel was accomplished by a strategic combination of > 500 independent tests performed on 170 unique reference material DNA samples, which included sequence variant and CNV accuracy, reproducibility, and specimen (blood, saliva, and buccal swab) controls. Among the accuracy controls were 32 samples from the 1000 Genomes Project that were selected based on their enrichment of sequence variants included in the pharmacogenetic panel (VarCover.org). Coupled with publicly available samples from the Genetic Testing Reference Materials Coordination Program (GeT-RM), accuracy validation material was available for the majority (77%) of interrogated sequence variants (100% with average allele frequencies > 0.1%), as well as additional structural alleles with unique copy number signatures (e.g., CYP2D6*5, *13, *36, *68; CYP2B6*29; and CYP2C19*36). Accuracy and reproducibility for both genotyping and copy number were > 99.9%, indicating that the optimized panel platforms were precise and robust. Importantly, multi-ethnic allele frequencies of the interrogated variants indicate that the vast majority of the general population carries at least one of these clinically relevant pharmacogenetic variants, supporting the implementation of this panel for pharmacogenetic research and/or clinical implementation programs.

Beta-lactam-induced immediate hypersensitivity reactions: A genome-wide association study of a deeply phenotyped cohort.

BACKGROUND: ß-lactam antibiotics are associated with a variety of immune-mediated or hypersensitivity reactions, including immediate (type I) reactions mediated by antigen-specific IgE. OBJECTIVE: We sought to identify genetic predisposing factors for immediate reactions to ß-lactam antibiotics. METHODS: Patients with a clinical history of immediate hypersensitivity reactions to either penicillins or cephalosporins, which were immunologically confirmed, were recruited from allergy clinics. A genome-wide association study was conducted on 662 patients (the discovery cohort) with a diagnosis of immediate hypersensitivity and the main finding was replicated in a cohort of 98 Spanish cases, recruited using the same diagnostic criteria as the discovery cohort. RESULTS: Genome-wide association study identified rs71542416 within the Class II HLA region as the top hit (P = 2 × 10-14); this was in linkage disequilibrium with HLA-DRB1∗10:01 (odds ratio, 2.93; P = 5.4 × 10-7) and HLA-DQA1∗01:05 (odds ratio, 2.93, P = 5.4 × 10-7). Haplotype analysis identified that HLA-DRB1∗10:01 was a risk factor even without the HLA-DQA1∗01:05 allele. The association with HLA-DRB1∗10:01 was replicated in another cohort, with the meta-analysis of the discovery and replication cohorts showing that HLA-DRB1∗10:01 increased the risk of immediate hypersensitivity at a genome-wide level (odds ratio, 2.96; P = 4.1 × 10-9). No association with HLA-DRB1∗10:01 was identified in 268 patients with delayed hypersensitivity reactions to ß-lactams. CONCLUSIONS: HLA-DRB1∗10:01 predisposed to immediate hypersensitivity reactions to penicillins. Further work to identify other predisposing HLA and non-HLA loci is required.

Genetic Risk Factors in Drug-Induced Liver Injury Due to Isoniazid-Containing Antituberculosis Drug Regimens.

Drug-induced liver injury (DILI) is a complication of treatment with antituberculosis (TB) drugs, especially in isoniazid (INH)-containing regimens. To investigate genetic risk factors, we performed a genomewide association study (GWAS) involving anti-TB DILI cases (55 Indian and 70 European) and controls (1,199 Indian and 10,397 European). Most cases were treated with a standard anti-TB drug regimen; all received INH. We imputed single nucleotide polymorphism and HLA genotypes and performed trans-ethnic meta-analysis on GWAS and candidate gene genotypes. GWAS found one significant association (rs117491755) in Europeans only. For HLA, HLA-B*52:01 was significant (meta-analysis odds ratio (OR) 2.67, 95% confidence interval (CI) 1.63-4.37, P = 9.4 × 10-5 ). For N-acetyltransferase 2 (NAT2), NAT2*5 frequency was lower in cases (OR 0.69, 95% CI 0.57-0.83, P = 0.01). NAT2*6 and NAT2*7 were more common, with homozygotes for NAT2*6 and/or NAT2*7 enriched among cases (OR 1.89, 95% CI 0.84-4.22, P = 0.004). We conclude HLA genotype makes a small contribution to TB drug-related DILI and that the NAT2 contribution is complex, but consistent with previous reports when differences in the metabolic effect of NAT2*5 compared with those of NAT2*6 and NAT2*7 are considered.

Polygenic architecture informs potential vulnerability to drug-induced liver injury.

Drug-induced liver injury (DILI) is a leading cause of termination in drug development programs and removal of drugs from the market; this is partially due to the inability to identify patients who are at risk1. In this study, we developed a polygenic risk score (PRS) for DILI by aggregating effects of numerous genome-wide loci identified from previous large-scale genome-wide association studies2. The PRS predicted the susceptibility to DILI in patients treated with fasiglifam, amoxicillin-clavulanate or flucloxacillin and in primary hepatocytes and stem cell-derived organoids from multiple donors treated with over ten different drugs. Pathway analysis highlighted processes previously implicated in DILI, including unfolded protein responses and oxidative stress. In silico screening identified compounds that elicit transcriptomic signatures present in hepatocytes from individuals with elevated PRS, supporting mechanistic links and suggesting a novel screen for safety of new drug candidates. This genetic-, cellular-, organoid- and human-scale evidence underscored the polygenic architecture underlying DILI vulnerability at the level of hepatocytes, thus facilitating future mechanistic studies. Moreover, the proposed 'polygenicity-in-a-dish' strategy might potentially inform designs of safer, more efficient and robust clinical trials.

Erratum to "Uterine Tumor Resembling Ovarian Sex Cord Tumors Type II with Vaginal Vault Recurrence".

[This corrects the article DOI: 10.1155/2019/5231219.].

Identifying the Potential Mechanism of Action of SNPs Associated With Breast Cancer Susceptibility With GVITamIN.

In the last decade, a large number of genome-wide association studies have uncovered many single-nucleotide polymorphisms (SNPs) that are associated with complex traits and confer susceptibility to diseases, such as cancer. However, so far only a few heritable traits with medium-to-high penetrance have been identified. The vast majority of the discovered variants only leads to disease in combination with other still unknown factors. Furthermore, while many studies aimed to link the effect of SNPs to changes in molecular phenotypes, the analysis has been often focused on testing associations between a single SNP and a transcript, hence disregarding the dysregulation of gene regulatory networks that has been shown to play an essential role in disease onset, notably in cancer. Here we take a systems biology approach and develop GVITamIN (Genetic VarIaTIoN functional analysis tool), a new statistical and computational approach to characterize the effect of a SNP on both genes and transcriptional regulatory programs. GVITamIN exploits a novel statistical approach to combine the usually small effect of disease-susceptibility SNPs, and reveals important potential oncogenic mechanisms, hence taking one step further in the direction of understanding the SNP mechanism of action. We apply GVITamIN on a breast cancer cohort and identify well-known cancer-related transcription factors, such as CTCF, LEF1, and FOXA1, as TFs dysregulated by breast cancer-associated SNPs. Furthermore, our results reveal that SNPs located on the RAD51B gene are significantly associated with an abnormal regulatory activity, suggesting a pivotal role for homologous recombination repair mechanisms in breast cancer.