Large-scale transcriptome-wide profiling of microRNAs in human placenta and maternal plasma at early to mid gestation

Melanie D. Smith, Katherine Pillman, Tanja Jankovic-Karasoulos, Dale McAninch, Qianhui Wan, K. Justinian Bogias, Dylan McCullough, Tina Bianco-Miotto, James Breen, Claire T. Roberts
doi: https://doi.org/10.1101/2020.08.19.20177873

Abstract

Background MicroRNAs (miRNAs) are increasingly seen as important regulators of placental development and opportunistic biomarker targets. Given the difficulty in obtaining samples from early gestation and subsequent paucity of the same, investigation of the role of miRNAs in early gestation human placenta has been limited. To address this, we generated miRNA profiles using 96 placentas from presumed normal pregnancies, across early gestation, in combination with matched profiles from maternal plasma. Placenta samples range from 6-23 weeks’ gestation, a time period that includes placenta from the early, relatively low but physiological (6-10 weeks’ gestation) oxygen environment, and later, physiologically normal oxygen environment (11-23 weeks’ gestation).

Results We identified 637 miRNAs with expression in 86 samples (after removing poor quality samples), showing a clear gestational age gradient from 6-23 weeks’ gestation. We identified 374 differentially expressed (DE) miRNAs between placentas from 6-10 weeks’ versus 11-23 weeks’ gestation. We see a clear gestational age group bias in miRNA clusters C19MC, C14MC, miR-17∼92 and paralogs, regions that also include many DE miRNAs. Proportional change in expression of placenta-specific miRNA clusters was reflected in maternal plasma.

Conclusion The presumed introduction of oxygenated maternal blood into the placenta (between ∼10-12 weeks’ gestation) changes the miRNA profile of the chorionic villus, particularly in placenta-specific miRNA clusters. Data presented here comprise a clinically important reference set for studying early placenta development and may underpin the generation of minimally invasive methods for monitoring placental health.

Background

The human placenta is a complex, rapidly developing and highly regulated organ shared between two genetically unique individuals, the mother and her fetus. Its function is vital for the transfer of nutrients, gases and wastes between maternal and fetal circulations, with dysfunction leading to placental insufficiencies and pregnancy complications [1,2]. Essential placental functions are mediated by the syncytiotrophoblast, large multinucleated cells covering the chorionic villi and providing a large surface area for bi-directional exchange of molecules between fetal and maternal circulations.

Early development of the human placenta occurs in a relatively low, but physiological, oxygen environment due to the occluding presence of extravillous cytotrophoblasts (EVTs) in the uterine spiral arterioles. The placenta transitions to that which is usually considered a physiologically normal, oxygen environment over several weeks from 10 weeks’ gestation when EVTs become dislodged permitting maternal blood to flow into the placental intervillous space [1,3]. This period of rising oxygen tension is a critical time for placental development, as an appropriate response to the accompanying burst of oxidative stress is crucial to the ongoing success of the pregnancy [3]. Robust characterisation of important regulatory mechanisms surrounding normal placental development can give insight into pregnancy health, and the potential development of pregnancy complications [2,4]. Identification of biomarkers of ectopic oxygenation or placental dysfunction leading to pregnancy complications is an important step in monitoring pregnancy health.

MicroRNAs (miRNAs) are an important class of RNA molecules that can help us to understand regulation of developmental pathways within human tissues. As such they are prime targets for biomarker research. miRNAs are short, approximately 20-24 nucleotide (nt) lengths of single-stranded, non-coding ribonucleic acid (RNA) originating in non-coding regions of the genome [5] and comprehensively reviewed in Bartel (2018) [6]. They most commonly function as post-transcriptional regulators of gene expression either by translational repression or targeted degradation by cleavage of mRNA transcripts in the cytoplasm [7,8]. miRNAs target mRNA regulatory sequences that are often located in the 3’ untranslated region (UTR) of the mRNA transcript with potential for multiple conserved target sites within a single mRNA [5,9]. Approximately half of known mammalian miRNAs are found within intergenic regions of the genome, with the remaining miRNAs mostly transcribed from the open reading frame of other genes, including long non-coding RNA (lncRNA) [10–12]. Adding to this complexity, miRNA may be transcribed as a single, or polycistronic unit [11].

The critical role of the placenta in successful pregnancy is clear and there is also evidence linking developmental programming of chronic adult disease such as heart disease, diabetes and obesity with the placental phenotype in utero [4,13,14]. However, there are currently no non-invasive tests to determine which women are likely to develop pregnancy complications in routine use, nor is there a comprehensive reference set of miRNA expression during early pregnancy that can be used for development of pregnancy health biomarkers. In this study, we performed high-throughput sequencing to generate robust miRNA expression profiles for 96 placentas and matched maternal plasma samples between 6-23 weeks’ gestation, creating a comprehensive set of miRNA expression profiles for presumably healthy placental function. We identified differential miRNA expression between 6-10 and 11-23 weeks’ gestation which could reflect physiological changes occurring in this important time period. Using placenta-associated miRNA clusters found on chromosome 14, 19, 13, 7 and X, we also investigated how maternal plasma reflects placental profiles to enable future biomarker discoveries.

Results

Characterising a comprehensive miRNA profile of chorionic villous tissue across early gestation to characterise the human placental miRNA profile across early to mid gestation, Illumina NextSeq 75 bp single-end read sequencing was performed on miRNA libraries from chorionic villous samples obtained from 96 singleton pregnancies following elective pregnancy terminations (44 female and 52 male bearing pregnancies). An average of 18.72 million reads were sequenced per sample (range ∼10-29 million reads; Additional file 1: Table S1). miRNA counts were generated by mapping to human miRBase V21 [15]. Libraries with low sequencing depth were removed, including samples PAC0131 (14,595 reads) and PAC0071 (removed due to low mapping after initial read count of 1,819,847 reads). Subsequent unsupervised clustering (Additional file 2: Figure S1) using principal component analysis (PCA) showed eight additional samples (PAC0041, PAC0039, PAC0045, PAC0008, PAC0035, PAC0024, PAC0034, PAC0006) that segregated from the remaining 86 samples, most likely due to the presence of adjacent non placental villous tissue. This hypothesis was tested using matched DNA methylation profiles which confirmed the presence of decidual tissue [16], leading to removal of these samples from subsequent analyses. After filtering, miRNAs were aligned to the human GRCh37 reference genome and counts produced from 4,665 miRBase annotations [15], identifying an initial set of 2,032 expressed miRNAs. To remove sequencing noise, miRNAs with <5 reads across all samples were removed. A total of 1,422 miRNAs were identified with non-zero expression means across all 86 samples, representing a standard known miRBase miRNA expression set for early to mid gestation human placenta. After accounting for sequencing batch effects, very low abundance miRNAs were discarded by converting to log2 CPM values and removing miRNAs with expression <2.5 CPM in ≥29 samples (see methods), leaving a total of 637 robust miRNAs for downstream profile analyses. An additional 588 candidate de novo miRNA sequences were predicted through miRDeep2 analysis.

Of the 637 identified miRNAs, the top ten miRNAs with the highest mean expression across all samples in placenta were miR-30d-5p, miR-125a-5p, miR-517a-3p, miR-199a-3p, miR-26b-5p, miR-26a-5p, let-7a-5p, miR-21-5p, miR-126-3p and miR-516b-5p (Table 1a). For a full list of miRNAs identified and profiled in this study see Additional File 1: Table S2. A principal component analysis (PCA) of all 637 miRNAs (Figure 1) indicates a clear gradient across early gestation. The PCA identified 23.2% variance accounted for in PC1, with the top 10 miRNAs contributing to PCA dimension one (PC1) being: miR-519c-3p, miR-328-3p, miR-519b-3p, miR-20a-5p, miR-19b-3p, miR-501-3p, miR-485-5p, miR-485-5p, miR-3605-3p and miR-515-5p (Additional file 1: Table S3). Interestingly both -3p and -5p arms of the miR-20a transcript, previously associated with placental angiogenesis [17], contributed highly to the gestational age change associated with PC1.

View this table:
Table 1: miRNA with highest expression in early gestation placenta.

Global miRNA expression was calculated using the average expression across early gestation (6-23 weeks’) and log2 transformed. (a) Top ten miRNA with highest average global expression. (b) Constitutively highly expressed miRNA with membership in the top ten by highest average expression if calculated independently for the 6-10 weeks’ and 11-23 weeks’ gestation groups.

Figure 1:
Figure 1: Unsupervised clustering of placental chorionic villous samples after batch correction.

The PCA plot indicates a clear gestational age gradient from early (left) to later gestation (right) in the first dimension. This trend reflects the gestational age signature identified in our miRNA expression data.

Differentially expressed miRNAs across early gestation

