Dynamic genome-wide gene expression and immune cell composition in the developing human placenta

Differential gene expression analysis pipelines and bioinformatic tools for the identification of specific biomarkers: A review

In recent years, the role of bioinformatics and computational biology together with omics techniques and transcriptomics has gained tremendous importance in biomedicine and healthcare, particularly for the identification of biomarkers for precision medicine and drug discovery. Differential gene expression (DGE) analysis is one of the most used techniques for RNA-sequencing (RNA-seq) data analysis. This tool, which is typically used in various RNA-seq data processing applications, allows the identification of differentially expressed genes across two or more sample sets. Functional enrichment analyses can then be performed to annotate and contextualize the resulting gene lists. These studies provide valuable information about disease-causing biological processes and can help in identifying molecular targets for novel therapies. This review focuses on differential gene expression (DGE) analysis pipelines and bioinformatic techniques commonly used to identify specific biomarkers and discuss the advantages and disadvantages of these techniques.

Placental homogeneity: Characterizing transcriptional variation among equine chorioallantoic locations

The proper function of the placenta is essential for the health and growth of the fetus and the mother. The placenta relies on dynamic gene expression for its correct and timely development and function. Although numerous studies have identified genes vital for placental functions, equine placental molecular research has primarily focused on single placental locations, in sharp contrast with the broader approach in human studies. Here, we hypothesized that the molecular differences across different regions of the equine placenta are negligible because of its diffuse placental type with a macroscopic homogenous distribution of villi across the placental surface. We compared the transcriptome and stereological findings of the body, pregnant horn, and non-pregnant horn within the equine chorioallantois. Our transcriptomic analysis indicates that the variation between regions of the placenta within individuals is less than the variation observed between individuals. A low number of differentially expressed genes (DEGs) (n = 8) was identified when comparing pregnant and non-pregnant horns within the same placenta, suggesting a remarkable molecular uniformity. A higher number of DEGs was identified when comparing each horn to the body (193 DEGs comparing pregnant horn with body and 207 DEGs comparing non-pregnant horn with body). Genes with a higher expression in the body were associated with processes such as extracellular matrix synthesis and remodeling, which is relevant for placental maturation and placenta-endometrial separation at term and implies asynchrony of these processes across locations. The stereological analysis showed no differences in microcotyledonary density, and width between the locations. However, we observed a greater chorioallantoic thickness in the body and pregnant horn compared to the non-pregnant horn. Overall, our findings reveal a uniform transcriptomic profile across the placental horns, alongside a more distinct gene expression pattern between the uterine body and horns. These regional differences in gene expression suggest a different pace in the placental maturation and detachment among the placental locations.

Revealing the molecular landscape of human placenta: a systematic review and meta-analysis of single-cell RNA sequencing studies

BACKGROUND

With increasing significance of developmental programming effects associated with placental dysfunction, more investigations are devoted to improving the characterization and understanding of placental signatures in health and disease. The placenta is a transitory but dynamic organ adapting to the shifting demands of fetal development and available resources of the maternal supply throughout pregnancy. Trophoblasts (cytotrophoblasts, syncytiotrophoblasts, and extravillous trophoblasts) are placental-specific cell types responsible for the main placental exchanges and adaptations. Transcriptomic studies with single-cell resolution have led to advances in understanding the placenta's role in health and disease. These studies, however, often show discrepancies in characterization of the different placental cell types.

OBJECTIVE AND RATIONALE

We aim to review the knowledge regarding placental structure and function gained from the use of single-cell RNA sequencing (scRNAseq), followed by comparing cell-type-specific genes, highlighting their similarities and differences. Moreover, we intend to identify consensus marker genes for the various trophoblast cell types across studies. Finally, we will discuss the contributions and potential applications of scRNAseq in studying pregnancy-related diseases.

SEARCH METHODS

We conducted a comprehensive systematic literature review to identify different cell types and their functions at the human maternal-fetal interface, focusing on all original scRNAseq studies on placentas published before March 2023 and published reviews (total of 28 studies identified) using PubMed search. Our approach involved curating cell types and subtypes that had previously been defined using scRNAseq and comparing the genes used as markers or identified as potential new markers. Next, we reanalyzed expression matrices from the six available scRNAseq raw datasets with cell annotations (four from first trimester and two at term), using Wilcoxon rank-sum tests to compare gene expression among studies and annotate trophoblast cell markers in both first trimester and term placentas. Furthermore, we integrated scRNAseq raw data available from 18 healthy first trimester and nine term placentas, and performed clustering and differential gene expression analysis. We further compared markers obtained with the analysis of annotated and raw datasets with the literature to obtain a common signature gene list for major placental cell types.

OUTCOMES

Variations in the sampling site, gestational age, fetal sex, and subsequent sequencing and analysis methods were observed between the studies. Although their proportions varied, the three trophoblast types were consistently identified across all scRNAseq studies, unlike other non-trophoblast cell types. Notably, no marker genes were shared by all studies for any of the investigated cell types. Moreover, most of the newly defined markers in one study were not observed in other studies. These discrepancies were confirmed by our analysis on trophoblast cell types, where hundreds of potential marker genes were identified in each study but with little overlap across studies. From 35 461 and 23 378 cells of high quality in the first trimester and term placentas, respectively, we obtained major placental cell types, including perivascular cells that previously had not been identified in the first trimester. Importantly, our meta-analysis provides marker genes for major placental cell types based on our extensive curation.

