Abanish Singh

The Characterization of Twenty Sequenced Human Genomes

We present the analysis of twenty human genomes to evaluate the prospects for identifying rare functional variants that contribute to a phenotype of interest. We sequenced at high coverage ten “case” genomes from individuals with severe hemophilia A and ten “control” genomes. We summarize the number of genetic variants emerging from a study of this magnitude, and provide a proof of concept for the identification of rare and highly-penetrant functional variants by confirming that the cause of hemophilia A is easily recognizable in this data set. We also show that the number of novel single nucleotide variants (SNVs) discovered per genome seems to stabilize at about 144,000 new variants per genome, after the first 15 individuals have been sequenced. Finally, we find that, on average, each genome carries 165 homozygous protein-truncating or stop loss variants in genes representing a diverse set of pathways.

Screening the human exome: a comparison of whole genome and whole transcriptome sequencing

BackgroundThere is considerable interest in the development of methods to efficiently identify all coding variants present in large sample sets of humans. There are three approaches possible: whole-genome sequencing, whole-exome sequencing using exon capture methods, and RNA-Seq. While whole-genome sequencing is the most complete, it remains sufficiently expensive that cost effective alternatives are important.ResultsHere we provide a systematic exploration of how well RNA-Seq can identify human coding variants by comparing variants identified through high coverage whole-genome sequencing to those identified by high coverage RNA-Seq in the same individual. This comparison allowed us to directly evaluate the sensitivity and specificity of RNA-Seq in identifying coding variants, and to evaluate how key parameters such as the degree of coverage and the expression levels of genes interact to influence performance. We find that although only 40% of exonic variants identified by whole genome sequencing were captured using RNA-Seq; this number rose to 81% when concentrating on genes known to be well-expressed in the source tissue. We also find that a high false positive rate can be problematic when working with RNA-Seq data, especially at higher levels of coverage.ConclusionsWe conclude that as long as a tissue relevant to the trait under study is available and suitable quality control screens are implemented, RNA-Seq is a fast and inexpensive alternative approach for finding coding variants in genes with sufficiently high expression levels.

SVA: software for annotating and visualizing sequenced human genomes

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. Contact: d.ge@duke.edu Supplementary information: Supplementary data are available at Bioinformatics online.

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.

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.

Gene by stress genome-wide interaction analysis and path analysis identify EBF1 as a cardiovascular and metabolic risk gene.

We performed gene–environment interaction genome-wide association analysis (G × E GWAS) to identify SNPs whose effects on metabolic traits are modified by chronic psychosocial stress in the Multi-Ethnic Study of Atherosclerosis (MESA). In Whites, the G × E GWAS for hip circumference identified five SNPs within the Early B-cell Factor 1 (EBF1) gene, all of which were in strong linkage disequilibrium. The gene-by-stress interaction (SNP × STRESS) term P-values were genome-wide significant (Ps=7.14E−09 to 2.33E−08, uncorrected; Ps=1.99E−07 to 5.18E−07, corrected for genomic control). The SNP-only (without interaction) model P-values (Ps=0.011–0.022) were not significant at the conventional genome-wide significance level. Further analysis of related phenotypes identified gene-by-stress interaction effects for waist circumference, body mass index (BMI), fasting glucose, type II diabetes status, and common carotid intimal–medial thickness (CCIMT), supporting a proposed model of gene-by-stress interaction that connects cardiovascular disease (CVD) risk factor endophenotypes such as central obesity and increased blood glucose or diabetes to CVD itself. Structural equation path analysis suggested that the path from chronic psychosocial stress to CCIMT via hip circumference and fasting glucose was larger (estimate=0.26, P=0.033, 95% CI=0.02–0.49) in the EBF1 rs4704963 CT/CC genotypes group than the same path in the TT group (estimate=0.004, P=0.34, 95% CI=−0.004–0.012). We replicated the association of the EBF1 SNPs and hip circumference in the Framingham Offspring Cohort (gene-by-stress term P-values=0.007–0.012) as well as identified similar path relationships. This observed and replicated interaction between psychosocial stress and variation in the EBF1 gene may provide a biological hypothesis for the complex relationship between psychosocial stress, central obesity, diabetes, and cardiovascular disease.

Socioeconomic indices as independent correlates of C-reactive protein in the national longitudinal study of adolescent health

Objectives To examine the association between socioeconomic status (SES) and C-reactive protein (CRP) to understand how SES may increase the risk of cardiovascular disease and thus identify targets for prevention measures. Methods Path models were used to examine direct and indirect associations of four indices of SES (objective early life built environment ratings, parental and participant education, and income) with CRP measured during early adulthood using data from the National Longitudinal Adolescent Health Study (n = 11,371; mean age = 29 years, range = 24–32 years; 53.8% women, 28.0% black participants). The present study examined potential mediation of the association of SES with CRP by way of body mass index (BMI), smoking, and alcohol consumption within white and black men and women. Results BMI was a mediator of the relation between parent education and CRP for white men (path coefficient [&ggr;] = −0.05, p < .001) and women (&ggr; = −0.05, p < .001). Smoking mediated the income-CRP (&ggr; = −0.01, p < .01) and the education-CRP (&ggr; = −0.07, p < .001) relation for white men. BMI mediated the relation between all measures of SES and CRP for white women (&ggr; values between −0.02 and −0.05; p values < .01). None of the risk factors mediated the SES-CRP relation in black participants. Conclusions These findings indicate that the association of SES with CRP is influenced by both the timing and type of SES measure examined. In addition, race and sex play a role in how potential mediators are involved with the SES-CRP relationship, such that BMI and smoking were mediators in white men, whereas BMI was the sole mediator in white women.

