The microRNA expression profile in porcine skeletal muscle is changed by constant heat stress

Heat stress has profound effects on animal performance and muscle function, and microRNAs (miRNAs) play a critical role in muscle development and stress responses. To characterize the changes in miRNAs in skeletal muscle responding to heat stress, the miRNA expression profiles of longissimus dorsi muscles of pigs raised under constant heat stress (30 °C; n = 8) or control temperature (22 °C; n = 8) for 21 days were analyzed by Illumina deep sequencing. A total of 58 differentially expressed miRNAs were identified with 30 down‐regulated and 28 up‐regulated, and 63 differentially expressed target genes were predicted by miRNA–mRNA joint analysis. GO and KEGG analyses showed that the genes regulated by differentially expressed miRNAs were enriched in glucose metabolism, cytoskeletal structure and function and stress response. Real‐time PCR showed that the mRNA levels of PDK4, HSP90 and DES were significantly increased, whereas those of SCD and LDHA significantly decreased by heat exposure. The protein levels of CALM1, DES and HIF1α were also significantly increased by constant heat. These results demonstrated that the change in miRNA expression in porcine longissimus dorsi muscle underlies the changes in muscle structure and metabolism in porcine skeletal muscle affected by constant heat stress.     

雌性藏猪不同生长发育阶段肌肉生长和脂肪沉积相关基因的表达模式

本研究利用包含140个与猪肌肉生长和脂肪沉积密切相关基因的Oligo功能分类基因芯片检测了藏猪在2、4、6和8月龄间背最长肌中这些基因的表达量变化,并在2月龄时与脂肪型的太湖猪和瘦肉型的长白猪进行比较.ANOVA分析结果表明:2-8月龄间藏猪分别有10和 7个基因的表达差异达极显著(P<0.01)和显著水平(P<0.05);2月龄时藏猪体重极显著低于长白猪(P<0.01)和显著低于太湖猪(P<0.05),而藏猪肌纤维面积却为最大,但品种间差异未达显著水平(P>0.05);2月龄时3个品种间分别有15和13个基因的表达差异达极显著 (P<0.01)和显著水平(P<0.05).STEM聚类分析结果表明:直线下降和上升是藏猪在2 -8月龄间最具代表性的基因表达模式(P<0.01).另外,5个差异表达基因的荧光定量RT -PC R验证结果与基因芯片结果的Person相关系数平均高达0.856±0.109.提示:藏猪在2-8月龄间骨骼肌生长发育强度较肌内脂肪合成沉积占优势,2月龄时藏猪脂肪酸合成相关基因的表达水平较其他两品种低,而脂肪酸β氧化和肌纤维生长相关基因的表达水平较高,与其在高原独特的自然生态环境和全放牧散养的饲喂方式下长期形成的品种特性相符 [动物学报 54(3):442-452,2008].     

Integrated analysis of miRNA and mRNA expression profiles in response to gut microbiota depletion in the abdomens of female Bactrocera dorsalis

Insect gut microbiota has been reported to participate in regulating host multiple biological processes including metabolism and reproduction. However, the corresponding molecular mechanisms remain largely unknown. Recent studies suggest that microRNAs (miRNAs) are involved in complex interactions between the gut microbiota and the host. Here, we used next-generation sequencing technology to characterize miRNA and mRNA expression profiles and construct the miRNA–gene regulatory network in response to gut microbiota depletion in the abdomens of female Bactrocera dorsalis. A total of 3016 differentially expressed genes (DEGs) and 18 differentially expressed miRNAs (DEMs) were identified. Based on the integrated analysis of miRNA and mRNA sequencing data, 229 negatively correlated miRNA–gene pairs were identified from the miRNA–mRNA network. Gene ontology enrichment analysis indicated that DEMs could target several genes involved in the metabolic process, oxidation–reduction process, oogenesis, and insulin signaling pathway. Finally, real-time quantitative polymerase chain reaction further verified the accuracy of RNA sequencing results. In conclusion, our study provides the profiles of miRNA and mRNA expressions under antibiotics treatment and provides an insight into the roles of miRNAs and their target genes in the interaction between the gut microbiota and its host.     

