Genome-wide identification and functional analysis of long non-coding RNAs and mRNAs in male mice testes at the onset of puberty after low dose lead exposure

Highlights: • A cluster of RNA transcripts are changed in testes after low dose lead exposure. • Several lncRNAs might be implicated in the MAPK signaling pathway. • The differentially expressed genes were also altered in Sertoli cells. Many researchers have studied the relationship between lead (Pb) and testis injury, but the underlying mechanisms are still unknown. The participation of long non-coding RNAs (lncRNAs) in biological processes has been proposed. To comprehensively gain insight into the molecular toxicity of Pb, expression patterns are analysed through RNA sequencing (RNA-seq) in male mice treated with 200 mg/L of Pb through the drinking water for 90 days at the onset of puberty. A total of 614 differentially expressed (DE) lncRNAs were included (p ≤ 0.05 and fold change ≥2), of which 288 were up-regulated, and 326 were down-regulated. A total of 2295 DE mRNAs (p ≤ 0.05 and fold change ≥2), including 1202 up-regulated and 1093 down-regulated ones, were found in the testes of Pb-exposed group. Functional analysis results showed that several lncRNAs might be implicated in the bio-pathway of mitogen-activated protein kinase (MAPK) signaling pathway. Finally, seven pairs of lncRNA-mRNA co-expression were established in mice testes and confirmed by RT-qPCR. Moreover, the DE genes were also altered in Sertoli cells. Therefore, our research might be helpful for future exploring the effects of Pb exposure on lncRNA in testis, as well as its function.