Quantifying single nucleotide variant detection sensitivity in exome sequencing

BACKGROUND: The targeted capture and sequencing of genomic regions has rapidly demonstrated its utility in genetic studies. Inherent in this technology is considerable heterogeneity of target coverage and this is expected to systematically impact our sensitivity to detect genuine polymorphisms. To fully interpret the polymorphisms identified in a genetic study it is often essential to both detect polymorphisms and to understand where and with what probability real polymorphisms may have been missed. RESULTS: Using down-sampling of 30 deeply sequenced exomes and a set of gold-standard single nucleotide variant (SNV) genotype calls for each sample, we developed an empirical model relating the read depth at a polymorphic site to the probability of calling the correct genotype at that site. We find that measured sensitivity in SNV detection is substantially worse than that predicted from the naive expectation of sampling from a binomial. This calibrated model allows us to produce single nucleotide resolution SNV sensitivity estimates which can be merged to give summary sensitivity measures for any arbitrary partition of the target sequences (nucleotide, exon, gene, pathway, exome). These metrics are directly comparable between platforms and can be combined between samples to give “power estimates” for an entire study. We estimate a local read depth of 13X is required to detect the alleles and genotype of a heterozygous SNV 95% of the time, but only 3X for a homozygous SNV. At a mean on-target read depth of 20X, commonly used for rare disease exome sequencing studies, we predict 5–15% of heterozygous and 1–4% of homozygous SNVs in the targeted regions will be missed. CONCLUSIONS: Non-reference alleles in the heterozygote state have a high chance of being missed when commonly applied read coverage thresholds are used despite the widely held assumption that there is good polymorphism detection at these coverage levels. Such alleles are likely to be of functional importance in population based studies of rare diseases, somatic mutations in cancer and explaining the “missing heritability” of quantitative traits

A meta-analysis of CAG (cytarabine, aclarubicin, G-CSF) regimen for the treatment of 1029 patients with acute myeloid leukemia and myelodysplastic syndrome

The regimen of cytarabine, aclarubicin and G-CSF (CAG) has been widely used in China and Japan for treatment of acute myeloid leukemia (AML) and myelodysplastic syndrome (MDS). We searched literature on CAG between 1995 and 2010 and performed a meta-analysis to determine its overall efficacy using a random-effects or fixed-effects model. Thirty five trials with a total of 1029 AML (n = 814) and MDS (n = 215) patients were included for analysis. The CR rate of AML (57.9%) was significantly higher than that of MDS (45.7%) (p < 0.01). No difference in CR was noted between the new (56.7%) and relapsed/refractory AML (60.1%) (p > 0.05). The CR rate was also significantly higher in patients with favorable (64.5%) and intermediate (69.6%) karyotypes than those with unfavorable one (29.5%) (p < 0.05). Remarkably, the CR rate of CAG was significantly higher than those of non-CAG regimens (odds ratio 2.43). CAG regimen was well tolerated, with cardiotoxicity in 2.3% and early death in 5.2% of the cases. In conclusion, CAG regimen was an effective and safe regimen for the treatment of AML, and may be more effective than non-CAG regimens. Randomized controlled trials are strongly recommended to evaluate its efficacy and safety in comparison with the current standard treatment

Identification of rare de novo epigenetic variations in congenital disorders