Expression profiling analysis for genes related to meat quality and carcass traits during postnatal development of backfat in two pig breeds

The competitive equilibrium of fatty acid biosynthesis and oxidation in vivo determines porcine subcutaneous fat thickness (SFT) and intramuscular fat (IMF) content. Obese and lean-type pig breeds show obvious differences in adipose deposition; however, the molecular mechanism underlying this phenotypic variation remains unclear. We used pathway-focused oligo microarray studies to examine the expression changes of 140 genes associated with meat quality and carcass traits in backfat at five growth stages (1–5 months) of Landrace (a leaner, Western breed) and Taihu pigs (a fatty, indigenous, Chinese breed). Variance analysis (ANOVA) revealed that differences in the expression of 25 genes in Landrace pigs were significant (FDR adjusted permutation, P<0.05) among 5 growth stages. Gene class test (GCT) indicated that a gene-group was very significant between 2 pig breeds across 5 growth stages (P _ErmineJ<0.01), which consisted of 23 genes encoding enzymes and regulatory proteins associated with lipid and steroid metabolism. These findings suggest that the distinct differences in fat deposition ability between Landrace and Taihu pigs may closely correlate with the expression changes of these genes. Clustering analysis revealed a very high level of significance (FDR adjusted, P<0.01) for 2 gene expression patterns in Landrace pigs and a high level of significance (FDR adjusted, P<0.05) for 2 gene expression patterns in Taihu pigs. Also, expression patterns of genes were more diversified in Taihu pigs than those in Landrace pigs, which suggests that the regulatory mechanism of micro-effect polygenes in adipocytes may be more complex in Taihu pigs than in Landrace pigs. Based on a dynamic Bayesian network (DBN) model, gene regulatory networks (GRNs) were reconstructed from time-series data for each pig breed. These two GRNs initially revealed the distinct differences in physiological and biochemical aspects of adipose metabolism between the two pig breeds; from these results, some potential key genes could be identified. Quantitative, real-time RT-PCR (QRT-PCR) was used to verify the microarray data for five modulated genes, and a good correlation between the two measures of expression was observed for both 2 pig breeds at different growth stages (R=0.874±0.071). These results highlight some possible candidate genes for porcine fat characteristics and provide some data on which to base further study of the molecular basis of adipose metabolism.     

MicroRNA expression and regulation in human, chimpanzee, and macaque brains

Among other factors, changes in gene expression on the human evolutionary lineage have been suggested to play an important role in the establishment of human-specific phenotypes. However, the molecular mechanisms underlying these expression changes are largely unknown. Here, we have explored the role of microRNA (miRNA) in the regulation of gene expression divergence among adult humans, chimpanzees, and rhesus macaques, in two brain regions: prefrontal cortex and cerebellum. Using a combination of high-throughput sequencing, miRNA microarrays, and Q-PCR, we have shown that up to 11% of the 325 expressed miRNA diverged significantly between humans and chimpanzees and up to 31% between humans and macaques. Measuring mRNA and protein expression in human and chimpanzee brains, we found a significant inverse relationship between the miRNA and the target genes expression divergence, explaining 2%-4% of mRNA and 4%-6% of protein expression differences. Notably, miRNA showing human-specific expression localize in neurons and target genes that are involved in neural functions. Enrichment in neural functions, as well as miRNA-driven regulation on the human evolutionary lineage, was further confirmed by experimental validation of predicted miRNA targets in two neuroblastoma cell lines. Finally, we identified a signature of positive selection in the upstream region of one of the five miRNA with human-specific expression, miR-34c-5p. This suggests that miR-34c-5p expression change took place after the split of the human and the Neanderthal lineages and had adaptive significance. Taken together these results indicate that changes in miRNA expression might have contributed to evolution of human cognitive functions.     

A Comparison of Brain Gene Expression Levels in Domesticated and Wild Animals

