RNA-Sequencing of Umbilical Cord Blood to Investigate Spontaneous Preterm Birth: A Pilot Study

Validation of reference genes for whole blood gene expression analysis in cord blood of preterm and full-term neonates and peripheral blood of healthy adults

Background

Preterm birth is the leading cause of neonatal morbidity and mortality, but research efforts in neonatology are complicated due to the unavailability of large volume blood samples. Whole blood assays can be used to overcome this problem by performing both functional and gene expression studies using small amounts of blood. Gene expression studies using RT-qPCR estimate mRNA-levels of target genes normalized to reference genes. The goal of this study was to identify and validate stable reference genes applicable to cord blood samples obtained from developing neonates of different gestational age groups as well as to adult peripheral blood samples. Eight reference gene candidates (ACTB, B2M, GAPDH, GUSB, HPRT, PPIB, RPLP0, RPL13) were analyzed using the three published software algorithms Bestkeeper, GeNorm and NormFinder.

Results

A normalization factor consisting of ACTB and PPIB allows for comparative expression analyses of neonatal samples from different gestational age groups. Normalization factors consisting of GAPDH and PPIB or ACTB and GAPDH are suitable when samples from preterm and full-term neonates and adults are compared. However, all candidate reference genes except RPLP0 exhibited significant intergroup gene expression variance and a higher gene expression towards an older age which resulted in a small but statistically significant systematic bias. Systematic analysis of RNA-seq data revealed new reference gene candidates with potentially superior stability.

Conclusions

The current study identified suitable normalization factors and proposed the use of the additional single gene RPLP0 to avoid systematic bias. This combination will enable comparative analyses not only between neonates of different gestational ages, but also between neonates and adults, as it facilitates more detailed investigations of developmental gene expression changes. The use of software algorithms did not prevent unintended systematic bias. This generally highlights the need for careful validation of such results to prevent false interpretation of potential age-dependent changes in gene expression. To identify the most stable reference genes in the future, RNA-seq based global approaches are recommended.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07801-0.

Comparative Analysis of Global Gene Expression and Complement Components Levels in Umbilical Cord Blood from Preterm and Term Neonates: Implications for Significant Downregulation of Immune Response Pathways related to Prematurity

Background: Preterm birth is the most frequent cause of neonatal death, but its aetiology remains unclear. It has been suggested that the imbalance of immunological mechanisms responsible for maintaining pregnancy is contributing to preterm birth pathogenesis. We aimed to investigate global gene expression and the levels of several complement system components in umbilical cord blood samples from preterm neonates and compare them to term newborns. We sought to examine how differentially expressed genes could affect various immune-related pathways that are believed to be crucial factors in preterm birth. Material and methods: We enrolled 27 preterm infants (<37 weeks GA) and 52 term infants (>37 weeks GA), from which umbilical cord blood samples were collected. From these samples, peripheral blood mononuclear cells were isolated and subsequent RNA isolation was performed. We used Affymetrix Human Gene 2.1 ST Array Strip for microarray experiment and DAVID resources for bioinformatics analysis of the obtained data. Concentrations of C2, C3a, C5/C5a, C9, FactorD, Properdin were measured in umbilical cord blood plasma samples using multiplex fluorescent bead-based immunoassays using Luminex technology. Results: The levels of C3a and C5/5a were significantly elevated in preterm neonates compared to term babies, whereas C9 concentration was evidently increased in term infants. The expression of 250 genes was upregulated at least 2-fold and 3781 genes were downregulated at least 2-fold in preterm neonates in comparison with term infants. Functional annotation analysis revealed that in preterm infants in comparison to term babies there was a significant downregulation of genes encoding several Toll-like receptors, interleukins and genes involved in major signalling pathways (e.g. NF-κB, MAPK, TNF, Notch, JAK) and vital cellular processes (e.g. intracellular signal transduction, protein ubiquitination, protein transport, RNA splicing, DNA-templated transcription). Conclusions: Preterm birth results in immediate and long-term complications. Our results indicate that infants born prematurely show significant differences in complement components concentration and a downregulation of over 3,000 genes, involved mainly in various immune-related pathways, including innate immune response, phagocytosis and TLR function, when compared to full-term babies. Further studies on larger cohorts are needed to elucidate the role of immunity in prematurity.