Skip to main content
medRxiv
  • Home
  • About
  • Submit
  • ALERTS / RSS
Advanced Search

3D genomic features across >50 diverse cell types reveal insights into the genomic architecture of childhood obesity

View ORCID ProfileKhanh B. Trang, Matthew C. Pahl, James A. Pippin, Chun Su, View ORCID ProfileSheridan H. Littleton, Prabhat Sharma, Nikhil N. Kulkarni, Louis R. Ghanem, Natalie A. Terry, Joan M. O’Brien, Yadav Wagley, Kurt D. Hankenson, Ashley Jermusyk, Jason W. Hoskins, Laufey T. Amundadottir, Mai Xu, View ORCID ProfileKevin M Brown, Stewart A. Anderson, Wenli Yang, Paul M. Titchenell, Patrick Seale, Laura Cook, View ORCID ProfileMegan K. Levings, View ORCID ProfileBabette S. Zemel, Alessandra Chesi, Andrew D. Wells, View ORCID ProfileStruan F.A. Grant
doi: https://doi.org/10.1101/2023.08.30.23294092
Khanh B. Trang
1Center for Spatial and Functional Genomics, The Children’s Hospital of Philadelphia, Philadelphia, PA, USA
2Division of Human Genetics, The Children’s Hospital of Philadelphia, Philadelphia, PA, USA
PhD
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Khanh B. Trang
Matthew C. Pahl
1Center for Spatial and Functional Genomics, The Children’s Hospital of Philadelphia, Philadelphia, PA, USA
2Division of Human Genetics, The Children’s Hospital of Philadelphia, Philadelphia, PA, USA
PhD
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
James A. Pippin
1Center for Spatial and Functional Genomics, The Children’s Hospital of Philadelphia, Philadelphia, PA, USA
2Division of Human Genetics, The Children’s Hospital of Philadelphia, Philadelphia, PA, USA
BA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Chun Su
1Center for Spatial and Functional Genomics, The Children’s Hospital of Philadelphia, Philadelphia, PA, USA
2Division of Human Genetics, The Children’s Hospital of Philadelphia, Philadelphia, PA, USA
PhD
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Sheridan H. Littleton
1Center for Spatial and Functional Genomics, The Children’s Hospital of Philadelphia, Philadelphia, PA, USA
2Division of Human Genetics, The Children’s Hospital of Philadelphia, Philadelphia, PA, USA
3Cell and Molecular Biology Graduate Group, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
4Department of Genetics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
BS
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Sheridan H. Littleton
Prabhat Sharma
1Center for Spatial and Functional Genomics, The Children’s Hospital of Philadelphia, Philadelphia, PA, USA
5Department of Pathology, The Children’s Hospital of Philadelphia, Philadelphia, PA, USA
PhD
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Nikhil N. Kulkarni
1Center for Spatial and Functional Genomics, The Children’s Hospital of Philadelphia, Philadelphia, PA, USA
5Department of Pathology, The Children’s Hospital of Philadelphia, Philadelphia, PA, USA
PhD
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Louis R. Ghanem
6Division of Gastroenterology, Hepatology, and Nutrition, Children’s Hospital of Philadelphia, PA, USA
MD, PhD
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Natalie A. Terry
6Division of Gastroenterology, Hepatology, and Nutrition, Children’s Hospital of Philadelphia, PA, USA
MD, PhD
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Joan M. O’Brien
7Scheie Eye Institute, Department of Ophthalmology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania, PA, USA
8Penn Medicine Center for Ophthalmic Genetics in Complex Disease
MD
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Yadav Wagley
9Department of Orthopedic Surgery University of Michigan Medical School Ann Arbor, MI, USA
PhD
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Kurt D. Hankenson
9Department of Orthopedic Surgery University of Michigan Medical School Ann Arbor, MI, USA
DVM, MS, PhD
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ashley Jermusyk
10Laboratory of Translational Genomics, Division of Cancer Epidemiology and Genetics, National Cancer Institute, Bethesda, MD, USA
PhD
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jason W. Hoskins
10Laboratory of Translational Genomics, Division of Cancer Epidemiology and Genetics, National Cancer Institute, Bethesda, MD, USA
PhD
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Laufey T. Amundadottir
10Laboratory of Translational Genomics, Division of Cancer Epidemiology and Genetics, National Cancer Institute, Bethesda, MD, USA
PhD
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Mai Xu
10Laboratory of Translational Genomics, Division of Cancer Epidemiology and Genetics, National Cancer Institute, Bethesda, MD, USA
PhD
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Kevin M Brown
10Laboratory of Translational Genomics, Division of Cancer Epidemiology and Genetics, National Cancer Institute, Bethesda, MD, USA
PhD
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Kevin M Brown
Stewart A. Anderson
11Department of Child and Adolescent Psychiatry, Children’s Hospital of Philadelphia, Philadelphia, PA, USA
12Department of Psychiatry, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
PhD
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Wenli Yang
13Institute for Diabetes, Obesity and Metabolism, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
14Department of Cell and Developmental Biology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
PhD
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Paul M. Titchenell
13Institute for Diabetes, Obesity and Metabolism, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
15Department of Physiology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
PhD
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Patrick Seale
13Institute for Diabetes, Obesity and Metabolism, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
14Department of Cell and Developmental Biology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
PhD
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Laura Cook
16Department of Microbiology and Immunology, University of Melbourne, Peter Doherty Institute for Infection and Immunity, Melbourne, VIC, Australia
17Department of Critical Care, Melbourne Medical School, University of Melbourne, Melbourne, VIC, Australia
18Division of Infectious Diseases, Department of Medicine, University of British Columbia, Vancouver, BC, Canada
PhD
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Megan K. Levings
19Department of Surgery, University of British Columbia, Vancouver, BC, Canada
20BC Children’s Hospital Research Institute, Vancouver, BC, Canada
21School of Biomedical Engineering, University of British Columbia, Vancouver, BC, Canada
PhD
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Megan K. Levings
Babette S. Zemel
6Division of Gastroenterology, Hepatology, and Nutrition, Children’s Hospital of Philadelphia, PA, USA
22Department of Pediatrics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
PhD
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Babette S. Zemel
Alessandra Chesi
1Center for Spatial and Functional Genomics, The Children’s Hospital of Philadelphia, Philadelphia, PA, USA
23Department of Pathology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
PhD
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Andrew D. Wells
1Center for Spatial and Functional Genomics, The Children’s Hospital of Philadelphia, Philadelphia, PA, USA
5Department of Pathology, The Children’s Hospital of Philadelphia, Philadelphia, PA, USA
23Department of Pathology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
24Institute for Immunology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
PhD
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Struan F.A. Grant
1Center for Spatial and Functional Genomics, The Children’s Hospital of Philadelphia, Philadelphia, PA, USA
2Division of Human Genetics, The Children’s Hospital of Philadelphia, Philadelphia, PA, USA
4Department of Genetics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
13Institute for Diabetes, Obesity and Metabolism, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
22Department of Pediatrics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
25Division Endocrinology and Diabetes, The Children’s Hospital of Philadelphia, Philadelphia, PA, USA
26Penn Neurodegeneration Genomics Center, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
PhD
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Struan F.A. Grant
  • For correspondence: grants{at}chop.edu
  • Abstract
  • Full Text
  • Info/History
  • Metrics
  • Supplementary material
  • Data/Code
  • Preview PDF
Loading

ABSTRACT

The prevalence of childhood obesity is increasing worldwide, along with the associated common comorbidities of type 2 diabetes and cardiovascular disease in later life. Motivated by evidence for a strong genetic component, our prior genome-wide association study (GWAS) efforts for childhood obesity revealed 19 independent signals for the trait; however, the mechanism of action of these loci remains to be elucidated. To molecularly characterize these childhood obesity loci we sought to determine the underlying causal variants and the corresponding effector genes within diverse cellular contexts. Integrating childhood obesity GWAS summary statistics with our existing 3D genomic datasets for 57 human cell types, consisting of high-resolution promoter-focused Capture-C/Hi-C, ATAC-seq, and RNA-seq, we applied stratified LD score regression and calculated the proportion of genome-wide SNP heritability attributable to cell type-specific features, revealing pancreatic alpha cell enrichment as the most statistically significant. Subsequent chromatin contact-based fine-mapping was carried out for genome-wide significant childhood obesity loci and their linkage disequilibrium proxies to implicate effector genes, yielded the most abundant number of candidate variants and target genes at the BDNF, ADCY3, TMEM18 and FTO loci in skeletal muscle myotubes and the pancreatic beta-cell line, EndoC-BH1. One novel implicated effector gene, ALKAL2 – an inflammation-responsive gene in nerve nociceptors – was observed at the key TMEM18 locus across multiple immune cell types. Interestingly, this observation was also supported through colocalization analysis using expression quantitative trait loci (eQTL) derived from the Genotype-Tissue Expression (GTEx) dataset, supporting an inflammatory and neurologic component to the pathogenesis of childhood obesity. Our comprehensive appraisal of 3D genomic datasets generated in a myriad of different cell types provides genomic insights into pediatric obesity pathogenesis.