Domestication has led to similar changes in morphology and behavior in several animal species, raising the question whether similarities between different domestication events also exist at the molecular level. We used mRNA sequencing to analyze genome-wide gene expression patterns in brain frontal cortex in three pairs of domesticated and wild species (dogs and wolves, pigs and wild boars, and domesticated and wild rabbits). We compared the expression differences with those between domesticated guinea pigs and a distant wild relative (Cavia aperea) as well as between two lines of rats selected for tameness or aggression towards humans. There were few gene expression differences between domesticated and wild dogs, pigs, and rabbits (30–75 genes (less than 1%) of expressed genes were differentially expressed), while guinea pigs and C. aperea differed more strongly. Almost no overlap was found between the genes with differential expression in the different domestication events. In addition, joint analyses of all domesticated and wild samples provided only suggestive evidence for the existence of a small group of genes that changed their expression in a similar fashion in different domesticated species. The most extreme of these shared expression changes include up-regulation in domesticates of SOX6 and PROM1, two modulators of brain development. There was almost no overlap between gene expression in domesticated animals and the tame and aggressive rats. However, two of the genes with the strongest expression differences between the rats (DLL3 and DHDH) were located in a genomic region associated with tameness and aggression, suggesting a role in influencing tameness. In summary, the majority of brain gene expression changes in domesticated animals are specific to the given domestication event, suggesting that the causative variants of behavioral domestication traits may likewise be different.     

Genome‐wide analysis reveals artificial selection on coat colour and reproductive traits in Chinese domestic pigs

Pigs from Asia and Europe were independently domesticated from c. 9000 years ago. During this period, strong artificial selection has led to dramatic phenotypic changes in domestic pigs. However, the genetic basis underlying these morphological and behavioural adaptations is relatively unknown, particularly for indigenous Chinese pigs. Here, we performed a genome‐wide analysis to screen 196 regions with selective sweep signals in Tongcheng pigs, which are a typical indigenous Chinese breed. Genes located in these regions have been found to be involved in lipid metabolism, melanocyte differentiation, neural development and other biological processes, which coincide with the evolutionary phenotypic changes in this breed. A synonymous substitution, c.669T>C, in ESR1, which colocalizes with a major quantitative trait locus for litter size, shows extreme differences in allele frequency between Tongcheng pigs and wild boars. Notably, the variant C allele in this locus exhibits high allele frequency in most Chinese populations, suggesting a consequence of positive selection. Five genes (PRM1, PRM2, TNP2, GPR149 and JMJD1C) related to reproductive traits were found to have high haplotype similarity in Chinese breeds. Two selected genes, MITF and EDNRB, are implied to shape the two‐end black colour trait in Tongcheng pig. Subsequent SNP microarray studies of five Chinese white‐spotted breeds displayed a concordant signature at both loci, suggesting that these two genes are responsible for colour variations in Chinese breeds. Utilizing massively parallel sequencing, we characterized the candidate sites that adapt to artificial and environmental selections during the Chinese pig domestication. This study provides fundamental proof for further research on the evolutionary adaptation of Chinese pigs.     

MicroRNA Buffering and Altered Variance of Gene Expression in Response to Salmonella Infection

One potential role of miRNAs is to buffer variation in gene expression, although conflicting results have been reported. To investigate the buffering role of miRNAs in response to Salmonella infection in pigs, we sequenced miRNA and mRNA in whole blood from 15 pig samples before and after Salmonella challenge. By analyzing inter-individual variation in gene expression patterns, we found that for moderately and lowly expressed genes, putative miRNA targets showed significantly lower expression variance compared with non-miRNA-targets. Expression variance between highly expressed miRNA targets and non-miRNA-targets was not significantly different. Further, miRNA targets demonstrated significantly reduced variance after challenge whereas non-miRNA-targets did not. RNA binding proteins (RBPs) are significantly enriched among the miRNA targets with dramatically reduced variance of expression after Salmonella challenge. Moreover, we found evidence that targets of young (less-conserved) miRNAs showed lower expression variance compared with targets of old (evolutionarily conserved) miRNAs. These findings point to the importance of a buffering effect of miRNAs for relatively lowly expressed genes, and suggest that the reduced expression variation of RBPs may play an important role in response to Salmonella infection.     

Population history and genomic signatures for high-altitude adaptation in Tibetan pigs

