Targeted gene expression profiling for accurate endometrial receptivity testing.

Expressional profiling of the endometrium enables the personalised timing of the window of implantation (WOI). This study presents and evaluates a novel analytical pipeline based on a TAC-seq (Targeted Allele Counting by sequencing) method for endometrial dating. The expressional profiles were clustered, and differential expression analysis was performed on the model development group, using 63 endometrial biopsies spanning over proliferative (PE, n = 18), early-secretory (ESE, n = 18), mid-secretory (MSE, n = 17) and late-secretory (LSE, n = 10) endometrial phases of the natural cycle. A quantitative predictor model was trained on the development group and validated on sequenced samples from healthy women, consisting of 52 paired samples taken from ESE and MSE phases and five LSE phase samples from 31 individuals. Finally, the developed test was applied to 44 MSE phase samples from a study group of patients diagnosed with recurrent implantation failure (RIF). In validation samples (n = 57), we detected displaced WOI in 1.8% of the samples from fertile women. In the RIF study group, we detected a significantly higher proportion of the samples with shifted WOI than in the validation set of samples from fertile women, 15.9% and 1.8% (p = 0.012), respectively. The developed model was evaluated with an average cross-validation accuracy of 98.8% and an accuracy of 98.2% in the validation group. The developed beREADY screening model enables sensitive and dynamic detection of selected transcriptome biomarkers, providing a quantitative and accurate prediction of endometrial receptivity status.

GWAS meta-analyses clarify the genetics of cervical phenotypes and inform risk stratification for cervical cancer.

Genome-wide association studies (GWAS) have successfully identified associations for cervical cancer, but the underlying mechanisms of cervical biology and pathology remain uncharacterised. Our GWAS meta-analyses fill this gap, as we characterise the genetic architecture of cervical phenotypes, including cervical ectropion, cervicitis, cervical dysplasia, as well as up to 9229 cases and 490 304 controls for cervical cancer from diverse ancestries. Leveraging the latest computational methods and gene expression data, we refine the association signals for cervical cancer and propose potential causal variants and genes at each locus. We prioritise PAX8/PAX8-AS1, LINC00339, CDC42, CLPTM1L, HLA-DRB1 and GSDMB as the most likely candidate genes for cervical cancer signals, providing insights into cervical cancer pathogenesis and supporting the involvement of reproductive tract development, immune response and cellular proliferation/apoptosis. We construct a genetic risk score (GRS) that is associated with cervical cancer [hazard ratios (HR) = 3.1 (1.7-5.6) for the top 15% vs lowest 15% of individuals], and with other HPV- and immune-system-related diagnoses in a phenome-wide association study analysis. Our results propose valuable leads for further functional studies and present a GRS for cervical cancer that allows additional risk stratification and could potentially be used to personalise the conventional screening strategies for groups more susceptible to cervical cancer.

Human endometrial cell-type-specific RNA sequencing provides new insights into the embryo-endometrium interplay.

