Utilization of geospatial distribution in the measurement of study cohort representativeness.

OBJECTIVE: The ability to apply results from a study to a broader population remains a primary objective in translational science. Distinct from intrinsic elements of scientific rigor, the extrinsic concept of generalization requires there be alignment between a study cohort and population in which results are expected to be applied. Widespread efforts have been made to quantify representativeness of study cohorts. These techniques, however, often consider the study and target cohorts as monolithic collections that can be directly compared. Overlooking known impacts to health from socio-demographic and environmental factors tied to individual's geographical location, and potentially obfuscating misalignment in underrepresented population subgroups. This manuscript introduces several measures to account for geographic information in the assessment of cohort representation. METHODS: Metrics were defined across two themes: First, measures of recruitment, to assess if a study cohort is drawn at an expected rate and in an expected geographical pattern with respect to individuals in a reference cohort. Second, measures of individual characteristics, to assess if the individuals in the study cohort accurately reflect the sociodemographic, clinical, and geographic diversity observed across a reference cohort while accounting for the geospatial proximity of individuals. RESULTS: As an empirical demonstration, methods are applied to an active clinical study examining asthma in Black/African American patients at a US Midwestern pediatric hospital. Results illustrate how areas of over- and under-recruitment can be identified and contextualized in light of study recruitment patterns at an individual-level, highlighting the ability to identify a subset of features for which the study cohort closely resembled the broader population. In addition they provide an opportunity to dive deeper into misalignments, to identify study cohort members that are in some way distinct from the communities for which they are expected to represent. CONCLUSION: Together, these metrics provide a comprehensive spatial assessment of a study cohort with respect to a broader target population. Such an approach offers researchers a toolset by which to target expected generalization of results derived from a given study.

ASCL2 reciprocally controls key trophoblast lineage decisions during hemochorial placenta development.

Invasive trophoblast cells are critical to spiral artery remodeling in hemochorial placentation. Insufficient trophoblast cell invasion and vascular remodeling can lead to pregnancy disorders including preeclampsia, preterm birth, and intrauterine growth restriction. Previous studies in mice identified achaete-scute homolog 2 (ASCL2) as essential to extraembryonic development. We hypothesized that ASCL2 is a critical and conserved regulator of invasive trophoblast cell lineage development. In contrast to the mouse, the rat possesses deep intrauterine trophoblast cell invasion and spiral artery remodeling similar to human placentation. In this study, we investigated invasive/extravillous trophoblast (EVT) cell differentiation using human trophoblast stem (TS) cells and a loss-of-function mutant Ascl2 rat model. ASCL2 transcripts are expressed in the EVT column and junctional zone, which represent tissue sources of invasive trophoblast progenitor cells within human and rat placentation sites, respectively. Differentiation of human TS cells into EVT cells resulted in significant up-regulation of ASCL2 and several other transcripts indicative of EVT cell differentiation. Disruption of ASCL2 impaired EVT cell differentiation, as indicated by cell morphology and transcript profiles. RNA sequencing analysis of ASCL2-deficient trophoblast cells identified both down-regulation of EVT cell-associated transcripts and up-regulation of syncytiotrophoblast-associated transcripts, indicative of dual activating and repressing functions. ASCL2 deficiency in the rat impacted placental morphogenesis, resulting in junctional zone dysgenesis and failed intrauterine trophoblast cell invasion. ASCL2 acts as a critical and conserved regulator of invasive trophoblast cell lineage development and a modulator of the syncytiotrophoblast lineage.

Intersection of regulatory pathways controlling hemostasis and hemochorial placentation.

