Identification of target genes and pathways associated with chicken microRNA miR-143

miRNA profiling of chicken follicles during follicular development

MicroRNAs (miRNAs) play a crucial role as transcription regulators in various aspects of follicular development, including steroidogenesis, ovulation, apoptosis, and gene regulation in poultry. However, there is a paucity of studies examining the specific impact of miRNAs on ovarian granulosa cells (GCs) across multiple grades in laying hens. Consequently, this study aims to investigate the roles of miRNAs in chicken GCs. By constructing miRNA expression profiles of GCs at 10 different time points, encompassing 4 pre-hierarchical, 5 preovulatory, and 1 postovulatory follicles stage, we identified highly expressed miRNAs involved in GC differentiation (miR-148a-3p, miR-143-3p), apoptosis (let7 family, miR-363-3p, miR-30c-5p, etc.), and autophagy (miR-128-3p, miR-21-5p). Furthermore, we discovered 48 developmentally dynamic miRNAs (DDMs) that target 295 dynamic differentially expressed genes (DDGs) associated with follicular development and selection (such as oocyte meiosis, progesterone-mediated oocyte maturation, Wnt signaling pathway, TGF-β signaling pathway) as well as follicular regression (including autophagy and cellular senescence). These findings contribute to a more comprehensive understanding of the intricate mechanisms underlying follicle recruitment, selection, and degeneration, aiming to enhance poultry's reproductive capacity.

Relationship between Functional miR-143/145 Cluster Variants and Susceptibility to Type 2 Diabetes Mellitus: A Preliminary Case-Control Study and Bioinformatics Analyses

Purpose: To investigate the link between two variants (rs4705342 and rs4705343) in the promoter of the miR-143/145 cluster with Type 2 diabetes mellitus (T2DM) risk. Methods:A total of 1200 subjects were genotyped using the ARMS-PCR method. Results: The rs4705342 variant enhanced the risk of T2DM under codominant CC (OR = 3.24; 95% CI: 1.89-5.60), recessive TT+TC (OR = 3.02; 95% CI: 1.77-5.17), and dominant TC+CC (OR = 1.35; 95% CI: 1.08-1.71) genetic models. Individuals carrying the C allele of rs4705342 conferred a 1.43 fold increased risk of T2DM. As regards rs4705343, decreased risk of T2DM was observed under codominant TC (OR = 0.53; 95% CI: 0.42-0.67), over-dominant TT+CC (OR = 0.51; 95% CI: 0.40-0.64), and dominant TC+CC (OR = 0.59; 95% CI: 0.48-0.75) models. Haplotype analysis of the variants showed a 1.941-fold increased risk of T2DM regarding the C T combination. Significant associations were noticed between different haplotypes and lipid indices of T2DM patients. There were no notable changes in p-values after adjustment for BMI. Computational analysis revealed that miR143 and/or miR145 target important genes involved in glucose and lipid metabolism. Conclusions: Functional miR-143/145 variants might influence the risk of T2DM. Hence, clarifying the precise regulatory mechanisms of gene expression in the development of T2DM will significantly guide researchers to find a novel target for therapeutic intervention.

Functional miR143/145 Cluster Variants and Haplotypes Are Associated with Chronic Kidney Disease: a Preliminary Case-Control Study and Computational Analyses

MiR-143/145 cluster is a novel transcriptional target of many signaling pathways, with variations within this cluster contributed to the risk of multiple diseases. To date, no data regarding the link between miR143/145 cluster polymorphisms and the risk of developing chronic kidney disease (CKD) has been reported. Hence, we aimed to examine such association in a population of Iranian ancestry. In this preliminary study, 276 CKD patients and 300 unrelated age and sex-matched healthy controls were recruited. Genotyping was performed by PCR-RFLP and allele-specific-PCR methods. Computational analyses were performed to predict the potential effects of the variants. Our findings indicated that rs41291957, rs12659504, and rs353292 polymorphisms were positively associated with CKD, while rs4705342 and rs4705343 polymorphisms demonstrated a significant negative association with the disease. Moreover, a significant association was observed between CC + TC and TT genotypes and CKD stages. We found that AACTT, AATTC, AATTT, GATTC, GATTT, and GGCTT haplotypes significantly enhanced the risk of CKD compared with the G_rs41291957AA_rs12659504C_rs353292T_rs4705342T_rs4705343 haplotype. Computational analysis showed that rs353292, rs4705342, and rs4705343 might alter the binding of the transcription factors in this gene cluster. We found that miR-143/145 cluster polymorphisms were associated with CKD risk in a sample of the Iranian population. Replicated studies on different ethnicities are necessary to investigate the association between these promoter variants and clinical outcomes.

