Association of Liver Injury From Specific Drugs, or Groups of Drugs, With Polymorphisms in HLA and Other Genes in a Genome-Wide Association Study.

BACKGROUND & AIMS: We performed a genome-wide association study (GWAS) to identify genetic risk factors for drug-induced liver injury (DILI) from licensed drugs without previously reported genetic risk factors. METHODS: We performed a GWAS of 862 persons with DILI and 10,588 population-matched controls. The first set of cases was recruited before May 2009 in Europe (n = 137) and the United States (n = 274). The second set of cases were identified from May 2009 through May 2013 from international collaborative studies performed in Europe, the United States, and South America. For the GWAS, we included only cases with patients of European ancestry associated with a particular drug (but not flucloxacillin or amoxicillin-clavulanate). We used DNA samples from all subjects to analyze HLA genes and single nucleotide polymorphisms. After the discovery analysis was concluded, we validated our findings using data from 283 European patients with diagnosis of DILI associated with various drugs. RESULTS: We associated DILI with rs114577328 (a proxy for A*33:01 a HLA class I allele; odds ratio [OR], 2.7; 95% confidence interval [CI], 1.9-3.8; P = 2.4 × 10-8) and with rs72631567 on chromosome 2 (OR, 2.0; 95% CI, 1.6-2.5; P = 9.7 × 10-9). The association with A*33:01 was mediated by large effects for terbinafine-, fenofibrate-, and ticlopidine-related DILI. The variant on chromosome 2 was associated with DILI from a variety of drugs. Further phenotypic analysis indicated that the association between DILI and A*33:01 was significant genome wide for cholestatic and mixed DILI, but not for hepatocellular DILI; the polymorphism on chromosome 2 was associated with cholestatic and mixed DILI as well as hepatocellular DILI. We identified an association between rs28521457 (within the lipopolysaccharide-responsive vesicle trafficking, beach and anchor containing gene) and only hepatocellular DILI (OR, 2.1; 95% CI, 1.6-2.7; P = 4.8 × 10-9). We did not associate any specific drug classes with genetic polymorphisms, except for statin-associated DILI, which was associated with rs116561224 on chromosome 18 (OR, 5.4; 95% CI, 3.0-9.5; P = 7.1 × 10-9). We validated the association between A*33:01 terbinafine- and sertraline-induced DILI. We could not validate the association between DILI and rs72631567, rs28521457, or rs116561224. CONCLUSIONS: In a GWAS of persons of European descent with DILI, we associated HLA-A*33:01 with DILI due to terbinafine and possibly fenofibrate and ticlopidine. We identified polymorphisms that appear to be associated with DILI from statins, as well as 2 non-drug-specific risk factors.

An Exome Sequencing Study to Assess the Role of Rare Genetic Variation in Pulmonary Fibrosis.

RATIONALE: Idiopathic pulmonary fibrosis (IPF) is an increasingly recognized, often fatal lung disease of unknown etiology. OBJECTIVES: The aim of this study was to use whole-exome sequencing to improve understanding of the genetic architecture of pulmonary fibrosis. METHODS: We performed a case-control exome-wide collapsing analysis including 262 unrelated individuals with pulmonary fibrosis clinically classified as IPF according to American Thoracic Society/European Respiratory Society/Japanese Respiratory Society/Latin American Thoracic Association guidelines (81.3%), usual interstitial pneumonia secondary to autoimmune conditions (11.5%), or fibrosing nonspecific interstitial pneumonia (7.2%). The majority (87%) of case subjects reported no family history of pulmonary fibrosis. MEASUREMENTS AND MAIN RESULTS: We searched 18,668 protein-coding genes for an excess of rare deleterious genetic variation using whole-exome sequence data from 262 case subjects with pulmonary fibrosis and 4,141 control subjects drawn from among a set of individuals of European ancestry. Comparing genetic variation across 18,668 protein-coding genes, we found a study-wide significant (P < 4.5 × 10-7) case enrichment of qualifying variants in TERT, RTEL1, and PARN. A model qualifying ultrarare, deleterious, nonsynonymous variants implicated TERT and RTEL1, and a model specifically qualifying loss-of-function variants implicated RTEL1 and PARN. A subanalysis of 186 case subjects with sporadic IPF confirmed TERT, RTEL1, and PARN as study-wide significant contributors to sporadic IPF. Collectively, 11.3% of case subjects with sporadic IPF carried a qualifying variant in one of these three genes compared with the 0.3% carrier rate observed among control subjects (odds ratio, 47.7; 95% confidence interval, 21.5-111.6; P = 5.5 × 10-22). CONCLUSIONS: We identified TERT, RTEL1, and PARN-three telomere-related genes previously implicated in familial pulmonary fibrosis-as significant contributors to sporadic IPF. These results support the idea that telomere dysfunction is involved in IPF pathogenesis.