Hemochorial placentation is characterized by the development of trophoblast cells specialized to interact with the uterine vascular bed. We utilized trophoblast stem (TS) cell and mutant rat models to investigate regulatory mechanisms controlling trophoblast cell development. TS cell differentiation was characterized by acquisition of transcript signatures indicative of an endothelial cell-like phenotype, which was highlighted by the expression of anticoagulation factors including tissue factor pathway inhibitor (TFPI). TFPI localized to invasive endovascular trophoblast cells of the rat placentation site. Disruption of TFPI in rat TS cells interfered with development of the endothelial cell-like endovascular trophoblast cell phenotype. Similarly, TFPI was expressed in human invasive/extravillous trophoblast (EVT) cells situated within first-trimester human placental tissues and following differentiation of human TS cells. TFPI was required for human TS cell differentiation to EVT cells. We next investigated the physiological relevance of TFPI at the placentation site. Genome-edited global TFPI loss-of-function rat models revealed critical roles for TFPI in embryonic development, resulting in homogeneous midgestation lethality prohibiting analysis of the role of TFPI as a regulator of the late-gestation wave of intrauterine trophoblast cell invasion. In vivo trophoblast-specific TFPI knockdown was compatible with pregnancy but had profound effects at the uterine-placental interface, including restriction of the depth of intrauterine trophoblast cell invasion while leading to the accumulation of natural killer cells and increased fibrin deposition. Collectively, the experimentation implicates TFPI as a conserved regulator of invasive/EVT cell development, uterine spiral artery remodeling, and hemostasis at the maternal-fetal interface.

Neonatal hyperoxia induces activated pulmonary cellular states and sex-dependent transcriptomic changes in a model of experimental bronchopulmonary dysplasia.

Hyperoxia disrupts lung development in mice and causes bronchopulmonary dysplasia (BPD) in neonates. To investigate sex-dependent molecular and cellular programming involved in hyperoxia, we surveyed the mouse lung using single cell RNA sequencing (scRNA-seq), and validated our findings in human neonatal lung cells in vitro. Hyperoxia-induced inflammation in alveolar type (AT) 2 cells gave rise to damage-associated transient progenitors (DATPs). It also induced a new subpopulation of AT1 cells with reduced expression of growth factors normally secreted by AT1 cells, but increased mitochondrial gene expression. Female alveolar epithelial cells had less EMT and pulmonary fibrosis signaling in hyperoxia. In the endothelium, expansion of Car4+ EC (Cap2) was seen in hyperoxia along with an emergent subpopulation of Cap2 with repressed VEGF signaling. This regenerative response was increased in females exposed to hyperoxia. Mesenchymal cells had inflammatory signatures in hyperoxia, with a new distal interstitial fibroblast subcluster characterized by repressed lipid biosynthesis and a transcriptomic signature resembling myofibroblasts. Hyperoxia-induced gene expression signatures in human neonatal fibroblasts and alveolar epithelial cells in vitro resembled mouse scRNA-seq data. These findings suggest that neonatal exposure to hyperoxia programs distinct sex-specific stem cell progenitor and cellular reparative responses that underpin lung remodeling in BPD.

Insurance denials and diagnostic rates in a pediatric genomic research cohort.

PURPOSE: This study aimed to assess the amount and types of clinical genetic testing denied by insurance and the rate of diagnostic and candidate genetic findings identified through research in patients who faced insurance denials. METHODS: Analysis consisted of review of insurance denials in 801 patients enrolled in a pediatric genomic research repository with either no previous genetic testing or previous negative genetic testing result identified through cross-referencing with insurance prior-authorizations in patient medical records. Patients and denials were also categorized by type of insurance coverage. Diagnostic findings and candidate genetic findings in these groups were determined through review of our internal variant database and patient charts. RESULTS: Of the 801 patients analyzed, 147 had insurance prior-authorization denials on record (18.3%). Exome sequencing and microarray were the most frequently denied genetic tests. Private insurance was significantly more likely to deny testing than public insurance (odds ratio = 2.03 [95% CI = 1.38-2.99] P = .0003). Of the 147 patients with insurance denials, 53.7% had at least 1 diagnostic or candidate finding and 10.9% specifically had a clinically diagnostic finding. Fifty percent of patients with clinically diagnostic results had immediate medical management changes (5.4% of all patients experiencing denials). CONCLUSION: Many patients face a major barrier to genetic testing in the form of lack of insurance coverage. A number of these patients have clinically diagnostic findings with medical management implications that would not have been identified without access to research testing. These findings support re-evaluation of insurance carriers' coverage policies.

The Impact of Prior Infection and Age on Antibody Persistence After Severe Acute Respiratory Syndrome Coronavirus 2 Messenger RNA Vaccine.

Determining the duration of immunity to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) vaccines is critical for informing the timing of booster immunization. Many genetic and environmental factors could influence both the magnitude and persistence of the antibody response. Here, we showed that SARS-CoV-2 infection before vaccination and age affected the decay of antibody responses to the SARS-CoV-2 messenger RNA vaccine.