Noncoding RNAs: the shot callers in tumor immune escape

Immunotherapy, designed to exploit the functions of the host immune system against tumors, has shown considerable potential against several malignancies. However, the utility of immunotherapy is heavily limited due to the low response rate and various side effects in the clinical setting. Immune escape of tumor cells may be a critical reason for such low response rates. Noncoding RNAs (ncRNAs) have been identified as key regulatory factors in tumors and the immune system. Consequently, ncRNAs show promise as targets to improve the efficacy of immunotherapy in tumors. However, the relationship between ncRNAs and tumor immune escape (TIE) has not yet been comprehensively summarized. In this review, we provide a detailed account of the current knowledge on ncRNAs associated with TIE and their potential roles in tumor growth and survival mechanisms. This review bridges the gap between ncRNAs and TIE and broadens our understanding of their relationship, providing new insights and strategies to improve immunotherapy response rates by specifically targeting the ncRNAs involved in TIE.

The role of microRNAs in the regulation of insulin signaling pathway with respect to metabolic and mitogenic cascades: A review

Insulin resistance (IR) is a shared pathological condition among type 2 diabetes, obesity, cardiovascular disease, and other metabolic disorders. It is growing significantly all over the world and consequently, a substantial effort is needed for developing the potential novel diagnostics and therapeutics. An insulin signaling pathway is tightly modulated by different mechanisms including the epigenetic modifications. Today, a deal of great attention has been shifted towards the regulatory role of noncoding RNAs on target proteins of the insulin signaling pathway. Noncoding RNAs are a major area of the epigenetics which control gene expression at the posttranscriptional levels and include a large class of microRNAs (miRNAs). With this in view, many studies have implicated the mediatory effects of miRNAs on the downstream metabolic and mitogenic proteins of the insulin signaling pathway. Since providing new biomarkers for the early diagnosis of IR and related metabolic traits are very significant, we intended to review the possible role of miRNAs in the regulation of the insulin signaling pathway, with a primary focus on the downstream target proteins of the metabolic and mitogenic cascades.

MicroRNA expression profiling of porcine mammary epithelial cells after challenge with Escherichia coli in vitro

Background

Coliform mastitis is a symptom of postpartum dysgalactia syndrome (PDS), a multifactorial infectious disease of sows. Our previous study showed gene expression profile change after bacterial challenge of porcine mammary epithelial cells (PMECs). These mRNA expression changes may be regulated through microRNAs (miRNAs) which play critical roles in biological processes. Therefore, miRNA expression profile was investigated in PMECs.

Results

PMECs were isolated from three lactating sows and challenged with heat-inactivated potential mastitis-causing pathogen Escherichia coli (E. coli) for 3 h and 24 h, in vitro. At 3 h post-challenge with E. coli, target gene prediction identified a critical role of miRNAs in regulation of host immune responses and homeostasis of PMECs mediated by affecting pathways including cytokine binding (miR-202, miR-3277, miR-4903); IL-10/PPAR signaling (miR-3277, miR-4317, miR-548); and NF-ĸB/TNFR2 signaling (miR-202, miR-2262, miR-885-3p). Target genes of miRNAs in PMECs at 24 h were significantly enriched in pathways associated with interferon signaling (miR-210, miR-23a, miR-1736) and protein ubiquitination (miR-125, miR-128, miR-1280).

Conclusions

This study provides first large-scale miRNA expression profiles and their predicted target genes in PMECs after contact with a potential mastitis-causing E. coli strain. Both, highly conserved miRNAs known from other species as well as novel miRNAs were identified in PMECs, representing candidate predictive biomarkers for PDS. Time-dependent pathogen clearance suggests an important role of PMECs in inflammatory response of the first cellular barrier of the porcine mammary gland.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4070-2) contains supplementary material, which is available to authorized users.

Enrichment and verification of differentially expressed miRNAs in bursa of Fabricius in two breeds of duck

Objective

The bursa of Fabricius (BF) is a central humoral immune organ belonging specifically to avians. Recent studies had suggested that miRNAs were active regulators involved in the immune processes. This study was to investigate the possible differences of the BF at miRNA level between two genetically disparate duck breeds.

Methods

