- Research article
- Open Access
Amniotic fluid cell-free transcriptome: a glimpse into fetal development and placental cellular dynamics during normal pregnancy
BMC Medical Genomics volume 13, Article number: 25 (2020)
The amniotic fluid (AF) cell-free transcriptome is modulated by physiologic and pathologic processes during pregnancy. AF gene expression changes with advancing gestation reflect fetal development and organ maturation; yet, defining normal expression and splicing patterns for biomarker discovery in obstetrics requires larger heterogeneous cohorts, evaluation of potential confounding factors, and novel analytical approaches.
Women with a normal pregnancy who had an AF sample collected during midtrimester (n = 30) or at term gestation (n = 68) were included. Expression profiling at exon level resolution was performed using Human Transcriptome Arrays. Differential expression was based on moderated t-test adjusted p < 0.05 and fold change > 1.25; for differential splicing, a splicing index > 2 and adjusted p < 0.05 were required. Functional profiling was used to interpret differentially expressed or spliced genes. The expression of tissue-specific and cell-type specific signatures defined by single-cell genomics was quantified and correlated with covariates. In-silico validation studies were performed using publicly available datasets.
1) 64,071 genes were detected in AF, with 11% of the coding and 6% of the non-coding genes being differentially expressed between midtrimester and term gestation. Expression changes were highly correlated with those previously reported (R > 0.79, p < 0.001) and featured increased expression of genes specific to the trachea, salivary glands, and lung and decreased expression of genes specific to the cardiac myocytes, uterus, and fetal liver, among others. 2) Single-cell RNA-seq signatures of the cytotrophoblast, Hofbauer cells, erythrocytes, monocytes, T and B cells, among others, showed complex patterns of modulation with gestation (adjusted p < 0.05). 3) In 17% of the genes detected, we found differential splicing with advancing gestation in genes related to brain development processes and immunity pathways, including some that were missed based on differential expression analysis alone.
This represents the largest AF transcriptomics study in normal pregnancy, reporting for the first time that single-cell genomic signatures can be tracked in the AF and display complex patterns of expression during gestation. We also demonstrate a role for alternative splicing in tissue-identity acquisition, organ development, and immune processes. The results herein may have implications for the development of fetal testing to assess placental function and fetal organ maturity.
Amniotic fluid (AF) provides mechanical protection to the developing fetus and has important nutritional and immunologic roles [1,2,3]. Given that AF is in direct contact with the placenta and fetal membranes, surrounding the fetus, while passing through several fetal cavities (e.g. gastrointestinal and respiratory tracts), it is expected that its molecular composition is both reflective of and contributes to fetal wellbeing [4,5,6].
For decades, AF sampled through transabdominal amniocentesis has been used to assess the fetal karyotype [7,8,9,10,11], lung maturity [12,13,14,15], and presence of inflammatory conditions, such as intra-amniotic infection and sterile intra-amniotic inflammation [16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82]. Although the ultimate goal is to develop fetal testing via non-invasive sampling, such as urine or maternal blood, the advantages of AF for biomarker discovery have often been recognized [83,84,85,86,87].
Among the omics platforms used to study AF (see Kamath-Rayne et al.  for a review), the analysis of cell-free mRNA (cfRNA) has the advantage of being easier to profile than its proteomics [45, 89, 90] and metabolomics [91,92,93,94,95,96] counterparts. The AF cfRNAs are thought to be contributed directly by the fetus and by apoptotic amniocytes  and have been shown to be altered by physiologic and pathologic factors such as gestational age [83, 85, 98], fetal sex , maternal obesity , genetic syndromes [100,101,102], and neonatal co-morbidities  (see Zwemer and Bianchi for a review ). Of all the factors previously reported to be reflected in the AF transcriptome, advancing gestation seems to have the most dramatic effect on the AF transcriptome considering the number of genes differentially expressed. AF gene expression changes with gestational age have been associated with cell types found in the intrauterine environment and with the development of multiple organ systems [78, 85, 98, 103]. These results point to the possible use of the AF transcriptome to complement fetal lung maturity evaluation [85, 98] when elective delivery prior to term is considered, and also to discover biomarkers for the ‘great obstetrical syndromes’ .
Currently, several limitations exist in defining a reference of gene expression and splicing patterns during gestation in normal pregnancy based on publicly available data. Chief among them is the rather limited sample size and use of measurement platforms that do not allow assessment of non-coding RNAs and splicing patterns. Moreover, while previous studies considered multiple functional databases and tissue-specific gene sets to interpret differentially expressed genes in AF [85, 98], detailed signatures of specific cell-types were not readily available. Recently, single-cell genomics studies enabled the assessment of the maternal-fetal crosstalk by quantifying mRNA signatures specific to sub-populations of placental cells [105,106,107]. These mRNA signatures were shown to be detectable by cell-free  and cellular  transcriptome analyses of maternal blood; however, to date, they have not been evaluated in AF.
Therefore, the current study aimed i) to evaluate the effect of gestational age on AF cfRNA expression and splicing while considering relevant potential covariates (fetal sex, maternal characteristics, mode of sample collection, and indication for amniocentesis), ii) to determine whether the molecular dynamics of specific placental cell populations, such as those recently described by single-cell RNA sequencing , can be tracked by AF cfRNA analysis, and iii) to assess the extent to which previously reported changes in the overall gene expression with gestational age, fetal sex, and maternal obesity can be replicated, given differences in populations and profiling techniques.
To conduct a prospective longitudinal study, we enrolled pregnant women attending the Center for Advanced Obstetrical Care and Research of the Perinatology Research Branch, Eunice Kennedy Shriver National Institute of Child Health and Human Development (NICHD), National Institutes of Health, U.S. Department of Health and Human Services; Wayne State University School of Medicine; and Hutzel Women’s Hospital of the Detroit Medical Center (Detroit, Michigan, USA). Based on this cohort, we designed a retrospective study to include 30 women who underwent transabdominal amniocentesis during the midtrimester at 16.4–24.0 weeks of gestation (median = 21.1 weeks) to assess the fetal karyotype or to rule out intra-amniotic inflammation/infection. In all cases tested, the karyotype was normal. The study also included 68 women at term not in labor (TNL) who had an AF sample collected either by transabdominal amniocentesis (to assess fetal lung maturity) or for research purposes during Cesarean section at 37.1–40.9 weeks of gestation (median = 39.0 weeks). All study participants included herein delivered at term and had a normal singleton pregnancy with appropriate fetal growth, according to the INTERGROWTH-21st birthweight standard , which was previously found to match the population of patients attending our research clinic . Five ml of AF were collected from each woman and processed according to the recommended protocol .
Starting with 5 ml of AF, we applied the Plasma/Serum RNA Purification Maxi Kit (#56200; Norgen Biotek Corp., Thorold, Ontario, Canada), including the optional DNAse treatment, according to the manufacturer’s protocol. Following RNA extraction, each sample was concentrated to a volume of 12 μl by using the RNA Clean & Concentrator-5 Kit (#R1015; Zymo Research, Irvine, California, USA). The concentrates were then quantified by UV spectrophotometry on a DropSense 96 system (PerkinElmer, Waltham, MA, USA) and then quality assessed on the Agilent 2200 TapeStation system (Agilent Technologies, Santa Clara, California, USA).
Ten ng of RNA was reverse-transcribed and amplified using the Affymetrix GeneChip® WT Pico Reagent Kit (Affymetrix, Inc., Santa Clara, California, USA), following the manufacturer’s suggested protocol. Briefly, 5.5 μg of sense-stranded cDNA was fragmented, labeled, and hybridized in a final volume of 200 μl to the Affymetrix GeneChip® Human Transcriptome Array 2.0 in an Affymetrix hybridization oven at 45 °C at 60 rpm for 16 h. Wash and stain steps were performed utilizing an Affymetrix GeneChip® Fluidics Station 450 and scanned on an Affymetrix GeneChip® Scanner 3000. Raw intensity data were generated from array images using Affymetrix GeneChip™ Command Console Software.
Raw gene expression data of more than 6.0 million microarray probes and 98 microarrays were preprocessed (background correction, normalization, and summarization) using the Robust Multi-array Average (RMA) approach  implemented in the oligo package . Expression summarization was obtained into one value per sample and transcript cluster (gene level data) for differential expression analysis as well as at the level of each exon or exon-exon junction (probe set level data) for differential splicing analysis. Transcript clusters were assigned to unique genes using annotation from the hta20transcriptcluster.db package of Bioconductor . Only genes expressed above background (p < 0.05) in at least 25% of the samples of either group (midtrimester or TNL) were retained for further differential expression and differential splicing analyses. Detection p-values were obtained using the Expression Console, version 1.4. (Affymetrix, Inc.).
Gene level log2 expression data were analyzed using linear models implemented in the limma  package in Bioconductor . The gestational-age effect (TNL versus midtrimester or early midtrimester versus late midtrimester) was assessed while adjusting only for covariates with significant effect on gene expression among those considered (fetal sex, maternal race, obesity, parity, smoking status, and mode of sample collection). The effect of fetal sex was assessed while adjusting for gestational age. Expression changes were deemed significant based on false discovery rate q-values < 0.05 and a minimum fold change of 1.25-fold, which are known to be rather conservative for this microarray platform .
Differential splicing was assessed in relation to gestational age, fetal sex, and obesity based on the splicing index (SI) method  implemented in the Transcriptome Analysis Console (TAC) version 4.0 (Affymetrix, Inc.) using netaffx_release_36 gene annotations. The splicing index represents the difference in average exon usage between groups, where exon usage is defined as exon level expression relative to overall gene abundance in a given sample. An adjusted p < 0.05 and SI > 2.0-fold for one exon/exon-exon junction of a gene was required to infer significant alternative splicing for that gene. A second differential splicing algorithm, also based on the SI concept, was applied: DiffSplice, implemented in the limma package , uses an F-test to assess whether log fold changes (between groups) differ among exons of the same gene. Furthermore, based on the SI and corresponding p-values, specific patterns of differential splicing were identified among the following: cassette exon, mutually exclusive exons, alternative 5′ sites, alternative 3′ sites, and intron retention (see Blencowe B. for a review ), implementing the Exon Event Estimation algorithm in TAC 4.0 software, also used to display differential splicing data for specific genes.
The lists of differentially expressed/spliced genes for each factor considered (e.g. gestational age) were tested for enrichment based on chromosomal location and membership in previously described functional categories and pathways as well as on specificity to tissues and cell types. The functional databases considered were the Developmental FunctionaL Annotation at Tufts (DFLAT) database  and the Curated Gene Sets (C2) collection from the Molecular Signatures Database (MSigDB) database . Tissue-specific genes were defined as those with median expression > 30 times higher in a given tissue than the median expression of all other tissues described in the Gene Atlas . This cut-off was chosen to enable the direct comparison of findings to previous reports . All enrichment analyses were based on a hypergeometric test (equivalent to a Fisher’s exact test) and accounted for multiple testing, with q < 0.05 being considered a significant result.
Analysis of single-cell RNA-Seq signature expression in AF
Log2 microarray expression data were transformed into Z-scores for each gene by subtracting the mean and dividing to the standard deviation calculated from the reference study group (e.g. midtrimester when assessing changes from midtrimester to term). The Z-scores in each sample were averaged over the set of genes previously defined as specific to a given population of cells defined by single-cell RNA-Seq analyses in Tsang et al.  The collection of gene sets included those for extravillous trophoblasts, cytotrophoblasts, the syncytiotrophoblast, decidual cells, dendritic cells, endothelial cells, erythrocytes, Hofbauer cells, stromal cells, vascular smooth muscle cells, B cells, T cells, and monocytes. Unlike averaging over the normalized expression of genes, as in Tsang et al. , the standardization of expression data ensures that genes contribute equally to the gene set summary [108, 121]. The average Z-score for each single-cell signature was compared between the TNL and midtrimester groups by using the Wilcoxon rank sum test as well as within the midtrimester group via linear models by correlating gene expression with gestational age (continuous).
Clinical characteristics of the study population
We profiled the cell-free transcriptome in AF samples collected during gestation at midtrimester (n = 30) and at term from women without labor (n = 68). The median gestational ages were 21 and 39 weeks in the midtrimester and TNL groups, respectively, at the time of sample collection. Women in the midtrimester group were more likely to be nulliparous (26.7%) compared to those in the TNL group (4.4%) (p = 0.003). There were no differences in maternal age, body mass index (BMI), smoking status, fetal sex, and gestational age at delivery between the midtrimester and TNL groups (Table 1).
Factors affecting the AF transcriptome in normal pregnancy
The Human Transcriptome Array 2.0 platform used in this study was designed to probe at exon-level resolution the expression of 44,699 protein coding and 22,829 non-protein coding transcript clusters, simply referred to herein as genes (67,528 total). Of these, 64,071 (95%) were deemed expressed (present) in at least 25% of the AF samples in either the midtrimester or the TNL group, and were retained for further differential expression and splicing analyses. An unsupervised principal component analysis representation of the genome-wide gene expression profiles based on the top 1000 most varying genes across all samples is shown in Fig. 1a. The samples in this figure are clustered by gestational age groups, which is suggestive  of large between-group differences. Moreover, the first principal component (PC1) was not only linearity correlated with gestational age overall (R = 0.96) but also within the subset of midtrimester samples alone (R = 0.72) (both, p < 0.001) (Fig. 1b).
We considered several maternal and fetal characteristics as well as the mode of AF sample collection and indication for amniocentesis to evaluate their effects on overall gene expression and splicing (Table 2). Overall, gene expression and splicing patterns were found to change dramatically with gestational age and modestly with fetal sex. There was also weak evidence of changes in maternal obesity, whereas other maternal characteristics (parity, smoking status, and race) and mode of AF sample collection (cesarean section versus transabdominal amniocentesis) and indication for amniocentesis did not have a significant effect on overall gene expression (Table 2). The effects of gestational age and fetal sex on gene expression and splicing are detailed in the next section.
Effect of gestational age on the AF cfRNA
Advancing gestational age from midtrimester to term was associated with expression changes in about 10% of the genes detected in AF (6194/64,071), representing 11% of the coding and 6% of the non-coding genes detected. Gene expression changes included both an increase in expression from midtrimester to term gestation (2776 genes) and a decrease (3418 genes) (q-value < 0.05 and fold change > 1.25) (Additional file 7: Table S1). In addition to assessing how the AF transcriptome changes from midtrimester to term, we explored for the first time how the transcriptional program changes with advancing gestation during the midtrimester, which was not feasible in previous studies due to sample size limitations. A differential expression analysis between samples collected at early (16.4–21.0 weeks) and late (21.1–24.0 weeks) midtrimester identified 413 differentially expressed genes (Additional file 8: Table S2).
In-silico validation of differential expression
To demonstrate the plausibility of AF differential expression from midtrimester to term, we conducted an in-silico analysis to determine whether previously reported findings support our data and vice versa. In the first analysis, we considered all 2719 genes reported by Hui et al.  to change with gestational age (term versus midtrimester) and they were also detected present in the samples in this study (regardless of statistical significance). We found a substantial agreement in terms of direction of change (92% match) and correlation of log fold changes (Spearman’s correlation, R = 79%, p < 0.0001) (Fig. 2a). When only the subset of 1332 genes that were significant in both studies were considered, the agreement in terms of direction of change reached 99%, while the correlation of log2 fold changes increased to 0.82 (p < 0.001). In a second in-silico experiment, we considered the genes detected in the current study (regardless of significance) that were differentially expressed (q-value < 0.05 fold change > 1.25) based on a re-analysis of the RNA-Seq data reported by Kamath-Rayne et al. . Our expression-change estimates for the TNL to midtrimester comparison were highly correlated with those of the 1234 genes that differed between late preterm and midtrimester (R = 0.83, p < 0.001, 97% direction of change agreement, Fig. 2b) and with those of the 1420 genes that differed between term and midtrimester (R = 0.79, p < 0.001, 97% direction of change agreement, Fig. 2c), obtained by re-analysis of the Kamath-Rayne et al.  dataset. Overall, these results demonstrate high cross-study reproducibility of gene-level differential-expression changes with gestation in amniotic fluid.
Although differential expression with advancing gestation was noted throughout the genome, five chromosomes (Chr1, Chr4, Chr6, Chr12, and Chr15) displayed slightly more differential expression than expected by chance (see Additional file 1: Figure S1) (q < 0.05, OR 1.2–1.4 for all).
Biological processes enrichment
A wide range of biological processes documented in the DFLAT database  were over-represented among the genes that change in expression from midterm to term gestation in the AF. Placental development, organ development (lung, liver, brain, heart, etc.), and immune-related pathways (positive regulation of immune response, activation of immune response, T-cell activation) are a few of more than 1500 biological processes related to gene expression changes (see Additional file 9: Table S3).
Canonical pathways enrichment
Similarly, about 250 pathways curated in the MSigDB collection were associated with gene expression changes with gestational age, such as the KEGG database systemic lupus erythematosus pathways and the Reactome database amyloids, cell cycle, developmental biology, immune system, cytokine signaling in immune system, and mRNA splicing pathways (q < 0.05) (Additional file 10: Table S4).
Tissue enrichment and signature analysis
Finally, we tested the association of mRNA modulation with advancing gestational age and defined sets based on the Gene Atlas  and found that most over-represented organs among genes with higher expression at term were the trachea, lung, salivary glands, tonsils, tongue, colon, bone marrow, skin, and fetal lung, among others listed in Additional file 11: Table S5 (q < 0.05). Although fetal skin was not represented in the Gene Atlas-based analysis, it was proposed by Hui et al.  that skin-specific transcripts identified in amniotic fluid are likely derived from the fetal skin. The most-enriched organs for genes with higher expression during midterm gestation were the small intestine, placenta, uterus, and specific cell types (e.g. CD105+ endothelial cells, cardiac myocytes), among others listed in Additional file 12: Table S6 (q < 0.05). These findings were also supported by an alternative analysis in which the expression signature of each tissue type (based on the average of the top 20 most-specific genes) was analyzed as a continuous response as a function of gestational age (Fig. 3a&b and Additional file 2: Figure S2). Of note, the expression of gene signatures for the trachea, salivary glands, and lungs increased while those for the cardiac myocytes and uterus decreased steadily throughout gestation; yet, more complex patterns emerged for the pituitary gland and fetal liver, whose expression signatures peaked and bottomed toward the end of the midtrimester, respectively (Fig. 3b).
Single-cell RNA-Seq signature analysis
The placenta emerged as one of the organs associated with both increase (Odds Ratio, OR = 2.1) (Additional file 11: Table S5) and decrease (OR = 3.1) (Additional file 12: Table S6) in gene expression from midtrimester to term gestation. We sought to further dissect placenta-specific gene expression by averaging the expression of genes specific to sub-populations of cells based on single-cell genomics studies. Among the 13 cell types identified by an RNA-Seq analysis of the placenta by Tsang et al. , the cytotrophoblast, monocyte, and syncytiotrophoblast expression increased monotonically from 16 weeks of gestation until term. By contrast, the Hofbauer cells, erythrocytes, vascular smooth muscle cells, B cells, T cells, and others showed more complex patterns (increase followed by decrease) (q < 0.05) (Fig. 4, Additional file 3: Figure S3). Of these single-cell signatures, that of the cytotrophoblast, defined as the average expression of the FAM3B, FOXO4, and MIR205HG genes, was the most highly modulated, being increased at term 1.8 SD relative to the midtrimester group, mostly due to the contribution of FAM3B gene expression.
Alternative splicing (AS) events associated with advancement from midtrimester to term gestation were identified in 17.5% (8566/48,820) of the genes detected and met the conditions for testing (see Methods) (exon q-value < 0.05 and SI > 2.0 or SI < − 2.0) (Additional file 13: Table S7). Of all tested genes, 25% of the coding genes and 4% of the non-coding genes displayed AS events associated with gestational-age difference from midtrimester to term gestation. Of note, 85% of all differential splicing results according to the SI method in TAC 4.0 software were also supported by the diffSplice algorithm (q < 0.05). Of approximately 54% of differentially spliced genes for which a particular type of AS event was identified by the Event Estimation algorithm, the most common types were the cassette exon (69%), alternative 5′ site (19%), alternative 3′ site (10%), and intron retention (2%) (Additional file 13: Table S7). For genes differentially expressed with gestational age, the AS analysis pinpointed the transcript isoforms likely responsible for overall gene expression changes: this is illustrated in Fig. 5 for the MUC7 (salivary gland-specific), SFTPD (lung-specific), and GKN1 (stomach-specific) genes. These genes were previously reported to be among those most differentially expressed with gestational age in AF based on 3′-end biased microarray platforms studies . For example, our data suggest that most of the increase in expression of MUC7 from midtrimester to term gestation can be explained by a short isoform of this gene (see transcript TR04000256, Fig. 5a).
The majority (76%) of differentially spliced genes were not differentially expressed, showcasing the importance of differential splicing in the study of the AF transcriptome. For instance, CNIH1 skipped and ZNF365 included a particular exon more frequently in the TNL group compared to the midtrimester group, yet neither gene met the criteria of being differentially expressed (Fig. 6). Functional profiling of differentially spliced genes identified about 800 DFLAT database biological processes as enriched that were not identified as such by analysis of differential expression, including 47 processes related to development (ear, central nervous system neurons, tongue, and spleen) as well as several immune-related processes (Additional file 14: Table S8). Among the MSigDB database pathways found to be enriched based on differential splicing, but not differential expression analysis, the REACTOME adaptive immune system, the BIOCARTA MEF2D pathway, and the KEGG MAPK pathway were among the most enriched (Additional file 15: Table S9). Tissues and cell types associated with gestational age by differential splicing analysis, but missed by differential expression analysis, were brain tissues (e.g. prefrontal cortex, globus pallidus, and cerebellum peduncles) and T cells (CD8+ and CD4+) (Table 3).
Effect of fetal sex on the AF transcriptome
AF mRNA expression changes with fetal sex were found for 252 genes, with 215 being decreased and 37 being increased in expression in pregnancies with a male compared to those with a female fetus (Additional file 16: Table S10). All genes with increased expression in AF samples of women with a male fetus compared to those with a female fetus, were located on the chromosome Y (q < 0.05, OR = 64, Additional file 4: Figure S4), with ribosomal protein S4, Y-linked 1 (RPS4Y1) being the most increased in pregnancies with a male fetus (8-fold increase) (Additional file 16: Table S10). This gene was also reported by Larrabee et al.  as being present in the AF of women with a male fetus but not in those with a female fetus. Functional profiling analysis revealed that fetal sex gene expression differences were associated with one PID database pathway (P 53 downstream pathway) as well as in particular tissues (tongue and tonsils) and cell types (bronchial epithelial cells, CD71+ early erythroid cells) (all q < 0.05). When differential splicing associated with fetal sex was assessed, we found 240 significant genes, including some deemed also differentially expressed (e.g. RPS4Y1) and others that were not differentially expressed, including a transcript cluster for TXLNGY (Additional file 5: Figure S5). Chromosomes Y (odds ratio [OR] = 70) and X (OR = 2.4) were significantly enriched in genes with differential splicing between pregnancies with a male or a female fetus (q < 0.05).
Effect of maternal characteristics and mode of sample collection
Motivated by the findings of Edlow et al.  regarding an association between the AF transcriptome and maternal obesity, we also compared gene expression between 38 obese (BMI > 30) and 33 lean (BMI < 25) women, while adjusting for fetal sex and gestational age at amniocentesis. No significant differences in expression and splicing were found with maternal obesity for individual genes. The same was true when limiting the analysis only to women sampled during the midtrimester, as in the original report by Edlow et al.  (Table 2). However, we found a weak but significant correlation of fold changes derived herein for the 182 genes reported to change with obesity by Edlow et al.  and present on the microarray platform herein (Spearman’s correlation 0.2, p = 0.02), with 62% of those genes matching in terms of the direction of change with obesity (Additional file 6: Figure S6).
Finally, we also examined the effect of other maternal characteristics, including ethnicity (African American versus Caucasian), parity (parous versus nulliparous), smoking status, as well as the mode of collection (cesarean delivery versus transabdominal amniocentesis) and indication for amniocentesis, while adjusting for variables that had a significant effect (gestational-age group and fetal sex). No significant differences in overall gene expression were found based on these analyses (all q > 0.1) (Table 2).
Amniotic fluid cfRNA analysis was proposed as a means to study real-time human fetal physiology and development [83, 85, 88, 99]. The results presented herein are in agreement with previous reports regarding differential gene expression from midtrimester to term gestation; yet, we have shown for the first time strong expression changes even during the midtrimester alone (before 21 weeks of gestation versus after). Tissue-specific mRNA expression patterns with gestation are found herein to be more complex than previously reported, owing to better coverage of the second-trimester gestational-age range. Leveraging previous single-cell genomics studies of the placenta, we also show for the first time that it is feasible to track signatures of placental single-cell populations by AF cfRNA analysis to assess the dynamic of crosstalk at the maternal-fetal interface. Finally, we present evidence that alternative splicing, a mechanism described to contribute to tissue-identity acquisition and organ development , is modulated in the AF with advancing gestational age and associated with the development of complex organ systems (e.g. brain).
Amniotic fluid differential expression with advancing gestation
Hui et al.  and Kamath-Rayne et al.  reported that AF cfRNA displays dramatic changes with gestational age using 3′-end biased microarrays and RNA-Seq, respectively. Using a more recent microarray platform that probes both the coding and non-coding genes at exon-level resolution, we found 6194 differentially expressed transcript clusters (3447 unique ENTREZ database annotated genes) to be differentially expressed between midtrimester and term samples (Additional file 7: Table S1). The agreement between results presented herein and these two previous studies was high in terms of direction of change and correlation of fold changes (Fig. 2), yet the number of genes identified herein was larger owing to an increased sample size and sensitivity of the microarray platform employed. Enrichment analyses based on differentially expressed genes identified similar tissues and organ systems related to the in utero environment, fetal development, and preparation for life outside the uterus [83, 98]. However, given better coverage of the gestational-age span during midtrimester gestation, we show that the expression of tissue-specific signatures can have a complex pattern of modulation, including the fetal liver and cardiac myocytes, among others (Fig. 3).
Single-cell RNA-Seq signature modulation with advancing gestation
Single-cell RNA-seq signatures of populations of placental cells were recently described as a vocabulary to interpret the molecular crosstalk at the maternal-fetal interface [105, 106, 124]. The ability to track the expression of these signatures throughout gestation  and their alterations with obstetrical disease [105, 124] were also reported. The current study provides new and original evidence that fetal-specific (syncytiotrophoblast, cytotrophoblasts, Hofbauer cells, and vascular smooth muscle cells) or maternal-fetal origin cell populations (erythrocytes, monocytes, B cells and T cells) change with advancing gestation in AF. The increase in the expression of the monocyte and syncytiotrophoblast signatures reported in Fig. 4 is in agreement with observations based on cell-free  transcriptome analysis in maternal circulation. Nonetheless, the sharp increase during midtrimester gestation, followed by a decrease at near-term gestation, for Hofbauer and other cell signatures shown in Fig. 3 are for the first time described herein. Importantly, the AF single-cell RNA signatures of monocytes, T cells, and B cells mirrored the abundance of these immune cells throughout normal pregnancy as quantified by flow cytometry .
Amniotic fluid differential splicing with advancing gestation
Through analysis of gene expression at the level of individual exons and exon-exon junctions, we identified AS events associated with gestational age and fetal sex. The exon junction arrays used in this study were previously compared to RNA-Seq for the purpose of differential splicing and found to have higher power when quantifying low-abundance transcripts as well as long non-coding RNAs that tend to be shorter than protein-coding gene counterparts . Of interest, we found more genes (17.5% of the genes detected) displaying differential splicing than differential expression (10%) with gestational age. This can be explained, in part, by the fact that there is a higher threshold to claim differential expression than differential splicing. For differential expression, changes need to concur for multiple exonic regions (Fig. 5), while for differential splicing, changes for even a single or a few exonic region are sufficient (Fig. 6).
Genes showing differential splicing were associated with processes and pathways related to development and immunity and were specific to certain complex organ systems that were not identified as enriched based on differential expression analysis alone (Table 3). AS is recognized as a fundamental process by which cells expand their transcriptomic diversity, and it is particularly widespread in the nervous system . Therefore, consideration of AS, as a means to assess the maturity or developmental stage of fetal organs, is important.
Strengths and limitations
This is the largest study of the AF transcriptome in which coding and non-coding gene expression was profiled at exon-level resolution in approximately 100 normal pregnancies. One of the strengths of the study is the simultaneous assessment of the effect of genetic (fetal sex and maternal ethnicity), physiologic (gestational age, parity), and maternal risk factors (smoking, obesity) on the transcriptome. Also, the observation that there are no significant effects related to the AF sample collection mode (cesarean delivery versus transabdominal amniocentesis) is an important addition to the literature. This suggests that samples collected by both modalities can be used to establish a gene expression reference. Although the sample size for this particular comparison (cesarean delivery versus transabdominal amniocentesis at term) was low, the reduced magnitude of such an effect was reassuring. Of note, the comparison between the two types of sample collection at term gestation was also cofounded by the indication for amniocentesis, since transabdominal collection was performed to assess fetal lung maturity while collection during cesarean delivery was done for research purposes. The use of multiple types of functional profiling approaches and single-cell signatures as a means to interpret differential expression and splicing results is also a strength. Finally, although no additional wet-lab confirmatory results were presented, the in silico validation of overall gene expression differences with gestational age and fetal sex, through correlations with previous reports, increases confidence in the novel results presented herein, such as differential splicing with advancing gestation and fetal gender. Among the limitations, we would also note a lower statistical power for some of the sub-analyses assessing the effects of maternal race and obesity in midtrimester samples.
We reported herein the largest AF cell-free transcriptomics study that catalogues physiologic adaptations with advancing gestation in normal pregnancy and surveys the effects of relevant maternal, fetal, and experimental covariates on the transcriptome. Our data show that AF mRNA profiles can be used to track placental function through single-cell specific signatures, as a readout of the maternal-fetal crosstalk during pregnancy. We also propose that alternative splicing evaluation should be a part of the future development of fetal testing to assess organ maturity; this information could be used to inform clinical management given the current debate about the usefulness of fetal lung maturity evaluation.
Availability of data and materials
The raw and processed gene expression data and relevant patient characteristics were deposited in the Gene Expression Omnibus (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE133824).
Davis LE, McLaren LC, Stewart JA, James CG, Levine MD, Skipper BJ. Immunological and microbiological studies of midtrimester amniotic fluid. Gynecol Obstet Investig. 1983;16(5):261–8.
Schmidt W. The amniotic fluid compartment: the fetal habitat. Adv Anat Embryol Cell Biol. 1992;127:1–100.
Underwood MA, Gilbert WM, Sherman MP. Amniotic fluid: not just fetal urine anymore. J Perinatol: official journal of the California Perinatal Association. 2005;25(5):341–8.
Cherry SH. Amniotic fluid analysis as an index of fetal health in utero. Med Times. 1967;95(7):713–7.
Barham KA. Amnioscopy, meconium and fetal well-being. J Obstet Gynaecol Br Commonw. 1969;76(5):412–8.
Schifrin BS, Guntes V, Gergely RC, Eden R, Roll K, Jacobs J. The role of real-time scanning in antenatal fetal surveillance. Am J Obstet Gynecol. 1981;140(5):525–30.
Jacobson CB, Barter RH. Intrauterine diagnosis and management of genetic defects. Am J Obstet Gynecol. 1967;99(6):796–807.
Valenti C, Schutta EJ, Kehaty T. Cytogenetic diagnosis of Down's syndrome in utero. JAMA : J Am Med Assoc. 1969;207(8):1513–5.
Santesson B, Akesson HO, Book JA, Brosset A. Karyotyping human amnitic-fluid cels. Lancet. 1969;2(7629):1067–8.
Lisgar F, Gertner M, Cherry S, Hsu LY, Hirschhorn K. Prenatal chromosome analysis. Nature. 1970;225(5229):280–1.
Gertner M, Hsu LY, Martin J, Hirshhorn K. The use of amniocentesis for prenatal genetic counseling. Bull N Y Acad Med. 1970;46(11):916–21.
Hobbins JC, Brock W, Speroff L, Anderson GG, Caldwell B. L-S ratio in predicting pulmonary maturity in utero. Obstet Gynecol. 1972;39(5):660–4.
Gluck L. The evaluation of fetal lung maturity. Calif Med. 1972;116(4):58–9.
Winn HN, Romero R, Roberts A, Liu H, Hobbins JC. Comparison of fetal lung maturation in preterm singleton and twin pregnancies. Am J Perinatol. 1992;9(5–6):326–8.
Palacio M, Bonet-Carne E, Cobo T, Perez-Moreno A, Sabria J, Richter J, et al. Prediction of neonatal respiratory morbidity by quantitative ultrasound lung texture analysis: a multicenter study. Am J Obstet Gynecol. 2017;217(2):196 e1–e14.
Romero R, Emamian M, Quintero R, Wan M, Hobbins JC, Mitchell MD. Amniotic fluid prostaglandin levels and intra-amniotic infections. Lancet. 1986;1(8494):1380.
Romero R, Emamian M, Quintero R, Wan M, Hobbins JC, Mazor M, et al. The value and limitations of the gram stain examination in the diagnosis of intraamniotic infection. Am J Obstet Gynecol. 1988;159(1):114–9.
Romero R, Emamian M, Quintero R, Wan M, Scioscia AL, Hobbins JC, et al. Diagnosis of intra-amniotic infection: the acridine orange stain. Am J Perinatol. 1989;6(1):41–5.
Romero R, Manogue KR, Mitchell MD, Wu YK, Oyarzun E, Hobbins JC, et al. Infection and labor. IV. Cachectin-tumor necrosis factor in the amniotic fluid of women with intraamniotic infection and preterm labor. Am J Obstet Gynecol. 1989;161(2):336–41.
Romero R, Sirtori M, Oyarzun E, Avila C, Mazor M, Callahan R, et al. Infection and labor. V. Prevalence, microbiology, and clinical significance of intraamniotic infection in women with preterm labor and intact membranes. Am J Obstet Gynecol. 1989;161(3):817–24.
Romero R, Avila C, Santhanam U, Sehgal PB. Amniotic fluid interleukin 6 in preterm labor. Association with infection. J Clin Invest. 1990;85(5):1392–400.
Romero R, Jimenez C, Lohda AK, Nores J, Hanaoka S, Avila C, et al. Amniotic fluid glucose concentration: a rapid and simple method for the detection of intraamniotic infection in preterm labor. Am J Obstet Gynecol. 1990;163(3):968–74.
Romero R, Quintero R, Nores J, Avila C, Mazor M, Hanaoka S, et al. Amniotic fluid white blood cell count: a rapid and simple test to diagnose microbial invasion of the amniotic cavity and predict preterm delivery. Am J Obstet Gynecol. 1991;165(4 Pt 1):821–30.
Baumann P, Romero R, Berry S, Gomez R, McFarlin B, Araneda H, et al. Evidence of participation of the soluble tumor necrosis factor receptor I in the host response to intrauterine infection in preterm labor. Am J Reprod Immunol. 1993;30(2–3):184–93.
Cherouny PH, Pankuch GA, Romero R, Botti JJ, Kuhn DC, Demers LM, et al. Neutrophil attractant/activating peptide-1/interleukin-8: association with histologic chorioamnionitis, preterm delivery, and bioactive amniotic fluid leukoattractants. Am J Obstet Gynecol. 1993;169(5):1299–303.
Romero R, Yoon BH, Kenney JS, Gomez R, Allison AC, Sehgal PB. Amniotic fluid interleukin-6 determinations are of diagnostic and prognostic value in preterm labor. Am J Reprod Immunol. 1993;30(2–3):167–83.
Gomez R, Romero R, Galasso M, Behnke E, Insunza A, Cotton DB. The value of amniotic fluid interleukin-6, white blood cell count, and gram stain in the diagnosis of microbial invasion of the amniotic cavity in patients at term. Am J Reprod Immunol. 1994;32(3):200–10.
Hsu CD, Meaddough E, Aversa K, Copel JA. The role of amniotic fluid L-selectin, GRO-alpha, and interleukin-8 in the pathogenesis of intraamniotic infection. Am J Obstet Gynecol. 1998;178(3):428–32.
Hsu CD, Meaddough E, Aversa K, Hong SF, Lu LC, Jones DC, et al. Elevated amniotic fluid levels of leukemia inhibitory factor, interleukin 6, and interleukin 8 in intra-amniotic infection. Am J Obstet Gynecol. 1998;179(5):1267–70.
Maymon E, Romero R, Pacora P, Gervasi MT, Edwin SS, Gomez R, et al. Matrilysin (matrix metalloproteinase 7) in parturition, premature rupture of membranes, and intrauterine infection. Am J Obstet Gynecol. 2000;182(6):1545–53.
Maymon E, Romero R, Pacora P, Gervasi MT, Gomez R, Edwin SS, et al. Evidence of in vivo differential bioavailability of the active forms of matrix metalloproteinases 9 and 2 in parturition, spontaneous rupture of membranes, and intra-amniotic infection. Am J Obstet Gynecol. 2000;183(4):887–94.
Maymon E, Romero R, Pacora P, Gomez R, Athayde N, Edwin S, et al. Human neutrophil collagenase (matrix metalloproteinase 8) in parturition, premature rupture of the membranes, and intrauterine infection. Am J Obstet Gynecol. 2000;183(1):94–9.
Pacora P, Maymon E, Gervasi MT, Gomez R, Edwin SS, Yoon BH, et al. Lactoferrin in intrauterine infection, human parturition, and rupture of fetal membranes. Am J Obstet Gynecol. 2000;183(4):904–10.
Yoon BH, Romero R, Kim M, Kim EC, Kim T, Park JS, et al. Clinical implications of detection of Ureaplasma urealyticum in the amniotic cavity with the polymerase chain reaction. Am J Obstet Gynecol. 2000;183(5):1130–7.
Hsu CD, Hong SF, Harirah H, Bahado-Singh R, Lu L. Amniotic fluid soluble fas levels in intra-amniotic infection. Obstet Gynecol. 2000;95(5):667–70.
Hsu CD, Aversa K, Meaddough E. The role of amniotic fluid interleukin-6, and cell adhesion molecules, intercellular adhesion molecule-1 and leukocyte adhesion molecule-1, in intra-amniotic infection. Am J Reprod Immunol. 2000;43(5):251–4.
Maymon E, Romero R, Chaiworapongsa T, Berman S, Conoscenti G, Gomez R, et al. Amniotic fluid matrix metalloproteinase-8 in preterm labor with intact membranes. Am J Obstet Gynecol. 2001;185(5):1149–55.
Yoon BH, Oh SY, Romero R, Shim SS, Han SY, Park JS, et al. An elevated amniotic fluid matrix metalloproteinase-8 level at the time of mid-trimester genetic amniocentesis is a risk factor for spontaneous preterm delivery. Am J Obstet Gynecol. 2001;185(5):1162–7.
Yoon BH, Romero R, Moon JB, Shim SS, Kim M, Kim G, et al. Clinical significance of intra-amniotic inflammation in patients with preterm labor and intact membranes. Am J Obstet Gynecol. 2001;185(5):1130–6.
Keelan JA, Yang J, Romero RJ, Chaiworapongsa T, Marvin KW, Sato TA, et al. Epithelial cell-derived neutrophil-activating peptide-78 is present in fetal membranes and amniotic fluid at increased concentrations with intra-amniotic infection and preterm delivery. Biol Reprod. 2004;70(1):253–9.
Espinoza J, Goncalves LF, Romero R, Nien JK, Stites S, Kim YM, et al. The prevalence and clinical significance of amniotic fluid 'sludge' in patients with preterm labor and intact membranes. Ultrasound Obstet Gynecol. 2005;25(4):346–52.
Kusanovic JP, Espinoza J, Romero R, Goncalves LF, Nien JK, Soto E, et al. Clinical significance of the presence of amniotic fluid 'sludge' in asymptomatic patients at high risk for spontaneous preterm delivery. Ultrasound Obstet Gynecol. 2007;30(5):706–14.
Lee SE, Romero R, Jung H, Park CW, Park JS, Yoon BH. The intensity of the fetal inflammatory response in intraamniotic inflammation with and without microbial invasion of the amniotic cavity. Am J Obstet Gynecol. 2007;197(3):294 e1–6.
Romero R, Espinoza J, Hassan S, Gotsch F, Kusanovic JP, Avila C, et al. Soluble receptor for advanced glycation end products (sRAGE) and endogenous secretory RAGE (esRAGE) in amniotic fluid: modulation by infection and inflammation. J Perinat Med. 2008;36(5):388–98.
Romero R, Espinoza J, Rogers WT, Moser A, Nien JK, Kusanovic JP, et al. Proteomic analysis of amniotic fluid to identify women with preterm labor and intra-amniotic inflammation/infection: the use of a novel computational method to analyze mass spectrometric profiling. J Maternal-Fetal Neonatal Med : Official J Eur Assoc Perinat Med, Fed Asia Oceania Perinat Soc, Int Soc Perinat Obstet. 2008;21(6):367–88.
Romero R, Schaudinn C, Kusanovic JP, Gorur A, Gotsch F, Webster P, et al. Detection of a microbial biofilm in intraamniotic infection. American journal of obstetrics and gynecology. 2008;198(1):135 e1–5.
Romero R, Chaiworapongsa T, Alpay Savasan Z, Xu Y, Hussein Y, Dong Z, et al. Damage-associated molecular patterns (DAMPs) in preterm labor with intact membranes and preterm PROM: a study of the alarmin HMGB1. J Matern Fetal Neonatal Med. 2011;24(12):1444–55.
Gervasi MT, Romero R, Bracalente G, Erez O, Dong Z, Hassan SS, et al. Midtrimester amniotic fluid concentrations of interleukin-6 and interferon-gamma-inducible protein-10: evidence for heterogeneity of intra-amniotic inflammation and associations with spontaneous early (<32 weeks) and late (>32 weeks) preterm delivery. J Perinat Med. 2012;40(4):329–43.
Gervasi MT, Romero R, Bracalente G, Chaiworapongsa T, Erez O, Dong Z, et al. Viral invasion of the amniotic cavity (VIAC) in the midtrimester of pregnancy. J Matern Fetal Neonatal Med. 2012;25(10):2002–13.
Romero R, Miranda J, Chaiworapongsa T, Chaemsaithong P, Gotsch F, Dong Z, et al. A novel molecular microbiologic technique for the rapid diagnosis of microbial invasion of the amniotic cavity and intra-amniotic infection in preterm labor with intact membranes. Am J Reprod Immunol. 2014;71(4):330–58.
Romero R, Miranda J, Chaiworapongsa T, Chaemsaithong P, Gotsch F, Dong Z, et al. Sterile intra-amniotic inflammation in asymptomatic patients with a sonographic short cervix: prevalence and clinical significance. J Matern Fetal Neonatal Med. 2014:1–17.
Romero R, Miranda J, Chaiworapongsa T, Korzeniewski SJ, Chaemsaithong P, Gotsch F, et al. Prevalence and clinical significance of sterile intra-amniotic inflammation in patients with preterm labor and intact membranes. Am J Reprod Immunol. 2014;72(5):458–74.
Romero R, Miranda J, Chaemsaithong P, Chaiworapongsa T, Kusanovic JP, Dong Z, et al. Sterile and microbial-associated intra-amniotic inflammation in preterm prelabor rupture of membranes. J Matern Fetal Neonatal Med. 2015;28(12):1394–409.
Romero R, Miranda J, Kusanovic JP, Chaiworapongsa T, Chaemsaithong P, Martinez A, et al. Clinical chorioamnionitis at term I: microbiology of the amniotic cavity using cultivation and molecular techniques. J Perinat Med. 2015;43(1):19–36.
Romero R, Grivel JC, Tarca AL, Chaemsaithong P, Xu Z, Fitzgerald W, et al. Evidence of perturbations of the cytokine network in preterm labor. Am J obstet Gynecol. 2015;213(6):836.e1–e18.
Chaemsaithong P, Romero R, Korzeniewski SJ, Martinez-Varea A, Dong Z, Yoon BH, et al. A point of care test for interleukin-6 in amniotic fluid in preterm prelabor rupture of membranes: a step toward the early treatment of acute intra-amniotic inflammation/infection. J Matern Fetal Neonatal Med. 2016;29(3):360–7.
Chaemsaithong P, Romero R, Korzeniewski SJ, Martinez-Varea A, Dong Z, Yoon BH, et al. A rapid interleukin-6 bedside test for the identification of intra-amniotic inflammation in preterm labor with intact membranes. J Matern Fetal Neonatal Med. 2016;29(3):349–59.
Romero R, Chaemsaithong P, Korzeniewski SJ, Kusanovic JP, Docheva N, Martinez-Varea A, et al. Clinical chorioamnionitis at term III: how well do clinical criteria perform in the identification of proven intra-amniotic infection? J Perinat Med. 2016;44(1):23–32.
Romero R, Chaemsaithong P, Korzeniewski SJ, Tarca AL, Bhatti G, Xu Z, et al. Clinical chorioamnionitis at term II: the intra-amniotic inflammatory response. J Perinat Med. 2016;44(1):5–22.
Yoneda N, Yoneda S, Niimi H, Ueno T, Hayashi S, Ito M, et al. Polymicrobial amniotic fluid infection with mycoplasma/Ureaplasma and other Bacteria induces severe intra-amniotic inflammation associated with poor perinatal prognosis in preterm labor. Am J Reprod Immunol. 2016;75(2):112–25.
Kim SM, Romero R, Lee J, Chaemsaithong P, Docheva N, Yoon BH. Gastric fluid versus amniotic fluid analysis for the identification of intra-amniotic infection due to Ureaplasma species. J Matern Fetal Neonatal Med. 2016;29(16):2579–87.
Son GH, You YA, Kwon EJ, Lee KY, Kim YJ. Comparative analysis of Midtrimester amniotic fluid cytokine levels to predict spontaneous very pre-term birth in patients with cervical insufficiency. Am J Reprod Immunol. 2016;75(2):155–61.
Maddipati KR, Romero R, Chaiworapongsa T, Chaemsaithong P, Zhou SL, Xu Z, et al. Lipidomic analysis of patients with microbial invasion of the amniotic cavity reveals up-regulation of leukotriene B4. FASEB J : Official Publication Fed Am Soc Exp Biol. 2016;30(10):3296–307.
Yoneda S, Shiozaki A, Yoneda N, Ito M, Shima T, Fukuda K, et al. Antibiotic therapy increases the risk of preterm birth in preterm labor without intra-amniotic microbes, but may prolong the gestation period in preterm labor with microbes, evaluated by rapid and high-sensitive PCR system. Am J Reprod Immunol. 2016;75(4):440–50.
Park JY, Romero R, Lee J, Chaemsaithong P, Chaiyasit N, Yoon BH. An elevated amniotic fluid prostaglandin F2alpha concentration is associated with intra-amniotic inflammation/infection, and clinical and histologic chorioamnionitis, as well as impending preterm delivery in patients with preterm labor and intact membranes. J Matern Fetal Neonatal Med. 2016;29(16):2563–72.
Romero R, Chaemsaithong P, Chaiyasit N, Docheva N, Dong Z, Kim CJ, et al. CXCL10 and IL-6: Markers of two different forms of intra-amniotic inflammation in preterm labor. Am J Reprod Immunol. 2017;78(1).
Chaiyasit N, Romero R, Chaemsaithong P, Docheva N, Bhatti G, Kusanovic JP, et al. Clinical chorioamnionitis at term VIII: a rapid MMP-8 test for the identification of intra-amniotic inflammation. J Perinat Med. 2017.
Gomez-Lopez N, Romero R, Xu Y, Leng Y, Garcia-Flores V, Miller D, et al. Are amniotic fluid neutrophils in women with intraamniotic infection and/or inflammation of fetal or maternal origin? Am J Obstet Gynecol. 2017;217(6):693 e1–e16.
Gomez-Lopez N, Romero R, Xu Y, Miller D, Unkel R, Shaman M, et al. Neutrophil extracellular traps in the amniotic cavity of women with intra-amniotic infection: a new mechanism of host defense. Reprod Sci. 2017;24(8):1139–53.
Martinez-Varea A, Romero R, Xu Y, Miller D, Ahmed AI, Chaemsaithong P, et al. Clinical chorioamnionitis at term VII: the amniotic fluid cellular immune response. J Perinat Med. 2017;45(5):523–38.
Maymon E, Romero R, Bhatti G, Chaemsaithong P, Gomez-Lopez N, Panaitescu B, et al. Chronic inflammatory lesions of the placenta are associated with an up-regulation of amniotic fluid CXCR3: a marker of allograft rejection. J Perinat Med. 2017.
Musilova I, Andrys C, Krejsek J, Drahosova M, Zednikova B, Pliskova L, et al. Amniotic fluid pentraxins: potential early markers for identifying intra-amniotic inflammatory complications in preterm pre-labor rupture of membranes. Am J Reprod Immunol. 2017.
Oh KJ, Kim SM, Hong JS, Maymon E, Erez O, Panaitescu B, et al. Twenty-four percent of patients with clinical chorioamnionitis in preterm gestations have no evidence of either culture-proven intraamniotic infection or intraamniotic inflammation. Am J Obstet Gynecol. 2017;216(6):604 e1–e11.
Tarca AL, Fitzgerald W, Chaemsaithong P, Xu Z, Hassan SS, Grivel JC, et al. The cytokine network in women with an asymptomatic short cervix and the risk of preterm delivery. Am J Reprod Immunol. 2017;78(3).
Rowlands S, Danielewski JA, Tabrizi SN, Walker SP, Garland SM. Microbial invasion of the amniotic cavity in midtrimester pregnancies using molecular microbiology. Am J Obstet Gynecol. 2017;217(1):71 e1–5.
Gomez-Lopez N, Romero R, Garcia-Flores V, Xu Y, Leng Y, Alhousseini A, et al. Amniotic fluid neutrophils can phagocytize bacteria: A mechanism for microbial killing in the amniotic cavity. Am J Reprod Immunol. 2017;78(4).
Gomez-Lopez N, Romero R, Panaitescu B, Leng Y, Xu Y, Tarca AL, Faro J, Pacora P, Hassan SS, Hsu CD. Inflammasome activation during spontaneous preterm labor with intra-amniotic infection or sterile intra-amniotic inflammation. Am J Reprod Immunol. 2018;80(5):e13049. https://doi.org/10.1111/aji.13049. Epub 2018 Sep 18.
Gomez-Lopez N, Romero R, Xu Y, Miller D, Leng Y, Panaitescu B, et al. The immunophenotype of amniotic fluid leukocytes in normal and complicated pregnancies. Am J Reprod Immunol. 2018;79(4):e12827.
Gomez-Lopez N, Romero R, Galaz J, Xu Y, Panaitescu B, Slutsky R, Motomura K, Gill N, Para R, Pacora P, Jung E, Hsu CD. Cellular immune responses in amniotic fluid of women with preterm labor and intra-amniotic infection or intra-amniotic inflammation. Am J Reprod Immunol. 2019;82(5):e13171. https://doi.org/10.1111/aji.13171. Epub 2019 Sep 3.
Romero R, Brody DT, Oyarzun E, Mazor M, Wu YK, Hobbins JC, et al. Infection and labor. III. Interleukin-1: a signal for the onset of parturition. Am J Obstet Gynecol. 1989;160(5 Pt 1):1117–23.
Espinoza J, Chaiworapongsa T, Romero R, Edwin S, Rathnasabapathy C, Gomez R, et al. Antimicrobial peptides in amniotic fluid: defensins, calprotectin and bacterial/permeability-increasing protein in patients with microbial invasion of the amniotic cavity, intra-amniotic inflammation, preterm labor and premature rupture of membranes. J Matern Fetal Neonatal Med. 2003;13(1):2–21.
Gomez-Lopez N, Romero R, Leng Y, Xu Y, Slutsky R, Levenson D, Pacora P, Jung E, Panaitescu B, Hsu CD. The origin of amniotic fluid monocytes/macrophages in women with intra-amniotic inflammation or infection. J Perinat Med. 2019;47(8):822–840. https://doi.org/10.1515/jpm-2019-0262.
Larrabee PB, Johnson KL, Lai C, Ordovas J, Cowan JM, Tantravahi U, et al. Global gene expression analysis of the living human fetus using cell-free messenger RNA in amniotic fluid. JAMA : the journal of the American Medical Association. 2005;293(7):836–42.
Cho CK, Shan SJ, Winsor EJ, Diamandis EP. Proteomics analysis of human amniotic fluid. Molecular & cellular proteomics : MCP. 2007;6(8):1406–15.
Kamath-Rayne BD, Du Y, Hughes M, Wagner EA, Muglia LJ, DeFranco EA, et al. Systems biology evaluation of cell-free amniotic fluid transcriptome of term and preterm infants to detect fetal maturity. BMC Med Genet. 2015;8:67.
Kang JH, Park HJ, Jung YW, Shim SH, Sung SR, Park JE, et al. Comparative Transcriptome analysis of cell-free fetal RNA from amniotic fluid and RNA from Amniocytes in uncomplicated pregnancies. PLoS One. 2015;10(7):e0132955.
Tarui T, Kim A, Flake A, McClain L, Stratigis JD, Fried I, et al. Amniotic fluid transcriptomics reflects novel disease mechanisms in fetuses with myelomeningocele. Am J Obstet Gynecol. 2017;217(5):587.e1–e10.
Kamath-Rayne BD, Smith HC, Muglia LJ, Morrow AL. Amniotic fluid: the use of high-dimensional biology to understand fetal well-being. Reprod Sci. 2014;21(1):6–19.
Bujold E, Romero R, Kusanovic JP, Erez O, Gotsch F, Chaiworapongsa T, et al. Proteomic profiling of amniotic fluid in preterm labor using two-dimensional liquid separation and mass spectrometry. J Matern Fetal Neonatal Med. 2008;21(10):697–713.
Govia RNM, Birse KD, Sepehri S, Khafipour E, Menticoglou SM, Burgener AD, et al. Amniotic fluid proteomic signatures of cervical insufficiency and their association with length of latency. Am J Reprod Immunol. 2018;80(5):e13030.
Romero R, Mazaki-Tovi S, Vaisbuch E, Kusanovic JP, Chaiworapongsa T, Gomez R, et al. Metabolomics in premature labor: a novel approach to identify patients at risk for preterm delivery. J Matern Fetal Neonatal Med. 2010;23(12):1344–59.
Baraldi E, Giordano G, Stocchero M, Moschino L, Zaramella P, Tran MR, et al. Untargeted Metabolomic analysis of amniotic fluid in the prediction of preterm delivery and Bronchopulmonary dysplasia. PLoS One. 2016;11(10):e0164211.
Wan J, Jiang F, Zhang J, Xu Q, Chen D, Yu B, et al. Amniotic fluid metabolomics and biochemistry analysis provides novel insights into the diet-regulated foetal growth in a pig model. Sci Rep. 2017;7:44782.
Bahado-Singh RO, Sonek J, McKenna D, Cool D, Aydas B, Turkoglu O, et al. Artificial intelligence and amniotic fluid multiomics analysis: the prediction of perinatal outcome in asymptomatic short cervix. Ultrasound Obstet Gynecol. 2018.
Iuculano A, Murgia F. Peddes C. Tronci L, Deiana M, et al. Metabolic characterization of amniotic fluids of fetuses with enlarged nuchal translucency. J Perinat Med: Santoru ML; 2019.
Shan J, Xie T, Xu J, Zhou H, Zhao X. Metabolomics of the amniotic fluid: is it a feasible approach to evaluate the safety of Chinese medicine during pregnancy? Journal of applied toxicology : JAT. 2019;39(1):163–71.
Zwemer LM, Bianchi DW. The amniotic fluid transcriptome as a guide to understanding fetal disease. Cold Spring Harb Perspect Med. 2015;5(4).
Hui L, Wick HC, Moise KJ Jr, Johnson A, Luks F, Haeri S, et al. Global gene expression analysis of amniotic fluid cell-free RNA from recipient twins with twin-twin transfusion syndrome. Prenat Diagn. 2013;33(9):873–83.
Edlow AG, Vora NL, Hui L, Wick HC, Cowan JM, Bianchi DW. Maternal obesity affects fetal neurodevelopmental and metabolic gene expression: a pilot study. PLoS One. 2014;9(2):e88661.
Slonim DK, Koide K, Johnson KL, Tantravahi U, Cowan JM, Jarrah Z, et al. Functional genomic analysis of amniotic fluid cell-free mRNA suggests that oxidative stress is significant in Down syndrome fetuses. Proc Natl Acad Sci U S A. 2009;106(23):9425–9.
Massingham LJ, Johnson KL, Scholl TM, Slonim DK, Wick HC, Bianchi DW. Amniotic fluid RNA gene expression profiling provides insights into the phenotype of turner syndrome. Hum Genet. 2014;133(9):1075–82.
Zwemer LM, Nolin SL, Okamoto PM, Eisenberg M, Wick HC, Bianchi DW. Global transcriptome dysregulation in second trimester fetuses with FMR1 expansions. Prenat Diagn. 2017;37(1):43–52.
Gomez-Lopez N, Romero R, Xu Y, Miller D, Arenas-Hernandez M, Garcia-Flores V, et al. Fetal T cell activation in the amniotic cavity during preterm labor: a potential mechanism for a subset of idiopathic preterm birth. J Immunol. 2019;203(7):1793–807.
Brosens I, Pijnenborg R, Vercruysse L, Romero R. The "great obstetrical syndromes" are associated with disorders of deep placentation. Am J Obstet Gynecol. 2011;204(3):193–201.
Tsang JCH, Vong JSL, Ji L, Poon LCY, Jiang P, Lui KO, et al. Integrative single-cell and cell-free plasma RNA transcriptomics elucidates placental cellular dynamics. Proc Natl Acad Sci U S A. 2017;114(37):E7786–E95.
Pavlicev M, Wagner GP, Chavan AR, Owens K, Maziarz J, Dunn-Fletcher C, et al. Single-cell transcriptomics of the human placenta: inferring the cell communication network of the maternal-fetal interface. Genome Res. 2017;27(3):349–61.
Vento-Tormo R, Efremova M, Botting RA, Turco MY, Vento-Tormo M, Meyer KB, et al. Single-cell reconstruction of the early maternal-fetal interface in humans. Nat. 2018;563(7731):347–53.
Tarca AL, Romero R, Xu Z, Gomez-Lopez N, Erez O, Hsu CD, et al. Targeted expression profiling by RNA-Seq improves detection of cellular dynamics during pregnancy and identifies a role for T cells in term parturition. Sci Rep. 2019;9(1):848.
Villar J, Cheikh Ismail L, Victora CG, Ohuma EO, Bertino E, Altman DG, et al. International standards for newborn weight, length, and head circumference by gestational age and sex: the newborn cross-sectional study of the INTERGROWTH-21st project. Lancet. 2014;384(9946):857–68.
Tarca AL, Hernandez-Andrade E, Ahn H, Garcia M, Xu Z, Korzeniewski SJ, et al. Single and serial fetal biometry to detect preterm and term small- and large-for-gestational-age neonates: a longitudinal cohort study. PLoS One. 2016;11(11):e0164161.
Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, et al. Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostat. 2003;4(2):249–64.
Carvalho BS, Irizarry RA. A framework for oligonucleotide microarray preprocessing. Bioinformatics. 2010;26(19):2363–7.
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5(10):R80.
Smyth GK. Limma: linear models for microarray data. In: Gentleman R, Carey VJ, Huber W, Irizarry RA, Dudoit S, editors. Bioinform Comput Biol Solut Using R and Bioconductor: Springer; 2012. p. 397–420.
Clark TA, Sugnet CW, Ares M Jr. Genomewide analysis of mRNA processing in yeast using splicing-specific microarrays. Science. 2002;296(5569):907–10.
Blencowe BJ. Alternative splicing: new insights from global analyses. Cell. 2006;126(1):37–47.
Wick HC, Drabkin H, Ngu H, Sackman M, Fournier C, Haggett J, et al. DFLAT: functional annotation for human development. BMC Bioinform. 2014;15:45.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. ProcNatlAcadSciUSA. 2005;102(43):15545–50.
Su AI, Wiltshire T, Batalov S, Lapp H, Ching KA, Block D, et al. A gene atlas of the mouse and human protein-encoding transcriptomes. Proc Natl Acad Sci U S A. 2004;101(16):6062–7.
Hui L, Slonim DK, Wick HC, Johnson KL, Bianchi DW. The amniotic fluid transcriptome: a source of novel information about human fetal development. Obstet Gynecol. 2012;119(1):111–8.
Hagenauer MH, Schulmann A, Li JZ, Vawter MP, Walsh DM, Thompson RC, et al. Inference of cell type content from human brain transcriptomic datasets illuminates the effects of age, manner of death, dissection, and psychiatric diagnosis. PLoS One. 2018;13(7):e0200003.
Tarca AL, Carey VJ, Chen XW, Romero R, Draghici S. Machine learning and its applications to biology. PLoS Comput Biol. 2007;3(6):e116.
Wang ET, Sandberg R, Luo S, Khrebtukova I, Zhang L, Mayr C, et al. Alternative isoform regulation in human tissue transcriptomes. Nature. 2008;456(7221):470–6.
Pique-Regi R, Romero R, Tarca AL, Sendler ED, Xu Y, Garcia-Flores V, et al. Single cell transcriptional signatures of the human placenta in term and preterm parturition. Elife. 2019;8.
Tarca AL, Romero R, Xu Z, Gomez-Lopez N, Erez O, Hsu CD, et al. Targeted expression profiling by RNA-Seq improves detection of cellular dynamics during pregnancy and identifies a role for T cells in term parturition. Sci Rep. 2018; accepted.
Nazarov PV, Muller A, Kaoma T, Nicot N, Maximo C, Birembaut P, et al. RNA sequencing and transcriptome arrays analyses show opposing results for alternative splicing in patient derived samples. BMC Genomics. 2017;18(1):443.
Raj B, Blencowe BJ. Alternative splicing in the mammalian nervous system: recent insights into mechanisms and functional roles. Neuron. 2015;87(1):14–27.
The authors express our heartfelt appreciation to the patients who made this work possible while attending the Center for Advanced Obstetrical Care and Research at the Perinatology Research Branch, Eunice Kennedy Shriver National Institute of Child Health and Development, National Institutes of Health, U.S. Department of Health and Human Services; Wayne State University; and the Detroit Medical Center, Detroit, Michigan. The authors are also grateful to Dan Lot for conducting the RNA extraction at the Applied Genomics Technology Center of Wayne State University in Detroit, Michigan; to Dr. Christopher Krebs for conducting the microarray experiments at the DNA Sequencing Core at the University of Michigan in Ann Arbor, Michigan; and to Maureen McGerty (Wayne State University) for proofreading the manuscript.
This research was supported, in part, by the Perinatology Research Branch, Division of Obstetrics and Maternal-Fetal Medicine, Division of Intramural Research, Eunice Kennedy Shriver National Institute of Child Health and Human Development, National Institutes of Health, U.S. Department of Health and Human Services (NICHD/NIH/DHHS); and, in part, with Federal funds from NICHD/NIH/DHHS under Contract No. HHSN275201300006C. ALT and NGL were also supported by the Wayne State University School of Medicine Perinatal Initiative.
Ethics approval and consent to participate
All patients provided written informed consent, prior to transabdominal amniocentesis and sample collection, for the use of biological specimens as well as clinical and ultrasound data for research purposes. The study protocols were approved by the Institutional Review Boards of Wayne State University and the Eunice Kennedy Shriver National Institute of Child Health and Human Development (NICHD), National Institutes of Health, U. S. Department of Health and Human Services.
Consent for publication
All authors meet the criteria of authorship and have approved the final version of the manuscript.
The authors report no conflicts of interest. Dr. Romero has contributed to this work as part of his official duties as an employee of the United States Federal Government.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Representation of chromosome-specific differential expression between midtrimester and term gestation groups. The outer circle shows the chromosomes, with significant enrichment being marked with * (q < 0.05). The inner circle shows the log2 fold change (term/midtrimester) of differentially expressed genes, positioned based on their genomic coordinates within each chromosome. Blue denotes increase and red denotes decrease in the term group. Values greater than 3.0, in absolute value, were truncated to 3.0 to enhance display
Changes in the expression of tissue-specific signatures with gestational age superposed to patient-specific values. For each tissue, the expression of the top 20 most-specific genes (based on the Gene Atlas dataset) was transformed into a Z-score and averaged in each AF sample (dots). A Robust Locally Weighted Regression and Smoothing Scatterplots (LOESS) model fit through the Z-scores as a function of gestational age is shown using lines. All differentially expressed tissue signatures (term versus midtrimester or linear correlation within the midtrimester group) are shown. AF, amniotic fluid
Changes in the expression of RNA Seq single-cell signatures with gestational age. For each single-cell signature, the expression the all specific genes (based on Tsang et al. ) was transformed into a Z-score and averaged in each AF sample (dots). A Robust Locally Weighted Regression and Smoothing Scatterplots (LOESS) model fit through the Z-scores as a function of gestational age is shown using lines. AF, amniotic fluid
Representation of chromosome-specific differential expression between pregnancies with a male those with a female fetus. The outer circle shows the chromosomes with significant enrichment being marked with * (q < 0.05). The inner circle shows the log2 fold change (male/female) of differentially expressed genes, positioned based on their genomic coordinates within each chromosome. Blue denotes increase and red denotes decrease in the males. Values greater than 3.0 in absolute value were truncated to 3.0
Example of differential splicing with fetal sex for TXLNGY gene. Details as shown in Fig. 5
Correlation of expression changes with maternal obesity between studies. Each dot represents a unique annotated gene. The y axis represents the log2 fold change (obese/lean) obtained in the current study. The x-axis represents the log2 expression of 182 genes detected as present in the current study among those reported as differentially expressed with obesity by Edlow et al. . R: Spearman’s correlation coefficient
Genes differentially expressed with advancing gestational age from midtrimester to term. The table includes the Affymetrix transcript cluster identifier (ID), gene symbol, gene name, ENTREZ database identifier, chromosome and strand information, genomic coordinates (start and stop), number of microarray probes used to measure the gene expression, assignment to coding or non-coding region, log2 expression change (term/midtrimester), nominal p-value, and adjusted p-value (q-value)
Genes differentially expressed with gestational age during midtrimester. Columns are the same as in Table S1, except that the comparison is between late (21.1–24.0 weeks) versus early (16.4–21.0 weeks) midtrimester gestation
Biological processes from DFLAT database associated with gestational age differential expression from midtrimester to term. Count: number of differentially expressed genes associated with the ontology. Size: total number of genes associated with the ontology that were tested; odds ratio of enrichment based on a Fisher’s exact test; q: adjusted p-value
MSigDB database canonical pathways associated with gestational age differential expression from midtrimester to term. Columns are as shown in the Table S3 legend
Gene Atlas-based tissues and cell types associated with increased expression in term vs midtrimester samples. Columns are as shown in the Table S3 legend
Gene Atlas-based tissues and cell types associated with increased expression in midtrimester vs term samples. Columns are as shown in the Table S3 legend
Genes differentially spliced with gestational age (term versus midtrimester). The table includes results for the single most-significant exon or exon junction based on splicing index analysis. The columns include the Affymetrix transcript cluster identifier (ID), gene symbol, gene name, probeset identifier, exon splicing index, corresponding p-value, adjusted p-value, name of the alternative splicing event type, and adjusted p-value based on the diffSplice algorithm
Biological processes from DFLAT database associated with gestational-age differential splicing but not with differential expression from midtrimester to term. Columns are as shown in the Table S3 legend
MSigDB database canonical pathways associated with gestational-age differential splicing but not with differential expression from midtrimester to term. Columns are as shown in the Table S3 legend
Genes differentially expressed with fetal sex. Columns are the same as in Table S1, except that the comparison is between pregnancies with a male versus a female fetus
About this article
Cite this article
Tarca, A.L., Romero, R., Pique-Regi, R. et al. Amniotic fluid cell-free transcriptome: a glimpse into fetal development and placental cellular dynamics during normal pregnancy. BMC Med Genomics 13, 25 (2020). https://doi.org/10.1186/s12920-020-0690-5
- Cell-free RNA
- Differential splicing
- Differential expression
- Single-cell genomic signature
- Tissue-specific signature
- Gestational age
- Maternal obesity
- Fetal sex