Blood DNA methylation at TXNIP and glycemic changes in response to weight-loss diet interventions: the POUNDS lost trial.

BACKGROUND: Thioredoxin Interacting Protein (TXNIP) functions as a master regulator for glucose homeostasis. Hypomethylation at the 5'-cytosine-phosphate-guanine-3' (CpG) site cg19693031 of TXNIP has been consistently related to islet dysfunction, hyperglycemia, and type 2 diabetes. DNA methylation (DNAm) may reveal the missing mechanistic link between obesity and type 2 diabetes. We hypothesize that baseline DNAm level at TXNIP in blood may be associated with glycemic traits and their changes in response to weight-loss diet interventions. METHODS: We included 639 adult participants with overweight or obesity, who participated in a 2-year randomized weight-loss diet intervention. Baseline blood DNAm levels were profiled by high-resolution methylC-capture sequencing. We defined the regional DNAm level of TXNIP as the average methylation level over CpGs within 500 bp of cg19693031. Generalized linear regression models were used for main analyses. RESULTS: We found that higher regional DNAm at TXNIP was significantly correlated with lower fasting glucose, HbA1c, and Homeostatic Model Assessment of Insulin Resistance (HOMA-IR) at baseline (P < 0.05 for all). Significant interactions were observed between dietary protein intake and DNAm on changes in insulin (P-interaction = 0.007) and HOMA-IR (P-interaction = 0.009) at 6 months. In participants with the highest tertile of regional DNAm at TXNIP, average protein (15%) intake was associated with a greater reduction in insulin (ß: -0.14; 95% CI: -0.24, -0.03; P = 0.011) and HOMA-IR (ß: -0.15; 95% CI: -0.26, -0.03; P = 0.014) than high protein (25%) intake, whereas no significant associations were found in those with the lower tertiles (P > 0.05). The interaction was attenuated to be non-significant at 2 years, presumably related to decreasing adherence to the diet intervention. CONCLUSIONS: Our data indicate that higher regional DNAm level at TXNIP was significantly associated with better fasting glucose, HbA1c, and HOMA-IR; and people with higher regional DNAm levels benefited more in insulin and HOMA-IR improvement by taking the average-protein weight-loss diet.

Genomic answers for children: Dynamic analyses of >1000 pediatric rare disease genomes.

PURPOSE: This study aimed to provide comprehensive diagnostic and candidate analyses in a pediatric rare disease cohort through the Genomic Answers for Kids program. METHODS: Extensive analyses of 960 families with suspected genetic disorders included short-read exome sequencing and short-read genome sequencing (srGS); PacBio HiFi long-read genome sequencing (HiFi-GS); variant calling for single nucleotide variants (SNV), structural variant (SV), and repeat variants; and machine-learning variant prioritization. Structured phenotypes, prioritized variants, and pedigrees were stored in PhenoTips database, with data sharing through controlled access the database of Genotypes and Phenotypes. RESULTS: Diagnostic rates ranged from 11% in patients with prior negative genetic testing to 34.5% in naive patients. Incorporating SVs from genome sequencing added up to 13% of new diagnoses in previously unsolved cases. HiFi-GS yielded increased discovery rate with >4-fold more rare coding SVs compared with srGS. Variants and genes of unknown significance remain the most common finding (58% of nondiagnostic cases). CONCLUSION: Computational prioritization is efficient for diagnostic SNVs. Thorough identification of non-SNVs remains challenging and is partly mitigated using HiFi-GS sequencing. Importantly, community research is supported by sharing real-time data to accelerate gene validation and by providing HiFi variant (SNV/SV) resources from >1000 human alleles to facilitate implementation of new sequencing platforms for rare disease diagnoses.

Large-scale analysis of circulating glutamate and adipose gene expression in relation to abdominal obesity.