Using Illumina next-generation sequencing, the miRNAs libraries of ducks were established.

Results

The results showed that there were 66 differentially expressed miRNAs and 28 novel miRNAs in bursa. A set of abundant miRNAs (i.e., let-7, miR-146a-5p, miR-21-5p, miR-17~92) which are involved in immunity and disease were detected and the predicted target genes of the novel miRNAs were associated with duck high anti-adversity ability. By gene ontology analysis and enriching KEGG pathway, the targets of differential expressed miRNAs were mainly involved in immunity and disease, supporting that there were differences in the BF immune functions between the two duck breeds. In addition, the metabolic pathway had the maximum enriched target genes and some enriched pathways that were related to cell cycle, protein synthesis, cell proliferation and apoptosis. It indicted that the difference of metabolism may be one of the reasons leading the immune difference between the BF of two duck breeds.

Conclusion

This data lists the main differences in the BF at miRNAs level between two genetically disparate duck breeds and lays a foundation to carry out molecular assisted breeding of poultry in the future.

An Integrated Analysis of MicroRNA and mRNA Expression Profiles to Identify RNA Expression Signatures in Lambskin Hair Follicles in Hu Sheep

Wave patterns in lambskin hair follicles are an important factor determining the quality of sheep’s wool. Hair follicles in lambskin from Hu sheep, a breed unique to China, have 3 types of waves, designated as large, medium, and small. The quality of wool from small wave follicles is excellent, while the quality of large waves is considered poor. Because no molecular and biological studies on hair follicles of these sheep have been conducted to date, the molecular mechanisms underlying the formation of different wave patterns is currently unknown. The aim of this article was to screen the candidate microRNAs (miRNA) and genes for the development of hair follicles in Hu sheep. Two-day-old Hu lambs were selected from full-sib individuals that showed large, medium, and small waves. Integrated analysis of microRNA and mRNA expression profiles employed high-throughout sequencing technology. Approximately 13, 24, and 18 differentially expressed miRNAs were found between small and large waves, small and medium waves, and medium and large waves, respectively. A total of 54, 190, and 81 differentially expressed genes were found between small and large waves, small and medium waves, and medium and large waves, respectively, by RNA sequencing (RNA-seq) analysis. Differentially expressed genes were classified using gene ontology and pathway analyses. They were found to be mainly involved in cell differentiation, proliferation, apoptosis, growth, immune response, and ion transport, and were associated with MAPK and the Notch signaling pathway. Reverse transcription-polymerase chain reaction (RT-PCR) analyses of differentially-expressed miRNA and genes were consistent with sequencing results. Integrated analysis of miRNA and mRNA expression indicated that, compared to small waves, large waves included 4 downregulated miRNAs that had regulatory effects on 8 upregulated genes and 3 upregulated miRNAs, which in turn influenced 13 downregulated genes. Compared to small waves, medium waves included 13 downregulated miRNAs that had regulatory effects on 64 upregulated genes and 4 upregulated miRNAs, which in turn had regulatory effects on 22 downregulated genes. Compared to medium waves, large waves consisted of 13 upregulated miRNAs that had regulatory effects on 48 downregulated genes. These differentially expressed miRNAs and genes may play a significant role in forming different patterns, and provide evidence for the molecular mechanisms underlying the formation of hair follicles of varying patterns.

Identification and Characterization of MicroRNAs in Snakehead Fish Cell Line upon Snakehead Fish Vesiculovirus Infection

MicroRNAs (miRNAs) play important roles in mediating multiple biological processes in eukaryotes and are being increasingly studied to evaluate their roles associated with cellular changes following viral infection. Snakehead fish Vesiculovirus (SHVV) has caused mass mortality in snakehead fish during the past few years. To identify specific miRNAs involved in SHVV infection, we performed microRNA deep sequencing on a snakehead fish cell line (SSN-1) with or without SHVV infection. A total of 205 known miRNAs were identified when they were aligned with the known zebrafish miRNAs, and nine novel miRNAs were identified using MiRDeep2 software. Eighteen and 143 of the 205 known miRNAs were differentially expressed at three and 24 h post-infection (poi), respectively. From the differentially-expressed miRNAs, five were randomly selected to validate their expression profiles using quantitative reverse transcription polymerase chain reaction (qRT-PCR), and their expression profiles were consistent with the microRNA sequencing results. In addition, the target gene prediction of the SHVV genome was performed for the differentially-expressed host miRNAs, and a total of 10 and 58 differentially-expressed miRNAs were predicted to bind to the SHVV genome at three and 24 h poi, respectively. The effects of three selected miRNAs (miR-130-5p, miR-214 and miR-216b) on SHVV multiplication were evaluated using their mimics and inhibitors via qRT-PCR and Western blotting. The results showed that all three miRNAs were able to inhibit the multiplication of SHVV; whereas the mechanisms underlying the SHVV multiplication inhibited by the specific miRNAs need to be further characterized in the future.

