Correction: High-resolution DNA size enrichment using a magnetic nano-platform and application in non-invasive prenatal testing.

Correction for 'High-resolution DNA size enrichment using a magnetic nano-platform and application in non-invasive prenatal testing' by Bo Zhang et al., Analyst, 2020, 145, 5733-5739, DOI: 10.1039/D0AN00813C.

High-resolution DNA size enrichment using a magnetic nano-platform and application in non-invasive prenatal testing.

Precise DNA sizing can boost sequencing efficiency, reduce cost, improve data quality, and even allow sequencing of low-input samples, while current pervasive DNA sizing approaches are incapable of differentiating DNA fragments under 200 bp with high resolution (<20 bp). In non-invasive prenatal testing (NIPT), the size distribution of cell-free fetal DNA in maternal plasma (main peak at 143 bp) is significantly different from that of maternal cell-free DNA (main peak at 166 bp). The current pervasive workflow of NIPT and DNA sizing is unable to take advantage of this 20 bp difference, resulting in sample rejection, test inaccuracy, and restricted clinical utility. Here we report a simple, automatable, high-resolution DNA size enrichment workflow, named MiniEnrich, on a magnetic nano-platform to exploit this 20 bp size difference and to enrich fetal DNA fragments from maternal blood. Two types of magnetic nanoparticles were developed, with one able to filter high-molecular-weight DNA with high resolution and the other able to recover the remaining DNA fragments under the size threshold of interest with >95% yield. Using this method, the average fetal fraction was increased from 13% to 20% after the enrichment, as measured by plasma DNA sequencing. This approach provides a new tool for high-resolution DNA size enrichment under 200 bp, which may improve NIPT accuracy by rescuing rejected non-reportable clinical samples, and enable NIPT earlier in pregnancy. It also has the potential to improve non-invasive screening for fetal monogenic disorders, differentiate tumor-related DNA in liquid biopsy and find more applications in autoimmune disease diagnosis.

G3viz: an R package to interactively visualize genetic mutation data using a lollipop-diagram.

