A fast algorithm for genome-wide haplotype pattern mining

<p>Abstract</p> <p>Background</p> <p>Identifying the genetic components of common diseases has long been an important area of research. Recently, genotyping technology has reached the level where it is cost effective to genotype single nucleotide polymorphism (SNP) markers covering the entire genome, in thousands of individuals, and analyse such data for markers associated with a diseases. The statistical power to detect association, however, is limited when markers are analysed one at a time. This can be alleviated by considering multiple markers simultaneously. The <it>Haplotype Pattern Mining </it>(HPM) method is a machine learning approach to do exactly this.</p> <p>Results</p> <p>We present a new, faster algorithm for the HPM method. The new approach use patterns of haplotype diversity in the genome: locally in the genome, the number of observed haplotypes is much smaller than the total number of possible haplotypes. We show that the new approach speeds up the HPM method with a factor of 2 on a genome-wide dataset with 5009 individuals typed in 491208 markers using default parameters and more if the pattern length is increased.</p> <p>Conclusion</p> <p>The new algorithm speeds up the HPM method and we show that it is feasible to apply HPM to whole genome association mapping with thousands of individuals and hundreds of thousands of markers.</p

Ash generation and distribution from the April-May 2010 eruption of Eyjafjallajökull, Iceland

The 39-day long eruption at the summit of Eyjafjallajökull volcano in April–May 2010 was of modest size but ash was widely dispersed. By combining data from ground surveys and remote sensing we show that the erupted material was 4.8±1.2·1011 kg (benmoreite and trachyte, dense rock equivalent volume 0.18±0.05 km3). About 20% was lava and water-transported tephra, 80% was airborne tephra (bulk volume 0.27 km3) transported by 3–10 km high plumes. The airborne tephra was mostly fine ash (diameter <1000 µm). At least 7·1010 kg (70 Tg) was very fine ash (<28 µm), several times more than previously estimated via satellite retrievals. About 50% of the tephra fell in Iceland with the remainder carried towards south and east, detected over ~7 million km2 in Europe and the North Atlantic. Of order 1010 kg (2%) are considered to have been transported longer than 600–700 km with <108 kg (<0.02%) reaching mainland Europe

Quantitative detection of myocardial ischaemia by stress echocardiography; a comparison with SPECT

<p>Abstract</p> <p>Aims</p> <p>Real-time perfusion (RTP) adenosine stress echocardiography (ASE) can be used to visually evaluate myocardial ischaemia. The RTP power modulation technique angio-mode (AM), provides images for off-line perfusion quantification using Qontrast^®software, generating values of peak signal intensity (A), myocardial blood flow velocity (β) and myocardial blood flow (Axβ). By comparing rest and stress values, their respective reserve values (A-r, β-r, Axβ-r) are generated. We evaluated myocardial ischaemia by RTP-ASE Qontrast^®quantification, compared to visual perfusion evaluation with ^99mTc-tetrofosmin single-photon emission computed tomography (SPECT).</p> <p>Methods and Results</p> <p>Patients admitted to SPECT underwent RTP-ASE (SONOS 5500) using AM during Sonovue^®infusion, before and throughout adenosine stress, also used for SPECT. Visual myocardial perfusion and wall motion analysis, and Qontrast^®quantification, were blindly compared to one another and to SPECT, at different time points off-line.</p> <p>We analyzed 201 coronary territories (left anterior descendent [LAD], left circumflex [LCx] and right coronary [RCA] artery territories) in 67 patients. SPECT showed ischaemia in 18 patients and 19 territories. Receiver operator characteristics and kappa values showed significant agreement with SPECT only for β-r and Axβ-r in all segments: area under the curve 0.678 and 0.665; P < 0.001 and < 0.01, respectively. The closest agreements were seen in the LAD territory: kappa 0.442 for both β-r and Axβ-r; P < 0.01. Visual evaluation of ischaemia showed good agreement with SPECT: accuracy 93%; kappa 0.67; P < 0.001; without non-interpretable territories.</p> <p>Conclusion</p> <p>In this agreement study with SPECT, RTP-ASE Qontrast^®quantification of myocardial ischaemia was less accurate and less feasible than visual evaluation and needs further development to be clinically useful.</p

Head to head comparisons of two modalities of perfusion adenosine stress echocardiography with simultaneous SPECT