Circulating levels of the amino acid glutamate are associated with central fat accumulation, yet the pathophysiology of this relationship remains unknown. We aimed to (i) refine and validate the association between circulating glutamate and abdominal obesity in a large twin cohort, and (ii) investigate whether transcriptomic profiles in adipose tissue could provide insight into the biological mechanisms underlying the association. First, in a cohort of 4665 individuals from the TwinsUK resource, we identified individuals with abdominal obesity and compared prevalence of the latter across circulating glutamate quintiles. Second, we used transcriptomic signatures generated from adipose tissue, both subcutaneous and visceral, to investigate associations with circulating glutamate levels. Individuals in the top circulating glutamate quintile had a sevenfold higher prevalence of abdominal obesity compared to those in the bottom quintile. The adipose tissue transcriptomic analyses identified GLUL, encoding Glutamate-Ammonia Ligase, as being associated with circulating glutamate and abdominal obesity, with pronounced signatures in the visceral depot. In conclusion, circulating glutamate is positively associated with the prevalence of abdominal obesity which relates to dysregulated GLUL expression specifically in visceral adipose tissue.

Humoral immune responses during SARS-CoV-2 mRNA vaccine administration in seropositive and seronegative individuals.

BACKGROUND: The global pandemic of coronavirus disease 2019 (COVID-19) is caused by infection with the SARS-CoV-2 virus. Currently, there are three approved vaccines against SARS-CoV-2 in the USA, including two based on messenger RNA (mRNA) technology that has demonstrated high vaccine efficacy. We sought to characterize humoral immune responses, at high resolution, during immunization with the BNT162b2 (Pfizer-BioNTech) vaccine in individuals with or without prior history of natural SARS-CoV-2 infection. METHODS: We determined antibody responses after each dose of the BNT162b2 SARS-CoV-2 vaccine in individuals who had no prior history of SARS-CoV-2 infection (seronegative) and individuals that had previous viral infection 30-60 days prior to first vaccination (seropositive). To do this, we used both an antibody isotype-specific multiplexed bead-based binding assays targeting multiple SARS-CoV-2 viral protein antigens and an assay that identified potential SARS-CoV-2 neutralizing antibody levels. Moreover, we mapped antibody epitope specificity after immunization using SARS-CoV-2 spike protein peptide arrays. RESULTS: Antibody levels were significantly higher after a single dose in seropositive individuals compared to seronegative individuals and were comparable to levels observed in seronegative individuals after two doses. While IgG was boosted by vaccination for both seronegative and seropositive individuals, only seronegative individuals had increased IgA or IgM antibody titers after primary immunization. We identified immunodominant peptides targeted on both SARS-CoV-2 spike S1 and S2 subunits after vaccination. CONCLUSION: These findings demonstrated the antibody responses to SARS-CoV-2 immunization in seropositive and seronegative individuals and provide support for the concept of using prior infection history as a guide for the consideration of future vaccination regimens. Moreover, we identified key epitopes on the SARS-CoV-2 spike protein that are targeted by antibodies after vaccination that could guide future vaccine and immune correlate development.

An epigenome-wide association study of total serum immunoglobulin E concentration.

Immunoglobulin E (IgE) is a central mediator of allergic (atopic) inflammation. Therapies directed against IgE can alleviate hay fever and allergic asthma. Genetic association studies have not yet identified novel therapeutic targets or pathways underlying IgE regulation. We therefore surveyed epigenetic associations between serum IgE concentrations and methylation at loci concentrated in CpG islands genome wide in 95 nuclear pedigrees, using DNA from peripheral blood leukocytes. We validated positive results in additional families and in subjects from the general population. Here we show replicated associations--with a meta-analysis false discovery rate less than 10(-4)--between IgE and low methylation at 36 loci. Genes annotated to these loci encode known eosinophil products, and also implicate phospholipid inflammatory mediators, specific transcription factors and mitochondrial proteins. We confirmed that methylation at these loci differed significantly in isolated eosinophils from subjects with and without asthma and high IgE levels. The top three loci accounted for 13% of IgE variation in the primary subject panel, explaining the tenfold higher variance found compared with that derived from large single-nucleotide polymorphism genome-wide association studies. This study identifies novel therapeutic targets and biomarkers for patient stratification for allergic diseases.

Dominant gut Prevotella copri in gastrectomised non-obese diabetic Goto-Kakizaki rats improves glucose homeostasis through enhanced FXR signalling.

