MiRNA Expression Profile of Human Subcutaneous Adipose and during Adipocyte Differentiation

Background

Potential regulators of adipogenesis include microRNAs (miRNAs), small non-coding RNAs that have been recently shown related to adiposity and differentially expressed in fat depots. However, to date no study is available, to our knowledge, regarding miRNAs expression profile during human adipogenesis. Thereby, the aim of this study was to investigate whether miRNA pattern in human fat cells and subcutaneous adipose tissue is associated to obesity and co-morbidities and whether miRNA expression profile in adipocytes is linked to adipogenesis.

Methodology/Principal Findings

We performed a global miRNA expression microarray of 723 human and 76 viral mature miRNAs in human adipocytes during differentiation and in subcutaneous fat samples from non-obese (n = 6) and obese with (n = 9) and without (n = 13) Type-2 Diabetes Mellitus (DM-2) women. Changes in adipogenesis-related miRNAs were then validated by RT-PCR. Fifty of 799 miRNAs (6.2%) significantly differed between fat cells from lean and obese subjects. Seventy miRNAs (8.8%) were highly and significantly up or down-regulated in mature adipocytes as compared to pre-adipocytes. Otherwise, 17 of these 799 miRNAs (2.1%) were correlated with anthropometrical (BMI) and/or metabolic (fasting glucose and/or triglycerides) parameters. We identified 11 miRNAs (1.4%) significantly deregulated in subcutaneous fat from obese subjects with and without DM-2. Interestingly, most of these changes were associated with miRNAs also significantly deregulated during adipocyte differentiation.

Conclusions/Significance

The remarkable inverse miRNA profile revealed for human pre-adipocytes and mature adipocytes hints at a closely crosstalk between miRNAs and adipogenesis. Such candidates may represent biomarkers and therapeutic targets for obesity and obesity-related complications.     

Differential expression of microRNAs in decidua-derived mesenchymal stem cells from patients with pre-eclampsia

Background

Mesenchymal stem cells (MSCs) at maternal-fetal interface are considered to play an important role in the pathogenesis of pre-eclampsia (PE). microRNAs (miRNAs) also have an important influence on differentiation, maturation, and functions of MSCs. Our aim in this study was to determine the differential expression of miRNAs in decidua-derived MSCs (dMSCs) from severe PE and normal pregnancies.

Results

miRNA expression profiles in dMSCs from five patients with severe PE and five healthy pregnant women were screened using microarray. Then, bioinformatic analysis of the microarray results was performed. Out of 179 differentially expressed miRNAs, 49 miRNAs had significant (p < 0.05) differential expression of ≥ 2.0-fold changes, including 21 up-regulated and 28 down-regulated. miRNA-Gene-network and miRNA-Gene ontology (GO) -network analyses were performed. Overall, 21 up-regulated and 15 down-regulated miRNAs showed high degrees in these analyses. Moreover, the significantly enriched signaling pathways and GOs were identified. The analyses revealed that pathways associated with cell proliferation, angiogenesis, and immune functions were highly regulated by the differentially expressed miRNAs, including Wnt signaling pathway, mitogen-activated protein kinase signaling pathway, transforming growth factor beta signaling pathway, T-cell receptor signaling pathway, and B cell receptor signaling pathway. Four miRNA predicted target genes, vascular endothelial growth factor A (VEGFA), indoleamine 2,3-dioxygenase, suppression of cytokine signaling 3, and serine/threonine protein phosphatase 2A 55 kDa regulatory subunit B α isoform (PPP2R2A) were all decreased in dMSCs from patients with PE. Furthermore, the physiological roles of miR-16 and miR-136 in the down-regulation of VEGFA and PPP2R2A, respectively, were confirmed through reporter assays.

Conclusions

These findings suggest that miRNAs in dMSCs may be important regulatory molecules in the development of PE.     

Age and sex differences in kidney microRNA expression during the life span of F344 rats

Background

Growing evidence suggests that epigenetic mechanisms of gene regulation may play a role in susceptibilities to specific toxicities and adverse drug reactions. MiRNAs in particular have been shown to be important regulators in cancer and other diseases and show promise as predictive biomarkers for diagnosis and prognosis. In this study, we characterized the global kidney miRNA expression profile in untreated male and female F344 rats throughout the life span. These findings were correlated with sex-specific susceptibilities to adverse renal events, such as male-biased renal fibrosis and inflammation in old age.

Methods