Expression of microRNA‑26b and identification of its target gene EphA2 in pituitary tissues in Yanbian cattle

microRNAs (miRNAs/miRs) are a class of single-stranded non-coding RNA molecules of 19–24 nucleotides (nt) in length. They are widely expressed in animals, plants, bacteria and viruses. Via specific mRNA complementary pairing of target genes, miRNAs are able to regulate the expression of mRNA levels or inhibit protein translation following transcription. miRNA expression has a time- and space specificity, and it is involved in cell proliferation and differentiation, apoptosis, development, tumor metastasis occurrence and other biological processes. miR-26b is an miRNA of 22 nt and is important in the regulation of cellular processes. With the advancement of molecular biology techniques in recent years, there have been extensive investigations into miR-26b. Numerous studies have observed that miR-26b is involved in early embryonic development, cell proliferation regulation, pituitary hormone secretion and other physiological activities. miRNAs are associated with the function of propagation. The present study used reverse transcription quantitative polymerase chain reaction to detect the relative expression levels of miR-26b in the pituitary tissue of Yanbian cattle at different developmental stages. The 2^−ΔΔCt method was used to calculate the relative gene expression levels. The miRNA target gene database TargetScan and RNA22 were used for prediction of the miR-26b target gene and selective recognition was also performed. The results demonstrated that miR-26b is expressed in the pituitary tissues of Yanbian cattle at 6 and 24 months of age. The relative expression levels of miR-26b in the pituitary tissues of 24-month-old Yanbian cattle were 2.41 times that of those in the six-month-old Yanbian cattle, demonstrating significant differences in the relative expression (P<0.01). The relative expression of the candidate target genes, EphA2 and miR-26b, exhibited the opposite expression pattern. The relative expression levels in the pituitary tissues of six-month-old Yanbian cattle were 3.34 times that of those in 24-month-old Yanbian cattle (P<0.01). There are miR-26b binding sites in the 3′-untranslated region (3′-UTR) of EphA2 in bovine, human, murine and other mammalian mRNAs, suggesting that the EphA2 gene may be a target gene of miR-26b. The results of a Luciferase reporter system assay revealed that miR-26b is able to suppress EphA2 expression at the transcription level. Following the site-directed mutagenesis of plasmid EphA2 3′-UTR pmirGLO-MUT- and miR-26b mimic-transfected HeLa cells, the dual-luciferase reporter gene assay revealed that there were three consecutive nucleotide mutations in the 3′-UTR, binding with the predicted seed region. This may have caused the miR-26b inhibition of luciferase activity to decrease from 60% in the wild-type to 26%, suggesting that miR-26b achieved its function via binding with the TACTTGAA sequence of the 3′-UTR in EphA2. In conclusion, the present study successfully assessed the expression pattern of miR-26b in the pituitary tissue of Yanbian cattle, and also confirmed that EphA2 was a target gene of miR-26b in Yanbian cattle in vitro. The present study provided the theoretical basis to further investigate the role of miR-26b in early embryonic development, pituitary hormone secretion and other reproductive functions.

Computational identification and characterization of miRNAs and their target genes from five cyprinidae fishes

MicroRNAs (miRNAs) are a kind of small single-strand RNA molecules with lengths of 18–25 nt, which do not encode any proteins. They play an essential role in gene expression regulation by binding to their target genes, leading to translational repression or transcript degradation. In this study, 23 miRNAs were predicted from five cyprinidae fishes by using a bioinformatics-based gene search based on blasting ESTs and GSS in NCBI, of which 21 miRNA genes have not been previously reported. To prove their validity, five of the computationally predicted miRNAs were verified by RTPCR, their transcripts were successfully detected, and, 46 potential target genes for these miRNAs were predicted, most target genes encode transcription factors, they are involved in signal transduction, metabolism and development processes.

Screening for characteristic microRNAs between pre-invasive and invasive stages of cervical cancer

