- Research article
- Open Access
A comprehensive microRNA expression profile of the backfat tissue from castrated and intact full-sib pair male pigs
BMC Genomics volume 15, Article number: 47 (2014)
It is widely known that castration has a significant effect on the accumulation of adipose tissue. microRNAs (miRNAs) are known to be involved in fat deposition and to be regulated by the androgen-induced androgen receptor (AR). However, there is little understanding of the relationship between miRNAs and fat deposition after castration. In this study, the high-throughput SOLiD sequencing approach was used to identify and characterize miRNA expression in backfat from intact and castrated full-sib male 23-week-old pigs. The patterns of adipogenesis and fat deposition were compared between castrated and intact male pigs.
A total of 366 unique miRNA genes were identified, comprising 174 known pre-miRNAs and 192 novel pre-miRNAs. One hundred and sixty-seven pre-miRNAs were common to both castrated (F3) and intact (F4) male pig small RNA libraries. The novel pre-miRNAs encoded 153 miRNAs/miRNA*s and 141 miRNAs/miRNA*s in the F3 and F4 libraries, respectively. One hundred and seventy-seven miRNAs, including 45 up- and 132 down-regulated, had more than 2-fold differential expression between the castrated and intact male pigs (p-value < 0.001). Thirty-five miRNAs were further selected, based on the expression abundance and differentiation between the two libraries, to predict their targets in KEGG pathways. KEGG pathway analyses suggested that miRNAs differentially expressed between the castrated and intact male pigs are involved in proliferation, apoptosis, differentiation, migration, adipose tissue development and other important biological processes. The expression patterns of eight arbitrarily selected miRNAs were validated by stem-loop reverse-transcription quantitative polymerase chain reaction. These data confirmed the expression tendency observed with SOLiD sequencing. miRNA isomiRs and mirtrons were also investigated in this study. Mirtrons are a recently described category of miRNA relying on splicing rather than processing by the microprocessor complex to generate the RNAi pathway. The functions of miRNAs important for regulating fat deposition were also investigated in this study.
This study expands the number of fat-deposition-related miRNAs in pig. The results also indicate that castration can significantly affect the expression patterns of fat-related miRNAs. The differentially expressed miRNAs may play important roles in fat deposition after castration.
microRNAs (miRNAs) are small regulatory molecules of between 18 to 24 nucleotides. They play important roles in many biological processes, including development, differentiation and metabolism [1, 2]. The miRNA regulation of target genes occurs mostly at the 3′-untranslated region of mRNAs. The transcriptional profiles of mouse 3T3-L1 cells (pre-adipocyte cells) have been monitored during differentiation using microarrays and of the 395 differentially expressed sequence tags, 71% of the differentially expressed genes could be regulated by miRNAs [3, 4]. The process of mature miRNA production contains three steps. First, miRNA genes yield pri-miRNAs via RNA polymerase II. Second, the pri-miRNAs are processed to generate pre-miRNAs by a complex containing Drosha. Finally, miRNA:miRNA* duplexes are generated under the effect of Dicer. The miRNA:miRNA* duplexes have two strands. One strand is more stable and is incorporated into an RNA-induced silencing complex (RISC), while the other strand (miRNA*) is usually discarded . However, Guo et al. reported two different fates for the miRNA* strand: it is either degraded as merely a carrier strand or it is expressed abundantly as a potential functional guide miRNA .
The involvement of miRNAs with lipid metabolism was first reported in 2003; a lack of miR-14 can increase the accumulation of diacylglycerol and triglyceride in Drosophila . In vivo and in vitro data show that miRNAs play an important role in lipogenesis [8–10]. For instance, the expression of miR-143 increased in differentiating adipocytes. Also the inhibition of miR-143 can effectively suppress the differentiation of adipocytes. This indicated a potential role of miR-143 in adipocyte differentiation . During differentiation of the mouse pre-adipocyte 3T3-L1 cell line, miR-17-92 was identified as a significantly up-regulated miRNA cluster at the clonal expansion stage . On the other hand, some miRNAs were also reported as inhibitors of adipocyte differentiation, such as miR-27b and let-7 [12, 13].
It is well known that castration of male pigs can decrease an unpleasant odor in pork. However, it also results in unwanted fat deposition [14, 15]. Mangelsdorf et al.  found that androgens influence gene transcription through the activation of the androgen receptor (AR), a ligand-activated transcription factor that binds specific DNA motifs in its target genes. Ribas et al. showed that androgen-induced AR can bind to the defined miR-21 promoter, miPPR-21, which demonstrated that AR can directly regulate the transcription of miR-21 mRNA [17, 18]. We hypothesized that decreased levels of AR may directly regulate the transcription of certain miRNAs after castration. These miRNAs may contribute to the fat deposition phenotype. We therefore compared the expression of miRNAs between castrated and intact male pigs from full-sib pairs, in order to identify novel fat-deposition-related and differentially expressed miRNAs after castration. The purpose of this research is to gain new insight into fat deposition-related miRNAs in pigs, which will improve our understanding of fat deposition after castration.
Overview of the fat-deposition-related miRNA transcription profile
A total of 19,699,505 and 19,502,567 reads ranging in size from 15 nt-35 nt were retrieved from the F3 and F4 libraries, respectively. The size distribution of the clean reads is shown in Figure 1. Interestingly the size distribution of the small RNAs was similar between the small RNA libraries of the castrated (F3) and intact (F4) male pigs. By aligning the clean reads against the pig genome sequences (Sscrofa10.2), 4,822,509 reads in the F3 library and 1,704,144 reads in the F4 library were matched to the pig genome. A read was assigned to a miRNA by blasting against the non-miRNA databases. The clean reads were annotated into different categories. The number of 21–24 nt sequences (89.06%) was significantly greater than that of shorter or longer sequences. Almost half of the sequences in the F3 and F4 libraries were 22 nt, which is consistent with the known specificity of Dicer processing and the features of mature miRNAs (Additional file 1: Table S3).
Identification of known and novel porcine miRNAs and Mirtrons
A total of 366 unique miRNA genes comprising 174 known pre-miRNAs and 192 novel pre-miRNAs were identified. One hundred and sixty-seven pre-miRNAs overlapped between the two libraries. One hundred and fifty-three and 148 known miRNAs genes were identified in the F3 and F4 libraries (Table 1), respectively, of which 73 and 77 were novel miRNA*s (Additional file 2: Table A, Table B). The novel miRNA*s constituted 1.17% of the total known expressed miRNAs. From this data set, 116 novel miRNAs and 39 novel miRNA*s were identified in the F3 library, corresponding to 116 novel miRNA genes. Of those 116 candidate novel miRNAs in the F3 library, 24 were conserved in other mammals and the other 92 were considered to be pig-specific (Additional file 2: Table C, Table D). In the F4 library, 116 novel miRNAs and 27 novel miRNA*s were detected, corresponding to 116 novel miRNA genes. Of those 116 candidate novel miRNAs in the F4 library, 13 were conserved and 103 were pig-specific (Additional file 2: Table E, Table F). Hairpin structures of the partial novel miRNA precursors are shown in Additional file 1: Table S4.
The chromosome locations of known and candidate novel pre-miRNAs were determined based on the pig draft genomic sequence (Sscrofa10.2) from Ensembl (http://www.ensembl.org/Sus_scrofa/Info/Index). All of these miRNAs were mapped on autosomes or the X chromosome (Additional file 2). About 54.3% (63/116) and 44.8% (52/116) of the novel pre-miRNAs in the F3 and F4 libraries, respectively, were located in intergenic regions. Interestingly some pre-miRNAs had multiple copies. Nine pre-miRNAs in the F3 library and 19 pre-miRNAs in the F4 library were mapped to two positions on the same chromosome. In addition, ssc-miR-F4-S60 has three matched loci (Chr4:104560738-104560800+; Chr4:105581425-105581487+; Chr4:105429841-105429903+). The genomic density of pre-miRNAs is shown in Figure 2. In the F3 library, the average density of pre-miRNAs located on the auto-chromosome and X chromosome ranged from 0.05 to 0.3 miRNAs per 1 Mbp. The shortest chromosome, 18, and the longest chromosome, 1, encoded 13 and 19 pre-miRNAs, respectively, corresponding to 0.24 and 0.06 genomic densities. In the F4 library, the average density was from 0.04 to 0.37. The shortest and the longest chromosomes encoded nine and 18 pre-miRNAs, respectively, corresponding to 0.17 and 0.06 genomic densities. Based on the dataset, 38 and 30 pre-miRNAs were found on the X chromosome in the F3 and F4 libraries, respectively, and the corresponding average genomic density was 0.3 and 0.24.
In this study, mirtrons were identified as special pre-miRNAs, which constitute a recently described category of miRNA that relies on splicing rather than processing by the microprocessor complex to generate the RNAi pathway. Mirtrons were originally recognized in the fly and worm, but they have also been reported in mammals [19, 20]. From our dataset, five mirtrons (ssc-miR-F3-S25, ssc-miR-F3-S47, ssc-miR-7142-3p, ssc-miR-F3-S63 and ssc-miR-F3-S69) were identified. However, they were only found in the F3 library at low levels of expression (< 50 reads). These mirtrons need to be further validated with RT-qPCR to confirm that they are not false-positives caused by specific expression from one single animal.
Differentially expressed miRNAs in backfat between castrated and intact male pigs
One hundred and seventy-seven miRNAs, including 45 up-regulated and 132 down-regulated (Figure 3; Additional file 3), were found to have more than 2-fold differential expression between the castrated and intact male pigs (p-value < 0.001). In the F3 library, ssc-miR-143, ssc-miR-145, ssc-miR-199b, ssc-miR-103, ssc-miR-191, ssc-miR-10a, ssc-miR-320a, ssc-miR-152, ssc-miR-23a and ssc-miR-23b were the dominant expressed miRNAs, with the number of reads ranging from 30,354 to 1,307,953. The reads of dominant miRNAs constituted 82.7% of the total reads. In the F4 library, the 10 most abundant miRNAs (ssc-miR-320a, ssc-miR-21, ssc-miR-191, ssc-miR-143, ssc-miR-145, ssc-miR-423-5p, ssc-miR-7134-3p, ssc-miR-152, ssc-miR-195 and ssc-miR-193a) had reads numbering from 30,883 to 111,160, which contributed to 58.15% of the total miRNAs. Interestingly, the ratio of known miRNAs sequence reads to total reads was 0.9934 and the ratio of novel miRNAs was 0.0064 in the castrated male pigs, which differed from that in the intact male pigs, where the ratios were 0.9424 and 0.0576, respectively.
Expression of miRNA clusters in castrated and intact male pigs
Usually, miRNAs located consecutively within 10 kb of each other along the genome are considered as an miRNA cluster . In this study, the miRNAs were categorized into different clusters based on their chromosome locations (miRNA-miRNA distance < 10 kb). Briefly, there were 29 miRNAs clusters detected in the F3 library, each of which had 2–5 miRNAs. Ten novel miRNAs were grouped into eight miRNA clusters. In the F4 library, there were 22 miRNA clusters, each of which had two to four miRNAs, and six novel miRNAs were distributed into five miRNA clusters (Additional file 4). Clustered miRNAs can be coordinately transcribed as a polycistronic transcript or expressed individually. Based on the obtained results more than half of the miRNA clusters had a similar expression pattern (Additional file 1: Table S5). As an example the ssc-miR-181a-1 cluster includes ssc-miR-181a-1 and ssc-miR-181b-1, two different miRNAs according to the miRbase19.0 database. In this study both of them were found to have a similar level of expression. Another polycistronic miR-17 cluster showed a differential expression pattern. It contains six mature miRNAs: ssc-miR-17, ssc-miR-18a, ssc-miR-19a, ssc-miR-20, ssc-miR-19b-1 and ssc-miR-92a-1; however, only ssc-miR-17, ssc-miR-18a and ssc-miR-19a were expressed in castrated male pigs (Additional file 1: Table S5).
Sequence variations between miRNAs of castrated and intact male pigs
Porcine miRNAs frequently exhibit variation from their “reference” miRBase sequences (which have the highest level of expression), producing multiple mature variants that we herein refer to as isomiRs [22–24]. Sequence variations of miRNA*s and isomiRs were also observed in this study. There were 2,131 miRNA isomiRs and 795 miRNA* isomiRs found in the F3 library, and different isomiRs for a given miRNA and miRNA* ranged from two to 41 and two to 22, respectively. In the F4 library, 1,904 miRNA isomiRs and 738 miRNA* isomiRs were detected. The observed isomiR number varied from two to 35. There were 38 miRNA genes in the F3 library and 42 in the F4 library that did not show any mature miRNA or miRNA* isomiRs. In some cases, the most abundant sequence did not correspond exactly to the current porcine miRBase 19.0 reference sequence. As an example, the size of ssc-miR-205 in our library was 23 bp (read number was 4453), while the size of ssc-miR-205 in miRBase is 22 bp (read number is 3477). ssc-miR-205* was identified in this study; however, it was not found in the miRBase database (Additional file 5).
Prediction of target genes and pathways analysis
To better understand the biological function of the 19 up-regulated and 16 down-regulated miRNAs (> 2-fold change and p-value < 0.001) (Additional file 1: Table S6) in the castrated male pigs compared with intact male pigs, their target genes were predicted using Targetscan and the NCBI Entrez database. The predicted target genes were further classified to identify pathways that were actively regulated by miRNAs in adipose tissues according to KEGG functional annotations (Tables 2 and 3). The predicted target genes of miRNAs that participate in signaling pathways, as shown by DAVID KEGG analysis, are shown in Additional files 6 and 7. It is worth noting that the targets of the majority of up-regulated miRNAs in the castrated pigs belonged to the MAPK (mitogen-activated protein kinase 1) signaling pathway, which is closely related to the inhibition of lipogenesis  (Table 2, Figure 4). Another pathway targeted by over-represented miRNAs is the Wnt signaling pathway. This pathway is known to be involved in adipocyte lineage commitment, adipogenesis metabolism, and also inhibits adipogenesis through beta-catenin-dependent and -independent mechanisms. The Wnt/beta-catenin signaling pathway was reported to inhibit porcine adipogenic differentiation potential . Furthermore, 20 putative target genes, regulated by down-regulated miRNAs, are involved in the GnRH signaling pathway. Some target genes are also components of the insulin signaling pathway, which is related to the biological function of adipose tissue. Another important pathway targeted by down-regulated miRNAs in adipose tissue is the TGF-beta signaling pathway. This is associated with cellular functions, such as proliferation, apoptosis, differentiation and migration . Based on our results, the TGFBR2 gene was the target of down-regulated ssc-miR-23b, ssc-miR-320a, ssc-miR-103, ssc-miR-142-5p and miR-21. TGFBR2 is part of the fat deposition-related TGF-beta signaling pathway and is known to inhibit adipose differentiation of preadipocyte cell lines and primary cultures . All putative target genes and their involvement in signaling pathways are shown in Additional file 8.
Validation of miRNA expression with Stem-loop RT-qPCR
Stem-loop qPCR was adopted to validate the sequencing data. A total of eight miRNAs were validated in the castrated and intact male pigs: ssc-miR-21, ssc-miR-30a, sssc-miR-27a, ssc-miR-143, ssc-miR-103, ssc-miR-142-5p ssc-miR-F3-C13 and ssc-miR-7134-3p. The expression levels of ssc-miR-21, ssc-miR-142-5p, ssc-miR-27a, ssc-miR-7134-3p and ssc-miR-103 were significantly higher in the intact male pigs than in the castrated male pigs, while the expression levels of ssc-miR-30a, ssc-miR-143 and ssc-miR-F3-C13 were lower in the intact male pigs (Figure 5). These data were consistent between RT-qPCR and SOLiD sequencing. The detailed data is shown in Additional file 1: Table S7.
The domestic pig is important economically for meat production and as a model for comparative genomics and biomedical studies [28, 29]. The castration of domestic pigs has an important effect on the accumulation of adipose tissues [14, 15]. The molecular mechanism behind this is still unclear. Recent studies have provided some clues to solve this problem since they show that miRNAs may play a regulatory role in white adipose tissue development in animals. To gain new insight into the role of miRNAs in fat deposition in castrated pigs, novel and differentially expressed miRNAs were identified by SOLiD sequencing. The regulatory pathways of putative target genes were also analyzed in the castrated and intact male pigs. It was found that there were significantly different miRNA expression patterns in backfat tissue between castrated and intact male pigs. The differences can be described according to two aspects. First is the varied expression of up-regulated and down-regulated genes. One hundred and seventy-seven miRNAs, including 45 up- and 132 down-regulated, had more than 2-fold differential expression between the castrated and intact male pigs (p-value < 0.001). Second is the specific expression of some miRNAs in each library. We speculate that the above two differences might be a reason for the different fat deposition patterns of castrated and intact male pigs. In conclusion, miRNAs may play an important role in backfat deposition in castrated male pigs.
ssc-miR-143 had the highest level of expression in pig backfat tissue among all expressed miRNAs. Moreover, it was expressed at a higher level in castrated male pigs than in intact pigs (4.96-fold). This is consistent with ssc-miR-143 being the most highly expressed miRNA (representing 20% of the total miRNA expression) in the backfat of 240 day-old obese Rongchang pigs . In addition, miR-143 was up-regulated after differentiation induction in both human preadipocytes and mouse 3T3-L1 cells. It was also reported that the introduction of antisense oligonucleotides against miR-143 inhibited adipocyte differentiation by up-regulating MAPK1, which thereby reduced fat deposition. The above results suggested that the miR-143 plays a role in prompting adipocyte differentiation [8, 9, 31]. Our results together with the above published findings show that miR-143 is directly related to fat deposition in pigs. A possible mechanism for this action is that miR-143 down-regulates MAPK1, and increases PPARγ transcriptional activity by mediating its phosphorylation [32, 33]. In addition, MAPK1 is also a predicted target gene of ssc-miR-204, ssc-miR-129-5p, ssc-miR-194a, ssc-miR-181a, ssc-miR-130a and ssc-miR-101 (Additional file 9). We have previously shown that MAPK1 expression is significantly decreased, while PPARγ expression is increased in the backfat of castrated pigs compared with intact male pigs , and both genes were involved in MAPK signaling, which is an important pathway in fat deposition. It is likely that these miRNAs promote the generation of pig adipose tissue by down-regulating MAPK1 and up-regulating PPARγ.
In contrast with miRNAs regulating the MAPK signaling pathway, five of the miRNAs involved in the TGF-beta signaling pathway for fat deposition (ssc-miR-23b, ssc-miR-320a, ssc-miR-103, ssc-miR-142-5p and ssc-miR-21) had decreased levels of expression in castrated male pigs. However, the significantly different expression of miR-21 between castrated and intact male pigs (0 vs 83,244 reads) was unexpected because it is inconsistent with an increase in miR-21 expression during adipogenic differentiation of hASCs . The expression pattern of the Let-7 gene family is another conflicting result since the expression of the Let-7 gene family was undetected in this study, while previously it accounted for 7.7% of the total miRNA expression in pig fat . The above inconsistences might be due to different fat-deposition patterns and different experimental designs.
We also observed a reduction in some miRNA levels in the backfat of castrated male pigs compared with intact pigs. As an example, ssc-miR-27a and ssc-miR-378 were decreased by 3.14-fold and 6.06-fold in the castrated male pigs, respectively. This result is consistent with a previous report showing that ssc-miR-27a negatively regulates adipocyte differentiation by inhibiting the expression of PPARγ (Additional file 9), and that miR-378 expression is also negatively correlated with cattle backfat thickness . The fat mass and obesity associated (FTO) gene was reported to be associated with fat deposition in Italian Duroc pigs  and its inactivation can protect from obesity . FTO was predicted by this research as a potential target of ssc-\miR-185, ssc-F3-C29, ssc-miR-101, ssc-miR-152 and ssc-miR-150, which were down-regulated in the castrated male pigs (Additional file 9). These five down-regulated miRNAs might play an important role in lipogenesis by targeting FTO; however, this needs experimental verification.
Levels of androgen and expression of the androgen receptor gene (AR) dramatically decrease after castration . In this study, we focused attention on ssc-miR-30a and ssc-miR-30e, which target the AR gene (Additional file 9). Levels of ssc-miR-30a were significantly increased in the backfat tissue in castrated male pigs compared with intact male pigs. This was consistent with a recent report showing that inhibition of miR-30a blocked adipogenesis, while over-expression of miR-30a stimulated human adipogenesis . From these observations, we suggest that the increased fat deposition in castrated pigs may be caused by the inhibitory effect of miR-30a on AR expression since expression of fat deposition-related genes will also be relaxed if the transcription of AR and androgen are decreased. However, whether miR-30a can cleave the AR mRNA needs to be further investigated. This rationale can be also applied to ssc-miR-30e, which exhibited a significant increase in the backfat tissue in castrated male pigs.
Interestingly, the expression of known miRNAs and miRNA*s was detected in the constructed libraries. Most miRNA*s were expressed at very low levels. This is consistent with the fact that the expression of mature miRNAs is always higher than that of the miRNA*s. This may be due to the gene structure and/or spatial and temporal regulation. However, there is an inconsistent result from the F4 library; the expression of ssc-miR-30e* (16,200 reads) was higher than that of ssc-miR-30e (3,652 reads), indicating that this miRNA*s is abundantly expressed, potentially as a functional guide miRNA . The expression of miR-30e* was significantly decreased in white adipose tissue of mice fed a high-fat diet . Our results therefore suggest that the high level of miR-30e* expression inhibited fat deposition in the intact male pigs. The ssc-miR-F3-C37 miRNA and miRNA* had the same expression level. The underlying reason for this may be that the 5′ ends of both miRNA and miRNA* have similar stabilities and, therefore, the same chance to combine with RISC to avoid degradation . Similar expression patterns of miRNA genes have also been identified in other species, such as chicken, fruit flies [41–43]. Above results can be explained by a specific molecular mechanism; however, the quality of the different constructed libraries cannot be excluded as an explanation of the results.
Several studies have shown that isomiRs may function in animals . Here, the length of ssc-miR-205 varied from 18 to 23 nucleotides. The length variation occurred at the 3′ end of miRNA, the 5′ end or both. This was mainly in the form of missing nucleotides and/or terminal additions of nucleotides and likely resulted from variations in the pre-miRNA secondary structures. These variations result in variable cleavage sites for Dicer or Drosha. On the other hand, it cannot be excluded that these isomiRs could have been generated by the digestion of miR-205. The 5′ end variation may affect the seed region, which is at the 2nd-8th position of a mature miRNA. The target of a miRNA may alter owing to change in the nucleotides in this region. These end-sequence variations are intriguing as they may allow miRNA variants to perform distinct roles by influencing miRNA/target mRNA hybrid duplex formation.
Stem-loop RT-qPCR was adopted to validate the sequencing data. Generally, the RT-qPCR validated the sequencing results and the results from the two different methods were consistent. However, there were few differences in expression levels, for example, the expression of ssc-miR-21, ssc-miR-142-5p and ssc-miR-7134-3p were not detected in the castrated male pigs according to the SOLiD sequencing, but a trace amount of expression was detected by RT-qPCR in the castrated male pigs. Such discrepancies in detection levels might be caused by the sensitivity and capability of the different methods.
This study significantly expands the cohort of fat-related miRNAs expressed in the pig and also identifies miRNAs whose expression patterns alter significantly in response to castration. In total, 232 novel miRNAs were identified in the backfat tissues from both castrated and intact male pigs. Approximately 48.4% (177/366) of the miRNAs were found to have more than 2-fold differential expression between intact and castrated male pigs, with p-values < 0.001. Putative target genes of these miRNAs and regulation pathways suggest that these miRNAs may play a regulatory role in the adipose accumulation caused by castration. This research provides novel information for the study of human obesity.
Sample preparation and RNA extraction
The animals used in this study were three pairs of male full sibs from Pietrain sire × Landrace dam crossed pigs. The paired male pigs had similar birth weight. One pig in each pair was castrated at the age of 1 week. The backfat tissues were collected from the 6th-7th rib of castrated and intact male pigs at 23 weeks from birth, and the backfat thickness is significantly different between the castrated and intact male pigs (p-value < 0.01). More details about phenotype information are described in Additional file 1: Table S1. Total RNA was extracted from each sample using Trizol reagent (Invitrogen, CA, USA) according to the manufacturer’s instructions. The quality of total RNA was checked using a Nano drop 2000 (Thermo, MA, USA). The RNA samples were stored at −80°C until use. The experimental procedures were approved by the animal welfare committee of the State Key Laboratory for Agro-biotechnology of China Agricultural University (Approval number XK257).
Small RNA library construction and SOLiD sequencing
Two small RNA libraries were constructed from the castrated and intact male pigs. The sequencing procedures were performed as follows. First, the RNA was purified to enrich the small RNA fraction using the flashPAGETM or PAGE (Ambion, Austin, USA). Second, adapters were hybridized and ligated to the purified small RNAs for sequencing. Third, reverse transcription and amplification were conducted using the SOLiD™ Small RNA Expression Kit (Applied Biosystems, CA, USA). Finally, the 105–150 base-pair PCR products were selected to construct sequencing libraries, which were sequenced using the SOLiD™ 2.0 system (SOLiD™ Small RNA Expression Kit, Applied Biosystems, CA, USA) (Additional file 10).
Data analysis and miRNA annotation
Small RNA sequences were produced by SOLiD. Low-quality reads were strictly trimmed from the raw reads. After removing adaptor sequences, 15–35 bp sequences were retrieved from the initial data set for further analysis. The 15–35 bp sequences that matched genomic sequences (Sscrofa10.2) were chosen (allowing for up to two end-nucleotide mismatches) and used to search against non-coding RNA data (rRNA, tRNA, etc.). After removing non-miRNA sequences, the remaining sequences were considered as candidate miRNA reads. Subsequently, precursors of the candidate sequences obtained from the short read library were searched for against the pig genome (Ensembl Trace Archive) to determine their chromosome locations. These miRNA precursor sequences were then used for fold-back secondary structure prediction. If a hairpin structure with a free energy of hybridization lower than −20 kcal/mol was predicted, the RNA sequence was then subjected to MiPred analysis. This predicts whether the input sequence is a genuine pre-miRNA-like hairpin sequence. The candidate sequences that matched known animal mature miRNAs in miRBase were searched against the porcine genome to identify the known porcine miRNAs. All unannotated matched genomic sequences were considered as novel candidate miRNAs or miRNA*s. Further, all novel miRNAs were classified into conserved and pig-specific miRNAs according to their sequence conservation among species. The conserved novel miRNAs were named with a suffix, for example, ssc-miR-F3-C or ssc-miR-F4-C. The pig-specific novel miRNAs were named, for example, as ssc-miR-F3-S or ssc-miR-F4-S.
The high-throughput SOLiD sequencing approach provides a useful means to estimate the expression profiles of miRNA genes by measuring the sequence reads. Raw gene counts were normalized by trimmed mean of M-values (TMM) normalization method in the edgeR package [45, 46]. This is a two-step procedure of TMM normalization between compared libraries. First, the normalization factors were calculated by edgeR for each library. Second, the normalized read counts were calculated by the following formula: normalized read counts = raw read counts/(library size × normalization factor) × 106. The obtained normalized read counts were then directly used for subsequent analyses. DEGseq packages were used to identify differentially expressed genes between two compared groups. P-values in the experiment without biological repeats were calculated using DEGseq software. It takes the technique variation into Poisson distribution without biological repeats. Scatter plots were generated based on log2normalized read counts. Fold change = log2(normalized read counts in F3/normalized read counts in F4).
Biological analysis software for small RNAs
For the primary analysis, necessary sequence data and software were retrieved from RNA2MAP (RNA_pipeline_0.4.0; http://solidsoftwaretools.com/gf/project/rna2map/). The pig genome sequence data was obtained from http://www.ensembl.org/Sus_scrofa/Info/Index. Pig mRNAs were obtained from the NCBI Nucleotide database (http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/) using the input “pig mRNA” AND “Sus scrofa [porgn:_txid9823]”. Known miRNAs/miRNA hairpins were obtained from miRBase Version 19.0 (http://microrna.sanger.ac.uk/). RNAfold (RNAfold: http://rna.tbi.univie.ac.at/cgi-bin/RNAfold.cgi) was downloaded for the secondary structure prediction. Targetscan 5.1 (http://www.targetscan.org/) and miRanda version 3.1 (http://www.microrna.org/microrna/getMirnaForm) were used for target gene prediction.
miRNA Identification by Stem-loop RT-qPCR
To validate miRNAs identified by SOLiD sequencing, the stem-loop reverse-transcription quantitative polymerase chain reaction (RT-qPCR) assay was used to specifically detect mature miRNAs . RT-qPCR includes two steps. First, the stem-loop RT primer is hybridized to a miRNA molecule and then reverse transcribed with a MultScribe reverse transcriptase. Second, the RT products are quantified. RT-qPCR was performed according to the manufacturer’s protocol using the mirVana™ miRNA isolation kit (Ambion, Austin, USA). The expression was validated in the same animals as used for the SOLiD experiment and analyzed individually. 5S was used as the endogenous control, and the mixture of every sample cDNA was used as a calibrator. The primers for miRNA amplification are shown in Additional file 1: Table S2. A BioRad CFX96 (BioRad, CA, USA) was used to perform the quantitative analysis using SYBRGreen PCR master mix (Applied Biosystems, CA, USA). The relative quantification of miRNA expression was calculated using the standard curve-based method for relative real-time PCR, which has been previously described .
Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116: 281-297. 10.1016/S0092-8674(04)00045-5.
Ambros V: MicroRNA pathway in files and worms: growth, death, fat, stress, and timing. Cell. 2003, 113: 673-676. 10.1016/S0092-8674(03)00428-8.
Hackl H, Burkard TR, Sturn A, Rubio R, Schleiffer A, Tian S, Quackenbush J, Eisenhaber F, Trajanoski Z: Molecular processes during fat cell development revealed by gene expression profiling and functional annotation. Genome Biol. 2005, 6: R108-10.1186/gb-2005-6-13-r108.
Carrington JC, Ambros V: Role of microRNAs in plant and animal development. Science. 2003, 301: 336-338. 10.1126/science.1085242.
Aguilar AL, Piskol R, Beitzinger M, Zhu JY, Kruspe D, Aszodi A, Moser M, Englert C, Meister G: The small RNA expression profile of the developing murine urinary and reproductive system. FEBS Lett. 2010, 584: 4426-4434. 10.1016/j.febslet.2010.09.050.
Guo L, Lu Z: The fate of miRNA* strand through evolutionary analysis: implication for degradation as merely carrier strand or potential regulatory molecule?. PLoS One. 2010, 5 (6): e11387-10.1371/journal.pone.0011387.
Xu P, Vernooy SY, Guo M, Hay BA: The Drosophila microRNA Mir-14 suppresses cell death and is required for normal fat metabolism. Curr Biol. 2003, 13: 790-795. 10.1016/S0960-9822(03)00250-1.
Esau C, Kang X, Peralta E, Hanson E, Marcusson EG, Ravichandran LV, Sun Y, Koo S, Perera RJ, Jain R, Dean NM, Freier SM, Bennett CF, Lollo B, Griffey R: MicroRNA-143 regulates adipocyte differentiation. J Biol Chem. 2004, 279: 52361-52365. 10.1074/jbc.C400438200.
Xie H, Lim B, Lodish HF: MicroRNAs induced during adipogenesis that accelerate fat cell development are downregulated in obesity. Diabetes. 2009, 58: 1050-1057. 10.2337/db08-1299.
Kim SY, Kim AY, Lee HW, Son YH, Lee GY, Lee JW, Lee YS, Kim JB: MiR-27a is a negative regulator of adipocyte differentiation via suppressing PPARg expression. Biochem Biophys Res Commun. 2010, 392: 323-328. 10.1016/j.bbrc.2010.01.012.
Wang Q, Li YC, Wang J, Kong J, Qi Y, Quigg RJ, Li X: MiR-17-92 cluster accelerates adipocyte differentiation by negatively regulating tumor-suppressor Rb2/p130. Proc Natl Acad Sci USA. 2008, 105: 2889-2894. 10.1073/pnas.0800178105.
Karbiener M, Fischer C, Nowitsch S, Opriessnig P, Papak C, Ailhaud G, Dani C, Amri EZ, Scheideler M: MicroRNA miR-27b impairs human adipocyte differentiation and targets PPARg. Biochem Biophys Res Commun. 2009, 390: 247-251. 10.1016/j.bbrc.2009.09.098.
Sun T, Fu M, Bookout AL, Kliewer SA, Mangelsdorf DJ: MicroRNA let-7 regulates 3T3-L1 adipogenesis. Mol Endocrinol. 2009, 23: 925-931. 10.1210/me.2008-0298.
Malmfors B, Nilsson R: Meat quality traits of boars in comparison with castrates and gilts. Swed J Agric. 1978, 8: 209-217.
Cronin GM, Dunshea FR, Butler KL, McCauley I, Barnett JL, Hemsworth PH: The effects of immuno- and surgical-castration on the behaviour and consequently growth of group-housed, male finisher pigs. Appl Anim Behav Sci. 2003, 81: 111-126. 10.1016/S0168-1591(02)00256-3.
Mangelsdorf DJ, Thummel C, Beato M, Herrlich P, Schütz G, Umesono K, Blumberg B, Kastner P, Mark M, Chambon P, Evans RM: The nuclear receptor superfamily: the second decade. Cell. 1995, 83: 835-839. 10.1016/0092-8674(95)90199-X.
Ribas J, Ni X, Haffner M, Wentzel EA, Salmasi AH, Chowdhury WH, Kudrolli TA, Yegnasubramanian S, Luo J, Rodriguez R, Mendell JT, Lupold SE: miR-21: An androgen receptor regulated microRNA which promotes hormone dependent and independent prostate cancer growth. Cancer Res. 2009, 69 (18): 7165-7169. 10.1158/0008-5472.CAN-09-1448.
Mammi C, Calanchini M, Antelmi A, Cinti F, Rosano GM, Lenzi A, Caprio M, Fabbri A: Androgens and adipose tissue in males: a complex and reciprocal interplay. Int J Endocrinol. 2012, 2012: 789653-
Berezikov E, Chung WJ, Willis J, Cuppen E, Lai EC: Mammalian mirtron genes. Mol Cell. 2007, 28: 328-336. 10.1016/j.molcel.2007.09.028.
Chan SP, Slack FJ: And now introducing Mammalian Mirtrons. Dev Cell. 2007, 13: 605-607. 10.1016/j.devcel.2007.10.010.
Yu J, Wang F, Yang GH, Wang FL, Ma YN, Du ZW, Zhang JW: Human microRNA clusters: genomic organization and expression profile in leukemia cell lines. Biochem Biophys Res Commun. 2006, 349 (1): 59-68. 10.1016/j.bbrc.2006.07.207.
Landgraf P, Rusu M, Sheridan R, Sewer A, Iovino N, Aravin A, Pfeffer S, Rice A, Kamphorst AO, Landthaler M, Lin C, Socci ND, Hermida L, Fulci V, Chiaretti S, Foà R, Schliwka J, Fuchs U, Novosel A, Müller RU, Schermer B, Bissels U, Inman J, Phan Q, Chien M, Weir DB, Choksi R, De Vita G, Frezzetti D, Trompeter HI, et al: A mammalian microRNA expression atlas based on small RNA library sequencing. Cell. 2007, 129: 1401-1414. 10.1016/j.cell.2007.04.040.
Reese TA, Xia J, Johnson LS, Zhou X, Zhang W, Virgin HW: Identification of novel microRNA-like molecules generated from herpesvirus and host tRNA transcripts. J Virol. 2010, 84: 10344-10353. 10.1128/JVI.00707-10.
Sdassi N, Silveri L, Laubier J, Tilly G, Costa J, Layani S, Vilotte JL, Le Provost F: Identification and characterization of new miRNAs cloned from normal mouse mammary gland. BMC Genomics. 2009, 10: 149-10.1186/1471-2164-10-149.
Wang M, Wang JJ, Li J, Park K, Qian X, Ma JX, Zhang SX: Pigment epithelium-derived factor suppresses adipogenesis via inhibition of the MAPK/ERK pathway in 3T3-L1 preadipocytes. Am J Physiol Endocrinol Metab. 2009, 297 (6): E1378-E1387. 10.1152/ajpendo.00252.2009.
Li HX, Lou X, Liu RX, Yang GS: Roles of Wnt/β-catenin signaling in adipogenic differentiation potential of adipose-derived mesenchymal stem cells. Mol Cell Endocrinol. 2008, 291 (1–2): 116-124.
Kim YJ, Hwang SJ, Bae YC, Jung JS: MiR-21 regulates adipogenic differentiation through the modulation of TGF-beta signaling in mesenchymal stem cells derived from human adipose tissue. Stem Cells. 2009, 27 (12): 3093-3102.
Lunney JK: Advances in swine biomedical model genomics. Int J Biol Sci. 2007, 3: 179-184.
Critser JK, Laughlin MH, Prather RS, Riley LK: Proceedings of the conference on swine in biomedical resarch. Ilar J. 2009, 50: 89-94. 10.1093/ilar.50.1.89.
Li G, Li Y, Li X, Ning X, Li M, Yang G: MicroRNA identity and abundance in developing swine adipose tissue as determined by Solexa sequencing. J Cell Biochem. 2011, 112: 1318-1328. 10.1002/jcb.23045.
Takanabe R, Ono K, Abe Y, Takaya T, Horie T, Wada H, Kita T, Satoh N, Shimatsu A, Hasegawa K: Up-regulated expression of microRNA-143 in association with obesity in adipose tissue of mice fed high-fat diet. Biochem Biophys Res Commun. 2008, 376: 728-732. 10.1016/j.bbrc.2008.09.050.
Hu E, Kim JB, Sarraf P, Bruce MN, Spiegelman N: Inhibition of adipogenesis Through MAP Kinase-Mediated phosphorylation of PPARgamma. Science. 1996, 274: 2100-2103. 10.1126/science.274.5295.2100.
Adams M, Reginato MJ, Shao D, Lazar MA, Chatterjee VK: Transcriptional activation by peroxisome proliferator-activated receptor gamma is inhibited by phosphorylation at a consensus mitogen-activated protein kinase site. J Biol Chem. 1997, 272: 5128-5132. 10.1074/jbc.272.8.5128.
Yao YC: Study on the growth and development of porcine skeletal muscles and adipose tissues regulated by male gonad. PhD Dissertation. 2009, China: College Animal Science and Technology, China Agricultural University
Jin W, Dodson MV, Moore SS, Basarab JA, Guan LL: Characterization of microRNA expression in bovine adipose tissues: a potential regulatory mechanism of subcutaneous adipose tissue development. BMC Mol Biol. 2010, 11: 29-10.1186/1471-2199-11-29.
Fischer J, Koch L, Emmerling C, Vierkotten J, Peters T, Brüning JC, Rüther U: Inactivation of the FTO gene protects from obesity. Nature. 2009, 458: 894-898. 10.1038/nature07848.
Fontanesi L, Scotti E, Buttazzoni L, Davoli R, Russo V: The porcine fat mass and obesity associated (FTO) gene is associated with fat deposition in Italian Duroc pigs. Anim Genet. 2009, 40 (1): 90-93. 10.1111/j.1365-2052.2008.01777.x.
Zaragosi LE, Wdziekonski B, Brigand KL, Villageois P, Mari B, Waldmann R, Dani C, Barbry P: Small RNA sequencing reveals miR-642a-3p as a novel adipocyte-specific microRNA and miR-30 as a key regulator of human adipogenesis. Genome Biol. 2011, 12 (7): R64-10.1186/gb-2011-12-7-r64.
Chartoumpekis DV, Zaravinos A, Ziros PG, Iskrenova RP, Psyrogiannis AI, Kyriazopoulou VE, Habeos IG: Differential expression of MicroRNAs in adipose tissue after long-term high-fat diet-induced obesity in mice. PLoS One. 2012, 7 (4): e34872-10.1371/journal.pone.0034872.
Winter J, Jung S, Keller S, Gregory RI, Diederichs S: Many roads to maturity: microRNA biogenesis pathways and their regulation. Nat Cell Biol. 2009, 11: 228-234. 10.1038/ncb0309-228.
Glazov EA, Cottee PA, Barris WC, Moore RJ, Dalrymple BP, Tizard ML: A microRNA catalog of the developing chicken embryo identified by a deep sequencing approach. Genome Res. 2008, 18: 957964-
Aravin AA, Lagos-Quintana M, Yalcin A, Zavolan M, Marks D, Snyder B, Gaasterland T, Meyer J, Tuschl T: The small RNA profile during Drosophila melanogaster development. Dev Cell. 2003, 5: 337-350. 10.1016/S1534-5807(03)00228-4.
Schwarz DS, Hutvágner G, Du T, Xu Z, Aronin N, Zamore PD: Asymmetry in the assembly of the RNAi enzyme complex. Cell. 2003, 115: 199-208. 10.1016/S0092-8674(03)00759-1.
Fernandez-Valverde SL, Taft RJ, Mattick JS: Dynamic isomiR regulation in Drosophila development. RNA. 2010, 16: 1881-1888. 10.1261/rna.2379610.
Robinson MD, Oshlack A: A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol. 2010, 11: R25-10.1186/gb-2010-11-3-r25.
Wang L, Feng Z, Wang X, Wang X, Zhang X: DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010, 70 (1): 136-138.
Chen C, Ridzon DA, Broomer AJ, Zhou Z, Lee DH, Nguyen JT, Barbisin M, Xu NL, Mahuvakar VR, Andersen MR, Lao KQ, Livak KJ, Guegler KJ: Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res. 2005, 33: e179-10.1093/nar/gni178.
Huang J, Ju Z, Li Q, Hou Q, Wang C, Li J, Li R, Wang L, Sun T, Hang S, Gao Y, Hou M, Zhong J: Solexa sequencing of novel and differentially expressed microRNAs in testicular and ovarian tissues in Holstein cattle. Int J Biol Sci. 2011, 7: 1016-1026.
Authors thank Prof. Changxin Wu for his helpful discussion in the experimental design. Thank Dr Huda Al-Kateb for English correction. And we also would like to thank the three anonymous reviewers, for their valuable comments and suggestions. This work was supported by National High Technology and Science Development Plan of China (2011AA100302), Natural Science Foundation of China (NSFC, Grant 31072002,31372275), Program for New Century Excellent Talents in University (NCET-11-0480) and Program for Changjiang Scholar and Innovation Research Team in University (IRT1191).
The sequence data from this study have been deposited at NCBI Sequence Read Archive database, and can be accessed with the number SRP029942.
The authors declare that they have no competing interests.
MYF designed the study. JMH, GL, JBZ and JYW collected the samples. YB and JMH generated the data. MYF, YB and JMH analyzed the data. YB and MYF contributed to the writing of the manuscript. CKL did the Figures. All authors have read and approved the final manuscript.
Ying Bai, Jin-Ming Huang contributed equally to this work.
Electronic supplementary material
Additional file 1: Table S1:Information of experimental male pigs; Table S2. Stem-loop primers for RT-PCR; Table S3. Numbers of miRNAs and tags in the castrated and intact male pigs; Table S4. Stem-loop structures of the partial novel microRNA precursors; Table S5. Expression of miRNA clusters in castrated and intact male pigs; Table S6. 35 differential expressed miRNAs for prediction of target genes and pathways analysis; Table S7. The expression log ratios obtained by RT-qPCR vs SOLiD sequencing. (XLS 194 KB)
Additional file 2:Table A. The known miRNAs in backfat in castrated male pig. Table B. The known miRNAs in backfat in intact male pig. Table C. The novel miRNAs conserved with ortholog of known mammalian miRNAs in castrated male pigs (F3); Table D. The candidate novel pig-specific miRNAs in backfat in castrated male pigs; Table E: The novel miRNAs conserved with ortholog of known mammalian miRNAs in intact male pigs (F4); Table F. The candidate novel pig-specific miRNAs in backfat in intact male pigs. (XLS 128 KB)
Additional file 9:Target prediction for selected differential expressed miRNAs. (A) AR was predicted as potential target of miR-30e and miR-30a; (B) PPARG was predicted as potential target of miR-27a and miR-27b; (C) MAPK1 was predicted as potential target of miR-143, miR-129-5p and miR-204; (D) ssc-miR-185, ssc-miR-150, ssc-F3-C29, ssc-miR-101, ssc-miR-152 were predicted target to FTO. (PDF 75 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Bai, Y., Huang, J., Liu, G. et al. A comprehensive microRNA expression profile of the backfat tissue from castrated and intact full-sib pair male pigs. BMC Genomics 15, 47 (2014) doi:10.1186/1471-2164-15-47
- Male pig
- Fat deposition