Both indirect maternal and direct fetal genetic effects reflect the observational relationship between higher birth weight and lower adult bone mass.

BACKGROUND: Birth weight is considered not only to undermine future growth, but also to induce lifelong diseases; the aim of this study is to explore the relationship between birth weight and adult bone mass. METHODS: We performed multivariable regression analyses to assess the association of birth weight with bone parameters measured by dual-energy X-ray absorptiometry (DXA) and by quantitative ultrasound (QUS), independently. We also implemented a systemic Mendelian randomization (MR) analysis to explore the causal association between them with both fetal-specific and maternal-specific instrumental variables. RESULTS: In the observational analyses, we found that higher birth weight could increase the adult bone area (lumbar spine, ß-coefficient= 0.17, P < 2.00 × 10-16; lateral spine, ß-coefficient = 0.02, P = 0.04), decrease bone mineral content-adjusted bone area (BMCadjArea) (lumbar spine, ß-coefficient= - 0.01, P = 2.27 × 10-14; lateral spine, ß-coefficient = - 0.05, P = 0.001), and decrease adult bone mineral density (BMD) (lumbar spine, ß-coefficient = - 0.04, P = 0.007; lateral spine; ß-coefficient = - 0.03, P = 0.02; heel, ß-coefficient = - 0.06, P < 2.00 × 10-16), and we observed that the effect of birth weight on bone size was larger than that on BMC. In MR analyses, the higher fetal-specific genetically determined birth weight was identified to be associated with higher bone area (lumbar spine; ß-coefficient = 0.15, P = 1.26 × 10-6, total hip, ß-coefficient = 0.15, P = 0.005; intertrochanteric area, ß-coefficient = 0.13, P = 0.0009; trochanter area, ß-coefficient = 0.11, P = 0.03) but lower BMD (lumbar spine, ß-coefficient = - 0.10, P = 0.01; lateral spine, ß-coefficient = - 0.12, P = 0.0003, and heel ß-coefficient = - 0.11, P = 3.33 × 10-13). In addition, we found that the higher maternal-specific genetically determined offspring birth weight was associated with lower offspring adult heel BMD (ß-coefficient = - 0.001, P = 0.04). CONCLUSIONS: The observational analyses suggested that higher birth weight was associated with the increased adult bone area but decreased BMD. By leveraging the genetic instrumental variables with maternal- and fetal-specific effects on birth weight, the observed relationship could be reflected by both the direct fetal and indirect maternal genetic effects.

General and abdominal obesity operate differently as influencing factors of fracture risk in old adults.

To infer the causality between obesity and fracture and the difference between general and abdominal obesity, a prospective study was performed in 456,921 participants, and 10,142 participants developed an incident fracture with follow-up period of 7.96 years. A U-shape relationship was observed between BMI and fracture, with the lowest risk of fracture in overweight participants. The obesity individuals had higher fracture risk when BMD was adjusted, and the protective effect of moderate-high BMI on fracture was mostly mediated by bone mineral density (BMD). However, for abdominal obesity, the higher WCadjBMI (linear) and HCadjBMI (J-shape) were found to be related to higher fracture risk, and less than 30% of the effect was mediated by BMD. By leveraging genetic instrumental variables, it provided additional evidences to support the aforementioned findings. In conclusion, keeping moderate-high BMI might be of benefit to old people in terms of fracture risk, whereas abdominal adiposity might increase risk of fracture.

Genomic analyses of 10,376 individuals in the Westlake BioBank for Chinese (WBBC) pilot project.