AIMS/HYPOTHESIS: Drug and surgical-based therapies in type 2 diabetes are associated with altered gut microbiota architecture. Here we investigated the role of the gut microbiome in improved glucose homeostasis following bariatric surgery. METHODS: We carried out gut microbiome analyses in gastrectomised (by vertical sleeve gastrectomy [VSG]) rats of the Goto-Kakizaki (GK) non-obese model of spontaneously occurring type 2 diabetes, followed by physiological studies in the GK rat. RESULTS: VSG in the GK rat led to permanent improvement of glucose tolerance associated with minor changes in the gut microbiome, mostly characterised by significant enrichment of caecal Prevotella copri. Gut microbiota enrichment with P. copri in GK rats through permissive antibiotic treatment, inoculation of gut microbiota isolated from gastrectomised GK rats, and direct inoculation of P. copri, resulted in significant improvement of glucose tolerance, independent of changes in body weight. Plasma bile acids were increased in GK rats following inoculation with P. copri and P. copri-enriched microbiota from VSG-treated rats; the inoculated GK rats then showed increased liver glycogen and upregulated expression of Fxr (also known as Nr1h4), Srebf1c, Chrebp (also known as Mlxipl) and Il10 and downregulated expression of Cyp7a1. CONCLUSIONS: Our data underline the impact of intestinal P. copri on improved glucose homeostasis through enhanced bile acid metabolism and farnesoid X receptor (FXR) signalling, which may represent a promising opportunity for novel type 2 diabetes therapeutics.

Global miRNA expression and correlation with mRNA levels in primary human bone cells.

MicroRNAs (miRNAs) are important post-transcriptional regulators that have recently introduced an additional level of intricacy to our understanding of gene regulation. The aim of this study was to investigate miRNA-mRNA interactions that may be relevant for bone metabolism by assessing correlations and interindividual variability in miRNA levels as well as global correlations between miRNA and mRNA levels in a large cohort of primary human osteoblasts (HOBs) obtained during orthopedic surgery in otherwise healthy individuals. We identified differential expression (DE) of 24 miRNAs, and found 9 miRNAs exhibiting DE between males and females. We identified hsa-miR-29b, hsa-miR-30c2, and hsa-miR-125b and their target genes as important modulators of bone metabolism. Further, we used an integrated analysis of global miRNA-mRNA correlations, mRNA-expression profiling, DE, bioinformatics analysis, and functional studies to identify novel target genes for miRNAs with the potential to regulate osteoblast differentiation and extracellular matrix production. Functional studies by overexpression and knockdown of miRNAs showed that, the differentially expressed miRNAs hsa-miR-29b, hsa-miR-30c2, and hsa-miR-125b target genes highly relevant to bone metabolism, e.g., collagen, type I, α1 (COL1A1), osteonectin (SPARC), Runt-related transcription factor 2 (RUNX2), osteocalcin (BGLAP), and frizzled-related protein (FRZB). These miRNAs orchestrate the activities of key regulators of osteoblast differentiation and extracellular matrix proteins by their convergent action on target genes and pathways to control the skeletal gene expression.

Human metabolic individuality in biomedical and pharmaceutical research.

Genome-wide association studies (GWAS) have identified many risk loci for complex diseases, but effect sizes are typically small and information on the underlying biological processes is often lacking. Associations with metabolic traits as functional intermediates can overcome these problems and potentially inform individualized therapy. Here we report a comprehensive analysis of genotype-dependent metabolic phenotypes using a GWAS with non-targeted metabolomics. We identified 37 genetic loci associated with blood metabolite concentrations, of which 25 show effect sizes that are unusually high for GWAS and account for 10-60% differences in metabolite levels per allele copy. Our associations provide new functional insights for many disease-related associations that have been reported in previous studies, including those for cardiovascular and kidney disorders, type 2 diabetes, cancer, gout, venous thromboembolism and Crohn's disease. The study advances our knowledge of the genetic basis of metabolic individuality in humans and generates many new hypotheses for biomedical and pharmaceutical research.

Phenotypic dissection of bone mineral density reveals skeletal site specificity and facilitates the identification of novel loci in the genetic regulation of bone mass attainment.