Question What are the causal variants and corresponding effector genes conferring pediatric obesity susceptibility in different cellular contexts?

Findings Our method of assessing 3D genomic data across a range of cell types revealed heritability enrichment of childhood obesity variants, particularly within pancreatic alpha cells. The mapping of putative causal variants to cis-regulatory elements revealed candidate effector genes for cell types spanning metabolic, neural, and immune systems.

Meaning We gain a systemic view of childhood obesity genomics by leveraging 3D techniques that implicate regulatory regions harboring causal variants, providing insights into the disease pathogenesis across different cellular systems.

INTRODUCTION

The prevalence of obesity has risen significantly worldwide1, especially among children and adolescents2. Obesity is associated with chronic diseases, such as diabetes, cardiovascular diseases, and certain cancers3–6, along with mechanical issues including osteoarthritis and sleep apnea7.

Modern lifestyle factors, including physical inactivity, excessive caloric intake, and socioeconomic inequity, along with disrupted sleep and microbiome, represent environmental risk factors for obesity pathogenesis. However, genetics also play a significant role, with the estimated heritability ranging from 40% to 70%8–10. Studies show that body weight and obesity remain stable from infancy to adulthood11–14, but variation between individuals does exisit15. Genome-wide association studies (GWAS) have improved our understanding of the genetic contribution to childhood obesity16–21. However, the functional consequences and molecular mechanisms of identified genetic variants in such GWAS efforts are yet to be fully elucidated. Efforts are now being made to predict target effector genes and explore potential drug targets using various computational and experimental approaches22–26, which subsequently warrant functional follow-up efforts.

With our extensive datasets generated on a range of different cell types, by combining 3D chromatin maps (Hi-C, Capture-C) with matched transcriptome (RNA-seq) and chromatin accessibility data (ATAC-seq), we investigated heritability patterns of pediatric obesity-associated variants and their gene-regulatory functions in a cell type-specific manner. This approach yielded 94 candidate causal variants mapped to their putative effector gene(s) and corresponding cell type(s) setting. In addition, using methods comparable to our prior efforts in other disease contexts27–34, we also uncovered new variant-to-gene combinations within specific novel cellular settings, most notably in immune cell types, which further confirmed the involvement of the immune system in the pathogenesis of obesity in the early stages of life.

METHODS

Data and resource

Datasets used in prior studies are listed in eTable 1. ATAC-seq, RNA-seq, Hi-C, and Capture-C library generation for each cell type is provided in their original published study and their pre-processing pipelines and tools can be found in eMethods.

Definition of cis-Regulatory Elements (cREs)

We intersected ATAC-seq open chromatin regions (OCRs) of each cell type with chromatin conformation capture data determined by Hi-C/Capture-C of the same cell type, and with promoters (−1,500/+500bp of TSS) defined by GENCODE v30.

Childhood obesity GWAS summary statistics

Data on childhood obesity from the EGG consortium was downloaded from www.egg-consortium.org. We used 8,566,179 European ancestry variants (consisting of 8,613 cases and 12,696 controls in stage 1,921 cases and 1,930 controls in stage 2), representing ∼55% of the total 15,504,218 variants observed across all ancestries in the original study35. The sumstats file was reformatted by munge_sumstats.py to standardize with the weighted variants from HapMap v3 within the LDSC baseline, which reduced the variants to 1,217,311 (7.8% of total).

Cell type specific partitioned heritability

We used LDSC (http://www.github.com/bulik/ldsc) v.1.0.1 with --h2 flag to estimate SNP-based heritability of childhood obesity within 4 defined sets of input genomic regions: (1) OCRs, (2) OCRs at gene promoters, (3) cREs, and (4) cREs with an expanded window of ±500 bp. The baseline model LD scores, plink filesets, allele frequencies and variants weights files for the European 1000 genomes project phase 3 in hg38 were downloaded from the provided link (https://alkesgroup.broadinstitute.org/LDSCORE/GRCh38/). The cREs of each cell type were used to create the annotation, which in turn were used to compute annotation-specific LD scores for each cell types cREs set.

Genetic loci included in variant-to-genes mapping

19 sentinel signals that achieved genome-wide significance in the trans-ancestral meta-analysis study35 were leveraged for our analyses. Proxies for each sentinel SNP were queried using TopLD36 and LDlinkR tool37 with the GRCh38 Genome assembly, 1000 Genomes phase 3 v5 variant set, European population, and LD threshold of r2>0.8, which resulted in 771 proxies, including the 21 SNPs from the 99% credible set of the original study (eTable 2).

GWAS-eQTL colocalization

The summary statistics for the European ancestry subset from the EGG consortium GWAS for childhood obesity was used. Common variants (MAF□≥□0.01) from the 1000 Genomes Project v3 samples were used as a reference panel. We used non-overlapped genomic windows of ±250,000 bases extended in both directions from the median genomic position of each of 19 sentinel loci as input. We used ColocQuiaL38 to test genome-wide colocalization of all possible variants included in each inputted window against GTEx v.8 eQTLs associations for all 49 tissues available from https://www.gtexportal.org/home/datasets. A conditional posterior probability of colocalization of 0.8 or greater was imposed.

RESULTS

Enrichment assessment of childhood obesity variants across cell types

To explore the enrichment of childhood obesity GWAS variants across cell types, we carried out Partitioned Linkage Disequilibrium Score Regression (LDSR)39 on all ATAC-seq-defined OCRs for each cell type (Fig. 1A). We observed that 41 of 57 cell types – including 22 metabolic, 21 immune, 7 neural cell types and 7 independent cell lines (eTable 1) – showed at least a degree of directional enrichment with the total set of OCRs (Fig. 1B.a). However, only four cell types – two pancreatic alpha and two pancreatic beta cell-based datasets – had statistically significant enrichments (P<0.05). These enrichments were less pronounced when focusing on promoter OCRs only (Fig. 1B.b). To further limit the LD enrichment assessment to just those OCRs that can putatively regulate gene expression via chromatin contacts with gene promoters, we used the putative cREs27,29. This reduced the number of cell types showing at least nominal enrichment (31 of 57), enlarged the dispersion of enrichment ranges across different cell types, increased the 95% confidence intervals (CI) of enrichments, and hence increased the P-value of the resulting regression score. cREs from pancreatic alpha cells derived from single-cell ATAC-seq were the only dataset that remained statically significant (Fig. 1B.c).

Figure 1:
  • Download figure
  • Open in new tab
Figure 1:

Partitioned Linkage Disequilibrium Score Regression analysis for open chromatin regions of all cell types. A. Bar-plot shows total number of open chromatin regions (OCRs) identified by ATAC-seq for each cell type on bulk cells - blue, or on single cell – red; and portion of OCRs that fall within promoter-interacting region (PIR-OCRs) identified by hi-C – green, or by capture-C – orange. B. 4 panels of dot-plots show heritability enrichment by LDSC analysis for each cell type, with standard error whiskers. Dots’ colors corresponding to p-values, dots with white asterisk are significant p-value<0.05, dots’ sizes corresponding to proportion of SNP contribute for heritability. Dash line at 1, i.e., no enrichment. a) Analysis done on whole OCRs set of each cell type (whiskers colors match with bulk/single cell from bar-plot A); b) On only OCRs that overlapped with promoters (whiskers’ colors match with bulk/single cell from bar-plot A); c) On the putative cREs of each cell type (whiskers’ colors match with hi-C/capture-C from bar-plot A); d) On the same PIR-OCRs as middle panel with their genomic positions expanded ±500 bases on both sides (whiskers’ colors match with hi-C/capture-C from bar-plot A).

The original reported LDSR method analyzed enrichment in the 500bp flanking regions of their regulatory categories39. However, when we expanded our analysis to the ±500bp window for our cREs, albeit incorporating more weighted variants into the enrichment (represented by larger dots in Fig. 1B.d), this resulted in a decrease in the number of cell types yielding at least nominal enrichment (26 cell types), the enrichment range across cell types, the 95% CI, and level of significance. The pancreatic alpha cell observation also dropped below the bar for significance with this expanded window definition.

Consistency and diversity of childhood obesity proxy variants mapped to cREs

Despite the enrichments above only being limited to just a small number of cell types, it is likely that individual loci have differing levels of contributions in various cellular contexts and could not be detected at the genome wide assessment scale. As such we elected to further explore the candidate effector genes that are directly affected by cREs harboring childhood obesity-associated variants by systematically mapping the genomic positions of the LD proxies onto each cell type’s cREs, resulting in several scenarios, as illustrated in Fig. 2A,B. The resulting 94 proxies clustered at 13 original loci (eTable 3). eFig. 1A outlines the number of signals at each locus included or excluded based on the criteria we defined for our regions of interest. The TMEM18 locus yielded the most variants through cREs mapping, with 46 proxies for the two lead independent variants, rs7579427 and rs62104180. The second most abundant locus was ADCY3, with 21 proxies for lead variant rs4077678 (Fig. 2C). The higher number of variants at one locus did not correlate with implicating more genes or cell types through mapping. The mapping frequency of various variants within a specific locus exhibited substantial differences.

