To complement our special issue on exome sequencing, Genome Biology asked several leaders in the field for their views on this new approach. Leslie G Biesecker (LGB), Jim C Mullikin (JM) and Kevin V Shianna (KVS) discuss the reasons for the popularity of exome sequencing and its contribution to genomics.
Human genetic variation contributes to differences in susceptibility to HIV-1 infection. To search for novel host resistance factors, we performed a genome-wide association study (GWAS) in hemophilia patients highly exposed to potentially contaminated factor VIII infusions.
Individuals with hemophilia A and a documented history of factor VIII infusions before the introduction of viral inactivation procedures (1979–1984) were recruited from 36 hemophilia treatment centers (HTCs), and their genome-wide genetic variants were compared with those from matched HIV-infected individuals. Homozygous carriers of known CCR5 resistance mutations were excluded. Single nucleotide polymorphisms (SNPs) and inferred copy number variants (CNVs) were tested using logistic regression. In addition, we performed a pathway enrichment analysis, a heritability analysis, and a search for epistatic interactions with CCR5 Δ32 heterozygosity.
A total of 560 HIV-uninfected cases were recruited: 36 (6.4%) were homozygous for CCR5 Δ32 or m303. After quality control and SNP imputation, we tested 1 081 435 SNPs and 3686 CNVs for association with HIV-1 serostatus in 431 cases and 765 HIV-infected controls. No SNP or CNV reached genome-wide significance. The additional analyses did not reveal any strong genetic effect.
Highly exposed, yet uninfected hemophiliacs form an ideal study group to investigate host resistance factors. Using a genome-wide approach, we did not detect any significant associations between SNPs and HIV-1 susceptibility, indicating that common genetic variants of major effect are unlikely to explain the observed resistance phenotype in this population.
Background & Aims
We studied the relationship between IL28B gene-related SNP rs12979860 and early viral kinetics (day 0–28) during peginterferon and ribavirin treatment, in 173 African Americans (AA) and 188 Caucasian Americans (CA) with HCV genotype 1.
We studied the relationship between IL28B 16 gene-related SNP rs12979860 and early viral kinetics (day 0–28) 17 during peginterferon and ribavirin treatment, in 171 African 18 Americans (AA) and 188 Caucasian Americans (CA) with HCV 19 genotype 1.
Compared to non-C/C genotypes, C/C was associated with greater declines in serum HCV RNA during phase 1 (day 0–2), phase 2 (day 7–28), and day 0–28 and higher response (undetected HCV RNA) rates at weeks 4 and 12 in AA and CA. A static phase and increases in HCV RNA from day 2 to 7 were more common in patients with non-C/C genotypes. C/C was also associated with higher week 24, 48, and 72 response rates in CA (p <0.01) but not in AA. At baseline, SNP genotype was the only independent predictor of phase 1; SNP genotype and phase 1 were independent predictors of phase 2 (p<0.001). There were no racial differences in HCV RNA declines during phase 1, day 2–7, phase 2, and day 0–28 with the same SNP genotype. AA with C/C and C/T genotypes had lower week 24, 48, and 72 (SVR) rates than did CA (p = 0.03). SNP C/C predicted higher SVR rates in AA and CA with high baseline HCV RNA (≥ 600,000 IU/ml), and in CA with ≥ 1 log10 IU/ml decrease in HCV RNA from day 0 to 28.
SNP rs12979860 is strongly associated with both phase 1 and phase 2 HCV RNA kinetics in AA and CA with HCV genotype 1.
IL28B gene; Hepatitis C virus; Viral kinetics; Single nucleotide polymorphism
Tibetan natives have lived on the Tibetan plateau (altitude ∼4,000 m) for at least 25,000 years, and as such they are adapted to life and reproduction in a hypoxic environment. Recent studies have identified two genetic loci, EGLN1 and EPAS1, that have undergone natural selection in Tibetans, and further demonstrated an association of EGLN1/EPAS1 genotype with hemoglobin concentration. Both genes encode major components of the hypoxia-inducible factor (HIF) transcriptional pathway, which coordinates an organism's response to hypoxia. Patients living at sea level with genetic disease of the HIF pathway have characteristic phenotypes at both the integrative-physiology and cellular level. We sought to test the hypothesis that natural selection to hypoxia within Tibetans results in related phenotypic differences. We compared Tibetans living at sea level with Han Chinese, who are Tibetans' most closely related major ethnic group. We found that Tibetans had a lower hemoglobin concentration, a higher pulmonary ventilation relative to metabolism, and blunted pulmonary vascular responses to both acute (minutes) and sustained (8 h) hypoxia. At the cellular level, the relative expression and hypoxic induction of HIF-regulated genes were significantly lower in peripheral blood lymphocytes from Tibetans compared with Han Chinese. Within the Tibetans, we found a significant correlation between both EPAS1 and EGLN1 genotype and the induction of erythropoietin by hypoxia. In conclusion, this study provides further evidence that Tibetans respond less vigorously to hypoxic challenge. This is evident at sea level and, at least in part, appears to arise from a hyporesponsive HIF transcriptional system.
hypoxic pulmonary vasoconstriction; EPAS1; EGLN1; erythropoietin; high altitude
Genetic variation in the IL28B region has been associated with sustained virological response (SVR) rates in chronic hepatitis C (CHC) patients treated with peginterferon-α and ribavirin. We hypothesized that IL28B polymorphism is associated with intrahepatic expression of interferon-stimulated genes (ISGs), known to influence treatment outcome. IL28B genotyping (rs12979860) and whole-genome RNA expression were performed using liver biopsies from 61 North American CHC patients. After correction for multiple testing (false discovery rate < 0.10), 164 transcripts were found to be differentially expressed by IL28B-type. The interferon signaling pathway was the most enriched canonical pathway differentially expressed by IL28B-type (p < 10−5), with most genes showing higher expression in livers of individuals carrying the poor-response IL28B-type. In 25 patients for which treatment response data were available, IL28B-type was associated with SVR (p = 0.0054). ISG expression was also associated with SVR; however, this was not independent of IL28B-type. Analysis of miR-122 expression in liver biopsies showed reduced miR-122 levels associated with poorer treatment outcome, independently of IL28B-type. No association was observed between IL28B-type and levels of liver IL28B or IL28A mRNA expression. IL28B protein sequence variants associated with rs12979860 were therefore investigated in vitro: no differences in ISG induction or inhibition of HCV replication were observed in Huh7.5 cells.
The good response IL28B variant was strongly associated with lower level ISG expression. The results suggest that IL28B genotype may explain the relationship between hepatic ISG expression and HCV treatment outcome, and this is independent of miR-122 expression. IL28B-type was not associated with intrahepatic IL28B mRNA expression in vivo. Further investigation of the precise molecular mechanism(s) by which IL28B genetic variation influences HCV outcomes is warranted.
interferon lambda; hepatitis C virus; gene expression
A single-nucleotide polymorphism (rs2395029) in the HCP5 gene associated with HLA-B*5701 is correlated with lower HIV-1 viral set point. The two allelic forms of coding region were ectopically expressed in TZM-bl cells for an effect on HIV-1 replication. No significant HIV-1 restriction was observed in the cells with infectivity assays throughout HIV-1 life cycle, suggesting that the association of HCP5 variant with viral control is likely due to HLA-B*5701-related effect or other functional variants in the haplotype or both.
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-α (peg-IFN)-related thrombocytopenia, neutropenia, and leukopenia.
1604/3070 patients in the IDEAL study consented to genetic testing. Trial inclusion criteria included a platelet (Pl) count ≥80 × 109/L and an absolute neutrophil count (ANC) ≥ 1500/mm3. 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).
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.
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.
GWAS; ITPA; Thrombocytopenia; Hepatitis C; Neutropenia; IL28B
Alternating hemiplegia of childhood (AHC) is a rare, severe neurodevelopmental syndrome characterized by recurrent hemiplegic episodes and distinct neurologic manifestations. AHC is usually a sporadic disorder with unknown etiology. Using exome sequencing of seven patients with AHC, and their unaffected parents, we identified de novo nonsynonymous mutations in ATP1A3 in all seven AHC patients. Subsequent sequence analysis of ATP1A3 in 98 additional patients revealed that 78% of AHC cases have a likely causal ATP1A3 mutation, including one inherited mutation in a familial case of AHC. Remarkably, six ATP1A3 mutations explain the majority of patients, including one observed in 36 patients. Unlike ATP1A3 mutations that cause rapid-onset-dystonia-parkinsonism, AHC-causing mutations revealed consistent reductions in ATPase activity without effects on protein expression. This work identifies de novo ATP1A3 mutations as the primary cause of AHC, and offers insight into disease pathophysiology by expanding the spectrum of phenotypes associated with mutations in this gene.
Background. A recent genome-wide association study reported a strong association with a single-nucleotide polymorphism (SNP) in the inosine triphosphate (ITPA) gene and hemolytic anemia in patients infected with hepatitis C virus (HCV) receiving pegylated interferon and ribavirin. We investigate these polymorphisms in a cohort of human immunodeficiency virus (HIV)/HCV–coinfected patients.
Methods. DNA was available for 161 patients with validated outcomes. We analyzed the association between the variants and week 4 hemoglobin reduction. Anemia over the course of therapy, ribavirin (RBV) dose reduction, serum RBV level, and rapid virological response (RVR) and sustained virological response (SVR) were also investigated. Using a candidate gene approach, ITPA variants rs1127354 and rs7270101 were tested using the ABI TaqMan kit. Multivariable models were used to identify predictors of anemia.
Results. A significant minority (33%) of patients were predicted to have reduced ITPase activity. The minor allele of each variant was associated with protection against week 4 anemia. In multivariable models only the genetic variants, creatinine, and zidovudine exposure remained significant. ITPase deficiency was not associated with RBV-dose reduction, RVR, or SVR.
Conclusions. This study confirms that polymorphisms in the ITPA gene are associated with protection from RBV-induced anemia in HIV/HCV-coinfected patients but not improved clinical outcomes.
We sequenced the genomes of ten unrelated individuals and identified heterozygous stop gain variants in protein-coding genes: we then sequenced their transcriptomes and assessed the expression levels of the stop gain alleles. An ANOVA showed statistically significant differences between their expression levels (p=4×10-16). This difference was almost entirely accounted for by whether the stop gain variant had a second, non-protein-truncating function in or near an alternate transcript: stop gains without alternate functions were generally not found in the cDNA (p=3×10-5). Additionally, stop gain variants in two intronless genes were not expressed, an unexpected outcome given previous studies. In this study, stop gain variants were either well expressed in all individuals or were never expressed. Our finding that stop gain variants were generally expressed only when they had an alternate function suggests that most naturally occurring stop gain variants in protein-coding genes are either not transcribed or have their transcripts destroyed.
Nonsense-mediated decay; whole-genome sequencing; RNA-Seq; premature termination codons
An association between carbamazepine-induced hypersensitivity and HLA-A*3101 has been reported in populations of both European and Asian descent. We aimed to investigate HLA-A*3101 and other common variants across the genome as markers for cutaneous adverse drug reactions (cADRs) attributed to lamotrigine and phenytoin.
Materials & methods
We recruited patients with lamotrigine-induced cADRs (n = 46) and patients with phenytoin-cADRs (n = 44) and the 1958 British birth cohort was used as a control (n = 1296). HLA-A*3101 was imputed from genome-wide association study data. We applied genome-wide association to study lamotrigine- and phenytoin-induced cADR, and total cADR cases combined.
Neither HLA-A*3101 nor any other genetic marker significantly predicted lamotrigine- or phenytoin-induced cADRs.
HLA-A*3101 does not appear to be a predictor for lamotrigine- and phenytoin-induced cADRs in Europeans. Our genome-wide association study results do not support the existence of a clinically relevant common variant for the development of lamotrigine- or phenytoin-induced cADRs. As a predictive marker, HLA-A*3101 appears to be specific for carbamazepine-induced cADRs.
epilepsy; GWAS; HLA-A*3101; hypersensitivity; lamotrigine; phenytoin
Mendelian analysis of disorders of immune regulation can provide insight into molecular pathways associated with host defense and immune tolerance.
We identified three families with a dominantly inherited complex of cold-induced urticaria, antibody deficiency, and susceptibility to infection and autoimmunity. Immunophenotyping methods included flow cytometry, analysis of serum immunoglobulins and autoantibodies, lymphocyte stimulation, and enzymatic assays. Genetic studies included linkage analysis, targeted Sanger sequencing, and next-generation whole-genome sequencing.
Cold urticaria occurred in all affected subjects. Other, variable manifestations included atopy, granulomatous rash, autoimmune thyroiditis, the presence of antinuclear antibodies, sinopulmonary infections, and common variable immunodeficiency. Levels of serum IgM and IgA and circulating natural killer cells and class-switched memory B cells were reduced. Linkage analysis showed a 7-Mb candidate interval on chromosome 16q in one family, overlapping by 3.5 Mb a disease-associated haplotype in a smaller family. This interval includes PLCG2, encoding phospholipase Cγ2 (PLCγ2), a signaling molecule expressed in B cells, natural killer cells, and mast cells. Sequencing of complementary DNA revealed heterozygous transcripts lacking exon 19 in two families and lacking exons 20 through 22 in a third family. Genomic sequencing identified three distinct in-frame deletions that cosegregated with disease. These deletions, located within a region encoding an autoinhibitory domain, result in protein products with constitutive phospholipase activity. PLCG2-expressing cells had diminished cellular signaling at 37°C but enhanced signaling at subphysiologic temperatures.
Genomic deletions in PLCG2 cause gain of PLCγ2 function, leading to signaling abnormalities in multiple leukocyte subsets and a phenotype encompassing both excessive and deficient immune function. (Funded by the National Institutes of Health Intramural Research Programs and others.)
There is considerable interest in the use of next-generation sequencing to help diagnose unidentified genetic conditions, but it is difficult to predict the success rate in a clinical setting that includes patients with a broad range of phenotypic presentations.
The authors present a pilot programme of whole-exome sequencing on 12 patients with unexplained and apparent genetic conditions, along with their unaffected parents. Unlike many previous studies, the authors did not seek patients with similar phenotypes, but rather enrolled any undiagnosed proband with an apparent genetic condition when predetermined criteria were met.
This undertaking resulted in a likely genetic diagnosis in 6 of the 12 probands, including the identification of apparently causal mutations in four genes known to cause Mendelian disease (TCF4, EFTUD2, SCN2A and SMAD4) and one gene related to known Mendelian disease genes (NGLY1). Of particular interest is that at the time of this study, EFTUD2 was not yet known as a Mendelian disease gene but was nominated as a likely cause based on the observation of de novo mutations in two unrelated probands. In a seventh case with multiple disparate clinical features, the authors were able to identify homozygous mutations in EFEMP1 as a likely cause for macular degeneration (though likely not for other features).
This study provides evidence that next-generation sequencing can have high success rates in a clinical setting, but also highlights key challenges. It further suggests that the presentation of known Mendelian conditions may be considerably broader than currently recognised.
Exome sequencing; unidentified genetic conditions; medical genetics; paediatrics; clinical genetics; complex traits; genetic screening/counselling; genetics; genome-wide; psychotic disorders (including schizophrenia); molecular genetics; gastroenterology; immunology (including allergy).
Understanding how human genetic variation impacts individual response to immunogens is fundamental for rational vaccine development. To explore host mechanisms involved in cellular immune responses to the MRKAd5 human immunodeficiency virus type 1 (HIV-1) gag/pol/nef vaccine tested in the Step trial, we performed a genome-wide association study of determinants of HIV-specific T cell responses, measured by interferon γ enzyme-linked immunospot assays. No human genetic variant reached genome-wide significance, but polymorphisms located in the major histocompatibility complex (MHC) region showed the strongest association with response to the HIV-1 Gag protein: HLA-B alleles known to be associated with differences in HIV-1 control were responsible for these associations. The implication of the same HLA alleles in vaccine-induced cellular immunity and in natural immune control is of relevance for vaccine design. Furthermore, our results demonstrate the importance of considering the host immunogenetic background in the analysis of immune responses to T cell vaccines.
To date, CCR5 variants remain the only human genetic factors to be confirmed to impact HIV-1 acquisition. However, protective CCR5 variants are largely absent in African populations, in which sporadic resistance to HIV-1 infection is still unexplained. Here we perform a genome-wide association study (GWAS) in a population of 1,532 individuals from Malawi, a country with high prevalence of HIV-1 infection, to investigate whether common genetic variants associate with HIV-1 susceptibility in Africans. Using single nucleotide polymorphisms (SNPs) present on the genome-wide chip, we also investigated previously reported associations with HIV-1 susceptibility or acquisition. Recruitment was coordinated by the Center for HIV/AIDS Vaccine Immunology at two sexually transmitted infection clinics. HIV status was determined by HIV rapid tests and nucleic acid testing.
After quality control, the population consisted of 848 high-risk seronegative and 531 HIV-1 seropositive individuals. Logistic regression testing in an additive genetic model was performed for SNPs that passed quality control. No single SNP yielded a significant P-value after correction for multiple testing. The study was sufficiently powered to detect markers with genotype relative risk ≥ 2.0 and minor allele frequencies ≥12%. This is the first GWAS of host determinants of HIV-1 susceptibility, performed in an African population. The absence of any significant association can have many possible explanations: rarer genetic variants or common variants with weaker effect could be responsible for the resistance phenotype; alternatively, resistance to HIV-1 infection might be due to non-genetic parameters or to complex interactions between genes, immunity and environment.
Human immunodeficiency virus (HIV-1); acquisition; resistance; Genome Wide Association Study (GWAS); Africa
Polymorphisms of IL28B gene are highly associated with sustained virological response (SVR) in patients with chronic hepatitis C treated with peginterferon and ribavirin. Quantitation of Interferon-γ Inducible Protein-10 (IP-10) may also differentiate antiviral response. We evaluated IP-10 levels in pretreatment serum from 115 non-responders and 157 sustained responders in the VIRAHEP-C cohort, including African Americans (AA) and Caucasian Americans (CA). Mean IP-10 was lower in sustained responders compared to non-responders (460 ± 37 pg/ml vs 697 ± 49 pg/ml, p<0.001), both in AA and CA. The positive predictive value of low IP-10 levels (<600 pg/ml) for SVR was 69% while the negative predictive value of high IP-10 levels (>600 pg/ml) was 67%. We assessed the combination of pretreatment IP-10 levels with IL28B genotype as predictors of treatment response. The IL28B polymorphism rs12979860 was tested in 210 participants. CC, CT, or TT genotypes were found in 30%, 49%, and 21%, respectively, with corresponding SVR rates of 87%, 50%, and 39% (p<0.0001). Serum IP-10 levels within the IL28B genotype groups provided additional information regarding the likelihood of SVR (p< 0.0001). CT carriers with low IP-10 had 64% SVR versus 24% with high IP-10. Similarly, a higher SVR rate was identified for TT and CC carriers with low versus high IP-10 (TT: 48% versus 20%, CC: 89% versus 79%). IL28B genotype and baseline IP-10 levels were additive but independent when predicting SVR in both AA and CA.
When IL28B genotype is combined with pretreatment serum IP-10 measurement, the predictive value for discrimination between SVR and non-response is significantly improved, especially in non-CC genotypes. This relationship warrants further investigation to elucidate mechanisms of antiviral response and prospective validation.
Peginterferon; Ribavirin; Race; Biomarker; Genotype
The authors that the number of activating and inhibitory KIR genes varies between individuals and plays a role in the regulation of immune mechanisms that determine 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.
There is marked intrinsic variation in the extent to which individuals are able to control HIV-1. We have identified a genetic copy number variable region (CNV) in humans that plays a significant role in the control of HIV-1. This CNV is located in the genomic region that encodes the killer cell immunoglobulin-like receptors (KIRs) and specifically affects the KIR3DS1 and KIR3DL1 genes, encoding two KIRs that interact with human leukocyte antigen B (HLA-B) ligands. KIRs are expressed on the surface of natural killer (NK) cells, which serve as important players in the innate immune response, and are involved in the recognition of infected and malignant cells through a loss or alteration in “self” ligands. We use both genetic association and functional evidence to show a strong interaction between KIR3DL1 and KIR3DS1, indicating that increasing gene counts for KIR3DL1 confer increasing levels of protection against HIV-1, but only in the presence of at least one copy of KIR3DS1. This effect was associated with a dramatic increase in the abundance of KIR3DS1+ NK cells in the peripheral blood, and strongly associated with a more robust capacity of peripheral NK cells to suppress HIV-1 replication in vitro. This work provides one of the few examples of an association between a relatively common CNV and a human complex trait.
Carbamazepine causes various forms of hypersensitivity reactions, ranging from maculopapular exanthema to severe blistering reactions. The HLA-B★1502 allele has been shown to be strongly correlated with carbamazepine-induced Stevens–Johnson syndrome and toxic epidermal necrolysis (SJS–TEN) in the Han Chinese and other Asian populations but not in European populations.
We performed a genomewide association study of samples obtained from 22 subjects with carbamazepine-induced hypersensitivity syndrome, 43 subjects with carbamazepine-induced maculopapular exanthema, and 3987 control subjects, all of European descent. We tested for an association between disease and HLA alleles through proxy single-nucleotide polymorphisms and imputation, confirming associations by high-resolution sequence-based HLA typing. We replicated the associations in samples from 145 subjects with carbamazepine-induced hypersensitivity reactions.
The HLA-A★3101 allele, which has a prevalence of 2 to 5% in Northern European populations, was significantly associated with the hypersensitivity syndrome (P = 3.5×10−8). An independent genomewide association study of samples from subjects with maculopapular exanthema also showed an association with the HLA-A★3101 allele (P = 1.1×10−6). Follow-up genotyping confirmed the variant as a risk factor for the hypersensitivity syndrome (odds ratio, 12.41; 95% confidence interval [CI], 1.27 to 121.03), maculopapular exanthema (odds ratio, 8.33; 95% CI, 3.59 to 19.36), and SJS–TEN (odds ratio, 25.93; 95% CI, 4.93 to 116.18).
The presence of the HLA-A★3101 allele was associated with carbamazepine-induced hypersensitivity reactions among subjects of Northern European ancestry. The presence of the allele increased the risk from 5.0% to 26.0%, whereas its absence reduced the risk from 5.0% to 3.8%. (Funded by the U.K. Department of Health and others.)
Summary: Here we present Sequence Variant Analyzer (SVA), a software tool that assigns a predicted biological function to variants identified in next-generation sequencing studies and provides a browser to visualize the variants in their genomic contexts. SVA also provides for flexible interaction with software implementing variant association tests allowing users to consider both the bioinformatic annotation of identified variants and the strength of their associations with studied traits. We illustrate the annotation features of SVA using two simple examples of sequenced genomes that harbor Mendelian mutations.
Availability and implementation: Freely available on the web at http://www.svaproject.org.
Supplementary information: Supplementary data are available at Bioinformatics online.
BACKGROUND & AIMS
In a genome-wide association study of patients being treated for chronic hepatitis C, 2 functional variants in ITPA that cause inosine triphosphatase (ITPase) deficiency were shown to protect against ribavirin (RBV)-induced hemolytic anemia during early stages of treatment. We aimed to replicate this finding in an independent cohort from the Study of Viral Resistance to Antiviral Therapy of Chronic Hepatitis C and to investigate the effects of these variants beyond week 4.
Genetic material was available from 318 patients. The ITPA variants, rs1127354 (exon 2, P32T) and rs7270101 (intron 2, splice altering), were genotyped and tested for association with hemoglobin (Hb) reduction at week 4. An ITPase deficiency variable was defined that combined both ITPA variants according to documented effect on ITPase activity. We investigated the impact of ITPA variants on Hb levels over the course of therapy and on the need for RBV dose reduction.
The final analysis included 304 patients with genotype 1 hepatitis C virus (167 white patients and 137 black patients). The polymorphisms rs1127354 and rs7270101 were associated with Hb reduction at week 4 (P = 3.1 × 10−13 and 1.3 × 10−3, respectively). The minor alleles of each variant protected against Hb reduction. Combining the variants into the ITPase deficiency variable strengthened the association (P = 2.4 × 10−18). The ITPase deficiency variable was associated with lower rates of anemia over the entire treatment period (48 weeks), as well as a lower rate of anemia-related RBV dose reduction (hazard ratio, 0.52; P = .0037). No association with sustained virological response was observed.
Two polymorphisms that cause ITPase deficiency are strongly associated with protection from RBV-induced hemolytic anemia and decrease the need for RBV dose reduction.
Pharmacogenomics; Genome-Wide Association Study; Polymorphism; Single Nucleotide Polymorphism; HCV; Adverse Event
We performed a whole-genome association study on HIV-1 viral load setpoint in an African American cohort (n=515), and an intronic SNP in the HLA-B gene showed one of the strongest associations. Using a subset of patients, we show that this SNP reflects the effect of the HLA-B*5703 allele, which shows a genome-wide significant association with HIV-1 VL setpoint (p=5.6×10−10). These analyses therefore confirm a member of the HLA-B*57 group of alleles as the most important common variant influencing viral load variation in African Americans, consistent with what is observed in individuals of European ancestry in which the most important common variant is HLA-B*5701.
HIV; viral load setpoint; host genetics; association study; HLA
Schizophrenia is a complex disorder, caused by both genetic and environmental factors and their interactions. Research on pathogenesis has traditionally focused on neurotransmitter systems in the brain, particularly those involving dopamine. Schizophrenia has been considered a separate disease for over a century, but in the absence of clear biological markers, diagnosis has historically been based on signs and symptoms. A fundamental message emerging from genome-wide association studies of copy number variations (CNVs) associated with the disease is that its genetic basis does not necessarily conform to classical nosological disease boundaries. Certain CNVs confer not only high relative risk of schizophrenia but also of other psychiatric disorders1–3. The structural variations associated with schizophrenia can involve several genes and the phenotypic syndromes, or the ‘genomic disorders’, have not yet been characterized4. Single nucleotide polymorphism (SNP)-based genome-wide association studies with the potential to implicate individual genes in complex diseases may reveal underlying biological pathways. Here we combined SNP data from several large genome-wide scans and followed up the most significant association signals. We found significant association with several markers spanning the major histocompatibility complex (MHC) region on chromosome 6p21.3-22.1, a marker located upstream of the neurogranin gene (NRGN) on 11q24.2 and a marker in intron four of transcription factor 4 (TCF4) on 18q21.2. Our findings implicating the MHC region are consistent with an immune component to schizophrenia risk, whereas the association with NRGN and TCF4 points to perturbation of pathways involved in brain development, memory and cognition.
Psychiatric disorders such as schizophrenia are commonly accompanied by cognitive impairments that are treatment resistant and crucial to functional outcome. There has been great interest in studying cognitive measures as endophenotypes for psychiatric disorders, with the hope that their genetic basis will be clearer. To investigate this, we performed a genome-wide association study involving 11 cognitive phenotypes from the Cambridge Neuropsychological Test Automated Battery. We showed these measures to be heritable by comparing the correlation in 100 monozygotic and 100 dizygotic twin pairs. The full battery was tested in ∼750 subjects, and for spatial and verbal recognition memory, we investigated a further 500 individuals to search for smaller genetic effects. We were unable to find any genome-wide significant associations with either SNPs or common copy number variants. Nor could we formally replicate any polymorphism that has been previously associated with cognition, although we found a weak signal of lower than expected P-values for variants in a set of 10 candidate genes. We additionally investigated SNPs in genomic loci that have been shown to harbor rare variants that associate with neuropsychiatric disorders, to see if they showed any suggestion of association when considered as a separate set. Only NRXN1 showed evidence of significant association with cognition. These results suggest that common genetic variation does not strongly influence cognition in healthy subjects and that cognitive measures do not represent a more tractable genetic trait than clinical endpoints such as schizophrenia. We discuss a possible role for rare variation in cognitive genomics.
This is a crucial transition time for human genetics in general, and for HIV host genetics in particular. After years of equivocal results from candidate gene analyses, several genome-wide association studies have been published that looked at plasma viral load or disease progression. Results from other studies that used various large-scale approaches (siRNA screens, transcriptome or proteome analysis, comparative genomics) have also shed new light on retroviral pathogenesis. However, most of the inter-individual variability in response to HIV-1 infection remains to be explained: genome resequencing and systems biology approaches are now required to progress toward a better understanding of the complex interactions between HIV-1 and its human host.