Certain human traits such as neurodevelopmental disorders (NDs) and congenital anomalies (CAs) are believed to be primarily genetic in origin. However, even after whole-genome sequencing (WGS), a substantial fraction of such disorders remain unexplained. We hypothesize that some cases of ND-CA are caused by aberrant DNA methylation leading to dysregulated genome function. Comparing DNA methylation profiles from 489 individuals with ND-CAs against 1534 controls, we identify epivariations as a frequent occurrence in the human genome. De novo epivariations are significantly enriched in cases, while RNAseq analysis shows that epivariations often have an impact on gene expression comparable to loss-of-function mutations. Additionally, we detect and replicate an enrichment of rare sequence mutations overlapping CTCF binding sites close to epivariations, providing a rationale for interpreting non-coding variation. We propose that epivariations contribute to the pathogenesis of some patients with unexplained ND-CAs, and as such likely have diagnostic relevance.The authors are grateful to the patients and families who participated in this study and to the collaborators who supported patient recruitment. This work was supported by NIH grant HG006696 and research grant 6-FY13-92 from the March of Dimes to A.J.S., grant HL098123 to B.D.G. and A.J.S., Gulbenkian Programme for Advanced Medical Education and the Portuguese Foundation for Science and Technology (SFRH/BDINT/51549/ 2011, PIC/IC/83026/2007, PIC/IC/83013/2007, SFRH/BD/90167/2012, Portugal) to P.M., F.L., and M.B., by the Northern Portugal Regional Operational Programme (NORTE 2020), under the Portugal 2020 Partnership Agreement, through the European Regional Development Fund (FEDER) (NORTE-01-0145-FEDER-000013) to P.M., a Beatriu de Pinos Postdoctoral Fellowship to R.S.J. (2011BP-A00515), and a Seaver Foundation fellowship to S.D.R. The views expressed are those of the authors and do not necessarily reflect those of the National Heart, Lung, and Blood Institute or the National Institutes of Health. Research reported in this paper was supported by the Office of Research Infrastructure of the National Institutes of Health under award number S10OD018522. This work was supported in part through the computational resources and staff expertise provided by Scientific Computing at the Icahn School of Medicine at Mount Sinai.The authors are grateful to the patients and families who participated in this study and to the collaborators who supported patient recruitment. This work was supported by NIH grant HG006696 and research grant 6-FY13-92 from the March of Dimes to A.J.S., grant HL098123 to B.D.G. and A.J.S., Gulbenkian Programme for Advanced Medical Education and the Portuguese Foundation for Science and Technology (SFRH/BDINT/51549/ 2011, PIC/IC/83026/2007, PIC/IC/83013/2007, SFRH/BD/90167/2012, Portugal) to P.M., F.L., and M.B., by the Northern Portugal Regional Operational Programme (NORTE 2020), under the Portugal 2020 Partnership Agreement, through the European Regional Development Fund (FEDER) (NORTE-01-0145-FEDER-000013) to P.M., a Beatriu de Pinos Postdoctoral Fellowship to R.S.J. (2011BP-A00515), and a Seaver Foundation fellowship to S.D.R. The views expressed are those of the authors and do not necessarily reflect those of the National Heart, Lung, and Blood Institute or the National Institutes of Health. Research reported in this paper was supported by the Office of Research Infrastructure of the National Institutes of Health under award number S10OD018522. This work was supported in part through the computational resources and staff expertise provided by Scientific Computing at the Icahn School of Medicine at Mount Sinai

β-defensin 1 expression in HCV infected liver/liver cancer: an important role in protecting HCV progression and liver cancer development

Abstract β-defensin family plays a role in host defense against viral infection, however its role in HCV infection is still unknown. In this study, we demonstrated that β-defensin 1 was significantly reduced in HCV-infected liver specimens. Treatment with interferon and ribavirin upregulated β-defensin-1, but not other β-defensin tested, with the extent and duration of upregulation associated with treatment response. We investigated β-defensin family expression in liver cancer in publicly available datasets and found that among all the β-defensins tested, only β-defensin 1 was significantly downregulated, suggesting β-defensin 1 plays a crucial role in liver cancer development. Further analysis identified E-cadherin as the top positive correlated gene, while hepatocyte growth factor-regulated tyrosine kinase substrate as the top negative correlated gene. Expression of two proteoglycans were also positively correlated with that of β-defensin 1. We have also identified small molecules as potential therapeutic agents to reverse β-defensin 1-associated gene signature. Furthermore, the downregulation of β-defensin 1 and E-cadherin, and upregulation of hepatocyte growth factor-regulated tyrosine kinase substrate, were further confirmed in liver cancer and adjacent normal tissue collected from in-house Chinese liver cancer patients. Together, our results suggest β-defensin 1 plays an important role in protecting HCV progression and liver cancer development