Features of Autoimmune Hepatitis in Patients With Drug-induced Liver Injury.

BACKGROUND & AIMS: Drug-induced liver injury (DILI) has features similar to those of other liver diseases including autoimmune hepatitis (AIH). We aimed to characterize the clinical and autoimmune features of liver injury caused by nitrofurantoin, minocycline, methyldopa, or hydralazine. METHODS: We analyzed data from 88 cases of DILI attributed to nitrofurantoin, minocycline, methyldopa, or hydralazine included in the Drug-Induced Liver Injury Network prospective study from 2004 through 2014. Sera were collected from patients at baseline and follow-up examination and tested for levels of immunoglobulin G (IgG), antibodies to nuclear antigen (ANA), smooth muscle (SMA), and soluble liver antigen (SLA). An autoimmune score was derived on the basis of increases in levels of IgG, ANA, SMA, and SLA (assigned values of 0, 1+, or 2+). AIH-associated HLA-DRB1*03:01 and HLA-DRB1*04:01 allele frequencies were compared with those of the general population (controls). RESULTS: Of the 88 cases, 80 were women (91%), 74% had hepatocellular injury, and 25% had severe injury. At the onset of DILI, 39% of cases had increased levels of IgG, 72% had increased levels of ANA, 60% had increased levels of SMA, and none had increases in SLA. A phenotype of autoimmunity (autoimmune score ≥2) was observed in 82% of cases attributed to nitrofurantoin and 73% of cases attributed to minocycline (73%) but only 55% of cases attributed to methyldopa and 43% of cases attributed to hydralazine (P = .16 for nitrofurantoin and minocycline vs methyldopa and hydralazine). We observed a decrease in numbers of serum samples positive for ANA (P = .01) or SMA (P < .001) and in autoimmune scores (P < .001) between DILI onset and follow-up. Similar percentages of patients with DILI had HLA-DRB1*03:01 (15%) and HLA-DRB1*04:01 (9%) as controls (12% and 9%, respectively). CONCLUSIONS: In analysis of data from the DILIN prospective study, we found that most cases of DILI attributed to nitrofurantoin or minocycline and about half of cases that were due to methyldopa and hydralazine have a phenotype of autoimmunity similar to AIH. These features decrease with recovery of the injury and are not associated with the typical HLA alleles found in patients with idiopathic AIH.

An exploratory genome-wide analysis of genetic risk for alcoholic hepatitis.

OBJECTIVES: To elucidate the genetic variability between heavy drinkers with and without alcoholic hepatitis (AH). MATERIALS AND METHODS: An exploratory genome-wide association study (GWAS; NCT02172898) was conducted comparing 90 AH cases with 93 heavy drinking matched controls without liver disease in order to identify variants or genes associated with risk for AH. Individuals were genotyped using the multi-ethnic genotyping array, after which the data underwent conventional quality control. Using bioinformatics tools, pathways associated with AH were explored on the basis of individual variants, and based on genes with a higher 'burden' of functional variation. RESULTS: Although no single variant reached genome-wide significance, an association signal was observed for PNPLA3 rs738409 (p = .01, OR 1.9, 95% CI 1.1-3.1), a common single nucleotide polymorphism that has been associated with a variety of liver-related pathologies including alcoholic cirrhosis. Using the improved gene set enrichment analysis for GWAS tool, it was shown that, based on the single variants' trait-association p-values, multiple pathways were associated with risk for AH with high confidence (false discovery rate [FDR] < 0.05), including several pathways involved in lymphocyte activation and chemokine signaling, which coincides with findings from other research groups. Several Tox Functions and Canonical Pathways were highlighted using Ingenuity Pathway Analysis, with an especially conspicuous role for pathways related to ethanol degradation, which is not surprising considering the phenotype of the genotyped individuals. CONCLUSION: This preliminary analysis suggests a role for PNPLA3 variation and several gene sets/pathways that may influence risk for AH among heavy drinkers.

