introduction to plink
play

Introduction to PLINK Scott Hazelhurst Sydney Brenner Institute for - PowerPoint PPT Presentation

May 24, 2023 •746 likes •958 views

Introduction Introduction to PLINK Scott Hazelhurst Sydney Brenner Institute for Molecular Bioscience and School of Electrical & Information Engineering University of the Witwatersrand Johannesburg 2014 Scott Hazelhurst Sydney Brenner

  1. Introduction Introduction to PLINK Scott Hazelhurst Sydney Brenner Institute for Molecular Bioscience and School of Electrical & Information Engineering University of the Witwatersrand Johannesburg 2014 Scott Hazelhurst Sydney Brenner Institute for Molecular Bioscience and School of Electrical & Information E Introduction to PLINK

  2. Introduction Data format Standard tool for manipulating genotype data vcftools PLINK/PSEQ Plink has multiple data formats Other tools for converting to/from other formats pngu.mgh.harvard.edu/~purcell/plink/ https://www.cog-genomics.org/plink2/ Scott Hazelhurst Sydney Brenner Institute for Molecular Bioscience and School of Electrical & Information Eng Introduction to PLINK

  3. Introduction PLINK in transition to PLINK 2 Current version of plink: 1.90b2 Previous version: 1.07 New version: Much faster Has more features Missing some features Data compatible Scott Hazelhurst Sydney Brenner Institute for Molecular Bioscience and School of Electrical & Information Eng Introduction to PLINK

  4. Introduction PLINK primarily aimed at genotype data SNPs “short” indels Some support for CNV A leading tool for GWAS, structure analysis – many other tools support format. Not appropriate for many SVs, or when great variability Scott Hazelhurst Sydney Brenner Institute for Molecular Bioscience and School of Electrical & Information Eng Introduction to PLINK

  5. Introduction PED format PED files with individuals information, MAP file with SNP information PED file One row per individual. Columns are Family ID, Individual ID Paternal ID, Maternal ID Sex (1=male; 2=female; other=unknown) Phenotype Missing: − 9 , 0; Control: 1; Case 2. (or QT) Pair of columns per SNP: different encodings possible Scott Hazelhurst Sydney Brenner Institute for Molecular Bioscience and School of Electrical & Information E Introduction to PLINK

  6. Introduction HCB181 1 0 0 1 1 2 2 2 2 2 2 1 2 2 2 2 2 HCB182 1 0 0 1 1 2 2 1 2 2 2 1 2 1 2 2 2 HCB183 1 0 0 1 2 2 2 1 2 2 2 1 2 1 1 2 2 HCB184 1 0 0 1 1 2 2 1 2 2 2 1 1 2 2 2 2 HCB185 1 0 0 1 1 2 2 1 2 2 2 2 2 2 2 2 2 HCB186 1 0 0 1 1 2 2 2 2 2 2 1 1 2 2 2 2 HCB187 1 0 0 1 1 2 2 2 2 2 2 1 2 1 2 2 2 HCB188 1 0 0 1 1 2 2 1 2 2 2 1 1 2 2 2 2 HCB189 1 0 0 1 1 2 2 2 2 2 2 2 2 2 2 2 2 HCB190 1 0 0 1 1 2 2 2 2 2 2 2 2 2 2 2 2 HCB191 1 0 0 1 2 1 2 2 2 2 2 1 2 1 2 2 2 HCB192 1 0 0 1 1 2 2 2 2 2 2 1 1 2 2 2 2 HCB193 1 0 0 1 1 1 2 2 2 2 2 2 2 2 2 2 2 HCB194 1 0 0 1 1 2 2 2 2 2 2 1 2 2 2 2 2 Scott Hazelhurst Sydney Brenner Institute for Molecular Bioscience and School of Electrical & Information Eng HCB195 1 0 0 1 1 2 2 2 2 2 2 2 2 2 2 2 2 Introduction to PLINK

  7. Introduction Can be used to model family studies AFAM 1 0 0 . . . AFAM 2 0 0 . . . AFAM 3 1 2 . . . AFAM 4 1 2 . . . AFAM 5 0 0 . . . AFAM 6 1 5 . . . AFAM 7 0 0 . . . AFAM 8 3 0 . . . AFAM 9 0 4 . . . Scott Hazelhurst Sydney Brenner Institute for Molecular Bioscience and School of Electrical & Information Eng Introduction to PLINK

  8. Introduction NB: Some commands/toolsx: Expect sex information by default --allow-no-sex / --must-have-sex Want phenotype data Scott Hazelhurst Sydney Brenner Institute for Molecular Bioscience and School of Electrical & Information Eng Introduction to PLINK

  9. Introduction MAP file MAP file has one row per SNP Chromosome: 1..26 (X, Y, XY, MT) SNP id genetic distance (morgans) base pair (which build!) Newer versions of PLINK have support for some non-human genomes Scott Hazelhurst Sydney Brenner Institute for Molecular Bioscience and School of Electrical & Information E Introduction to PLINK

  10. Introduction 1 rs3094315 0 742429 1 rs3131972 0 742584 1 rs12562034 0 758311 1 rs12124819 0 766409 1 rs11240777 0 788822 1 rs6681049 0 789870 1 rs4970383 0 828418 1 rs4475691 0 836671 1 rs7537756 0 844113 1 rs13302982 0 851671 1 rs1110052 0 863421 1 rs2272756 0 871896 Scott Hazelhurst Sydney Brenner Institute for Molecular Bioscience and School of Electrical & Information Eng Introduction to PLINK

  11. Introduction Binary PED format Faster, more compact FAM file: one row per individual – identification information (first 6 columns of PED file). Human readable BIM file: one row per SNP. MAP file + two variants for that SNP. Human readable. BED file: one row per individual – genotype information (rest of the columns of the PED file). Not human readable Don’t confuse with UCSC BED format for genomic data – can have both in a study Scott Hazelhurst Sydney Brenner Institute for Molecular Bioscience and School of Electrical & Information Eng Introduction to PLINK

  12. Introduction BIM file 1 rs2185539 0 566875 C C 1 rs11510103 0 567753 A A 1 rs11240767 0 728951 C C 1 rs3131972 0 752721 G G 1 rs3131969 0 754182 G G 1 rs1048488 0 760912 T T 1 rs12562034 0 768448 A G 1 rs12124819 0 776546 A A 1 rs4040617 0 779322 A A 1 rs2905036 0 792480 T T 1 rs4245756 0 799463 C C 1 rs12086311 0 808769 G G Scott Hazelhurst Sydney Brenner Institute for Molecular Bioscience and School of Electrical & Information Eng Introduction to PLINK

  13. Introduction Other formats: transposed long format Not commonly used – typically when you need to import from another format. May be easy to write a script that does the conversion. Scott Hazelhurst Sydney Brenner Institute for Molecular Bioscience and School of Electrical & Information Eng Introduction to PLINK

  14. Introduction Transposed data tped/tfam files. tped: one row per SNP with SNP info followed by genotype of each individual; tfam: info about individuals Scott Hazelhurst Sydney Brenner Institute for Molecular Bioscience and School of Electrical & Information Eng Introduction to PLINK

  15. Introduction Long format Very inefficient – but may be useful in conversion MAP file FAM file LGEN file containing genotypes LGEN family ID, individual ID SNP ID two alleles Scott Hazelhurst Sydney Brenner Institute for Molecular Bioscience and School of Electrical & Information E Introduction to PLINK

  16. Introduction A 1 rs123 A C A 2 rs28782 C G A 3 rs919878 T T A 2 rs123 A C B 7 rs123 A C B 8 rs123 A C B 9 rs28782 C T Scott Hazelhurst Sydney Brenner Institute for Molecular Bioscience and School of Electrical & Information Eng Introduction to PLINK

  17. Introduction Phenotype/Cluster file FID IID PHE FID IID PHE FID IID PHE FID IID PHE Scott Hazelhurst Sydney Brenner Institute for Molecular Bioscience and School of Electrical & Information Eng Introduction to PLINK

  18. Introduction Can have multiple phenotypes FID IID PHE1 PHE2 PHE3 FID IID PHE1 PHE2 PHE3 FID IID PHE1 PHE2 PHE3 FID IID PHE1 PHE2 PHE3 Scott Hazelhurst Sydney Brenner Institute for Molecular Bioscience and School of Electrical & Information Eng Introduction to PLINK

  19. Introduction Tri-allelic alleles PLINK can represent tri-allelic alleles Only very limited ability to analyse them Same SNP may appear several times in the MAP or BIM file Usually filter out tri-allelic alleles Often an issue when merging data sets Scott Hazelhurst Sydney Brenner Institute for Molecular Bioscience and School of Electrical & Information Eng Introduction to PLINK

  20. Introduction Strandedness Different chips or experiments may record a SNP using different strand A C T G See when merge data — appears to be multi-allelic PLINK reports apparently multi-allelic SNPs You can flip them – create a new data set Try merge again – if really multi-allelic should work Filter out remaining May incorrectly flip a few Scott Hazelhurst Sydney Brenner Institute for Molecular Bioscience and School of Electrical & Information Eng Introduction to PLINK