STUDY QUESTION: Which genes regulate receptivity in the epithelial and stromal cellular compartments of the human endometrium, and which molecules are interacting in the implantation process between the blastocyst and the endometrial cells? SUMMARY ANSWER: A set of receptivity-specific genes in the endometrial epithelial and stromal cells was identified, and the role of galectins (LGALS1 and LGALS3), integrin ß1 (ITGB1), basigin (BSG) and osteopontin (SPP1) in embryo-endometrium dialogue among many other protein-protein interactions were highlighted. WHAT IS KNOWN ALREADY: The molecular dialogue taking place between the human embryo and the endometrium is poorly understood due to ethical and technical reasons, leaving human embryo implantation mostly uncharted. STUDY DESIGN SIZE DURATION: Paired pre-receptive and receptive phase endometrial tissue samples from 16 healthy women were used for RNA sequencing. Trophectoderm RNA sequences were from blastocysts. PARTICIPANTS/MATERIALS SETTING METHODS: Cell-type-specific RNA-seq analysis of freshly isolated endometrial epithelial and stromal cells using fluorescence-activated cell sorting (FACS) from 16 paired pre-receptive and receptive tissue samples was performed. Endometrial transcriptome data were further combined in silico with trophectodermal gene expression data from 466 single cells originating from 17 blastocysts to characterize the first steps of embryo implantation. We constructed a protein-protein interaction network between endometrial epithelial and embryonal trophectodermal cells, and between endometrial stromal and trophectodermal cells, thereby focusing on the very first phases of embryo implantation, and highlighting the molecules likely to be involved in the embryo apposition, attachment and invasion. MAIN RESULTS AND THE ROLE OF CHANCE: In total, 499 epithelial and 581 stromal genes were up-regulated in the receptive phase endometria when compared to pre-receptive samples. The constructed protein-protein interactions identified a complex network of 558 prioritized protein-protein interactions between trophectodermal, epithelial and stromal cells, which were grouped into clusters based on the function of the involved molecules. The role of galectins (LGALS1 and LGALS3), integrin ß1 (ITGB1), basigin (BSG) and osteopontin (SPP1) in the embryo implantation process were highlighted. LARGE SCALE DATA: RNA-seq data are available at www.ncbi.nlm.nih.gov/geo under accession number GSE97929. LIMITATIONS REASONS FOR CAUTION: Providing a static snap-shot of a dynamic process and the nature of prediction analysis is limited to the known interactions available in databases. Furthermore, the cell sorting technique used separated enriched epithelial cells and stromal cells but did not separate luminal from glandular epithelium. Also, the use of biopsies taken from non-pregnant women and using spare IVF embryos (due to ethical considerations) might miss some of the critical interactions characteristic of natural conception only. WIDER IMPLICATIONS OF THE FINDINGS: The findings of our study provide new insights into the molecular embryo-endometrium interplay in the first steps of implantation process in humans. Knowledge about the endometrial cell-type-specific molecules that coordinate successful implantation is vital for understanding human reproduction and the underlying causes of implantation failure and infertility. Our study results provide a useful resource for future reproductive research, allowing the exploration of unknown mechanisms of implantation. We envision that those studies will help to improve the understanding of the complex embryo implantation process, and hopefully generate new prognostic and diagnostic biomarkers and therapeutic approaches to target both infertility and fertility, in the form of new contraceptives. STUDY FUNDING/COMPETING INTERESTS: This research was funded by the Estonian Research Council (grant PRG1076); Horizon 2020 innovation grant (ERIN, grant no. EU952516); Enterprise Estonia (grant EU48695); the EU-FP7 Marie Curie Industry-Academia Partnerships and Pathways (IAPP, grant SARM, EU324509); Spanish Ministry of Economy, Industry and Competitiveness (MINECO) and European Regional Development Fund (FEDER) (grants RYC-2016-21199, ENDORE SAF2017-87526-R, and Endo-Map PID2021-127280OB-100); Programa Operativo FEDER Andalucía (B-CTS-500-UGR18; A-CTS-614-UGR20), Junta de Andalucía (PAIDI P20_00158); Margarita Salas program for the Requalification of the Spanish University system (UJAR01MS); the Knut and Alice Wallenberg Foundation (KAW 2015.0096); Swedish Research Council (2012-2844); and Sigrid Jusélius Foundation; Academy of Finland. A.S.-L. is funded by the Spanish Ministry of Science, Innovation and Universities (PRE2018-085440). K.G.-D. has received consulting fees and/or honoraria from RemovAid AS, Norway Bayer, MSD, Gedeon Richter, Mithra, Exeltis, MedinCell, Natural cycles, Exelgyn, Vifor, Organon, Campus Pharma and HRA-Pharma and NIH support to the institution; D.B. is an employee of IGENOMIX. The rest of the authors declare no conflict of interest.

Genome-wide association study identifies five risk loci for pernicious anemia.