Kidney miRNA expression was examined in F344 rats at 2, 5, 6, 8, 15, 21, 78, and 104 weeks of age in both sexes using Agilent miRNA microarrays. Differential expression was determined using filtering criteria of ≥1.5 fold change and ANOVA or pairwise t-test (FDR <5%) to determine significant age and sex effects, respectively. Pathway analysis software was used to investigate the possible roles of these target genes in age- and sex-specific differences.

Results

Three hundred eleven miRNAs were found to be expressed in at least one age and sex. Filtering criteria revealed 174 differentially expressed miRNAs in the kidney; 173 and 34 miRNAs exhibiting age and sex effects, respectively. Principal component analysis revealed age effects predominated over sex effects, with 2-week miRNA expression being much different from other ages. No significant sexually dimorphic miRNA expression was observed from 5 to 8 weeks, while the most differential expression (13 miRNAs) was observed at 21 weeks. Potential target genes of these differentially expressed miRNAs were identified.

Conclusions

The expression of 56% of detected renal miRNAs was found to vary significantly with age and/or sex during the life span of F344 rats. Pathway analysis suggested that 2-week-expressed miRNAs may be related to organ and cellular development and proliferation pathways. Male-biased miRNA expression at older ages correlated with male-biased renal fibrosis and mononuclear cell infiltration. These miRNAs showed high representation in renal inflammation and nephritis pathways, and included miR-214, miR-130b, miR-150, miR-223, miR-142-5p, miR-185, and miR-296*. Analysis of kidney miRNA expression throughout the rat life span will improve the use of current and future renal biomarkers and inform our assessments of kidney injury and disease.

Electronic supplementary material

The online version of this article (doi:10.1186/s13293-014-0019-1) contains supplementary material, which is available to authorized users.     

MicroRNA Profiling of BRCA1/2 Mutation-Carrying and Non-Mutation-Carrying High-Grade Serous Carcinomas of Ovary

Background

MicroRNAs (miRNA) are 20∼25 nucleotide non-coding RNAs that inhibit the translation of targeted mRNA, and they have been implicated in the development of human malignancies. High grade serous ovarian carcinomas, the most common and lethal subtype of ovarian cancer, can occur sporadically or in the setting of BRCA1/2 syndromes. Little is known regarding the miRNA expression profiles of high grade serous carcinoma in relation to BRCA1/2 status, and compared to normal tubal epithelium, the putative tissue of origin for high grade serous carcinomas.

Methodology/Principal Findings

Global miRNA expression profiling was performed on a series of 33 high grade serous carcinomas, characterized with respect to BRCA1/2 status (mutation, epigenetic silencing with loss of expression or normal), and with clinical follow-up, together with 2 low grade serous carcinomas, 2 serous borderline tumors, and 3 normal fallopian tube samples, using miRNA microarrays (328 human miRNA). Unsupervised hierarchical clustering based on miRNA expression profiles showed no clear separation between the groups of carcinomas with different BRCA1/2 status. There were relatively few miRNAs that were differentially expressed between the genotypic subgroups. Comparison of 33 high grade serous carcinomas to 3 normal fallopian tube samples identified several dysregulated miRNAs (false discovery rate <5%), including miR-422b and miR-34c. Quantitative RT-PCR analysis performed on selected miRNAs confirmed the pattern of differential expression shown by microarray analysis. Prognostically, lower level miR-422b and miR-34c in high grade serous carcinomas were both associated with decreased disease-specific survival by Kaplan-Meier analysis (p<0.05).

Conclusions/Significance

High grade serous ovarian carcinomas with and without BRCA1/2 abnormalities demonstrate very similar miRNA expression profiles. High grade serous carcinomas as a group exhibit significant miRNA dysregulation in comparison to tubal epithelium and the levels of miR-34c and miR-422b appear to be prognostically important.     

Profiling the circulating miRNAs in mice exposed to gram-positive and gram-negative bacteria by Illumina small RNA deep sequencing

Background

We profiled the expression of circulating microRNAs (miRNAs) in mice using Illumina small RNA deep sequencing in order to identify the miRNAs that may potentially be used as biomarkers to distinguish between gram-negative and gram-positive bacterial infections.

Results