<p>Abstract</p> <p>Background</p> <p>Real-time perfusion (RTP) contrast echocardiography can be used during adenosine stress echocardiography (ASE) to evaluate myocardial ischemia. We compared two different types of RTP power modulation techniques, angiomode (AM) and high-resolution grayscale (HR), with ^99mTc-tetrofosmin single-photon emission computed tomography (SPECT) for the detection of myocardial ischemia.</p> <p>Methods</p> <p>Patients with known or suspected coronary artery disease (CAD), admitted to SPECT, were prospectively invited to participate. Patients underwent RTP imaging (SONOS 5500) using AM and HR during Sonovue^®infusion, before and throughout the adenosine stress, also used for SPECT. Analysis of myocardial perfusion and wall motion by RTP-ASE were done for AM and HR at different time points, blinded to one another and to SPECT. Each segment was attributed to one of the three main coronary vessel areas of interest.</p> <p>Results</p> <p>In 50 patients, 150 coronary areas were analyzed by SPECT and RTP-ASE AM and HR. SPECT showed evidence of ischemia in 13 out of 50 patients. There was no significant difference between AM and HR in detecting ischemia (p = 0.08). The agreement for AM and HR, compared to SPECT, was 93% and 96%, with Kappa values of 0.67 and 0.75, respectively (p < 0.001).</p> <p>Conclusion</p> <p>There was no significant difference between AM and HR in correctly detecting myocardial ischemia as judged by SPECT. This suggests that different types of RTP modalities give comparable data during RTP-ASE in patients with known or suspected CAD.</p

A Flexible and Accurate Genotype Imputation Method for the Next Generation of Genome-Wide Association Studies

Genotype imputation methods are now being widely used in the analysis of genome-wide association studies. Most imputation analyses to date have used the HapMap as a reference dataset, but new reference panels (such as controls genotyped on multiple SNP chips and densely typed samples from the 1,000 Genomes Project) will soon allow a broader range of SNPs to be imputed with higher accuracy, thereby increasing power. We describe a genotype imputation method (IMPUTE version 2) that is designed to address the challenges presented by these new datasets. The main innovation of our approach is a flexible modelling framework that increases accuracy and combines information across multiple reference panels while remaining computationally feasible. We find that IMPUTE v2 attains higher accuracy than other methods when the HapMap provides the sole reference panel, but that the size of the panel constrains the improvements that can be made. We also find that imputation accuracy can be greatly enhanced by expanding the reference panel to contain thousands of chromosomes and that IMPUTE v2 outperforms other methods in this setting at both rare and common SNPs, with overall error rates that are 15%–20% lower than those of the closest competing method. One particularly challenging aspect of next-generation association studies is to integrate information across multiple reference panels genotyped on different sets of SNPs; we show that our approach to this problem has practical advantages over other suggested solutions

GWAS Meets Microarray: Are the Results of Genome-Wide Association Studies and Gene-Expression Profiling Consistent? Prostate Cancer as an Example

Genome-wide association studies (GWASs) and global profiling of gene expression (microarrays) are two major technological breakthroughs that allow hypothesis-free identification of candidate genes associated with tumorigenesis. It is not obvious whether there is a consistency between the candidate genes identified by GWAS (GWAS genes) and those identified by profiling gene expression (microarray genes).We used the Cancer Genetic Markers Susceptibility database to retrieve single nucleotide polymorphisms from candidate genes for prostate cancer. In addition, we conducted a large meta-analysis of gene expression data in normal prostate and prostate tumor tissue. We identified 13,905 genes that were interrogated by both GWASs and microarrays. On the basis of P values from GWASs, we selected 1,649 most significantly associated genes for functional annotation by the Database for Annotation, Visualization and Integrated Discovery. We also conducted functional annotation analysis using same number of the top genes identified in the meta-analysis of the gene expression data. We found that genes involved in cell adhesion were overrepresented among both the GWAS and microarray genes.We conclude that the results of these analyses suggest that combining GWAS and microarray data would be a more effective approach than analyzing individual datasets and can help to refine the identification of candidate genes and functions associated with tumor development

Planck intermediate results. XLI. A map of lensing-induced B-modes

The secondary cosmic microwave background (CMB) $B$-modes stem from the post-decoupling distortion of the polarization $E$-modes due to the gravitational lensing effect of large-scale structures. These lensing-induced $B$-modes constitute both a valuable probe of the dark matter distribution and an important contaminant for the extraction of the primary CMB $B$-modes from inflation. Planck provides accurate nearly all-sky measurements of both the polarization $E$-modes and the integrated mass distribution via the reconstruction of the CMB lensing potential. By combining these two data products, we have produced an all-sky template map of the lensing-induced $B$-modes using a real-space algorithm that minimizes the impact of sky masks. The cross-correlation of this template with an observed (primordial and secondary) $B$-mode map can be used to measure the lensing $B$-mode power spectrum at multipoles up to $2000$. In particular, when cross-correlating with the $B$-mode contribution directly derived from the Planck polarization maps, we obtain lensing-induced $B$-mode power spectrum measurement at a significance level of $12\,\sigma$, which agrees with the theoretical expectation derived from the Planck best-fit $\Lambda$CDM model. This unique nearly all-sky secondary $B$-mode template, which includes the lensing-induced information from intermediate to small ($10\lesssim \ell\lesssim 1000$) angular scales, is delivered as part of the Planck 2015 public data release. It will be particularly useful for experiments searching for primordial $B$-modes, such as BICEP2/Keck Array or LiteBIRD, since it will enable an estimate to be made of the lensing-induced contribution to the measured total CMB $B$-modes.Comment: 20 pages, 12 figures; Accepted for publication in A&A; The B-mode map is part of the PR2-2015 Cosmology Products; available as Lensing Products in the Planck Legacy Archive http://pla.esac.esa.int/pla/#cosmology; and described in the 'Explanatory Supplement' https://wiki.cosmos.esa.int/planckpla2015/index.php/Specially_processed_maps#2015_Lensing-induced_B-mode_ma