The aim of the present study was to investigate the characteristic microRNAs (miRNAs) expressed during the pre-invasive and invasive stages of cervical cancer. A gene expression profile (GSE7803) containing 21 invasive squamous cell cervical carcinoma samples, 10 normal squamous cervical epithelium samples and seven high-grade squamous intraepithelial cervical lesion samples, was obtained from the Gene Expression Omnibus. Differentially expressed genes (DEGs) were identified using significance analysis of microarray software, and a Gene Ontology (GO) enrichment analysis was conducted using the Database for Annotation, Visualization and Integrated Discovery. The miRNAs that interacted with the identified DEGs were selected, based on the TarBase v5.0 database. Regulatory networks were constructed from these selected miRNAs along with their corresponding target genes among the DEGs. The regulatory networks were visualized using Cytoscape. A total of 1,160 and 756 DEGs were identified in the pre-invasive and invasive stages of cervical cancer, respectively. The results of the GO enrichment demonstrated that the DEGs were predominantly involved in the immune response and the cell cycle, in the pre-invasive and invasive stages, respectively. Furthermore, a total of 18 and 26 characteristic miRNAs were screened in the pre-invasive and invasive stages, respectively. These miRNAs may be potential biomarkers and targets for the diagnosis and treatment of the different stages of cervical cancer.

Oxysterol-binding protein-related protein 8 (ORP8) increases sensitivity of hepatocellular carcinoma cells to Fas-mediated apoptosis

Human hepatoma (HCC) has been reported to be strongly resistant to Fas-mediated apoptosis. However, the underlying mechanisms are poorly understood. In this study the function of oxysterol-binding protein-related protein 8 (ORP8) in human hepatoma cells apoptosis was assessed. We found that ORP8 is down-regulated, whereas miR-143, which controls ORP8 expression, is up-regulated in clinical HCC tissues as compared with liver tissue from healthy subjects. ORP8 overexpression triggered apoptosis in primary HCC cells and cell lines, which coincided with a relocation of cytoplasmic Fas to the cell plasma membrane and FasL up-regulation. Co-culture of HepG2 cells or primary HCC cells with Jurkat T-cells or T-cells, respectively, provided further evidence that ORP8 increases HCC cell sensitivity to Fas-mediated apoptosis. ORP8-induced Fas translocation is p53-dependent, and FasL was induced upon ORP8 overexpression via the endoplasmic reticulum stress response. Moreover, ORP8 overexpression and miR-143 inhibition markedly inhibited tumor growth in a HepG2 cell xenograft model. These results indicate that ORP8 induces HCC cell apoptosis through the Fas/FasL pathway. The role of ORP8 in Fas translocation to the plasma membrane and its down-regulation by miR-143 offer a putative mechanistic explanation for HCC resistance to apoptosis. ORP8 may be a potential target for HCC therapy.

Lessons from miR-143/145: the importance of cell-type localization of miRNAs

miR-143 and miR-145 are co-expressed microRNAs (miRNAs) that have been extensively studied as potential tumor suppressors. These miRNAs are highly expressed in the colon and are consistently reported as being downregulated in colorectal and other cancers. Through regulation of multiple targets, they elicit potent effects on cancer cell growth and tumorigenesis. Importantly, a recent discovery demonstrates that miR-143 and miR-145 are not expressed in colonic epithelial cells; rather, these two miRNAs are highly expressed in mesenchymal cells such as fibroblasts and smooth muscle cells. The expression patterns of miR-143 and miR-145 and other miRNAs were initially determined from tissue level data without consideration that multiple different cell types, each with their own unique miRNA expression patterns, make up each tissue. Herein, we discuss the early reports on the identification of dysregulated miR-143 and miR-145 expression in colorectal cancer and how lack of consideration of cellular composition of normal tissue led to the misconception that these miRNAs are downregulated in cancer. We evaluate mechanistic data from miR-143/145 studies in context of their cell type-restricted expression pattern and the potential of these miRNAs to be considered tumor suppressors. Further, we examine other examples of miRNAs being investigated in inappropriate cell types modulating pathways in a non-biological fashion. Our review highlights the importance of determining the cellular expression pattern of each miRNA, so that downstream studies are conducted in the appropriate cell type.

The miRNAome of the postpartum dairy cow liver in negative energy balance

Background