HLA-DRB1*16: 01-DQB1*05: 02 is a novel genetic risk factor for flupirtine-induced liver injury.

OBJECTIVE: Flupirtine is a nonopioid analgesic with regulatory approval in a number of European countries. Because of the risk of serious liver injury, its use is now limited to short-term pain management. We aimed to identify genetic risk factors for flupirtine-related drug-induced liver injury (DILI) as these are unknown. MATERIALS AND METHODS: Six flupirtine-related DILI patients from Germany were included in a genome-wide association study (GWAS) involving a further 614 European cases of DILI because of other drugs and 10,588 population controls. DILI was diagnosed by causality assessment and expert review. Human leucocyte antigen (HLA) and single nucleotide polymorphism genotypes were imputed from the GWAS data, with direct HLA typing performed on selected cases to validate HLA predictions. Four replication cases that were unavailable for the GWAS were genotyped by direct HLA typing, yielding an overall total of 10 flupirtine DILI cases. RESULTS: In the six flupirtine DILI cases included in the GWAS, we found a significant enrichment of the DRB1*16:01-DQB1*05:02 haplotype compared with the controls (minor allele frequency cases 0.25 and minor allele frequency controls 0.013; P=1.4 × 10(-5)). We estimated an odds ratio for haplotype carriers of 18.7 (95% confidence interval 2.5-140.5, P=0.002) using population-specific HLA control data. The result was replicated in four additional cases, also with a haplotype frequency of 0.25. In the combined cohort (six GWAS plus four replication cases), the haplotype was also significant (odds ratio 18.7, 95% confidence interval 4.31-81.42, P=6.7 × 10(-5)). CONCLUSION: We identified a novel HLA class II association for DILI, confirming the important contribution of HLA genotype towards the risk of DILI generally.

A nonsense mutation in IKBKB causes combined immunodeficiency.

Identification of the molecular etiologies of primary immunodeficiencies has led to important insights into the development and function of the immune system. We report here the cause of combined immunodeficiency in 4 patients from 2 different consanguineous Qatari families with similar clinical and immunologic phenotypes. The patients presented at an early age with fungal, viral, and bacterial infections and hypogammaglobulinemia. Although their B- and T-cell numbers were normal, they had low regulatory T-cell and NK-cell numbers. Moreover, patients' T cells were mostly CD45RA(+)-naive cells and were defective in activation after T-cell receptor stimulation. All patients contained the same homozygous nonsense mutation in IKBKB (R286X), revealed by whole-exome sequencing with undetectable IKKß and severely decreased NEMO proteins. Mutant IKKß(R286X) was unable to complex with IKKα/NEMO. Immortalized patient B cells displayed impaired IκBα phosphorylation and NFκB nuclear translocation. These data indicate that mutated IKBKB is the likely cause of immunodeficiency in these 4 patients.

ITPA gene variants protect against anaemia in patients treated for chronic hepatitis C.

Chronic infection with the hepatitis C virus (HCV) affects 170 million people worldwide and is an important cause of liver-related morbidity and mortality. The standard of care therapy combines pegylated interferon (pegIFN) alpha and ribavirin (RBV), and is associated with a range of treatment-limiting adverse effects. One of the most important of these is RBV-induced haemolytic anaemia, which affects most patients and is severe enough to require dose modification in up to 15% of patients. Here we show that genetic variants leading to inosine triphosphatase deficiency, a condition not thought to be clinically important, protect against haemolytic anaemia in hepatitis-C-infected patients receiving RBV.