Figure 2:
  • Download figure
  • Open in new tab
Figure 2:

Mapping 771 proxies to the open chromatin regions of each cell type A. Schematic shows 4 different ways OCRs can be subset due to overlapping with chromatin contact sites and gene promoters, colors of sets of OCRs applied for the areas in B panel. B. Venn diagram shows how 771 proxies mapped onto the OCRs: a. Blue area: 758 proxies were located within contact regions of at least one cell type regardless of chromatin state; b. Red area: 417 proxies were located within contact regions marked as open by overlapping with OCR; c. Yellow area: If we only considered open chromatin regions, 178 proxies were included; d. Green area: To focus on just those variants residing within open chromatin and contacting promoter regions in any cell type, we overlapped the genomic positions of these proxies with each cell type’s cRE set, yielding 90 variants (3 from the 99% credible set) directly contacting open gene promoters (eTable 3), with 10 located within a promoter of one gene but contacting another different gene promoter. There were an additional 4 variants located within gene promoters but in chromatin contact with promoter(s) of nearby transcript(s) of the same gene. C. Bar-plot shows number of proxies, cell types and target genes mapped at each locus. D. The upSet plot shows the degree of overlap across cell types of the variants; ranked from the most common variant (red) – rs61888800 from BDNF locus, a well-known 5’ untranslated region variant of this gene that is associated with anti-depression and therapeutic response81,82 – appeared in 39 cell types, to the group of variants (grey) which appeared in only one cell type.

Inspecting individual variants regardless of their locus, we found that 28 of 94 proxies appeared in cREs across multiple cell types, with another 66 observed in just one cell type (Fig. 2D). 45 variants of these 66 just contacted one gene promoter, such as at the GPR1 and TFAP2B loci (eFig. 2).

Overall, the number of cell types in which a variant was observed in open chromatin correlated with the number of genes contacted via chromatin loops (eFig. 3A). However, we also observed that some variants found in cREs in multiple cell types were more selective with respect to their candidate effector genes (eFig. 3B-red arrow), or conversely, more selective across given cell types but implicated multiple genes (eFig. 3B-blue arrow). eFig. 4 outlines our observations at the TMEM18 locus – an example locus involved in both scenarios.

Figure 3:
  • Download figure
  • Open in new tab
Figure 3:

Profiles of 111 implicated genes by 94 proxies through cREs of each cell type Main panel: Bubble plot show corresponding expression level (size) and number of variants (color) target each implicated gene of each cell type. Squares represent genes with variants at their promoters. Circles represent genes with variants contacted through chromatin loops. Some genes were implicated by both types, these “double implications” are represented as diamond shapes, and were identified across several cell types: two cell types (plasmacytoid dendritic cells and pre-differentiated adipocytes) for ADCY3 gene, and five for BDNF (human embryonic stem cells - hESC, differentiated human fetal osteoblast cells - hFOB_Diff, neural progenitor cells derived from induced pluripotent stem cells - NPC_iPSC, PANC-1, and NCIH716 cell lines) Genes with expression undetected in our arrays are shown as triangles. Top panel: bar-plot shows numbers of cell types each gene was implicated within, color-coded by which systems the cell types belong to. Right panel: bar-plot shows numbers of genes implicated by the variants with each cell type.

Figure 4:
  • Download figure
  • Open in new tab
Figure 4:

Colocalization of target effector genes with eQTLs A. Venn diagram shows the overlaps between sets of genes yielded by ColocQuiaL and the variant-to-gene mapping process. B. Circos plot of the 10 loci demonstrates the differences in the ranges of associations between the two approaches, with long-ranged chromatin contacts between obesity variants and target genes displayed as orange links and short-range eQTLs colocalizations as green links. Two SNPs – rs35796073, and rs35142762 within the TMEM18 locus, in linkage disequilibrium with rs7579427 – were estimated with high probability (cond.PP.H4=0.78) of colocalizing with the expression of ALKAL2 gene in subcutaneous adipose tissue. These pairs of SNP-gene were also identified by our variant-to-gene mapping approach in natural killer cells, plasmacytoid dendritic cells, unstimulated PBMC naïve CD4 T cells and astrocytes. The rs7132908 variant at the FAIM2 locus colocalized with the expression of AQP6 in thyroid tissue and with ASIC1 in prostate tissue, not only with high cond.PP.H4 but also with high individual SNP causal probability (SNP.PP.H4 > 0.95). rs7132908 was the second most consistent observation in our variant-to-gene mapping, namely across 25 different cell types (Figure 2C) and all three systems plus the other independent cell lines. The pair of rs7132908-contacting-AQP6 was observed in 15 different cell types −8 metabolic and 4 neural cell types, and 3 independent cell lines. The pair of rs7132908-contacting-ASIC1 was observed in 11 different cell types - 8 metabolic and 2 neural cell types, and plasmacytoid dendritic cells. The other eQTL signals that overlapped with our variant-to-gene mapping results were: BDNF at the METTL15 locus with its promoter physically contacted by rs11030197 in 4 cell types and its expression significantly colocalized (cond.PP.H4=0.82) in tibial artery; ADCY9 at its locus with its promoter physically contacted by rs2531995 in natural killer cells and its expression significantly colocalized in skin tissue (“Skin_Not_Sun_Exposed_Suprapubic”, cond.PP.H4=0.97). And ADCY3 in the C panel. C. ColocQuiaL estimated that these SNPs highly colocalize with the expression of ADCY3 in 11 different tissues, where the overlapping with the 16 cell types is represented, color-coded by the proxies rs numbers.

Implicated genes cluster at loci strongly associated with childhood obesity consistently across multiple cell types

Mapping the variants across all the cell types resulted in a total of 111 implicated childhood obesity candidate effector genes (Table 1). Among these, 45 genes were specific to just one cell type (eFig. 5A), including 13 in myotubes and 7 in natural killer cells. Conversely and notably, BDNF appeared across 42 different cell types. Across the metabolic, neural, and immune systems and seven other cell lines, there were 9 genes consistently implicated in all four categories (top panel Fig. 3 – red stars, eFig. 5B: “all”), while 5 genes were consistently implicated in metabolic, neural, and immune systems (top panel Fig. 3 – blue stars, eFig. 5B: “all_main”). Two genes, ADCY3 and BDNF, had variants both at their promoters and contacted variants in cREs via chromatin loops (eFig. 6).

View this table:
  • View inline
  • View popup
Table 1 PubMed-query known functions for 111 genes implicated by obesity variants.

At the TMEM18 locus on chr 2p25.3, a highly significant human obesity locus that has long been associated with both adult and childhood obesity, we obeserved differing degrees of evidence for 16 genes, but noted that rs6548240, rs35796073, and rs35142762 consistently contacted the SH3YL1, ACP1, and ALKAL2 promoters across multiple cell types (Fig. 2D-third and fourth column).

At the chr 2p23 locus, ADCY3 yielded the most contacts (i.e. many proxies contacting the same gene via chromatin loops), suggesting this locus acts as a regulatory hub. However, we observed a similar composition in cell types for four other genes: DNAJC27, DNAJC27-AS1 (both previously implicated in obesity and/or diabetes traits40), AC013267.1, and SNORD14 (RF00016). ITSN2, NCOA1, and EFR3B were three genes within this locus that were only implicated in immune cell types. NCOA1 encodes a prominent meta-inflammation factor41 known to reduce adipogenesis and shift the energy balance between white and brown fat, and its absence known to induce obesity42.

CALCR was the most frequently implicated gene at its locus, supported by 20 cell types across all systems. While within the BDNF locus, METTL15 and KIF18A – two non-cell-type-specific genes - plus some lncRNA genes, were contacted by childhood obesity-associated proxies within the same multiple cell types as BNDF, again suggesting the presence of a regulatory hub.

At the FAIM2 locus on chr 12q13.12, we observed known genes associated with obesity, eating patterns, and diabetes-related traits, including ASIC1, AQP2, AQP5, AQP6, RACGAP1, and AC025154.2 (AQP5-AS1) along with FAIM2 (Table 1). These genes were harbored within PIR-OCRs of astrocytes, neural progenitors, hypothalamic neurons, and multiple metabolic cell types. Plasmacytoid and CD1c+ conventional dendritic cells were the only two immune cell types that harbored such proxies within their cREs, implicating ASIC1, PRPF40B, RPL35AP28, TMBIM6, and LSM6P2 at the FAIM2 locus.

The independent ADCY9 and FTO loci are both located on chromosome 16. Genes at the ADCY9 locus were only implicated in a subset of immune cell types. Interestingly, genes at the FTO locus were only implicated in Hi-C datasets (as opposed to Capture C), including 6 metabolic cell types and astrocytes. Most genes at the FTO locus were implicated in skeletal myotubes, differentiated osteoblasts, and astrocytes, namely FTO and IRX3; while IRX5, CRNDE, and AC106738.1 were also implicated in adipocytes and hepatocytes.