Selective Autophagy of Mitochondria on a Ubiquitin-Endoplasmic-Reticulum Platform

The dynamics and coordination between autophagy machinery and selective receptors during mitophagy are unknown. Also unknown is whether mitophagy depends on pre-existing membranes or is triggered on the surface of damaged mitochondria. Using a ubiquitin-dependent mitophagy inducer, the lactone ivermectin, we have combined genetic and imaging experiments to address these questions. Ubiquitination of mitochondrial fragments is required the earliest, followed by auto-phosphorylation of TBK1. Next, early essential autophagy proteins FIP200 and ATG13 act at different steps, whereas ULK1 and ULK2 are dispensable. Receptors act temporally and mechanistically upstream of ATG13 but downstream of FIP200. The VPS34 complex functions at the omegasome step. ATG13 and optineurin target mitochondria in a discontinuous oscillatory way, suggesting multiple initiation events. Targeted ubiquitinated mitochondria are cradled by endoplasmic reticulum (ER) strands even without functional autophagy machinery and mitophagy adaptors. We propose that damaged mitochondria are ubiquitinated and dynamically encased in ER strands, providing platforms for formation of the mitophagosomes

Role of BRCA gene dysfunction in breast and ovarian cancer predisposition

Tumor suppressor genes that perform apparently generic cellular functions nonetheless cause tissue-specific syndromes in the human population when they are mutated in the germline. The two major hereditary breast/ovarian cancer predisposition genes, BRCA1 and BRCA2, appear to participate in a common pathway that is involved in the control of homologous recombination and in the maintenance of genomic integrity. How might such functions translate into the specific suppression of cancers of the breast and ovarian epithelia? Recent advances in the study of BRCA1 and BRCA2, discussed herein, have provided new opportunities to address this question

Genetic Determinants of Cardiovascular Events among Women with Migraine: A Genome-Wide Association Study

Migraine is associated with an increased risk for cardiovascular disease (CVD). Both migraine and CVD are highly heritable. However, the genetic liability for CVD among migraineurs is unclear.We performed a genome-wide association study for incident CVD events during 12 years of follow-up among 5,122 migraineurs participating in the population-based Women's Genome Health Study. Migraine was self-reported and CVD events were confirmed after medical records review. We calculated odds ratios (OR) and 95% confidence intervals (CI) and considered a genome-wide p-value <5×10(-8) as significant.Among the 5,122 women with migraine 164 incident CVD events occurred during follow-up. No SNP was associated with major CVD, ischemic stroke, myocardial infarction, or CVD death at the genome-wide level; however, five SNPs showed association with p<5×10(-6). Among migraineurs with aura rs7698623 in MEPE (OR = 6.37; 95% CI 3.15-12.90; p = 2.7×10(-7)) and rs4975709 in IRX4 (OR = 5.06; 95% CI 2.66-9.62; p = 7.7×10(-7)) appeared to be associated with ischemic stroke, rs2143678 located close to MDF1 with major CVD (OR = 3.05; 95% CI 1.98-4.69; p = 4.3×10(-7)), and the intergenic rs1406961 with CVD death (OR = 12.33; 95% CI 4.62-32.87; p = 5.2×10(-7)). Further, rs1047964 in BACE1 appeared to be associated with CVD death among women with any migraine (OR = 4.67; 95% CI 2.53-8.62; p = 8.0×10(-7)).Our results provide some suggestion for an association of five SNPs with CVD events among women with migraine; none of the results was genome-wide significant. Four associations appeared among migraineurs with aura, two of those with ischemic stroke. Although our population is among the largest with migraine and incident CVD information, these results must be treated with caution, given the limited number of CVD events among women with migraine and the low minor allele frequencies for three of the SNPs. Our results await independent replication and should be considered hypothesis generating for future research