Assessment of myometrial transcriptome changes associated with spontaneous human labour by high-throughput RNA-seq

The transition of the human uterus from a quiescent to a contractile state takes place over a number of weeks. On such biological time scales, cellular phenotype is modified by changes in the transcriptome, which in turn is under the control of the underlying endocrine, paracrine, and biophysical processes resulting from the ongoing pregnancy. In this study, we characterize the transition of the human myometrial transcriptome at term from not in labour (NIL) to in labour (LAB) using high throughput RNA sequencing (RNA-seq). RNA was isolated from the myometrium of uterine biopsies from patients at term who were not in labour (n = 5) and at term in spontaneous labour (n = 5) without augmentation. A total of 143.6 million separate reads were sequenced, achieving, on average, ∼13 times coverage of the expressed human transcriptome per sample. Principal component analysis indicated that the NIL and LAB transcriptomes could be distinguished as two distinct clusters. A comparison of the NIL and LAB groups, using three different statistical approaches (baySeq, edgeR, and DESeq), demonstrated an overlap of 764 differentially expressed genes. A comparison with currently available microarray data revealed only a partial overlap in differentially expressed genes. We conclude that the described RNA-seq data sets represent the first fully annotated catalogue of expressed mRNAs in human myometrium. When considered together, the full expression repertoire and the differentially expressed gene sets should provide an excellent resource for formulating new hypotheses of physiological function, as well as the discovery of novel therapeutic targets.