Pernicious anemia is a rare condition characterized by vitamin B12 deficiency anemia due to lack of intrinsic factor, often caused by autoimmune gastritis. Patients with pernicious anemia have a higher incidence of other autoimmune disorders, such as type 1 diabetes, vitiligo, and autoimmune thyroid issues. Therefore, the disease has a clear autoimmune basis, although the genetic susceptibility factors have thus far remained poorly studied. We conduct a genome-wide association study meta-analysis in 2166 cases and 659,516 European controls from population-based biobanks and identify genome-wide significant signals in or near the PTPN22 (rs6679677, p = 1.91 × 10-24, OR = 1.63), PNPT1 (rs12616502, p = 3.14 × 10-8, OR = 1.70), HLA-DQB1 (rs28414666, p = 1.40 × 10-16, OR = 1.38), IL2RA (rs2476491, p = 1.90 × 10-8, OR = 1.22) and AIRE (rs74203920, p = 2.33 × 10-9, OR = 1.83) genes, thus providing robust associations between pernicious anemia and genetic risk factors.

In vitro fertilization does not increase the incidence of de novo copy number alterations in fetal and placental lineages.

Although chromosomal instability (CIN) is a common phenomenon in cleavage-stage embryogenesis following in vitro fertilization (IVF)1-3, its rate in naturally conceived human embryos is unknown. CIN leads to mosaic embryos that contain a combination of genetically normal and abnormal cells, and is significantly higher in in vitro-produced preimplantation embryos as compared to in vivo-conceived preimplantation embryos4. Even though embryos with CIN-derived complex aneuploidies may arrest between the cleavage and blastocyst stages of embryogenesis5,6, a high number of embryos containing abnormal cells can pass this strong selection barrier7,8. However, neither the prevalence nor extent of CIN during prenatal development and at birth, following IVF treatment, is well understood. Here we profiled the genomic landscape of fetal and placental tissues postpartum from both IVF and naturally conceived children, to investigate the prevalence and persistence of large genetic aberrations that probably arose from IVF-related CIN. We demonstrate that CIN is not preserved at later stages of prenatal development, and that de novo numerical aberrations or large structural DNA imbalances occur at similar rates in IVF and naturally conceived live-born neonates. Our findings affirm that human IVF treatment has no detrimental effect on the chromosomal constitution of fetal and placental lineages.

MUC20 expression marks the receptive phase of the human endometrium.

RESEARCH QUESTION: How does mucin MUC20 expression change during the menstrual cycle in different cell types of human endometrium? DESIGN: Study involved examination of MUC20 expression in two previously published RNA-seq datasets in whole endometrial tissue (n = 10), sorted endometrial epithelial (n = 44) or stromal (n = 42) cell samples. RNA-Seq results were validated by quantitative reverse transcription polymerase chain reaction (qRT-PCR) in whole tissue (n = 10), sorted epithelial (n = 17) and stromal (n = 17) cell samples. MUC20 protein localization and expression were analysed in human endometrium by immunohistochemical analysis of intact endometrial tissue (n = 6) and also Western blot of cultured stromal and epithelial cells (n = 2). RESULTS: MUC20 is differentially expressed in the endometrium between the pre-receptive and receptive phases. We show that MUC20 is predominantly expressed by epithelial cells of the receptive endometrium, both at the mRNA (RNA-Seq, P = 0.005; qRT-PCR, P = 0.039) and protein levels (Western blot; immunohistochemistry, P = 0.029). CONCLUSION: Our results indicate MUC20 as a novel marker of mid-secretory endometrial biology. We propose a model of MUC20 function in the hepatocyte growth factor (HGF)-activated mesenchymal-epithelial transition (MET) receptor signalling specifically in the receptive phase. Further investigations should reveal the precise function of MUC20 in human endometrium and the possible connection between MUC20 and HGF-activated MET receptor signalling. MUC20 could potentially be included in the list of endometrial receptivity markers after further clinical validation.

Hyperglycosylated hCG activates LH/hCG-receptor with lower activity than hCG.