Population pharmacokinetic modeling of plasma and intracellular ribavirin concentrations in patients with chronic hepatitis C virus infection.

Ribavirin, a guanosine analog, is a broad-spectrum antiviral agent. Ribavirin has been a fundamental component of the treatment of hepatitis C virus (HCV) infection for decades, but there is a very limited understanding of the clinical pharmacology of this drug. Furthermore, it is associated with a major dose-limiting toxicity, hemolytic anemia. Ribavirin undergoes intracellular phosphorylation by host enzymes to ribavirin monophosphate (RMP), ribavirin diphosphate (RDP), and ribavirin triphosphate (RTP). The intracellular forms have been associated with antiviral and toxic effects in vitro, but the kinetics of these phosphorylated moieties have not been fully elucidated in vivo. We developed a model to characterize the plasma pharmacokinetics of ribavirin and the difference between intracellular phosphorylation kinetics in red cells (nonnucleated) and in peripheral blood mononuclear cells (nucleated). A time-independent two-compartment model with first-order absorption described the plasma data well. The cellular phosphorylation kinetics was described by a one-compartment model for RMP, with the formation rate driven by plasma concentrations and the first-order degradation rate. RDP and RTP rapidly reached equilibrium with RMP. Concomitant telaprevir use, inosine triphosphatase genetics, creatinine clearance, weight, and sex were significant covariates. The terminal ribavirin half-life in plasma and phosphorylated anabolites in cells was approximately 224 h. We found no evidence of time-dependent kinetics. These data provide a foundation for uncovering concentration-effect associations for ribavirin and determining the optimal dose and duration of this drug for use in combination with newer direct-acting HCV agents. (This study has been registered at ClinicalTrials.gov under registration no. NCT01097395.).

Prioritizing genetic variants for causality on the basis of preferential linkage disequilibrium.

To date, the widely used genome-wide association studies (GWASs) of the human genome have reported thousands of variants that are significantly associated with various human traits. However, in the vast majority of these cases, the causal variants responsible for the observed associations remain unknown. In order to facilitate the identification of causal variants, we designed a simple computational method called the "preferential linkage disequilibrium (LD)" approach, which follows the variants discovered by GWASs to pinpoint the causal variants, even if they are rare compared with the discovery variants. The approach is based on the hypothesis that the GWAS-discovered variant is better at tagging the causal variants than are most other variants evaluated in the original GWAS. Applying the preferential LD approach to the GWAS signals of five human traits for which the causal variants are already known, we successfully placed the known causal variants among the top ten candidates in the majority of these cases. Application of this method to additional GWASs, including those of hepatitis C virus treatment response, plasma levels of clotting factors, and late-onset Alzheimer disease, has led to the identification of a number of promising candidate causal variants. This method represents a useful tool for delineating causal variants by bringing together GWAS signals and the rapidly accumulating variant data from next-generation sequencing.

Prevalence of genetic variants of keratins 8 and 18 in patients with drug-induced liver injury.

BACKGROUND: Keratin 8 and 18 (K8/K18) cytoskeletal proteins protect hepatocytes from undergoing apoptosis and their mutations predispose to adverse outcomes in acute liver failure (ALF). All known K8/K18 variants occur at relatively non-conserved residues and do not cause keratin cytoskeleton reorganization, whereas epidermal keratin-conserved residue mutations disrupt the keratin cytoskeleton and cause severe skin disease. The aim of our study was to identify keratin variants in idiosyncratic drug-induced liver injury (DILI). METHODS: Genomic DNA was isolated from 800 patients enrolled in an ongoing US multicenter study, with DILI attributed to a wide range of drugs. Specific K8/K18 exonic regions were PCR-amplified and screened by denaturing HPLC followed by DNA sequencing. The functional impact of keratin variants was assessed using cell transfection and immune staining. RESULTS: Heterozygous and compound amino acid-altering K8/K18 variants were identified in 86 DILI patients and non-coding variants in 15 subjects. Five novel amino acid-altering (K8 Lys393Arg, K8 Ala351Val, K8 Ala358Val, K8 Ile346Val, K18 Asp89His) and two non-coding variants were observed. Several variants segregated with specific ethnic backgrounds but were found at similar frequencies in DILI subjects and ethnically matched population controls. Notably, variants in highly conserved residues of K8 Lys393Arg (ezetimibe/simvastatin-related) and K18 Asp89His (isoniazid-related) were found in patients with fatal DILI. These novel variants also led to keratin network disruption in transfected cells. CONCLUSIONS: Novel K8/K18 cytoskeleton-disrupting variants were identified in two patients and segregated with fatal DILI. Other non-cytoskeleton-disrupting keratin variants did not preferentially associate with DILI.

