Distal mediator-enriched, placental transcriptome-wide analyses illustrate the Developmental Origins of Health and Disease

As the master regulator of the intrauterine environment, the placenta is core to the Developmental Origins of Health and Disease (DOHaD) but is understudied in relation to tissue-specific gene and trait regulation. We performed distal mediator-enriched transcriptome-wide association studies (TWAS) for 40 health traits across 5 physiological categories, using gene expression models trained with multi-omic data from the Extremely Low Gestational Age Newborn Study (N = 272). At P < 2.5 x 10-6, we detected 248 gene-trait associations (GTAs) across 176 genes, mostly for metabolic and neonatal traits and enriched for cell growth and immunological pathways. Of these GTAs, 89 showed significant mediation through genetic variants distal to the gene, identifying potential targets for functional validation. Functional validation of a mediator gene (EPS15) in human placenta-derived JEG-3 trophoblasts resulted in increased expression of its predicted targets, SPATA13 and FAM214A, both associated with the trait of waist-hip ratio in TWAS. These results illustrate the profound health impacts of placental genetic and genomic regulation in developmental programming across the life course.


MAIN
The placenta serves as the master regulator of the intrauterine environment via nutrient transfer, 2 metabolism, gas exchange, neuroendocrine signaling, growth hormone production, and immunologic 3 control 1-5 . Due to strong influences on postnatal health, the placenta is central to the Developmental 4 Origins of Health and Disease (DOHaD) hypothesis, which purports that the in utero experience has 5 lifelong impacts on child health by altering developmental programming and influencing risk of common, 6 noncommunicable health conditions 6 . For example, placental biology has been linked to neuropsychiatric, 7 developmental, and metabolic diseases or health traits (collectively referred to as traits) that manifest 8 throughout the life course, either early-or later-in-life (Figure 1) [7][8][9][10] . Despite its long-lasting influences on 9 health, the placenta has not been well-studied in large consortia studies of multi-tissue gene 10 regulation 11, 12 . Studying regulatory mechanisms in the placenta underlying biological processes in 11 developmental programming will provide novel insight into health and disease etiology. 12 13 The complex interplay between genetics and placental transcriptomics and epigenomics has strong 14 effects on gene expression that may explain variation in gene-trait associations (GTAs). Quantitative trait 15 loci (QTL) analyses have identified a strong influence of cis-genetic variants on both placental gene 16 expression and DNA methylation 13,14 . Furthermore, there is growing evidence that the placental 17 epigenome influences gene regulation, often distally (more than 1-3 Megabases away in the genome) 15 , 18 and that placental DNA methylation and microRNA (miRNA) expression are associated with health traits 19 in children [16][17][18] . Dysfunction of transcription factor regulation in the placenta has also shown profound 20 effects on childhood traits [19][20][21][22] . Although combining genetics, transcriptomics, and epigenomics lends 21 insight into the influence of placental genomics on complex traits 23 , genome-wide screens for GTAs that 22 integrate different molecular profiles and generate functional hypotheses require more sophisticated 23 computational methods. To this end, advances in transcriptome-wide association studies (TWAS) have allowed for integration of 26 genome-wide association studies (GWAS) and eQTL datasets to boost power in identifying GTAs, 27 specific to a relevant tissue [24][25][26] . However, traditional methods for TWAS largely overlook genetic variants 28 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted April 14, 2021. ; https://doi.org/10.1101/2021.04.12.21255170 doi: medRxiv preprint distal to genes of interest, ostensibly mediated through regulatory biomarkers (e.g., transcription factors, 1 miRNAs, and DNA methylation sites) 27,28 . Not only may these distal biomarkers explain a significant 2 portion of both gene expression heritability and trait heritability on the tissue-specific expression level 29-32 , 3 they may also influence tissue-specific trait associations for individual genes. Due to the strong interplay 4 of regulatory elements in placental gene regulation, we sought to systematically characterize portions of 5 gene expression that are influenced by these distal regulatory elements. 6 7 Here, we investigate three broad questions: (1) which genes show associations between their placental 8 genetically-regulated expression (GReX) and various traits across the life course, (2) which traits along 9 the life course can be explained by placental GReX, in aggregate, and (3) which transcription factors, 10 miRNAs, or CpG sites potentially regulate trait-associated genes in the placenta (Figure 1). We 11 leveraged gene expression, CpG methylation, and miRNA expression data from fetal-side placenta tissue 12 from the Extremely Low Gestational Age Newborn (ELGAN) Cohort Study 33 . We trained predictive models 13 of gene expression, enriched for distal SNPs using MOSTWAS, a recent TWAS extension that integrates 14 multi-omic data 34 . Re-analyzing 40 GWAS of European-ancestry subjects from large consortia 35-39 , we 15 performed a series of TWAS for non-communicable health traits and disorders that may be influenced by 16 the placenta to identify GTAs and functional hypotheses for regulation (Figure 2). To our knowledge, this 17 is the first distal mediator-enriched TWAS of health traits that integrates placental multi-omics. Results 18 from our analysis can be explored at our Shiny R app, the ELGAN DOHaD Atlas: https://elgan-19 twas.shinyapps.io/dohad/. 20 23 From large consortia 35-39 , we curated GWAS summary statistics from subjects of European ancestry for 24 40 complex, non-communicable traits and disorders across five health categories to systematically 25 identify potential links to genetically-regulated placental expression ( CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted April 14, 2021. ; https://doi.org/10.1101/2021.04.12.21255170 doi: medRxiv preprint (Supplemental Table 1). These five categories of traits have been linked previously to placental biology 1 and morphology 7-10 . 2 3 To assess the percent variance explained by genetics in each trait and the genetic associations shared 4 between traits, we estimated the SNP heritability (ℎ 2 ) and genetic correlation ( ) of these traits, 5 respectively (Supplemental Figure S1 and S2). Of the 40 traits, 37 showed significantly positive SNP 6 heritability and 18 with ĥ 2 > 0.10 (Supplemental Figure S1, Supplemental Table S1), with the largest 7 heritability for childhood BMI (ĥ 2 = 0.69, = 0.064). As expected, we observed strong, statistically 8 significant genetic correlations between traits of similar categories (i.e., between neuropsychiatric traits or 9 between metabolic traits) (Supplemental Figure S2; Supplemental Gene expression prediction models 18 To train predictive models of placental expression, the first step of our TWAS (Figure 2A), we leveraged 19 MOSTWAS 34 , a recent extension that includes distal variants in transcriptomic prediction. As large 20 proportions of total heritable gene expression are explained by distal-eQTLs local to regulatory 21 hotspots 27,28,30 , MOSTWAS uses data-driven approaches to either identify mediating regulatory 22 biomarkers or distal-eQTLs mediated through local regulatory biomarkers to increase predictive power for 23 gene expression and power to detect GTAs (Supplemental Figure S3) 34  Using genotypes (from umbilical cord blood) 40 , mRNA expression, CpG methylation, and miRNA 1 expression data (from fetal-side placenta) 23 from the ELGAN Study 33 for 272 infants born pre-term, we 2 built genetic models to predict RNA expression levels for genes in the fetal placenta (demographic 3 summary in Supplemental Table S3). Out of a total of 12,020 genes expressed across all samples in 4 ELGAN, we successfully built significant models for 2,994 genes, such that SNP-based expression 5 heritability is significantly positive (nominal < 0.05) and five-fold cross-validation (CV) adjusted 2 ≥ 6 0.01 ( Figure 3A [Step 1]); only these 2,994 models are used in subsequent TWAS steps. Mean SNP 7 heritability for these genes was 0.39 (25% quantile = 0.253, 75% quantile = 0.511), and mean CV 2 was  Table S4). 12 13 Placental transcriptome-wide association studies 14 Overall associations and permutation tests 15 We integrated GWAS summary statistics for 40 traits from European-ancestry subjects with placental 16 gene expression using our predictive models. Using the weighted burden test 25,43 , we detected 932 GTAs 17 (spanning 686 unique genes) at < 2.5 × 10 −6 (corresponding to | | > 4.56), a transcriptome-wide 18 significance threshold consistent with previous TWAS 25,31 ( Figure 3A [Step 2], Supplemental Data). As 19 many of these loci carry significant signal because of strong trait-associated GWAS architecture, we 20 employed Gusev et al's permutation test to assess how much signal is added by the SNP-expression 21 weights and confidently conclude that integration of expression data significantly refines association with 22 the trait 25 . At FDR-adjusted permutation < 0.05 and spanning 176 unique genes, we detected 248 such 23 GTAs, of which 11 were found in autoimmune/autoreactive disorders, 136 in body size/metabolic traits,   176397). We only detected a significant association between PSG8 and fetal birthweight ( = −7.77). 10 Similarly, of the 6 childhood BMI-associated genes identified by Peng et al, only 1 had a significant model 11 in ELGAN and showed no association with the trait; there were no overlaps with childhood obesity-12 associated genes from Peng et al 10 . We hypothesize that minimal overlap with susceptibility genes 13 identified by Peng et al is due to differing eQTL architectures in the datasets and different inclusion 14 criteria for significant gene expression models 10,34,44,45 . 15 16 We conducted over-representation analysis for biological process, molecular function, and PANTHER 17 gene pathway ontologies for TWAS-detected susceptibility genes (Supplemental Figure S10, 18 Supplemental Table S6) 46 . Overall, considering all 176 TWAS-identified genes, we observed 19 enrichments for nucleic acid binding and immune or cell growth signaling pathways (e.g., B-cell/T-cell 20 activation and EGF receptor, interleukin, PDGF, and Ras signaling pathways). By trait, we found related 21 pathways (sphingolipid biosynthesis, cell motility, etc) for TWAS genes for metabolic and morphological 22 traits (e.g., BMI and childhood BMI); for most traits, we were underpowered to detect ontology 23 enrichments (Supplemental Table S6). We also assessed the overlap of TWAS genes with GWAS 24 signals. A total of 112 TWAS genes did not overlap with GWAS loci ( < 5 × 10 −8 ) within a 500 kilobase 25 interval around any SNPs (local and distal) included in predictive models ( Table 2). CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted April 14, 2021. ; https://doi.org/10.1101/2021.04.12.21255170 doi: medRxiv preprint To assess how, on the whole, genetically regulated placental expression explains trait variance, we 1 computed trait heritability on the placental expression level (ℎ 2 ) using all examined and all TWAS-2 prioritized susceptibility genes using RHOGE, an linkage disequilibrium (LD) score regression 3 approach 25,47 . Overall, we found 3/14 neonatal traits (childhood BMI, total puberty growth, and pubertal 4 growth start) with significant ĥ 2 > 0 (FDR-adjusted < 0.05 for jack-knife test of significance) 31 ; none of 5 the 26 traits outside the neonatal category were appreciably explained by placental GReX. Figure 3D   6 shows that mean ĥ 2 is higher in neonatal traits than other groups. A comparison of the number of 7 GWAS-significant SNPs and TWAS-significant genes also shows that neonatal/childhood traits are 8 enriched for placental TWAS associations, even though significant genome-wide GWAS architecture 9 cannot be inferred for these traits (Supplemental Figure S11). Not only do these results highlight the 10 power advantage of properly aligned tissue-specific TWAS compared to GWAS, they suggest that 11 placental GReX affects neonatal traits more profoundly, as a significantly larger proportion of neonatal 12 traits showed significant heritability on the placental GReX level than later-in-life traits. 13 14 Similarly, using RHOGE 31 , we assessed genetic correlations ( ) between traits at the level of placental 15 GReX (Supplemental Figure S12). We found several known correlations, such as between cholesterol 16 and triglycerides (̂= 0.99, = 1.44 × 10 −118 ) and childhood BMI and adult BMI (̂= 0.55, = 17 3.67 × 10 −8 ). Interestingly, we found correlations between traits across categories: IQ and diastolic blood 18 pressure (̂= −0.55, = 2.44 × 10 −5 ) and age of asthma diagnosis and glucose levels (̂= 0.86, = 19 3.05 × 10 −6 ); these traits have been linked in morphological analyses of the placenta, but our results 20 suggest possible gene regulatory contributions 48 . Overall, these correlations may suggest shared genetic 21 pathways for these pairs of traits or for etiologic antecedents of these traits; these shared pathways could 22 be either at the susceptibility genes or through shared distal loci, mediated by transcription factors, 23 miRNAs, or CpG methylation sites. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted April 14, 2021. ; https://doi.org/10.1101/2021.04.12.21255170 doi: medRxiv preprint TWAS-significant genes with overlapping genetic loci, FOCUS estimates posterior inclusion probabilities 1 (PIP) in a credible set of genes that explains the association signal at the locus. We found 8 such 2 overlaps and estimated a 90% credible set of genes explaining the signal for each locus (Supplemental 3   Table S8). For example, we identified 3 genes associated with triglycerides at the 12q24. 13  We also noticed that ERP29 and RPL6 (Ribosomal Protein L6, OMIM: 603703) were identified in GTAs 10 with multiple traits, leading us to examine potential horizontally pleiotropic genes. Of the 176 TWAS-11 prioritized genes, we identified 50 genes associated with multiple traits, many of which are genetically 12 correlated ( Table 3). Nine genes showed more than 3 GTAs across different categories. For example, CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) intrauterine growth restriction due to involvement with endothelial vascularization 55 , potentially suggesting 1 that CMTM4 has a more direct effect in utero, which mediates its associations with body fat percentage 2 and hypertension. 3 4 We further studied the 9 genes with 3 or more distinct GTAs across different categories ( Figure 4A). 5 Using UK Biobank 35 GWAS summary statistics, we conducted TWAS for a variety of traits across 8 6 groups, defined generally around ICD code blocks ( Figure 4A, Supplemental Figure S13); here, we 7 grouped metabolic and cardiovascular traits into one category for ease of analysis. At FDR-adjusted < We next examined whether placental GReX of these 9 genes correlate with fundamental traits at birth. 25 We imputed expression into individual-level ELGAN genotypes ( = 729) (Online Methods). Controlling 26 for race, sex, gestational duration, inflammation of the chorion, and maternal age, as described before 27 and in Online Methods 23 , we tested for associations for 6 representative traits measured at birth or at 24 28 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