Given the importance of oxygen tension in orchestrating the growth and differentiation of the placenta between 6–23 weeks’ gestation, we separated the 86 samples in our cohort by their gestation relative to the presumed introduction of oxygenated maternal blood flow into the intervillous space. Using prior knowledge we designated that ≤10 weeks’ (6-10 weeks’) gestation be considered pre- and >10 weeks’ (11-23 weeks’) gestation be considered post-initiation of maternal blood flow into the intervillous space [3]. This enabled us to identify changes in miRNA expression across this time-point, and to independently assess highly expressed miRNA in each group. When miRNAs from all 86 samples were analysed according to expression prior to and post 10 weeks’ gestation, after accounting for potentially confounding factors such as maternal age, maternal BMI, maternal smoking status, we identified 13 constitutively expressed miRNAs with consistent high expression (log2 CPM 14.35 - 15.60) in both gestation groups (Table 1b). Notably, four of the top twenty most highly expressed miRNAs are from the placenta-specific chromosome 19 cluster, previously identified as being highly expressed in placenta [18–20].

After identifying the highly expressed miRNAs in placenta, we then investigated the change in expression that occurs after 10 weeks’ gestation once maternal blood flow into the placenta is initiated and presumably when oxygen tension begins to rise. Differential expression analysis identified 374 significantly (FDR <0.05) different miRNAs, with 163 down-regulated and 211 up-regulated in the 11-23 weeks’ gestation placenta compared to 6-10 weeks’ gestation (Figure 2a; Additional file 1: Table S4). Strongly downregulated miRNAs (negative log2FC) include miR-4483, miR-129-5p, miR-124-3p and miR-122-5p, and highest up-regulated (positive log2FC) include miR-4645-3p, miR-195-3p, miR-137, miR-139-3p, miR-6715b-3p, miR-3927-3p and miR-1269b. Differentially expressed miRNAs had predominantly low expression (Figure 2b) with an average log2 expression >10 only observed in miR-9-5p (log2FC of −1.43, Figure 2c) and let-7a-5p (log2FC of 1.22, Figure 2d). We see that for both let-7 and miR-9, the 5’ and 3’ arms follow a similar expression trajectory. These mature fragments were originally transcribed as a single transcript before further processing (Figure 2e).

Figure 2:
Figure 2: Differential expression analysis between 6-10 weeks’ and 11-23 weeks’ gestation placenta.

(a) Volcano plot of differential expression. Linear regression identified 26 down-regulated with a log2FC <-1 (red) and 26 up-regulated with a log2FC >1 (green) miRNA between gestation groups. (b) MA plot of log2 fold change as a function of log2 average expression. Box plots of average expression as a function of gestational age for miR-9-5p/3p (c) and let-7a-5p/3p (d) show a clear changing pattern of expression from early to mid gestation. (e) Schematic of a representative miRNA hairpin structure indicates the origin of the 5’ and 3’ arms of the mature miRNA transcript.

Conserved placenta-associated miRNA clusters display variable expression between early and mid gestation

In vertebrate species miRNAs are significantly enriched in clusters, with ∼50% of miRNA clusters being the result of random or non-local duplications. The placenta-associated clusters analysed herein included chromosome 19 miRNA cluster commonly known as C19MC (19q13.41) [21], chromosome 14 cluster (C14MC; 14q.32) [22], and chromosome 13 cluster (miR-17∼92; 13q31.3) with its paralogs the chromosome X cluster miR-106a∼363 (Xq26.2), and chromosome 7 cluster miR-106b∼25 (7q22.1) [23]. We identified the expression of 48 mature miRNAs from C19MC, 77 from C14MC and 20 from miR-17∼92 and paralogs. Of the 48 mature miRNAs identified in the C19MC cluster, 34 miRNA had significant (FDR <0.05) differential expression, that were down regulated in 11-23 weeks’ compared to 6-10 weeks’ gestation placenta (Additional file 1: Table S5). C19MC cluster expression, as a proportion of total miRNA transcripts, was reduced after 10 weeks’ gestation (Figure 3a). Interestingly, while miRNAs in this cluster have been demonstrated to be transcribed by RNA-Pol II as a single polycistron [24], the individual members display a range of expression levels (Figure 3b).

Figure 3:
Figure 3: The proportion of C19MC transcripts decreased after 10 weeks’ gestation in placenta.

The proportion of combined C19MC member expression as a function of total miRNA expression is decreased in 11-23 weeks’ compared to 6-10 weeks’ gestation (violin plot) (a) with expression of individual cluster members tending to decrease across early to mid gestation (b).

The C14MC cluster, that contains 77 mature miRNAs, had 56 miRNAs with significant (FDR <0.05) differential expression. 55 of these miRNAs were up-regulated, while only one miRNA (miR-410-5p; FDR 1.99 x 10×3) was down-regulated in the 11-23 weeks’ gestation group (Additional file 1: Table S6). C14MC-member transcript expression, as a proportion of total expression, increased across early to mid gestation in placenta, further differing from C19MC finding by being more centered around the group mean (Figure 4a). Similar to C19MC cluster members, we found a range of expression levels amongst the C14MC members (Figure 4b) that are also thought to be transcribed as a single polycistronic transcript, suggesting more complex regulation of primary miRNA transcripts processing [25]. The miRNA miR-412-5p, previously observed as highly expressed in first trimester placenta [26], is notable in that it appears to have a bimodal distribution of expression (Figure 4c) that, upon further investigation, was not found to be related to sample characteristics such as fetal sex, maternal age, maternal BMI or maternal smoking status.

Figure 4:
Figure 4: The proportion of C14MC transcripts increased in 11-23 weeks’ compared to 6-10 weeks’ gestation placenta.

The proportion of C14MC member expression as a function of total miRNA expression is increased after 10 weeks’ gestation compared to 6-10 weeks’ gestation (violin plot) (a). Individual cluster members’ expression is increased across early gestation in 55 of the 56 miRNAs (b). miR-412-5p appears to have a bimodal distribution which is unexplained by fetal sex, maternal age or maternal smoking status (c).

Of the 20 mature miRNAs in the miR-17∼92 cluster and its paralogs, 17 had significant (FDR <0.05) differential expression with 15 being down- and 2 being up-regulated in the 11-23 weeks’ gestation group compared to the 6-10 weeks’ (Additional file 1: Table S7). Investigation into the proportion of placenta expression originating from the miR-17∼92 cluster and its paralogs revealed a pattern similar to that of the C19MC member transcripts with the proportion of transcripts identified in placenta tissue as a function of total transcripts decreasing across early gestation (Figure 5a, b, c). We identified a variation in expression patterns between each of the three paralogs (Figure 5d), with chromosome 13 expression decreasing across early to mid gestation.

Figure 5:
Figure 5: The proportion of miR-17∼92 (Chr 13) and paralogs 106b∼25 (Chr 7) and 106a∼363 (Chr X) transcripts decreased in 11-23 weeks’ compared to 6-10 weeks’ gestation placenta.

Violin plot showing decreased proportion of miR-17∼92 member expression as a function of total miRNA expression after 10 weeks’ gestation (a, b, c) with a general trend of decreased expression of individual cluster members across early gestation (d). Although chromosome 13 and chromosome X appear correlated, expression of miR-17∼92 (chromosome 13) and its paralogs miR-106a∼363 (chromosome 7) and miR-106b∼25 (X chromosome) does not appear to follow a common expression pattern in placental tissue as we would expect if all 3 clusters are identically regulated by C-Myc.

Whilst we find significant (FDR <0.05) differential expression in all of the above miRNA clusters, the log2FC between 6-10 weeks’ and 11-23 weeks’ gestation groups is small (C19MC log2FC −1.02 to −0.19; C14MC log2FC 0.23 to 0.94 and log2FC 0.67; miR-17∼92 log2FC −0.2 to −1.13 and log2FC 0.84 to 0.89), with changes in the proportional expression of each cluster as a function of total miRNA expression suggesting an important regulatory mechanism for these clusters in placenta.

Characterising a comprehensive profile of placental and endogenous miRNA in maternal plasma across early gestation

To characterise the miRNA profile across early to mid gestation in maternal plasma, we sequenced 96 matched maternal plasma samples (hereafter referred to as plasma) from the same women sampled at the time of pregnancy termination using the same sequencing protocol as detailed for the chorionic villous tissue. miRNA sequencing resulted in approximately 17.77 million 75 bp single-end reads per sample in plasma (range 9,986,347-41,344,361 reads: Additional file 3: Table S8) across two sequencing runs. After initial quality control and sequence analysis, we identified a number of samples with aberrant expression of miR-451a and miR-16 indicating potential for red blood cell miRNA contamination. We performed a subsequent analysis using the delta-delta CT method (ΔΔCt) to calculate the relative fold change of gene expression of samples using quantitative polymerase chain reaction (qPCR). The ΔΔCt value of miR-23a minus miR-451 is examined, as miR-23a is known to be stable in samples affected by haemolysis whilst miR-451 is known to vary with contamination. This method identified 11 samples (PAC0033, PAC0034, PAC0041, PAC0048, PAC0050, PAC0051, PAC0054, PAC0056, PAC0062, PAC0084 and PAC0105) as more likely to contain haemolysis contamination, a common issue in plasma research. Despite potential issues involving haemolysis, we aimed to establish the presence of placenta specific miRNAs within maternal plasma by including all available samples regardless of whether they had matched placenta data or not.