Kinetic differences in the induction of interferon stimulated genes by interferon-α and interleukin 28B are altered by infection with hepatitis C virus.

UNLABELLED: Several genome-wide association studies (GWAS) have identified a genetic polymorphism associated with the gene locus for interleukin 28B (IL28B), a type III interferon (IFN), as a major predictor of clinical outcome in hepatitis C. Antiviral effects of the type III IFN family have previously been shown against several viruses, including hepatitis C virus (HCV), and resemble the function of type I IFN including utilization of the intracellular Janus kinase signal transducer and activator of transcription (JAK-STAT) pathway. Effects unique to IL28B that would distinguish it from IFN-α are not well defined. By analyzing the transcriptomes of primary human hepatocytes (PHH) treated with IFN-α or IL28B, we sought to identify functional differences between IFN-α and IL28B to better understand the roles of these cytokines in the innate immune response. Although our data did not reveal distinct gene signatures, we detected striking kinetic differences between IFN-α and IL28B stimulation for interferon stimulated genes (ISGs). While gene induction was rapid and peaked at 8 hours of stimulation with IFN-α in PHH, IL28B produced a slower, but more sustained increase in gene expression. We confirmed these findings in the human hepatoma cell line Huh7.5.1. Interestingly, in HCV-infected cells the rapid response after stimulation with IFN-α was blunted, and the induction pattern resembled that caused by IL28B. CONCLUSION: The kinetics of gene induction are fundamentally different for stimulations with either IFN-α or IL28B in hepatocytes, suggesting distinct roles of these cytokines within the immune response. Furthermore, the observed differences are substantially altered by infection with HCV.

Toward predicting drug-induced liver injury: parallel computational approaches to identify multidrug resistance protein 4 and bile salt export pump inhibitors.

Drug-induced liver injury (DILI) is an important cause of drug toxicity. Inhibition of multidrug resistance protein 4 (MRP4), in addition to bile salt export pump (BSEP), might be a risk factor for the development of cholestatic DILI. Recently, we demonstrated that inhibition of MRP4, in addition to BSEP, may be a risk factor for the development of cholestatic DILI. Here, we aimed to develop computational models to delineate molecular features underlying MRP4 and BSEP inhibition. Models were developed using 257 BSEP and 86 MRP4 inhibitors and noninhibitors in the training set. Models were externally validated and used to predict the affinity of compounds toward BSEP and MRP4 in the DrugBank database. Compounds with a score above the median fingerprint threshold were considered to have significant inhibitory effects on MRP4 and BSEP. Common feature pharmacophore models were developed for MRP4 and BSEP with LigandScout software using a training set of nine well characterized MRP4 inhibitors and nine potent BSEP inhibitors. Bayesian models for BSEP and MRP4 inhibition/noninhibition were developed with cross-validated receiver operator curve values greater than 0.8 for the test sets, indicating robust models with acceptable false positive and false negative prediction rates. Both MRP4 and BSEP inhibitor pharmacophore models were characterized by hydrophobic and hydrogen-bond acceptor features, albeit in distinct spatial arrangements. Similar molecular features between MRP4 and BSEP inhibitors may partially explain why various drugs have affinity for both transporters. The Bayesian (BSEP, MRP4) and pharmacophore (MRP4, BSEP) models demonstrated significant classification accuracy and predictability.

Genetic variation in IL28B predicts hepatitis C treatment-induced viral clearance.