Recombinant-specific gram-negative pathogen Escherichia coli (Xen14) and gram-positive pathogen Staphylococcus aureus (Xen29) were used to induce bacterial infection in mice at a concentration of 1 × 10⁸ bacteria/100 μL of phosphate buffered saline (PBS). Small RNA libraries generated from the serum of mice after exposure to PBS, Xen14, Xen29, and Xen14 + Xen29 via the routes of subcutaneous injection (I), cut wound (C), or under grafted skin (S) were analyzed using an Illumina HiSeq2000 Sequencer. Following exposure to gram-negative bacteria alone, no differentially expressed miRNA was found in the injection, cut, or skin graft models. Exposure to mixed bacteria induced a similar expression pattern of the circulating miRNAs to that induced by gram-positive bacterial infection. Upon gram-positive bacterial infection, 9 miRNAs (mir-193b-3p, mir-133a-1-3p, mir-133a-2-3p, mir-133a-1-5p, mir-133b-3p, mir-434-3p, mir-127-3p, mir-676-3p, mir-215-5p) showed upregulation greater than 4-fold with a p-value < 0.01. Among them, mir-193b-3p, mir-133a-1-3p, and mir-133a-2-3p presented the most common miRNA targets expressed in the mice exposed to gram-positive bacterial infection.

Conclusions

This study identified mir-193b-3p, mir-133a-1-3p, and mir-133a-2-3p as potential circulating miRNAs for gram-positive bacterial infections.

Electronic supplementary material

The online version of this article (doi:10.1186/s12929-014-0106-y) contains supplementary material, which is available to authorized users.     

Profiling circulating microRNA expression in a mouse model of nerve allotransplantation

Background

The lack of noninvasive biomarkers of rejection remains a challenge in the accurate monitoring of deeply buried nerve allografts and precludes optimization of therapeutic intervention. This study aimed to establish the expression profile of circulating microRNAs (miRNAs) during nerve allotransplantation with or without immunosuppression.

Results

Balb/c mice were randomized into 3 experimental groups, that is, (1) untreated isograft (Balb/c → Balb/c), (2) untreated allograft (C57BL/6 → Balb/c), and (3) allograft (C57BL/6 → Balb/c) with FK506 immunosuppression. A 1-cm Balb/c or C57BL/6 donor sciatic nerve graft was transplanted into sciatic nerve gaps created in recipient mice. At 1, 3, 7, 10, and 14 d after nerve transplantation, nerve grafts, whole blood, and sera were obtained for miRNA expression analysis with an miRNA array and subsequent validation with quantitative real-time PCR (qRT-PCR). Three circulating miRNAs (miR-320, miR-762, and miR-423-5p) were identified in the whole blood and serum of the mice receiving an allograft with FK506 immunosuppression, within 2 weeks after nerve allotransplantation. However, these 3 circulating miRNAs were not expressed in the nerve grafts. The expression of all these 3 upregulated circulating miRNAs significantly decreased at 2, 4, and 6 d after discontinuation of FK506 immunosuppression. In the nerve graft, miR-125-3b and miR-672 were significantly upregulated in the mice that received an allograft with FK506 only at 7 d after nerve allotransplantation.

Conclusions

We identified the circulating miR-320, miR-762, and miR-423-5p as potential biomarkers for monitoring the immunosuppression status of the nerve allograft. However, further research is required to investigate the mechanism behind the dysregulation of these markers and to evaluate their prognostic value in nerve allotransplantation.     

Rare Circulating MicroRNAs as Biomarkers of Colorectal Neoplasia

Background

MicroRNAs (miRNAs) are regulatory RNAs, stable in circulation, and implicated in colorectal cancer (CRC) etiology and progression. Therefore they are promising as early detection biomarkers of colorectal neoplasia. However, many circulating miRNAs are highly expressed in blood cells, and therefore may not be specific to colorectal neoplasia.

Methods

We selected 7 miRNA candidates with previously reported elevated expression in adenoma tissue but low expression in blood cells (“rare” miRNAs), 2 previously proposed as adenoma biomarkers, and 3 implicated in CRC. We conducted a colonoscopy-based case-control study including 48 polyp-free controls, 43 advanced adenomas, 73 non-advanced adenomas, and 8 CRC cases. miRNAs from plasma were quantified by qRT-PCR. Correlations between miRNA expression levels, adjusted for age and sex, were assessed. We used polytomous logistic regression to estimate odds ratios (ORs) and 95% confidence intervals quantifying the association between expression levels of miRNAs and case groups. We also conducted nonparametric receiver operating characteristic (ROC) analyses and estimated area under the curve (AUC).

Results

