Abstract
To conduct differential gene expression analysis, ovaries from the cattle tick Rhipicephalus microplus were dissected at three distinct developmental stages (preingurgitated, immature ingurgitated, and mature ingurgitated). Additionally, undissected intact mature males and complete ingurgitated female ticks without ovaries (carcasses) were also collected to serve as reference samples for analysis. To perform total RNA purification, tissue from ten individuals representing each of the five previously described conditions was pooled. mRNA was isolated from the purified total RNA using the oligo (dT) method. Following fragmentation, double stranded cDNA was synthesized and ligated to sequencing adapters. Suitable-sized fragments were subsequently used for PCR amplification. Libraries were analyzed and quantified using an Agilent 2100 Bioanalyzer and an ABI StepOnePlus Real-Time PCR System. A total of 45.64 Gb bases were sequenced using the Illumina HiSeq sequencing platform. After assembling the samples and correcting for abundance, we obtained 82,877 unigenes. The total length, average length, N50, and GC content of the unigenes were 89,754,828 bp,1,082 bp,2,068 bp and 49.04 % respectively. For functional annotation, the unigenes were aligned with 7 functional databases. The number of unigenes identified in the functional databases were as follows: 32,518 (NR:39.24 %), 10,259 (NT:12.38 %), 23,624 (Swissprot:28.50 %), 22,203 (KOG:26.79 %), 25,072 (KEGG:30.25 %), 17,435(GO:21.04 %), and 23,220 (InterPro:28.02 %). Unigene candidate coding regions (CDS) among the unigenes were predicted using TransDecoder software and 42,143 CDS were detected. We also detected 10,522 simple sequence repeats (SSRs) distributed on 8,126 unigenes, and predicted 4,672 transcription factors (TF) coding unigenes. Our data can be used to identify genes that are important for male and female tick and arachnid reproduction and tick general physiology.