Alignment and filtering were performed using the same protocol as for chorionic villous tissue. The number of unique miRNA species identified per sample varied, ranging from 1022 miRNA species in PAC0077 to 355 miRNA species in PAC0010. Prior to filtering we identified 1775 unique mature miRNAs across all plasma samples, with a core set of 790 high confidence miRNAs present after filtering. Of these, 531 were also present in placenta samples, 259 unique to plasma, and 106 miRNAs unique to placenta (Additional file 3: Table S9). Unsupervised clustering (PCA) of all 790 high confidence miRNAs identified in plasma shows no gestational age gradient for maternal plasma from 6-23 weeks’ gestation (Additional file 2: Figure S2).

Given the variability of miRNA in plasma, we hypothesised that detection of placental miRNAs in plasma may be restricted to miRNA present in high concentrations in matched placenta. We therefore investigated whether ten miRNAs (Table 1a) with the highest abundance in placenta were present in plasma and found that all ten miRNA were at high abundance levels in matched plasma (Figure 6, Table 2a). Many of the most highly expressed plasma miRNA are red blood cell (RBC) associated.

View this table:
Table 2: miRNA by average expression in maternal plasma.

Average miRNA expression was calculated using the average normalised CPM expression across early to mid gestation (6-23 weeks’) and log2 transformed. (a) All top 10 miRNA with highest abundance in placenta were also detected at high abundance in plasma. Two miRNA (miR-517a-3p and miR-516b-5p) originate from the C19MC. (b) Top 10 miRNA by abundance in maternal plasma. RBC-associated miRNA identified herein have been previously identified in plasma samples including adult plasma and umbilical cord plasma but have not previously been profiled extensively in maternal plasma in early to mid gestation.

Figure 6:
Figure 6: Abundance of maternal plasma miRNAs, identified as highly expressed in chorionic villous tissue, from 6-23 weeks’ gestation.

All ten miRNAs identified as the most highly expressed in chorionic villous tissue were identified as highly abundant in matched plasma samples. Abundance of the C19MC miRNA members, miR-516b-5p and miR-517a-3p, appear to increase across early to mid gestation but these were not statistically significant (FDR <0.05).

Of the 790 miRNAs identified in plasma, the top ten miRNAs with the highest global expression were miR-16-5p, miR-486-5p, miR-92a-3p, let-7a-5p, let-7b-5p, miR-451a, let-7f-5p, miR-25-3p, let-7i-5p, miR-93-5p (Table 2b). These were also all highly expressed in placental tissue (average global expression of 8.95-14.47 Log2 normalised CPM). miR-141, miR-200b, miR-139 and miR-184, previously identified in third trimester plasma, were also identified at high abundance in our early gestation plasma, with miR-141-3p abundance increasing across early gestation, as shown in previous studies [27].

We investigated plasma levels of the 21 miRNAs identified as both differentially and highly expressed in the matched placentas (Table 3). Whilst all 21 placental miRNAs were also detected in plasma, their origin is not necessarily placenta-specific, with cellular enrichment linked to various systemic compartments [28]. For example, five of these miRNAs have cellular enrichment in immune response, another three miRNAs in pluripotent and embryonic stem cells, and one miRNA is enriched in extra-embryonic stem cells. The latter may originate from the fetal membranes or other extra-embryonic tissue such as chorionic villi.

View this table:
Table 3: The abundance in plasma of miRNA differentially expressed in matched placenta.

Average miRNA expression was calculated using the average CPM expression across early to mid gestation (6-23 weeks’) and log2 transformed.

We further investigated the abundance of plasma miRNAs from three miRNA clusters previously identified as highly placenta specific or placenta associated [21–23]. We found 44 mature miRNA from the C19MC cluster, 48 miRNA from the C14MC cluster and 20 miRNA from the miR-17∼92 cluster and its paralogs in plasma. These cluster miRNAs have previously been the target of biomarker research in pregnancy health but to our knowledge these are not being used clinically.

The change in relative proportion of expression by miRNA cluster members in plasma across the 10-11 weeks’ gestation threshold is difficult to detect with C14MC and C19MC members each representing less than 1% of the miRNA abundance quantified in plasma. As such the proportional change identified in placenta is not entirely reflected in plasma. Abundance levels of the miR-17∼92 cluster and its paralogs represent approximately 0.03% of the total miRNA identified in plasma.

Discussion

In this study we have assessed the miRNA profile of 96 human placenta and matched plasma samples in early to mid gestation (6-23 weeks’), generating a placental tissue dataset using high-throughput sequencing that is substantially larger than previous reports that used microarray and qPCR profiling. We derive robust and comprehensive profiles that define miRNA expression at high-resolution, applying stringent bioinformatic protocols for the analysis of sequencing data, akin to large atlas studies.

Differentially expressed miRNAs reflect the dynamically changing utero-placental environment during early to mid gestation

In order to determine the potential association between placental miRNA expression and pregnancy health, particularly the expression of highly abundant miRNAs across early gestation, we examined more closely a period in early gestation involving a steep rise in oxygen tension in the intervillous space (IVS). Rising oxygen tension occurs between 10-12 weeks’ gestation, at onset of maternal blood flow into the IVS, resulting in an oxidative burst in placenta [3]. Previous studies have addressed miRNA response to hypoxia and the introduction of oxygen [31–35], including how ectopic oxygen tension may contribute to diseases of pregnancy [31,36–38]; and reviewed in [39]. The number and extensive gestation range of our samples enabled us to analyse the differential expression of miRNAs prior to, during and post the critical 10-11 weeks’ gestational time point. This has potential to lead to identification of placental mechanisms for the remediation of oxidative stress at this crucial time in pregnancy.

That the placenta develops initially in a condition of low (but physiological) oxygen tension has been known for some time [3] and it has been suggested that this is important for the protection of the early fetus from the effects of oxidative stress such as tissue damage from oxygen free radicals [40]. miRNAs are uniquely poised to play a significant role in the mediation of gene expression throughout the early change of oxygen tension in the placenta by employing post transcriptional regulatory mechanisms to rapidly and precisely fine tune the placental environment in response to oxidative stress. We identified upregulation of let-7b/e5p, miR-21, miR-23a/b-3p, miR-24-3p and miR-199a-5p, also associated with HIF1α, in the 6-10 weeks’ gestation group. All of these have previously been shown to be upregulated by hypoxia. We also found lower abundance of let-7a and miR-101-3p in our 6-10 weeks’ gestation group, both of which have previously been shown to be downregulated by hypoxia [41–45]; and reviewed in [46]. miRNAs identified in this study likely play an important role in the mitigation and regulation of hypoxia in the placenta. This may prove important to our understanding of diseases of pregnancy such as preeclampsia and early pregnancy loss.

The importance of miRNA function in the early placenta is not limited to the mitigation of oxidative stress but also in cellular proliferation, trophoblast invasion and cellular differentiation, which are important to the placenta and also to the success of the pregnancy. let-7a-5p expression in first trimester placental explants has been implicated in a reduction in cytotrophoblast proliferation. In our differential expression analysis we identified let-7a-5p as up-regulated in the 11-23 weeks’ gestation compared to the 6-10 weeks’ gestation group. Our detailed week-by-week expression data suggests an incremental increase across early to mid gestation which together with previous term data is consistent with let-7a-5p’s proposed role in regulation of placental growth by the reduction of cytotrophoblast proliferation [47]. The miR-34 family is known as a tumour suppressor miRNA family due to their synergistic effects with p53 [48]. Consistent with this, we identified up-regulation of all six miR-34 family members reflecting a reduction in cellular proliferation as the pregnancy progresses and placental growth trajectory begins to slow. Similarly, miR-378a-5p was previously found to enhance cell survival and promote trophoblast migration and invasion [49]. Our data showed significant down-regulation of the miR-34 family in 11-23 weeks’ gestation compared to the 6-10 weeks’ gestation placenta potentially linking down-regulation of these miRNAs with the reduced trophoblast migration as pregnancy progresses. miR-200 family members are also known to have an established role in placental function, in particular in the cellular transition from epithelial to mesenchymal (EMT) phenotype which is a prerequisite for extravillous trophoblast invasion of the maternal decidua. miR-200 family members play a well-known role in EMT through a feedback loop with the ZEB family of transcription factors. We identified seven members of the miR-200 family, six of which are down-regulated in the 11-23 weeks’ gestation group. The miR-200 family is known to be highly expressed in epithelial cell types [50] and thus its down-regulation is consistent with development of the vascularised mesenchymal core in chorionic villi in the second trimester.