Download Presentation
Download Policy: The content available on the website is offered to you 'AS IS' for your personal information and use only. It cannot be commercialized, licensed, or distributed on other websites without prior consent from the author. To download a presentation, simply click this link. If you encounter any difficulties during the download process, it's possible that the publisher has removed the file from their server.

Recommend

The PLINK example GWAS analysed by PLINK and Sib-pair David Duffy Genetic Epidemiology

The PLINK example GWAS analysed by PLINK and Sib-pair David Duffy Genetic Epidemiology

The PLINK example GWAS analysed by PLINK and Sib-pair David Duffy Genetic Epidemiology Laboratory Introduction Overview of development of Sib-pair PLINK v. Sib-pair Overview of Sib-pair An extensible platform for genetic data

687 views • 51 slides
Computer Practical Exercise using GCTA (with PLINK and R) Overview Purpose This exercise

Computer Practical Exercise using GCTA (with PLINK and R) Overview Purpose This exercise

Computer Practical Exercise using GCTA (with PLINK and R) Overview Purpose This exercise repeats the FaST-LMM analysis from the previous exercise using the program GCTA. In addition, we use GCTA to estimate the heritability accounted for by

325 views • 7 slides
Management of Risks Related to Long Waves Impact and Oil Spills European Joint Master in Water

