Next Article in Journal
Longitudinal Study of Cognitive Functioning in Adults with Juvenile Idiopathic Arthritis
Next Article in Special Issue
Small Nucleolar Derived RNAs as Regulators of Human Cancer
Previous Article in Journal
Adipose-Derived Stem Cell Exosomes as a Novel Anti-Inflammatory Agent and the Current Therapeutic Targets for Rheumatoid Arthritis
Previous Article in Special Issue
Extraction-Free Absolute Quantification of Circulating miRNAs by Chip-Based Digital PCR
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Maternal Body Mass Index Is Associated with Profile Variation in Circulating MicroRNAs at First Trimester of Pregnancy

by
Kathrine Thibeault
1,
Cécilia Légaré
1,
Véronique Desgagné
1,2,
Frédérique White
3,
Andrée-Anne Clément
1,
Michelle S. Scott
1,
Pierre-Étienne Jacques
3,4,
Renée Guérin
1,2,
Patrice Perron
4,5,
Marie-France Hivert
5,6,7 and
Luigi Bouchard
1,2,4,*
1
Department of Biochemistry and Functional Genomics, Faculty of Medicine and Health Sciences (FMHS), Université de Sherbrooke, Sherbrooke, QC J1K 2R1, Canada
2
Clinical Department of Laboratory Medicine, Centre Intégré Universitaire de Santé et de Services Sociaux (CIUSSS) du Saguenay-Lac-Saint-Jean—Hôpital de Chicoutimi, Saguenay, QC G7H 5H6, Canada
3
Département de Biologie, Faculté des Science, Université de Sherbrooke, Sherbrooke, QC J1K 2R1, Canada
4
Department of Medicine, FMHS, Université de Sherbrooke, Sherbrooke, QC J1K 2R1, Canada
5
Centre de Recherche du Centre Hospitalier Universitaire de Sherbrooke (CR-CHUS), Université de Sherbrooke, Sherbrooke, QC J1H 5N4, Canada
6
Department of Population Medicine, Harvard Medical School, Harvard Pilgrim Health Care Institute, Boston, MA 02115, USA
7
Diabetes Unit, Massachusetts General Hospital, Boston, MA 02114, USA
*
Author to whom correspondence should be addressed.
Biomedicines 2022, 10(7), 1726; https://doi.org/10.3390/biomedicines10071726
Submission received: 2 June 2022 / Revised: 13 July 2022 / Accepted: 14 July 2022 / Published: 18 July 2022
(This article belongs to the Special Issue microRNAs in Health and Disease)

Abstract

:
Many women enter pregnancy with overweight and obesity, which are associated with complications for both the expectant mother and her child. MicroRNAs (miRNAs) are short non-coding RNAs that regulate many biological processes, including energy metabolism. Our study aimed to identify first trimester plasmatic miRNAs associated with maternal body mass index (BMI) in early pregnancy. We sequenced a total of 658 plasma samples collected between the 4th and 16th week of pregnancy from two independent prospective birth cohorts (Gen3G and 3D). In each cohort, we assessed associations between early pregnancy maternal BMI and plasmatic miRNAs using DESeq2 R package, adjusting for sequencing run and lane, gestational age, maternal age at the first trimester of pregnancy and parity. A total of 38 miRNAs were associated (FDR q < 0.05) with BMI in the Gen3G cohort and were replicated (direction and magnitude of the fold change) in the 3D cohort, including 22 with a nominal p-value < 0.05. Some of these miRNAs were enriched in fatty acid metabolism-related pathways. We identified first trimester plasmatic miRNAs associated with maternal BMI. These miRNAs potentially regulate fatty acid metabolism-related pathways, supporting the hypothesis of their potential contribution to energy metabolism regulation in early pregnancy.

1. Introduction

Overweight and obesity (OW/O) are worldwide health issues [1] increasingly affecting women of reproductive age [2,3] and are considered the most common health conditions in pregnancy [4]. The prevalence of obesity among women of reproductive age is 17.8% in Canada [5] and 39.7% in the United States [6]. Maternal OW/O can have serious short- and long-term consequences on the health of the mother and her child. OW/O in pregnancy are associated with an increased risk of gestational diabetes mellitus (GDM), gestational hypertension, preeclampsia (PE) and cesarean delivery [4]. Offspring exposed to maternal OW/O are more likely to suffer from macrosomia and neonatal hyperinsulinemic hypoglycemia [7]. Possibly through fetal metabolic programming (based on the Developmental Origin of Health and Disease-DOHaD), these offspring are also at increased risk of childhood obesity, type 2 diabetes, cardiovascular disease, asthma and neurodevelopmental disorders [8,9]. A better understanding of the pathophysiological mechanisms leading to pregnancy complications associated with maternal OW/O could lead to strategies to prevent to maternal morbidity and possibly lower their consequences on the exposed offspring.
MicroRNAs (miRNAs) are short non-coding, single-stranded, RNA molecules of 19 to 25 nucleotides. They regulate many biological processes by targeting messenger RNAs (mRNAs), leading to a decrease in protein synthesis [10]. Three miRNA clusters are predominantly expressed in the placenta: chromosome 19 miRNA (C19MC), chromosome 14 miRNA (C14MC) and miR-371-3 miRNA clusters. The abundance of miRNAs from C19MC gradually increases in maternity from conception until the end of pregnancy, whereas miRNAs from the C14MC cluster are more abundant at the beginning of pregnancy and their blood levels decrease throughout pregnancy [11,12]. These placental miRNAs may participate in the regulation of maternal physiology by intercellular communication through placental extracellular vesicle secretion in maternal circulation [13]. Indeed, miRNAs are suspected of having an important role in pregnancy, its maintenance as well as adaptation to its very specific physiologic needs [10,11]. Accordingly, dysregulation of the C19MC and C14MC has been associated with pregnancy complications, including PE, intrauterine growth restriction and insulin sensitivity regulation [13,14,15,16,17].
Until now, only a few studies have investigated miRNAs in pregnancy complicated by OW/O [18,19,20,21]. However, none of these studies were done by next-generation sequencing or in plasma in the first trimester of pregnancy. We hypothesized that maternal BMI at the first trimester of pregnancy is associated with variations in circulating levels of miRNAs. Therefore, our objectives were to identify plasmatic miRNAs associated with maternal BMI and the metabolic pathways they potentially regulate.

2. Materials and Methods

2.1. Discovery Cohort: Genetics of Glucose Regulation in Gestation and Growth (Gen3G) Cohort