Locus of emotion influences psychophysiological reactions to music

It is now widely accepted that the perception of emotional expression in music can be vastly different from the feelings evoked by it. However, less understood is how the locus of emotion affects the experience of music, that is how the act of perceiving the emotion in music compares with the act of assessing the emotion induced in the listener by the music. In the current study, we compared these two emotion loci based on the psychophysiological response of 40 participants listening to 32 musical excerpts taken from movie soundtracks. Facial electromyography, skin conductance, respiration and heart rate were continuously measured while participants were required to assess either the emotion expressed by, or the emotion they felt in response to the music. Using linear mixed effects models, we found a higher mean response in psychophysiological measures for the “perceived” than the “felt” task. This result suggested that the focus on one’s self distracts from the music, leading to weaker bodily reactions during the “felt” task. In contrast, paying attention to the expression of the music and consequently to changes in timbre, loudness and harmonic progression enhances bodily reactions. This study has methodological implications for emotion induction research using psychophysiology and the conceptualization of emotion loci. Firstly, different tasks can elicit different psychophysiological responses to the same stimulus and secondly, both tasks elicit bodily responses to music. The latter finding questions the possibility of a listener taking on a purely cognitive mode when evaluating emotion expression

The effects of β-glucan on human immune and cancer cells

Non-prescriptional use of medicinal herbs among cancer patients is common around the world. The alleged anti-cancer effects of most herbal extracts are mainly based on studies derived from in vitro or in vivo animal experiments. The current information suggests that these herbal extracts exert their biological effect either through cytotoxic or immunomodulatory mechanisms. One of the active compounds responsible for the immune effects of herbal products is in the form of complex polysaccharides known as β-glucans. β-glucans are ubiquitously found in both bacterial or fungal cell walls and have been implicated in the initiation of anti-microbial immune response. Based on in vitro studies, β-glucans act on several immune receptors including Dectin-1, complement receptor (CR3) and TLR-2/6 and trigger a group of immune cells including macrophages, neutrophils, monocytes, natural killer cells and dendritic cells. As a consequence, both innate and adaptive response can be modulated by β-glucans and they can also enhance opsonic and non-opsonic phagocytosis. In animal studies, after oral administration, the specific backbone 1→3 linear β-glycosidic chain of β-glucans cannot be digested. Most β-glucans enter the proximal small intestine and some are captured by the macrophages. They are internalized and fragmented within the cells, then transported by the macrophages to the marrow and endothelial reticular system. The small β-glucans fragments are eventually released by the macrophages and taken up by other immune cells leading to various immune responses. However, β-glucans of different sizes and branching patterns may have significantly variable immune potency. Careful selection of appropriate β-glucans is essential if we wish to investigate the effects of β-glucans clinically. So far, no good quality clinical trial data is available on assessing the effectiveness of purified β-glucans among cancer patients. Future effort should direct at performing well-designed clinical trials to verify the actual clinical efficacy of β-glucans or β-glucans containing compounds

Pooled sequencing of 531 genes in inflammatory bowel disease identifies an associated rare variant in BTNL2 and implicates other immune related genes.

The contribution of rare coding sequence variants to genetic susceptibility in complex disorders is an important but unresolved question. Most studies thus far have investigated a limited number of genes from regions which contain common disease associated variants. Here we investigate this in inflammatory bowel disease by sequencing the exons and proximal promoters of 531 genes selected from both genome-wide association studies and pathway analysis in pooled DNA panels from 474 cases of Crohn's disease and 480 controls. 80 variants with evidence of association in the sequencing experiment or with potential functional significance were selected for follow up genotyping in 6,507 IBD cases and 3,064 population controls. The top 5 disease associated variants were genotyped in an extension panel of 3,662 IBD cases and 3,639 controls, and tested for association in a combined analysis of 10,147 IBD cases and 7,008 controls. A rare coding variant p.G454C in the BTNL2 gene within the major histocompatibility complex was significantly associated with increased risk for IBD (p = 9.65x10-10, OR = 2.3[95% CI = 1.75-3.04]), but was independent of the known common associated CD and UC variants at this locus. Rare (T) or decreased risk (IL12B p.V298F, and NICN p.H191R) of IBD. These results provide additional insights into the involvement of the inhibition of T cell activation in the development of both sub-phenotypes of inflammatory bowel disease. We suggest that although rare coding variants may make a modest overall contribution to complex disease susceptibility, they can inform our understanding of the molecular pathways that contribute to pathogenesis

