-
Notifications
You must be signed in to change notification settings - Fork 9
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
Move rasi's last lecture after Gavin's
- Loading branch information
Showing
15 changed files
with
86 additions
and
86 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
File renamed without changes.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -1,4 +1,30 @@ | ||
# Lecture 14: Working with biological sequences and annotations using `Bioconductor` | ||
# Lecture 15: Introduction to Sequencing Data Analysis | ||
|
||
Now that we have a basic grasp of concepts surrounding data management, manipulation, and visualization, we're ready to start focusing on some of the more specialized data encountered in computational biology research. Sequencing of nucleic acids is almost ubiquitous in biological research. In this lecture, we will introduce some common resources for depositing and retrieving sequence data generated by consortium efforts and independent laboratories. We will introduce concepts and practical steps of querying, inspecting, and visualizing sequence data. Then, we will cover the types of genomic variation and common tools used to predict these from sequencing data. | ||
|
||
We will learn how to work with DNA sequences and transcript annotations. Open [lecture14.ipynb](./lecture14.ipynb) in VSCode. Make sure to select the `kernel` for `R` so that you can execute `R` code. You should have already set this up following the software installation instructions [here](../../software/README.md). | ||
This lecture focuses on concepts surrounding genome sequence data and their associated workflows. This lecture will include demonstrations and student exercises. We will dive into details of sequencing data and formats, as well as outputs for specific sequencing analysis commands. There will also be materials included as a resource for your future reference. | ||
|
||
## Learning objectives | ||
|
||
- Identify common databases and file formats used for sequence data | ||
- Describe the steps involved in processing and analyzing sequence data to predict different types of genomic variants | ||
- Recognize common tools (databases and software) used to assess variation in genomic data | ||
|
||
## Class materials | ||
|
||
Outline of [content from the slides](MCB536_lecture15_IntroSeqData.pdf): | ||
|
||
1. Sequence data | ||
- Databases and online resources for sequence data | ||
- Learn the common sequence data file formats | ||
|
||
2. Tools for sequencing data | ||
- Tools to query, inspect, visualize an aligned sequence file (demo + exercise) | ||
- Learn the contents of sequence data files (demo + exercise) | ||
- Learn to generate sequencing metrics and to process sequence data (demo + exercise) | ||
|
||
## Before the class | ||
|
||
1. Please be able to locate files for in-class exercises. Data and examples shown in `lecture 15` and [lecture 16](../lecture16/) are available on Fred Hutch filesystem at `/fh/fast/subramaniam_a/tfcb` and on [DropBox](https://www.dropbox.com/sh/zoitjnobgp7l7c2/AABBIpTQcNA4lWYOFnV5dlMKa?dl=0). For both `lecture 15` and `lecture 16`, you will need to download these files onto your laptop. | ||
|
||
2. Please install [Integrative Genomics Viewer (IGV)](https://software.broadinstitute.org/software/igv/). |
File renamed without changes.
File renamed without changes.
File renamed without changes.
File renamed without changes.
File renamed without changes.
File renamed without changes.
File renamed without changes.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -1,30 +1,65 @@ | ||
# Lecture 15: Introduction to Sequencing Data Analysis | ||
# Lecture 16: Genomic data in R | ||
|
||
Now that we have a basic grasp of concepts surrounding data management, manipulation, and visualization, we're ready to start focusing on some of the more specialized data encountered in computational biology research. Sequencing of nucleic acids is almost ubiquitous in biological research. In this lecture, we will introduce some common resources for depositing and retrieving sequence data generated by consortium efforts and independent laboratories. We will introduce concepts and practical steps of querying, inspecting, and visualizing sequence data. Then, we will cover the types of genomic variation and common tools used to predict these from sequencing data. | ||
This lecture will unite the last lecture's content on genomic analysis with our previous coding in R. The packages we'll use this week are from [Bioconductor](http://bioconductor.org), a collection of software specifically designed for genomic analysis in R. | ||
|
||
This lecture focuses on concepts surrounding genome sequence data and their associated workflows. This lecture will include demonstrations and student exercises. We will dive into details of sequencing data and formats, as well as outputs for specific sequencing analysis commands. There will also be materials included as a resource for your future reference. | ||
## Lecture Notes | ||
[Lecture slides](./MCB536_lecture16_GenomicDataInR.pdf) | ||
|
||
## Learning objectives | ||
|
||
- Identify common databases and file formats used for sequence data | ||
- Describe the steps involved in processing and analyzing sequence data to predict different types of genomic variants | ||
- Recognize common tools (databases and software) used to assess variation in genomic data | ||
Genome variant analysis (Background) | ||
- Types of genomic variation | ||
- Tools to predict genomic variations | ||
- Learn the common file formats for variation data | ||
- Databases and online resources for human variation data | ||
|
||
## Class materials | ||
Genomic Data (hands-on tutorials) | ||
- Use Bioconductor packages to work with genomic data in R | ||
- Load, inspect, and query genomic data (BED/SEG, BAM, VCF files) | ||
- Identify and annotate genomic variants | ||
|
||
Outline of [content from the slides](MCB536_lecture15_IntroSeqData.pdf): | ||
## Before the class | ||
|
||
1. Sequence data | ||
- Databases and online resources for sequence data | ||
- Learn the common sequence data file formats | ||
We will be working through some tutorials directly on your laptop using R Studio. | ||
|
||
2. Tools for sequencing data | ||
- Tools to query, inspect, visualize an aligned sequence file (demo + exercise) | ||
- Learn the contents of sequence data files (demo + exercise) | ||
- Learn to generate sequencing metrics and to process sequence data (demo + exercise) | ||
### 1. Install the R packages | ||
|
||
## Before the class | ||
- Tutorial is tested for R-4.0.3 | ||
- You should run [this script in VSCode](../../software/genomic_data.R) to ensure all Bioconductor packages are installed. | ||
``` | ||
## start R session ## | ||
R | ||
## run this command within R session ## | ||
source("../../software/genomic_data.R") | ||
``` | ||
- This script will install the following packages: | ||
- `Rsamtools`: querying BAM files | ||
- `VariantAnnotation`: reading VCF files | ||
- `GenomicRanges`: manipulating genomic data | ||
- `plyranges`: fast & easy tool for mannipulating GRanges | ||
### 2. Class materials: R Markdown files containing the tutorials | ||
- If you have not done so already, update your local copy of the class repository from GitHub. You should have a directory (`lecture16`) containing the following three RMarkdown tutorials: | ||
- [Lecture16_GenomicData.Rmd](Lecture16_GenomicData.Rmd): store genomic data as objects, assess genomic ranges, apply operations on genomic data | ||
- [Lecture16_Rsamtools.Rmd](Lecture16_Rsamtools.Rmd): load and query sequencing data; compute “pile-up” statistics at genomic loci to identify genomic variants | ||
- [Lecture16_VariantCalls.Rmd](Lecture16_VariantCalls.Rmd): load and assess variant (vcf) data | ||
### 3. Install R Extension in VSCode | ||
- `Extensions` (on left panel) > Type in search bar: `"R Extension"` > Select `R Extension for Visual Studio Code` by Yuki Ueda | ||
- The extension page should look something like this: https://marketplace.visualstudio.com/items?itemName=Ikuyadeu.r | ||
### 4. Install `pandoc` | ||
1. Please be able to locate files for in-class exercises. Data and examples shown in `lecture 15` and [lecture 16](../lecture16/) are available on Fred Hutch filesystem at `/fh/fast/subramaniam_a/tfcb` and on [DropBox](https://www.dropbox.com/sh/zoitjnobgp7l7c2/AABBIpTQcNA4lWYOFnV5dlMKa?dl=0). For both `lecture 15` and `lecture 16`, you will need to download these files onto your laptop. | ||
- To knit R Markdown files, you'll need the R Extension as well as `pandoc`. | ||
- Install `pandoc` outside of VScode by downloading the installer here: https://pandoc.org/installing.html | ||
### 5. Class materials: Genomic and sequencing data for the tutorials | ||
2. Please install [Integrative Genomics Viewer (IGV)](https://software.broadinstitute.org/software/igv/). | ||
- Please download all data files found in [this folder](https://www.dropbox.com/sh/zoitjnobgp7l7c2/AABBIpTQcNA4lWYOFnV5dlMKa?dl=0) and add them to your `lecture16` directory. The files should have the following filenames: | ||
- `BRCA.genome_wide_snp_6_broad_Level_3_scna.seg` | ||
- `BRCA_IDC_cfDNA.bam` | ||
- `BRCA_IDC_cfDNA.bam.bai` | ||
- `GIAB_highconf_v.3.3.2.vcf.gz` (if this file was automatically uncompressed on your computer, resulting in a file named `GIAB_highconf_v.3.3.2.vcf`, look in your Trash folder to find the original file ending in `gz`) | ||
- `GIAB_highconf_v.3.3.2.vcf.gz.tbi` |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -1,65 +1,4 @@ | ||
# Lecture 16: Genomic data in R | ||
# Lecture 14: Working with biological sequences and annotations using `Bioconductor` | ||
|
||
This lecture will unite the last lecture's content on genomic analysis with our previous coding in R. The packages we'll use this week are from [Bioconductor](http://bioconductor.org), a collection of software specifically designed for genomic analysis in R. | ||
|
||
## Lecture Notes | ||
[Lecture slides](./MCB536_lecture16_GenomicDataInR.pdf) | ||
|
||
## Learning objectives | ||
|
||
Genome variant analysis (Background) | ||
- Types of genomic variation | ||
- Tools to predict genomic variations | ||
- Learn the common file formats for variation data | ||
- Databases and online resources for human variation data | ||
|
||
Genomic Data (hands-on tutorials) | ||
- Use Bioconductor packages to work with genomic data in R | ||
- Load, inspect, and query genomic data (BED/SEG, BAM, VCF files) | ||
- Identify and annotate genomic variants | ||
|
||
## Before the class | ||
|
||
We will be working through some tutorials directly on your laptop using R Studio. | ||
|
||
### 1. Install the R packages | ||
|
||
- Tutorial is tested for R-4.0.3 | ||
- You should run [this script in VSCode](../../software/genomic_data.R) to ensure all Bioconductor packages are installed. | ||
``` | ||
## start R session ## | ||
R | ||
## run this command within R session ## | ||
source("../../software/genomic_data.R") | ||
``` | ||
- This script will install the following packages: | ||
- `Rsamtools`: querying BAM files | ||
- `VariantAnnotation`: reading VCF files | ||
- `GenomicRanges`: manipulating genomic data | ||
- `plyranges`: fast & easy tool for mannipulating GRanges | ||
### 2. Class materials: R Markdown files containing the tutorials | ||
- If you have not done so already, update your local copy of the class repository from GitHub. You should have a directory (`lecture16`) containing the following three RMarkdown tutorials: | ||
- [Lecture16_GenomicData.Rmd](Lecture16_GenomicData.Rmd): store genomic data as objects, assess genomic ranges, apply operations on genomic data | ||
- [Lecture16_Rsamtools.Rmd](Lecture16_Rsamtools.Rmd): load and query sequencing data; compute “pile-up” statistics at genomic loci to identify genomic variants | ||
- [Lecture16_VariantCalls.Rmd](Lecture16_VariantCalls.Rmd): load and assess variant (vcf) data | ||
### 3. Install R Extension in VSCode | ||
- `Extensions` (on left panel) > Type in search bar: `"R Extension"` > Select `R Extension for Visual Studio Code` by Yuki Ueda | ||
- The extension page should look something like this: https://marketplace.visualstudio.com/items?itemName=Ikuyadeu.r | ||
### 4. Install `pandoc` | ||
- To knit R Markdown files, you'll need the R Extension as well as `pandoc`. | ||
- Install `pandoc` outside of VScode by downloading the installer here: https://pandoc.org/installing.html | ||
### 5. Class materials: Genomic and sequencing data for the tutorials | ||
- Please download all data files found in [this folder](https://www.dropbox.com/sh/zoitjnobgp7l7c2/AABBIpTQcNA4lWYOFnV5dlMKa?dl=0) and add them to your `lecture16` directory. The files should have the following filenames: | ||
- `BRCA.genome_wide_snp_6_broad_Level_3_scna.seg` | ||
- `BRCA_IDC_cfDNA.bam` | ||
- `BRCA_IDC_cfDNA.bam.bai` | ||
- `GIAB_highconf_v.3.3.2.vcf.gz` (if this file was automatically uncompressed on your computer, resulting in a file named `GIAB_highconf_v.3.3.2.vcf`, look in your Trash folder to find the original file ending in `gz`) | ||
- `GIAB_highconf_v.3.3.2.vcf.gz.tbi` | ||
We will learn how to work with DNA sequences and transcript annotations. Open [lecture14.ipynb](./lecture14.ipynb) in VSCode. Make sure to select the `kernel` for `R` so that you can execute `R` code. You should have already set this up following the software installation instructions [here](../../software/README.md). |
File renamed without changes.
File renamed without changes.
File renamed without changes.