Management of Risks Related to Long Waves Impact and Oil Spills European Joint Master in Water

Management of Risks Related to Long Waves Impact and Oil Spills European Joint Master in Water and Coastal Management EU-COMET2 Nikolay Plink ICM Department Russian State Hydrometeorological University St. Petersburg plink@rshu.ru

469 views • 34 slides
INTRODUCTION INTRODUCTION INTRODUCTION INTRODUCTION INTRODUCTION INTRODUCTION INTRODUCTION

INTRODUCTION INTRODUCTION INTRODUCTION INTRODUCTION INTRODUCTION INTRODUCTION INTRODUCTION

TRAINING COURSE ON SOCIAL PROTECTION & FORMALIZATION TRINIDAD AND TOBAGO MARCH 15, 2017 INTRODUCTION INTRODUCTION INTRODUCTION INTRODUCTION INTRODUCTION INTRODUCTION INTRODUCTION INTRODUCTION Design of the NIS Assistance from the

547 views • 23 slides
Differentiating tide-influenced and tide-dominated deltas: Devonian Baltic Basin Piret

Differentiating tide-influenced and tide-dominated deltas: Devonian Baltic Basin Piret

Differentiating tide-influenced and tide-dominated deltas: Devonian Baltic Basin Piret Plink-Bjorklund MAIN POINTS: - facies, architecture & stratigraphy of tide-influenced and dominated deltas - the nature of tidal signatures in both

624 views • 28 slides
INTRODUCTION Introduction 2/42 INTRODUCTION Alternations I am giving bread to the

INTRODUCTION Introduction 2/42 INTRODUCTION Alternations I am giving bread to the

Patterns of variation in the expression of case and agreement Andrs Brny Leiden University Centre for Linguistics 9 November 2019, BaSIS Workshop a.barany@hum.leidenuniv.nl INTRODUCTION Introduction 2/42 INTRODUCTION