SUMMARY: The lollipop-diagram is one of the widely used graphical representations to visualize and explore translational effects of genetic mutations in cancer genomics. However, an easy-to-use lollipop-diagram tool with full functionality is still lacking. Here, we introduce g3viz, an R package that enables researchers to explore genetic mutation data using a lollipop-diagram in a web browser. With a few lines of R code, users can interactively visualize data details, annotate findings and export resultant diagrams in high-quality figures. Because of usefulness and usability, g3viz can be generally exploited by researchers with different levels of bioinformatics skills and programming experience. AVAILABILITY AND IMPLEMENTATION: The R package is freely available under the MIT license from CRAN (http://cran.r-project.org/web/packages/g3viz). The g3lollipop JavaScript package is freely available under MIT license at GitHub (https://github.com/g3viz/g3lollipop.js). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Differential Serum Cytokine Profiles in Patients with Chronic Hepatitis B, C, and Hepatocellular Carcinoma.

Cytokines play an important role in the pathogenesis of cirrhosis and hepatocellular carcinoma (HCC), most cases of which are related to either hepatitis B virus (HBV) or hepatitis C virus (HCV). Prior studies have examined differences in individual cytokine levels in patients with chronic liver disease, but comprehensive cytokine profiling data across different clinical characteristics are lacking. We examined serum cytokine profiles of 411 patients with HCC (n = 102: 32% HBV, 54% HCV, 14% non-viral) and without HCC (n = 309: 39% HBV, 39% HCV, 22% non-viral). Multiplex analysis (Luminex 200 IS) was used to measure serum levels of 51 common cytokines. Random forest machine learning was used to obtain receiver operator characteristic curves and to determine individual cytokine importance using Z scores of mean fluorescence intensity for individual cytokines. Among HCC and non-HCC patients, cytokine profiles differed between HBV and HCV patients (area under curve (AUC) 0.82 for HCC, 0.90 for non-HCC). Cytokine profiles did not distinguish cirrhotic HBV patients with and without HCC (AUC 0.503) or HCV patients with and without HCC (AUC 0.63). In conclusion, patients with HBV or HCV infection, with or without HCC, have distinctly different cytokine profiles, suggesting potential differences in disease pathogenesis and/or disease characteristics.

Soluble intercellular adhesion molecule-1 is associated with hepatocellular carcinoma risk: multiplex analysis of serum markers.

Individualized assessment of hepatocellular carcinoma (HCC) risk in chronic liver disease remains challenging. Serum biomarkers including cytokines may offer helpful adjuncts to standard parameters for risk prediction. Our aim was to identify markers associated with increased HCC incidence. This was a prospective cohort study of 282 patients with both viral or non-viral chronic liver disease. Baseline serum cytokines and other markers were measured in multiplex with a commercially-available Luminex-based system. Patients were followed until death or HCC diagnosis. We performed Lasso-based survival analysis to determine parameters associated with HCC development. Cytokine mean florescence intensity (MFI) was the primary predictor and HCC development the primary outcome. 25 patients developed HCC with total follow-up of 1,363 person-years. Parameters associated with increased HCC incidence were cirrhosis, hepatic decompensation, and soluble serum intercellular adhesion molecule 1 (sICAM-1) MFI. No other molecules increased predictive power for HCC incidence. On univariate analysis, the parameters associated with HCC incidence in patients with cirrhosis were age, antiviral treatment, and high sICAM-1 MFI; on multivariate analysis, sICAM-1 remained associated with HCC development (adjusted HR = 2.75). On unbiased screening of serum cytokines and other markers in a diverse cohort, baseline sICAM-1 MFI is associated with HCC incidence.

Big Data Analytics for Genomic Medicine.

Genomic medicine attempts to build individualized strategies for diagnostic or therapeutic decision-making by utilizing patients' genomic information. Big Data analytics uncovers hidden patterns, unknown correlations, and other insights through examining large-scale various data sets. While integration and manipulation of diverse genomic data and comprehensive electronic health records (EHRs) on a Big Data infrastructure exhibit challenges, they also provide a feasible opportunity to develop an efficient and effective approach to identify clinically actionable genetic variants for individualized diagnosis and therapy. In this paper, we review the challenges of manipulating large-scale next-generation sequencing (NGS) data and diverse clinical data derived from the EHRs for genomic medicine. We introduce possible solutions for different challenges in manipulating, managing, and analyzing genomic and clinical data to implement genomic medicine. Additionally, we also present a practical Big Data toolset for identifying clinically actionable genetic variants using high-throughput NGS data and EHRs.

Baseline interpatient hepatitis B viral diversity differentiates HBsAg outcomes in patients treated with tenofovir disoproxil fumarate.

BACKGROUND AND AIMS: HBsAg loss is a desired, but rare, treatment-induced clinical endpoint in chronic hepatitis B (CHB). Few studies have evaluated viral factors contributing to HBsAg loss. METHODS: This study evaluated baseline interpatient sequence diversity across the HBV genome in tenofovir disoproxil fumarate-treated patients who lost HBsAg and compared it to that of control patients with high HBsAg levels throughout therapy. Twenty-one HBeAg+ patients (14 genotype (GT) A and 7 GT D) who achieved HBsAg loss and 27 controls (17 GT A and 10 GT D), were analyzed. Population sequencing was performed on baseline samples and pairwise genetic distances were calculated for 17 overlapping regions across the HBV genome as a measure of interpatient viral diversity. RESULTS: Overall, viral diversity was up to 10-fold higher across GT D patients compared to GT A patients throughout the HBV genome. Within the pol/RT and HBs genes, interpatient viral diversity was significantly lower among HBsAg loss patients for both GT A and D, with the difference driven largely by a reduction in diversity in the small S gene. Conversely, interpatient viral diversity was generally higher in HBsAg loss patients across the HBx gene regulatory elements and precore region. CONCLUSION: In HBsAg loss patients, less interpatient viral diversity was observed within structural-coding regions while specific regions across the HBx and precore genes encoding nonstructural regulatory elements generally displayed higher interpatient viral diversity. These distinct patterns may reflect different responses to adaptive pressure for HBV genomic structural and nonstructural elements.

Genome-wide association study to characterize serum bilirubin elevations in patients with HCV treated with GS-9256, an HCV NS3 serine protease inhibitor.

BACKGROUND: Protease inhibitors for the treatment of HCV can cause mild and reversible elevations of unconjugated bilirubin. We sought to characterize genetic determinants of bilirubin elevations using a genome-wide approach among patients with genotype 1 HCV who received combination therapy that included GS-9256, a novel potent inhibitor of HCV NS3 serine protease, as part of a Phase IIb trial. METHODS: Of the 200 patients sampled, 176 had confirmed European ancestry and were included in the analysis. Infinium HumanOmni5BeadChip (Illumina, Inc., San Diego, CA, USA) was used for genotyping. A categorical analysis of low (grade 0-1) versus high (grade 2-4) bilirubin toxicity grade and a quantitative trait locus mapping of peak bilirubin concentrations was performed. RESULTS: A total of 4,466,809 genetic markers were analysed. No single variant showed a statistically significant association with observed bilirubin elevations in this patient population. In a targeted analysis of single nucleotide polymorphisms in genes known to be involved in bilirubin transport, no significant differences in allele frequency between high and low bilirubin toxicity grade were observed. CONCLUSIONS: These results indicate that risk for bilirubin elevation in patients receiving GS-9256 is unlikely to be strongly influenced by common genetic variants with large effects. The current study cannot rule out a role for common variants of weak effect, or a more complex model, including multiple contributing factors, such as rare variants and as yet unidentified environmental influences.

Mutation in TECPR2 reveals a role for autophagy in hereditary spastic paraparesis.

We studied five individuals from three Jewish Bukharian families affected by an apparently autosomal-recessive form of hereditary spastic paraparesis accompanied by severe intellectual disability, fluctuating central hypoventilation, gastresophageal reflux disease, wake apnea, areflexia, and unique dysmorphic features. Exome sequencing identified one homozygous variant shared among all affected individuals and absent in controls: a 1 bp frameshift TECPR2 deletion leading to a premature stop codon and predicting significant degradation of the protein. TECPR2 has been reported as a positive regulator of autophagy. We thus examined the autophagy-related fate of two key autophagic proteins, SQSTM1 (p62) and MAP1LC3B (LC3), in skin fibroblasts of an affected individual, as compared to a healthy control, and found that both protein levels were decreased and that there was a more pronounced decrease in the lipidated form of LC3 (LC3II). siRNA knockdown of TECPR2 showed similar changes, consistent with aberrant autophagy. Our results are strengthened by the fact that autophagy dysfunction has been implicated in a number of other neurodegenerative diseases. The discovered TECPR2 mutation implicates autophagy, a central intracellular mechanism, in spastic paraparesis.

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.

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.

Using ERDS to infer copy-number variants in high-coverage genomes.

Although there are many methods available for inferring copy-number variants (CNVs) from next-generation sequence data, there remains a need for a system that is computationally efficient but that retains good sensitivity and specificity across all types of CNVs. Here, we introduce a new method, estimation by read depth with single-nucleotide variants (ERDS), and use various approaches to compare its performance to other methods. We found that for common CNVs and high-coverage genomes, ERDS performs as well as the best method currently available (Genome STRiP), whereas for rare CNVs and high-coverage genomes, ERDS performs better than any available method. Importantly, ERDS accommodates both unique and highly amplified regions of the genome and does so without requiring separate alignments for calling CNVs and other variants. These comparisons show that for genomes sequenced at high coverage, ERDS provides a computationally convenient method that calls CNVs as well as or better than any currently available method.

Heritability of strabismus: genetic influence is specific to eso-deviation and independent of refractive error.

Strabismus represents a complex oculomotor disorder characterized by the deviation of one or both eyes and poor vision. A more sophisticated understanding of the genetic liability of strabismus is required to guide searches for associated molecular variants. In this classical twin study of 1,462 twin pairs, we examined the relative influence of genes and environment in comitant strabismus, and the degree to which these influences can be explained by factors in common with refractive error. Participants were examined for the presence of latent ('phoria') and manifest ('tropia') strabismus using cover-uncover and alternate cover tests. Two phenotypes were distinguished: eso-deviation (esophoria and esotropia) and exo-deviation (exophoria and exotropia). Structural equation modeling was subsequently employed to partition the observed phenotypic variation in the twin data into specific variance components. The prevalence of eso-deviation and exo-deviation was 8.6% and 20.7%, respectively. For eso-deviation, the polychoric correlation was significantly greater in monozygotic (MZ) (r = 0.65) compared to dizygotic (DZ) twin pairs (r = 0.33), suggesting a genetic role (p = .003). There was no significant difference in polychoric correlation between MZ (r = 0.55) and DZ twin pairs (r = 0.53) for exo-deviation (p = .86), implying that genetic factors do not play a significant role in the etiology of exo-deviation. The heritability of an eso-deviation was 0.64 (95% CI 0.50-0.75). The additive genetic correlation for eso-deviation and refractive error was 0.13 and the bivariate heritability (i.e., shared variance) was less than 1%, suggesting negligible shared genetic effect. This study documents a substantial heritability of 64% for eso-deviation, yet no corresponding heritability for exo-deviation, suggesting that the genetic contribution to strabismus may be specific to eso-deviation. Future studies are now needed to identify the genes associated with eso-deviation and unravel their mechanisms of action.

Exome sequencing followed by large-scale genotyping fails to identify single rare variants of large effect in idiopathic generalized epilepsy.

Idiopathic generalized epilepsy (IGE) is a complex disease with high heritability, but little is known about its genetic architecture. Rare copy-number variants have been found to explain nearly 3% of individuals with IGE; however, it remains unclear whether variants with moderate effect size and frequencies below what are reliably detected with genome-wide association studies contribute significantly to disease risk. In this study, we compare the exome sequences of 118 individuals with IGE and 242 controls of European ancestry by using next-generation sequencing. The exome-sequenced epilepsy cases include study subjects with two forms of IGE, including juvenile myoclonic epilepsy (n = 93) and absence epilepsy (n = 25). However, our discovery strategy did not assume common genetic control between the subtypes of IGE considered. In the sequence data, as expected, no variants were significantly associated with the IGE phenotype or more specific IGE diagnoses. We then selected 3,897 candidate epilepsy-susceptibility variants from the sequence data and genotyped them in a larger set of 878 individuals with IGE and 1,830 controls. Again, no variant achieved statistical significance. However, 1,935 variants were observed exclusively in cases either as heterozygous or homozygous genotypes. It is likely that this set of variants includes real risk factors. The lack of significant association evidence of single variants with disease in this two-stage approach emphasizes the high genetic heterogeneity of epilepsy disorders, suggests that the impact of any individual single-nucleotide variant in this disease is small, and indicates that gene-based approaches might be more successful for future sequencing studies of epilepsy predisposition.

Exome sequencing followed by large-scale genotyping suggests a limited role for moderately rare risk factors of strong effect in schizophrenia.

Schizophrenia is a severe psychiatric disorder with strong heritability and marked heterogeneity in symptoms, course, and treatment response. There is strong interest in identifying genetic risk factors that can help to elucidate the pathophysiology and that might result in the development of improved treatments. Linkage and genome-wide association studies (GWASs) suggest that the genetic basis of schizophrenia is heterogeneous. However, it remains unclear whether the underlying genetic variants are mostly moderately rare and can be identified by the genotyping of variants observed in sequenced cases in large follow-up cohorts or whether they will typically be much rarer and therefore more effectively identified by gene-based methods that seek to combine candidate variants. Here, we consider 166 persons who have schizophrenia or schizoaffective disorder and who have had either their genomes or their exomes sequenced to high coverage. From these data, we selected 5,155 variants that were further evaluated in an independent cohort of 2,617 cases and 1,800 controls. No single variant showed a study-wide significant association in the initial or follow-up cohorts. However, we identified a number of case-specific variants, some of which might be real risk factors for schizophrenia, and these can be readily interrogated in other data sets. Our results indicate that schizophrenia risk is unlikely to be predominantly influenced by variants just outside the range detectable by GWASs. Rather, multiple rarer genetic variants must contribute substantially to the predisposition to schizophrenia, suggesting that both very large sample sizes and gene-based association tests will be required for securely identifying genetic risk factors.

Personal receptor repertoires: olfaction as a model.

BACKGROUND: Information on nucleotide diversity along completely sequenced human genomes has increased tremendously over the last few years. This makes it possible to reassess the diversity status of distinct receptor proteins in different human individuals. To this end, we focused on the complete inventory of human olfactory receptor coding regions as a model for personal receptor repertoires. RESULTS: By performing data-mining from public and private sources we scored genetic variations in 413 intact OR loci, for which one or more individuals had an intact open reading frame. Using 1000 Genomes Project haplotypes, we identified a total of 4069 full-length polypeptide variants encoded by these OR loci, average of ~10 per locus, constituting a lower limit for the effective human OR repertoire. Each individual is found to harbor as many as 600 OR allelic variants, ~50% higher than the locus count. Because OR neuronal expression is allelically excluded, this has direct effect on smell perception diversity of the species. We further identified 244 OR segregating pseudogenes (SPGs), loci showing both intact and pseudogene forms in the population, twenty-six of which are annotatively "resurrected" from a pseudogene status in the reference genome. Using a custom SNP microarray we validated 150 SPGs in a cohort of 468 individuals, with every individual genome averaging 36 disrupted sequence variations, 15 in homozygote form. Finally, we generated a multi-source compendium of 63 OR loci harboring deletion Copy Number Variations (CNVs). Our combined data suggest that 271 of the 413 intact OR loci (66%) are affected by nonfunctional SNPs/indels and/or CNVs. CONCLUSIONS: These results portray a case of unusually high genetic diversity, and suggest that individual humans have a highly personalized inventory of functional olfactory receptors, a conclusion that might apply to other receptor multigene families.

Urinary norepinephrine and epinephrine excretion rates are heritable, but not associated with office and ambulatory blood pressure.

Genetic and environmental contributions to urinary excretion rates of norepinephrine (U(NE)V) and epinephrine (U(E)V) and their association with blood pressure (BP) were investigated in 91 African American (mean age, 17.3±2.6 years) and 101 European American (mean age, 18.7±3.4 years) mono- and di-zygotic twins. Genetic modeling was performed using Mx software. U(NE)V (1.9±1.3 µg h(-1)) and U(E)V (0.2±0.2 µg h(-1)) were highly correlated (r=0.81, P<0.001). Significant heritabilities for U(NE)V (0.68) and U(E)V (0.74) without ethnic and gender effects were observed. The genetic correlation between U(NE)V and U(E)V was 0.86. There was no clear pattern of correlations for U(NE)V and U(E)V with BP measures in European Americans, but African Americans showed some inverse correlations of moderate size. Measurements of U(NE)V and U(E)V provide a viable method for the study of sympathetic tone and are substantially heritable.

Genome-wide association study in Han Chinese identifies four new susceptibility loci for coronary artery disease.

We performed a meta-analysis of 2 genome-wide association studies of coronary artery disease comprising 1,515 cases and 5,019 controls followed by replication studies in 15,460 cases and 11,472 controls, all of Chinese Han ancestry. We identify four new loci for coronary artery disease that reached the threshold of genome-wide significance (P < 5 × 10(-8)). These loci mapped in or near TTC32-WDR35, GUCY1A3, C6orf10-BTNL2 and ATP2B1. We also replicated four loci previously identified in European populations (in or near PHACTR1, TCF21, CDKN2A-CDKN2B and C12orf51). These findings provide new insights into pathways contributing to the susceptibility for coronary artery disease in the Chinese Han population.

The association of genetic variants with hepatic steatosis in patients with genotype 1 chronic hepatitis C infection.

BACKGROUND: Single-nucleotide polymorphisms (SNPs) in the IL28B and PNPLA3 gene regions have been associated with hepatic steatosis in genotype 1 (G1) chronic HCV infection but their clinical impacts remain to be determined. AIM: We sought to validate these associations and to explore their impact on treatment response to peginterferon and ribavirin therapy. METHODS: A total of 972 G1 HCV-infected Caucasian patients were genotyped for the SNPs rs12979860 (IL28B) and rs2896019 (PNPLA3). Multivariable analysis tested IL28B and PNPLA3 for association with the presence of any steatosis (>0 %); clinically significant steatosis (>5 %); steatosis severity (grade 0-3/4); and the interacting associations of the SNPs and hepatic steatosis to sustained viral response (SVR). RESULTS: IL28B and PNPLA3 polymorphisms were associated with the presence of any steatosis (rs12979860, p = 1.87 × 10(-7); rs2896019, p = 7.56 × 10(-4)); clinically significant steatosis (rs12979860, p = 1.82 × 10(-3); rs2896019, p = 1.27 × 10(-4)); and steatosis severity (rs12979860, p = 2.05 × 10(-8); rs2896019, p = 2.62 × 10(-6)). Obesity, hypertriglyceridemia, hyperglycemia, liver fibrosis, and liver inflammation were all independently associated with worse steatosis. Hepatic steatosis was associated with lower SVR, and this effect was attenuated by IL28B. PNPLA3 had no independent association with SVR. CONCLUSIONS: IL28B and PNPLA3 are associated with hepatic steatosis prevalence and severity in Caucasians with G1 HCV, suggesting differing potential genetic risk pathways to steatosis. IL28B attenuates the association between steatosis and SVR. Remediable metabolic risk factors remain important, independently of these polymorphisms, and remain key therapeutic goals to achieve better outcomes for patients with HCV-associated hepatic steatosis.

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.