Chronic infection with hepatitis C virus (HCV) affects 170 million people worldwide and is the leading cause of cirrhosis in North America. Although the recommended treatment for chronic infection involves a 48-week course of peginterferon-alpha-2b (PegIFN-alpha-2b) or -alpha-2a (PegIFN-alpha-2a) combined with ribavirin (RBV), it is well known that many patients will not be cured by treatment, and that patients of European ancestry have a significantly higher probability of being cured than patients of African ancestry. In addition to limited efficacy, treatment is often poorly tolerated because of side effects that prevent some patients from completing therapy. For these reasons, identification of the determinants of response to treatment is a high priority. Here we report that a genetic polymorphism near the IL28B gene, encoding interferon-lambda-3 (IFN-lambda-3), is associated with an approximately twofold change in response to treatment, both among patients of European ancestry (P = 1.06 x 10(-25)) and African-Americans (P = 2.06 x 10(-3)). Because the genotype leading to better response is in substantially greater frequency in European than African populations, this genetic polymorphism also explains approximately half of the difference in response rates between African-Americans and patients of European ancestry.

Genetic basis of drug-induced liver injury: present and future.

There is considerable evidence that susceptibility to idiosyncratic drug-induced liver injury (DILI) is genetically determined. Though genetic associations with DILI have been reported since the 1980s, the development of genome-wide association studies has enabled genetic risk factors for DILI, in common with other diseases, to be detected and confirmed more confidently. Human leukocyte antigen (HLA) genotype has been demonstrated to be a strong risk factor for development of DILI with a range of drugs and the underlying mechanism, probably involving presentation of a drug-peptide complex to T cells is increasingly well understood. However, specific HLA alleles are not associated with all forms of DILI and non-HLA genetic risk factors, especially those relating to drug disposition, also appear to contribute. For some drugs, there is evidence of a dual role for HLA and drug metabolism genes. Though the associations with non-HLA genes have been less well replicated than the HLA associations, there is increasing evidence that drug metabolism genes such as NAT2 and UGT2B7 contribute to some forms of DILI. Translating current genetic findings on DILI susceptibility to the clinic has been relatively slow, but some progress is now being made. In the future, DNA sequencing may lead to the identification of rare variants that contribute to DILI. Developments in the related area of epigenomics and in the development of improved models for DILI by use of genetically defined induced pluripotent stem cells should improve understanding of the biology of DILI and inform drug development.

Risk factors for development of cholestatic drug-induced liver injury: inhibition of hepatic basolateral bile acid transporters multidrug resistance-associated proteins 3 and 4.

Impaired hepatic bile acid export may contribute to development of cholestatic drug-induced liver injury (DILI). The multidrug resistance-associated proteins (MRP) 3 and 4 are postulated to be compensatory hepatic basolateral bile acid efflux transporters when biliary excretion by the bile salt export pump (BSEP) is impaired. BSEP inhibition is a risk factor for cholestatic DILI. This study aimed to characterize the relationship between MRP3, MRP4, and BSEP inhibition and cholestatic potential of drugs. The inhibitory effect of 88 drugs (100 µM) on MRP3- and MRP4-mediated substrate transport was measured in membrane vesicles. Drugs selected for investigation included 50 BSEP non-inhibitors (24 non-cholestatic; 26 cholestatic) and 38 BSEP inhibitors (16 non-cholestatic; 22 cholestatic). MRP4 inhibition was associated with an increased risk of cholestatic potential among BSEP non-inhibitors. In this group, for each 1% increase in MRP4 inhibition, the odds of the drug being cholestatic increased by 3.1%. Using an inhibition cutoff of 21%, which predicted a 50% chance of cholestasis, 62% of cholestatic drugs inhibited MRP4 (P < 0.05); in contrast, only 17% of non-cholestatic drugs were MRP4 inhibitors. Among BSEP inhibitors, MRP4 inhibition did not provide additional predictive value of cholestatic potential; almost all BSEP inhibitors were also MRP4 inhibitors. Inclusion of pharmacokinetic predictor variables (e.g., maximal unbound concentration in plasma) in addition to percent MRP4 inhibition in logistic regression models did not improve cholestasis prediction. Association of cholestasis with percent MRP3 inhibition was not statistically significant, regardless of BSEP-inhibition status. Inhibition of MRP4, in addition to BSEP, may be a risk factor for the development of cholestatic DILI.

