title: Sequencing bias: comparison of different protocols of microRNA library construction.
abstract: 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.
authors.affiliation: Beijing Genomics Institute, Shenzhen 518000, China.
authors.forename: Geng
authors.initials: G
authors.lastname: Tian
authors.forename: XuYang
authors.initials: X
authors.lastname: Yin
authors.forename: Hong
authors.initials: H
authors.lastname: Luo
authors.forename: XiaoHong
authors.initials: X
authors.lastname: Xu
authors.forename: Lars
authors.initials: L
authors.lastname: Bolund
authors.forename: XiuQing
authors.initials: X
authors.lastname: Zhang
authors.forename: Shang-Quang
authors.initials: SQ
authors.lastname: Gan
authors.forename: Ning
authors.initials: N
authors.lastname: Li
mesh_terms: 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
chemical_list: D035683:MicroRNAs
references.citation: Methods. 2008 Jan;44(1):3-12
references.pmid: 18158127
references.citation: Methods. 2007 Oct;43(2):110-7
references.pmid: 17889797
references.citation: BMC Genomics. 2009 Feb 05;10:65
references.pmid: 19196471
references.citation: J Mol Biol. 1999 May 21;288(5):911-40
references.pmid: 10329189
references.citation: Haematologica. 2008 Jul;93(7):1009-16
references.pmid: 18508790
references.citation: Nucleic Acids Res. 2005 Nov 27;33(20):e179
references.pmid: 16314309
references.citation: Genome Res. 2004 Jun;14(6):1188-90
references.pmid: 15173120
references.citation: Nat Cell Biol. 2009 Mar;11(3):228-34
references.pmid: 19255566
references.citation: Curr Protoc Mol Biol. 2008 Oct;Chapter 3:Unit3.15
references.pmid: 18972386
references.citation: Nucleic Acids Res. 2003 Jul 1;31(13):3406-15
references.pmid: 12824337
references.citation: Neoplasia. 2009 Feb;11(2):167-76
references.pmid: 19177201
references.citation: Curr Protoc Mol Biol. 2005 Nov;Chapter 26:Unit 26.4
references.pmid: 18265364
references.citation: Exp Biol Med (Maywood). 2008 Feb;233(2):131-8
references.pmid: 18222968
references.citation: Cell. 2007 Jun 29;129(7):1401-14
references.pmid: 17604727
references.citation: Plant Cell. 2005 May;17(5):1397-411
references.pmid: 15805478
references.citation: Stem Cells. 2008 Oct;26(10):2496-505
references.pmid: 18583537
references.citation: RNA. 2003 Mar;9(3):277-9
references.pmid: 12592000
references.citation: Nucleic Acids Res. 2006 Jan 1;34(Database issue):D140-4
references.pmid: 16381832
references.citation: Nature. 2005 Mar 17;434(7031):338-45
references.pmid: 15735639