Screening the human exome: a comparison of whole genome and whole transcriptome sequencing.
BACKGROUND
There is considerable interest in the development of methods to efficiently identify all coding variants present in large sample sets of humans. There are three approaches possible: whole-genome sequencing, whole-exome sequencing using exon capture methods, and RNA-Seq. While whole-genome sequencing is the most complete, it remains sufficiently expensive that cost effective alternatives are important.


RESULTS
Here we provide a systematic exploration of how well RNA-Seq can identify human coding variants by comparing variants identified through high coverage whole-genome sequencing to those identified by high coverage RNA-Seq in the same individual. This comparison allowed us to directly evaluate the sensitivity and specificity of RNA-Seq in identifying coding variants, and to evaluate how key parameters such as the degree of coverage and the expression levels of genes interact to influence performance. We find that although only 40% of exonic variants identified by whole genome sequencing were captured using RNA-Seq; this number rose to 81% when concentrating on genes known to be well-expressed in the source tissue. We also find that a high false positive rate can be problematic when working with RNA-Seq data, especially at higher levels of coverage.


CONCLUSIONS
We conclude that as long as a tissue relevant to the trait under study is available and suitable quality control screens are implemented, RNA-Seq is a fast and inexpensive alternative approach for finding coding variants in genes with sufficiently high expression levels.
Cirulli
Elizabeth T
ET
Center for Human Genome Variation, Duke University School of Medicine, Box 91009, Durham, NC 27708, USA. etc3@duke.edu
Singh
Abanish
A
Shianna
Kevin V
KV
Ge
Dongliang
D
Smith
Jason P
JP
Maia
Jessica M
JM
Heinzen
Erin L
EL
Goedert
James J
JJ
Goldstein
David B
DB
D001483:Base Sequence; D030541:Databases, Genetic; D005091:Exons / Q000235:genetics*; D020869:Gene Expression Profiling; D005786:Gene Expression Regulation; D015894:Genome, Human / Q000235:genetics*; D006801:Humans; D007963:Leukocytes, Mononuclear / Q000378:metabolism; D020641:Polymorphism, Single Nucleotide / Q000235:genetics; D016415:Sequence Alignment; D017422:Sequence Analysis, DNA / Q000379:methods*; D012689:Sequence Homology, Nucleic Acid