We initiate the Westlake BioBank for Chinese (WBBC) pilot project with 4,535 whole-genome sequencing (WGS) individuals and 5,841 high-density genotyping individuals, and identify 81.5 million SNPs and INDELs, of which 38.5% are absent in dbSNP Build 151. We provide a population-specific reference panel and an online imputation server ( https://wbbc.westlake.edu.cn/ ) which could yield substantial improvement of imputation performance in Chinese population, especially for low-frequency and rare variants. By analyzing the singleton density of the WGS data, we find selection signatures in SNX29, DNAH1 and WDR1 genes, and the derived alleles of the alcohol metabolism genes (ADH1A and ADH1B) emerge around 7,000 years ago and tend to be more common from 4,000 years ago in East Asia. Genetic evidence supports the corresponding geographical boundaries of the Qinling-Huaihe Line and Nanling Mountains, which separate the Han Chinese into subgroups, and we reveal that North Han was more homogeneous than South Han.

Integrative analysis of genomic and epigenomic data reveal underlying superenhancer-mediated microRNA regulatory network for human bone mineral density.

Bone mineral density (BMD) is a highly heritable complex trait and is a key indicator for diagnosis and treatment for osteoporosis. In the last decade, numerous susceptibility loci for BMD and fracture have been identified by genome-wide association studies (GWAS); however, fine mapping of these loci is challengeable. Here, we proposed a new long-range fine-mapping approach that combined superenhancers (SEs) and microRNAs (miRNAs) data, which were two important factors in control of cell identity and specific differentiation, with the GWAS summary datasets in cell-type-restricted way. Genome-wide SE-based analysis found that the BMD-related variants were significantly enriched in the osteoblast SE regions, indicative of potential long-range effects of such SNPs. With the SNP-mapped SEs (mSEs), 13 accessible long-range mSE-interacted miRNAs (mSE-miRNAs) were identified by integrating osteoblast Hi-C and ATAC-seq data, including three known bone-related miRNAs (miR-132-3p, miR-212-3p and miR-125b-5p). The putative targets of the two newly identified mSE-miRNAs (miR-548aj-3p and miR-190a-3p) were found largely enriched in osteogenic-related pathway and processes, suggesting that these mSE-miRNAs could be functional in the regulation of osteoblast differentiation. Furthermore, we identified 54 genes with the long-range 'mSE-miRNA' approach, and 24 of them were previously reported to be related to skeletal development. Besides, enrichment analysis found that these genes were specifically enriched in the post-transcriptional regulation and bone formation processes. This study provided a new insight into the approach of fine-mapping of GWAS loci. A tool was provided for the genome-wide SE-based analysis and the detection of long-range osteoblast-restricted mSE-miRNAs (https://github.com/Zheng-Lab-Westlake/Osteo-Fine-Mapp-SNP2SE2miRNA).

Comprehensive assessment of the association between genes on JAK-STAT pathway (IFIH1, TYK2, IL-10) and systemic lupus erythematosus: a meta-analysis.

Previous studies have reported that genes relating to JAK-STAT pathway (IFIH1, TYK2 and IL-10) conferred the susceptibility to SLE. In this study, we performed a meta-analysis (including 43 studies) to evaluate the association between IFIH1 (9288 patients and 24,040 controls), TYK2 (4928 patients and 11,536 controls), IL-10 (3623 patients and 4907 controls) polymorphisms and systemic lupus erythematosus (SLE) in a comprehensive way. We found that IFIH1 rs1990760_T allele was associated with risk of SLE in overall population under three models (allelic: P = 2.56 × 10-11, OR 1.135, 95% CI 1.094-1.179, dominant: P = 1.8 × 10-8, OR 1.203, 95% CI 1.128-1.284, recessive: P = 2.6 × 10-7, OR 1.163, 95% CI 1.098-1.231). A strong association had been observed between TYK2 polymorphism rs2304256_C allele and SLE in Europeans (P = 5.82 × 10-5, OR 1.434, 95% CI 1.203-1.710). When coming to overall population, TYK2 rs2304256_C showed a significant association with SLE under recessive model (P = 8.05 × 10-3, OR 1.314, 95% CI 1.074-1.608). However, the other two SNPs (rs12720270, rs280519) of TYK2 were not significant. The results also indicated an association between IL-10 rs1800896_G allele and SLE in Asians under recessive model (P = 4.65 × 10-3, OR 2.623, 95% CI 1.346-5.115), while, IL-10 rs1800896_G had a trend of association with SLE in European population in dominant model (P = 1.21 × 10-2, OR 1.375, 95% CI 1.072-1.764). In addition, we found IL-10 rs1800896 GG homozygote might be associated with increased susceptibility to SLE (GG vs AA, P = 4.65 × 10-3, OR 1.539, 95% CI 1.142-2.072). We concluded that IFIH1 rs1990760_T and TYK2 rs2304256_C alleles were significantly associated with SLE, and IL-10 rs1800896 GG homozygote might have an enhancement effect on SLE risk.