WIDER IMPLICATIONS

This review and meta-analysis emphasizes the need for establishing a consensus for annotating placental cell types from scRNAseq data. The marker genes identified here can be deployed for defining human placental cell types, thereby facilitating and improving the reproducibility of trophoblast cell annotation.

High-throughput mRNA-seq atlas of human placenta shows vast transcriptome remodeling from first to third trimester†

The placenta, composed of chorionic villi, changes dramatically across gestation. Understanding differences in ongoing pregnancies are essential to identify the role of chorionic villi at specific times in gestation and develop biomarkers and prognostic indicators of maternal-fetal health. The normative mRNA profile is established using next-generation sequencing of 124 first trimester and 43 third trimester human placentas from ongoing healthy pregnancies. Stably expressed genes (SEGs) not different between trimesters and with low variability are identified. Differential expression analysis of first versus third trimester adjusted for fetal sex is performed, followed by a subanalysis with 23 matched pregnancies to control for subject variability using the same genetic and environmental background. Placenta expresses 14,979 polyadenylated genes above sequencing noise (transcripts per million > 0.66), with 10.7% SEGs across gestation. Differentially expressed genes (DEGs) account for 86.7% of genes in the full cohort [false discovery rate (FDR) < 0.05]. Fold changes highly correlate between the full cohort and subanalysis (Pearson = 0.98). At stricter thresholds (FDR < 0.001, fold change > 1.5), there remains 50.1% DEGs (3353 upregulated in first and 4155 upregulated in third trimester). This is the largest mRNA atlas of healthy human placenta across gestation, controlling for genetic and environmental factors, demonstrating substantial changes from first to third trimester in chorionic villi. Specific differences and SEGs may be used to understand the specific role of the chorionic villi throughout gestation and develop first trimester biomarkers of placental health that transpire across gestation, which can be used for future development of biomarkers for maternal-fetal health.

Genome-wide identification of transcriptional enhancers during human placental development and association with function, differentiation, and disease†

The placenta is a dynamic organ that must perform a remarkable variety of functions during its relatively short existence in order to support a developing fetus. These functions include nutrient delivery, gas exchange, waste removal, hormone production, and immune barrier protection. Proper placenta development and function are critical for healthy pregnancy outcomes, but the underlying genomic regulatory events that control this process remain largely unknown. We hypothesized that mapping sites of transcriptional enhancer activity and associated changes in gene expression across gestation in human placenta tissue would identify genomic loci and predicted transcription factor activity related to critical placenta functions. We used a suite of genomic assays [i.e., RNA-sequencing (RNA-seq), Precision run-on-sequencing (PRO-seq), and Chromatin immunoprecipitation-sequencing (ChIP-seq)] and computational pipelines to identify a set of >20 000 enhancers that are active at various time points in gestation. Changes in the activity of these enhancers correlate with changes in gene expression. In addition, some of these enhancers encode risk for adverse pregnancy outcomes. We further show that integrating enhancer activity, transcription factor motif analysis, and transcription factor expression can identify distinct sets of transcription factors predicted to be more active either in early pregnancy or at term. Knockdown of selected identified transcription factors in a trophoblast stem cell culture model altered the expression of key placental marker genes. These observations provide a framework for future mechanistic studies of individual enhancer-transcription factor-target gene interactions and have the potential to inform genetic risk prediction for adverse pregnancy outcomes.

Placental single cell transcriptomics: Opportunities for endocrine disrupting chemical toxicology

The placenta performs essential biologic functions for fetal development throughout pregnancy. Placental dysfunction is at the root of multiple adverse birth outcomes such as intrauterine growth restriction, preeclampsia, and preterm birth. Exposure to endocrine disrupting chemicals during pregnancy can cause placental dysfunction, and many prior human studies have examined molecular changes in bulk placental tissues. Placenta-specific cell types, including cytotrophoblasts, syncytiotrophoblasts, extravillous trophoblasts, and placental resident macrophage Hofbauer cells play unique roles in placental development, structure, and function. Toxicant-induced changes in relative abundance and/or impairment of these cell types likely contribute to placental pathogenesis. Although gene expression insights gained from bulk placental tissue RNA-sequencing data are useful, their interpretation is limited because bulk analysis can mask the effects of a chemical on individual populations of placental cells. Cutting-edge single cell RNA-sequencing technologies are enabling the investigation of placental cell-type specific responses to endocrine disrupting chemicals. Moreover, in situ bioinformatic cell deconvolution enables the estimation of cell type proportions in bulk placental tissue gene expression data. These emerging technologies have tremendous potential to provide novel mechanistic insights in a complex heterogeneous tissue with implications for toxicant contributions to adverse pregnancy outcomes.