miRNAs with high expression levels were statistically significantly correlated with one another. No miRNAs were significantly associated with non-advanced or advanced adenomas. Strong (ORs >5) and significant associations with CRC were observed for 6 miRNA candidates, with corresponding AUCs significantly >0.5.

Conclusions

These candidate miRNAs, assayed by qRT-PCR, are probably unsuitable as blood-based adenoma biomarkers. Strong associations between miRNAs and CRC were observed, but primarily with miRNAs highly expressed in blood cells. These results suggest that rare miRNAs will require new detection methods to serve as circulating biomarkers of adenomas.     

Identification of Restricted Subsets of Mature microRNA Abnormally Expressed in Inactive Colonic Mucosa of Patients with Inflammatory Bowel Disease

Background

Ulcerative Colitis (UC) and Crohn''s Disease (CD) are two chronic Inflammatory Bowel Diseases (IBD) affecting the intestinal mucosa. Current understanding of IBD pathogenesis points out the interplay of genetic events and environmental cues in the dysregulated immune response. We hypothesized that dysregulated microRNA (miRNA) expression may contribute to IBD pathogenesis. miRNAs are small, non-coding RNAs which prevent protein synthesis through translational suppression or mRNAs degradation, and regulate several physiological processes.

Methodology/Findings

Expression of mature miRNAs was studied by Q-PCR in inactive colonic mucosa of patients with UC (8), CD (8) and expressed relative to that observed in healthy controls (10). Only miRNAs with highly altered expression (>5 or <0.2 -fold relative to control) were considered when Q-PCR data were analyzed. Two subsets of 14 (UC) and 23 (CD) miRNAs with highly altered expression (5.2->100 -fold and 0.05–0.19 -fold for over- and under- expression, respectively; 0.001<p≤0.05) were identified in quiescent colonic mucosa, 8 being commonly dysregulated in non-inflamed UC and CD (mir-26a,-29a,-29b,-30c,-126*,-127-3p,-196a,-324-3p). Several miRNA genes with dysregulated expression co-localize with acknowledged IBD-susceptibility loci while others, (eg. clustered on 14q32.31), map on chromosomal regions not previously recognized as IBD-susceptibility loci. In addition, in silico clustering analysis identified 5 miRNAs (mir-26a,-29b,-126*,-127-3p,-324-3p) that share coordinated dysregulation of expression both in quiescent and in inflamed colonic mucosa of IBD patients. Six miRNAs displayed significantly distinct alteration of expression in non-inflamed colonic biopsies of UC and CD patients (mir-196b,-199a-3p,-199b-5p,-320a,-150,-223).

Conclusions/Significance

Our study supports miRNAs as crucial players in the onset and/or relapse of inflammation from quiescent mucosal tissues in IBD patients. It allows speculating a role for miRNAs as contributors to IBD susceptibility and suggests that some of the miRNA with altered expression in the quiescent mucosa of IBD patients may define miRNA signatures for UC and CD and help develop new diagnostic biomarkers.     

Ultra-deep sequencing reveals the microRNA expression pattern of the human stomach

Background

While microRNAs (miRNAs) play important roles in tissue differentiation and in maintaining basal physiology, little is known about the miRNA expression levels in stomach tissue. Alterations in the miRNA profile can lead to cell deregulation, which can induce neoplasia.

Methodology/Principal Findings

A small RNA library of stomach tissue was sequenced using high-throughput SOLiD sequencing technology. We obtained 261,274 quality reads with perfect matches to the human miRnome, and 42% of known miRNAs were identified. Digital Gene Expression profiling (DGE) was performed based on read abundance and showed that fifteen miRNAs were highly expressed in gastric tissue. Subsequently, the expression of these miRNAs was validated in 10 healthy individuals by RT-PCR showed a significant correlation of 83.97% (P<0.05). Six miRNAs showed a low variable pattern of expression (miR-29b, miR-29c, miR-19b, miR-31, miR-148a, miR-451) and could be considered part of the expression pattern of the healthy gastric tissue.

Conclusions/Significance

This study aimed to validate normal miRNA profiles of human gastric tissue to establish a reference profile for healthy individuals. Determining the regulatory processes acting in the stomach will be important in the fight against gastric cancer, which is the second-leading cause of cancer mortality worldwide.     

The role of bronchial epithelial cells in the pathogenesis of COPD in Z-alpha-1 antitrypsin deficiency

Background