The most implicated cell types by two sets of analyses

EndoC-BH1 and myotubes are the two cell types in which we implicated the most effector genes, with 38 and 42, respectively – Fig. 3 side panel. This phenomenon is likely proportional in the case of myotubes, given the large number of cREs identified by overlapped Hi-C contact data and ATAC-seq open regions (Fig. 1A), but not for EndoC-BH1. Albeit harboring an average number of cREs compared to other cell types, EndoC-BH1 cells were consistently among the top-ranked heritability estimates for the childhood obesity variants resulting from the EGG consortium GWAS (Fig. 1) and harbored a significant number of implicated genes by the mapping of proxies. Interestingly, the pancreatic alpha cell type – shown above to be the most significant for heritability estimate by LDSC – revealed only 6 implicated genes contacted by the defined proxies, namely BDNF and five lncRNA genes.

Pathway analysis

Of the 111 implicated genes in total, PubMed query revealed functional studies for 66 genes. The remaining were principally lncRNA and miRNA genes with currently undefined functions (Table 1). To investigate how our implicated genes could confer obesity risk, we performed several pathway analyses keeping them either separated for each cell type or pooling into the respective metabolic, neural, or immune system sets. eFig. 7 shows simple Gene Ontology (GO) biological process terms enrichment results.

Leveraging the availability of our expression data generated via RNA-seq (available for 46 of 57 cell types), we performed pathway analysis. Given that our gene sets from the variant-to-gene process was stringently mapped, the sparse enrichment from normal direct analyses is not ideal for exploring obesity genetic etiology. Thus, we incorporated two methods from the pathfindR package43 and our customized SPIA (details in eMethods). The result of 60 enriched KEGG terms is shown in eFig. 8 (eTable 4), with 13 genes in 14 cell types for pathfindR and 39 enriched KEGG terms shown in eFig. 9 (eTable 5), with 10 genes in 42 cell types for customized SPIA. There were 20 overlapping pathways between the two approaches (yellow rows in eTable 4&5) including many signaling pathways such as the GnRH (hsa04912), cAMP (hsa04024), HIF-1 (hsa04066), Glucagon (hsa04922), Relaxin (hsa04926), Apelin (hsa04371), and Phospholipase D (hsa04072) signaling pathways. They were all driven by one or more of these 5 genes: ADCY3, ADCY9, CREBBP, MMP2, and NCOA1. Interestingly, we observed the involvement of natural killer cells in nearly all the enriched KEGG terms from pathfindR due to the high expression of the two adenylyl cyclase encoded genes, ADCY3 and ADCY9, along with CREBBP. The SPIA approach disregarded the aquaporin genes (given they appear so frequently in so many pathways that involve cellular channels) but highlighted the central role of BDNF which single-handedly drove four signaling pathways: the Ras, Neurotrophin, PI3K-Akt, and MAPK signaling pathways. This also revealed the role of TRAP1 in neurodegeneration.

These two approaches did not discount the role of FAIM2 and CALCR. However, their absence was mainly due to the content of the current KEGG database. On the other hand, these approaches accentuated the role of the MMP2 gene at the FTO locus in skeletal myotubes, given its consistency within the GnRH signaling pathway (eFig. 10), which is in line with previous studies linking its expression with obesity44–46.

Supportive evidence by colocalization of target effector genes with eQTLs

The GTEx consortium has characterized thousands of eQTLs, albeit in heterogeneous bulk tissues47. To assess how many observed gene-SNP pairs agreed with our physical variant-to-gene mapping approach in our multiple separate cellular settings, we performed colocalization analysis using ColocQuiaL38.

282 genes were reported to be associated with the variants within 13 loci from our variant-to-genes analysis. We found 114 colocalizations for ten of our loci that had high conditioned posterior probabilities (cond.PP.H4.abf ≥ 0.8), involving 44 genes and 41 tissues among the eQTLs. We extracted the posterior probabilities for each SNP within each colocalization and selected the 95% credible set as the likely causal variants (complete list in eTable 6). Despite sensitivity differences and varying cellular settings, when compared with our variant-to-gene mapping results, colocalization analysis yielded consistent identification for 21 pairs of SNP-gene interactions when considering the analyses across all our cell types, composed of 20 SNPs and 7 genes. Details of these SNP-gene pairs are shown in Fig. 4A and B.

Of these 20 SNPs, 15 were at the ADCY3 locus, in LD with sentinel variant rs4077678, and all implicated ADCY3 as the effector gene in 29 cell types – 15 metabolic, 6 immune, 4 neural cell types and 4 independent cell lines (Fig. 4C). Indeed, missense mutations have been previously reported for this gene in the context of obesity48,49 while another member of this gene family, ADCY5, has also been extensively implicated in metabolic traits50.

Predicting transcription factors (TFs) binding disruption at implicated genes contributing to obesity risk

TFs regulate gene expression by binding to DNA motifs at enhancers and silencers, where any disruption by a SNP can potentially cause dysregulation of a target gene. Thus, we used motifbreakR (R package) to predict such possible events at the loci identified by our variant-to-gene mapping. Each variant was predicted to disrupt the binding of several different TFs, thus requiring further literature cross-examination to select the most probable effects. For example, rs7132908 (consistently contacting FAIM2 in 25 cell types) was predicted to disrupt the binding of 12 different transcription factors. Among them, SREBF1 (eFig. 11A) was the only TF that concurred with evidence that it regulates AQP2 and FAIM2 at the same enhancer51. The full prediction list can be found in eTable 7.

To narrow down the list of putative TF binding sites at each variant position, we leveraged the ATAC-seq footprint analysis using the RGT suite52. The final set of Motif-Predicted Binding Sites (MPBS) within each cell type ATAC-seq footprints was used to overlap with the genomic locations of the OCRs, and then overlapped with our obesity variants, resulting in annotated 29 variants. Mosaic plot in eFig. 11B shows the number and proportions of variants predicted by motifbreakR and/or overlapped with MPBS. Insignificant P-value from Fisher’s exact test indicated the independence of the two analyses. Only seven variants were found within the cREs for the same TF motifs predicted to be disrupted by motifbreakR (eFig. 11C). eFig. 12 outlines the seven variants that motifbreakR and ATAC-seq footprint analysis agreed on the TF bindings they might disrupt.

DISCUSSION

Given the challenge of uncovering the underlying molecular mechanisms driving such a multifactorial disease as obesity, our approach leveraging GWAS summary statistics, RNA-seq, ATAC-seq, and promoter Capture C / Hi-C offers new insights. This is particularly true as it is becoming increasingly evident that multiple effector genes can operate in a temporal fashion at a given locus depending on cell state, including at the FTO locus53. Our approach offers an opportunity to implicate relevant cis-regulatory regions across different cell types contributing to the genetic etiology of the disease. By assigning GWAS signals to candidate causal variants and corresponding putative effector genes via open chromatin and chromatin contact information, we enhanced the fine-mapping process with an experimental genomic perspective to yield new insights into the biological pathways influencing childhood obesity.

LD score regression is a valuable method that estimates the relationship between linkage disequilibrium score and the summary statistics of GWAS SNPs to quantify the separate contributions of polygenic effects and various confounding factors that produce SNP-based heritability of disease. The general positive heritability enrichment across our open chromatin features spanning multiple cell types (Fig. 1B.a) reinforces the notion that obesity etiology involves many systems in our body.

While obesity has long been known to be a risk factor for pancreatitis and pancreatic cancer, the significant enrichment of pancreatic alpha and beta cell related 3D genomic features for childhood obesity GWAS signals demonstrates the bidirectional relationship between obesity and the pancreas; indeed, it is well established that insulin has obesogenic properties.

Moreover, the comorbidity of obesity and diabetes (either causal or a result of the overlap between SNPs associated with these two diseases) is tangible. When focusing on genetic annotation of the cREs only, the association with obesity became more diverse across cell types, especially in metabolic cells. Interestingly, the lack of enrichment (only 8 of 57 cell types yielded no degree of enrichment) of obesity SNPs heritability in open gene promoters (Fig. 1B.b) reveals that cRE regions harboring obesity SNPs are more involved in gene regulation than disruption, and therefore potentially contributing more weight to the manifestation of the disease.

Of course, we should factor in the effective sample sizes of the GWAS efforts that are wide-ranging (2,000-24,000 – given that the N for each variant is different within a single dataset, thus contributing to the weights and P-value of each SNP when the algorithm calculates the genome-wide heritability), which could result in noise and negative enrichment observed in the analysis – a methodology limitation of partial linkage regression that has been extensively discussed in the field54. Thus, it is crucial to interpret the enrichment (or lack thereof) of disease variants in a certain cellular setting with an ad hoc biological context.

From mapping the common proxies of 19 independent sentinel SNPs that were genome-wide significantly associated with childhood obesity to putative effector genes through chromatin contacting cREs, one striking finding was the several potential “hubs” of putatively core effector genes, whose occurrence spread across three human physiological systems. With the data available from so many cell types, our approach connected new candidate causal variants to known obesity-related genes and new implications of cell modality for previously known associations.

