Endometrial receptivity revisited: endometrial transcriptome adjusted for tissue cellular heterogeneity
Authors: Suhorutshenko Marina , Kukushkina Viktorija , Velthut-Meikas Agne , Altmäe Signe , Peters Maire , Mägi Reedik , Krjutškov Kaarel , Koel Mariann , Martinez-Blanch Juan Fco. , Codoner Francisco M. , Vilella Felipe , Simon Carlos , Salumets Andres , Laisk Triin 
Published: July 2018

AbstractSTUDY QUESTIONDoes cellular composition of the endometrial biopsy affect the gene expression profile of endometrial whole-tissue samples?SUMMARY ANSWERThe differences in epithelial and stromal cell proportions in endome-trial biopsies modify whole-tissue gene expression profiles, and also affect the results of differential expression analysis.WHAT IS ALREADY KNOWNEach cell type has its unique gene expression profile. The proportions of epithelial and stromal cells vary in endometrial tissue during the menstrual cycle, along with individual and technical variation due to the way and tools used to obtain the tissue biopsy.STUDY DESIGN, SIZE, DURATIONUsing cell-population specific transcriptome data and computational deconvolution approach, we estimated the epithelial and stromal cell proportions in whole-tissue biopsies taken during early secretory and mid-secretory phases. The estimated cellular proportions were used as covariates in whole-tissue differential gene expression analysis. Endometrial transcriptomes before and after deconvolution were compared and analysed in biological context.PARTICIPANTS/MATERIAL, SETTING, METHODSPaired early- and mid-secretory endometrial biopsies were obtained from thirty-five healthy, regularly cycling, fertile volunteers, aged 23 to 36 years, and analysed by RNA sequencing. Differential gene expression analysis was performed using two approaches. In one of them, computational deconvolution was applied as an intermediate step to adjust for epithelial and stromal cells’ proportions in endometrial biopsy. The results were then compared to conventional differential expression analysis.MAIN RESULTS AND THE ROLE OF CHANCEThe estimated average proportions of stromal and epithelial cells in early secretory phase were 65% and 35%, and during mid-secre-tory phase 46% and 54%, respectively, that correlated well with the results of histological evaluation (r=0.88, p=1.1×10−6). Endometrial tissue transcriptomic analysis showed that approximately 26% of transcripts (n=946) differentially expressed in receptive endometrium in cell-type unadjusted analysis also remain differentially expressed after adjustment for biopsy cellular composition. However, the other 74% (n=2,645) become statistically non-significant after adjustment for biopsy cellular composition, underlining the impact of tissue heterogeneity on differential expression analysis. The results suggest new mechanisms involved in endometrial maturation involving genes like LINC01320, SLC8A1 and GGTA1P, described for the first time in context of endometrial receptivity.LIMITATIONS, REASONS FOR CAUTIONOnly dominant endometrial cell types were considered in gene expression profile deconvolution; however, other less frequent endometrial cell types also contribute to the whole-tissue gene expression profile.WIDER IMPLICATIONS OF THE FINDINGSThe better understanding of molecular processes during transition from pre-receptive to receptive endometrium serves to improve the effectiveness and personalization of assisted reproduction protocols. Biopsy cellular composition should be taken into account in future endometrial ‘omics’ studies, where tissue heterogeneity could potentially influence the results.TRIAL REGISTRATION NON/A