Proportional change in expression miRNA clusters play a role in placenta gene regulation during development

miRNA clusters are known to play an important role in the placenta, with these regions thought to have co-evolved and to act cooperatively to repress target genes [51]. Interestingly, while miRNA clusters are typically transcribed as a single polycistronic transcript, we identified large intra-sample variances in the expression levels of cluster members (Additional file 2, Figure S3) hinting at additional post transcriptional regulatory mechanisms. There are a number of miRNA clusters expressed either exclusively, or preferentially in placenta [18,52,53]. In particular, the C19MC, C14MC and miR-17∼92 clusters have been shown to regulate phenotypic and functional diversity in frequently used cell lines and isolated primary trophoblasts, and in the immune system, by which they may play a role in immune tolerance to paternal antigens [19,54,55].

C19MC is expressed almost exclusively in the placenta but also in certain tumors and undifferentiated cells [21,56]. C19MC is a large placenta specific imprinted cluster located on chromosome 19q13.41 and mono-allelically expressed from the paternally inherited chromosome [57]. This cluster has been difficult to study as its expression is both placenta and primate specific [57] and hence lacks an ortholog in mouse [22]. Importantly, C19MC members are expressed at lower than trophoblast levels in JEG-3 and JAR cells, two cell lines derived from choriocarcinoma [58], and are not expressed in the commonly used HTR8/SVneo cell line [53,56] highlighting the need for primary placental tissue in research. Consistent with our findings, previous research has shown that C19MC constitutes a large proportion of total miRNA transcripts (15% at term), with our data additionally showing a gestational decline from 27% at 6-10 weeks’ gestation to 23% at 11-23 weeks’ gestation [20,56]. In addition, there was a significant (FDR < 0.05) down-regulation in 34 C19MC cluster members in 11-23 weeks’ placenta compared to the 6-10 weeks’ gestation. Interestingly, a previous report showed that 46 of the 47 miRNAs from the C19MC cluster were up-regulated in third trimester villous trophoblasts compared to first trimester [19]. It is possible that this miRNA cluster exhibits a dynamic gestational expression pattern but perhaps more likely is the possibility that differences between studies can be explained by differences in expression between assessment of trophoblast versus whole chorionic villous tissue. The latter is more physiologically relevant and fits our question.

A recent miR-517a/c in situ hybridisation analysis localised expression from these C19MC miRs to cytotrophoblast and syncytiotrophoblast in chorionic villi, proximal cytotrophoblasts in anchoring villi, decreased expression in distal cytotrophoblasts with further reduction in extra-villous cytotrophoblasts. This finding highlights variable expression of the C19MC cluster in different placental cell types, and shows loss of expression of this cluster as trophoblasts differentiate to extravillous phenotypes post EMT [59]. Another report concurs showing higher expression of C19MC members in villous trophoblasts compared to extravillous trophoblasts suggestive of a role for the C19MC cluster miRNA in attenuation of EVT migration through the direct targeting of mRNA transcripts related to cellular movement [60]. Mouillet et al. (2015) suggested a role for C19MC in cellular differentiation or the maintenance of pluripotency [61]. Further investigation into the proportional change in expression from pregnancies complicated by placental pathology may help to identify miRNA species with biomarker potential in plasma. However, low expression levels in some maternal plasma samples have frustrated our efforts in this regard.

C14MC, located at the imprinted DLK1-DIO3 domain on human chromosome 14q32 is the largest known human miRNA cluster comprising 52 miRNA genes expressed from the maternally inherited chromosome [53]. The C14MC cluster is found exclusively in eutherian mammals and has been suggested to be essential to the evolution of this lineage, with expression in humans demonstrating a strong bias towards brain, placenta and some embryonic tissues [52]. Our analysis identified expression from 84 mature miRNAs in the C14MC cluster with the proportion of member transcripts increased in the 11-23 weeks’ gestation group compared to 6-10 weeks’ gestation. Previous studies investigating C14MC expression, in both whole villous tissue and in primary cytotrophoblasts, have reported a decrease in expression from first to third trimester. Using whole villous tissue, Gu et al. (2013) identified 11 miRNA from this cluster down-regulated in third compared to first trimester, and proposed a link between both C19MC and C14MC in immune suppression and innate/adaptive immune response in the mother. In primary cytotrophoblast, Morales-Prieto et al. (2012) identified 34 miRNA from this cluster which were also down-regulated in third compared to first trimester [53,62]. Whilst the discrepancy between our data and Morales-Prieto et al. (2013) can be explained by the use of different biological source (single cell type versus whole tissue), that between our data and Gu et al. (2013) is not so clear but may be due to differences in technology, the number of cluster members identified, or that our samples range up to 23 weeks’ gestation and do not include term. It is plausible that rather than a steady decline from early gestation to term, the expression of C14MC members increases across early gestation, before decreasing at term. Consistent with a recent report by Wommack et al. (2018) who profiled circulating miRNAs in maternal plasma and found an inverse relationship between placenta-specific clusters C19MC and C14MC we found this in our tissue samples.

Finally, in the analysis of placenta-associated clusters we investigated the expression of miR-17∼92 and paralogs miR-106a∼363 and miR-106b∼25, identifying 19 mature miRNAs. The proportion of total expression from miR-17∼92 members and paralogs are less than 1% per cluster, with all three clusters showing an overall decrease in the proportion of total expression across early to mid gestation. Using the miRNA sequences provided by Kumar et al. (2013), we determined the 5’ or 3’ origin of the transcripts by cross-checking the given sequences against miRBase v21[54,63], and were thus able to directly compare their microarray data with our sequencing data (detailed in Additional File 3: Table S11).

Plasma miRNA sequencing offers a potential source of non-invasive biomarkers for detecting placenta health

In addition to establishing the placental miRNA profile across early to mid gestation and, given the potential for placental health biomarker discovery using our access to matched maternal plasma, we also aimed to analyse miRNA profiles in maternal plasma sampled in women at the same time as placenta. The identification of placenta derived, or pregnancy associated miRNA in plasma offers a potential window into minimally invasive monitoring of pregnancy progression.

Working with miRNA in plasma is known to be difficult due to potential inclusion of red blood cells or other cellular debris [64–66]. Red blood cell contamination may mask any signal we hoped to detect from the placenta as haemolysis may perturb miRNA species and count distributions due to overloading of red blood cell miRNAs [64,66,67]. In our study we see evidence that 11 of our samples were likely affected by haemolysis. This common difficulty with plasma data, led us to assess the detection of miRNAs of interest revealing that the 10 most abundant placental miRNAs were reflected by high abundance levels in matched plasma. Furthermore, these plasma miRNAs in our study show strong concordance with previous miRNA plasma studies [65,68]. Differential expression and clear gestational age gradient clustering pattern in placenta miRNAs between 6-10 and 11-23 weeks’ gestation were not particularly reflected by maternal plasma, indicating that plasma does not reflect the changes in the miRNA profile in placenta across early to mid gestation. miRNAs that were abundant in placenta were not always found to be abundant in plasma. In particular, miRNAs enriched in Pluripotent stem cell/embryonic stem cell were of low abundance in plasma compared to placenta. With the exception of placenta-specific miRNAs, cellular enrichment could be linked to various other systemic compartments. For example, five miRNAs have cellular enrichment in immune response cells, whilst another three miRNAs have enrichment in pluripotent and embryonic stem cells, and one miRNA is enriched in extra-embryonic stem cells which may originate from the fetal membranes or other extra-embryonic tissue such as chorionic villi. Of particular interest, miR-99a-5p enrichment includes enrichment in mast cells which are involved in inflammatory response and allergic reactions. We see an apparent increase in abundance of miR-99a-5p across early to mid gestation, and some additional increase in male-compared with female-bearing pregnancies which may be indicative of the increased maternal immune response to the male conceptus [69].

Complementary to the investigation of miRNA clusters in placenta, we also investigated the expression of three highly placenta-associated miRNA clusters and confirmed the abundance of these cluster miRNAs in plasma, with the proportion of cluster member expression revealing biomarker potential in plasma. One pathway for the transfer of placental molecules into the maternal circulation is via exosomes. It has previously been reported that members of the C19MC are the predominant miRNA species expressed in exosomes released from primary human trophoblasts [56]. We were able to confirm the presence of 39 C19MC, 48 C14MC and 20 miR-17∼92 miRNA in matched plasma, with their presence showing promise for the use of miRNA as biomarkers for placental health. Perhaps more importantly, week-by-week examination of miRNAs across early gestation, as described in this study, are key elements for developing pregnancy surveillance measures in the future, especially when combined with disease-specific profiles associated with preeclampsia, gestational diabetes and preterm birth. Our investigation here of miRNAs as biomarkers of placental, and by proxy pregnancy health have identified many areas for further investigation. With the implementation of machine learning algorithms, cell free DNA work and placenta specific analysis standards, much more can be expected from this field.