Alpha-1 antitrypsin is the main inhibitor of neutrophil elastase in the lung. Although it is principally synthesized by hepatocytes, alpha-1 antitrypsin is also secreted by bronchial epithelial cells. Gene mutations can lead to alpha-1 antitrypsin deficiency, with the Z variant being the most clinically relevant due to its propensity to polymerize. The ability of bronchial epithelial cells to produce Z-variant protein and its polymers is unknown.We investigated the expression, accumulation, and secretion of Z-alpha-1 antitrypsin and its polymers in cultures of transfected cells and in cells originating from alpha-1 antitrypsin-deficient patients.

Methods

Experiments using a conformation-specific antibody were carried out on M- and Z-variant–transfected 16HBE cells and on bronchial biopsies and ex vivo bronchial epithelial cells from Z and M homozygous patients. In addition, the effect of an inflammatory stimulus on Z-variant polymer formation, elicited by Oncostatin M, was investigated. Comparisons of groups were performed using t-test or ANOVA. Non-normally distributed data were assessed by Mann–Whitney U test or the Kruskal-Wallis test, where appropriate. A P value of < 0.05 was considered to be significant.

Results

Alpha-1 antitrypsin polymers were found at a higher concentration in the culture medium of ex vivo bronchial epithelial cells from Z-variant homozygotes, compared with M-variant homozygotes (P < 0.01), and detected in the bronchial epithelial cells and submucosa of patient biopsies. Oncostatin M significantly increased the expression of alpha-1 antitrypsin mRNA and protein (P < 0.05), and the presence of Z-variant polymers in ex vivo cells (P < 0.01).

Conclusions

Polymers of Z-alpha-1 antitrypsin form in bronchial epithelial cells, suggesting that these cells may be involved in the pathogenesis of lung emphysema and in bronchial epithelial cell dysfunction.

Electronic supplementary material

The online version of this article (doi:10.1186/s12931-014-0112-3) contains supplementary material, which is available to authorized users.     

Circulating microRNAs have a sex-specific association with metabolic syndrome

Background

The microRNAs let-7 g and miR-221 have been demonstrated to be related to the glucose metabolism. This study assessed the serum levels of these two microRNAs in subjects with and without metabolic syndrome (MetS).

Results

The serum microRNA levels were detected in 102 subjects aged 40 to 80 years who were recruited from the general population. The status of MetS was defined by the Adult Treatment Panel III (ATP III) criteria modified for Asians. Subjects with histories of cardiovascular diseases or who were receiving treatment with hypoglycemic or lipid-lowering agents were excluded. The levels of both circulating microRNAs (let-7 g and miR-221) were higher in subjects with MetS (p = 0.004 and p = 0.01, respectively). The sex-specific analysis showed that the difference was more prominent in women (for both miRNAs, p < 0.05 in women and p > 0.1 in men). In the female subjects, increased expression of both microRNAs was associated with an increased number of MetS risk components (p = 0.002 for let-7 g and p = 0.022 for miR-221). Moreover, the elevation of serum let-7 g was significantly associated with a low level of high-density lipoprotein cholesterol (p = 0.022) and high blood pressure (p = 0.023). In contrast, the miR-221 level was not associated with any individual MetS risk component.

Conclusions

The circulating levels of let-7 g and miR-221 displayed a female-specific elevation in individuals with metabolic syndrome.     

Inhibition of LPS-induced airway neutrophilic inflammation in healthy volunteers with an oral CXCR2 antagonist

Background

Inhaled lipopolysaccharide (LPS) induces a dose-dependent, acute neutrophilic response in the airways of healthy volunteers that can be quantified in induced sputum. Chemokines, such as CXCL1 and CXCL8, play an important role in neutrophilic inflammation in the lung through the activation of CXCR2 and small molecule antagonists of these receptors have now been developed. We investigated the effect of AZD8309, a CXCR2 antagonist, compared with placebo on LPS-induced inflammation measured in sputum of healthy volunteers.

Methods

Twenty healthy subjects were randomized in a double-blind placebo-controlled, cross-over study. AZD8309 (300 mg) or placebo was dosed twice daily orally for 3 days prior to challenge with inhaled LPS and induced sputum was collected 6 h later.

Results

Treatment with AZD8309 showed a mean 77% reduction in total sputum cells (p < 0.001) and 79% reduction in sputum neutrophils (p < 0.05) compared with placebo after LPS challenge. There was also a reduction in neutrophil elastase activity (p < 0.05) and CXCL1 (p < 0.05) and trends for reductions in sputum macrophages (47%), leukotriene B₄ (39%) and CXCL8 (52%).