While human chorionic gonadotropin (hCG) appears to have an essential role in early pregnancy, it is controversial whether the hyperglycosylated form of hCG (hCG-h), which is the major hCG isoform during the first 4-5 weeks of pregnancy, is able to activate LH/hCG receptor (LHCGR). To address this, we utilized different extensively characterized hCG and hCGß reference reagents, cell culture- and urine-derived hCG-h preparations, and an in vitro reporter system for LHCGR activation. The WHO hCG reference reagent (99/688) was found to activate LHCGR with an EC50-value of 3.3 ±â€¯0.6 pmol/L (n = 9). All three studied hCG-h preparations were also able to activate LHCGR, but with a lower potency (EC50-values between 7.1 ±â€¯0.5 and 14 ±â€¯3 pmol/L, n = 5-11, for all P < 0.05 as compared to the hCG reference). The activities of commercial urinary hCG (Pregnyl) and recombinant hCG (Ovitrelle) preparations were intermediate between those of the hCG reference and the hCG-h. These results strongly suggest that the hCG-h is functionally similar to hCG, although it has lower potency for LHCGR activation. Whether this explains the reduced proportion of hCG-h to hCG reported in patients developing early onset pre-eclampsia or those having early pregnancy loss remains to be determined.

TAC-seq: targeted DNA and RNA sequencing for precise biomarker molecule counting.

Targeted next-generation sequencing (NGS) methods have become essential in medical research and diagnostics. In addition to NGS sensitivity and high-throughput capacity, precise biomolecule counting based on unique molecular identifier (UMI) has potential to increase biomolecule detection accuracy. Although UMIs are widely used in basic research its introduction to clinical assays is still in progress. Here, we present a robust and cost-effective TAC-seq (Targeted Allele Counting by sequencing) method that uses UMIs to estimate the original molecule counts of mRNAs, microRNAs, and cell-free DNA. We applied TAC-seq in three different clinical applications and compared the results with standard NGS. RNA samples extracted from human endometrial biopsies were analyzed using previously described 57 mRNA-based receptivity biomarkers and 49 selected microRNAs at different expression levels. Cell-free DNA aneuploidy testing was based on cell line (47,XX, +21) genomic DNA. TAC-seq mRNA profiling showed identical clustering results to transcriptome RNA sequencing, and microRNA detection demonstrated significant reduction in amplification bias, allowing to determine minor expression changes between different samples that remained undetermined by standard NGS. The mimicking experiment for cell-free DNA fetal aneuploidy analysis showed that TAC-seq can be applied to count highly fragmented DNA, detecting significant (p = 7.6 × 10-4) excess of chromosome 21 molecules at 10% fetal fraction level. Based on three proof-of-principle applications we demonstrate that TAC-seq is an accurate and highly potential biomarker profiling method for advanced medical research and diagnostics.

Endometrial receptivity revisited: endometrial transcriptome adjusted for tissue cellular heterogeneity.

STUDY QUESTION: Does cellular composition of the endometrial biopsy affect the gene expression profile of endometrial whole-tissue samples? SUMMARY ANSWER: The differences in epithelial and stromal cell proportions in endometrial biopsies modify the whole-tissue gene expression profiles and affect the results of differential expression analyses. WHAT IS ALREADY KNOWN: Each 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 method and tools used to obtain the tissue biopsy. STUDY DESIGN, SIZE, DURATION: Using 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, METHODS: Paired early- and mid-secretory endometrial biopsies were obtained from 35 healthy, regularly cycling, fertile volunteers, aged 23-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 the proportions of epithelial and stromal cells in the endometrial biopsy. The results were then compared to conventional differential expression analysis. Ten paired endometrial samples were analysed with qPCR to validate the results. MAIN RESULTS AND THE ROLE OF CHANCE: The estimated average proportions of stromal and epithelial cells in early secretory phase were 65% and 35%, and during mid-secretory phase, 46% and 54%, respectively, correlating 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 = 2645) 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. LARGE-SCALE DATA: The RNA-seq data presented in this study is deposited in the Gene Expression Omnibus database with accession number GSE98386. LIMITATIONS REASONS FOR CAUTION: Only 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 FINDINGS: The 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. STUDY FUNDING/COMPETING INTEREST(S): This study was funded by: Estonian Ministry of Education and Research (grant IUT34-16); Enterprise Estonia (EU48695); the EU-FP7 Eurostars program (NOTED, EU41564); the EU-FP7 Marie Curie Industry-Academia Partnerships and Pathways (SARM, EU324509); Horizon 2020 innovation program (WIDENLIFE, EU692065); MSCA-RISE-2015 project MOMENDO (No 691058) and the Miguel Servet Program Type I of Instituto de Salud Carlos III (CP13/00038); Spanish Ministry of Economy, Industry and Competitiveness (MINECO) and European Regional Development Fund (FEDER): grants RYC-2016-21199 and ENDORE SAF2017-87526. Authors confirm no competing interests.