Conclusion

Appropriate placental growth and development is essential for pregnancy success and miRNAs play roles in mediating these. We investigated placental and maternal plasma miRNAs using matched samples taken after elective termination of otherwise normal pregnancies and found 374 significantly (FDR <0.05) differentially expressed miRNA between 6-10 and 11-23 weeks’ gestation placenta.

The research conducted herein is the most comprehensive attempt to accurately profile the miRNA landscape of the human placenta, in part, because of the use of miRNA sequencing rather than microarray or qPCR. This dataset represents unprecedented access to a large number of human placenta samples from 6-23 weeks of gestation. Together these data are an important reference set for miRNA expression of early to mid gestation placental development and function, at time points rarely seen in previous placental sequencing work. These provide an important resource for the placental biology field. Further investigation using maternal plasma in early gestation samples from pregnancies with known outcomes will provide further insights into the search for miRNAs as biomarkers of pregnancy health and disease.

Methods

Sample collection

Placental chorionic villous tissue samples were obtained with informed, written consent from women undergoing elective terminations of otherwise healthy pregnancies. Samples of peripheral blood (6-9 mL) were collected from the same women into standard EDTA blood tubes at the time of termination and stored on ice until processed. Whole blood underwent centrifugation at 800 x g for 15 minutes at 4°C before plasma removal and then spun for a further 15 minutes. All samples were stored at −80°C until further processing. Prior to termination, accurate gestational age was determined using transvaginal ultrasonography. Termination samples were collected from the Pregnancy Advisory Centre (PAC), Woodville, South Australia. All placenta samples were assessed for gross morphology before inclusion in the study.

RNA extraction and library preparation

All placenta tissue was collected post termination and processed as soon as possible after collection and placed into RNALater (Thermo Fisher) within 15 minutes. Total RNA was isolated using a modified protocol of the RNeasy mini plus kit (Qiagen, Hilden, Germany) whereby the RW1 buffer, a proprietary component of the RNeasy Kit that eliminates small RNAs, is replaced with 100% ethanol. The sample is then subjected to size fractionation for small RNAs before sequencing. For plasma, miRNA was isolated from 200 μL plasma samples using the QIAGEN miRNA serum/plasma kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. All samples were stored at −80°C until further processing. miRNA real-time qPCR for assessment of haemolysis was conducted by QIAGEN Genomic Services (Qiagen, Hilden, Germany). 2ml RNA was reverse transcribed in 10 ml reactions using miRCURY LNA Kit (QIAGEN version5). Each RT was performed including an artificial RNA spike-in (UniSp6). cDNA was diluted 50 x and assayed in 10 ml PCR reactions according to the protocol for miRCURY LNA miRNA PCR; each miRNA was assayed once by qPCR using assays for miR-23a, miR-30c, miR-103, miR-142-3p, and miR-451. In addition to these miRNA assays, the RNA spike-ins were assayed. The amplification was performed in a LightCycler® 480 Real-Time PCR System (Roche) in 384 well plates. The amplification curves were analysed using in-house software, both for determination of the Cq (by the 2nd derivative method) and for melting curve analysis.

All placenta samples were genotyped for fetal sex using high resolution melt curve analysis of the gene that encodes amelogenin. Amelogenin is found on both the X (AMELX) and Y (AMELY) chromosomes, the X allele features a 3 bp deletion in exon 3 allowing the identification of samples with only X chromosomes (female) or X and Y chromosomes (male) prior to sequencing. Forward primer 5’-CCCTGGGCTCTGTAAAGAATAGTG-3’, reverse primer 5’- ATCAGAGCTTAAACTGGGAAGCTG-3. qPCR was performed using SSOFast EvaGreen Supermix (Bio-Rad, CA), primers at 250 nM final concentration, 5 ng of DNA per reaction, on a Bio-Rad CFX384 Real-Time PCR System. Cycling conditions: initial denaturation 98°C 30s, 40 cycles of 98°C for 5s and 60°C for 5s. High resolution melt curve analysis was performed from 65°C to 85°C with a 0.2°C increment every 10s. Melt curve between 65°C and 70°C was analysed using Bio-Rad Precision melt software (Bio-Rad, CA) to identify sex genotypes. Library preparation and sequencing was performed by Qiagen (Valencia, CA) using the QIAseq miRNA Library Kit and QIAseq miRNA 48 Index IL kits as per manufacturer’s instructions. Amplified cDNA libraries underwent single-end sequencing by synthesis (Illumina 1.9).

High-throughput sequencing analysis