A potential application of this association could be to fine-tune the effect of a drug toward controlling appetite. An example of bringing new aspects to the old is for the signal within the FTO locus that contacted IRX3 and IRX5: previous studies have suggested these obesogenic effects operate in adipocytes55, brain56, or pancreas57; here we confirmed this association in adipocytes and uncover the presence of distal chromatin contacts in myotubes for the first time.

Besides the above-mentioned genes with known associations with obesity, we discovered newly implicated genes. For example, the LRRIQ3 gene at the TNNI3K locus had its open promoter contacted by two SNPs, rs1040070 and rs10493544, in NTERA2 cells only. The published studies58,59 that associated LRRIQ3 with major depressive disorder and opioid usage acknowledged the overlapping promoter of this gene, albeit in the opposite direction, with a run-through transcript of FPGT-TNNI3K – previously shown to be associated with BMI in European60 and Korean populations61.

It is apparent that not all the implicated genes we report would contribute equally to the susceptibility of obesity pathogenesis. Each locus comprises genes whose functions are obviously related to obesity or similar traits like BMI, fat weight, etc., while other genes are not so directly obvious in their relation to these traits.

It is encouraging that for implicated genes within these multi-cell-type loci across different physiological systems we could find previous associations to the corresponding cell types or systems. Examples are the two aforementioned genes at the TMEM18 locus (SH3YL1 and ACP1)62–66 with the broad spectrum of their functions, HEPACAM2 implicated in the NCIH716 cell line at the CALCR locus67,68, and LRRIQ3 in the NTERA2 cell line at theTNNI3K locus69.

Chronic inflammation is an essential characteristic of obesity pathogenesis. Adipose tissue-resident immune cells have been observed, leading to an increased focus in recent years on their potential contribution to metabolic dysfunction. On the other hand, neurological or psychological conditions, such as stress, induce the secretion of both glucocorticoids (increase motivation for food) and insulin (promotes food intake and obesity). Pleasure feeding then reduces activity in the stress-response network, reinforcing the feeding habit. It has been shown that voluntary behaviors, stimulated by external or internal stressors or pleasurable feelings, memories, and habits, can override the basic homeostatic controls of energy balance70. The potential link between the immune system and metabolic disease, and moreover, through the neural system, was tangible in our findings.

Two of the three SNPs which ranked the third most consistent in our variant-to-gene mapping (Fig. 2C) – rs35796073 and rs35142762 – contacted the ALKAL2 promoter (supported by GTEx evidence to colocalize with ALKAL2 expression). The anaplastic lymphoma kinase (encoded by ALK gene) is a receptor tyrosine kinase, belongs to the insulin receptor family, and has been reported to promote nerve cell growth and differentiation71,72. Despite ALKAL2 (ALK and LTK ligand 2) being studied principally in the context of immunity, a recent study using the EGCUT biobank GWAS identified ALK as a candidate thinness gene and genetic deletion showed that its expression in hypothalamic neurons acts as a negative regulator in controlling energy expenditure via sympathetic control of adipose tissue lipolysis73. ALKAL2 – encoding a high-affinity agonist of ALK/LTK receptors – which has been reported to enhance expression in response to inflammatory pain in nociceptors74,75 - has been recently implicated as a novel candidate gene for childhood BMI by transcriptome-wide association study76, and achieved genome-wide significance in a GWAS study contrasting persistent healthy thinness with severe early-onset obesity using the STILTS and SCOOP cohorts77. The finding that overexpression of ALKAL2 could potentiate neuroblastoma progression in the absence of ALK mutation78 echoes the relationship between ADCY3 and MC4R79, where a peripheral gene, ADCY3, can regulate/impair the function of a core gene, i.e. MC4R, within the energy-regulating melanocortin signaling pathway80.

Our approach implicates putative target genes based on a mechanism of regulation for these variants to alter gene expression – through regulator TF(s) that bind to these contact sites. A potential limitation of the predictions from motifbreakR and matching TF motifs to ATAC-seq footprint by the RGT toolkit is that they were both based on the position probability matrixes of Jaspar and Hocomoco, which come from public motif databases. The ATAC-seq footprint analysis also carries sequence bias that can lead to false positive discovery. Thus, our attempt to call such regulators by predicting TF binding disruption can only serve as nominations – but warrant further functional follow up.

Another limitation of this work is the diversity in data quality among different samples, since different datasets were sampled and collected at different time points, from different patients, using different protocols, with libraries sequenced at different depths and qualities, and initially preprocessed with different pipelines and parameters. Thus, it is crucial to keep in mind that the discrepancy in data points might have resulted from variations in data quality. Importantly, any association discovered must be validated functionally before effector genes of the genetic variants can be leveraged to develop new therapies. Their putative function(s) must be characterized, together with the mechanism whereby the given variant’s alleles differentially affect the expression of the targeted genes. The next step is to explore how the target genes affect the trait of interest more directly.

Our results have provided a set of leads for future exploratory experiments in specific cellular settings in order to further expand our knowledge of childhood obesity genomics and hence equip us with more effective means to overcome the burden of this systematic disease.

CONCLUSION

Our approach of combining RNA-seq, ATAC-seq, and promoter Capture C/Hi-C datasets with GWAS summary statistics offers a systemic view of the multi-cellular nature of childhood obesity, shedding light on potential regulatory regions and effector genes. By leveraging physical properties, such as open chromatin status and chromatin contacts, we enhanced the fine-mapping process and gained new insights into the biological pathways influencing the disease. Although further functional validation is required, our findings provide valuable leads together with their cellular contexts for future research and the development of more effective strategies to address the burden of childhood obesity.

Data Availability

All data produced in the present study are available upon reasonable request to the authors

ACKNOWLEDGMENTS

This work was supported by National Institutes of Health awards R01 HD056465, R01 DK122586 and UM1 DK126194, and the Daniel B. Burke Endowed Chair for Diabetes Research. Given the use of de-identified datasets and biospecimens was not considered human subjects research, ethical oversight was waived by the Institutional Review Board of the Children’s Hospital of Philadelphia

Footnotes

  • Abstract updated

