Sequencing bias: comparison of different protocols of microRNA library construction.
BACKGROUND
MicroRNAs(miRNAs) are 18-25 nt small RNAs playing critical roles in many biological processes. The majority of known miRNAs were discovered by conventional cloning and a Sanger sequencing approach. The next-generation sequencing (NGS) technologies enable in-depth characterization of the global repertoire of miRNAs, and different protocols for miRNA library construction have been developed. However, the possible bias between the relative expression levels and sequences introduced by different protocols of library preparation have rarely been explored.


RESULTS
We assessed three different miRNA library preparation protocols, SOLiD, Illumina versions 1 and 1.5, using cloning or SBS sequencing of total RNA samples extracted from skeletal muscles from Hu sheep and Dorper sheep, and then validated 9 miRNAs by qRT-PCR. Our results show that SBS sequencing data highly correlate with Illumina cloning data. The SOLiD data, when compared to Illumina's, indicate more dispersed distribution of length, higher frequency variation for nucleotides near the 3'- and 5'-ends, higher frequency occurrence for reads containing end secondary structure (ESS), and higher frequency for reads that do not map to known miRNAs. qRT-PCR results showed the best correlation with SOLiD cloning data. Fold difference of Hu sheep and Dorper sheep between qRT-PCR result and SBS sequencing data correlated well (r = 0.937), and fold difference of miR-1 and miR-206 among SOLiD cloning data, qRT-PCR and SBS sequencing data was similar.


CONCLUSIONS
The sequencing depth can influence the quantitative measurement of miRNA abundance, but the discrepancy caused by it was not statistically significant as high correlation was observed between Illumina cloning and SBS sequencing data. Bias of length distribution, sequence variation, and ESS was observed between data obtained with the different protocols. SOLiD cloning data differ from Illumina cloning data mainly because of distinct methods of adapter ligation. The good correlation between qRT-PCR result and SOLiD data might be due to the similarities of the hybridization-based methods. The fold difference analysis indicated that methods based on hybridization may be superior for quantitative measurement of miRNA abundance. Because of the genome sequence of the sheep is not available, our data may not explain how the entire miRNA bias in the natural miRNAs in sheep or other mammal miRNA expression, unbiased artificially synthesized miRNA will help on evaluating the methodology of miRNA library preparation.
Tian
Geng
G
Beijing Genomics Institute, Shenzhen 518000, China.
Yin
XuYang
X
Luo
Hong
H
Xu
XiaoHong
X
Bolund
Lars
L
Zhang
XiuQing
X
Gan
Shang-Quang
SQ
Li
Ning
N
D000818:Animals; D015723:Gene Library; D035683:MicroRNAs / Q000235:genetics / Q000378:metabolism*; D009690:Nucleic Acid Conformation; D020133:Reverse Transcriptase Polymerase Chain Reaction; D017423:Sequence Analysis, RNA / Q000379:methods*; D012756:Sheep
D035683:MicroRNAs
Methods. 2008 Jan;44(1):3-12
18158127
Methods. 2007 Oct;43(2):110-7
17889797
BMC Genomics. 2009 Feb 05;10:65
19196471
J Mol Biol. 1999 May 21;288(5):911-40
10329189
Haematologica. 2008 Jul;93(7):1009-16
18508790
Nucleic Acids Res. 2005 Nov 27;33(20):e179
16314309
Genome Res. 2004 Jun;14(6):1188-90
15173120
Nat Cell Biol. 2009 Mar;11(3):228-34
19255566
Curr Protoc Mol Biol. 2008 Oct;Chapter 3:Unit3.15
18972386
Nucleic Acids Res. 2003 Jul 1;31(13):3406-15
12824337
Neoplasia. 2009 Feb;11(2):167-76
19177201
Curr Protoc Mol Biol. 2005 Nov;Chapter 26:Unit 26.4
18265364
Exp Biol Med (Maywood). 2008 Feb;233(2):131-8
18222968
Cell. 2007 Jun 29;129(7):1401-14
17604727
Plant Cell. 2005 May;17(5):1397-411
15805478
Stem Cells. 2008 Oct;26(10):2496-505
18583537
RNA. 2003 Mar;9(3):277-9
12592000
Nucleic Acids Res. 2006 Jan 1;34(Database issue):D140-4
16381832
Nature. 2005 Mar 17;434(7031):338-45
15735639