Alignment and analysis of miRNA sequencing data were processed from raw FASTQ files using the bcbio-nextgen pipeline [70]. Briefly, adapter detection and trimming were performed using Atropos [71]. Alignment was performed using STAR [72] using the human reference genome build GRCh37 [73] (https://www.ncbi.nlm.nih.gov/grc). Quality control metrics were assessed using FastQC [74] (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) to check for per base sequence quality, sequence length distribution and duplication levels, and summarised using multiQC [75]. miRNA were clustered and collapsed using SeqCluster [76], and individual miRNAs detected using SeqBuster [70]. Annotation was performed using the miraligner [70] with miRBase version 21.0 [15,63] (http://www.mirbase.org) as the reference database.

Differential expression analysis

All profile and expression analyses were conducted in the R statistical environment (v.3.3.2), using the edgeR (v.3.16.5) and limma [77] (v.3.30.11) R/Bioconductor packages. edgeR [78] was used to filter miRNA with low expression and normalise for library composition bias. Batch effects, introduced by sequencing samples in multiple sequencing runs, were evident after normalisation and corrected using a limma batch correction method [77]. All samples were then normalised using the Trimmed Mean of M values (TMM).

Sample-weights and log transformation was performed using limma package [79] with the voom function used to estimate the mean-variance relationship between individual observations and then applied to the normalised log-counts data. Differential expression analysis including, moderated F-statistic evaluation, adjusted p-value estimation and Log Fold Change (LFC) analysis were performed using a moderated t-test [80] with Benjamini-Hochberg (BH) multiple hypothesis test corrections [81]. After adjustment, expression of miRNAs was considered significantly different at FDR ≤ 0.05. The workflow, source code and input files associated with this research available at (https://github.com/mxhp75/earlyPlacentamiRNA_SeqProfile.git).

Data Availability

The dataset(s) supporting the conclusion of this article are available in NCBI's Gene Expression Omnibus [82] and are accessible through GEO Series accession number GSE151362

https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE151362

https://github.com/mxhp75/earlyPlacentamiRNA_SeqProfile.git

Declarations

Ethics approval and consent to participate

Ethics approval for the collection of placenta tissue and blood from women undergoing elective pregnancy termination between 6–23 weeks’ gestation was provided under HREC/16/TQEH/33, by The Queen Elizabeth Hospital Human Research Ethics Committee (TQEH/LMH/MH).

Consent for publication

Availability of data and materials

The dataset(s) supporting the conclusion of this article are available in NCBI’s Gene Expression Omnibus [82] and are accessible through GEO Series accession number GSE151362 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE151362).

Competing interests

All authors declare no competing interests.

Funding

This research is supported by NIH NICHD R01 HD089685-01 Maternal molecular profiles reflect placental function and development across gestation PI Roberts, an Australian Government Research Training Program (RTP) Scholarship awarded to MDS, a National Health and Medical Research Council Investigator Grant (GNT1174971) awarded to CTR and a Matthew Flinders Professorial Fellowship awarded to CTR and funded by Flinders University. JB is supported by the James & Diana Ramsay Foundation.

Authors contributions

CTR created the concept. MS, CTR, JB and TBM conceived and developed experimental plans. TJK, DMcA and DMcC performed experiments. MS and JB analysed the sequencing data with help from KB and QW. Manuscript written by MS, JB and CTR with help from TBM and KP. All authors read and approved the final version of the manuscript.

Author information

1Robinson Research Institute, 2Adelaide Medical School, 3Bioinformatics Hub, 4The University of South Australia, 5School of Agriculture Food and Wine; University of Adelaide, 6South Australian Health & Medical Research Institute (SAHMRI), 7Flinders Health and Medical Research Institute, Flinders University

Additional file

Additional_File_1.xlxs

Table S1. Read statistics for 96 placenta samples miRNA-Seq data. Table S2. miRNA expression data from 86 placenta chorionic villous samples. Expression values calculated as log2 CPM across 86 samples passing filters (all_mean), across samples of 6-10 weeks’ gestation (≤10_mean) and for samples 11-23 weeks’ gestation (>10_mean). Relative ranks for mean expression is also given as all_rank, ≤10_rank and >10_rank respectively. Table S3. Principal component analysis results showing the contribution of each miRNA across PC1 and PC2. Table S4. Differential expression analysis of data from 86 placenta chorionic villous samples identified 374 significant (FDR <0.05) different miRNA between 6-10 weeks’ and 11-23 weeks’ gestation groups. Table 5. Differential expression analysis of data from 86 placenta chorionic villous samples revealed 34 C19MC cluster miRNAs down-regulated in the 6-10 weeks’ compared to the 11-23 weeks’ gestation groups. Table S6. Differential expression analysis of data from 86 placenta chorionic villous samples revealed 56 C14MC cluster miRNAs with significant (FDR <0.05) differential expression between the 6-10 weeks’ compared to the 11-23 weeks’ gestation. Table S7. Differential expression analysis of data from 86 placenta chorionic villous samples revealed 17 miR∼17-92 cluster miRNAs with significant (FDR <0.05) differential expression between the 6-10 weeks’ and 11-23 weeks’ gestation groups.

Additional_File_2.docx

Figure S1. Using principal component analyses of the miRNA-Seq data from 94 placenta samples, we identified 8 samples (named) that segregated from the remaining 86 samples. This separation on PC1 is most likely due to the presence of surrounding non-villous tissue. Figure S2. The PCA plot of batch corrected data shows no gestational age gradient for maternal plasma from 6-23 weeks’ gestation. The four samples in the upper right quadrant were sequenced in the same sequencing run and may represent a technical artefact in line with previous difficulties in plasma sequencing but have not been removed from further analysis. Figure S3. Heatmap of C19MC cluster members plotted from 5’ (left) to 3’ (right), and from early (top) to later (bottom) gestation clearly identifies different expression levels for the polycistronic members of this cluster.

Additional_File_3.xlxs

Table S8. Read statistics for 96 plasma samples miRNA-Seq data. Table S9. miRNA abundance data from 96 maternal plasma samples. Expression values calculated as log2 CPM across 96 samples. miRNA ranked 1:790 by abundance. Table S10. Differential expression analysis of data from 86 placenta chorionic villous samples using gestational age (6-23 weeks’) as a continuous variable. Table S11. We compared our T1 vs T2 differential expression data with those of Kumar et al. (2013) who compared expression of miR-17∼92, and its paralogs miR-106a∼363 and miR-106b∼25 clusters and identified a down-regulation of cluster members post syncytialisation.

Acknowledgements

We wish to acknowledge the generosity of the women who donated their placenta for our research. Without them, this research would not be possible. We also acknowledge valuable input from QIAGEN Genomic Services.

Footnotes

  • Colours changed in the heatmaps (Figure 3b and Figure 4b) to add clarity around the gestational age annotation.

References

  1. 1.↵
    Gude NM, Roberts CT, Kalionis B, King RG. Growth and function of the normal human placenta. Thromb Res. 2004;114:397–407.
    OpenUrlCrossRefPubMedWeb of Science
  2. 2.↵
    Sibai B, Dekker G, Kupferminc M. Pre-eclampsia. Lancet. Elsevier; 2005. p. 785–99.
  3. 3.↵
    Jauniaux E, Watson AL, Hempstock J, Bao YP, Skepper JN, Burton GJ. Onset of maternal arterial blood flow and placental oxidative stress: A possible factor in human early pregnancy failure. Am J Pathol. American Society for Investigative Pathology; 2000;157:2111–22.
    OpenUrl
  4. 4.↵
    Henriksen T, Clausen T. The fetal origins hypothesis: Placental insufficiency and inheritance versus maternal malnutrition in well-nourished populations. Acta Obstet Gynecol Scand. Wiley/Blackwell (10.1111); 2002;81:112–4.
    OpenUrl
  5. 5.↵
    Bartel DP. MicroRNAs: Genomics, Biogenesis, Mechanism, and Function [Internet]. Cell. 2004. p. 281–97. Available from: http://linkinghub.elsevier.com/retrieve/pii/S0092867404000455
  6. 6.↵
    Bartel DP. Metazoan MicroRNAs. Cell. 2018;173:20–51.
    OpenUrlCrossRefPubMed
  7. 7.↵
    Winter J, Jung S, Keller S, Gregory RI, Diederichs S. Many roads to maturity: MicroRNA biogenesis pathways and their regulation [Internet]. Nature Cell Biology. Nature Publishing Group; 2009. p. 228–34. Available from: http://www.nature.com/articles/ncb0309-228
  8. 8.↵
    Ameres SL, Zamore PD. Diversifying microRNA sequence and function [Internet]. Nature Reviews Molecular Cell Biology. Nature Publishing Group; 2013. p. 475–88. Available from: http://www.nature.com/articles/nrm3611
  9. 9.↵
    Bartel DP. MicroRNAs: Target Recognition and Regulatory Functions. Cell. 2009;136:215–33.
    OpenUrlCrossRefPubMedWeb of Science
  10. 10.↵
    Saini HK, Griffiths-Jones S, Enright AJ. Genomic analysis of human microRNA transcripts. Proceedings of the National Academy of Sciences. National Academy of Sciences; 2007;104:17719–24.
    OpenUrl
  11. 11.↵
    Kim VN, Han J, Siomi MC. Biogenesis of small RNAs in animals [Internet]. Nature Reviews Molecular Cell Biology. 2009. p. 126–39. Available from: http://dx.doi.org/10.1038/nrm2632
  12. 12.↵
    Dykes IM, Emanueli C. Transcriptional and Post-transcriptional Gene Regulation by Long Non-coding RNA [Internet]. Genomics, Proteomics and Bioinformatics. Elsevier; 2017. p. 177–86. Available from: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=PMC5487525
  13. 13.↵
    Morley R. Fetal origins of adult disease [Internet]. Seminars in Fetal and Neonatal Medicine. Elsevier; 2006. p. 73–8. Available from: http://dx.doi.org/10.1016/j.siny.2005.11.001
  14. 14.↵
    Burton GJ, Fowden AL, Thornburg KL. Placental Origins of Chronic Disease. Physiol Rev. American Physiological Society Bethesda, MD; 2016;96:1509–65.
    OpenUrl
  15. 15.↵
    Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enright AJ. miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acids Res. Oxford University Press; 2006;34:D140–4.
    OpenUrlCrossRefPubMedWeb of Science
  16. 16.↵
    Wan Q, Yiner-Lee Leemaqz S, Pederson SM, McCullough D, McAninch DC, Jankovic-Karasoulos T, et al. Quality control measures for placental sample purity in DNA methylation array analyses. Placenta. 2019;88:8–11.
    OpenUrl
  17. 17.↵
    Wang W, Feng L, Zhang H, Hachy S, Satohisa S, Laurent LC, et al. Preeclampsia Up-Regulates Angiogenesis-Associated MicroRNA (i.e., miR-17, −20a, and −20b) That Target Ephrin-B2 and EPHB4 in Human Placenta. J Clin Endocrinol Metab. 2012;97:E1051–9.
    OpenUrlCrossRefPubMedWeb of Science
  18. 18.↵
    Luo S-S, Ishibashi O, Ishikawa G, Ishikawa T, Katayama A, Mishima T, et al. Human Villous Trophoblasts Express and Secrete Placenta-Specific MicroRNAs into Maternal Circulation via Exosomes1. Biol Reprod. Oxford University Press; 2009;81:717–29.
    OpenUrlCrossRefPubMedWeb of Science
  19. 19.↵
    Morales-Prieto DM, Chaiwangyen W, Ospina-Prieto S, Schneider U, Herrmann J, Gruhn B, et al. MicroRNA expression profiles of trophoblastic cells. Placenta. 2012;33:725–34.
    OpenUrlCrossRefPubMedWeb of Science
  20. 20.↵
    Carreras-Badosa G, Bonmat A, Ortega FJ, Mercader JM, Guindo-Martnez M, Torrents D, et al. Dysregulation of placental miRNA in maternal obesity is associated with pre-and postnatal growth. J Clin Endocrinol Metab. 2017;102:2584–94.
    OpenUrlPubMed
  21. 21.↵
    Bentwich I, Avniel A, Karov Y, Aharonov R, Gilad S, Barad O, et al. Identification of hundreds of conserved and nonconserved human microRNAs. Nat Genet. Nature Publishing Group; 2005;37:766–70.
    OpenUrlCrossRefPubMedWeb of Science
  22. 22.↵
    Morales-Prieto DM, Ospina-Prieto S, Schmidt A, Chaiwangyen W, Markert UR. Elsevier Trophoblast Research Award Lecture: Origin, evolution and future of placenta miRNAs. Placenta [Internet]. 2014;35. Available from: http://dx.doi.org/10.1016/j.placenta.2013.11.017
  23. 23.↵
    Olive V, Jiang I, He L. Mir-17-92, a cluster of miRNAs in the midst of the cancer network. Int J Biochem Cell Biol. 2010;42:1348–54.
    OpenUrlCrossRefPubMedWeb of Science
  24. 24.↵
    Bortolin-Cavaille ML, Dance M, Weber M, Cavaille J. C19MC microRNAs are processed from introns of large Pol-II, non-protein-coding transcripts. Nucleic Acids Res [Internet]. 2009; Available from: http://dx.doi.org/10.1093/nar/gkp205
  25. 25.↵
    Seitz H, Royo H, Bortolin ML, Lin SP, Ferguson-Smith AC, Cavaillé J. A large imprinted microRNA gene cluster at the mouse Dlk1-Gtl2 domain. Genome Res. 2004;14:1741–8.
    OpenUrlAbstract/FREE Full Text
  26. 26.↵
    Schmidt JK, Block LN, Golos TG. Defining the rhesus macaque placental miRNAome: Conservation of expression of placental miRNA clusters between the macaque and human. Placenta. Elsevier Ltd; 2018;65:55–64.
    OpenUrl
  27. 27.↵
    Chim SSC, Shing TKF, Hung ECW, Leung T-Y, Lau T-K, Chiu RWK, et al. Detection and characterization of placental microRNAs in maternal plasma. Clin Chem. Clinical Chemistry; 2008;54:482–90.
    OpenUrl
  28. 28.↵
    Zhao C, Sun X, Li L. Biogenesis and function of extracellular miRNAs [Internet]. ExRNA. 2019. Available from: http://dx.doi.org/10.1186/s41544-019-0039-4
  29. 29.
    Ludwig N, Leidinger P, Becker K, Backes C, Fehlmann T, Pallasch C, et al. Distribution of miRNA expression across human tissues. Nucleic Acids Res. 2016;44:3865–77.
    OpenUrlCrossRefPubMed
  30. 30.
    Lizio M, Harshbarger J, Abugessaisa I, Noguchi S, Kondo A, Severin J, et al. Update of the FANTOM web resource: High resolution transcriptome of diverse cell types in mammals. Nucleic Acids Res. Oxford University Press; 2017;45:D737–43.
    OpenUrlCrossRefPubMed
  31. 31.↵
    Mouillet JF, Chu T, Hubel CA, Nelson DM, Parks WT, Sadovsky Y. The levels of hypoxia-regulated microRNAs in plasma of pregnant women with fetal growth restriction. Placenta. W.B. Saunders; 2010;31:781–4.
    OpenUrl
  32. 32.
    Mouillet JF, Chu T, Nelson DM, Mishima T, Sadovsky Y. MiR-205 silences MED1 in hypoxic primary human trophoblasts. The FASEB Journal. Federation of American Societies for Experimental Biology; 2010;24:2030–9.
    OpenUrl
  33. 33.
    Shen J, Xia W, Khotskaya YB, Huo L, Nakanishi K, Lim S-O, et al. EGFR modulates microRNA maturation in response to hypoxia through phosphorylation of AGO2. Nature. 2013;497:383–7.
    OpenUrlCrossRefPubMedWeb of Science
  34. 34.
    Camps C, Saini HK, Mole DR, Choudhry H, Reczko M, Guerra-Assunção J, et al. Integrated analysis of microRNA and mRNA expression and association with HIF binding reveals the complexity of microRNA expression regulation under hypoxia. Mol Cancer. BioMed Central; 2014;13:28.
    OpenUrl
  35. 35.↵
    Mori A, Nishi H, Sasaki T, Nagamitsu Y, Kawaguchi R, Okamoto A, et al. HLA-G expression is regulated by miR-365 in trophoblasts under hypoxic conditions. Placenta. 2016;45:37–41.
    OpenUrl
  36. 36.↵
    Chelbi ST, Vaiman D. Genetic and epigenetic factors contribute to the onset of preeclampsia. Mol Cell Endocrinol. Elsevier; 2008;282:120–9.
    OpenUrlCrossRefPubMedWeb of Science
  37. 37.
    1. Zenclussen AC,
    Cronqvist T, Saljé K, Familari M, Guller S, Schneider H, Gardiner C, et al. Syncytiotrophoblast Vesicles Show Altered micro-RNA and Haemoglobin Content after Ex-vivo Perfusion of Placentas with Haemoglobin to Mimic Preeclampsia. Zenclussen AC, editor. PLoS One. Public Library of Science; 2014;9:e90020.
    OpenUrlPubMed
  38. 38.↵
    Li Q, Pan Z, Wang X, Gao Z, Ren C, Yang W. miR-125b-1-3p inhibits trophoblast cell invasion by targeting sphingosine-1-phosphate receptor 1 in preeclampsia. Biochem Biophys Res Commun. 2014;453:57–63.
    OpenUrlCrossRefPubMed
  39. 39.↵
    Poston L, Raijmakers MTM. Trophoblast Oxidative Stress, Antioxidants and Pregnancy Outcome—A Review. Placenta. W.B. Saunders; 2004;25:S72–8.
    OpenUrl
  40. 40.↵
    Burton GJ, Hempstock J, Jauniaux E. Oxygen, early embryonic metabolism and free radical-mediated embryopathies. Reprod Biomed Online. 2003;6:84–96.
    OpenUrlCrossRefPubMed
  41. 41.↵
    Hua Z, Lv Q, Ye W, Wong C-KA, Cai G, Gu D, et al. MiRNA-directed regulation of VEGF and other angiogenic factors under hypoxia. PLoS One. 2006;1:e116.
    OpenUrlCrossRefPubMed
  42. 42.
    Hebert C, Norris K, Scheper MA, Nikitakis N, Sauk JJ. High mobility group A2 is a target for miRNA-98 in head and neck squamous cell carcinoma. Mol Cancer. 2007;6:5.
    OpenUrlCrossRefPubMed
  43. 43.
    Kulshreshtha R, Ferracin M, Wojcik SE, Garzon R, Alder H, Agosto-Perez FJ, et al. A microRNA signature of hypoxia. Mol Cell Biol. 2007;27:1859–67.
    OpenUrlAbstract/FREE Full Text
  44. 44.
    Rane S, He M, Sayed D, Vashistha H, Malhotra A, Sadoshima J, et al. Downregulation of MiR-199a derepresses hypoxia-inducible factor-1α and sirtuin 1 and recapitulates hypoxia preconditioning in cardiac myocytes. Circ Res. 2009;104:879–86.
    OpenUrlAbstract/FREE Full Text
  45. 45.↵
    Gonsalves CS, Kalra VK. Hypoxia-mediated expression of 5-lipoxygenase-activating protein involves HIF-1alpha and NF-kappaB and microRNAs 135a and 199a-5p. J Immunol. 2010;184:3878–88.
    OpenUrlAbstract/FREE Full Text
  46. 46.↵
    Nallamshetty S, Chan SY, Loscalzo J. Hypoxia: a master regulator of microRNA biogenesis and activity. Free Radic Biol Med. 2013;64:20–30.
    OpenUrlCrossRefPubMed
  47. 47.↵
    Forbes K, Farrokhnia F, Aplin JD, Westwood M. Dicer-dependent miRNAs provide an endogenous restraint on cytotrophoblast proliferation. Placenta [Internet]. 2012; Available from: http://dx.doi.org/10.1016/j.placenta.2012.03.006
  48. 48.↵
    Zhang L, Liao Y, Tang L. MicroRNA-34 family: a potential tumor suppressor and therapeutic candidate in cancer. J Exp Clin Cancer Res. 2019;38:53.
    OpenUrl
  49. 49.↵
    Luo L, Ye G, Nadeem L, Fu G, Yang BB, Honarparvar E, et al. MicroRNA-378a-5p promotes trophoblast cell survival, migration and invasion by targeting Nodal. J Cell Sci. 2012;125:3124–32.
    OpenUrlAbstract/FREE Full Text
  50. 50.↵
    Humphries B, Yang C. The microRNA-200 family: Small molecules with novel roles in cancer development, progression and therapy. Oncotarget. Impact Journals LLC; 2015;6:6472–98.
    OpenUrl
  51. 51.↵
    Wang Y, Luo J, Zhang H, Lu J. MicroRNAs in the Same Clusters Evolve to Coordinately Regulate Functionally Related Genes. Mol Biol Evol. Oxford University Press; 2016;33:2232–47.
    OpenUrlCrossRefPubMed
  52. 52.↵
    Glazov EA, McWilliam S, Barris WC, Dalrymple BP. Origin, evolution, and biological role of miRNA cluster in DLK-DIO3 genomic region in placental mammals. Mol Biol Evol. 2008;25:939–48.
    OpenUrlCrossRefPubMedWeb of Science
  53. 53.↵
    Morales-Prieto DM, Ospina-Prieto S, Chaiwangyen W, Schoenleben M, Markert UR. Pregnancy-associated miRNA-clusters. J Reprod Immunol. Elsevier Ireland Ltd; 2013;97:51–61.
    OpenUrlCrossRefPubMed
  54. 54.↵
    Kumar P, Luo Y, Tudela C, Alexander JM, Mendelson CR. The c-Myc-Regulated MicroRNA-17 92 (miR-17 92) and miR-106a 363 Clusters Target hCYP19A1 and hGCM1 To Inhibit Human Trophoblast Differentiation. Mol Cell Biol. American Society for Microbiology Journals; 2013;33:1782–96.
    OpenUrlAbstract/FREE Full Text
  55. 55.↵
    Kuo G, Wu CY, Yang HY. MiR-17-92 cluster and immunity [Internet]. Journal of the Formosan Medical Association. 2019. Available from: http://dx.doi.org/10.1016/j.jfma.2018.04.013
  56. 56.↵
    Donker RB, Mouillet JF, Chu T, Hubel CA, Stolz DB, Morelli AE, et al. The expression profile of C19MC microRNAs in primary human trophoblast cells and exosomes. Mol Hum Reprod. Narnia; 2012;18:417–24.
    OpenUrlCrossRefPubMedWeb of Science
  57. 57.↵
    Noguer-Dance M, Abu-Amero S, Al-Khtib M, Lefèvre A, Coullin P, Moore GE, et al. The primate-specific microRNA gene cluster (C19MC) is imprinted in the placenta. Hum Mol Genet. 2010;19:3566–82.
    OpenUrlCrossRefPubMedWeb of Science
  58. 58.↵
    Lee CQE, Gardner L, Turco M, Zhao N, Murray MJ, Coleman N, et al. What Is Trophoblast? A Combination of Criteria Define Human First-Trimester Trophoblast. Stem Cell Reports. 2016;6:257–72.
    OpenUrl
  59. 59.↵
    Mong EF, Yang Y, Akat KM, Canfield J, VanWye J, Lockhart J, et al. Chromosome 19 microRNA cluster enhances cell reprogramming by inhibiting epithelial-to-mesenchymal transition. Sci Rep. 2020;10:3029.
    OpenUrl
  60. 60.↵
    Xie L, Mouillet JF, Chu T, Parks WT, Sadovsky E, Knöfler M, et al. C19MC microRNAs regulate the migration of human trophoblasts. Endocrinology. 2014;155:4975–85.
    OpenUrlCrossRefPubMed
  61. 61.↵
    Mouillet J-FJ-F, Ouyang Y, Coyne CB, Sadovsky Y. MicroRNAs in placental health and disease. Am J Obstet Gynecol. 2015;213:S163–72.
    OpenUrlCrossRefPubMed
  62. 62.↵
    Gu Y, Sun J, Groome LJ, Wang Y. Differential miRNA expression profiles between the first and third trimester human placentas. AJP: Endocrinology and Metabolism. 2013;304:E836–43.
    OpenUrl
  63. 63.↵
    Kozomara A, Griffiths-Jones S. MiRBase: Annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res [Internet]. 2014;42. Available from: http://dx.doi.org/10.1093/nar/gkt1181
  64. 64.↵
    Pritchard CC, Cheng HH, Tewari M. MicroRNA profiling: Approaches and considerations [Internet]. Nature Reviews Genetics. NIH Public Access; 2012. p. 358–69. Available from: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=PMC4517822
  65. 65.↵
    Blondal T, Jensby Nielsen S, Baker A, Andreasen D, Mouritzen P, Wrang Teilum M, et al. Assessing sample and miRNA profile quality in serum and plasma or other biofluids [Internet]. Methods. Academic Press; 2013. p. 164–9. Available from: https://www.sciencedirect.com/science/article/pii/S1046202312002551
  66. 66.↵
    1. Kiechl S
    Cheng HH, Yi HS, Kim Y, Kroh EM, Chien JW, Eaton KD, et al. Plasma Processing Conditions Substantially Influence Circulating microRNA Biomarker Levels. Kiechl S, editor. PLoS One. Public Library of Science; 2013;8:e64795.
    OpenUrlCrossRefPubMed
  67. 67.↵
    1. Pfeffer S
    Kirschner MB, Kao SC, Edelman JJ, Armstrong NJ, Vallely MP, van Zandwijk N, et al. Haemolysis during Sample Preparation Alters microRNA Content of Plasma. Pfeffer S, editor. PLoS One. 2011;6:e24145.
    OpenUrlCrossRefPubMed
  68. 68.↵
    Brennan GP, Vitsios DM, Casey S, Looney A-M, Hallberg B, Henshall DC, et al. RNA-sequencing analysis of umbilical cord plasma microRNAs from healthy newborns. PLoS One. 2018;13:e0207952.
    OpenUrl
  69. 69.↵
    Clifton VL. Review: Sex and the human placenta: mediating differential strategies of fetal growth and survival. Placenta. Elsevier; 2010;31 Suppl:S33–9.
    OpenUrlCrossRefPubMedWeb of Science
  70. 70.↵
    Pantano L, Estivill X, Martí E. SeqBuster, a bioinformatic tool for the processing and analysis of small RNAs datasets, reveals ubiquitous miRNA modifications in human embryonic cells. Nucleic Acids Res [Internet]. 2009;38. Available from: https://watermark.silverchair.com/gkp1127.pdf?token=AQECAHi208BE49Ooan9kkhW_Ercy7Dm3ZL_9Cf3qfKAc485ysgAAAbQwggGwBgkqhkiG9w0BBwagggGhMIIBnQIBADCCAZYGCSqGSIb3DQEHATAeBglghkgBZQMEAS4wEQQM1TUUngRW8zyuVolBAgEQgIIBZ7gj8AvZ9tPmZnDNEGCVulK4NZusRxNE3WbfR5NrHrMCKeH
  71. 71.↵
    Didion JP, Martin M, Collins FS. Atropos: specific, sensitive, and speedy trimming of sequencing reads. PeerJ. PeerJ Inc.; 2017;5:e3720.
    OpenUrl
  72. 72.↵
    Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. 2013;29:15–21.
    OpenUrl
  73. 73.↵
    Church DM, Schneider VA, Graves T, Auger K, Cunningham F, Bouk N, et al. Modernizing reference genome assemblies. PLoS Biol. 2011;9:e1001091.
    OpenUrlCrossRefPubMed
  74. 74.↵
    Andrews S. FastQC: a quality control tool for high throughput sequence data [Internet]. 2010. Available from: http://www.bioinformatics.babraham.ac.uk/projects/fastqc/
  75. 75.↵
    Ewels P, Magnusson M, Lundin S, Käller M. MultiQC: Summarize analysis results for multiple tools and samples in a single report. Bioinformatics. 2016;32:3047–8.
    OpenUrlCrossRefPubMed
  76. 76.↵
    Pantano L, Estivill X, Martí E. A non-biased framework for the annotation and classification of the non-mirna small RNA transcriptome. Bioinformatics. 2011;27:3202–3.
    OpenUrlCrossRefPubMedWeb of Science
  77. 77.↵
    Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43:e47.
    OpenUrlCrossRefPubMed
  78. 78.↵
    Robinson MD, McCarthy DJ, Smyth GK. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2009;26:139–40.
    OpenUrlCrossRefPubMedWeb of Science
  79. 79.↵
    Law CW, Chen Y, Shi W, Smyth GK. Voom: Precision weights unlock linear model analysis tools for RNA-seq read counts. Genome Biol. BioMed Central; 2014;15:R29.
    OpenUrl
  80. 80.↵
    Smyth GK. Linear Models and Empirical Bayes Methods for Assessing Differential Expression in Microarray Experiments. Stat Appl Genet Mol Biol. 2004;3:1–25.
    OpenUrlCrossRefPubMed
  81. 81.↵
    Benjamini Y, Hochberg Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J R Stat Soc Series B Stat Methodol. 1995;57:289–300.
    OpenUrlCrossRef
  82. 82.↵
    Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. Narnia; 2002;30:207–10.
    OpenUrlCrossRefPubMedWeb of Science
Back to top
Posted October 13, 2020.
Download PDF

Supplementary Material

Data/Code
Email
Large-scale transcriptome-wide profiling of microRNAs in human placenta and maternal plasma at early to mid gestation
Melanie D. Smith, Katherine Pillman, Tanja Jankovic-Karasoulos, Dale McAninch, Qianhui Wan, K. Justinian Bogias, Dylan McCullough, Tina Bianco-Miotto, James Breen, Claire T. Roberts
medRxiv 2020.08.19.20177873; doi: https://doi.org/10.1101/2020.08.19.20177873
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area