The Placenta: A Maternofetal Interface

The placenta is the gatekeeper between the mother and the fetus. Over the first trimester of pregnancy, the fetus is nourished by uterine gland secretions in a process known as histiotrophic nutrition. During the second trimester of pregnancy, placentation has evolved to the point at which nutrients are delivered to the placenta via maternal blood (hemotrophic nutrition). Over gestation, the placenta must adapt to these variable nutrient supplies, to alterations in maternal physiology and blood flow, and to dynamic changes in fetal growth rates. Numerous questions remain about the mechanisms used to transport nutrients to the fetus and the maternal and fetal determinants of this process. Growing data highlight the ability of the placenta to regulate this process. As new technologies and omics approaches are utilized to study this maternofetal interface, greater insight into this unique organ and its impact on fetal development and long-term health has been obtained.

High-throughput mRNA-seq atlas of human placenta shows vast transcriptome remodeling from first to third trimester

Background

The placenta, composed of chorionic villi, changes dramatically across gestation. Understanding differences in ongoing pregnancies are essential to identify the role of chorionic villi at specific times in gestation and develop biomarkers and prognostic indicators of maternal- fetal health.

Methods

The normative mRNA profile is established using next-generation sequencing of 124 first trimester and 43 third trimester human placentas from ongoing healthy pregnancies. Stably expressed genes not different between trimesters and with low variability are identified. Differential expression analysis of first versus third trimester adjusted for fetal sex is performed, followed by a subanalysis with 23 matched pregnancies to control for subject variability using the same genetic and environmental background.

Results

Placenta expresses 14,979 mRNAs above sequencing noise (TPM>0.66), with 1,545 stably expressed genes across gestation. Differentially expressed genes account for 86.7% of genes in the full cohort (FDR<0.05). Fold changes highly correlate between the full cohort and subanalysis (Pearson = 0.98). At stricter thresholds (FDR<0.001, fold change>1.5), there are 6,941 differentially expressed protein coding genes (3,206 upregulated in first and 3,735 upregulated in third trimester).

Conclusion

This is the largest mRNA atlas of healthy human placenta across gestation, controlling for genetic and environmental factors, demonstrating substantial changes from first to third trimester in chorionic villi. Specific differences and stably expressed genes may be used to understand the specific role of the chorionic villi throughout gestation and develop first trimester biomarkers of placental health that transpire across gestation, which can be used for future development of biomarkers in maternal-fetal disease.

Placental cell type deconvolution reveals that cell proportions drive preeclampsia gene expression differences

The placenta mediates adverse pregnancy outcomes, including preeclampsia, which is characterized by gestational hypertension and proteinuria. Placental cell type heterogeneity in preeclampsia is not well-understood and limits mechanistic interpretation of bulk gene expression measures. We generated single-cell RNA-sequencing samples for integration with existing data to create the largest deconvolution reference of 19 fetal and 8 maternal cell types from placental villous tissue (n = 9 biological replicates) at term (n = 40,494 cells). We deconvoluted eight published microarray case-control studies of preeclampsia (n = 173 controls, 157 cases). Preeclampsia was associated with excess extravillous trophoblasts and fewer mesenchymal and Hofbauer cells. Adjustment for cellular composition reduced preeclampsia-associated differentially expressed genes (log2 fold-change cutoff = 0.1, FDR < 0.05) from 1154 to 0, whereas downregulation of mitochondrial biogenesis, aerobic respiration, and ribosome biogenesis were robust to cell type adjustment, suggesting direct changes to these pathways. Cellular composition mediated a substantial proportion of the association between preeclampsia and FLT1 (37.8%, 95% CI [27.5%, 48.8%]), LEP (34.5%, 95% CI [26.0%, 44.9%]), and ENG (34.5%, 95% CI [25.0%, 45.3%]) overexpression. Our findings indicate substantial placental cellular heterogeneity in preeclampsia contributes to previously observed bulk gene expression differences. This deconvolution reference lays the groundwork for cellular heterogeneity-aware investigation into placental dysfunction and adverse birth outcomes.

Identifying novel regulators of placental development using time-series transcriptome data

The placenta serves as a connection between the mother and the fetus during pregnancy, providing the fetus with oxygen, nutrients, and growth hormones. However, the regulatory mechanisms and dynamic gene interaction networks underlying early placental development are understudied. Here, we generated RNA-sequencing data from mouse fetal placenta at embryonic days 7.5, 8.5, and 9.5 to identify genes with timepoint-specific expression, then inferred gene interaction networks to analyze highly connected network modules. We determined that timepoint-specific gene network modules were associated with distinct developmental processes, and with similar expression profiles to specific human placental cell populations. From each module, we identified hub genes and their direct neighboring genes, which were predicted to govern placental functions. We confirmed that four novel candidate regulators identified through our analyses regulate cell migration in the HTR-8/SVneo cell line. Overall, we predicted several novel regulators of placental development expressed in specific placental cell types using network analysis of bulk RNA-sequencing data. Our findings and analysis approaches will be valuable for future studies investigating the transcriptional landscape of early development.