A Functional Polymorphism in the 5HTR2C Gene Associated with Stress Responses Also Predicts Incident Cardiovascular Events

Previously we have shown that a functional nonsynonymous single nucleotide polymorphism (rs6318) of the 5HTR2C gene located on the X-chromosome is associated with hypothalamic-pituitary-adrenal axis response to a stress recall task, and with endophenotypes associated with cardiovascular disease (CVD). These findings suggest that individuals carrying the rs6318 Ser23 C allele will be at higher risk for CVD compared to Cys23 G allele carriers. The present study examined allelic variation in rs6318 as a predictor of coronary artery disease (CAD) severity and a composite endpoint of all-cause mortality or myocardial infarction (MI) among Caucasian participants consecutively recruited through the cardiac catheterization laboratory at Duke University Hospital (Durham, NC) as part of the CATHGEN biorepository. Study population consisted of 6,126 Caucasian participants (4,036 [65.9%] males and 2,090 [34.1%] females). A total of 1,769 events occurred (1,544 deaths and 225 MIs; median follow-up time =  5.3 years, interquartile range  = 3.3–8.2). Unadjusted Cox time-to-event regression models showed, compared to Cys23 G carriers, males hemizygous for Ser23 C and females homozygous for Ser23C were at increased risk for the composite endpoint of all-cause death or MI: Hazard Ratio (HR)  = 1.47, 95% confidence interval (CI)  = 1.17, 1.84, p  = .0008. Adjusting for age, rs6318 genotype was not related to body mass index, diabetes, hypertension, dyslipidemia, smoking history, number of diseased coronary arteries, or left ventricular ejection fraction in either males or females. After adjustment for these covariates the estimate for the two Ser23 C groups was modestly attenuated, but remained statistically significant: HR  = 1.38, 95% CI = 1.10, 1.73, p = .005. These findings suggest that this functional polymorphism of the 5HTR2C gene is associated with increased risk for CVD mortality and morbidity, but this association is apparently not explained by the association of rs6318 with traditional risk factors or conventional markers of atherosclerotic disease.

A Putatively Functional Polymorphism in the HTR2C Gene is Associated with Depressive Symptoms in White Females Reporting Significant Life Stress

Psychosocial stress is well known to be positively associated with subsequent depressive symptoms. Cortisol response to stress may be one of a number of biological mechanisms that links psychological stress to depressive symptoms, although the precise causal pathway remains unclear. Activity of the x-linked serotonin 5-HTR2C receptor has also been shown to be associated with depression and with clinical response to antidepressant medications. We recently demonstrated that variation in a single nucleotide polymorphism on the HTR2C gene, rs6318 (Ser23Cys), is associated with different cortisol release and short-term changes in affect in response to a series of stress tasks in the laboratory. Based on this observation, we decided to examine whether rs6318 might moderate the association between psychosocial stress and subsequent depressive symptoms. In the present study we use cross-sectional data from a large population-based sample of young adult White men (N = 2,366) and White women (N = 2,712) in the United States to test this moderation hypothesis. Specifically, we hypothesized that the association between self-reported stressful life events and depressive symptoms would be stronger among homozygous Ser23 C females and hemizygous Ser23 C males than among Cys23 G carriers. In separate within-sex analyses a genotype-by-life stress interaction was observed for women (p = .022) but not for men (p = .471). Homozygous Ser23 C women who reported high levels of life stress had depressive symptom scores that were about 0.3 standard deviations higher than female Cys23 G carriers with similarly high stress levels. In contrast, no appreciable difference in depressive symptoms was observed between genotypes at lower levels of stress. Our findings support prior work that suggests a functional SNP on the HTR2C gene may confer an increased risk for depressive symptoms in White women with a history of significant life stress.

Chronic family stress moderates the association between a TOMM40 variant and triglyceride levels in two independent Caucasian samples

TOMM40 SNP rs157580 has been associated with triglyceride levels in genome-wide association studies (GWAS). Chronic caregiving stress moderates the association between triglyceride levels and a nearby SNP rs439401 that is associated with triglyceride levels in GWAS. Here, we report data from two independent Caucasian samples (242 U.S. women and men; 466 Danish men) testing the hypothesis that chronic family stress also moderates the association between rs157580 and triglyceride levels. The interaction of rs157580 and family stress in predicting triglyceride levels was statistically significant in the U.S. sample (p=0.004) and marginally significant (p=0.075) in the Danish sample. The G allele of rs157580 was associated with increased triglyceride levels among family stressed cases in both samples compared with A/A cases, but not among controls. Chronic family stress moderates the association of rs157580 variants with triglyceride levels and should be taken into account for disease risk assessment and potential intervention.