Background

The Tibetan pig is one of domestic animals indigenous to the Qinghai-Tibet Plateau. Several geographically isolated pig populations are distributed throughout the Plateau. It remained an open question if these populations have experienced different demographic histories and have evolved independent adaptive loci for the harsh environment of the Plateau. To address these questions, we herein investigated ~ 40,000 genetic variants across the pig genome in a broad panel of 678 individuals from 5 Tibetan geographic populations and 34 lowland breeds.

Results

Using a series of population genetic analyses, we show that Tibetan pig populations have marked genetic differentiations. Tibetan pigs appear to be 3 independent populations corresponding to the Tibetan, Gansu and Sichuan & Yunnan locations. Each population is more genetically similar to its geographic neighbors than to any of the other Tibetan populations. By applying a locus-specific branch length test, we identified both population-specific and -shared candidate genes under selection in Tibetan pigs. These genes, such as PLA2G12A, RGCC, C9ORF3, GRIN2B, GRID1 and EPAS1, are involved in high-altitude physiology including angiogenesis, pulmonary hypertension, oxygen intake, defense response and erythropoiesis. A majority of these genes have not been implicated in previous studies of highlanders and high-altitude animals.

Conclusion

Tibetan pig populations have experienced substantial genetic differentiation. Historically, Tibetan pigs likely had admixture with neighboring lowland breeds. During the long history of colonization in the Plateau, Tibetan pigs have developed a complex biological adaptation mechanism that could be different from that of Tibetans and other animals. Different Tibetan pig populations appear to have both distinct and convergent adaptive loci for the harsh environment of the Plateau.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-834) contains supplementary material, which is available to authorized users.     

Comparison of stomach microRNA transcriptomes of Tibetan and Yorkshire pigs by deep sequencing

MiRNAs regulate the expression of target genes in diverse cellular processes and hence play important roles in different physiological processes, yet little is known about the stomach microRNAome (miRNAome) of the Tibetan pig. The objective of this experiment was to investigate differentially expressed stomach miRNAs participating in digestion. Firstly, we isolated total RNA by Trizol reagent from three Tibetan and three Yorkshire purebred pigs stomach samples at 90-day-old. Secondly, a comprehensive analysis of Tibetan and Yorkshire pig stomach miRNAomes was performed by small RNA sequencing in the Illumina HiSeq 2000 system. Finally, SYBR Green Real-time RT-PCR was performed to validate the differentially expressed miRNAs. We identified 318 unique miRNAs, 260 were co-expressed in both libraries, 17 and 31 miRNAs were specifically expressed in Tibetan and Yorkshire pigs respectively. Fifty six differentially expressed miRNAs were identified by the identifying differentially expressed genes 6 (IDEG6). Kyoto encyclopedia of genes and genomes analysis revealed that some of the differentially expressed miRNAs were associated with protein and fat digestion. Two differentially expressed miRNAs (miR-214-3p and ssc-un39) participating in the digestion of lipid were identified. Additionally, qRT-PCR results suggested that a higher expression of miR-214-3p in the Tibetan pig stomach could lead to relatively lower expression of calcium-dependent phospholipase A2, which is an enzyme important for the digestion of glycerol phospholipid. This study has delineated the different stomach miRNAs expression patterns of Tibetan and Yorkshire pigs, which would help explain the regulatory mechanisms of miRNAs in digestion of Tibetan pigs, and contribute to utilize a the unique digestion merits of Tibetan pig in future porcine hybridization breeding.     

Expression profiling analysis for genes related to meat quality and carcass traits during postnatal development of backfat in two pig breeds

