Integrated analysis of microRNA expression profiles and function network in mice testes after low dose lead exposure from early puberty

There is evidence suggesting the participation of non-coding RNAs in male reproductive dysfunction induced by lead, and the significance of microRNAs has been highlighted recently because of their essential roles in gene regulatory networks. To comprehensively understand the functions of miRNA and the regulatory networks, RNA sequencing was carried out to obtain miRNA expression profiles in mice testes exposed to low dose Pb for 90 days at the onset of puberty. In total, 44 differentially expressed miRNAs with 26 up-regulated and 18 down-regulated were identified between 200 mg/L Pb group and control group (p < 0.05). Enrichment analysis confirmed that the target genes of DE miRNAs might participate in the metabolism of testicular cells. Furthermore, a miRNA-mRNA co-expression network consisting of 19 miRNAs and 106 mRNAs and a competing endogenous RNA network of lncRNA-miRNA-mRNA including 179 genes were established. Finally, the expressions of 4 miRNAs (mmu-miR-451a, mmu-miR-133a-3p, mmu-miR-1a-3p and mmu-miR-486a-3p) and 4 mRNAs (Gramd1b, Tcf7l2, Mov10 and Srcin1) involved in regulatory networks were verified by RT-qPCR. In conclusion, our research might provide targets for the mechanism studies of miRNAs in reproductive toxicity of Pb.