11
Investigating mediators of distal SNP to gene relationships 12 An advantage of MOSTWAS's methodology is in functional hypothesis generation by identifying potential 13 mediators that affect TWAS-identified genes. Using the distal-SNPs added-last test from MOSTWAS 34 , 14 we interrogated distal loci incorporated into expression models for trait associations, beyond the  Table S10), namely catabolic and metabolic processes, response 5 to lipids, and multiple nucleic acid-binding processes 46 .  Table S11). We found a significant association between predicted EPS15 and FAM214A expressions 14 (effect size -0.24, FDR-adjusted = 0.019). In addition, we detected a significant association between 15 predicted NFKBIA (NF-Kappa-B Inhibitor Alpha, OMIM: 164008) and HNRNPU (Heterogeneous Nuclear 16 Ribonucleoprotein U, OMIM: 602869) (effect size -0.26, FDR-adjusted = 1.9 × 10 −4 ). We also 17 considered an Egger regression-based Mendelian randomization framework 66 in RICHS to estimate the 18 causal effects of TFs on the associated TWAS genes (Methods) using, as instrumental variables, cis-19 SNPs correlated to the TF and uncorrelated with the TWAS genes. We estimated significant causal 20 effects for two TF-TWAS gene pairs ( Figure 4C, Supplemental We also examined CpG methylation sites MOSTWAS marked as potential mediators for expression of 24 TWAS genes for overlap with cis-regulatory elements in the placenta from the ENCODE Project Phase 25 II 12 , identifying 34 CpG sites (mediating 29 distinct TWAS genes) that fall in cis-regulatory regions 26 (Supplemental Table S13). Interestingly, one CpG site mediating (cg15733049, Chromosome 27 1:2334974) FAM214A is found in low-DNase activity sites in placenta samples taken at various 28 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted April 14, 2021. ; https://doi.org/10.1101/2021.04.12.21255170 doi: medRxiv preprint timepoints; additionally, cg15733049 is local to EPS15, the transcription factor predicted to mediate 1 genetic regulation of FAM214A. Furthermore, expression of LARS2, a TWAS gene for BMI, is mediated 2 by cg04097236 (found within ELOVL2), a CpG site found in low DNase or high H3K27 activity regions; 3 LARS2 houses multiple GWAS risk SNPs for type 2 diabetes 67 and has shown BMI TWAS associations in 4 other tissues 25,31 . Results from these external datasets add more evidence that these mediators play a 5 role in gene regulation of these TWAS-identified genes. 6 7 In-vitro assays of transcription factor activity 8 Based on our computational results, we experimentally investigated whether the inverse relationship 9 between TF EPS15 and its two prioritized target TWAS genes, SPATA13 and FAM214A, is supported in The placenta has historically been understudied in large multi-tissue consortia efforts that study tissue-27 specific regulatory mechanisms 11,12 . To address this gap, we systematically categorized placental gene-28 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted April 14, 2021. ; https://doi.org/10.1101/2021.04.12.21255170 doi: medRxiv preprint trait associations relevant to the DOHaD hypothesis using distal mediator-enriched TWAS and deployed 1 these results at the ELGAN DOHaD Atlas (https://elgan-twas.shinyapps.io/dohad/). By integrating multi-2 omic data from the ELGAN Study 33 with 40 GWAS, we detected 176 unique genes (enriched for cell 3 growth and immune pathways) with transcriptome-wide significant associations, with the majority of GTAs 4 linked to metabolic and neonatal/childhood traits. Many of these TWAS-identified genes, especially those 5 with neonatal GTAs, showed multiple GTAs across trait categories (9 genes with 3 or more GTAs). We 6 examined phenome-wide GTAs for these 9 genes in UKBB and found enrichments for traits affecting in 7 immune and circulatory system (e.g., immune cell, erythrocyte, and platelet counts). We followed up with 8 selected at-birth traits in ELGAN and found associations with neonatal body size and infant cognitive 9 development. Furthermore, we could only estimate significantly positive placental GReX-mediated 10 heritability for neonatal traits but not for later-in-life traits. These results suggest that placental expression, 11 mediated by fetal genetics, is most likely to have large effects on early life traits, but these effects may 12 carry over later-in-life or as etiologic antecedents for complex traits. We conclude with limitations of this study and future directions. First, although TWAS is unlikely to be 2 subject to reverse-causality (trait cannot affect expression, independent of genetics), instances of 3 horizontal SNP pleiotropy, where SNPs influence the trait and expression independently, were not 4 examined here. Second, the ELGAN Study gathered molecular data from infants born extremely pre-5 term. If unmeasured confounders affect both prematurity and a trait of interest, GTAs could be subject to 6 backdoor collider confounding 74,75 . However, significant TWAS genes did not show associations for 7 gestational duration, suggesting minimal bias from this collider effect. An interesting future endeavor 8 could include negative controls to account for unmeasured confounders in predictive models 76 to allow for 9 more generalizability of predictive models. Fourth, though we did scan neonatal traits in ELGAN using 10 individual-level genotypes, the sample size is small; larger GWAS with longitudinal traits could allow for 11 rigorous Mendelian randomization studies 77 that investigate relationships between traits across the life 12 course, in the context of placental regulation. Lastly, due to small sample sizes of other ancestry groups 13 in ELGAN, we could only credibly impute expression into samples from European ancestry and our 14 TWAS only considers GWAS in populations of European ancestry. We emphasize acquisition of larger 15 genetic and genomic datasets from understudied and underserved populations, especially related to 16 early-in-life traits.