The competitive equilibrium of fatty acid biosynthesis and oxidation in vivo determines porcine sub-cutaneous fat thickness(SFT) and intramuscular fat(IMF) content.Obese and lean-type pig breeds show obvious differences in adipose deposition;however, the molecular mechanism underlying this phenotypic variation remains unclear.We used pathway-focused oligo microarray studies to examine the expression changes of 140 genes associated with meat quality and carcass traits in backfat at five growth stages(1―5 months) of Landrace(a leaner, Western breed) and Taihu pigs(a fatty, indigenous, Chinese breed).Variance analysis(ANOVA) revealed that differences in the expression of 25 genes in Landrace pigs were significant(FDR adjusted permutation, P<0.05) among 5 growth stages.Gene class test(GCT) indicated that a gene-group was very significant between 2 pig breeds across 5 growth stages(PErmineJ<0.01), which consisted of 23 genes encoding enzymes and regulatory proteins associ-ated with lipid and steroid metabolism.These findings suggest that the distinct differences in fat deposition ability between Landrace and Taihu pigs may closely correlate with the expression changes of these genes.Clustering analysis revealed a very high level of significance(FDR adjusted, P<0.01) for 2 gene expression patterns in Landrace pigs and a high level of significance(FDR adjusted, P<0.05) for 2 gene expression patterns in Taihu pigs.Also, expression patterns of genes were more diversified in Taihu pigs than those in Landrace pigs, which suggests that the regulatory mechanism of micro-effect polygenes in adipocytes may be more complex in Taihu pigs than in Landrace pigs.Based on a dy-namic Bayesian network(DBN) model, gene regulatory networks(GRNs) were reconstructed from time-series data for each pig breed.These two GRNs initially revealed the distinct differences in physiological and biochemical aspects of adipose metabolism between the two pig breeds;from these results, some potential key genes could be identified.Quantitative, real-time RT-PCR(QRT-PCR) was used to verify the microarray data for five modulated genes, and a good correlation between the two measures of expression was observed for both 2 pig breeds at different growth stages(R=0.874±0.071).These results highlight some possible candidate genes for porcine fat characteristics and provide some data on which to base further study of the molecular basis of adipose metabolism.     

Comparison of microRNAs in the intramuscular adipose tissue from Jinhua and Landrace pigs

MicroRNA (miRNA) is critically involved in lipogenesis occurring in various body parts of humans and animals. In this study, to further investigate the role and distribution of miRNA in porcine intramuscular adipose tissue, small RNAs were extracted from Jinhua and Landrace pigs to identify the expression profiles of miRNAs. miRNA expression profiles revealed that 558 miRNAs including 287 known and 271 novel miRNAs were identified, and 220 of them showed differential expression in the pigs. Gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis suggested that the target genes of the differentially expressed miRNAs were involved in fatty metabolism. In conclusion, the current study reveals the active participation of miRNAs in the regulation of adipogenesis in the intramuscular adipose tissue of Jinhua and Landrace pigs.     

Lipid metabolism related gene-expression profiling in liver,skeletal muscle and adipose tissue in crossbred Duroc and Pietrain Pigs

Body-weight differences in animals may be ascribed to genetic and environmental factors. Here we utilized two divergent porcine genotypes, the highly muscled, leaner Pietrian × Yorkshire pigs and less muscled, fatter Duroc × Yorkshire growing pigs (75–110 kg), to examine the role of genetic background on expression of genes associated with anabolic (Fatty acid synthase, FAS; glucose transporter 4, GLUT-4; stearoyl CoA desaturase, SCD; Sterol regulatory binding protein-1, SREBP-1; leptin) and catabolic lipid metabolism (Carnitine palmitoyltransferase-1B, CPT-1B; acyl-CoA dehydrogenase, ACDH) in adipose tissue (AT), liver (L) and skeletal muscle (SKM). Pietrain pigs had lower mRNA abundance for FAS, SREBP-1, SCD and leptin in AT and L, but higher mRNA abundance for L ACDH and SKM ACDH and CPT-1B than Durocs. Duroc pigs exhibited higher expression of FAS, SREBP-1, SCD, leptin in AT and FAS in L and lower expression of ACDH and CPT-1B in L SKM. GLUT-4 expression did not differ in SKM between the two genotypes. Feeding of a beta adrenergic agonist (Paylean) for 52 days lowered expression of lipid anabolic and enhanced lipid catabolic genes expressions similarly in both genotypes. Overall, the lipid metabolism genes differential expression patterns documented here showed that in Pietrain pigs mRNA abundances of synthesis genes were lower and of catabolic genes were higher than in Duroc pigs.