Negative energy balance (NEB) is an altered metabolic state in high yielding cows that occurs during the first few weeks postpartum when energy demands for lactation and maintenance exceed the energy supply from dietary intake. NEB can, in turn, lead to metabolic disorders and to reduced fertility. Alterations in the expression of more than 700 hepatic genes have previously been reported in a study of NEB in postpartum dairy cows. miRNAs (microRNA) are known to mediate many alterations in gene expression post transcriptionally. To study the hepatic miRNA content of postpartum dairy cows, including their overall abundance and differential expression, in mild NEB (MNEB) and severe NEB (SNEB), short read RNA sequencing was carried out. To identify putative targets of differentially expressed miRNAs among differentially expressed hepatic genes reported previously in dairy cows in SNEB computational target identification was employed.

Results

Our results indicate that the dairy cow liver expresses 53 miRNAs at a lower threshold of 10 reads per million. Of these, 10 miRNAs accounted for greater than 95% of the miRNAome (miRNA content). Of the highly expressed miRNAs, miR-122 constitutes 75% followed by miR-192 and miR-3596. Five out of thirteen let-7 miRNA family members are also among the highly expressed miRNAs. miR-143, down-regulated in SNEB, was found to have 4 putative up-regulated gene targets associated with SNEB including LRP2 (low density lipoprotein receptor-related protein 2), involved in lipid metabolism and up-regulated in SNEB.

Conclusions

This is the first liver miRNA-seq profiling study of moderate yielding dairy cows in the early postpartum period. Tissue specific miR-122 and liver enriched miR-192 are two of the most abundant miRNAs in the postpartum dairy cow liver. miR-143 is significantly down-regulated in SNEB and putative targets of miRNA-143 which are up-regulated in SNEB, include a gene involved in lipid metabolism.

Misregulated RNA processing in amyotrophic lateral sclerosis

Amyotrophic lateral sclerosis (ALS) research is undergoing an era of unprecedented discoveries with the identification of new genes as major genetic causes of this disease. These discoveries reinforce the genetic, clinical and pathological overlap between ALS and frontotemporal lobar degeneration (FTLD). Common causes of these diseases include mutations in the RNA/DNA-binding proteins, TDP-43 and FUS/TLS and most recently, hexanucleotide expansions in the C9orf72 gene, discoveries that highlight the overlapping pathogenic mechanisms that trigger ALS and FTLD. TDP-43 and FUS/TLS, both of which participate in several steps of RNA processing, are abnormally aggregated and mislocalized in ALS and FTLD, while the expansion in the C9orf72 pre-mRNA strongly suggests sequestration of one or more RNA binding proteins in pathologic RNA foci. Hence, ALS and FTLD converge in pathogenic pathways disrupting the regulation of RNA processing. This article is part of a Special Issue entitled RNA-Binding Proteins.

OLFML3 expression is decreased during prenatal muscle development and regulated by microRNA-155 in pigs

The Olfactomedin-like 3 (OLFML3) gene has matrix-related function involved in embryonic development. MicroRNA-155 (miR-155), 21- to 23-nucleotides (nt) noncoding RNA, regulated myogenesis by target mRNA. Our LongSAGE analysis suggested that OLFML3 gene was differently expressed during muscle development in pig. In this study, we cloned the porcine OLFML3 gene and detected its tissues distribution in adult Tongcheng pigs and dynamical expression in developmental skeletal muscle (12 prenatal and 10 postnatal stages) from Landrace (lean-type) and Tongcheng (obese-type) pigs. Subsequently, we analyzed the interaction between OLFML3 and miR-155. The OLFML3 was abundantly expressed in liver and pancreas, moderately in lung, small intestine and placenta, and weakly in other tissues and postnatal muscle. There were different dynamical expression patterns between Landrace and Tongcheng pigs during prenatal skeletal muscle development. The OLFML3 was down-regulated (33-50 days post coitus, dpc), subsequently up-regulated (50-70 dpc), and then down-regulated (70-100 dpc) in Landrace pigs, while in Tongcheng pigs, it was down-regulated (33-50 dpc), subsequently up-regulated (50-55 dpc) and then down-regulated (55-100 dpc). There was higher expression in Tongcheng than Landrace in prenatal muscle from 33 to 60 dpc, and opposite situation from 65 to 100 dpc. Dual luciferase assay and real time PCR documented that OLFML3 expression was regulated by miR-155 at mRNA level. Our research indicated that OLFML3 gene may affect prenatal skeletal muscle development and was regulated by miR-155. These finding will help understanding biological function and expression regulation of OLFML3 gene in mammal animals.