Dynamic whole-transcriptome landscape of acute bilirubin encephalopathy in newborns

Abstract

Hyperbilirubinemia in newborns may progress to acute bilirubin encephalopathy (ABE), posing short- and long-term health risks. Despite extensive research identifying numerous mRNAs, lncRNAs, circRNAs, and miRNAs associated with brain injury, their roles in neonatal bilirubin-induced brain injury remain elusive. This study employed whole-transcriptome sequencing to ascertain the differentially expressed (DE) RNA profiles in a newborn ABE rat model, followed by bioinformatic analysis. A time-series competing endogenous RNA (ceRNA) regulatory network was established, and the expression trends of 9 arbitrarily chosen RNAs were verified through quantitative real-time polymerase chain reaction(qRT-PCR). In comparison with the control group, we identified 595, 888, and 1448 DE mRNAs; 22, 37, and 37 DE miRNAs; 1945, 1869, and 1997 DE lncRNAs; and 31, 28, and 36 DE circRNAs at 6 h, 12 h, and 24 h, respectively. Predominantly, these DERNAs contribute to biological functions and pathways associated with inflammation, immunity, metabolism, cell death, and neurodevelopmental regulation. Moreover, we constructed ceRNA networks of DE lncRNA/circRNA-DE miRNA-DE mRNA based on time series. The qRT-PCR expression trends for the selected 9 RNAs were generally similar to the RNA-seq outcomes. This investigation uniquely delineated the temporal expression patterns of mRNA and non-coding RNA in ABE, establishing ceRNA networks and identifying potential molecular mechanisms underlying bilirubin-induced hippocampal damage. Nonetheless, further studies are warranted to corroborate these findings in humans.