1.37k views • 42 slides
Introduction Introduction Introduction Introduction Outline Motivation Failures

Introduction Introduction Introduction Introduction Outline Motivation Failures

Introduction Introduction Introduction Introduction Outline Motivation Failures Definition and concepts Process and product properties Principles SE Approaches Motivation Software and the economy The economies of

1.03k views • 49 slides
Introduction Introduction Introduction Nationwide Cause for Concern 1

Introduction Introduction Introduction Nationwide Cause for Concern 1

Introduction Introduction Introduction Nationwide Cause for Concern 1 https://www.cdc.gov/std/stats17/infographic.htm Introduction Epidemiology in Western Colorado CPHE Surveillance Data Report, May 2018 Chlamydia Epidemiology in Western

977 views • 39 slides
DAQ introduction Purpose of this talk : (1) Introduction for those who have not been in every

DAQ introduction Purpose of this talk : (1) Introduction for those who have not been in every

DAQ introduction Purpose of this talk : (1) Introduction for those who have not been in every meeting (2) Some ideas of diagrams for section 7.1 (Introduction) of TP. Giles Barr 25 Jan 2018, Pembroke College On DUNE, the reality is there is a

474 views • 10 slides
INTRODUCTION cf. Schneider, Chapter 1 INTRODUCTION THEY ARE OUT TO GET YOU INTRODUCTION WHAT

INTRODUCTION cf. Schneider, Chapter 1 INTRODUCTION THEY ARE OUT TO GET YOU INTRODUCTION WHAT

ADVANCED COMPUTER SECURITY INTRODUCTION cf. Schneider, Chapter 1 INTRODUCTION THEY ARE OUT TO GET YOU INTRODUCTION WHAT IS SECURITY? A systems security policies describe What the system is supposed to do Store and provide access

479 views • 16 slides
Overview Introduction to SMIL Introduction to W3C and XML Introduction to SMIL

Overview Introduction to SMIL Introduction to W3C and XML Introduction to SMIL

Overview Introduction to SMIL Introduction to W3C and XML Introduction to SMIL Writing a SMIL file Adding Media Objects Martin Halvey Clipping media files January 2008 Timing and Synchronising Layout World Wide

387 views • 14 slides
14 Introduction Introduction Bad guys can put malware into Example: hosts

14 Introduction Introduction Bad guys can put malware into Example: hosts

Introduction Introduction source Encapsulation ISO/OSI reference model message M application segment H t H t M transport network datagram H n H n H t M presentation: allow applications to frame link H l H n H t M interpret

583 views • 4 slides
Previous Lecture Introduction to SMIL 2 Introduction to W3C and XML Introduction to SMIL

Previous Lecture Introduction to SMIL 2 Introduction to W3C and XML Introduction to SMIL

Previous Lecture Introduction to SMIL 2 Introduction to W3C and XML Introduction to SMIL Writing a SMIL file Adding Media Objects Martin Halvey Clipping media files January 2008 Timing and Synchronising Layout

322 views • 8 slides
Lecture 1 Andreas Habegger Introduction Zynq Introduction Zynq Introduction Zynq PS vs. PL

Lecture 1 Andreas Habegger Introduction Zynq Introduction Zynq Introduction Zynq PS vs. PL

Introduction Zynq Lecture 1 Andreas Habegger Introduction Zynq Introduction Zynq Introduction Zynq PS vs. PL Processing System Processor Peripherals Data Buses AXI Bus Conclusion BTE5380 - Embedded Systems Octobre 2014 Andreas Habegger

1.81k views • 30 slides
2018.06 01 SMILE5 Introduction S E 5 02 Alpha Cloud M I L 03 Company Introduction 04

2018.06 01 SMILE5 Introduction S E 5 02 Alpha Cloud M I L 03 Company Introduction 04

2018.06 01 SMILE5 Introduction S E 5 02 Alpha Cloud M I L 03 Company Introduction 04 Project References S E 5 PART.1 SMILE5 Introduction M I L Introduction STORION-SMILE5 Residential Energy Storage System Inverter 5kW Output Input