We selected participants for the discovery step of the study from the Gen3G prospective pregnancy and birth cohort [22]. Briefly, we recruited women in the first trimester of pregnancy (between the 4th and the 16th week), and we followed them until delivery. An oral glucose tolerance test (OGTT—75 g) was performed between the 24th and 28th week of pregnancy for 854 women. For this study, our selection criteria were: women of European descent, 18 years old and older, not taking any medication that influences glycemia, free from pre-gestational diabetes, having singleton pregnancy as well as the availability of a plasma sample at the first trimester of pregnancy (500 µL), anthropometric measures (e.g., maternal BMI), follow-ups of the offspring at 3 and 5 years old, and genetics (Mother: Infinium MEGAEX Array, Illumina; Offspring: whole genome sequencing) and epigenetics (EPIC array) data. A total of 444 women fulfilled these criteria. Before the study began, consent was obtained from all participating women, and all protocols were approved by the Centre Hospitalier Universitaire de Sherbrooke (CHUS) ethics committee.

2.2. Anthropometric Measurements in Gen3G

BMI is used to classify individuals by dividing their weight (kg) by their height (m) squared. This measurement was taken in the first trimester of pregnancy (between the 4th and 16th week) [23], and a description of the BMI measurement has been published [22]. Briefly, weight was measured in kg on a calibrated electronic scale, and height was measured in meters with a wall stadiometer (without shoes). The BMI value was calculated from these data.

2.3. Replication Cohort: Design, Develop, Discover (3D) Cohort Study

A replication analysis was run on 226 participants selected from the 3D prospective birth cohort [24]. Briefly, the participants were recruited during their first trimester of pregnancy in nine different sites across the province of Québec, Canada. For this study, women of European descent with plasma samples (500 µL) at the 1st trimester of pregnancy and at least fasting and the 2 h post-OGTT glucose levels measured between the 24th and 28th week of pregnancy were included. Exclusion criteria included pre-existing diabetes, GDM diagnosed at the first trimester of pregnancy, chronic hypertension, or gestational hypertension diagnosed at the first trimester of pregnancy (Supplementary Figure S1). Biological specimens (e.g., plasma) and anthropometric measurements were collected (measured by research staff). BMI was measured between the 8th and 14th week of pregnancy during the first trimester of pregnancy. Women without BMI data were also excluded. All women gave their informed consent before the beginning of the study according to the Helsinki declaration, and the protocols have been accepted by the ethical committees

2.4. RNA Extraction