REFERENCES

  1. 1.↵
    Worldwide trends in body-mass index, underweight, overweight, and obesity from 1975 to 2016: a pooled analysis of 2416 population-based measurement studies in 128·9 million children, adolescents, and adults. Lancet. Dec 16 2017;390(10113):2627–2642. doi:10.1016/s0140-6736(17)32129-3
    OpenUrlCrossRef
  2. 2.↵
    National Health and Nutrition Examination Survey 2017–March 2020 Prepandemic Data Files Development of Files and Prevalence Estimates for Selected Health Outcomes, doi:10.15620/cdc:106273 (2021). https://stacks.cdc.gov/view/cdc/106273
    OpenUrlCrossRef
  3. 3.↵
    Health Effects of Overweight and Obesity in 195 Countries over 25 Years. New England Journal of Medicine. 2017;377(1):13–27. doi:10.1056/NEJMoa1614362
    OpenUrlCrossRefPubMed
  4. 4.
    Singh GM, Danaei G, Farzadfar F, et al. The Age-Specific Quantitative Effects of Metabolic Risk Factors on Cardiovascular Diseases and Diabetes: A Pooled Analysis. PLOS ONE. 2013;8(7):e65174. doi:10.1371/journal.pone.0065174
    OpenUrlCrossRefPubMed
  5. 5.
    The Emerging Risk Factors C. Separate and combined associations of body-mass index and abdominal adiposity with cardiovascular disease: collaborative analysis of 58 prospective studies. The Lancet. 2011/03/26/ 2011;377(9771):1085–1095. doi:10.1016/S0140-6736(11)60105-0
    OpenUrlCrossRefPubMedWeb of Science
  6. 6.↵
    Lauby-Secretan B, Scoccianti C, Loomis D, Grosse Y, Bianchini F, Straif K. Body Fatness and Cancer — Viewpoint of the IARC Working Group. New England Journal of Medicine. 2016;375(8):794–798. doi:10.1056/NEJMsr1606602
    OpenUrlCrossRefPubMed
  7. 7.↵
    Fontaine KR, Barofsky I. Obesity and health-related quality of life. Obes Rev. Aug 2001;2(3):173–82. doi:10.1046/j.1467-789x.2001.00032.x
    OpenUrlCrossRefPubMed
  8. 8.↵
    Loos RJF, Yeo GSH. The genetics of obesity: from discovery to biology. Nature Reviews Genetics. 2022/02/01 2022;23(2):120–133. doi:10.1038/s41576-021-00414-z
    OpenUrlCrossRef
  9. 9.
    Maes HH, Neale MC, Eaves LJ. Genetic and environmental factors in relative body weight and human adiposity. Behavior genetics. Jul 1997;27(4):325–51.
    OpenUrlCrossRefPubMedWeb of Science
  10. 10.↵
    Elks CE, Den Hoed M, Zhao JH, et al. Variability in the heritability of body mass index: a systematic review and meta-regression. Frontiers in endocrinology. 2012;3:29.
    OpenUrl
  11. 11.↵
    Demerath EW, Choh AC, Czerwinski SA, et al. Genetic and environmental influences on infant weight and weight change: the Fels Longitudinal Study. Am J Hum Biol. Sep-Oct 2007;19(5):692–702. doi:10.1002/ajhb.20660
    OpenUrlCrossRefPubMedWeb of Science
  12. 12.
    Dubois L, Girard M, Girard A, Tremblay R, Boivin M, Pérusse D. Genetic and environmental influences on body size in early childhood: a twin birth-cohort study. Twin Res Hum Genet. Jun 2007;10(3):479–85. doi:10.1375/twin.10.3.479
    OpenUrlCrossRefPubMed
  13. 13.
    Wardle J, Carnell S, Haworth CM, Plomin R. Evidence for a strong genetic influence on childhood adiposity despite the force of the obesogenic environment. Am J Clin Nutr. Feb 2008;87(2):398–404. doi:10.1093/ajcn/87.2.398
    OpenUrlAbstract/FREE Full Text
  14. 14.↵
    Bouchard C. Childhood obesity: are genetic differences involved? The American Journal of Clinical Nutrition. 2009;89(5):1494S–1501S. doi:10.3945/ajcn.2009.27113C
    OpenUrlAbstract/FREE Full Text
  15. 15.↵
    Littleton SH, Berkowitz RI, Grant SFA. Genetic Determinants of Childhood Obesity. Mol Diagn Ther. Dec 2020;24(6):653–663. doi:10.1007/s40291-020-00496-1
    OpenUrlCrossRef
  16. 16.↵
    Vogelezang S, Bradfield JP, Ahluwalia TS, et al. Novel loci for childhood body mass index and shared heritability with adult cardiometabolic traits. PLOS Genetics. 2020;16(10):e1008718. doi:10.1371/journal.pgen.1008718
    OpenUrlCrossRef
  17. 17.
    Yaghootkar H, Zhang Y, Spracklen CN, et al. Genetic Studies of Leptin Concentrations Implicate Leptin in the Regulation of Early Adiposity. Diabetes. Dec 2020;69(12):2806–2818. doi:10.2337/db20-0070
    OpenUrlAbstract/FREE Full Text
  18. 18.
    Couto Alves A, De Silva NMG, Karhunen V, et al. GWAS on longitudinal growth traits reveals different genetic factors influencing infant, child, and adult BMI. Sci Adv. Sep 2019;5(9):eaaw3095. doi:10.1126/sciadv.aaw3095
    OpenUrlFREE Full Text
  19. 19.
    Fu J, Wang Y, Li G, et al. Childhood sleep duration modifies the polygenic risk for obesity in youth through leptin pathway: the Beijing Child and Adolescent Metabolic Syndrome cohort study. Int J Obes (Lond). Aug 2019;43(8):1556–1567. doi:10.1038/s41366-019-0405-1
    OpenUrlCrossRef
  20. 20.
    Turcot V, Lu Y, Highland HM, et al. Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity. Nat Genet. Jan 2018;50(1):26–41. doi:10.1038/s41588-017-0011-x
    OpenUrlCrossRefPubMed
  21. 21.↵
    Littleton SH, Grant SFA. Strategies to identify causal common genetic variants and corresponding effector genes for paediatric obesity. Pediatr Obes. Dec 2022;17(12):e12968. doi:10.1111/ijpo.12968
    OpenUrlCrossRef
  22. 22.↵
    Yu F, Cato LD, Weng C, et al. Variant to function mapping at single-cell resolution through network propagation. Nature Biotechnology. 2022/06/06 2022;doi:10.1038/s41587-022-01341-y
    OpenUrlCrossRef
  23. 23.
    Avsec Ž, Agarwal V, Visentin D, et al. Effective gene expression prediction from sequence by integrating long-range interactions. Nature Methods. 2021/10/01 2021;18(10):1196–1203. doi:10.1038/s41592-021-01252-x
    OpenUrlCrossRefPubMed
  24. 24.
    Gazal S, Weissbrod O, Hormozdiari F, et al. Combining SNP-to-gene linking strategies to identify disease genes and assess disease omnigenicity. Nature Genetics. 2022/06/01 2022;54(6):827–836. doi:10.1038/s41588-022-01087-y
    OpenUrlCrossRef
  25. 25.
    Zhou J, Troyanskaya OG. Predicting effects of noncoding variants with deep learning– based sequence model. Nature Methods. 2015/10/01 2015;12(10):931–934. doi:10.1038/nmeth.3547
    OpenUrlCrossRefPubMed
  26. 26.↵
    Nasser J, Bergman DT, Fulco CP, et al. Genome-wide enhancer maps link risk variants to disease genes. Nature. 2021/05/01 2021;593(7858):238–243. doi:10.1038/s41586-021-03446-x
    OpenUrlCrossRef
  27. 27.↵
    Chesi A, Wagley Y, Johnson ME, et al. Genome-scale Capture C promoter interactions implicate effector genes at GWAS loci for bone mineral density. OriginalPaper. Nature Communications. 2019-03-19 2019;10(1):1–11. doi:doi:10.1038/s41467-019-09302-x
    OpenUrlCrossRef
  28. 28.
    Su C, Johnson ME, Torres A, et al. Mapping effector genes at lupus GWAS loci using promoter Capture-C in follicular helper T cells. Nature Communications. 2020/07/03 2020;11(1):3294. doi:10.1038/s41467-020-17089-5
    OpenUrlCrossRefPubMed
  29. 29.↵
    Pahl MC, Doege CA, Hodge KM, et al. Cis-regulatory architecture of human ESC-derived hypothalamic neuron differentiation aids in variant-to-gene mapping of relevant complex traits. OriginalPaper. Nature Communications. 2021-11-19 2021;12(1):1–12. doi:doi:10.1038/s41467-021-27001-4
    OpenUrlCrossRef
  30. 30.
    Cousminer DL, Wagley Y, Pippin JA, et al. Genome-wide association study implicates novel loci and reveals candidate effector genes for longitudinal pediatric bone accrual. Genome Biol. Jan 4 2021;22(1):1. doi:10.1186/s13059-020-02207-9
    OpenUrlCrossRefPubMed
  31. 31.
    Vujkovic M, Ramdas S, Lorenz KM, et al. A multiancestry genome-wide association study of unexplained chronic ALT elevation as a proxy for nonalcoholic fatty liver disease with histological and radiological validation. Nat Genet. Jun 2022;54(6):761–771. doi:10.1038/s41588-022-01078-z
    OpenUrlCrossRefPubMed
  32. 32.
    Pahl MC, Le Coz C, Su C, et al. Implicating effector genes at COVID-19 GWAS loci using promoter-focused Capture-C in disease-relevant immune cell types. Genome Biology. 2022/06/03 2022;23(1):125. doi:10.1186/s13059-022-02691-1
    OpenUrlCrossRef
  33. 33.
    Su C, Gao L, May CL, et al. 3D chromatin maps of the human pancreas reveal lineage-specific regulatory architecture of T2D risk. Cell Metabolism. 6 September 2022 2022;34(9)1409. doi:10.1016/j.cmet.2022.08.014
    OpenUrlCrossRef
  34. 34.↵
    Palermo J, Chesi A, Zimmerman A, et al. Variant-to-gene mapping followed by cross-species genetic screening identifies GPI-anchor biosynthesis as a regulator of sleep. Sci Adv. Jan 6 2023;9(1):eabq0844. doi:10.1126/sciadv.abq0844
    OpenUrlCrossRef
  35. 35.↵
    Bradfield JP, Vogelezang S, Felix JF, et al. A trans-ancestral meta-analysis of genome-wide association studies reveals loci associated with childhood obesity. Human Molecular Genetics. 2019;28(19):3327–3338. doi:10.1093/hmg/ddz161
    OpenUrlCrossRef
  36. 36.↵
    Huang L, Rosen JD, Sun Q, et al. TOP-LD: A tool to explore linkage disequilibrium with TOPMed whole-genome sequence data. The American Journal of Human Genetics. 2022/06/02/ 2022;109(6):1175–1181. doi:10.1016/j.ajhg.2022.04.006
    OpenUrlCrossRef
  37. 37.↵
    Myers TA, Chanock SJ, Machiela MJ. LDlinkR: An R Package for Rapidly Calculating Linkage Disequilibrium Statistics in Diverse Populations. Front Genet. 2020;11:157. doi:10.3389/fgene.2020.00157
    OpenUrlCrossRefPubMed
  38. 38.↵
    Chen BY, Bone WP, Lorenz K, Levin M, Ritchie MD, Voight BF. ColocQuiaL: a QTL-GWAS colocalization pipeline. Bioinformatics. 2022;38(18):4409–4411. doi:10.1093/bioinformatics/btac512
    OpenUrlCrossRef
  39. 39.↵
    Finucane HK, Bulik-Sullivan B, Gusev A, et al. Partitioning heritability by functional annotation using genome-wide association summary statistics. OriginalPaper. Nature Genetics. 2015-09-28 2015;47(11):1228–1235. doi:doi:10.1038/ng.3404
    OpenUrlCrossRefPubMed
  40. 40.↵
    Cherian PT, Al-Khairi I, Sriraman D, et al. Increased Circulation and Adipose Tissue Levels of DNAJC27/RBJ in Obesity and Type 2-Diabetes. Front Endocrinol (Lausanne). 2018;9:423. doi:10.3389/fendo.2018.00423
    OpenUrlCrossRef
  41. 41.↵
    Rollins DA, Coppo M, Rogatsky I. Minireview: nuclear receptor coregulators of the p160 family: insights into inflammation and metabolism. Mol Endocrinol. Apr 2015;29(4):502–17. doi:10.1210/me.2015-1005
    OpenUrlCrossRefPubMed
  42. 42.↵
    Mohsen G A-M, Abu-Taweel GM, Rajagopal R, et al. Betulinic acid lowers lipid accumulation in adipocytes through enhanced NCoA1–PPARγ interaction. Journal of Infection and Public Health. 2019/09/01/ 2019;12(5):726–732. doi:10.1016/j.jiph.2019.05.011
    OpenUrlCrossRef
  43. 43.↵
    Ulgen E, Ozisik O, Sezerman OU. pathfindR: An R Package for Comprehensive Identification of Enriched Pathways in Omics Data Through Active Subnetworks. Methods. Frontiers in Genetics. 2019-September-25 2019;10 doi:10.3389/fgene.2019.00858
    OpenUrlCrossRefPubMed
  44. 44.↵
    Derosa G, Ferrari I, D’Angelo A, et al. Matrix metalloproteinase-2 and −9 levels in obese patients. Endothelium. Jul-Aug 2008;15(4):219–24. doi:10.1080/10623320802228815
    OpenUrlCrossRefPubMedWeb of Science
  45. 45.
    Sezgin SBA, Bayoglu B, Ersoz F, et al. Downregulation of MMP-2 and MMP-9 genes in obesity patients and their relation with obesity-related phenotypes. Turkish Journal of Biochemistry. 2022;47(4):425–433. doi:doi:10.1515/tjb-2021-0124
    OpenUrlCrossRef
  46. 46.↵
    Nonino CB, Noronha NY, de Araújo Ferreira-Julio M, et al. Differential Expression of MMP2 and TIMP2 in Peripheral Blood Mononuclear Cells After Roux-en-Y Gastric Bypass. Brief Research Report. Frontiers in Nutrition. 2021-October-13 2021;8doi:10.3389/fnut.2021.628759
    OpenUrlCrossRef
  47. 47.↵
    Mohammadi Pejman 5 6 Park YoSon 11 Parsana Princy 12 Segrè Ayellet V. 1 Strober Benjamin J. 9 Zappala Zachary 7 8 GCLaAFBAACSEDJRHYJB, P. 19 Volpi Simona 19 NpmAAGPKSLARLNCMHMRASJ, 16 PSLBMEBPA, 137 NCFNCR. Genetic effects on gene expression across human tissues. Nature. 2017;550(7675):204–213.
    OpenUrlCrossRefPubMedWeb of Science
  48. 48.↵
    Grarup N, Moltke I, Andersen MK, et al. Loss-of-function variants in ADCY3 increase risk of obesity and type 2 diabetes. Nature genetics. 2018;50(2):172–174. doi:10.1038/s41588-017-0022-7
    OpenUrlCrossRefPubMed
  49. 49.↵
    Stergiakouli E, Gaillard R, Tavaré JM, et al. Genome-wide association study of height-adjusted BMI in childhood identifies functional variant in ADCY3. Obesity (Silver Spring). Oct 2014;22(10):2252–9. doi:10.1002/oby.20840
    OpenUrlCrossRefPubMed
  50. 50.↵
    Sinnott-Armstrong N, Sousa IS, Laber S, et al. A regulatory variant at 3q21.1 confers an increased pleiotropic risk for hyperglycemia and altered bone mineral density. Cell Metabolism. 2021/03/02/ 2021;33(3):615–628.e13. doi:10.1016/j.cmet.2021.01.001
    OpenUrlCrossRefPubMed
  51. 51.↵
    Kikuchi H, Jung HJ, Raghuram V, et al. Bayesian identification of candidate transcription factors for the regulation of Aqp2 gene expression. Am J Physiol Renal Physiol. Sep 1 2021;321(3):F389–f401. doi:10.1152/ajprenal.00204.2021
    OpenUrlCrossRef
  52. 52.↵
    Li Z, Schulz MH, Look T, Begemann M, Zenke M, Costa IG. Identification of transcription factor binding sites using ATAC-seq. Genome Biology. 2019/02/26 2019;20(1):45. doi:10.1186/s13059-019-1642-2
    OpenUrlCrossRefPubMed
  53. 53.↵
    Sobreira DR, Joslin AC, Zhang Q, et al. Extensive pleiotropism and allelic heterogeneity mediate metabolic effects of IRX3 and IRX5. Science (New York, NY). 2021;372(6546):1085–1091. doi:10.1126/science.abf1008
    OpenUrlAbstract/FREE Full Text
  54. 54.↵
    Steinsaltz D, Dahl A, Wachter KW. On Negative Heritability and Negative Estimates of Heritability. Genetics. 2020;215(2):343–357. doi:10.1534/genetics.120.303161
    OpenUrlAbstract/FREE Full Text
  55. 55.↵
    Claussnitzer M, Hui CC, Kellis M. FTO Obesity Variant and Adipocyte Browning in Humans. N Engl J Med. Jan 14 2016;374(2):192–3. doi:10.1056/NEJMc1513316
    OpenUrlCrossRefPubMed
  56. 56.↵
    Smemo S, Tena JJ, Kim K-H, et al. Obesity-associated variants within FTO form long-range functional connections with IRX3. Nature. 2014/03/01 2014;507(7492):371–375. doi:10.1038/nature13138
    OpenUrlCrossRefPubMedWeb of Science
  57. 57.↵
    Ragvin A, Moro E, Fredman D, et al. Long-range gene regulation links genomic type 2 diabetes and obesity risk regions to HHEX, SOX4, and IRX3. Proc Natl Acad Sci U S A. Jan 12 2010;107(2):775–80. doi:10.1073/pnas.0911591107
    OpenUrlAbstract/FREE Full Text
  58. 58.↵
    Johnston KJA, Adams MJ, Nicholl BI, et al. Identification of novel common variants associated with chronic pain using conditional false discovery rate analysis with major depressive disorder and assessment of pleiotropic effects of LRFN5. Transl Psychiatry. Nov 20 2019;9(1):310. doi:10.1038/s41398-019-0613-4
    OpenUrlCrossRef
  59. 59.↵
    Sanchez-Roige S, Fontanillas P, Jennings MV, et al. Genome-wide association study of problematic opioid prescription use in 132,113 23andMe research participants of European ancestry. Mol Psychiatry. Nov 2021;26(11):6209–6217. doi:10.1038/s41380-021-01335-3
    OpenUrlCrossRef
  60. 60.↵
    Graff M, Ngwa JS, Workalemahu T, et al. Genome-wide analysis of BMI in adolescents and young adults reveals additional insight into the effects of genetic loci over the life course. Hum Mol Genet. Sep 1 2013;22(17):3597–607. doi:10.1093/hmg/ddt205
    OpenUrlCrossRefPubMed
  61. 61.↵
    Lee JS, Cheong HS, Shin HD. BMI prediction within a Korean population. PeerJ. 2017;5:e3510. doi:10.7717/peerj.3510
    OpenUrlCrossRef
  62. 62.↵
    Fernandes SJ, Morikawa H, Ewing E, et al. Non-parametric combination analysis of multiple data types enables detection of novel regulatory mechanisms in T cells of multiple sclerosis patients. Sci Rep. Aug 19 2019;9(1):11996. doi:10.1038/s41598-019-48493-7
    OpenUrlCrossRef
  63. 63.
    Blessing AM, Ganesan S, Rajapakshe K, et al. Identification of a Novel Coregulator, SH3YL1, That Interacts With the Androgen Receptor N-Terminus. Mol Endocrinol. Oct 2015;29(10):1426–39. doi:10.1210/me.2015-1079
    OpenUrlCrossRef
  64. 64.
    Kobayashi M, Harada K, Negishi M, Katoh H. Dock4 forms a complex with SH3YL1 and regulates cancer cell migration. Cell Signal. May 2014;26(5):1082–8. doi:10.1016/j.cellsig.2014.01.027
    OpenUrlCrossRefPubMed
  65. 65.
    Choi GS, Min HS, Cha JJ, et al. SH3YL1 protein as a novel biomarker for diabetic nephropathy in type 2 diabetes mellitus. Nutr Metab Cardiovasc Dis. Feb 8 2021;31(2):498–505. doi:10.1016/j.numecd.2020.09.024
    OpenUrlCrossRef
  66. 66.↵
    Gaynor SC, Monson ET, Gaine ME, et al. Male-specific association of the 2p25 region with suicide attempt in bipolar disorder. J Psychiatr Res. Feb 2020;121:151–158. doi:10.1016/j.jpsychires.2019.11.009
    OpenUrlCrossRef
  67. 67.↵
    Wu Z, Liu Z, Ge W, et al. Analysis of potential genes and pathways associated with the colorectal normal mucosa-adenoma-carcinoma sequence. Cancer Med. Jun 2018;7(6):2555–2566. doi:10.1002/cam4.1484
    OpenUrlCrossRef
  68. 68.↵
    Huang Z, Yang Q, Huang Z. Identification of Critical Genes and Five Prognostic Biomarkers Associated with Colorectal Cancer. Med Sci Monit. Jul 5 2018;24:4625–4633. doi:10.12659/msm.907224
    OpenUrlCrossRef
  69. 69.↵
    Pleasure SJ, Lee VM. NTera 2 cells: a human cell line which displays characteristics expected of a human committed neuronal progenitor cell. J Neurosci Res. Aug 15 1993;35(6):585–602. doi:10.1002/jnr.490350603
    OpenUrlCrossRefPubMedWeb of Science
  70. 70.↵
    Dallman MF. Stress-induced obesity and the emotional nervous system. Trends Endocrinol Metab. Mar 2010;21(3):159–65. doi:10.1016/j.tem.2009.10.004
    OpenUrlCrossRefPubMedWeb of Science
  71. 71.↵
    Iwahara T, Fujimoto J, Wen D, et al. Molecular characterization of ALK, a receptor tyrosine kinase expressed specifically in the nervous system. Oncogene. Jan 30 1997;14(4):439–49. doi:10.1038/sj.onc.1200849
    OpenUrlCrossRefPubMedWeb of Science
  72. 72.↵
    Motegi A, Fujimoto J, Kotani M, Sakuraba H, Yamamoto T. ALK receptor tyrosine kinase promotes cell growth and neurite outgrowth. J Cell Sci. Jul 1 2004;117(Pt 15):3319–29. doi:10.1242/jcs.01183
    OpenUrlAbstract/FREE Full Text
  73. 73.↵
    Orthofer M, Valsesia A, Mägi R, et al. Identification of ALK in Thinness. Cell. 2020/06/11/ 2020;181(6):1246–1262.e22. doi:10.1016/j.cell.2020.04.034
    OpenUrlCrossRef
  74. 74.↵
    Defaye M, Iftinca MC, Gadotti VM, et al. The neuronal tyrosine kinase receptor ligand ALKAL2 mediates persistent pain. The Journal of Clinical Investigation. 06/15/ 2022;132(12)doi:10.1172/JCI154317
    OpenUrlCrossRef
  75. 75.↵
    Sun W, Yang S, Wu S, et al. Transcriptome analysis reveals dysregulation of inflammatory and neuronal function in dorsal root ganglion of paclitaxel-induced peripheral neuropathy rats. Mol Pain. May 24 2022:17448069221106167. doi:10.1177/17448069221106167
    OpenUrlCrossRef
  76. 76.↵
    Yao S, Wu H, Ding J-M, et al. Transcriptome-wide association study identifies multiple genes associated with childhood body mass index. International Journal of Obesity. 2021/05/01 2021;45(5):1105–1113. doi:10.1038/s41366-021-00780-y
    OpenUrlCrossRef
  77. 77.↵
    Riveros-McKay FA-O, Mistry V, Bounds RA-O, et al. Genetic architecture of human thinness compared to severe obesity. 2019;(1553–7404 (Electronic))
  78. 78.↵
    Borenäs M, Umapathy G, Lai W-Y, et al. ALK ligand ALKAL2 potentiates MYCN-driven neuroblastoma in the absence of ALK mutation. The EMBO Journal. 2021;40(3):e105784. doi:10.15252/embj.2020105784
    OpenUrlCrossRef
  79. 79.↵
    Siljee JE, Wang Y, Bernard AA, et al. Subcellular localization of MC4R with ADCY3 at neuronal primary cilia underlies a common pathway for genetic predisposition to obesity. Nat Genet. Feb 2018;50(2):180–185. doi:10.1038/s41588-017-0020-9
    OpenUrlCrossRefPubMed
  80. 80.↵
    Timshel PN, Thompson JJ, Pers TH. Genetic mapping of etiologic brain cell types for obesity. eLife. 2020/09/21 2020;9:e55851. doi:10.7554/eLife.55851
    OpenUrlCrossRefPubMed
  81. 81.↵
    Licinio J, Dong C Fau - Wong M-L, Wong ML. Novel sequence variations in the brain-derived neurotrophic factor gene and association with major depression and antidepressant treatment response. (1538–3636 (Electronic))
  82. 82.↵
    Colle R, Deflesselle E, Martin S, et al. BDNF/TRKB/P75NTR polymorphisms and their consequences on antidepressant efficacy in depressed patients. (1744–8042 (Electronic))