Conclusions

AZD8309 inhibited LPS-induced inflammation measured in induced sputum of normal volunteers, indicating that this treatment may be useful in the treatment of neutrophilic diseases of the airways, such as COPD, severe asthma and cystic fibrosis.

Trial registration

{"type":"clinical-trial","attrs":{"text":"NCT00860821","term_id":"NCT00860821"}}NCT00860821.     

Roxithromycin treatment inhibits TGF-β1-induced activation of ERK and AKT and down-regulation of Caveolin-1 in rat airway smooth muscle cells

Background

Roxithromycin (RXM) has been widely used in asthma treatment; however, the mechanism has not been fully understood. The aim of our study was to investigate the underlying mechanism of RXM treatment in mediating the effect of transforming growth factor (TGF)-β1 on airway smooth muscle cells (ASMCs) proliferation and caveolinn-1 expression.

Methods

Firstly, the rat ovalbumin (OVA) model was built according to the previous papers. Rat ASMCs were prepared and cultured, and then TGF-β1 production in ASMCs was measured by enzyme-linked immunosorbent assay (ELISA). Moreover, the proliferation of ASMCs was determined using cell counting kit (CCK-8) assay. Additionally, the expressions of caveolin-1, phosphorylated-ERK1/2 (p-ERK1/2) and phosphorylated–AKT (p-AKT) in ASMCs treated with or without PD98059 (an ERK1/2 inhibitor), wortannin (a PI3K inhibitor), β-cyclodextrin (β-CD) and RXM were measured by Western blot. Finally, data were evaluated using t–test or one-way ANOVA, and then a P value < 0.05 was set as a threshold.

Results

Compared with normal control, TGF-β1 secretion was significantly increased in asthmatic ASMCs; meanwhile, TGF-β1 promoted ASMCs proliferation (P < 0.05). However, ASMCs proliferation was remarkably inhibited by RXM, β-CD, PD98059 and wortmannin (P < 0.05). Moreover, the expressions of p-ERK1/2 and p-AKT were increased and peaked at 20 min after TGF-β1 stimulation, and then suppressed by RXM. Further, caveolin-1 level was down-regulated by TGF-β1 and up-regulated by inhibitors and RXM.

Conclusion

Our findings demonstrate that RXM treatment inhibits TGF-β1-induced activation of ERK and AKT and down-regulation of caveolin-1, which may be the potential mechanism of RXM protection from chronic inflammatory diseases, including bronchial asthma.     

Temperature during early development has long-term effects on microRNA expression in Atlantic cod

Background

Environmental temperature has serious implications in life cycle of aquatic ectotherms. Understanding the molecular mechanisms of temperature acclimation and adaptation of marine organisms is of the uttermost importance for ecology, fisheries, and aquaculture, as it allows modeling the effects of global warming on population dynamics. Regulatory molecules are major modulators of acclimation and adaptation; among them, microRNAs (miRNAs) are versatile and substantial contributors to regulatory networks of development and adaptive plasticity. However, their role in thermal plasticity is poorly known. We have asked whether the temperature and its shift during the early ontogeny (embryonic and larval development) affect the miRNA repertoire of Atlantic cod (Gadus morhua), and if thermal experience has long-term consequences in the miRNA profile.

Results

We characterized miRNA during different developmental stages and in juvenile tissues using next generation sequencing. We identified 389 putative miRNA precursor loci, 120 novel precursor miRNAs, and 281 mature miRNAs. Some miRNAs showed stage- or tissue-enriched expression and miRNAs, such as the miR-17 ~ 92 cluster, myomiRs (miR-206), neuromiRs (miR-9, miR-124), miR-130b, and miR-430 showed differential expression in different temperature regimes. Long-term effect of embryonic incubation temperature was revealed on expression of some miRNAs in juvenile pituitary (miR-449), gonad (miR-27c, miR-30c, and miR-200a), and liver (let-7 h, miR-7a, miR-22, miR-34c, miR-132a, miR-192, miR-221, miR-451, miR-2188, and miR-7550), but not in brain. Some of differentially expressed miRNAs in the liver were confirmed using LNA-based rt-qPCR. The effect of temperature on methylation status of selected miRNA promoter regions was mostly inconclusive.

Conclusions

Temperature elevation by several degrees during embryonic and larval developmental stages significantly alters the miRNA profile, both short-term and long-term. Our results suggest that a further rise in seas temperature might affect life history of Atlantic cod.

Electronic supplementary material

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