5 %). In order to extract biological processes and molecular functions statistically over-represented in SO libraries, we performed a hyper-geometrical test between GO terms from the SO library and those from the AO library, which represents the natural physiological conditions. The p-values were then adjusted
using Bonferroni’s correction. In order to perform a functional enrichment analysis of the unigenes extracted from the SSH, we used the FatiGO web tool [39] against the SO library. With respect to the GO analysis, four different levels of description (3, 4, 6, and 9) were chosen for the biological processes. Quantitative expression by Real-Time RT-PCR Gene expression quantification was performed in whole animal, ovaries, and immune tissues Aurora Kinase inhibitor (hemocytes and hematopoietic organs pooled) Topoisomerase inhibitor of asymbiotic and symbiotic females. RNA extractions For the whole animal condition,
each individual was crushed with pestle and mortar in liquid nitrogen. Total RNA extraction was performed from about 30 mg of powder with TRIzol® reagent according to the manufacturer’s instructions (Invitrogen). For ovaries and immune tissues, total RNA extractions were performed from 25 and 50 females respectively with RNeasy Mini Kit according to the manufacturer’s instructions (QIAGEN). Real-Time RT-PCR First-strand cDNA was synthesized with the SuperScript III kit (Invitrogen) in accordance with manufacturer’s instructions, starting from 1 µg of total RNA using random hexamer primers. For whole animal samples, 0.2 µg of 5 individual extractions were pooled in 1 µg. Three biological replicates of each sample (whole animals, ovaries, and immune tissues) were used. For each gene, DNA Damage inhibitor primer pairs were designed with the Real-time PCR function of PerlPrimer [40]. The Tm and the length of each primer pair were fixed at 60°C and 18-22 bp, respectively.
Primers used for quantitative PCR are summarized in Additional File 1. Quantitative RT-PCR was performed using LightCycler LC480 system (Roche) as follows: 10 min at 95°C, 45 times [10 sec at 95°C, 10 sec at 60°C, 20 sec at 72°C]. A melting curve (65°C to 97°C) was recorded at the end of each reaction in order to check that the PCR product was unique. The reaction mixture consisted of 1.25 µL of each primer (10 µM), 5 µL of Fast SYBR-Green Master Mix (Roche) and 2.5 µL of diluted cDNA (corresponding to 12.5 ng of cDNA). Standard curves were plotted using 4 dilutions (125 ng, 25 ng, 5 ng, 1.25 ng) of pooled cDNAs from whole animals and ovaries. Efficiency of the PCR reaction was calculated. Expression data for each gene were estimated using the efficiency of the primer pair and the crossing point [41]. All gene expressions were normalized by the geometric mean of the expression level of the L8-ribosomal (RbL8) and Elongation Factor 2 (EF2) reference genes. Normalization and statistical pair-wise comparisons have been determined using REST [42].