Back to top
PreviousNext
Posted February 05, 2024.
Download PDF

Supplementary Material

Data/Code
Email

Thank you for your interest in spreading the word about medRxiv.

NOTE: Your email address is requested solely to identify you as the sender of this article.

Enter multiple addresses on separate lines or separate them with commas.
3D genomic features across >50 diverse cell types reveal insights into the genomic architecture of childhood obesity
(Your Name) has forwarded a page to you from medRxiv
(Your Name) thought you would like to see this page from the medRxiv website.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
3D genomic features across >50 diverse cell types reveal insights into the genomic architecture of childhood obesity
Khanh B. Trang, Matthew C. Pahl, James A. Pippin, Chun Su, Sheridan H. Littleton, Prabhat Sharma, Nikhil N. Kulkarni, Louis R. Ghanem, Natalie A. Terry, Joan M. O’Brien, Yadav Wagley, Kurt D. Hankenson, Ashley Jermusyk, Jason W. Hoskins, Laufey T. Amundadottir, Mai Xu, Kevin M Brown, Stewart A. Anderson, Wenli Yang, Paul M. Titchenell, Patrick Seale, Laura Cook, Megan K. Levings, Babette S. Zemel, Alessandra Chesi, Andrew D. Wells, Struan F.A. Grant
medRxiv 2023.08.30.23294092; doi: https://doi.org/10.1101/2023.08.30.23294092
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools
3D genomic features across >50 diverse cell types reveal insights into the genomic architecture of childhood obesity
Khanh B. Trang, Matthew C. Pahl, James A. Pippin, Chun Su, Sheridan H. Littleton, Prabhat Sharma, Nikhil N. Kulkarni, Louis R. Ghanem, Natalie A. Terry, Joan M. O’Brien, Yadav Wagley, Kurt D. Hankenson, Ashley Jermusyk, Jason W. Hoskins, Laufey T. Amundadottir, Mai Xu, Kevin M Brown, Stewart A. Anderson, Wenli Yang, Paul M. Titchenell, Patrick Seale, Laura Cook, Megan K. Levings, Babette S. Zemel, Alessandra Chesi, Andrew D. Wells, Struan F.A. Grant
medRxiv 2023.08.30.23294092; doi: https://doi.org/10.1101/2023.08.30.23294092

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Subject Area

  • Endocrinology (including Diabetes Mellitus and Metabolic Disease)