Meta-signature of human endometrial receptivity: a meta-analysis and validation study of transcriptomic biomarkers.

Previous transcriptome studies of the human endometrium have revealed hundreds of simultaneously up- and down-regulated genes that are involved in endometrial receptivity. However, the overlap between the studies is relatively small, and we are still searching for potential diagnostic biomarkers. Here we perform a meta-analysis of endometrial-receptivity associated genes on 164 endometrial samples (76 from 'pre-receptive' and 88 from mid-secretory, 'receptive' phase endometria) using a robust rank aggregation (RRA) method, followed by enrichment analysis, and regulatory microRNA prediction. We identify a meta-signature of endometrial receptivity involving 57 mRNA genes as putative receptivity markers, where 39 of these we confirm experimentally using RNA-sequencing method in two separate datasets. The meta-signature genes highlight the importance of immune responses, the complement cascade pathway and the involvement of exosomes in mid-secretory endometrial functions. Bioinformatic prediction identifies 348 microRNAs that could regulate 30 endometrial-receptivity associated genes, and we confirm experimentally the decreased expression of 19 microRNAs with 11 corresponding up-regulated meta-signature genes in our validation experiments. The 57 identified meta-signature genes and involved pathways, together with their regulatory microRNAs could serve as promising and sought-after biomarkers of endometrial receptivity, fertility and infertility.

Optimizing bone morphogenic protein 4-mediated human embryonic stem cell differentiation into trophoblast-like cells using fibroblast growth factor 2 and transforming growth factor-ß/activin/nodal signalling inhibition.

Several studies have demonstrated that human embryonic stem cells (hESC) can be differentiated into trophoblast-like cells if exposed to bone morphogenic protein 4 (BMP4) and/or inhibitors of fibroblast growth factor 2 (FGF2) and the transforming growth factor beta (TGF-ß)/activin/nodal signalling pathways. The goal of this study was to investigate how the inhibitors of these pathways improve the efficiency of hESC differentiation when compared with basic BMP4 treatment. RNA sequencing was used to analyse the effects of all possible inhibitor combinations on the differentiation of hESC into trophoblast-like cells over 12 days. Genes differentially expressed compared with untreated cells were identified at seven time points. Additionally, expression of total human chorionic gonadotrophin (HCG) and its hyperglycosylated form (HCG-H) were determined by immunoassay from cell culture media. We showed that FGF2 inhibition with BMP4 activation up-regulates syncytiotrophoblast-specific genes (CGA, CGB and LGALS16), induces several molecular pathways involved in embryo implantation and triggers HCG-H production. In contrast, inhibition of the TGF-ß/activin/nodal pathway decreases the ability of hESC to form trophoblast-like cells. Information about the conditions needed for hESC differentiation toward trophoblast-like cells helps us to find an optimal model for studying the early development of human trophoblasts in normal and in complicated pregnancy.

Globin mRNA reduction for whole-blood transcriptome sequencing.

The transcriptome analysis of whole-blood RNA by sequencing holds promise for the identification and tracking of biomarkers; however, the high globin mRNA (gmRNA) content of erythrocytes hampers whole-blood and buffy coat analyses. We introduce a novel gmRNA locking assay (GlobinLock, GL) as a robust and simple gmRNA reduction tool to preserve RNA quality, save time and cost. GL consists of a pair of gmRNA-specific oligonucleotides in RNA initial denaturation buffer that is effective immediately after RNA denaturation and adds only ten minutes of incubation to the whole cDNA synthesis procedure when compared to non-blood RNA analysis. We show that GL is fully effective not only for human samples but also for mouse and rat, and so far incompletely studied cow, dog and zebrafish.