25
Data acquisition and quality control 26 Genotypes and multi-omic (mRNA, miRNA, and CpG methylation) data were collected from umbilical cord 27 blood and the fetal side of the placenta of subjects enrolled in the ELGAN 33 Study, as described in 28 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) function was used from the sva package to adjust for batch effects from sample plate and cell-type 7 heterogeneity 84 . mRNA and miRNA were aligned to the GENCODE Release 31 reference transcriptome 8 and quantified using Salmon 85 and the HTG EdgeSeq System 86 . We upper-quartile normalized 9 distributional differences between lanes 87 and used RUVSeq and limma to estimate and remove 10 unwanted variation 88,89 . Overall, we considered 846,233 CpG sites, 12,020 genes, and 1,898 miRNAs. 11 We downloaded quality-controlled genotypes and obtained normalized RNA-seq data for RICHS data for 12 validation of gene expression models 42 . Summary statistics were downloaded from the following  Table 1). Genomic coordinates were transformed to the hg38 reference genome using 16 liftOver 90,91 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted April 14, 2021. ; https://doi.org/10.1101/2021.04.12.21255170 doi: medRxiv preprint Gene expression models 1 We used MOSTWAS to train predictive models of gene expression from germline genetics, including 2 distal variants that were either close to associated mediators (transcription factors, miRNAs, CpG sites) or 3 had large indirect effects on gene expression 34 (Supplemental Figure S1, Supplemental Methods). 4 Briefly, MOSTWAS contains two methods of predicting expression: (1) mediator-enriched TWAS 5 (MeTWAS) and (2) distal-eQTL prioritization via mediation analysis. For MeTWAS, we first identified 6 mediators strongly associated with genes through correlation analyses between all genes of interest and 7 a set of distal mediators (FDR-adjusted < 0.05). We then trained local predictive models (using SNPs  19 The association between predicted expression and traits was assessed in GWAS summary statistics 20 using the weighted burden test and 1000 Genomes Project CEU population as an LD reference 25,34,43,94 21 with Bonferroni-corrected significance threshold of < 2.5 × 10 −6 . We only consider GWAS of subjects 22 from European ancestry, as ELGAN data does not have a large enough sample size of non-Europeans to 23 accurately map distal-eQTLs. In individual level data from ELGAN, we multiplied the genotype matrix by 24 the SNP-gene weights to construct imputed expression in ELGAN; for samples used in model training, we 25 used the cross-validated predicted expression. We tested the significance of expression-trait associations 26 conditional on SNP-trait effects at a locus using the permutation test from Gusev et al 25 . We also tested 27 the trait association at distal variants using the added last test from MOSTWAS 34 . Briefly, we computed 28 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted April 14, 2021. ; https://doi.org/10.1101/2021.04.12.21255170 doi: medRxiv preprint the weighted Z-score at the distal loci, conditional on the weighted Z-score at the local locus and test this 1 using the same null distribution assumptions as in the weighted burden test from Gusev et al 25  sample correlations and inferred ℎ 2 using ordinary least squares. We employed RHOGE 31 to estimate 12 and test for significant genetic correlations between traits at the predicted expression level (details in 13 Supplemental Methods).
14 15 16 For 9 genes with 3 or more associations across covariates. These covariates have been previously used in placental genomic studies of neonatal 28 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. In-vitro functional assays 13 Cell culture and treatment CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted April 14, 2021. ; https://doi.org/10.1101/2021.04.12.21255170 doi: medRxiv preprint manufacturer's protocol. RNA was quantified using a NanoDrop 1000 spectrophotometer (Thermo 1 Scientific, Waltham, MA). RNA was then converted to cDNA, the next step toward analyzing gene 2 expression. Next, mRNA expression was measured for EPS15, SPATA13, and FAM214A using real-time  Post-hoc pairwise t-tests (3 degrees of freedom for biological and technical duplicate) were utilized to 10 investigate direct comparisons within sample groups. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted April 14, 2021. ; https://doi.org/10.1101/2021.04.12.21255170 doi: medRxiv preprint CJM, RCF, and HPS collected the data. AB, TMO, RH, RCF, and HPS interpreted results. AB, RCF, and 1 HPS wrote the paper. All authors read and edited the paper. 28 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted April 14, 2021.