DNA methylation signature of chronic low-grade inflammation and its role in cardio-respiratory diseases

We performed a multi-ethnic Epigenome Wide Association study on 22,774 individuals to describe the DNA methylation signature of chronic low-grade inflammation as measured by C-Reactive protein (CRP). We find 1,511 independent differentially methylated loci associated with CRP. These CpG sites show correlation structures across chromosomes, and are primarily situated in euchromatin, depleted in CpG islands. These genomic loci are predominantly situated in transcription factor binding sites and genomic enhancer regions. Mendelian randomization analysis suggests altered CpG methylation is a consequence of increased blood CRP levels. Mediation analysis reveals obesity and smoking as important underlying driving factors for changed CpG methylation. Finally, we find that an activated CpG signature significantly increases the risk for cardiometabolic diseases and COPD

Trans-ancestry genome-wide association study identifies 12 genetic loci influencing blood pressure and implicates a role for DNA methylation

We carried out a trans-ancestry genome-wide association and replication study of blood pressure phenotypes among up to 320,251 individuals of East Asian, European and South Asian ancestry. We find genetic variants at 12 new loci to be associated with blood pressure (P = 3.9 &times; 10-11 to 5.0 &times; 10-21). The sentinel blood pressure SNPs are enriched for association with DNA methylation at multiple nearby CpG sites, suggesting that, at some of the loci identified, DNA methylation may lie on the regulatory pathway linking sequence variation to blood pressure. The sentinel SNPs at the 12 new loci point to genes involved in vascular smooth muscle (IGFBP3, KCNK3, PDE3A and PRDM6) and renal (ARHGAP24, OSR1, SLC22A7 and TBX2) function. The new and known genetic variants predict increased left ventricular mass, circulating levels of NT-proBNP, and cardiovascular and all-cause mortality (P = 0.04 to 8.6 &times; 10-6). Our results provide new evidence for the role of DNA methylation in blood pressure regulation

A Large-Scale Multi-ancestry Genome-wide Study Accounting for Smoking Behavior Identifies Multiple Significant Loci for Blood Pressure

Genome-wide association analysis advanced understanding of blood pressure (BP), a major risk factor for vascular conditions such as coronary heart disease and stroke. Accounting for smoking behavior may help identify BP loci and extend our knowledge of its genetic architecture. We performed genome-wide association meta-analyses of systolic and diastolic BP incorporating gene-smoking interactions in 610,091 individuals. Stage 1 analysis examined similar to 18.8 million SNPs and small insertion/deletion variants in 129,913 individuals from four ancestries (European, African, Asian, and Hispanic) with follow-up analysis of promising variants in 480,178 additional individuals from five ancestries. We identified 15 loci that were genome-wide significant (p <5 x 10(-8)) in stage 1 and formally replicated in stage 2. A combined stage 1 and 2 meta-analysis identified 66 additional genome-wide significant loci (13, 35, and 18 loci in European, African, and trans-ancestry, respectively). A total of 56 known BP loci were also identified by our results (p <5 x 10(-8)). Of the newly identified loci, ten showed significant interaction with smoking status, but none of them were replicated in stage 2. Several loci were identified in African ancestry, highlighting the importance of genetic studies in diverse populations. The identified loci show strong evidence for regulatory features and support shared pathophysiology with cardiometabolic and addiction traits. They also highlight a role in BP regulation for biological candidates such as modulators of vascular structure and function (CDKN1B, BCAR1-CFDP1, PXDN, EEA1), ciliopathies (SDCCAG8, RPGRIP1L), telomere maintenance (TNKS, PINX1, AKTIP), and central dopaminergic signaling MSRA, EBF2).Peer reviewe