期刊
BMC GENOMICS
卷 14, 期 -, 页码 -出版社
BMC
DOI: 10.1186/1471-2164-14-536
关键词
Next-generation sequencing; Mapping bias; Drosophila melanogaster; Drosophila simulans; DGRP; Allelic imbalance; Genomics; Gene expression; Illumina
资金
- National Science Foundation [MCB-1021398, 0903629]
- National Institutes of Health [1F32-GM089009-01A1]
- Direct For Biological Sciences [1021398] Funding Source: National Science Foundation
- Div Of Molecular and Cellular Bioscience [1021398] Funding Source: National Science Foundation
Background: RNA-seq can be used to measure allele-specific expression (ASE) by assigning sequence reads to individual alleles; however, relative ASE is systematically biased when sequence reads are aligned to a single reference genome. Aligning sequence reads to both parental genomes can eliminate this bias, but this approach is not always practical, especially for non-model organisms. To improve accuracy of ASE measured using a single reference genome, we identified properties of differentiating sites responsible for biased measures of relative ASE. Results: We found that clusters of differentiating sites prevented sequence reads from an alternate allele from aligning to the reference genome, causing a bias in relative ASE favoring the reference allele. This bias increased with greater sequence divergence between alleles. Increasing the number of mismatches allowed when aligning sequence reads to the reference genome and restricting analysis to genomic regions with fewer differentiating sites than the number of mismatches allowed almost completely eliminated this systematic bias. Accuracy of allelic abundance was increased further by excluding differentiating sites within sequence reads that could not be aligned uniquely within the genome (imperfect mappability) and reads that overlapped one or more insertions or deletions (indels) between alleles. Conclusions: After aligning sequence reads to a single reference genome, excluding differentiating sites with at least as many neighboring differentiating sites as the number of mismatches allowed, imperfect mappability, and/or an indel(s) nearby resulted in measures of allelic abundance comparable to those derived from aligning sequence reads to both parental genomes.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据