22.
Iglesias 28 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. 28 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. 28 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

21
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted April 14, 2021. ; https://doi.org/10.1101/2021.04.12.21255170 doi: medRxiv preprint Figure 1: Overview of the placenta and the DOHaD Hypothesis. The placenta facilitates many important functions in utero, including nutrient transfer, metabolism, gas exchange, neuroendocrine signaling, growth hormone production, and immunologic control. As such, it is known as a master regulator of the intrauterine environment and is core to the Developmental Origins of Health and Disease (DOHaD) hypothesis. Placental genomic regulation is purported to be influenced by both genetic and environmental factors and affects placental developmental programming. In turn, this programming has been shown to have profound impacts on a variety of disorders and traits, both early-and later-in-life.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted April 14, 2021. ; https://doi.org/10.1101/2021.04.12.21255170 doi: medRxiv preprint . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review)  Kernel density plots of in-(through cross-validation in ELGAN, red) and out-sample (external validation in RICHS, blue) McNemar's adjusted 2 between predicted and observed expression. Dotted and solid lines represent the mean and median of the respective distribution, respectively. (C) Bar graph of numbers of TWAS GTAs at overall TWAS < 2.5 × 10 −6 and permutation FDR-adjusted < 0.05 (X-axis) across traits (Y-axis). The total number of GTAs per trait are labeled, colored by the category of each trait. The bar is broken down by numbers of GTAs with (orange) and without (green) significant distal expression-mediated associations, as indicated by FDR-adjusted < 0.05 for the distal-SNPs added-last test. (D) Box plot of expression-mediated trait heritability (Y-axis) by category (X-axis), with labels if ĥ 2 is significantly greater than 0 using jack-knife test of significance.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted April 14, 2021. ; https://doi.org/10.1101/2021.04.12.21255170 doi: medRxiv preprint Figure 4: Computational follow-up analyses of TWAS-prioritized genes. (A) Boxplot of -log10 FDR-adjusted P-value of multi-trait scans of GTAs in UKBB, grouped by 8 ICD code blocks across 9 genes with multiple TWAS GTAs across different trait categories. The red dotted line represents FDR-adjusted P = 0.05. (B) Forest plot of GTA association estimates and 95% FDR-adjusted confidence intervals for 6 neonatal traits in ELGAN for 9 genes with multiple TWAS GTAs across categories. The red line shows a null effect size of 0, and associations are colored blue for associations at FDR-adjusted P < 0.05. (C) Follow-up GBAT and Mendelian randomization (MR) analysis results using RICHS data. On the left, effect size and 95% adjusted confidence intervals from GBAT (X-axis) between GReX of TF-encoding genes and TWAS gene associations (pairs given on Y-axis). On the right, MR effect size and 95% adjusted confidence interval (X-axis) of TF-gene on TWAS gene (pairs on Y-axis). The red line shows a null effect size of 0, and associations are colored blue for associations at FDR-adjusted P < 0.05.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted April 14, 2021. ; https://doi.org/10.1101/2021.04.12.21255170 doi: medRxiv preprint . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review)  . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted April 14, 2021. ; https://doi.org/10.1101/2021.04.12.21255170 doi: medRxiv preprint Table 3: Susceptibility genes associated with multiple traits. TWAS gene, location (chromosome:start-end), and associated trait are provided with genetic correlations between traits at SNP level are provided if significant FDRadjusted < 0.05. Adult-onset asthma, Body fat percentage, Autism spectrum disorder AOA/BFP: 0.16 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review)