Heritability of bone mineral density (BMD) varies across skeletal sites, reflecting different relative contributions of genetic and environmental influences. To quantify the degree to which common genetic variants tag and environmental factors influence BMD, at different sites, we estimated the genetic (rg) and residual (re) correlations between BMD measured at the upper limbs (UL-BMD), lower limbs (LL-BMD) and skull (SK-BMD), using total-body DXA scans of ∼ 4,890 participants recruited by the Avon Longitudinal Study of Parents and their Children (ALSPAC). Point estimates of rg indicated that appendicular sites have a greater proportion of shared genetic architecture (LL-/UL-BMD rg = 0.78) between them, than with the skull (UL-/SK-BMD rg = 0.58 and LL-/SK-BMD rg = 0.43). Likewise, the residual correlation between BMD at appendicular sites (r(e) = 0.55) was higher than the residual correlation between SK-BMD and BMD at appendicular sites (r(e) = 0.20-0.24). To explore the basis for the observed differences in rg and re, genome-wide association meta-analyses were performed (n ∼ 9,395), combining data from ALSPAC and the Generation R Study identifying 15 independent signals from 13 loci associated at genome-wide significant level across different skeletal regions. Results suggested that previously identified BMD-associated variants may exert site-specific effects (i.e. differ in the strength of their association and magnitude of effect across different skeletal sites). In particular, variants at CPED1 exerted a larger influence on SK-BMD and UL-BMD when compared to LL-BMD (P = 2.01 × 10(-37)), whilst variants at WNT16 influenced UL-BMD to a greater degree when compared to SK- and LL-BMD (P = 2.31 × 10(-14)). In addition, we report a novel association between RIN3 (previously associated with Paget's disease) and LL-BMD (rs754388: ß = 0.13, SE = 0.02, P = 1.4 × 10(-10)). Our results suggest that BMD at different skeletal sites is under a mixture of shared and specific genetic and environmental influences. Allowing for these differences by performing genome-wide association at different skeletal sites may help uncover new genetic influences on BMD.

Higher chylomicron remnants and LDL particle numbers associate with CD36 SNPs and DNA methylation sites that reduce CD36.

Cluster of differentiation 36 (CD36) variants influence fasting lipids and risk of metabolic syndrome, but their impact on postprandial lipids, an independent risk factor for cardiovascular disease, is unclear. We determined the effects of SNPs within a ∼410 kb region encompassing CD36 and its proximal and distal promoters on chylomicron (CM) remnants and LDL particles at fasting and at 3.5 and 6 h following a high-fat meal (Genetics of Lipid Lowering Drugs and Diet Network study, n = 1,117). Five promoter variants associated with CMs, four with delayed TG clearance and five with LDL particle number. To assess mechanisms underlying the associations, we queried expression quantitative trait loci, DNA methylation, and ChIP-seq datasets for adipose and heart tissues that function in postprandial lipid clearance. Several SNPs that associated with higher serum lipids correlated with lower adipose and heart CD36 mRNA and aligned to active motifs for PPARγ, a major CD36 regulator. The SNPs also associated with DNA methylation sites that related to reduced CD36 mRNA and higher serum lipids, but mixed-model analyses indicated that the SNPs and methylation independently influence CD36 mRNA. The findings support contributions of CD36 SNPs that reduce adipose and heart CD36 RNA expression to inter-individual variability of postprandial lipid metabolism and document changes in CD36 DNA methylation that influence both CD36 expression and lipids.

Global analysis of DNA methylation variation in adipose tissue from twins reveals links to disease-associated variants in distal regulatory elements.

Epigenetic modifications such as DNA methylation play a key role in gene regulation and disease susceptibility. However, little is known about the genome-wide frequency, localization, and function of methylation variation and how it is regulated by genetic and environmental factors. We utilized the Multiple Tissue Human Expression Resource (MuTHER) and generated Illumina 450K adipose methylome data from 648 twins. We found that individual CpGs had low variance and that variability was suppressed in promoters. We noted that DNA methylation variation was highly heritable (h(2)median = 0.34) and that shared environmental effects correlated with metabolic phenotype-associated CpGs. Analysis of methylation quantitative-trait loci (metQTL) revealed that 28% of CpGs were associated with nearby SNPs, and when overlapping them with adipose expression quantitative-trait loci (eQTL) from the same individuals, we found that 6% of the loci played a role in regulating both gene expression and DNA methylation. These associations were bidirectional, but there were pronounced negative associations for promoter CpGs. Integration of metQTL with adipose reference epigenomes and disease associations revealed significant enrichment of metQTL overlapping metabolic-trait or disease loci in enhancers (the strongest effects were for high-density lipoprotein cholesterol and body mass index [BMI]). We followed up with the BMI SNP rs713586, a cg01884057 metQTL that overlaps an enhancer upstream of ADCY3, and used bisulphite sequencing to refine this region. Our results showed widespread population invariability yet sequence dependence on adipose DNA methylation but that incorporating maps of regulatory elements aid in linking CpG variation to gene regulation and disease risk in a tissue-dependent manner.

Genetic determinants of trabecular and cortical volumetric bone mineral densities and bone microstructure.

Most previous genetic epidemiology studies within the field of osteoporosis have focused on the genetics of the complex trait areal bone mineral density (aBMD), not being able to differentiate genetic determinants of cortical volumetric BMD (vBMD), trabecular vBMD, and bone microstructural traits. The objective of this study was to separately identify genetic determinants of these bone traits as analysed by peripheral quantitative computed tomography (pQCT). Separate GWA meta-analyses for cortical and trabecular vBMDs were performed. The cortical vBMD GWA meta-analysis (n = 5,878) followed by replication (n = 1,052) identified genetic variants in four separate loci reaching genome-wide significance (RANKL, rs1021188, p = 3.6×10⁻¹4; LOC285735, rs271170, p = 2.7×10⁻¹²; OPG, rs7839059, p = 1.2×10⁻¹°; and ESR1/C6orf97, rs6909279, p = 1.1×10⁻9). The trabecular vBMD GWA meta-analysis (n = 2,500) followed by replication (n = 1,022) identified one locus reaching genome-wide significance (FMN2/GREM2, rs9287237, p = 1.9×10⁻9). High-resolution pQCT analyses, giving information about bone microstructure, were available in a subset of the GOOD cohort (n = 729). rs1021188 was significantly associated with cortical porosity while rs9287237 was significantly associated with trabecular bone fraction. The genetic variant in the FMN2/GREM2 locus was associated with fracture risk in the MrOS Sweden cohort (HR per extra T allele 0.75, 95% confidence interval 0.60-0.93) and GREM2 expression in human osteoblasts. In conclusion, five genetic loci associated with trabecular or cortical vBMD were identified. Two of these (FMN2/GREM2 and LOC285735) are novel bone-related loci, while the other three have previously been reported to be associated with aBMD. The genetic variants associated with cortical and trabecular bone parameters differed, underscoring the complexity of the genetics of bone parameters. We propose that a genetic variant in the RANKL locus influences cortical vBMD, at least partly, via effects on cortical porosity, and that a genetic variant in the FMN2/GREM2 locus influences GREM2 expression in osteoblasts and thereby trabecular number and thickness as well as fracture risk.

Genetic risk factors for decreased bone mineral accretion in children with asthma receiving multiple oral corticosteroid bursts.

BACKGROUND: Long-term intermittent oral corticosteroid (OCS) use in children with asthma leads to significant decreases in bone mineral accretion (BMA). OBJECTIVE: We aimed to identify genetic factors influencing OCS dose effects on BMA in children with asthma. METHODS: We first performed a gene-by-OCS interaction genome-wide association study (GWAS) of BMA in 489 white participants in the Childhood Asthma Management Program trial who took short-term oral prednisone bursts when they experienced acute asthma exacerbations. We selected the top-ranked 2000 single nucleotide polymorphisms (SNPs) in the GWAS and determined whether these SNPs also had cis-regulatory effects on dexamethasone-induced gene expression in osteoblasts. RESULTS: We identified 2 SNPs (rs9896933 and rs2074439) associated with decreased BMA and related to the tubulin γ pathway. The rs9896933 variant met the criteria for genome-wide significance (P = 3.15 × 10(-8) in the GWAS) and is located on the intron of tubulin folding cofactor D (TBCD) gene. The rs2074439 variant (P = 2.74 × 10(-4) in the GWAS) showed strong cis-regulatory effects on dexamethasone-induced tubulin γ gene expression in osteoblasts (P = 8.64 × 10(-4)). Interestingly, we found that BMA worsened with increasing prednisone dose as the number of mutant alleles of the 2 SNPs increased. CONCLUSIONS: We have identified 2 novel tubulin γ pathway SNPs, rs9896933 and rs2074439, showing independent interactive effects with cumulative corticosteroid dose on BMA in children with asthma receiving multiple OCS bursts.