Copy number variation of KIR genes influences HIV-1 control.

A genome-wide screen for large structural variants showed that a copy number variant (CNV) in the region encoding killer cell immunoglobulin-like receptors (KIR) associates with HIV-1 control as measured by plasma viral load at set point in individuals of European ancestry. This CNV encompasses the KIR3DL1-KIR3DS1 locus, encoding receptors that interact with specific HLA-Bw4 molecules to regulate the activation of lymphocyte subsets including natural killer (NK) cells. We quantified the number of copies of KIR3DS1 and KIR3DL1 in a large HIV-1 positive cohort, and showed that an increase in KIR3DS1 count associates with a lower viral set point if its putative ligand is present (p = 0.00028), as does an increase in KIR3DL1 count in the presence of KIR3DS1 and appropriate ligands for both receptors (p = 0.0015). We further provide functional data that demonstrate that NK cells from individuals with multiple copies of KIR3DL1, in the presence of KIR3DS1 and the appropriate ligands, inhibit HIV-1 replication more robustly, and associated with a significant expansion in the frequency of KIR3DS1+, but not KIR3DL1+, NK cells in their peripheral blood. Our results suggest that the relative amounts of these activating and inhibitory KIR play a role in regulating the peripheral expansion of highly antiviral KIR3DS1+ NK cells, which may determine differences in HIV-1 control following infection.

A program for real-time surveillance of SARS-CoV-2 genetics.

The COVID-19 pandemic brought forth an urgent need for widespread genomic surveillance for rapid detection and monitoring of emerging SARS-CoV-2 variants. It necessitated design, development, and deployment of a nationwide infrastructure designed for sequestration, consolidation, and characterization of patient samples that disseminates de-identified information to public authorities in tight turnaround times. Here, we describe our development of such an infrastructure, which sequenced 594,832 high coverage SARS-CoV-2 genomes from isolates we collected in the United States (U.S.) from March 13th 2020 to July 3rd 2023. Our sequencing protocol ('Virseq') utilizes wet and dry lab procedures to generate mutation-resistant sequencing of the entire SARS-CoV-2 genome, capturing all major lineages. We also characterize 379 clinically relevant SARS-CoV-2 multi-strain co-infections and ensure robust detection of emerging lineages via simulation. The modular infrastructure, sequencing, and analysis capabilities we describe support the U.S. Centers for Disease Control and Prevention national surveillance program and serve as a model for rapid response to emerging pandemics at a national scale.

Rare deletions at 16p13.11 predispose to a diverse spectrum of sporadic epilepsy syndromes.

Deletions at 16p13.11 are associated with schizophrenia, mental retardation, and most recently idiopathic generalized epilepsy. To evaluate the role of 16p13.11 deletions, as well as other structural variation, in epilepsy disorders, we used genome-wide screens to identify copy number variation in 3812 patients with a diverse spectrum of epilepsy syndromes and in 1299 neurologically-normal controls. Large deletions (> 100 kb) at 16p13.11 were observed in 23 patients, whereas no control had a deletion greater than 16 kb. Patients, even those with identically sized 16p13.11 deletions, presented with highly variable epilepsy phenotypes. For a subset of patients with a 16p13.11 deletion, we show a consistent reduction of expression for included genes, suggesting that haploinsufficiency might contribute to pathogenicity. We also investigated another possible mechanism of pathogenicity by using hybridization-based capture and next-generation sequencing of the homologous chromosome for ten 16p13.11-deletion patients to look for unmasked recessive mutations. Follow-up genotyping of suggestive polymorphisms failed to identify any convincing recessive-acting mutations in the homologous interval corresponding to the deletion. The observation that two of the 16p13.11 deletions were larger than 2 Mb in size led us to screen for other large deletions. We found 12 additional genomic regions harboring deletions > 2 Mb in epilepsy patients, and none in controls. Additional evaluation is needed to characterize the role of these exceedingly large, non-locus-specific deletions in epilepsy. Collectively, these data implicate 16p13.11 and possibly other large deletions as risk factors for a wide range of epilepsy disorders, and they appear to point toward haploinsufficiency as a contributor to the pathogenicity of deletions.