We extracted total RNA with the MirVana PARIS kit (ThermoFisher Scientific, Waltham, United States, catalog #AM1556) from 500 µL of plasma collected between the 4th and 16th week of pregnancy following the standard protocol and eluted in 75 µL of nuclease-free water, in random order. We concentrated our RNA samples following the protocol established by Burgos et al. [25]. In brief, RNA was mixed with 35 µL of cold (4 °C) 7 M ammonium acetate solution (ThermoFisher Scientific, Waltham, United States, catalog #02002268) and mixed with 420 µL of chilled (−20 °C) absolute ethanol (Commercial Alcohols, ON, Canada; catalog #P006EAAN). RNA was precipitated overnight at −20 °C and centrifuged at 16,000× g at 4 °C for 30 min. The RNA pellet was washed twice with 200 µL of 80% ethanol and then centrifuged at 16,000× g, 4 °C for 5 min. The RNA pellet was dried for 30 min at room temperature and resuspended in 5 µL of nuclease-free water.

2.5. Library Preparation

We used the Truseq Small RNA Sample Prep kit (Illumina, BC, Canada, catalog #RS-200-0012) for library preparation. Concentrated RNA samples (5 µL) were randomly treated following the standard protocol adapted by Burgos et al. [25]. Briefly, half of the reagents were used for ligation of RNA at the 3′ and 5′ ends, reverse transcription, indexing (1–48, one index per sample) and PCR amplification (15 cycles), to maintain an optimal ratio between RNA and reagents. The libraries were purified by migration on a Novex polyacrylamide TBE Gel, 6% (ThermoFisher Scientific, Waltham, United States, catalog #EC265BOX) by selecting bands between 145–160 bp, eluted in 300 µL of nuclease-free water, and incubated overnight at room temperature and 500 RPM on an incubating microplate shaker (VWR, ON, Canada, catalog #12620-930). Then, the libraries were concentrated by precipitation following a standardized procedure (including a 30 min incubation of the precipitation mix at −80 °C), and the cDNA pellet was suspended in 25 μL of 10 mM Tris-HCI pH 8.5 buffer.

2.6. Library Quality Control and Sequencing

The libraries were sequenced at the McGill University and Génome Québec Innovation Centre (Montréal, Canada), either on a HiSeq 2500 or HiSeq 4000 platform (50 cycles, with 7 cycles indexing read) for the Gen3G samples. Twelve samples were extracted twice and sequenced on both the HiSeq 2500 and HiSeq 4000 platforms that we leveraged during our QC and normalization process to take into account potential technical and batch effects from the different sequencing platforms. Overall, Pearson correlation coefficients between miRNA levels measured on the 2 platforms for these 12 samples were ≥0.94. The data from both platforms were then combined for processing and statistical analysis but adjusted for run and lane as normally recommended [12]. Libraries were quantified by qPCR, equimolarly pooled (HiSeq 2500: 12 libraries with different indexes per lane at a molarity of 7 pM; HiSeq 4000: 20 libraries with different indexes per lane at a molarity of 10 pM), denatured and clustered on single-read Illumina flowcells (catalog # GD-401-3001 and catalog GD-410-1001) according to the manufacturer’s standard protocol.
The 3D replication study samples were sequenced on the NovaSeq 6000 platform at the McGill University and Génome Québec Innovation Centre (Montréal, Canada). Libraries were prepared following the procedure applied to Gen3G samples and quantified using the Kapa Illumina GA with Revised Primers-SYBR Fast Universal kit (Kapa Biosystems, Wilmington, United States), normalized, equimolarly pooled (48 libraries with different indexes per lane at a molarity of 225 pM), denatured and clustered on an Illumina NovaSeq S1 flowcell following the Xp protocol from the manufacturer’s recommendations. The run was performed for 100 cycles in single-end mode.

2.7. Bioinformatics Analysis

We first applied the extracellular RNA processing toolkit (exceRpt) pipeline (version 4.6) [26]. Briefly, exceRpt uses FASTX-Toolkit and FastQC to assess sequencing data quality after removing the adapters and poor quality (Phred score < 20 for 80% or more of the read) sequences. The remaining reads were then aligned to the human genome (GRCh37) and miRbase [27] (version 21) using STAR [28] (version 2.4.2a) alignment algorithm. After performing data visualization of the raw read counts, 9 outliers were excluded from the Gen3G cohort, and 3 outliers were excluded from the 3D cohort.

2.8. Statistical Analysis

Since all the participant characteristics (Gen3G and 3D cohorts) were not normally distributed based on a Shapiro-Wilk test, non-parametric Mann-Whitney U tests were applied to compare the two cohorts. Association between plasmatic miRNA levels and maternal BMI at the first trimester of pregnancy was assessed using the default parameters in DESeq2 R package [29]. The duplicate samples (n = 12) were combined using the collapseReplicates function from DESeq2 package. The statistical model was adjusted for sequencing run and lane, gestational age, maternal age at the first trimester of pregnancy, and parity. A sensitivity analysis was also performed by adding the fetal sex as a covariate to the analysis model. The results remain overall unchanged. MiRNAs were considered significantly associated with maternal BMI with a false discovery rate (FDR) adjusted q-value < 0.05, where the fold change represents the change in the miRNAs normalized read counts per unit of BMI. The EnhancedVolcano [30] package was used for the creation of the volcano plot. For replication, criteria were first the direction of the associations, the fold changes, and finally the nominal p-values (one-tailed). The statistical analyses were all done with R version 4.0.3 in R studio version 1.4.1103.

2.9. Biological Pathway Analysis

The potential biological function of the maternal BMI-associated miRNAs was evaluated using miRPath v.3 software from DIANA tools [31]. In brief, Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathway analysis was done on all the miRNAs that were significantly (q-values < 0.05) associated with maternal BMI at the first trimester of pregnancy and replicated for their fold change value in the 3D cohort. Tarbase 7.0 was selected for miRNA-mRNA interactions as it considers those validated experimentally. The pathway union parameter was used to merge the results with the Fisher exact test enrichment analysis method. The default settings in miRPath v.3 and FDR correction were applied, and the q-value threshold was set to 0.05.

3. Results

3.1. Participant’s Characteristics

Table 1 presents the characteristic of the selected Gen3G and 3D cohorts’ participants. Women from Gen3G were on average 28.48 ± 4.26 years old, had a mean BMI of 25.95 ± 5.98 kg/m2 at the first trimester of pregnancy (mean gestational age: 9.63 ± 2.26 weeks), and had a mean parity of 0.69 ± 0.91 (53.56% nulliparous). On average, women from the 3D cohort were slightly older (30.60 ± 3.90 years old, p < 0.001) and had a greater gestational age in the 1st trimester visit (11.89 ± 1.52 weeks, p < 0.001) and a lower parity (0.47 ± 0.66, p = 0.01; 60.99% nulliparous), whereas their mean BMI (25.39 ± 5.91 kg/m2, p = 0.22) was similar to that of the women from the Gen3G cohort.

3.2. Association between miRNA Levels and Maternal BMI at the First Trimester of Pregnancy

In Gen3G, 2170 miRNAs were expressed in the plasma of pregnant women. A miRNA was considered detected when it has >1 read per participant. A total of 61 miRNAs (FDR-adjusted q-values < 0.05) were associated with maternal BMI in the first trimester of pregnancy (Figure 1). The list of significant miRNAs with their mean normalized read counts, fold changes (FC), unadjusted p-values, and FDR-adjusted q-values are shown in Supplementary Table S1. Among these 61 miRNAs, higher BMI was associated with lower circulating levels for 48 miRNAs and with higher circulating levels for 13 miRNAs. Interestingly, 28 of these miRNAs (46%) are encoded by the C19MC, 1 miRNA by the C14MC, and 3 miRNAs by the miR-371-3 miRNAs cluster. Table 2 shows the top 10 miRNAs with the most significant associations with maternal BMI; for these 10 miRNAs, higher BMI was associated with lower circulating levels.

3.3. Replication of miRNAs Associated with Maternal BMI at the First Trimester of Pregnancy in the 3D Cohort

We conducted replication analyses for the 61 identified miRNAs (with q < 0.05) in 3D cohort (see Supplementary Table S1). Interestingly, we achieved full replication based on the strength of the associations (one tailed p-value < 0.05) as well as the direction and magnitude of the fold changes between Gen3G and 3D cohorts for 22 (36%) miRNAs (Table 3). Moreover, 8 (80%) out of the 10 most strongly associated miRNAs in Gen3G cohort were replicated in 3D (p < 0.05). Sixteen (16) additional miRNAs (total 38) had a similar fold change (direction and magnitude) between the Gen3G and 3D cohort without reaching our a priori selected statistical significance threshold (Supplementary Table S1). These miRNAs were included in the pathway analysis below.

3.4. Metabolic Pathway Analysis of miRNAs Associated with Maternal BMI at the First Trimester of Pregnancy

To assess the potential biological role of the 38 miRNAs associated with maternal BMI at the first trimester of pregnancy, we did a metabolic pathway analysis using miRPath v3 [32]. Seven KEGG pathways from the union pathway analysis were enriched with targets of these miRNAs. Figure 2 shows these targeted metabolic pathways and their FDR-adjusted q-value. Interestingly, the top 2 pathways were related to fatty acid metabolism (p = 1 × 10−325, 6 miRNAs) and fatty acids biosynthesis (p < 1 × 10−325, 4 miRNAs).

4. Discussion

Many women enter pregnancy with OW/O, which are associated with complications, such as GDM and pre-eclampsia. In the current study, we sequenced 658 plasma samples from two independent birth cohorts and identified over 20 miRNAs circulating in maternal plasma at first trimester that are associated with maternal BMI in early pregnancy. To the best of our knowledge, this is the largest study using next generation sequencing to identify miRNAs at the first trimester of pregnancy associated with maternal BMI.
Only a few studies have associated maternal OW/O in pregnancy with plasmatic miRNA profile dysregulation [15,16,18,19]. Interestingly, four of our top 10 miRNAs, all replicated in 3D, were previously associated with maternal obesity in cord blood [18]. Jing et al. reported that cord blood hsa-miR-1323, hsa-miR-516b-5p, hsa-miR-516a-5p, and hsa-miR-520a-3p levels were positively associated with maternal OW/O [18], which is contrary to what we have observed in the current study. This apparent discrepancy could be explained by the difference in sample origin (fetal whole cord blood vs. maternal plasma) and the timing during pregnancy (1st trimester vs. at delivery). These four miRNAs are all encoded in the C19MC, the expression of which is known to increase throughout pregnancy. Overall, our results suggest that the impact of OW/O might be dependent on the timing of collection and the tissue tested (maternal plasma vs. cord blood). In addition to these four miRNAs, our study also identified novel miRNAs associated with maternal BMI, of which many are also encoded by the C19MC. None of our miRNAs were reported in the other three studies [15,16,19].
One fascinating hypothesis is that the placenta secretes miRNAs into maternal blood contributing to feto-maternal communication and metabolic adaptation throughout pregnancy [32]. Indeed, more than 600 miRNAs dynamically expressed in the human placenta have been reported so far [10]. From these, 127 miRNAs are encoded into the C14MC, C19MC and miR-371-3 clusters. [10,11] Interestingly, 17 of the 22 (77%) replicated miRNAs are encoded by the C19MC. More broadly, 21 out of the 38 replicated miRNAs (55%) are from the C19MC, 1 from the C14MC, and 2 from the miR-371-3. C19MC is a large miRNA cluster specific to primates that encodes 46 miRNAs genes, 58 mature miRNAs and is maternally imprinted (paternally expressed) [10]. These miRNAs are also strongly expressed in trophoblasts, suggesting an important role in pregnancy and embryonic and fetal development [10]. In a previous study, we also identified miRNAs specific to pregnancy (pregnant vs. non-pregnant women) and varying between the 4th and the 16th week of pregnancy, using the same Gen3G microtranscriptomic dataset [12]. Accordingly, 5 more miRNAs (in addition to the 8 from the C19MC) of our replicated miRNAs were found upregulated in pregnancy in that previous study. This provides novel support of the roles of plasma miRNAs in pregnancy. Specifically, maternal BMI levels in pregnancy may influence the expression of miRNAs that are more abundant in pregnancy and may thus have a possible role in metabolic adaptation and development.
Maternal OW/O are also risk factors for GDM and PE and consequently the miRNAs we have identified could also contribute to the development of these complications. In another previous study from our group and using the same Gen3G microtranscriptomic dataset, we identified miRNAs in the first trimester of pregnancy associated with and predictive of insulin sensitivity in the second trimester of pregnancy. Among our replicated miRNAs, 15 were also associated with insulin sensitivity [17]. Similar analyses were also conducted to identify miRNAs associated and predictive of GDM still using this microtranscriptomics and Gen3G datasets. Among the replicated miRNAs, 16 were also associated with GDM [33]. Overall, 13 were associated with maternal obesity, insulin sensitivity and GDM. Moreover, higher circulating levels of serum hsa-miR-1323 has been associated with GDM [34], hsa-miR-517-5p and hsa-miR-520a-5p were reported to be downregulated in PE, which is consistent with our results with maternal BMI [16,35]. Finally, lower circulating levels of hsa-miR-526b-5p were detected in patients with metabolic syndrome – characterized by central excess adiposity—a direction of association that is concordant with our results [36]. Through their mechanism of action, the miRNAs that we found associated with maternal BMI might potentially have a role to play in the pathophysiology of how excess weight in early gestation may lead to pregnancy complications. However, further studies are needed to disentangle their specific roles and the causal relationship.
Our pathway analysis showed that the miRNAs associated with maternal BMI regulate a few metabolic pathways, of which lipid metabolism is of high interest for its biologic relevance. Overall, 10 miRNAs were linked to fatty acid biosynthesis (5 miRNAs), fatty acid metabolism (10 miRNAs), or fatty acid elongation (5 miRNAs) pathways. They were either up—(n = 5) or down—(n = 5) regulated in relation to maternal BMI. Obesity, inside and outside of pregnancy, is associated with dyslipidemia, among several other metabolic complications [37]. Early pregnancy is characterized by a lipogenic profile which is believed to favor energy storage that will be later used to meet the metabolic demand of both the mother and her growing fetus during pregnancy [38]. By targeting fatty acid biosynthesis and metabolism, our BMI-associated miRNAs could play an early role in the regulation of the lipid metabolism pathways in pregnancy or dysregulation in response to maternal OW/O. Further studies are also needed to understand the roles of plasma miRNAs in metabolic adaptation to pregnancy and their links with maternal OW/O.

Strengths and Limitations

To the best of our knowledge, this is the largest study investigating circulating miRNAs associated with maternal BMI in the first trimester of pregnancy using next-generation sequencing. Among the strengths, our study was conducted in two large, well-characterized prospective pregnancy and birth cohorts and on plasma samples collected early in pregnancy, from the 4th to 16th weeks. This allowed us to assess obesity-related miRNA dysregulation months before pregnancy complications, such as GDM and PE, develop. Finally, we applied next generation sequencing technology, which allowed us to fully profile and quantify plasma miRNAs. This technology is a robust and sensitive approach for miRNA quantification [39]. Finally, although our study was by design associative, the independent replication significantly improves the robustness of the results we are reporting.
Our study also has some limitations. First, BMI is an easy and convenient assessment of OW/O but remains a surrogate measure of fat mass. Nevertheless, BMI remains recommended by the WHO to assess OW/O. Also, our BMIs were measured at the time plasma samples were collected. Overall, BMI in the first trimester of pregnancy is considered an acceptable assessment of that before pregnancy as weight is relatively stable in the first weeks. Also, it seems clear that functional studies are needed to confirm the role of these miRNAs associated during pregnancy (metabolic adaptation, fetal development and in pregnancy complications.)

5. Conclusions

We have identified 22 plasma miRNAs associated with maternal BMI in the first trimester of pregnancy. Identified miRNAs were enriched in biological pathways related to fatty acids and lipids that are implicated in the pathophysiology of OW/O and mostly encoded by the C19MC, which is mainly expressed by the placenta. These results could provide new insights into the understanding of the effect of OW/O on the development of different complications in pregnancy such as GDM and PE.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biomedicines10071726/s1, Figure S1: Selection of participants for the 3D cohort. This diagram is showing the different inclusion and exclusion criteria of our selected participants in the 3D cohort. Abbreviations: OGTT: oral glucose tolerance test, GDM: gestational diabetes mellitus, BMI: body mass index. Table S1: List of miRNAs significantly associated with maternal BMI at the first trimester of pregnancy in the Gen3G cohort and replication results obtained in the 3D cohort for miRNAs significantly associated with maternal BMI at 1st trimester of pregnancy in the Gen3G cohort.

Author Contributions

L.B. designed the study with the contribution of V.D., M.S.S., P.-É.J., P.P., M.-F.H. and R.G.; C.L. and V.D. performed data collection with K.T.; C.L., A.-A.C. and F.W. worked on bioinformatics analysis; K.T. performed statistical analyses; K.T. wrote the manuscript with the collaboration of V.D. and L.B.; P.P., M.-F.H., R.G. and L.B. supervised all steps of the study. All authors have read and agreed to the published version of the manuscript.

Funding

K.T., C.L. and A.-A.C. were supported by a doctoral research award from Fonds de la recherche du Québec en santé (FRQS). V.D. received funding from Diabète Québec. L.B., M.S.S. and P.É.J. are research scholars from the FRQS and a member of the CR-CHUS, a FRQS-funded Research Center. This study was supported by the CIHR (Grant #IGH-155183), the Fondation de ma vie of the CIUSSS du Saguenay–Lac-St-Jean—Hôpital Universitaire de Chicoutimi and Diabète Québec. The Gen3G birth cohort recruitment was supported by operating grants from the FRQS (Grant #20697), the CIHR (Grant #MOP 115071), Diabète Québec, and the Canadian Diabetes Association (CDA; Grant #OG-3-08-2622-JA).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Ethics Committee of CIUSSS de l’Estrie-CHUS (protocol code #MP-31-2019-3059 and February 7 2019)” for studies involving humans.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The datasets used and/or analyzed during the current study are available from the corresponding author upon reasonable request.

Acknowledgments

We thank the McGill University and Génome Québec Innovation Centre (Montréal, Canada) staff for their work in libraries quality control and sequencing. We are also grateful to Calcul Québec (Montréal, Canada) and Compute Canada (Toronto, Canada) for their support in this research. Finally, we thank all the participants of the Gen3G and 3D cohorts.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Caballero, B. The Global Epidemic of Obesity: An Overview. Epidemiol. Rev. 2007, 29, 1–5. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Valsamakis, G.; Kyriazi, E.; Mouslech, Z.; Siristatidis, C.; Mastorakos, G. Effect of Maternal Obesity on Pregnancy Outcomes and Long-Term Metabolic Consequences. Hormones 2015, 14, 345–357.e3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Tauqeer, Z.; Gomez, G.; Stanford, F.C. Obesity in Women: Insights for the Clinician. J. Womens Health 2018, 27, 444–457. [Google Scholar] [CrossRef] [PubMed]
  4. Catalano, P.M.; Shankar, K. Obesity and Pregnancy: Mechanisms of Short Term and Long Term Adverse Consequences for Mother and Child. BMJ Br. Med. J. Online 2017, 356, j1. [Google Scholar] [CrossRef]
  5. Berger, H.; Melamed, N.; Murray-Davis, B.; Hasan, H.; Mawjee, K.; Barrett, J.; McDonald, S.D.; Geary, M.; Ray, J.G. Prevalence of Pre-Pregnancy Diabetes, Obesity, and Hypertension in Canada. J. Obstet. Gynaecol. Can. 2019, 41, 1579–1588.e2. [Google Scholar] [CrossRef]
  6. Overweight & Obesity Statistics | NIDDK. Available online: https://www.niddk.nih.gov/health-information/health-statistics/overweight-obesity (accessed on 21 March 2022).
  7. Catalano, P.M.; McIntyre, H.D.; Cruickshank, J.K.; McCance, D.R.; Dyer, A.R.; Metzger, B.E.; Lowe, L.P.; Trimble, E.R.; Coustan, D.R.; Hadden, D.R.; et al. The Hyperglycemia and Adverse Pregnancy Outcome Study: Associations of GDM and Obesity with Pregnancy Outcomes. Diabetes Care 2012, 35, 780–786. [Google Scholar] [CrossRef] [Green Version]
  8. Kelly, A.C.; Powell, T.L.; Jansson, T. Placental Function in Maternal Obesity. Clin. Sci. 2020, 134, 961–984. [Google Scholar] [CrossRef] [Green Version]
  9. Godfrey, K.M.; Reynolds, R.M.; Prescott, S.L.; Nyirenda, M.; Jaddoe, V.W.V.; Eriksson, J.G.; Broekman, B.F.P. Influence of Maternal Obesity on the Long-Term Health of Offspring. Lancet Diabetes Endocrinol. 2017, 5, 53–64. [Google Scholar] [CrossRef] [Green Version]
  10. Poirier, C.; Desgagné, V.; Guérin, R.; Bouchard, L. MicroRNAs in Pregnancy and Gestational Diabetes Mellitus: Emerging Role in Maternal Metabolic Regulation. Curr. Diab. Rep. 2017, 17, 35. [Google Scholar] [CrossRef]
  11. Morales-Prieto, D.M.; Ospina-Prieto, S.; Chaiwangyen, W.; Schoenleben, M.; Markert, U.R. Pregnancy-Associated MiRNA-Clusters. J. Reprod. Immunol. 2013, 97, 51–61. [Google Scholar] [CrossRef]
  12. Légaré, C.; Clément, A.-A.; Desgagné, V.; Thibeault, K.; White, F.; Guay, S.-P.; Arsenault, B.J.; Scott, M.S.; Jacques, P.-É.; Perron, P.; et al. Human Plasma Pregnancy-Associated MiRNAs and Their Temporal Variation within the First Trimester of Pregnancy. Reprod. Biol. Endocrinol. 2022, 20, 14. [Google Scholar] [CrossRef]
  13. Morales-Prieto, D.M.; Favaro, R.R.; Markert, U.R. Placental MiRNAs in Feto-Maternal Communication Mediated by Extracellular Vesicles. Placenta 2020, 102, 27–33. [Google Scholar] [CrossRef] [PubMed]
  14. Awamleh, Z.; Gloor, G.B.; Han, V.K.M. Placental MicroRNAs in Pregnancies with Early Onset Intrauterine Growth Restriction and Preeclampsia: Potential Impact on Gene Expression and Pathophysiology. BMC Med. Genom. 2019, 12, 91. [Google Scholar] [CrossRef] [Green Version]
  15. Hromadnikova, I.; Kotlabova, K.; Ondrackova, M.; Pirkova, P.; Kestlerova, A.; Novotna, V.; Hympanova, L.; Krofta, L. Expression Profile of C19MC MicroRNAs in Placental Tissue in Pregnancy-Related Complications. DNA Cell Biol. 2015, 34, 437–457. [Google Scholar] [CrossRef]
  16. Hromadnikova, I.; Dvorakova, L.; Kotlabova, K.; Krofta, L. The Prediction of Gestational Hypertension, Preeclampsia and Fetal Growth Restriction via the First Trimester Screening of Plasma Exosomal C19MC MicroRNAs. Int. J. Mol. Sci. 2019, 20, 2972. [Google Scholar] [CrossRef] [PubMed]
  17. Légaré, C.; Desgagné, V.; Poirier, C.; Thibeault, K.; White, F.; Clément, A.-A.; Scott, M.S.; Jacques, P.-É.; Perron, P.; Guérin, R.; et al. First Trimester Plasma MicroRNAs Levels Predict Matsuda Index-Estimated Insulin Sensitivity between 24th and 29th Week of Pregnancy. BMJ Open Diabetes Res. Care 2022, 10, e002703. [Google Scholar] [CrossRef]
  18. Jing, J.; Wang, Y.; Quan, Y.; Wang, Z.; Liu, Y.; Ding, Z. Maternal Obesity Alters C19MC MicroRNAs Expression Profile in Fetal Umbilical Cord Blood. Nutr. Metab. 2020, 17, 52. [Google Scholar] [CrossRef] [PubMed]
  19. Tsamou, M.; Martens, D.S.; Winckelmans, E.; Madhloum, N.; Cox, B.; Gyselaers, W.; Nawrot, T.S.; Vrijens, K. Mother’s Pre-Pregnancy BMI and Placental Candidate MiRNAs: Findings from the ENVIRONAGE Birth Cohort. Sci. Rep. 2017, 7, 5548. [Google Scholar] [CrossRef]
  20. Enquobahrie, D.A.; Wander, P.L.; Tadesse, M.G.; Qiu, C.; Holzman, C.; Williams, M.A. Maternal Pre-Pregnancy Body Mass Index and Circulating MicroRNAs in Pregnancy. Obes. Res. Clin. Pract. 2017, 11, 464–474. [Google Scholar] [CrossRef]
  21. Carreras-Badosa, G.; Bonmatí, A.; Ortega, F.-J.; Mercader, J.-M.; Guindo-Martínez, M.; Torrents, D.; Prats-Puig, A.; Martinez-Calcerrada, J.-M.; de Zegher, F.; Ibáñez, L.; et al. Dysregulation of Placental MiRNA in Maternal Obesity Is Associated with Pre- and Postnatal Growth. J. Clin. Endocrinol. Metab. 2017, 102, 2584–2594. [Google Scholar] [CrossRef] [Green Version]
  22. Guillemette, L.; Allard, C.; Lacroix, M.; Patenaude, J.; Battista, M.-C.; Doyon, M.; Moreau, J.; Ménard, J.; Bouchard, L.; Ardilouze, J.-L.; et al. Genetics of Glucose Regulation in Gestation and Growth (Gen3G): A Prospective Prebirth Cohort of Mother–Child Pairs in Sherbrooke, Canada. BMJ Open 2016, 6, e010031. [Google Scholar] [CrossRef]
  23. Krukowski, R.A.; West, D.S.; DiCarlo, M.; Shankar, K.; Cleves, M.A.; Saylors, M.E.; Andres, A. Are Early First Trimester Weights Valid Proxies for Preconception Weight? BMC Pregnancy Childbirth 2016, 16, 357. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Fraser, W.D.; Shapiro, G.D.; Audibert, F.; Dubois, L.; Pasquier, J.; Julien, P.; Bérard, A.; Muckle, G.; Trasler, J.; Tremblay, R.E.; et al. 3D Cohort Study: The Integrated Research Network in Perinatology of Quebec and Eastern Ontario. Paediatr. Perinat. Epidemiol. 2016, 30, 623–632. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Burgos, K.L.; Javaherian, A.; Bomprezzi, R.; Ghaffari, L.; Rhodes, S.; Courtright, A.; Tembe, W.; Kim, S.; Metpally, R.; Van Keuren-Jensen, K. Identification of Extracellular MiRNA in Human Cerebrospinal Fluid by Next-Generation Sequencing. RNA 2013, 19, 712–722. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Rozowsky, J.; Kitchen, R.; Park, J.J.; Galeev, T.R.; Diao, J.; Warrell, J.; Thistlethwaite, W.; Subramanian, S.L.; Milosavljevic, A.; Gerstein, M. ExceRpt: A Comprehensive Analytic Platform for Extracellular RNA Profiling. Cell Syst. 2019, 8, 352–357.e3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Kozomara, A.; Griffiths-Jones, S. MiRBase: Annotating High Confidence MicroRNAs Using Deep Sequencing Data. Nucleic Acids Res. 2014, 42, D68–D73. [Google Scholar] [CrossRef] [Green Version]
  28. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast Universal RNA-Seq Aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef]
  29. Love, M.I.; Huber, W.; Anders, S. Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [Green Version]
  30. Blighe, K.; Rana, S.; Lewis, M. EnhancedVolcano: Publication-Ready Volcano Plots with Enhanced Colouring and Labeling. Available online: https://github.com/kevinblighe/EnhancedVolcano (accessed on 16 February 2022).
  31. Vlachos, I.S.; Zagganas, K.; Paraskevopoulou, M.D.; Georgakilas, G.; Karagkouni, D.; Vergoulis, T.; Dalamagas, T.; Hatzigeorgiou, A.G. DIANA-MiRPath v3.0: Deciphering MicroRNA Function with Experimental Support. Nucleic Acids Res. 2015, 43, W460–W466. [Google Scholar] [CrossRef]
  32. Luo, S.-S.; Ishibashi, O.; Ishikawa, G.; Ishikawa, T.; Katayama, A.; Mishima, T.; Takizawa, T.; Shigihara, T.; Goto, T.; Izumi, A.; et al. Human Villous Trophoblasts Express and Secrete Placenta-Specific MicroRNAs into Maternal Circulation via Exosomes. Biol. Reprod. 2009, 81, 717–729. [Google Scholar] [CrossRef] [Green Version]
  33. Légaré, C.; Desgagné, V.; Thibeault, K.; White, F.; Clément, A.-A.; Poirier, C.; Luo, Z.-C.; Scott, M.; Jacques, P.-É.; Perron, P.; et al. First Trimester Plasma MicroRNA Levels Predict Risk of Developing Gestational Diabetes Mellitus. Front. Endocrinol. 2022, 13. [Google Scholar]
  34. Liu, L.; Zhang, J.; Liu, Y. MicroRNA-1323 Serves as a Biomarker in Gestational Diabetes Mellitus and Aggravates High Glucose-Induced Inhibition of Trophoblast Cell Viability by Suppressing TP53INP1. Exp. Ther. Med. 2021, 21, 230. [Google Scholar] [CrossRef] [PubMed]
  35. Fu, J.-Y.; Xiao, Y.-P.; Ren, C.-L.; Guo, Y.-W.; Qu, D.-H.; Zhang, J.-H.; Zhu, Y.-J. Up-Regulation of MiR-517-5p Inhibits ERK/MMP-2 Pathway: Potential Role in Preeclampsia. Eur. Rev. Med. Pharmacol. Sci. 2018, 22, 6599–6608. [Google Scholar] [CrossRef] [PubMed]
  36. Liu, G.; Lei, Y.; Luo, S.; Huang, Z.; Chen, C.; Wang, K.; Yang, P.; Huang, X. MicroRNA Expression Profile and Identification of Novel MicroRNA Biomarkers for Metabolic Syndrome. Bioengineered 2021, 12, 3864–3872. [Google Scholar] [CrossRef] [PubMed]
  37. Klop, B.; Elte, J.W.F.; Castro Cabezas, M. Dyslipidemia in Obesity: Mechanisms and Potential Targets. Nutrients 2013, 5, 1218–1240. [Google Scholar] [CrossRef] [Green Version]
  38. Herrera, E. Lipid Metabolism in Pregnancy and Its Consequences in the Fetus and Newborn. Endocrine 2002, 19, 43–55. [Google Scholar] [CrossRef]
  39. Tam, S.; de Borja, R.; Tsao, M.-S.; McPherson, J.D. Robust Global MicroRNA Expression Profiling Using Next-Generation Sequencing Technologies. Lab. Investig. J. Tech. Methods Pathol. 2014, 94, 350–358. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Plasmatic miRNAs associated with maternal BMI at 1st trimester of pregnancy. This volcano plot shows the miRNAs associated with maternal BMI. Each point represents a single miRNA. The vertical dotted line represents a log2 fold change of 0, and the horizontal dotted line represents the FDR-adjusted q-value threshold of < 0.05. The model was adjusted for sequencing lane and run, maternal and gestational age at 1st trimester, and parity. The fold change represents the change in miRNA abundance for an increase of one unit of maternal BMI in 1st trimester of pregnancy. Abbreviations: FDR: false discovery rate; NS: not significant.
Figure 1. Plasmatic miRNAs associated with maternal BMI at 1st trimester of pregnancy. This volcano plot shows the miRNAs associated with maternal BMI. Each point represents a single miRNA. The vertical dotted line represents a log2 fold change of 0, and the horizontal dotted line represents the FDR-adjusted q-value threshold of < 0.05. The model was adjusted for sequencing lane and run, maternal and gestational age at 1st trimester, and parity. The fold change represents the change in miRNA abundance for an increase of one unit of maternal BMI in 1st trimester of pregnancy. Abbreviations: FDR: false discovery rate; NS: not significant.
Biomedicines 10 01726 g001
Figure 2. KEGG pathways targeted by miRNAs associated with maternal BMI at the first trimester of pregnancy. This bar graph is showing the KEGG pathways enriched with targets of miRNAs associated with maternal BMI. Each bar represents one pathway, ranked according to its level of significance (FDR adjusted q-value). The number of miRNAs involved in the pathway is shown directly in the bar. Abbreviations: ECM: extracellular matrix; FDR: false discovery rate; KEGG: Kyoto Encyclopedia of Genes and Genomes.
Figure 2. KEGG pathways targeted by miRNAs associated with maternal BMI at the first trimester of pregnancy. This bar graph is showing the KEGG pathways enriched with targets of miRNAs associated with maternal BMI. Each bar represents one pathway, ranked according to its level of significance (FDR adjusted q-value). The number of miRNAs involved in the pathway is shown directly in the bar. Abbreviations: ECM: extracellular matrix; FDR: false discovery rate; KEGG: Kyoto Encyclopedia of Genes and Genomes.
Biomedicines 10 01726 g002
Table 1. Maternal characteristics at the first trimester of pregnancy in Gen3G (n = 435) and 3D (n = 223) cohorts.
Table 1. Maternal characteristics at the first trimester of pregnancy in Gen3G (n = 435) and 3D (n = 223) cohorts.
CharacteristicsGen3G Cohort3D CohortComparison
between Cohorts
(p-Value) *
Mean ± SDRangeMean ± SDRange
BMI (kg/m2)25.95 ± 5.9816.10–54.1025.39 ± 5.9116.80–48.500.22
Age (years)28.48 ± 4.2618–4730.60 ± 3.9020–42<0.001
Gestational age (weeks)9.63 ± 2.264.09–16.3011.89 ± 1.527.71–16.43<0.001
Parity(% nulliparous)0.69 ± 0.91 (53.56)0–60.47 ± 0.66 (60.99)0–30.01
* Comparisons between Gen3G and 3D cohorts were performed using a Mann-Whitney U test. Abbreviations: 3D: Design, Develop, Discover birth cohort; BMI: body mass index; Gen3G: Genetics of Glucose regulation in Gestation and Growth birth cohort; SD: standard deviation.
Table 2. Top 10 miRNAs most significantly associated with maternal BMI at first trimester of pregnancy in the Gen3G cohort.
Table 2. Top 10 miRNAs most significantly associated with maternal BMI at first trimester of pregnancy in the Gen3G cohort.
miRNADESeq2 Normalized Read Count (Mean ± SD)Fold Change *Unadjusted p-ValueFDR-Adjusted q-Value
hsa-miR-1323 a146.39 ± 230.600.9572.79 × 10109.60 × 10−8
hsa-miR-516b-5p a101.50 ± 150.350.9581.79 × 10109.60 × 10−8
hsa-miR-371a-5p b11.02 ± 17.090.9412.13 × 10109.60 × 10−8
hsa-miR-525-5p a9.07 ± 15.310.9521.63 × 10−94.21 × 10−7
hsa-miR-516a-5p a31.25 ± 53.810.9594.55 × 10−97.83 × 10−7
hsa-miR-524-5p a8.64 ± 14.480.9523.90 × 10−97.83 × 10−7
hsa-miR-518e-5p a43.96 ± 62.520.9628.92 × 10−91.22 × 10−6
hsa-miR-520a-3p a86.75 ± 132.370.9609.41 × 10−91.22 × 10−6
hsa-miR-518e-3p a7.99 ± 14.220.9562.02 × 10−82.09 × 10−6
hsa-miR-520d-5p a5.92 ± 10.350.9551.98 × 10−82.09 × 10−6
Model adjusted for sequencing run and lane, gestational age, maternal age at the first trimester of pregnancy, and parity. * Fold changes represent the change in miRNA abundance for each increase of one unit of maternal BMI at 1st trimester of pregnancy. a miRNAs from C19MC. b miRNAs from miR-371-3 miRNAs cluster. Abbreviations: FDR: false discovery rate; SD: standard deviation.
Table 3. miRNAs significantly associated with maternal BMI at 1st trimester of pregnancy in Gen3G cohort and fully replicated in 3D cohort.
Table 3. miRNAs significantly associated with maternal BMI at 1st trimester of pregnancy in Gen3G cohort and fully replicated in 3D cohort.
Gen3G3D
miRNADESeq2 Normalized Read Count (Mean ± SD)Fold Change *Unadjusted p-ValueFDR-Adjusted q-ValueDESeq2 Normalized Read Count (Mean ± SD)Fold Change *Unadjusted p-Value
hsa-miR-1323 a146.39 ± 230.600.9572.79 × 10109.60 × 10−8676.45 ± 581.880.9662.56 × 10−5
hsa-miR-516b-5p a101.50 ± 150.350.9581.79 × 10109.60 × 10−8327.56 ± 247.260.9662.92 × 10−5
hsa-miR-525-5p a9.07 ± 15.310.9521.63 × 10−94.21 × 10−779.89 ± 71.160.9820.04471
hsa-miR-516a-5p a31.25 ± 53.810.9594.55 × 10−97.83 × 10−7150.92 ± 135.760.9740.00141
hsa-miR-518e-5p a43.96 ± 62.520.9628.92 × 10−91.22 × 10−699.07 ± 88.520.9810.02118
hsa-miR-520a-3p a86.75 ± 132.370.9609.41 × 10−91.22 × 10−6220.91 ± 243.900.9820.04009
hsa-miR-518e-3p a7.99 ± 14.220.9562.02 × 10−82.09 × 10−629.60 ± 25.750.9710.00189
hsa-miR-512-3p a287.04 ± 575.480.9633.63 × 10−83.41 × 10−6771.74 ± 982.110.9820.02747
hsa-miR-1283 a57.10 ± 87.820.9624.17 × 10−83.59 × 10−6203.52 ± 173.840.9770.00635
hsa-miR-517a-3p a20.19 ± 37.780.9605.17 × 10−84.11 × 10−676.56 ± 73.880.9800.01599
hsa-miR-526b-5p a17.12 ± 25.900.9622.19 × 10−71.62 × 10−539.62 ± 34.830.9770.01056
hsa-miR-517-5p a16.90 ± 29.750.9629.56 × 10−76.17 × 10−570.67 ± 67.450.9810.02191
hsa-miR-519c-3p a11.08 ± 20.080.9655.11 × 10−60.0002636.44 ± 32.830.9660.00034
hsa-miR-519d-5p a6.12 ± 9.510.9611.53 × 10−50.0006310.75 ± 11.770.9690.01205
hsa-miR-515-5p a10.50 ± 19.120.9671.50 × 10−50.0006326.70 ± 26.080.9680.00076
hsa-miR-27b-3p29,699.76 ± 20,347.781.0084.98 × 10−50.0017255,214.27 ± 41,098.161.0110.00110
hsa-miR-885-5p13.50 ± 17.921.0364.98 × 10−50.0017240.63 ± 78.981.0340.00588
hsa-miR-520a-5p a3.06 ± 5.660.9630.000320.0083027.63 ± 27.090.9650.00136
hsa-miR-3751938.28 ± 3939.320.9750.000400.009862900.08 ± 2878.580.9820.03889
hsa-miR-520d-3p a9.80 ± 15.840.9720.000910.0187710.77 ± 13.370.9660.00638
hsa-miR-5921.64 ± 3.001.0420.001620.031032.43 ± 4.891.0580.00400
hsa-miR-21-5p85,256.17 ± 61,108.551.0050.002240.0381182,095.57 ± 63,195.511.0080.02738
Models adjusted for sequencing run and lane, gestational age, maternal age at the first trimester of pregnancy, and parity. * Fold changes represent the change in miRNA abundance for each increase of one unit of maternal BMI at 1st trimester of pregnancy. a miRNAs from C19MC. Abbreviations: 3D: Design, Develop, Discover birth cohort; FDR: false discovery rate; Gen3G: Genetics of Glucose regulation in Gestation and Growth birth cohort; SD: standard deviation.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Thibeault, K.; Légaré, C.; Desgagné, V.; White, F.; Clément, A.-A.; Scott, M.S.; Jacques, P.-É.; Guérin, R.; Perron, P.; Hivert, M.-F.; et al. Maternal Body Mass Index Is Associated with Profile Variation in Circulating MicroRNAs at First Trimester of Pregnancy. Biomedicines 2022, 10, 1726. https://doi.org/10.3390/biomedicines10071726

AMA Style

Thibeault K, Légaré C, Desgagné V, White F, Clément A-A, Scott MS, Jacques P-É, Guérin R, Perron P, Hivert M-F, et al. Maternal Body Mass Index Is Associated with Profile Variation in Circulating MicroRNAs at First Trimester of Pregnancy. Biomedicines. 2022; 10(7):1726. https://doi.org/10.3390/biomedicines10071726

Chicago/Turabian Style

Thibeault, Kathrine, Cécilia Légaré, Véronique Desgagné, Frédérique White, Andrée-Anne Clément, Michelle S. Scott, Pierre-Étienne Jacques, Renée Guérin, Patrice Perron, Marie-France Hivert, and et al. 2022. "Maternal Body Mass Index Is Associated with Profile Variation in Circulating MicroRNAs at First Trimester of Pregnancy" Biomedicines 10, no. 7: 1726. https://doi.org/10.3390/biomedicines10071726

APA Style

Thibeault, K., Légaré, C., Desgagné, V., White, F., Clément, A. -A., Scott, M. S., Jacques, P. -É., Guérin, R., Perron, P., Hivert, M. -F., & Bouchard, L. (2022). Maternal Body Mass Index Is Associated with Profile Variation in Circulating MicroRNAs at First Trimester of Pregnancy. Biomedicines, 10(7), 1726. https://doi.org/10.3390/biomedicines10071726

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop