Biol425 2013

From QiuLab
Revision as of 18:34, 30 March 2013 by imported>Weigang (→‎April 3: Transcriptome with R (Part 1))
Jump to navigation Jump to search
Computational Molecular Biology (BIOL 425/790.49, Spring 2013)
Instructors: Weigang Qiu (Associate Professor of Biology) & Che Martin (Assistant)
Room:1000G HN (10th Floor, North Building, Computer Science Department, Linux Lab FAQ)
Hours: Wednesdays 10:10 am-12:40 pm
Office Hours: Room 839 HN; Wednesdays 5-7pm or by appointment
Contacts: Dr Qiu: weigang@genectr.hunter.cuny.edu, 212-772-5296; Che: cmartin@gc.cuny.edu, 917-684-0864

General Information

Course Description

  • Background: Biomedical research is becoming a high-throughput science. As a result, information technology plays an increasingly important role in biomedical discovery. Bioinformatics is a new interdisciplinary field formed between molecular biology and computer science.
  • Contents: This course will introduce both bioinformatics theories and practices. Topics include: database searching, sequence alignment, molecular phylogenetics, structure prediction, and microarray analysis. The course is held in a UNIX-based instructional lab specifically configured for bioinformatics applications.
  • Problem-based Learning (PBL): For each session, students will work in groups to solve a set of bioinformatics problems. Instructor will serve as the facilitator rather than a lecturer. Evaluation of student performance include both active participation in the classroom work as well as quality of assignments (see #Grading Policy).
  • Learning Goals: After competing the course, students should be able to perform most common bioinformatics analysis in a biomedical research setting. Specifically, students will be able to
    • Approach biological questions evolutionarily ("Tree-thinking")
    • Evaluate and interpret computational results statistically ("Statistical-thinking")
    • Formulate informatics questions quantitatively and precisely ("Abstraction")
    • Design efficient procedures to solve problems ("Algorithm-thinking")
    • Manipulate high-volume textual data using UNIX tools, Perl/BioPerl, R, and Relational Database ("Data Visualization")
  • Pre-requisites: This 3-credit course is designed for upper-level undergraduates and graduate students. Prior experiences in the UNIX Operating System and at least one programming language are required. Hunter pre-requisites are CSCI132 (Practical Unix and Perl Programming) and BIOL300 (Biochemistry) or BIOL302 (Molecular Genetics), or permission by the instructor. Warning: This is a programming-based bioinformatics course. Working knowledge of UNIX and Perl is required for successful completion of the course.
  • Textbook: Krane & Raymer (2003). Fundamental Concepts of Bioinformatics. Pearson Education, Inc. (ISBN 0-8053-4633-3)
  • Academic Honesty: Hunter College regards acts of academic dishonesty (e.g., plagiarism, cheating on examinations, obtaining unfair advantage, and falsification of records and official documents) as serious offenses against the values of intellectual honesty. The College is committed to enforcing the CUNY Policy on Academic Integrity and will pursue cases of academic dishonesty according to the Hunter College Academic Integrity Procedures.

Grading Policy

  • Treat assignments as take-home exams. Student performance will be evaluated by weekly assignments and projects. While these are take-home projects and students are allowed to work in groups and answers to some of the questions are provided in the back of the textbook, students are expected to compose the final short answers, computer commands, and code independently. There are virtually an unlimited number of ways to solve a computational problem, as are ways and personal styles to implement an algorithm. Writings and blocks of codes that are virtually exact copies between individual students will be investigated as possible cases of plagiarism (e.g., copies from the Internet, text book, or each other). In such a case, the instructor will hold closed-door exams for involved individuals. Zero credits will be given to ALL involved individuals if the instructor considers there is enough evidence for plagiarism. To avoid being investigated for plagiarism, Do NOT copy from others or let others copy your work.
  • Submit assignments in Printed Hard Copies. Email attachments will NOT be accepted. Each assignment will be graded based on timeliness (10%), whether executable or having major errors (50%), algorithm efficiency (10%), and readability in programming styles (30%, see #Assignment Expectations).
  • The grading scheme
    • Assignments (100 pts): 10 exercises.
    • Mid-term (50 pts).
    • Final Project (50 pts)
    • Classroom performance (50 pts): Active engagement in classroom exercises and discussions
    • Attendance (50 pts): 1 unexcused absences = 40; 2 absences = 30; More than 2 = 0.

Assignment Expectations

  • Use a programming editor (e.g., vi or emacs) so you could have features like automatic syntax highlighting, indentation, and matching of quotes and parenthesis.
  • All PERL code must begin with "use strict; and use warnings;" statements. For each assignment, unless otherwise stated, I would like the full text of the source code. Since you cannot print using the text editor in the lab (even if you are connected from home), you must copy and paste the code into a word processor or a local text editor. If you are using a word processor, change the font to a fixed-width/monospace font. On Windows, this is usually Courier.
  • Also, unless otherwise stated, both the input and the output of the program must be submitted as well. This should also be in fixed-width font, and you should label it in such a way so that I know it is the program's input/output. This is so that I know that you've run the program, what data you have used, and what the program produced. If you are working from the lab, one option is to email the code to yourself, change the font, and then print it somewhere else as there is no printer in the lab.
  • Recommended Style
  • Bad Style

Course Schedule (All Wednesdays)

January 30. Browse Genome & Transcriptome Files with Unix Tools

  • Course Overview
  • Learning Goal: The power of Unix text filters
  • In-Class Exercises:
  1. You will need these 2 files for the following questions: /data/biocs/b/bio425/data/GBB.1con, GBB.seq
    1. What is a genome? What does a bacterial genome typically consist of? Explain the following terms: chromosome, plasmids, and contig
    2. Specify the FASTA file format
    3. What is the size of the Borrelia burgdorferi (the Lyme disease pathogen) B31 genome in terms of
      1. Number of replicons: From your home directory, run:
        grep -c "^>" ../../bio425/data/GBB.1con [Answer: N=22 replicons] 
      2. Number of genes: From your home directory, run:
        grep -c "^>" ../../bio425/data/GBB.seq [Answer: N=1,738 genes] 
      3. Number of bases: First filter out FASTA headers using "grep -v" and then remove newline characters using "tr -d":
        grep -v "^>" ../../bio425/data/GBB.1con | tr -d '\n' | wc -m [Answer: N=1,519,856 bases] 
  2. Based on the file /data/biocs/b/bio425/data/ge-breast-cancer-cell-lines.dat, answer the following questions:
    1. What is a transcriptome?
    2. How many unique genes are represented in this chip?
      cut -f1 ../../bio425/data/ge-breast-cancer-cell-lines.dat | grep -vc "^Description" [Answer: N=18,900 genes]
    3. How many cell lines in the file?
      grep "^Description" ../../bio425/data/ge-breast-cancer-cell-lines.dat| cut -f2- | wc -w [Answer: N=59 cells]
    4. Extract gene expression values at 3 breast cancer clinical markers: ERBB2, ESR1, and PGR
      grep -wP "ERBB2|ESR1|PGR" ../../bio425/data/ge-breast-cancer-cell-lines.dat [You should see 3 rows of gene expression values]
Assignment #1 (Final Version; Q4 file name corrected)
Unix Text Filters
  1. Display the absolute path of your home directory
  2. List files in your home directory in long format & ordered by their time stamps
  3. List files and directories in the "/data/biocs/b/bio425/" directory from your home directory
  4. Count the number of plasmids in the B. burgdorferi genome using the file "/data/biocs/b/bio425/GBB.1con"
  5. Show the first five lines of the file "GBB.seq" & save it to a file with arbitrary name
  6. Show your last ten commands using "history"
Read Chapter 1

February 6. Central Dogma with PERL (Part 1)

  • Learning goals:
  1. DNA Replication and Transcription
  2. String manipulations using BASH & PERL
  • In-Class Exercises:
  1. Find out base composition (%A, %T, %C, and %G) of a DNA sequence
    1. In the folder /data/biocs/b/bio425/data/GBB.1con-splitted, the whole B. burgdorferi genome has been splitted into files each of which containing a single replicon. Write a short BASH script to count how many bases in each replicon.
      for f in ../../bio425/data/GBB.1con-splitted/*.fas; do echo -ne "$f\t"; grep -v "^>" $f | tr -d "\n" | wc -m; done
    2. Write a PERL script to count GC% given a FASTA sequence.
      Script posted at /data/biocs/b/bio425/scripts/gc-counter.pl
    3. Run the GC-counter for each replicon using a BASH loop.
      for f in ../../bio425/data/GBB.1con-splitted/*.fas; do perl /data/biocs/b/bio425/scripts/gc-counter.pl $f; done
  2. Central Dogma
    1. Illustrate the Central Dogma using lines representing DNA/RNA/Protein. Indicate directionality of all parts.
      DNA replication, Transcription & Translation: all 5' to 3'
    2. Define cDNA. What are some of the common uses of cDNA?
      Identify alternative-splicing variants (isoforms); Quantify gene expression levels (RNA-SEQ)
  3. DNA replication
    1. What is the implied directionality of a DNA sequence if none is specified?
      5' to 3'
    2. What is the complement strand of this DNA sequence gatactaatgaagtat in 5' to 3' direction?
      atacttcattagtatc
    3. Write a PERL script to output the complementary strand of the ospC gene /data/biocs/b/bio425/data/ospC.seq in 5' to 3' direction. Save the result in a FASTA file ospC.revcom.fas
      Script posted at /data/biocs/b/bio425/scripts/revcom.pl
Assignment #2 (Final)
  1. Modify the gc-counter code /data/biocs/b/bio425/scripts/gc-counter.pl to (1) count only a,t,c,g's (not other letters, e.g., amino acids); see below (2) count bases in both upper-case AND lower-case letters; if ($base =~ /[ATCG]/i) { $base =~ tr/atcg/ATCG/; $count{$base}++; $total_valid++} (3) output GC% (i.e., percentage of G+C among the total bases; formated to no more than 2 decimal points, by using the printf function). my $gc_percent=($count{"G"}+$count{"C"})/$total_valid; printf "%.2f\n", $gc_percent;Show results for the file /data/biocs/b/bio425/data/cp9.seq. 23.69%
  2. Add comments to the revcom script /data/biocs/b/bio425/scripts/revcom.pl to explain what each statement does.
  3. Chapter 1. Problems (pg.31-32):
    1. 1.4. Would you expect the amino acids containing hydroxyl groups (-OH) hydrophobic or hydrophilic? Why? Hydrophilic, because it forms hydrogen bond with the water molecule
    2. 1.6. How does a cDNA library differ from a genomic library? cDNA library does not contain intronic sequences while genomic DNA does.

February 13. Central Dogma with PERL (Part 2)

  • Learning goals: DNA translation & ORF Identification
  • In-Class Exercises:
  1. Six-frame translation
    1. Define ORF and reading frame. How many possible reading frames are there for a given sequence?
      ORF: Hypothetical, computer-predicted protein-coding sequences. An ORF needs experimental data (e.g., mRNA) to be verified as a true gene. There are six possible reading frames for a given DNA sequence: three (+1, +2, +3) on one strand and another three (-1, -2, -3) on the reverse-complement strand. 
    2. How does a cell normally choose only one reading frame out of other possibilities for any given parts of a genome?
      In a prokaryotic cell, transcription factors (TFs) bind to a particular sequence (called "promoters") at the 5'-upstream of an ORF.  TFs recruits RNA polymerase to initiative the transcription of a gene. After the transcription, ribosomes bind to another particular 5'-upstream sequence (SD sequence, normally "AGGAG") to start the translation of a protein.
    3. Using a DNA genetic code table to manually translate the following sequence in all possible reading frames: gatactaatgaagtat . Which reading frame do you think is the most likely one?
      +1: DTNEV; +2: ILMKY; +3:Y**S; -1:ILH*Y; -2:YFISI; -3:TSLV. Frames +3 and -1 could be eliminated.
  2. ORF identification in a bacterial genome
    1. Pick a plasmid sequence from /data/biocs/b/bio425/data/GBB.1con-splitted
      Run ls -lh ../../bio425/data/GBB.1con-splitted/ and choose a small plasmid file, e.g., Borrelia_burgdorferi_4041_cp9_plasmid_C.fas
    2. Identify ORFs on the plasmid by running the program /data/biocs/b/bio425/bin/long-orfs
      /data/biocs/b/bio425/bin/long-orfs /data/biocs/b/bio425/data/GBB.1con-splitted/Borrelia_burgdorferi_4041_cp9_plasmid_C.fas cp9.coord
    3. Extract ORF sequences using the program /data/biocs/b/bio425/bin/extract
      /data/biocs/b/bio425/bin/extract /data/biocs/b/bio425/data/GBB.1con-splitted/Borrelia_burgdorferi_4041_cp9_plasmid_C.fas cp9.coord
Assignment #3 (Final Version)
  1. (10 pts) Write a PERL program to perform 6-frame translation given a FASTA sequence (output: a FASTA file with 6 protein sequences; single-letter amino acid code; * for stop codon). Code posted as /data/biocs/b/bio425/scripts/6-frame.pl. Show results for the input file /data/biocs/b/bio425/data/mystery_seq1.fas. How many long (>=200 bases) ORFs do you see? What are their reading frames? Two long ORFs: one in the "+1" frame (L=107 aa) and another in the "-3" frame (L=177 aa).
  2. (5 pts) Run the "long-orf" program on the above mystery sequence using the following 2 options: "ATG" as start codon and a minimum ORF length of 200 bases. Save the coordinate file as "mystery_seq1.coord". ../../bio425/bin/long-orfs -A ATG -g 200 ../../bio425/data/mystery_seq1.fas mystery_seq1.coord
  3. (10 pts) Write a PERL program that does exactly the same as the program /data/biocs/b/bio425/bin/extract. Input files: (1) "mystery_seq1.fas" (2) "mystery_seq1.coord". Output: A single FASTA file with multiple ORF sequences. Code posted as /data/biocs/b/bio425/scripts/extract.pl.

February 20 (No Class)

  • Monday Schedule

February 27. Genomics with BioPerl (Part 1)

  • Learning goals: Introduction to Objects
  • In-Class Exercises:
  1. Construct a sequence object with a hash
    1. Select and write down at least five features (e.g., strain, gene) associated with a DNA sequence based on this GenBank record. Include a sequence label and the sequence itself as two of the features.
    2. Write a PERL program named "seq_object.pl", which captures these five features with a hash by using feature names as keys and features as values. Example: my %seq_obj=("strain"=>"SD91, "gene"=>"ospC"); Print the hash with Data::Dumper.
      my $seq_hash=("id"=>"SD1_ospC", "seq"=>"tgtaataattcaggaaaaga", "organism"=>"Borrelia burgdorferi", "strain"=>"SD1", "gene"=>"ospC"); print Dumper(\%seq_hash)
    3. Create a scalar variable called "$ref_seq_obj", which is a reference to the hash %seq_obj. Print the hash reference with Data::Dumper.
      my $ref_seq_obj=\%seq_obj; print Dumper($ref_seq_obj)
    4. Name some of the differences between a sequence and a sequence object.
      Sequence refers to a string of nucleotides or amino acids. A sequence is usually represented by a string variable. A sequence object includes biological features associated with a sequence besides the sequence itself. A sequence object is commonly represented by a hash reference in Perl.
  2. Construct a sequence object with Bio::Seq
    1. Modify your program by calling the Bioperl Bio::Seq module:
      use lib '/data/biocs/b/bio425/bioperl-live'; use Bio::Seq;
    2. Convert your hash-based sequence object into a Bio::Seq object:
      my $seqobj = Bio::Seq->new( -display_id => "SD1_ospC", -seq =>"tgtaataattcaggaaaaga" );
    3. Print the Bio::Seq object with Data::Dumper
      print Dumper($seqobj);
  3. Apply Bio::Seq methods
    1. Find the names of methods for the following operations on the Bio::Seq object: reverse complement, length, id, seq, translate, sub-sequence by running
      perldoc /data/biocs/b/bio425/bioperl-live/Bio/Seq.pm
    2. Add the above methods into your script & print all results
      my $seq_rev=$seqobj->revcom()->seq(); my $eq_length=$seqobj->length(); my $seq_id=$seqobj->display_id(); my $seq_string=$seq_obj->seq(); my $seq_translate=$seqobj->translate()->seq(); my $subseq = $seqobj->subseq(1,10)
  4. Discussion Questions
    1. What is Object-oriented Programming?
    2. What is an Object and what is a Method?
      A Perl object is a "blessed" (strictly defined) reference to a hash. A method is a subroutine or function associated with an object. An object and its associated methods together make a class (e.g., the Bio::Seq class).
    3. What are the advantages and disadvantages of Object-Oriented programming in comparison with traditional programming based on subroutines and functions?
      Advantages: Simpler and more readable code. Easier maintenance. Disadvantages: Code dependent on installation of external libraries (e.g., BioPerl). An extra layer of abstraction, e.g.,  the concept of objects and methods.
Assignment #4 (Final version)
  • (10 pts) Write a Bio::Seq-based PERL program to perform 6-frame translation of this sequence "/data/biocs/b/bio425/data/mystery_seq1.fas". [Hints: Reuse the file-reading code from your original script (or copy from the 6-frame.pl code). You need to construct two Bio::Seq objects: one for the original sequence, the other for its reverse complement. You may either use the following methods: revcom(), trunc(), seq(), translate(), or (a simpler but more advanced way) call the translate() method with a "-frame=>" argument. For each method, read the documents on usage. Pay particular attention to the input object (which has to be a Bio::Seq object here), arguments (if any), and the return value (whether a number, a string, or a reference). Code posted as /data/biocs/b/bio425/scripts/6-frame-bioperl.pl.

March 6: Genomics with BioPerl (Part 2)

  • Learning goal: File I/O with BioPerl
  • In-Class Exercises:
  1. Read a FASTA file with Bio::SeqIO
    1. Find out how to read a FASTA file using Bio::SeqIO by running
      perldoc /data/biocs/b/bio425/bioperl-live/Bio/SeqIO.pm
    2. Write a PERL script fasta-io.pl, which uses Bio::SeqIO to read a multi-FASTA file (input) and print the name and length of each sequence (output). Test your script on this file ../../bio425/data/Borrelia_osp.dna.fasta
  2. Output a FASTA file with Bio::SeqIO
    1. Find out how to output a FASTA file using Bio::SeqIO by (again) running
      perldoc /data/biocs/b/bio425/bioperl-live/Bio/SeqIO.pm
    2. Modify the PERL script fasta-io.pl to read a multi-FASTA file (input) and print a translation of each sequence (output). Test your script on this file ../../bio425/data/Borrelia_osp.dna.fasta
  3. Start working on assignment 5 (see below) if you have time.
Assignment #5 (Final Version)
  • (10 pts) Write a BioPerl version of extract.pl by using the Bio::SeqIO (for both reading and outputing FASTA files) and Bio::Seq modules. Input files: (1) "mystery_seq1.fas" (2) "mystery_seq1.coord". Output: A single FASTA file with translated protein sequences. [Hints: Use regular "open" function to read the COORD file (you may copy code from "/data/biocs/b/bio425/scripts/extract.pl" for this part).] Code posted as /data/biocs/b/bio425/scripts/extract-bioperl-version.pl

March 13: Mid-Term Review

  • In-class Exercise: Write a BASH script named as "get-long-orfs.bash" to output long ORFs and their translated protein sequences given a genome sequence. Your script will combine commands and scripts you have previously composed into a single-command utility.
    • Input file: Locate whole plasmid sequences in the "/data/biocs/b/bio425/data/GBB.1con-splitted/" folder:
      • Group 1: use "Borrelia_burgdorferi_4075_lp54_plasmid_A.fas"
      • Group 2: use "Borrelia_burgdorferi_4091_cp26_plasmid_B.fas"
      • Group 3: use "Borrelia_burgdorferi_4041_cp9_plasmid_C.fas"
      • Group 4: use "Borrelia_burgdorferi_4076_lp17_plasmid_D.fas"
      • Group 5: use "Borrelia_burgdorferi_4005_lp25_plasmid_E.fas"
    • Output files: two files (2 bonus points if your output file names are generated dynamically, not hard-coded)
      • "plasmid_X.nuc": ORF sequences in FASTA
      • "plasmid_X.pep": Protein sequences in FASTA
filename=$1; base=$(basename $filename ".fas"); coordfile=$base.coord; nucfile=$base.nuc; pepfile=$base.pep; /data/biocs/b/bio425/bin/long-orfs -A ATG -g 200 $filename $coordfile; perl extract.pl $filename $coordfile > $nucfile; perl extract-bioperl-version.pl $filename $coordfile > $pepfile
  • Sample Questions:
  1. Define elements of a prokaryotic gene: ORF, promoter, RBS-binding sites
  2. Construct a single "for" loop to translate a sequence in six frames
  3. Print the reverse-complement sequence given a Bio::Seq object
  • Mid-term will cover the following contents (Review all in-class exercises and assignments)
    • Biology: genome, transcriptome, 6-frame translation, gene structure
    • Linux file system & text tools: ls [-lart]; less; grep [-ivcwP]; cat; cut [-fd]; head; tail; bash script (for loop); tr; wc [-lwm]
    • Gene-finding using GLIMMER: long-orfs, extract
    • Basic Perl: file I/O; loops (for; foreach; while); regex; using hash for counts and translation
    • Object-oriented BioPerl: the use Bio::SeqIO & Bio::Seq

March 20. Midterm Practicum


March 27 (No Class)

Spring Break


April 3: Transcriptome with R (Part 1)

  • Learning goal: Introduction to R
  • R resources
  1. Forget about Excel! Instead, use R for statistical analysis and visualization of large data sets.
  2. Download and install R to your own computer from: R website
  3. Download and install R studio to your own computer from: R Studio website
  4. A nice tutorial: SimpleR
  • In-Class Exercises: Part 1.
  1. Set a working directory
    1. Using a Linux terminal, make a directory called R-work
    2. Start R-Studio and find working directory with getwd()
    3. Change working directory with setwd("~/R-work")
  2. The basic R data structure: Vector
    1. Construct a character vector with c.vect=c("red", "blue", "green", "orange", "yellow", "cyan")
    2. Construct a numerical vector with n.vect=c(2, 3, 6, 10, 1, 0, 1, 1, 10, 30)
    3. Construct vectors using functions
      1. n.seq1=1:20
      2. n.seq2=seq(from=1, to=20, by=2)
      3. n.rep1=rep(1:4, times=2)
      4. n.rep2=rep(1:4, times=2, each=2)
    4. Use built-in help of R functions: ?seq or help(seq)
    5. Find out data class using the class function: class(c.vect)
  3. Access vector elements
    1. A single value: n.vect[2]
    2. A subset: n.vect[3:6]
    3. An arbitrary subset: n.vect[c(3,7)]
    4. Exclusion: n.vect[-8]
    5. Conditional: n.vect[n.vect>5]
  4. Apply functions
    1. Length: length(n.vect)
    2. Range: range(n.vect); min(n.vect); max(n.vect)
    3. Descriptive statistics: sum(n.vect); var(n.vect); sd(n.vect)
    4. Sort: order(n.vect). How would you get an ordered vector of n.vect? [Hint: use ?order to find the return values]
    5. Arithmetics: n.vect +10; n.vect *10; n.vect / 10
  5. Quit an R session
    1. Click on the "History" tab to review all your commands
    2. Save your commands into a file by opening a new "R Script" and save it as vector.R
    3. Save all your data to a default file .RData and all your commands to a default file ".Rhistory" with save.image()
    4. Quit R-Studio with q()

  • In-Class Exercises: Part 2
  1. Start a new R studio session and set working directory to ~/R-work
  2. Generate a vector of 100 normal deviates using x.ran=rnorm(100)
  3. Visualize the data distribution using hist(x.ran)
  4. Explore help file using ?help; example(help)
  5. How to break into smaller bins? How to add color? How to change x- and y-axis labels? Change the main title?
  6. Test if the data points are normally distributed.
  7. Save a copy of your plot using "Export"->"Save Plot as PDF"

  • In-Class Exercises: Part 3
  1. Using a Linux terminal, make a copy of breast-cancer transcriptome file /data/biocs/b/bio425/data/ge.dat in your R working directory
  2. Start a new R studio session and set working directory to ~/R-work
  3. Read the transcriptome file using ge=read.table("ge.dat", header=T, row.names=1, sep="\t")
    1. What is the data class of ge?
    2. Access data frame. What do the following commands return? ge[1,]; ge[1:5,]; ge[,1]; ge[,1:5]; ge[1:5, 1:5]
    3. Apply functions: head(ge); tail(ge); dim(ge)
  4. Save a copy of an object: ge.df=ge.
  5. Transform the transcriptome into a numerical matrix for subsequent operations: ge=as.matrix(ge)
  6. Test if the expression values are normally distributed.
  7. Save and quit the R session

Assignment #6 (Not Final)
  • (0 pts) Install R and R Studio to your home/laptop computer
  • (5 pts) Make a histogram of the "ge.dat". Are the data points normally distributed? Show all R commands and the final plot.
  • (5 pts) Formulate and describe three questions the ge transcriptome data set may help answer for the research and treatment of breast cancers. [Examples (feel free to re-use and refine): Are all breast cancer cells similar or different in gene-expression levels? What genes are associated with early or late clinical stages? What genes are associated with drug/treatment resistance? Could gene expression be used for early diagnosis? For predicting drug responses (as in personalized medicine)?

April 10: Transcriptome with R (Part 2)

  • Learning goal: Introduction to R
  • In-Class Exercises:
  1. R-Studio
    1. TBA
    2. TBA
  2. The c() constructor
    1. TBA
    2. TBA
  3. R plots
Assignment #6 (TBA)
  • (5 pts) Guided Reading. Read this paper and answer the following questions
  1. Describe the main purpose of the study
  2. List at least 3 genomic features of cancer cells the study generated
  3. What microarray platform did the study use? How was the transcriptome data set used for predicting drug effects?

April 17: Transcriptome with R (Part 3)

  • Learning goal: Introduction to R
  • In-Class Exercises:
  1. R-Studio
    1. TBA
    2. TBA
  2. The c() constructor
    1. TBA
    2. TBA
  3. R plots
Assignment #6 (TBA)
  • (10 pts) TBA

April 24: Molecular Phylogenetics (Part 1)

  • Learning goal: Introduction to R
  • In-Class Exercises:
  1. R-Studio
    1. TBA
    2. TBA
  2. The c() constructor
    1. TBA
    2. TBA
  3. R plots
Assignment #6 (TBA)
  • (10 pts) TBA

May 1: Molecular Phylogenetics (Part 2)

  • Learning goal: Introduction to R
  • In-Class Exercises:
  1. R-Studio
    1. TBA
    2. TBA
  2. The c() constructor
    1. TBA
    2. TBA
  3. R plots
Assignment #6 (TBA)
  • (10 pts) TBA

May 8: Genome Database with SQL (Part 1)

  • Learning goal: Introduction to R
  • In-Class Exercises:
  1. R-Studio
    1. TBA
    2. TBA
  2. The c() constructor
    1. TBA
    2. TBA
  3. R plots
Assignment #6 (TBA)
  • (10 pts) TBA

May 15: Genome Database with SQL (Part 2)

  • Learning goal: Introduction to R
  • In-Class Exercises:
  1. R-Studio
    1. TBA
    2. TBA
  2. The c() constructor
    1. TBA
    2. TBA
  3. R plots
Assignment #6 (TBA)
  • (10 pts) TBA

May 22: Final Project Due

Useful Links

Unix Tutorials

Perl Help

  • Professor Stewart Weiss has taught CSCI132, a UNIX and Perl class. His slides go into much greater detail and are an invaluable resource. They can be found on his course page here.
  • Perl documentation at perldoc.perl.org. Besides that, running the perldoc command before either a function (with the -f option ie, perldoc -f substr) or a perl module (ie, perldoc Bio::Seq) can get you similar results without having to leave the terminal.

Bioperl

SQL

R Project

  • Install location and instructions for Windows
  • Install location and instructions for Mac OS X
  • Install R-Studio
  • For users of Ubuntu/Debian:
sudo apt-get install r-base-core

Utilities

Other Resources

© Weigang Qiu, Hunter College, Last Update ~~