Limited contribution of common genetic variants to risk for liver injury due to a variety of drugs.

BACKGROUND AND AIMS: Drug-induced liver injury (DILI) is a serious adverse drug event that is suspected to have a heritable component. We carried out a genome-wide association study of 783 individuals of European ancestry who experienced DILI due to more than 200 implicated drugs. METHODS: DILI patients from the US-based Drug-Induced Liver Injury Network (n=401) and three international registries (n=382) were genotyped with the Illumina 1Mduo BeadChip and compared with population controls (n=3001). Potential associations were tested in 307 independent Drug-Induced Liver Injury Network cases. RESULTS: After accounting for known major histocompatibility complex risk alleles for flucloxacillin-DILI and amoxicillin/clavulanate-DILI, there were no genome-wide significant associations, including in the major histocompatibility complex region. Stratification of DILI cases according to clinical phenotypes (injury type, latency, age of onset) also did not show significant associations. An analysis of hepatocellular DILI (n=285) restricted to 193 single-nucleotide polymorphisms previously associated with autoimmune disease showed a trend association for rs7574865, in the vicinity of signal transducer and activator of transcription 4 (STAT4) (P=4.5×10(-4)). This association was replicated in an independent cohort of 168 hepatocellular DILI cases (P=0.011 and 1.5×10(-5) for combined cohorts). No significant associations were found with stratification by other clinical or demographic variables. CONCLUSION: Although not significant at the genome-wide level, the association between hepatocellular DILI and STAT4 is consistent with the emerging role of the immune system in DILI. However, the lack of genome-wide association study findings supports the idea that strong genetic determinants of DILI may be largely drug-specific or may reflect rare genetic variations, which were not assessed in our study.

Genome-wide association study of interferon-related cytopenia in chronic hepatitis C patients.

BACKGROUND & AIMS: Interferon-alfa (IFN)-related cytopenias are common and may be dose-limiting. We performed a genome wide association study on a well-characterized genotype 1 HCV cohort to identify genetic determinants of peginterferon-α (pegIFN)-related thrombocytopenia, neutropenia, and leukopenia. METHODS: 1604/3070 patients in the IDEAL study consented to genetic testing. Trial inclusion criteria included a platelet (Pl) count ≥80×10(9)/L and an absolute neutrophil count (ANC) ≥1500/mm(3). Samples were genotyped using the Illumina Human610-quad BeadChip. The primary analyses focused on the genetic determinants of quantitative change in cell counts (Pl, ANC, lymphocytes, monocytes, eosinophils, and basophils) at week 4 in patients >80% adherent to therapy (n=1294). RESULTS: 6 SNPs on chromosome 20 were positively associated with Pl reduction (top SNP rs965469, p=10(-10)). These tag SNPs are in high linkage disequilibrium with 2 functional variants in the ITPA gene, rs1127354 and rs7270101, that cause ITPase deficiency and protect against ribavirin (RBV)-induced hemolytic anemia (HA). rs1127354 and rs7270101 showed strong independent associations with Pl reduction (p=10(-12), p=10(-7)) and entirely explained the genome-wide significant associations. We believe this is an example of an indirect genetic association due to a reactive thrombocytosis to RBV-induced anemia: Hb decline was inversely correlated with Pl reduction (r=-0.28, p=10(-17)) and Hb change largely attenuated the association between the ITPA variants and Pl reduction in regression models. No common genetic variants were associated with pegIFN-induced neutropenia or leucopenia. CONCLUSIONS: Two ITPA variants were associated with thrombocytopenia; this was largely explained by a thrombocytotic response to RBV-induced HA attenuating IFN-related thrombocytopenia. No genetic determinants of pegIFN-induced neutropenia were identified.