|miRNA Expression Profile Analysis in Kidney of Different Porcine Breeds.|
|Jump to Full Text|
|PMID: 23372853 Owner: NLM Status: In-Data-Review|
|microRNAs (miRNAs) are important post-transcriptional regulators in eukaryotes that target mRNAs repressing their expression. The uncertain process of pig domestication, with different origin focuses, and the selection process that commercial breeds suffered, have generated a wide spectrum of breeds with clear genetic and phenotypic variability. The aim of this work was to define the miRNAs expression profile in kidney of several porcine breeds. Small RNA libraries from kidney were elaborated and high-throughput sequenced with the 454 Genome Sequencer FLX (Roche). Pigs used were classified into three groups: the European origin group (Iberian breed and European Wild Boar ancestor), European commercial breeds (Landrace, Large White and Piétrain breeds) and breeds with Asian origin (Meishan and Vietnamese breeds). A total of 229 miRNAs were described in the pig kidney miRNA profile, including 110 miRNAs out of the 257 previously described pig miRNAs and 119 orthologous miRNAs. The most expressed miRNAs in pig kidney microRNAome were Hsa-miR-200b-3p, Ssc-miR-125b and Ssc-miR-23b. Moreover, 5 novel porcine miRNAs and 3 orthologous miRNAs could be validated through RT-qPCR. miRNA sequence variation was determined in 116 miRNAs, evidencing the presence of isomiRs. 125 miRNAs were differentially expressed between breed groups. The identification of breed-specific miRNAs, which could be potentially associated to certain phenotypes, is becoming a new tool for the study of the genetic variability underlying complex traits and furthermore, it adds a new layer of complexity to the interesting process of pig evolution.|
|Oriol Timoneda; Ingrid Balcells; Jose Ignacio Núñez; Raquel Egea; Gonzalo Vera; Anna Castelló; Anna Tomàs; Armand Sánchez|
Related Documents :
|20542143 - Development of a reverse transcription-loop-mediated isothermal amplification (rt-lamp)...
10903033 - Study of cellular responses to polymeric biomaterials using the differential display me...
16888033 - The antineoplastic agent bryostatin-1 differentially regulates ifn-gamma receptor subun...
|Type: Journal Article Date: 2013-01-25|
|Title: PloS one Volume: 8 ISSN: 1932-6203 ISO Abbreviation: PLoS ONE Publication Date: 2013|
|Created Date: 2013-02-01 Completed Date: - Revised Date: -|
Medline Journal Info:
|Nlm Unique ID: 101285081 Medline TA: PLoS One Country: United States|
|Languages: eng Pagination: e55402 Citation Subset: IM|
|Departament de Genètica Animal, Centre de Recerca en AgriGenòmica (CRAG), Universitat Autònoma de Barcelona, Campus UAB, Bellaterra, Spain.|
|APA/MLA Format Download EndNote Download BibTex|
Journal ID (nlm-ta): PLoS One
Journal ID (iso-abbrev): PLoS ONE
Journal ID (publisher-id): plos
Journal ID (pmc): plosone
Publisher: Public Library of Science, San Francisco, USA
Copyright: 2013 Timoneda et al
Received Day: 26 Month: 10 Year: 2012
Accepted Day: 23 Month: 12 Year: 2012
collection publication date: Year: 2013
Electronic publication date: Day: 25 Month: 1 Year: 2013
Volume: 8 Issue: 1
E-location ID: e55402
PubMed Id: 23372853
Publisher Id: PONE-D-12-33324
|miRNA Expression Profile Analysis in Kidney of Different Porcine Breeds Alternate Title:miRNA Expression Profile Analysis in Pig Kidney|
|Jose Ignacio Núñez2|
1Departament de Genètica Animal, Centre de Recerca en AgriGenòmica (CRAG), Universitat Autònoma de Barcelona, Campus UAB, Bellaterra, Spain
2Centre de Recerca en Sanitat Animal (CReSA), UAB, Bellaterra, Spain
3Program Infection and Immunity, FISIB, Bunyola, Spain
|University of Iowa, United States of America
|Correspondence: * E-mail: email@example.com
[conflict] Competing Interests: The authors have declared that no competing interests exist.
Contributed by footnote: Bioinformatic analysis: RE GV. Conceived and designed the experiments: AS JIN. Performed the experiments: OT IB AC AT. Analyzed the data: OT. Wrote the paper: OT.
MicroRNAs (miRNAs) are small regulatory RNAs that play important roles in the regulation of gene expression , . These single-stranded non-coding RNAs that are approximately 22 nucleotides long are involved in post-transcriptional regulation mechanisms acting mainly through down-regulation of target messenger RNAs (mRNAs) in a wide range of biological and pathological processes , . Thus, miRNAs inhibit gene expression by blocking protein translation or inducing mRNA degradation –. mRNAs can be regulated by several miRNAs, and each miRNA can target hundreds of mRNAs in different binding sites , . The recent evidence about the genetic variability at both 5′ and 3′ ends of the mature miRNA sequence generates a large spectrum of sequence variants called miRNA isoforms or isomiRs –. Mechanisms involved in isomiRs generation and their biological relevance has increased the complexity of molecular mechanisms related to regulation of mRNA expression mediated by miRNAs.
Nowadays, a large number of miRNAs have been reported in animals, plants and viruses, with up to 18,226 entries in the miRBase database (v18, November 2011, URL: http://www.mirbase.org/) –. In mammalian genomes, the number of encoded miRNAs has been predicted to be up to 1000 miRNAs, comprising about 3% of all protein-coding genes . In mammals, nearly all miRNAs are conserved in closely related species , and may have homologous in distant species, suggesting that miRNA functions could also be conserved throughout the evolution of animal lineages . Several studies showed that variability in miRNA sequences has been lost over time since miRNAs have been described in unicellular eukaryotes, showing its deeper evolutionary history among eukaryotes .
The pig is an important livestock species, not only for its production in industry, but also as a suitable animal model for comparative genomics and biomedical studies . However, the number of porcine miRNAs available in public databases is still poor, with only 257 porcine miRNAs described in the pig genome compared to the completely sequenced human (1,921), mice (1,157), bovine (676), poultry (544), equine (360) or canine (289) genomes (miRBase v18). The description of the microRNAome (miRNAome) held by different tissues under different pathophysiological states has become of great interest in the last years –. Determining altered patterns of miRNA expression related to disease or specific treatments or conditions would be useful in order to identify differentially expressed miRNAs that could be used as novel biomarkers –. In pigs, several miRNAomes have been described in different tissues such as muscle, fat, heart, liver, thymus, intestine and testes –, or using in vitro cells models such as porcine PK-15 cells (derived from porcine kidney epithelial cells) and porcine dendritic cells –. However, the porcine kidney miRNAome has not yet been described although it is an essential organ involved in functions like blood filtering, gluconeogenesis and in the secretion of important hormones like erythropoietin, renine and vitamin D.
The process of pig domestication has been very complex. Independent geographical origins have been described  which have generated multiple phenotypically different breeds , . In particular, high differences are found regarding reproductive and meat production and quality traits between Asian and European breeds. In addition, pig has become an important production animal for human, and, consequently, it has been strongly selected for traits of economical interest. As a result of this high selection, many commercial breeds have been generated. Different consuming necessities in pig industry have expanded the large phenotypic variability due to different selected traits. Genetic diversity between pig breeds has been further studied – showing variability in gene expression and elucidating the artificial selection performed by porcine industry during pig domestication . However, the genetic variability held at miRNA sequences, which could imply changes at post-transcriptional level, has not yet been explored although it may lead to great phenotypic differences and/or pathological disorders.
Studies with mitochondrial DNA (mtDNA) assume significant differences between European and Asian pigs prior to domestication and several studies describe the phylogenetic relationships between European and Asian pig breeds, determining different groups into European diversity and other groups into Asian diversity , , . Other studies confirmed some introgression events of Asian pig breeds in most European commercial breeds and evidence that local Spanish breeds have never been introgressed with Asian neither with European commercial populations –. In this sense, studying miRNAs differentially expressed between European and Asian pig breeds could contribute to understand their phenotypical differences.
In the present work, the miRNAome of several pig breeds belonging to the pure European and Asian branches as well as European commercial breeds that have been introgressed with Asian genes has been undertaken. The main goals of this study were to describe the pig kidney miRNAome through high throughput sequencing (HTS) technology, to study the miRNA sequence variability among porcine breeds and to determinate new porcine miRNAs.
Kidney samples were collected from twenty pigs belonging to six different breeds and the European ancestor (Wild boar). Animals were classified into three main groups based on breed origin; (1) European breeds: Wild boar (WB, n = 2) and Iberian (IB, n = 4), (2) Asian breeds: Meishan (ME, n = 3) and Vietnamese (VT, n = 3), and (3) European commercial breeds with strong influence from Asian breeds: Landrace (LD, n = 3), Large White (LW, n = 2) and Piétrain (PT, n = 3). All samples were taken from slaughterhouse (PRIMAYOR, Mollerussa, Spain) under veterinary supervision and they were immediately snap-frozen in liquid nitrogen and stored at −80°C until use. Total RNA was isolated using TRIzol® reagent (Invitrogen, Carlsbad, USA) following the manufacturer's recommendations. Total RNA was quantified using ND 1000 Nanodrop® Spectrophotometer (Thermo Scientific, Wilmington, USA) and its quality was assessed on an Agilent 2100 Bioanalyzer using the RNA 6000 Nano kits (Agilent Technologies, Santa Clara, USA). RNA quality threshold was set to RNA Integrity Number (RIN) above seven. For RT-qPCR validations, additional kidney samples were added (n = 29: 2 WB, 4 IB, 4 LD, 4 LW, 4 PT, 4 ME and 7 VT).
Small RNA fraction (15–30 nt) was excised and isolated from denaturing 12.5% polyacrilamide gel electrophoresis (PAGE) using miSpikeTM (IDT®, Coralville, USA) as internal size marker. For each breed, 60 μg of total RNA were loaded on separate gels to avoid cross-contamination. Gels were stained with GelStar® Acid Nucleic Gel Stain (Lonza, Basel, Switzerland) for UV visualization. Excised small RNA fraction were purified using Performa® DTR gel filtration cartridges (EdgeBio, Gaithersburg, USA) and, then, isolated small RNAs were pooled by breed to construct small RNA libraries. Briefly, 3′ and 5′ linkers from miRCatTM kit (IDT, Coralville, USA) were ligated at both ends of the small RNAs in two separated reactions using a T4 RNA ligase without ATP (Fermentas, Germany) and T4 RNA ligase with ATP (Ambion, Austin, USA), respectively. Between 3′ and 5′ primer ligations, the 60 nt RNAs were purified by PAGE to eliminate unligated products. Then, linked products were used to perform a reverse transcription reaction using the SuperScriptTM III Reverse Transcriptase kit (InvitrogenTM, Carlsbad, USA) and the cDNA obtained was amplified with the Platinum® Taq DNA Polimerase High Fidelity kit (InvitrogenTM, Carlsbad, USA). PCRs were done with primers complementary to 3′ and 5′ linkers and, in addition, they included multiplex identifiers at the 5′ end (a five nucleotide sequence tag) to allow differentiation between libraries (Table S1). Obtained PCR products were purified using the QIAquick PCR Purification Kit (Qiagen®, Germany). Libraries were quantified with QubitTM fluorometer, Quant-ITTM (InvitrogenTM, Carlsbad, USA) and prepared to a 108 DNA molecules/μL concentration for sequencing by 454 Genome Sequencer FLX Titanium (Roche, Germany) following the manufacturer's protocol at DNA sequencing facilities at CRAG (Bellaterra, Spain).
Primers were trimmed from sequences and only those sequences between 15 and 29 nucleotides and with total number of sequences ≥3 were kept for further analysis. Sequences were compared to all available miRNA sequences (miRBase v18) using local Blast. Parameters were set to 100% identity and up to 4 mismatches allowed at the end of the sequences to assume variability on 3′ and 5′ ends .
Differences in miRNA expression between breed groups (European breeds: WB and IB, Asian breeds: ME and VT, and European commercial breeds: LD, LW and PT) were assessed using the miRNAs annotated through miRBase database. Total number of sequences obtained for each miRNA was normalised by library size (in counts per thousand) and, then, averaged by group. Fold changes (FC) between groups were calculated using normalised data. According to sequence count, nine miRNAs highly expressed (Hsa-miR-200b-3p, Hsa-miR-200c-3p, Ssc-miR-126, Ssc-miR-126*, Ssc-miR-99a, Ssc-miR-532-5p, Ssc-miR-92a, Ssc-miR-26a and Bta-miR-193b, n>100) and four miRNAs lowly expressed (Ssc-miR-423-5p, Ssc-miR-29c, Ssc-miR-486 and Ssc-let-7f, n<100) in the kidney miRNAome were selected to measure their expression levels by RT-qPCR. Hsa-let-7a, Hsa-miR-25 and Hsa-miR-93 were used as reference miRNAs based on literature , . Primers used were designed by Exiqon® (Denmark). For low expressed miRNAs and Bta-miR-193b primers were designed following the methodology suggested by Balcells et al  (Table 1), where it is described that miR-specific quantitative RT-qPCR with DNA primers is a highly specific, sensitive and accurate method for microRNA quantification.
Reverse transcription (RT) reactions were performed in duplicate with the universal cDNA synthesis kit (Exiqon, Denmark) using 300 ng of total RNA following the manufacturer's instructions. Non template controls (NTC) and minus poly(A) polymerase controls for each tissue were included. qPCRs were performed in duplicate using miRCURYTM LNATM Universal RT microRNA PCR kit (Exiqon, Vedbaek, Denmark) on an 7900HT Sequence Detection System (Applied Biosystems, Warrington, UK). Each qPCR was done in 20 µL total volume including 10 µL SYBR Green master mix (Exiqon, Vedbaek, Denmark), 0.5 µL of each LNA primer set (Exiqon, Denmark) and 5 µl of a 1∶20 dilution of the cDNA. To assess qPCR efficiency, standard curves with 10-fold serial dilutions of a pool of equal amounts of cDNA from all samples were included in each assay. Thermal profile was set as follows: 95°C for 10 min and 40 cycles at 95°C for 15 sec and 60°C for 60 sec. NTC and minus poly(A) polymerase controls were included. Melting curve analysis was included at the end of the qPCR to detect unspecific amplifications.
Quantities from each sample were obtained from the calibration (standard) curve added in each RT-qPCR reaction. GeNorm v.3.5 software  was used to examine the stability of the reference miRNAs (M<1.5) and to obtain a normalization factor (NF). The quantity obtained from each miRNA and sample was normalised by the NF and FCs were calculated in relation to the lowest normalised value. FCs were log2 transformed and expression data were analysed by a two-way analysis of variance (ANOVA) with the General Linear Models procedure of the Statistical Package for the Social Scientists (IBM® SPSS® Statistics 19; IBM Corporation, Armonk, USA) including the RT and breed as fixed factors. Significance threshold was set at α<0.05. Scheffe test was used to determine significant differential expression between breed groups.
All primer trimmed sequences were mapped in the pig genome retrieved from Ensembl Genome Browser (Ensembl release 62, ftp://ftp.ensembl.org/pub/release-62/fasta/sus_scrofa/dna/, Sscrofa9.62) considering 100% of alignment and identity (perfect match). Sequences that mapped only once in the genome and with a length from 19 to 23 nt were selected. Among these, only sequences with unknown annotation at Ensembl genome browser and with copy number (CN) higher than 2 were considered. They were clustered taking into account only the position in the genome. Hence, sequences positioned in the same region were grouped and the sequence with higher CN was selected as the reference sequence for each cluster. Fifteen clusters were considered novel miRNAs to be validated. Flanking regions (50 nt) of the selected reference sequences were used to predict pre-miRNA folding structure using MFold software  following the guidelines reported by Ambros et al. . At the end, 8 clusters were selected for RT-qPCR validation. RT reactions were performed in duplicate using total RNA as previously described by Balcells et al . Briefly, 600 ng of total RNA in a final volume of 20 μL including 2 μL of 10x poly(A) polymerase buffer, 0.1 mM of ATP, 0.1 mM of each dNTP, 1 μM of RT-primer, 200 U of M-MuLV Reverse Transcriptase (New England Biolabs, USA) and 2 U of poly(A) polimerase (New England Biolabs, USA) was incubated at 42°C for 1 hour and at 95°C for 5 minutes for enzyme inactivation. Non template controls (NTC) and minus poly(A) polymerase controls for each tissue were included.
DNA primers for each miRNA were designed following the methodology suggested by Balcells et al  (Table 1). qPCR reactions were performed in duplicate in 20 μL final volume including 10 μL FastStart Universal SYBR Green Master (Roche, Germany), 125–500 nM of each primer (Table 1) and 5 μL of a 1∶100 dilution of the cDNA. Standard curves were generated by 10 fold serial dilutions of a pool of all cDNAs in order to calculate the qPCR efficiency. For Cluster 2 (Cl-2), standard curve was performed by using 2-fold serial dilutions. qPCR settings and data analysis was performed as explained above. Hsa-let-7a, Hsa-miR-25 and Hsa-miR-103 were used as reference miRNAs , .
DIANA – microT v3.0 web server ,  was used to identify in silico potential mRNA targets for differentially expressed miRNAs. Porcine genes are not included in the current version of DIANA – microT v3.0 and predictions were based on the human mRNA: miRNA interactions assuming sequence conservation. In silico functional annotation of putative mRNA target genes for each miRNA were analyzed with WEB-based Gene Set Analysis Toolkit (WebGestalt, ). Predicted miRNA targets were functionally annotated through the biological process information supported by Gene Ontology (GO, ) and the pathways in which they were involved were described by using the Kyoto Encyclopedia of Genes and Genomes database (KEGG, , ). Over or under represented functional categories were identified with hypergeometric test corrected by the multiple test adjustment proposed by Benjamini & Hochberg . Significant threshold was set at α<0.05.
Small RNA libraries from kidney of six porcine breeds (Iberian, Landrace, Large White, Piétrain, Meishan and Vietnamese) and the European ancestor (Wild Boar) were sequenced in a high throughput 454 GS FLX sequencer (Roche). A total of 115,305 reads (corresponding to 17,913 unique sequences) containing short inserts of length ranging from 15 to 29 nt (corresponding to miRNA size) were obtained (Table 2). The total number of reads obtained per breed ranged from 13,704 to 23,815 counts; however, the Iberian breed yielded a much lower number of total counts (4,022). If we consider the number of unique sequences obtained per breed, as an average, 17% of total counts corresponded to unique sequences except for Iberian and Vietnamese breeds that had 39% and 38% of unique sequences, respectively.
Alignment of sequences to miRBase database revealed that 77.7% (89,594) of the total sequences yielded a positive match to a known miRNA. In total, 229 different miRNAs were found to be expressed in porcine kidney. Among them, 110 miRNAs corresponded to previously described porcine miRNAs and 119 to orthologous miRNAs. The percentage of detected miRNAs per breed varied from 55 to 83%, being the Iberian breed the library with less detected miRNAs and the Large White the most, in agreement with the total number of sequences obtained in these libraries. The proportion of porcine miRNAs and orthologous miRNAs was similar in all libraries, around 50% both. Sixty-two miRNAs had a CN below 3 and were excluded from the analysis. Overall, the porcine kidney miRNA profile was estimated to be of 167 miRNAs (Table S2). The most expressed miRNAs (CN>350, 0.30%) in porcine kidney are listed in Table 3. The most abundant miRNA was Hsa-miR-200b-3p (27,097, representing 23.50% of all porcine kidney miRNAs), followed by Ssc-miR-125b (8,809; 7.64%), Ssc-miR-23b (5,412; 4.69%), Ssc-miR-126 (5,274; 4.57%) and Ssc-miR-23a (5,156; 4.47%). The miRNA expression pattern was similar between breeds; however, some differences were found. As an example, Ssc-miR-125b was the most expressed miRNA in Iberian breed, followed by Ssc-miR-99a and Hsa-miR-200b-3p, and Ssc-miR-192 and Hsa-miR-200c-3p were the fourth and fifth most expressed miRNAs in Large White breed, respectively. In addition, Hsa-miR-200c-3p appeared as the fifth most expressed miRNA in Meishan breed.
All miRNA sequences showed variability at both ends although it was higher in the 3′ end. In some cases, nucleotide variability could be also detected inside the miRNA sequence. In general terms, for each miRNA there was one sequence more expressed called reference miRNA sequence and, then, there was a number of sequence variations with less number of copies (Figure 1). Nucleotide variation was found in 70% of the miRNAs expressed in kidney and only a small proportion of kidney miRNAs did not show isomiRs. Nevertheless, those miRNAs without sequence variations were low expressed, from 3 to 38 total reads. Overall, 116 miRNAs presented isomiRs with an average rate of 10 isomiRs per miRNA. The miRNA with more variants was Hsa-miR-200b-3p with 123 variants, followed by Ssc-miR-23b and Ssc-miR-125b, with 59 and 51 variants, respectively. These three miRNAs are, in fact, the most expressed miRNAs in porcine kidney (Table 3). In this sense, variants expression followed two main patterns: according to those miRNAs with more than 1,000 total reads, they were distributed in those miRNAs with a strong predominant isomiR, such as Hsa-miR-200b-3p, Ssc-miR-125b, Ssc-miR-23b, Ssc-miR-23a, Ssc-miR-192, Ssc-miR-10b, Ssc-miR-126* and Ssc-miR-10a, and those miRNAs where there is not a really strong predominant isomiR, like Ssc-miR-126, Ssc-miR-99a, Hsa-miR-200c-3p, Ssc-miR-30d and Ssc-miR-125a (Table S3). Furthermore, the most expressed isoform was, in most cases, the same in all breeds; it only differed between breeds in Ssc-miR-99a and Has-miR-200c-3p.
Interestingly, reference miRNA sequence for each miRNA was not always the same as the one described in miRBase database. This was the case of 72 miRNAs (Table S4). Sequence differences were mainly found at 3′ end of the miRNA (in 90% of cases) and consisted, basically, in the addition or deletion of one or two nucleotides. In Table 4 are shown those most expressed miRNA sequences (CN>100) differing from the miRBase reference sequence. These differences in miRNA sequences are not influenced by the miRNA abundance, as total reads in these 72 miRNAs include high expressed miRNAs (Ssc-miR-23a, 5,156 reads) and low expressed miRNAs (Mmu-miR-5100, 3 reads). Moreover, counts of these 72 miRNAs represent the 21% of total counts of kidney miRNA profile, a considerable amount.
Furthermore, the most abundant reference sequence for a specific miRNA varied also between breeds. We found 76 miRNAs for which the most abundant sequence in each breed was not the same (data not shown). Nonetheless, differences in isomiR predominance between breeds did not follow any pattern related to breed origin. In addition, for some miRNAs, the described sequence in miRBase was expressed at a low level or even could not be detected in the porcine kidney samples analysed in this study. This was the case of Bta-miR-193b, Ssc-miR-374a, Ssc-miR-145, Ssc-miR-362, and Hsa-miR-324-3p, that had a CN of 21, 13, 3, 0 and 0, respectively.
Differential miRNA expression analysis was performed between breed groups: EU, EA and AS. miRNAs were considered differentially expressed (DE) when fold change difference between groups was greater than 1.5 or when a miRNA was not expressed in one group although p-values were not significant (p-value<0.05). In this sense, 42 miRNAs (25%) were equally expressed among analysed groups whereas 125 miRNAs (75%) were DE between breed groups (Table S5). About 77% of the described miRNAs were present in all groups. Focusing on miRNAs that were differential expressed from one group to the remaining groups, 26 miRNAs were up regulated in EU, 25 in EA and 7 in AS. Conversely, 16 miRNAs were down regulated only in EU, 9 in EA and 31 in AS (Table S5). Looking at the most expressed miRNAs, Ssc-miR-125b and Ssc-miR-99a were up regulated in EU, Hsa-miR-200c-3p and Ssc-miR-192 were up regulated in EA and Hsa-miR-200b-3p, Ssc-miR-23a and Ssc-miR-23b were up regulated in AS. Several miRNAs displayed high fold changes (FC≥5) in at least one comparison (Table S5); however, their CN tended to be low, suggesting that FCs were overestimated. Additionally, looking at normalised counts for each breed (Table S6), there were some miRNAs with low frequency which their expression was higher in some breed, such as Hsa-miR-193b-5p, Hsa-miR-140-5p, Ssc-miR-19b, Ssc-miR-423-5p and Ssc-miR-24 in Iberian breed, Ssc-miR-151-5p, Hsa-miR-4454, Ssc-miR-199a-3p and Ssc-miR-374b-5p in Landrace breed, Ssc-miR-486 in Piétrain breed and Ssc-let-7f and Hsa-let-7i-5p in Meishan breed. Thirteen miRNAs displaying expression differences between breeds by HTS were selected to be validated through RT-qPCR (Table 5). For RT-qPCR validations, additional kidney samples were added (n = 29: 2 WB, 4 IB, 4 LD, 4 LW, 4 PT, 4 ME and 7 VT). Hsa-miR-200b-3p, Hsa-miR-200c-3p and Bta-miR-193b were chosen to be orthologous miRNAs and thus, not described in pig yet. Ssc-miR-99a, Bta-miR-193b and Ssc-miR-423-5p were selected to be up regulated in EU, Hsa-miR-200b-3p, Ssc-miR-532-5p and Ssc-let-7f to be up regulated in AS and, Hsa-miR-200c-3p, Ssc-miR-92a, Ssc-miR-26a, Ssc-miR-486 and Ssc-miR-29c to be up regulated in EA. Finally, Ssc-miR-126 and Ssc-miR-126* were also included to be suggestively up regulated in EA. These miRNAs were mapped in the pig genome sequence (Sscrofa 9.62) and miRNA folding was predicted for Hsa-miR-200b-3p and Hsa-miR-200c-3p using MFold software  following the guidelines described by Ambros et al. . Unfortunately, Bta-miR-193b could not be mapped in the pig genome and, consequently, its miRNA folding could not be predicted. RT-qPCR efficiencies obtained were high (ranging from 90% to 110%) and standard curves correlations were at least of 0.99, allowing us to perform a good gene expression quantification . Eight miRNAs (Hsa-miR-200b-3p, Hsa-miR-200c-3p, Ssc-miR-126, Ssc-miR-126*, Ssc-miR-99a, Bta-miR-193b, Ssc-miR-486 and Ssc-let-7f) were differentially expressed in at least one comparison by RT-qPCR (p-value<0.05) (Table 6). Relative quantification from RT-qPCR data determined that Hsa-miR-200b-3p, Bta-miR-193b and Ssc-let-7f were up regulated in EU while Ssc-miR-126, Ssc-miR-126* and Ssc-miR-99a were down regulated in AS. On the other hand, Hsa-miR-200c-3p and Ssc-miR-486 were down and up regulated in EA, respectively. Correlation between HTS and qPCR data was medium to low (0.52>r2>0.003). Therefore, not all the results obtained with both techniques are in agreement.
Unannotated sequences (486 unique sequences) that mapped in a unique region within the pig genome (Sscrofa 9.62) were clustered taking into account their position in the genome. After filtering by CN>2, 15 clusters remained for further analysis (Table S7). A compatible pre-miRNA folding structure could be successfully predicted in silico for 8 clusters (Cl-2, Cl-5, Cl-15, Cl-16, Cl-24, Cl-25, Cl-29 and Cl-38, Figure 2). Then, a validation protocol through RT-qPCR was designed for them. Five clusters (Cl-2, Cl-5, Cl-25, Cl-29 and Cl-38) were successfully amplified and, therefore, they have been considered as novel porcine miRNAs (Table 5). RT-qPCR efficiencies were high (ranging from 90% to 110%) and standard curves correlations were at least of 0.99. In addition, expression data showed that 4 clusters (Cl-5, Cl-25, Cl-29 and Cl-38) were differentially expressed (p-value<0.05) between breed groups (Table 6). Cl-5, Cl-29 and Cl-38 were down regulated in AS whereas Cl-25 was up regulated in EU.
To propose a nomenclature for the novel miRNAs, they were aligned to miRBase database accepting internal mismatches and some variation at the ends (Table 6). Cl-5 had a positive match with Hsa-miR-194-5p (p-value = 0.0004) with two nucleotide deletions at the ends, and also with Ssc-miR-194 (p-value = 0.002) although in this alignment one internal mismatch was detected in addition of the two nucleotide deletions at the ends. Cl-25 had homology with Hsa-miR-551a and Cl-29 did not match to any orthologous miRNA in miRBase database. Cl-38 matched with Bta-miR-1468 (p-value = 0.0003) with one nucleotide deletion in the 5′ end. Cl-2 was not differentially expressed through RT-qPCR but it matched with Hsa-miR-31-3p (p-value = 0.001) having two nucleotide deletion at the ends.
Target genes for the eight differentially expressed miRNAs (Hsa-miR-200b-3p, Hsa-miR-200c-3p, Ssc-miR-126, Ssc-miR-126*, Ssc-miR-99a, Bta-miR-193b, Ssc-miR-486 and Ssc-let-7f) were predicted in silico. A total of 1,884 target genes were identified (Table S8) and they were functionally analyzed through KEGG pathways database where 96 biological pathways were determined (data not shown). Interestingly, Hsa-miR-200b-3p and Hsa-miR-200c-3p, from the same miRNA family, practically share the same pathways, related to cellular processes, signal transduction and some biological processes, like immune, nervous, circulatory and endocrine systems. They have also been related to cancer pathways, like renal cell carcinoma, among others. Ssc-miR-126, Ssc-miR-126*, Ssc-miR-99a and Ssc-let-7f have also been described in some cancer pathways, emphasizing the importance of miRNAs in pathways related with cell cycle, growth and death. No significant related pathways were obtained from targets associated to Bta-miR-193b and Ssc-miR-486.
The present study represents the first work that describes the porcine kidney microRNAome in 6 different breeds and the European Ancestor, the Wild Boar, to determine differentially expressed miRNAs related to breed origin. A total of 110 out of 257 known porcine miRNAs in miRBase database have been determined in this study by high throughput sequencing. It is known that some miRNAs could be tissue specific  and this could explain why this study has only described the 43,28% of the porcine miRNAs reported up to date which is in agreement with other studies , . Interestingly, among all described kidney miRNAs, 119 mammalian orthologous miRNAs were identified. Assuming high miRNA conservation in mammals , , the high proportion of orthologous miRNAs (52%) suggests that there are still many miRNAs not described in pigs yet.
miRNA expression profile in kidney revealed that the most expressed miRNAs were Hsa-miR-200b-3p, Ssc-miR-125b and Ssc-miR-23b. Interestingly, Hsa-miR-200b-3p is the most expressed miRNA in all breeds except for Iberian breed, which Ssc-miR-125b (567 reads), Ssc-miR-99a (339 reads) and Hsa-miR-200b-3p (239 reads) are the first, the second and the third most expressed miRNAs, respectively. This variation of the expression in Iberian breed could be explained by many reasons. Firstly, some bias could be assumed because of the low number of reads obtained in this library in comparison with the other libraries (Table 2). However, the total reads of each miRNA in Iberian breed evidence that there is some tendency in its expression. Secondly, this change in the most expressed miRNA could also be explained as Iberian breed was exclusively originated from European Wild Boar and not mixed with any other ancestor neither being in contact with other commercial breeds in Europe.
Importantly, the most expressed miRNA in kidney was Hsa-miR-200b-3p which has not already been described in pig. Its sequence could not be found in the previous pig genome sequence (v9.2) but it could be now annotated in the last version (Sscrofa10, ftp://ftp.ensembl.org/pub/release-67/fasta/sus_scrofa/dna/) probably due to it was located on a gap in the pig genome sequence that has been solved in the new version. Its expression was significantly higher than other profiled miRNAs, as it had been detected 27,097 times corresponding to the 23.50% of the total reads. This expression pattern in miRNAs profile has been defined in more studies , where a miRNA is predominant over the rest meaning that it could have a significant role in the expressed tissue. It is hypothesised that it could governs or be implicated on the major constitutive functions carried out by this tissue.
The multiple sequences obtained for each specific miRNA were aligned to study the variability in the 3′ and 5′ miRNA ends. Most of the identified miRNAs (70%) showed nucleotide polymorphisms and variations in sequence length, according to previous studies , . One example in this study is Ssc-miR-23b, where the length varies from 18 to 25 nucleotides (Figure 1). This variability in the 5′ and 3′ ends could be favoured by processes prior to the high throughput sequencing sample preparation steps , being explained by the cleavage of the same miRNA precursor at different nucleotides during Drosha and Dicer processing. Another reason of this variability could be the nucleotide replacement and additions or deletions in the 5′ and 3′ miRNA ends. All this post-transcriptional mechanisms, called RNA editing, generate the isomiRs, defined as sequence variations in mature miRNAs and commonly reported in miRNA studies , . Mechanisms of miRNA sequence diversification are extensively reviewed by Berezikov , suggesting a functional meaning to all these processes. Furthermore, different patterns described on isomiR expression suggest different strategies in the regulation of the expression. Assuming the expression data obtained in this study, there are miRNAs where the main isomiR could represent the major regulation function and there are miRNAs which their regulation function could be assumed by more than one isomiR (Table S3).
Considering the proportion of unique sequences obtained per breed in reference to the total sequences, Iberian and Vietnamese breeds obtained the highest proportion (38% and 39%, respectively) showing that, if the number of miRNAs was not increased in these libraries, what was increased was the number of described isomiRs. The percentage of unique sequences in the other libraries and also in the whole data (all libraries) was between 15 and 19%. Thus, the variability in miRNA sequences seemed to be higher in European breeds like Iberian breed and in Asian breeds such as Vietnamese breed. In contrast, when European commercial breeds were created through the European and Asian breeds introgression, their variability could be reduced due to their strong selection in production traits, as it is indicated if we consider the number of isomiRs in each European commercial breed studied.
Looking at the most expressed isomiR by miRNA in each breed we can determine that largely it was the same isomiR in all studied breeds. However, in 76 miRNAs, like Ssc-miR-99a, the most expressed isomiR varies depending on the breed, being more expressed one isomiR (U0009355: [gene: 5′-AACCCGTAGATCCGATCTTGTG-3′]) in European and Asian breeds, and other isomiR (U0009353: [gene: 5′-AACCCGTAGATCCGATCTTGTGA-3′]) in European commercial breeds. In other cases, there is an agreement in all breeds for the most expressed isomiR except in one breed, like in Ssc-miR-21 and Ssc-miR-100 in Iberian breed, Ssc-miR-152 in Piétrain breed, Has-miR-200c-3p in Vietnamese breed or Hsa-let-7d-5p in Landrace breed. There are also few cases where the most expressed isomiR is different in almost each breed, such as in Ssc-miR-199a* and Ssc-miR-423-5p. All these different situations suggest that the genetic variability degree between breeds is high and it could play an important role in the post-transcriptional regulation mechanisms leading to perform different functions.
Another variation factor was that there were some differences between the most expressed isomiR in each miRNA and its concordance with the described sequence in miRBase database, diverging them in 72 miRNAs found in this study. As example, the most representative sequence for Ssc-miR-23a in this study was [gene: 5′-ATCACATTGCCAGGGATTTCCA-3′], found 3,116 times and it corresponded to Bta-miR-23a. The sequence described in miRBase database for Ssc-miR-23a ([gene: 5′-ATCACATTGCCAGGGATTTCC-3′]) was found only 622 times (Table 4). It also suggests that there is not a fixed predominant isomiR, but there is variability in isomiR expression given not only by many factors like age, tissue or disease , but also by other factors like breed .
Despite of all miRNA homologies found, this study was also focused on describing novel miRNAs. The strict methodology used proposed 8 novel miRNAs to be validated through RT-qPCR and finally 5 miRNAs were successfully validated. These confirmed novel porcine miRNAs were mapped once in the pig genome (Table S7) and their pre-miRNA folding were successfully predicted (Figure 2). Moreover, three out of the thirteen selected miRNAs from the kidney miRNAome to be amplified through RT-qPCR were orthologous (Hsa-miR-200b-3p, Hsa-miR-200c-3p and Bta-miR-193b, Table 5) and, therefore, they were also confirmed as new porcine miRNAs. Significant differential expression regarding breed groups was obtained by RT-qPCR in eight miRNAs: Hsa-miR-200b-3p, Hsa-miR-200c-3p, Ssc-miR-126, Ssc-miR-126*, Ssc-miR-99a, Bta-miR-193b, Ssc-miR-486 and Ssc-let-7f. However, expression results from RT-qPCR were not in agreement with the differential expression study from high throughput sequencing data in some cases. These differences could be explained by many factors, like some possible bias, the fact that not each miRNA is amplified in the same proportion with high throughput sequencing ,  and the presence of isomiRs, which make difficult to perform an accurate filtering of sequences. Other factors to be considered are some diversity and complexity in the miRNAs nomenclature between species listed in the miRNAs databases , which make difficult to perform a proper sequence classification. In this sense, cluster analysis methodology for novel miRNA detection was chosen to avoid problems with miRNA nomenclature. It is also important to consider the specificity in qPCR amplification, where primers designed in each miRNA were exclusive for the most expressed isomiR in each miRNA, amplifying only that one. Overall, the relevant topic of next generation sequencing about if it can be treated as expression data taking into account its large amount of reads, is still discussed . The slightly different expression results between high throughput sequencing and RT-qPCR data reveal a change of expression tendency in some miRNAs evidencing the “semi-quantitative” nature of high throughput sequencing methods. Low correlations (from 0.52 to 0.003) were determined between high throughput sequencing by 454 (Roche, Germany) and RT-qPCR data. In this sense, the bias in this study could be explained by the lower number of sequences obtained. For instance, Hsa-miR-200b-3p appears in a higher expression in Asian breeds in sequencing data while it is more expressed in European breeds according to real time RT-qPCR data. Furthermore, Ssc-miR-126* was more represented in European commercial breeds, whereas in RT-qPCR it was up regulated in European breeds. The same case happened with Ssc-miR-126, a conclusive fact taking into account that both miRNAs are transcribed together. There are similar cases with low expressed miRNAs (according to sequence count) such as Ssc-let-7f, which was more expressed in Asian breeds but appears up regulated in European breeds in real time RT-qPCR data.
Analysing the target pathways of differentially expressed miRNAs between breed groups, it was found that Hsa-mir-200b-3p, the most expressed miRNA in this study and up regulated in European breeds, was related in several kidney diseases like tubulointerstitial fibrosis  or hypertensive glomerulosclerosis , elucidating the importance of this miRNA in kidney physiology. Targets of Hsa-miR-200b-3p were involved in several pathways like renal cell carcinoma, associated to some oncogenes such as MET, or tumors suppressors like VHL, FH and BHD. Interestingly, miRNA targets were related with the vascular smooth muscle contraction and calcium signaling pathways. Both pathways are important in muscular growth and development in which Hsa-miR-200b-3p, being more expressed in European pig breeds, could contribute to a greater efficiency. Furthermore, European commercial breeds are precisely selected for production traits and strongly characterized for animal production. Kidney and muscle are strongly related because kidney is in charge of secreting erythropoietin, the hormone responsible to activate erythropoiesis and also to transport oxygen to the muscles and other tissues, and thus, favoring muscular growth. Hsa-miR-200c-3p, which belongs to the same family that Hsa-miR-200b-3p, is also differentially expressed being up regulated in Asian breeds. Some studies have related it to some mechanisms of cancer development in many organs, such as pancreas, bladder, ovaries, prostate or breast. –. Functional analysis related targets of Hsa-miR-200c-3p with neoplasias, specifically in renal cell carcinoma. Furthermore, Hsa-miR-200c-3p is also involved in some pathways related to reproduction, like oocyte meiosis or progesterone-mediated oocyte maturation. Breeds with an Asian origin have better phenotypes for reproduction traits and, interestingly, Asian breeds have Hsa-miR-200c-3p up regulated suggesting its involvement in the regulation of the reproduction traits. However, Ssc-miR-99a and Hsa-miR-200b-3p targets are also related to these reproduction pathways, both up regulated miRNAs in European breeds which could play an antagonistic role in European breeds. Some targets of Ssc-miR-126, a down regulated miRNA in Asian breeds, participate in cellular processes like focal adhesion and regulation of actin cytoskeleton, while Ssc-miR-126*, also down regulated in Asian breeds, regulate targets implicated with cellular processes like adherens, gap and tight junction. Moreover, both miRNAs also have a strong relation with many cancer pathways, which were related in all miRNAs. Ssc-let-7f is an ubiquitous miRNA from let-7 family, related to many pathways involved in wide range of physiological functions (cellular, metabolic and environmental information processes) and also to many diseases, like renal cell carcinoma , .
This study opens a wide field about the miRNA specificity, which is not only by tissue, age or species level, but also breed factor is crucial to manage the phenotypic changes and, therefore, the genetic variability. The role of miRNAs in gene expression regulation is much more complex than it was thought initially, being necessary to study in depth which roles isomiRs play and to uncover all factors that make the post-transcriptional regulation through miRNAs so complex.
Adaptors used for the construction of each library.1: Forward adaptor, 2: Reverse adaptor. The 5 nt code used for each library is in bold in the adaptor sequence. Reverse adaptor was used in all breeds libraries.
Click here for additional data file (pone.0055402.s001.doc)
Complete miRNA profile in healthy kidney of 6 porcine breeds and the European ancestor from 454 GS FLX run data. IB: Iberian breed, WB: European Wild Boar, LD: Landrace breed, LW: Large White breed, PT: Piétrain breed, ME: Meishan breed, VT: Vietnamese breed. miRNA name represents the most expressed sequence in the cluster. Bta: Bos taurus, Dre: Danio rerio, Eca: Equus caballus, Hsa: Homo sapiens, Mdo: Monodelphis domestica, Mmu: Mus musculus, Rno: Ratus norvegicus, Sha: Sarcophilus harrisii, Ssc: Sus scrofa.
Click here for additional data file (pone.0055402.s002.doc)
Summary of the isomiR distribution in the most expressed miRNAs (CN>1,000). miRNA name represents the most expressed sequence in the cluster. Hsa: Homo sapiens, Ssc: Sus scrofa.
Click here for additional data file (pone.0055402.s003.doc)
Described miRNAs from HTS presenting variation between the reference miRNA sequence and the miRBase described miRNA sequence. Bta: Bos taurus, Eca: Ecuus caballus, Hsa: Homo sapiens, Mdo: Monodelphis domestica, Mmu: Mus musculus, Rno: Rattus norvegicus, Sha: Sarcophilus harrisii, Ssc: Sus scrofa. Marked in bold the nucleotide variation between pair sequences.
Click here for additional data file (pone.0055402.s004.doc)
List of differentially expressed miRNAs (Fold Change >1.5times) between porcine breed groups. EU: European breeds; EA: European commercial breeds; AS: Asian breeds. Bta: Bos taurus, Dre: Danio rerio, Eca: Equus caballus, Hsa: Homo sapiens, Mdo: Monodelphis domestica, Mmu: Mus musculus, Rno: Ratus norvegicus, Sha: Sarcophilus harrisii, Ssc: Sus scrofa. Positive and negative signs indicate that the level of gene expression is higher for the first or the second group of the test, respectively.
Click here for additional data file (pone.0055402.s005.doc)
Normalised 454 GS FLX run data by library in counts per thousand. IB: Iberian breed, WB: European Wild Boar, LD: Landrace breed, LW: Large White breed, PT: Piétrain breed, ME: MeiShan breed, VT: Vietnamese breed. miRNA name represents the most expressed sequence in the cluster. Bta: Bos taurus, Dre: Danio rerio, Eca: Equus caballus, Hsa: Homo sapiens, Mdo: Monodelphis domestica, Mmu: Mus musculus, Rno: Ratus norvegicus, Sha: Sarcophilus harrisii, Ssc: Sus scrofa.
Click here for additional data file (pone.0055402.s006.doc)
Resulting clusters from sequence analysis for novel miRNAs discovery.1: chromosome: start position: end position: strand. Sequences mapped at pig genome sequence (Sscrofa 9.62).
Click here for additional data file (pone.0055402.s007.doc)
Putative target genes of the eight differentially expressed miRNAs analised by qRT-PCR. Bta: Bos taurus, Hsa: Homo sapiens, Ssc: Sus scrofa. Potential mRNA target genes for differentially expressed miRNAs predicted in silico with DIANA – microT v3.0 web server.
Click here for additional data file (pone.0055402.s008.doc)
|1.||Ambros V, (Year: 2008) The evolution of our thinking about microRNAs. Nature Medicine14: 1036–1040|
|2.||Ruvkun G, (Year: 2008) The perfect storm of tiny RNAs. Nature Medicine14: 1041–1045|
|3.||Bartel DP, (Year: 2009) MicroRNAs: Target recognition and regulatory functions. Cell136: 215–23319167326|
|4.||Krol J,, Loedige I,, Filipowicz W, (Year: 2010) The widespread regulation of microRNA biogenesis, function and decay. Nature Reviews. Genetics11: 597–61020661255|
|5.||Bartel DP, (Year: 2004) MicroRNAs: Genomics, biogenesis, mechanism, and function. Cell116: 281–29714744438|
|6.||Guo H,, Ingolia NT,, Weissman JS,, Bartel DP, (Year: 2010) Mammalian microRNAs predominantly act to decrease target mRNA levels. Nature466: 835–84020703300|
|7.||Huntzinger E,, Izaurralde E, (Year: 2011) Gene silencing by microRNAs: Contributions of translational repression and mRNA decay. Nature Reviews. Genetics12: 99–11021245828|
|8.||Ambros V, (Year: 2004) The functions of animal microRNAs. Nature Medicine431: 350–355|
|9.||Lewis BP,, Burge CB,, Bartel DP, (Year: 2005) Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell120: 15–2015652477|
|10.||Morin RD,, O'Connor MD,, Griffith M,, Kuchenbauer F,, Delaney A,, et al. (Year: 2008) Application of massively parallel sequencing to microRNA profiling and discovery in human embryonic stem cells. Genome Research18: 610–62118285502|
|11.||Saal S,, Harvey SJ, (Year: 2009) MicroRNAs and the kidney: Coming of age. Current Opinion in Nephrology and Hypertension18: 317–32319424061|
|12.||Li JY,, Yong TY,, Michael MZ,, Gleadle JM, (Year: 2010) Review: The role of microRNAs in kidney disease. Nephrology (Carlton, Vic.)15: 599–608|
|13.||Kozomara A,, Griffiths-Jones S, (Year: 2011) miRBase: Integrating microRNA annotation and deep-sequencing data. Nucleic Acids Research39: D152–721037258|
|14.||Griffiths-Jones S,, Saini HK,, van Dongen S,, Enright AJ, (Year: 2008) miRBase: Tools for microRNA genomics. Nucleic Acids Research36: D154–817991681|
|15.||Griffiths-Jones S,, Grocock RJ,, van Dongen S,, Bateman A,, Enright AJ, (Year: 2006) miRBase: MicroRNA sequences, targets and gene nomenclature. Nucleic Acids Research34: D140–416381832|
|16.||Berezikov E, (Year: 2011) Evolution of microRNA diversity and regulation in animals. Nature Reviews. Genetics12: 846–86022094948|
|17.||Lim LP,, Glasner ME,, Yekta S,, Burge CB,, Bartel DP, (Year: 2003) Vertebrate microRNA genes. Science (New York, N.Y.)299: 1540|
|18.||Tarver JE,, Donoghue PC,, Peterson KJ, (Year: 2012) Do miRNAs have a deep evolutionary history?BioEssays: News and Reviews in Molecular, Cellular and Developmental Biology34: 857–866|
|19.||Lunney JK, (Year: 2007) Advances in swine biomedical model genomics. International Journal of Biological Sciences3: 179–18417384736|
|20.||Glazov EA,, Cottee PA,, Barris WC,, Moore RJ,, Dalrymple BP,, et al. (Year: 2008) A microRNA catalog of the developing chicken embryo identified by a deep sequencing approach. Genome Research18: 957–96418469162|
|21.||Huang TH,, Zhu MJ,, Li XY,, Zhao SH, (Year: 2008) Discovery of porcine microRNAs and profiling from skeletal muscle tissues during development. PloS One3: e322518795099|
|22.||Ribeiro-dos-Santos A,, Khayat AS,, Silva A,, Alencar DO,, Lobato J,, et al. (Year: 2010) Ultra-deep sequencing reveals the microRNA expression pattern of the human stomach. PloS One5: e1320520949028|
|23.||Shao NY,, Hu HY,, Yan Z,, Xu Y,, Hu H,, et al. (Year: 2010) Comprehensive survey of human brain microRNA by deep sequencing. BMC Genomics11: 40920591156|
|24.||Xie SS,, Li XY,, Liu T,, Cao JH,, Zhong Q,, et al. (Year: 2011) Discovery of porcine microRNAs in multiple tissues by a solexa deep sequencing approach. PloS One6: e1623521283541|
|25.||Chen X,, Ba Y,, Ma L,, Cai X,, Yin Y,, et al. (Year: 2008) Characterization of microRNAs in serum: A novel class of biomarkers for diagnosis of cancer and other diseases. Cell Research18: 997–100618766170|
|26.||Chow TF,, Youssef YM,, Lianidou E,, Romaschin AD,, Honey RJ,, et al. (Year: 2010) Differential expression profiling of microRNAs and their potential involvement in renal cell carcinoma pathogenesis. Clinical Biochemistry43: 150–15819646430|
|27.||Godwin JG,, Ge X,, Stephan K,, Jurisch A,, Tullius SG,, et al. (Year: 2010) Identification of a microRNA signature of renal ischemia reperfusion injury. Proceedings of the National Academy of Sciences of the United States of America107: 14339–1434420651252|
|28.||Patel N,, Sauter ER, (Year: 2011) Body fluid micro(mi)RNAs as biomarkers for human cancer. Journal of Nucleic Acids Investigation2: e1|
|29.||Fendler A,, Stephan C,, Yousef GM,, Jung K, (Year: 2011) MicroRNAs as regulators of signal transduction in urological tumors. Clinical Chemistry57: 954–96821632885|
|30.||Catto JW,, Alcaraz A,, Bjartell AS,, De Vere White R,, Evans CP,, et al. (Year: 2011) MicroRNA in prostate, bladder, and kidney cancer: A systematic review. European Urology59: 671–68121296484|
|31.||Nielsen M,, Hansen JH,, Hedegaard J,, Nielsen RO,, Panitz F,, et al. (Year: 2010) MicroRNA identity and abundance in porcine skeletal muscles determined by deep sequencing. Animal Genetics41: 159–16819917043|
|32.||Reddy AM,, Zheng Y,, Jagadeeswaran G,, Macmil SL,, Graham WB,, et al. (Year: 2009) Cloning, characterization and expression analysis of porcine microRNAs. BMC Genomics10: 6519196471|
|33.||Cho IS,, Kim J,, Seo HY,, Lim do H,, Hong JS,, et al. (Year: 2010) Cloning and characterization of microRNAs from porcine skeletal muscle and adipose tissue. Molecular Biology Reports37: 3567–357420180025|
|34.||Sharbati S,, Friedlander MR,, Sharbati J,, Hoeke L,, Chen W,, et al. (Year: 2010) Deciphering the porcine intestinal microRNA transcriptome. BMC Genomics11: 27520433717|
|35.||Lian C, Sun B, Niu S, Yang R, Liu B, et al. . (2012) A comparative profile of the microRNA transcriptome in immature and mature porcine testes using solexa deep sequencing. The FEBS Journal.|
|36.||Wu YQ,, Chen DJ,, He HB,, Chen DS,, Chen LL,, et al. (Year: 2012) Pseudorabies virus infected porcine epithelial cell line generates a diverse set of host microRNAs and a special cluster of viral microRNAs. PloS One7: e3098822292087|
|37.||Anselmo A,, Flori L,, Jaffrezic F,, Rutigliano T,, Cecere M,, et al. (Year: 2011) Co-expression of host and viral microRNAs in porcine dendritic cells infected by the pseudorabies virus. PloS One6: e1737421408164|
|38.||Li S,, Ran XQ,, Xu L,, Wang JF, (Year: 2011) microRNA and mRNA expression profiling analysis of dichlorvos cytotoxicity in porcine kidney epithelial PK15 cells. DNA and Cell Biology30: 1073–108321702683|
|39.||Larson G,, Dobney K,, Albarella U,, Fang M,, Matisoo-Smith E,, et al. (Year: 2005) Worldwide phylogeography of wild boar reveals multiple centers of pig domestication. Science (New York, N.Y.)307: 1618–1621|
|40.||Larson G,, Albarella U,, Dobney K,, Rowley-Conwy P,, Schibler J,, et al. (Year: 2007) Ancient DNA, pig domestication, and the spread of the neolithic into europe. Proceedings of the National Academy of Sciences of the United States of America104: 15276–1528117855556|
|41.||Larson G,, Liu R,, Zhao X,, Yuan J,, Fuller D,, et al. (Year: 2010) Patterns of east asian pig domestication, migration, and turnover revealed by modern and ancient DNA. Proceedings of the National Academy of Sciences of the United States of America107: 7686–769120404179|
|42.||Amaral AJ,, Ferretti L,, Megens HJ,, Crooijmans RP,, Nie H,, et al. (Year: 2011) Genome-wide footprints of pig domestication and selection revealed through massive parallel sequencing of pooled DNA. PloS One6: e1478221483733|
|43.||Ojeda A,, Ramos-Onsins SE,, Marletta D,, Huang LS,, Folch JM,, et al. (Year: 2011) Evolutionary study of a potential selection target region in the pig. Heredity106: 330–33820502482|
|44.||Giuffra E,, Kijas JM,, Amarger V,, Carlborg O,, Jeon JT,, et al. (Year: 2000) The origin of the domestic pig: Independent domestication and subsequent introgression. Genetics154: 1785–179110747069|
|45.||Gao Y,, Zhang YH,, Jiang H,, Xiao SQ,, Wang S,, et al. (Year: 2011) Detection of differentially expressed genes in the longissimus dorsi of northeastern indigenous and large white pigs. Genetics and Molecular Research: GMR10: 779–79121563072|
|46.||Fang M,, Andersson L, (Year: 2006) Mitochondrial diversity in european and chinese pigs is consistent with population expansions that occurred prior to domestication. Proceedings.Biological Sciences/the Royal Society273: 1803–181016790414|
|47.||Clop A,, Amills M,, Noguera JL,, Fernandez A,, Capote J,, et al. (Year: 2004) Estimating the frequency of asian cytochrome B haplotypes in standard european and local spanish pig breeds. Genetics, Selection, Evolution: GSE36: 97–104|
|48.||Alves E,, Ovilo C,, Rodriguez MC,, Silio L, (Year: 2003) Mitochondrial DNA sequence variation and phylogenetic relationships among iberian pigs and other domestic and wild pig populations. Animal Genetics34: 319–32414510666|
|49.||Kim KI,, Lee JH,, Li K,, Zhang YP,, Lee SS,, et al. (Year: 2002) Phylogenetic relationships of asian and european pig breeds determined by mitochondrial DNA D-loop sequence polymorphism. Animal Genetics33: 19–2511849133|
|50.||Timoneda O,, Balcells I,, Cordoba S,, Castello A,, Sanchez A, (Year: 2012) Determination of reference microRNAs for relative quantification in porcine tissues. PloS One7: e4441322970213|
|51.||Peltier HJ,, Latham GJ, (Year: 2008) Normalization of microRNA expression levels in quantitative RT-PCR assays: Identification of suitable reference RNA targets in normal and cancerous human solid tissues. RNA (New York, N.Y.)14: 844–852|
|52.||Balcells I,, Cirera S,, Busk PK, (Year: 2011) Specific and sensitive quantitative RT-PCR of miRNAs with DNA primers. BMC Biotechnology11: 7021702990|
|53.||Vandesompele J,, De Preter K,, Pattyn F,, Poppe B,, Van Roy N,, et al. (Year: 2002) Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biology3: RESEARCH003412184808|
|54.||Zuker M, (Year: 2003) Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Research31: 3406–341512824337|
|55.||Ambros V,, Bartel B,, Bartel DP,, Burge CB,, Carrington JC,, et al. (Year: 2003) A uniform system for microRNA annotation. RNA (New York, N.Y.)9: 277–279|
|56.||Maragkakis M,, Alexiou P,, Papadopoulos GL,, Reczko M,, Dalamagas T,, et al. (Year: 2009) Accurate microRNA target prediction correlates with protein repression levels. BMC Bioinformatics10: 29519765283|
|57.||Maragkakis M,, Reczko M,, Simossis VA,, Alexiou P,, Papadopoulos GL,, et al. (Year: 2009) DIANA-microT web server: Elucidating microRNA functions through target prediction. Nucleic Acids Research37: W273–619406924|
|58.||Zhang B,, Kirov S,, Snoddy J, (Year: 2005) WebGestalt: An integrated system for exploring gene sets in various biological contexts. Nucleic Acids Research33: W741–815980575|
|59.||Ashburner M,, Ball CA,, Blake JA,, Botstein D,, Butler H,, et al. (Year: 2000) Gene ontology: Tool for the unification of biology. the gene ontology consortium. Nature Genetics25: 25–2910802651|
|60.||Kanehisa M,, Goto S, (Year: 2000) KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 28: 27–30|
|61.||Ogata H,, Goto S,, Sato K,, Fugibuchi W,, Bono H,, et al. (Year: 1999) KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 27: 29–34|
|62.||Benjamini Y,, Hochberg Y, (Year: 1995) Controlling the false discovery rate – a practical and powerful approach to multiple testing rid C-4219-2008. J. R. Stat. Soc. Ser. B Stat. Methodol. 57: 283–300|
|63.||Bustin SA,, Benes V,, Garson JA,, Hellemans J,, Huggett J,, et al. (Year: 2009) The MIQE guidelines: Minimum information for publication of quantitative real-time PCR experiments. Clinical Chemistry55: 611–62219246619|
|64.||Lagos-Quintana M,, Rauhut R,, Yalcin A,, Meyer J,, Lendeckel W,, et al. (Year: 2002) Identification of tissue-specific microRNAs from mouse. Current Biology: CB12: 735–73912007417|
|65.||Li HY, Xi QY, Xiong YY, Liu XL, Cheng X, et al. . (2012) Identification and comparison of microRNAs from skeletal muscle and adipose tissues from two porcine breeds. Animal Genetics.|
|66.||Lau NC,, Lim LP,, Weinstein EG,, Bartel DP, (Year: 2001) An abundant class of tiny RNAs with probable regulatory roles in caenorhabditis elegans. Science (New York, N.Y.)294: 858–862|
|67.||Zhou Y, Tang X, Song Q, Ji Y, Wang H, et al. . (2012) Identification and characterization of pig embryo MicroRNAs by solexa sequencing. Reproduction in Domestic Animals = Zuchthygiene.|
|68.||Landgraf P,, Rusu M,, Sheridan R,, Sewer A,, Iovino N,, et al. (Year: 2007) A mammalian microRNA expression atlas based on small RNA library sequencing. Cell129: 1401–141417604727|
|69.||Lee LW,, Zhang S,, Etheridge A,, Ma L,, Martin D,, et al. (Year: 2010) Complexity of the microRNA repertoire revealed by next-generation sequencing. RNA (New York, N.Y.)16: 2170–2180|
|70.||Fernandez-Valverde SL,, Taft RJ,, Mattick JS, (Year: 2010) Dynamic isomiR regulation in drosophila development. RNA (New York, N.Y.)16: 1881–1888|
|71.||Li SC,, Liao YL,, Ho MR,, Tsai KW,, Lai CH,, et al. (Year: 2012) miRNA arm selection and isomiR distribution in gastric cancer. BMC Genomics13 Suppl 1S1322369582|
|72.||Dohm JC,, Lottaz C,, Borodina T,, Himmelbauer H, (Year: 2008) Substantial biases in ultra-short read data sets from high-throughput DNA sequencing. Nucleic Acids Research36: e10518660515|
|73.||Taub MA,, Corrada Bravo H,, Irizarry RA, (Year: 2010) Overcoming bias and systematic errors in next generation sequencing data. Genome Medicine2: 8721144010|
|74.||Malone JH,, Oliver B, (Year: 2011) Microarrays, deep sequencing and the true measure of the transcriptome. BMC Biology9: 3421627854|
|75.||Oba S,, Kumano S,, Suzuki E,, Nishimatsu H,, Takahashi M,, et al. (Year: 2010) miR-200b precursor can ameliorate renal tubulointerstitial fibrosis. PloS One5: e1361421049046|
|76.||Li Y, VandenBoom TG,2nd, Kong D, Wang Z, Ali S, et al (Year: 2009) Up-regulation of miR-200 and let-7 by natural agents leads to the reversal of epithelial-to-mesenchymal transition in gemcitabine-resistant pancreatic cancer cells. Cancer Research69: 6704–671219654291|
|77.||Han Y,, Chen J,, Zhao X,, Liang C,, Wang Y,, et al. (Year: 2011) MicroRNA expression signatures of bladder cancer revealed by deep sequencing. PloS One6: e1828621464941|
|78.||Marchini S,, Cavalieri D,, Fruscio R,, Calura E,, Garavaglia D,, et al. (Year: 2011) Association between miR-200c and the survival of patients with stage I epithelial ovarian cancer: A retrospective study of two independent tumour tissue collections. The Lancet Oncology12: 273–28521345725|
|79.||Tang X,, Tang X,, Gal J,, Kyprianou N,, Zhu H,, et al. (Year: 2011) Detection of microRNAs in prostate cancer cells by microRNA array. Methods in Molecular Biology (Clifton, N.J.)732: 69–88|
|80.||Bockmeyer CL, Christgen M, Muller M, Fischer S, Ahrens P, et al. . (2011) MicroRNA profiles of healthy basal and luminal mammary epithelial cells are distinct and reflected in different breast cancer subtypes. Breast Cancer Research and Treatment.|
|81.||Gottardo F,, Liu CG,, Ferracin M,, Calin GA,, Fassan M,, et al. (Year: 2007) Micro-RNA profiling in kidney and bladder cancers. Urologic Oncology25: 387–39217826655|
|82.||White NM,, Bui A,, Mejia-Guerrero S,, Chao J,, Soosaipillai A,, et al. (Year: 2010) Dysregulation of kallikrein-related peptidases in renal cell carcinoma: Potential targets of miRNAs. Biological Chemistry391: 411–42320180642|
[Figure ID: pone-0055402-g001]
Figure 1 Sequence variability in Ssc-miR-23b.
The number of total counts for each sequence is indicated. First sequence belongs to the precursor miRNA. Red bases match with precursor miRNA and belong to mature miRNA. Blue bases match also with precursor miRNA but do not belong to mature miRNA sequence. Sequence corresponding to the annotated miRNA in miRBase database is marked with an asterisk.
[Figure ID: pone-0055402-g002]
Figure 2 Predicted miRNA folding for cluster candidates to novel miRNAs using MFold software1.
1: . Orange line points out the miRNA position into the pre-miRNA structure.
Table 1 Primers used for the RT-qPCR validation design.
|miRNA/Cluster||Forward primer (5′-3′)||Reverse primer (5′-3′)||Primer conc. (nM each)|
|Bta-miR-193b||[gene: AGAACTGGCCCACAAAGTC]||[gene: TCCAGTTTTTTTTTTTTTTTAGCGG]||250|
|Ssc-miR-423-5p||[gene: TGAGGGGCAGAGAGCGA]||[gene: GTCCAGTTTTTTTTTTTTTTTAAAGTC]||250|
|Ssc-miR-29c||[gene: GCAGTAGCACCATTTGAAATC]||[gene: GGTCCAGTTTTTTTTTTTTTTTAACC]||250|
|Ssc-miR-486||[gene: GCAGTCCTGTACTGAGCTG]||[gene: GGTCCAGTTTTTTTTTTTTTTTCTCG]||250|
|Ssc-let-7f||[gene: CGCAGTGAGGTAGTAGATTG]||[gene: AGGTCCAGTTTTTTTTTTTTTTTAACT]||250|
|Cl-2||[gene: CAGCTGCTATGCCAACA]||[gene: CCAGTTTTTTTTTTTTTTTGGCAA]||250|
|Cl-5||[gene: GCTGTAACAGCAACTCCA]||[gene: TCCAGTTTTTTTTTTTTTTTCCAC]||500|
|Cl-15||[gene: GCGACCCACTCTTGGT]||[gene: TCCAGTTTTTTTTTTTTTTTCATGG]||No amplification|
|Cl-16||[gene: CAGTTGGTGACCAGGTG]||[gene: AGTTTTTTTTTTTTTTCTCCCTGA]||No amplification|
|Cl-24||[gene: GCTGCATTTCCTGGCTG]||[gene: GGTCCAGTTTTTTTTTTTTTTTAATTAAG]||No amplification|
|Cl-25||[gene: CAGCTGGTGTTGTGAATC]||[gene: GTTTTTTTTTTTTTTTCGGCCT]||125|
|Cl-29||[gene: GGTTGGTGTACACTGGAA]||[gene: GTCCAGTTTTTTTTTTTTTTTAGCT]||250|
|Cl-38||[gene: GTCTCCGTTTGCCTGT]||[gene: CCAGTTTTTTTTTTTTTTTCAGCA]||250|
Table 2 Summary of sequence alignment with miRBase database1.
|Library||Item||Total||miRBase homology*||Sus scrofa homology*||Other species homology*||No homology in miRBase|
|Iberian (IB)||Total counts||4,022||2,748||1,897||851||1,274|
|Wild Boar (WB)||Total counts||18,391||14,787||7,146||7,641||3,604|
|Landrace (LD)||Total counts||13,704||11,131||4,897||6,234||2,573|
|Large White (LW)||Total counts||23,815||19,728||10,609||9,119||4,087|
|Piétrain (PT)||Total counts||20,215||16,971||9,890||7,081||3,244|
|Meishan (ME)||Total counts||17,407||14,990||5,380||9,610||2,417|
|Vietnamese (VT)||Total counts||17,751||9,239||5,105||4,134||8,512|
1: miRBase database (v18, November 2011).
*Mismatches or gaps have not been allowed into the sequence alignment. Up to 4 mismatches allowed in the ends of the sequences.
Table 3 Most abundant miRNA profile in pig kidney (CN>350).
|miRNA name||Total counts (%)||IsomiRs||IB||WB||LD||LW||PT||ME||VT|
IB: Iberian breed, WB: European Wild Boar, LD: Landrace breed, LW: Large White breed, PT: Piétrain breed, ME: Meishan breed, VT: Vietnamese breed.
miRNA name represents the most expressed miRNA sequence in the cluster.
Underlined miRNAs correspond to miRNAs that had not been previously described.
φ: miRNA clusters whose predominant isomiR is different from that described in miRBase database (v18, November 2011).
Table 4 Most expressed miRNAs diverging between the reference miRNA sequence and the miRBase1 described miRNA sequence.
|miRNA name||Total counts||IsomiRs||Most expressed isomiR sequence (5′-3′)||Counts||Described miRNA sequence in miRBase database (5′-3′)||Counts|
|Ssc-miR-23a||5,156||51||[gene: ATCACATTGCCAGGGATTTCCA]||3,116||[gene: ATCACATTGCCAGGGATTTCC]||622|
|Ssc-miR-99a||3,781||32||[gene: AACCCGTAGATCCGATCTTGTGA]||1,124||[gene: AACCCGTAGATCCGATCTTGTG]||1,085|
|Ssc-miR-30d||1,977||31||[gene: TGTAAACATCCCCGACTGGAAGC]||711||[gene: TGTAAACATCCCCGACTGGAAGCT]||190|
|Ssc-miR-125a||1,369||29||[gene: TCCCTGAGACCCTTTAACCTGT]||479||[gene: TCCCTGAGACCCTTTAACCTGTG]||326|
|Bta-miR-193b||473||26||[gene: AACTGGCCCACAAAGTCCCGCT]||196||[gene: AACTGGCCCACAAAGTCCCGCTTT]||21|
|Ssc-miR-30e-5p||461||14||[gene: TGTAAACATCCTTGACTGGAAGC]||268||[gene: TGTAAACATCCTTGACTGGAAGCT]||74|
|Ssc-miR-139-5p||329||8||[gene: TCTACAGTGCACGTGTCTCCAGT]||216||[gene: TCTACAGTGCACGTGTCTCCAG]||85|
|Ssc-miR-30a-5p||314||10||[gene: TGTAAACATCCTCGACTGGAAGC]||198||[gene: TGTAAACATCCTCGACTGGAAGCT]||15|
|Ssc-miR-362||303||13||[gene: AATCCTTGGAACCTAGGTGTGAGT]||133||[gene: AATCCTTGGAACCTAGGTGTGAGTG]||0|
|Hsa-miR-29c-5p||260||20||[gene: GACCGATTTCTCCTGGTGTTCA]||37||[gene: TGACCGATTTCTCCTGGTGTTC]||22|
|Ssc-miR-374a||253||6||[gene: TTATAATACAACCTGATAAGTGT]||117||[gene: TTATAATACAACCTGATAAGTG]||13|
|Ssc-miR-145||222||12||[gene: GTCCAGTTTTCCCAGGAATCCCT]||91||[gene: GTCCAGTTTTCCCAGGAATCCCTT]||3|
|Hsa-miR-324-3p||211||22||[gene: ACTGCCCCAGGTGCTGCTGGT]||45||[gene: ACTGCCCCAGGTGCTGCTGG]||0|
|Ssc-miR-21||210||4||[gene: TAGCTTATCAGACTGATGTTGAC]||105||[gene: TAGCTTATCAGACTGATGTTGA]||77|
|Hsa-miR-874||195||8||[gene: CTGCCCTGGCCCGAGGGACCGAC]||56||[gene: CTGCCCTGGCCCGAGGGACCGA]||42|
|Ssc-miR-218b||185||7||[gene: TTGTGCTTGATCTAACCATGT]||56||[gene: TTGTGCTTGATCTAACCATGTG]||27|
|Ssc-miR-191||183||8||[gene: CAACGGAATCCCAAAAGCAGCT]||63||[gene: CAACGGAATCCCAAAAGCAGCTG]||45|
|Hsa-miR-150-5p||174||12||[gene: GTTCTCCCAACCCTTGTACCAGT]||47||[gene: TCTCCCAACCCTTGTACCAGTG]||14|
|Ssc-miR-22-3p||106||5||[gene: AAGCTGCCAGTTGAAGAAC]||43||[gene: AAGCTGCCAGTTGAAGAACTGT]||28|
1: miRBase database (v18, November 2011).
Bta: Bos taurus, Hsa: Homo sapiens, Ssc: Sus scrofa.
Marked in bold the nucleotide variation between pair sequences.
Table 5 Selected miRNA and putative novel miRNAs (Clusters) for further validation through RT-qPCR.
|miRNA/Cluster||Total counts||IsomiRs||Sequence (5′-3′)|
|Cl-2 (Hsa-miR-31-3p)1||4||1||[gene: CTGCTATGCCAACATATTGCCA]|
|Cl-5 (Hsa-miR-194-5p)1||12||2||[gene: CTGTAACAGCAACTCCATGTGGAA]|
|Cl-15 (Hsa-miR-551a)1||9||2||[gene: GCGACCCACTCTTGGTTTCCATG]|
|Cl-25 (Hsa-miR-138-5p)1||7||3||[gene: CAGCTGGTGTTGTGAATCAGGCCG]|
|Cl-38 (Bta-miR-1468)1||21||3||[gene: TCTCCGTTTGCCTGTTTTGCTGA]|
1: miRBase (v18) homology (e-value ≤1e-03). Internal mismatches were accepted.
Bta: Bos taurus, Hsa: Homo sapiens, Ssc: Sus scrofa.
Table 6 Fold Change (FC) comparison from HTS differential expression study and RT-qPCR data.
|miRNA/Cluster||EU2vs. EA3||EU vs. AS4||EA vs. AS||ANOVA group factor p-value||EU vs. EA||EU vs. AS||EA vs. AS||ANOVA group factor p-value|
1: miRBase (v18) homology (e-value≤1e-03). Internal mismatches were accepted.
2: EU: European breeds. 3: EA: European commercial breeds. 4: AS: Asian breeds.
For HTS data, fold changes from sequence counting between breed groups were calculated from normalised data in counts per thousand for each library and averaged per groups. Positive and negative signs indicate that the level of gene expression is higher for the first or the second group of the test, respectively. Analysis of variance (ANOVA) including breed as fixed factor was performed. Significance was set at P<0.05. Fold changes for clusters could not be calculated due to their low total counts.
For qPCR data expression study, the quantity obtained of each miRNA in each sample was normalised by the Normalization Factor and corrected in relation to the lowest normalised value. Analysis of variance (ANOVA) including the RT and breed as fixed factors was performed and fold change from Least Squares Means (LSM) between breed groups were calculated. Significance was set at P<0.05. Scheffe test determined whether there was significant differential expression between breed groups (§: suggestive p-value<0.1, *: significant with p-value<0.05, **: significant with p-value<0.01, ***: significant with p-value<0.001).
Previous Document: Variable Strength of Forest Stand Attributes and Weather Conditions on the Questing Activity of Ixod...
Next Document: Molecular vibration-sensing component in human olfaction.