Subject Areas
All Articles
  • Addiction Medicine (430)
  • Allergy and Immunology (755)
  • Anesthesia (221)
  • Cardiovascular Medicine (3288)
  • Dentistry and Oral Medicine (364)
  • Dermatology (277)
  • Emergency Medicine (479)
  • Endocrinology (including Diabetes Mellitus and Metabolic Disease) (1169)
  • Epidemiology (13359)
  • Forensic Medicine (19)
  • Gastroenterology (898)
  • Genetic and Genomic Medicine (5147)
  • Geriatric Medicine (481)
  • Health Economics (782)
  • Health Informatics (3264)
  • Health Policy (1140)
  • Health Systems and Quality Improvement (1190)
  • Hematology (429)
  • HIV/AIDS (1017)
  • Infectious Diseases (except HIV/AIDS) (14622)
  • Intensive Care and Critical Care Medicine (912)
  • Medical Education (476)
  • Medical Ethics (127)
  • Nephrology (522)
  • Neurology (4919)
  • Nursing (262)
  • Nutrition (727)
  • Obstetrics and Gynecology (882)
  • Occupational and Environmental Health (795)
  • Oncology (2519)
  • Ophthalmology (723)
  • Orthopedics (280)
  • Otolaryngology (347)
  • Pain Medicine (323)
  • Palliative Medicine (90)
  • Pathology (543)
  • Pediatrics (1299)
  • Pharmacology and Therapeutics (550)
  • Primary Care Research (556)
  • Psychiatry and Clinical Psychology (4205)
  • Public and Global Health (7499)
  • Radiology and Imaging (1704)
  • Rehabilitation Medicine and Physical Therapy (1011)
  • Respiratory Medicine (980)
  • Rheumatology (479)
  • Sexual and Reproductive Health (497)
  • Sports Medicine (424)
  • Surgery (547)
  • Toxicology (72)
  • Transplantation (235)
  • Urology (205)