1.06k views • 25 slides
Introduction to CICS Course introduction Course introduction What is CICS? What is an

Introduction to CICS Course introduction Course introduction What is CICS? What is an

Introduction to CICS Course introduction Course introduction What is CICS? What is an application server? Why use an application server? Course introduction Services provided by CICS How CICS applications are defined Scaling CICS to meet

1.84k views • 146 slides
PhaRmacodynamic Effects of Switching thErapy in PCI patients with high on Treatment platelet

PhaRmacodynamic Effects of Switching thErapy in PCI patients with high on Treatment platelet

PhaRmacodynamic Effects of Switching thErapy in PCI patients with high on Treatment platelet reactivity and genotype variation: high Clopidogrel dose versus Prasugrel (RESET GENE Study). (ClinicalTrials.gov NCT01465828) GENNARO SARDELLA MD,

385 views • 22 slides
Computations with Markers Paulino Prez 1 Jos Crossa 1 1 ColPos-Mxico 2 CIMMyT-Mxico June,

Computations with Markers Paulino Prez 1 Jos Crossa 1 1 ColPos-Mxico 2 CIMMyT-Mxico June,

Computations with Markers Paulino Prez 1 Jos Crossa 1 1 ColPos-Mxico 2 CIMMyT-Mxico June, 2015. CIMMYT, Mxico-SAGPDB Computations with Markers 1/20 Contents Genomic relationship matrix 1 Examples 2 Big Data! 3 CIMMYT,

491 views • 20 slides
Familial Pulmonary Fibrosis And Role of Genetic Testing Saturday, November 3, 2018 UCSF CME

Familial Pulmonary Fibrosis And Role of Genetic Testing Saturday, November 3, 2018 UCSF CME

Genetics of Pulmonary Fibrosis: Time for Genetic Screening? Familial Pulmonary Fibrosis And Role of Genetic Testing Saturday, November 3, 2018 UCSF CME Program: ILD The Changing Landscape San Francisco Marriott Union Square Christine Kim

492 views • 31 slides
Detection of copy number alterations and loss of heterozygosity Control-FREEC tutorial Valentina

Detection of copy number alterations and loss of heterozygosity Control-FREEC tutorial Valentina

Complete Genome analysis: Detection of copy number alterations and loss of heterozygosity Control-FREEC tutorial Valentina BOEVA Institut Curie, INSERM, Mines ParisTech Workshop outlines Motivation for copy number detection in cancer

1.32k views • 82 slides
Using RNA-Seq to reveal expression & diversity in maize

Using RNA-Seq to reveal expression & diversity in maize

Using RNA-Seq to reveal expression & diversity in maize C Robin Buell Michigan State University Department of Plant Biology Purdue 8 th InternaGonal

811 views • 35 slides
Todd A. Johnson RIKEN Center for Genomic Medicine Tokyo Medical & Dental University R User

Todd A. Johnson RIKEN Center for Genomic Medicine Tokyo Medical & Dental University R User

hzAnalyzer: Detection, quantification, and visualization of contiguous homozygosity in human populations from high-density genotyping datasets using R and Java Todd A. Johnson RIKEN Center for Genomic Medicine Tokyo Medical & Dental

464 views • 23 slides
Predic'ng 'ssue-specific effects of rare gene'c variants Farhan Damani Biological Data Sciences

Predic'ng 'ssue-specific effects of rare gene'c variants Farhan Damani Biological Data Sciences

Predic'ng 'ssue-specific effects of rare gene'c variants Farhan Damani Biological Data Sciences 2016 Goal: develop a framework to predict 'ssue- specific regulatory effects of rare variants Rare variants are abundant and poten'ally

659 views • 38 slides
management law Dr. Sndor Flp PhD, associate professor NUPS, Department of Water and

management law Dr. Sndor Flp PhD, associate professor NUPS, Department of Water and

Environmental law, water management law Dr. Sndor Flp PhD, associate professor NUPS, Department of Water and Environmental Policy 1 The system of environmental la law 1a. A (medium level) branch of law (4) note: number in bracket signs the

1.51k views • 101 slides

More recommend