Dissection of Cellular Communication between Human Primary Osteoblasts and Bone Marrow Mesenchymal Stem Cells in Osteoarthritis at Single-Cell Resolution.

Background and Objectives: Osteoblasts are derived from bone marrow mesenchymal stem cells (BMMSCs) and play important role in bone remodeling. While our previous studies have investigated the cell subtypes and heterogeneity in osteoblasts and BMMSCs separately, cell-to-cell communications between osteoblasts and BMMSCs in vivo in humans have not been characterized. The aim of this study was to investigate the cellular communication between human primary osteoblasts and bone marrow mesenchymal stem cells. Methods and Results: To investigate the cell-to-cell communications between osteoblasts and BMMSCs and identify new cell subtypes, we performed a systematic integration analysis with our single-cell RNA sequencing (scRNA-seq) transcriptomes data from BMMSCs and osteoblasts. We successfully identified a novel preosteoblasts subtype which highly expressed ATF3, CCL2, CXCL2 and IRF1. Biological functional annotations of the transcriptomes suggested that the novel preosteoblasts subtype may inhibit osteoblasts differentiation, maintain cells to a less differentiated status and recruit osteoclasts. Ligand-receptor interaction analysis showed strong interaction between mature osteoblasts and BMMSCs. Meanwhile, we found FZD1 was highly expressed in BMMSCs of osteogenic differentiation direction. WIF1 and SFRP4, which were highly expressed in mature osteoblasts were reported to inhibit osteogenic differentiation. We speculated that WIF1 and sFRP4 expressed in mature osteoblasts inhibited the binding of FZD1 to Wnt ligand in BMMSCs, thereby further inhibiting osteogenic differentiation of BMMSCs. Conclusions: Our study provided a more systematic and comprehensive understanding of the heterogeneity of osteogenic cells. At the single cell level, this study provided insights into the cell-to-cell communications between BMMSCs and osteoblasts and mature osteoblasts may mediate negative feedback regulation of osteogenesis process.

Genome-wide meta-analysis of alcohol use disorder in East Asians.

Alcohol use disorder (AUD) is a leading cause of death and disability worldwide. Genome-wide association studies (GWAS) have identified ~30 AUD risk genes in European populations, but many fewer in East Asians. We conducted GWAS and genome-wide meta-analysis of AUD in 13,551 subjects with East Asian ancestry, using published summary data and newly genotyped data from five cohorts: (1) electronic health record (EHR)-diagnosed AUD in the Million Veteran Program (MVP) sample; (2) DSM-IV diagnosed alcohol dependence (AD) in a Han Chinese-GSA (array) cohort; (3) AD in a Han Chinese-Cyto (array) cohort; and (4) two AD Thai cohorts. The MVP and Thai samples included newly genotyped subjects from ongoing recruitment. In total, 2254 cases and 11,297 controls were analyzed. An AUD polygenic risk score was analyzed in an independent sample with 4464 East Asians (Genetic Epidemiology Research in Adult Health and Aging (GERA)). Phenotypes from survey data and ICD-9-CM diagnoses were tested for association with the AUD PRS. Two risk loci were detected: the well-known functional variant rs1229984 in ADH1B and rs3782886 in BRAP (near the ALDH2 gene locus) are the lead variants. AUD PRS was significantly associated with days per week of alcohol consumption (beta = 0.43, SE = 0.067, p = 2.47 × 10-10) and nominally associated with pack years of smoking (beta = 0.09, SE = 0.05, p = 4.52 × 10-2) and ever vs. never smoking (beta = 0.06, SE = 0.02, p = 1.14 × 10-2). This is the largest GWAS of AUD in East Asians to date. Building on previous findings, we were able to analyze pleiotropy, but did not identify any new risk regions, underscoring the importance of recruiting additional East Asian subjects for alcohol GWAS.

Identification of novel pleiotropic gene for bone mineral density and lean mass using the cFDR method.

Bone mineral density (BMD) and whole-body lean mass (WBLM) are two important phenotypes of osteoporosis and sarcopenia. Previous studies have shown that BMD and lean mass were phenotypically and genetically correlated. To identify the novel common genetic factors shared between BMD and WBLM, we performed the conditional false discovery rate (cFDR) analysis using summary data of the genome-wide association study of femoral neck BMD (n = 53,236) and WBLM (n = 38,292) from the Genetic Factors for Osteoporosis Consortium (GEFOS). We identified eight pleiotropic Single Nucleotide Polymorphism (SNPs) (PLCL1 rs11684176 and rs2880389, JAZF1 rs198, ADAMTSL3 rs10906982, RFTN2/MARS2 rs7340470, SH3GL3 rs1896797, ST7L rs10776755, ANKRD44/SF3B1 rs11888760) significantly associated with femoral neck BMD and WBLM (ccFDR < 0.05). Bayesian fine-mapping analysis showed that rs11888760, rs198, and rs1896797 were the possible functional variants in the ANKRD44/SF3B1, JAZF1i, and SH3GL3 loci, respectively. Functional annotation suggested that rs11888760 was likely to comprise a DNA regulatory element and linked to the expression of RFTN2 and PLCL1. PLCL1 showed differential expression in laryngeal posterior cricoarytenoid muscle between rats of 6 months and 30 months of age. Our findings, together with PLCL1's potential functional relevance to bone and skeletal muscle function, suggested that rs11888760 was the possible pleiotropic functional variants appearing to coregulate both bone and muscle metabolism through regulating the expression of PLCL1. The findings enhanced our knowledge of genetic associations between BMD and lean mass and provide a rationale for subsequent functional studies of the implicated genes in the pathophysiology of diseases, such as osteoporosis and sarcopenia.

Identification of novel variants associated with osteoporosis, type 2 diabetes and potentially pleiotropic loci using pleiotropic cFDR method.

AIMS: Clinical and epidemiological findings point to an association between type 2 diabetes (T2D) and osteoporosis. Genome-wide association studies (GWASs) have been fruitful in identifying some loci potentially associated with osteoporosis and T2D respectively. However, the total genetic variance for each of these two diseases and the shared genetic determination between them are largely unknown. The aim of this study was to identify novel genetic variants for osteoporosis and/or T2D. METHODS: First, using a pleiotropic conditional false discovery rate (cFDR) method, we analyzed two GWAS summary data of femoral neck bone mineral density (FN_BMD, n = 53,236) and T2D (n = 159,208) to identify novel shared genetic loci. FN_BMD is an important risk factor for osteoporosis. Next, to explore the potential functions of the identified potential pleiotropic genes, differential expression analysis was performed for them in monocytes and peripheral blood mononuclear cells (PBMCs) as these cells are relevant to the etiology of osteoporosis and/or T2D. Further, weighted gene co-expression analysis (WGCNA) was conducted to identify functional connections between novel pleiotropic genes and known osteoporosis/T2D susceptibility genes by using transcriptomic expression datasets in bone biopsies (E-MEXP-1618) and pancreatic islets (GSE50397). Finally, multi-trait fine mapping for the detected pleiotropic risk loci were conducted to identify the SNPs that have the highest probability of being causal for both FN_BMD and T2D. RESULTS: We identified 27 significant SNPs with cFDR<0.05 for FN_BMD and 61 SNPs for T2D respectively. Four loci, rs7068487 (PLEKHA1), rs10885421 (TCF7L2), rs944082 (GNG12-AS1 (WLS)) and rs2065929 (PIFO||PGCP1), were found to be potentially pleiotropic and shared between FN_BMD and T2D (ccFDR<0.05). PLEKHA1 was found differentially expressed in circulating monocytes between high and low BMD subjects, and PBMCs between diabetic and non-diabetic conditions. WGCNA showed that PLEKHA1 and TCF7L2 were interconnected with multiple osteoporosis and T2D associated genes in bone biopsy and pancreatic islets, such as JAG, EN1 and CPE. Fine mapping showed that rs11200594 was a potentially causal variant in the locus of PLEKHA1. rs11200594 is also an eQTL of PLEKHA1 in multiple tissue (e.g. peripheral blood cells, adipose and ovary) and is in strong LD with a number of functional variants. CONCLUSIONS: Four potential pleiotropic loci were identified for shared genetic determination of osteoporosis and T2D. Our study highlights PLEKHA1 as an important potentially pleiotropic gene. The findings may help us gain a better understanding of the shared genetic determination between these two important disorders.

Impact of female cigarette smoking on circulating B cells in vivo: the suppressed ICOSLG, TCF3, and VCAM1 gene functional network may inhibit normal cell function.

As pivotal immune guardians, B cells were found to be directly associated with the onset and development of many smoking-induced diseases. However, the in vivo molecular response of B cells underlying the female cigarette smoking remains unknown. Using the genome-wide Affymetrix HG-133A GeneChip microarray, we firstly compared the gene expression profiles of peripheral circulating B cells between 39 smoking and 40 non-smoking healthy US white women. A total of 125 differential expressed genes were identified in our study, and 75.2% of them were down-regulated in smokers. We further obtained genotypes of 702 single nucleotide polymorphisms in those promising genes and assessed their associations with smoking status. Using a novel multicriteria evaluation model integrating information from microarray and the association studies, several genes were further revealed to play important roles in the response of smoking, including ICOSLG (CD275, inducible T-cell co-stimulator ligand), TCF3 (E2A immunoglobulin enhancer binding factors E12/E47), VCAM1 (CD106, vascular cell adhesion molecule 1), CCR1 (CD191, chemokine C-C motif receptor 1) and IL13 (interleukin 13). The differential expression of ICOSLG (p = 0.0130) and TCF3 (p = 0.0125) genes between the two groups were confirmed by real-time reverse transcription PCR experiment. Our findings support the functional importance of the identified genes in response to the smoking stimulus. This is the first in vivo genome-wide expression study on B cells at today's context of high prevalence rate of smoking for women. Our results highlight the potential usage of integrated analyses for unveiling the novel pathogenesis mechanism and emphasized the significance of B cells in the etiology of smoking-induced disease.

Genome-wide association study suggested copy number variation may be associated with body mass index in the Chinese population.

Obesity is a major public health problem characterized with high body mass index (BMI). Copy number variations (CNVs) have been identified to be associated with complex human diseases. The effect of CNVs on obesity is unknown. In this study, we explored the association of CNVs with BMI in 597 Chinese Han subjects using Affymetrix GeneChip Human Mapping 500K Array Set. We found that one CNV at 10q11.22 (from 46.36 Mb to 46.56 Mb) was associated with BMI (the raw P=0.011). The CNV contributed 1.6% of BMI variation, and it covered one important obesity gene-pancreatic polypeptide receptor 1(PPYR1). It was reported that PPYR1 was a key regulator of energy homeostasis. Our findings suggested that CNV might be potentially important for the BMI variation. In addition, our study suggested that CNV might be used as a genetic marker to locate genes associated with BMI in Chinese population.

Genome-wide copy-number-variation study identified a susceptibility gene, UGT2B17, for osteoporosis.

Osteoporosis, a highly heritable disease, is characterized mainly by low bone-mineral density (BMD), poor bone geometry, and/or osteoporotic fractures (OF). Copy-number variation (CNV) has been shown to be associated with complex human diseases. The contribution of CNV to osteoporosis has not been determined yet. We conducted case-control genome-wide CNV analyses, using the Affymetrix 500K Array Set, in 700 elderly Chinese individuals comprising 350 cases with homogeneous hip OF and 350 matched controls. We constructed a genomic map containing 727 CNV regions in Chinese individuals. We found that CNV 4q13.2 was strongly associated with OF (p = 2.0 x 10(-4), Bonferroni-corrected p = 0.02, odds ratio = 1.73). Validation experiments using PCR and electrophoresis, as well as real-time PCR, further identified a deletion variant of UGT2B17 in CNV 4q13.2. Importantly, the association between CNV of UGT2B17 and OF was successfully replicated in an independent Chinese sample containing 399 cases with hip OF and 400 controls. We further examined this CNV's relevance to major risk factors for OF (i.e., hip BMD and femoral-neck bone geometry) in both Chinese (689 subjects) and white (1000 subjects) samples and found consistently significant results (p = 5.0 x 10(-4) -0.021). Because UGT2B17 encodes an enzyme catabolizing steroid hormones, we measured the concentrations of serum testosterone and estradiol for 236 young Chinese males and assessed their UGT2B17 copy number. Subjects without UGT2B17 had significantly higher concentrations of testosterone and estradiol. Our findings suggest the important contribution of CNV of UGT2B17 to the pathogenesis of osteoporosis.

Receiver-operating characteristic analyses of body mass index, waist circumference and waist-to-hip ratio for obesity: Screening in young adults in central south of China.

BACKGROUND & AIMS: Recent data suggest that current obesity diagnostic criterion based on body mass index (BMI) above 30 in Caucasians may not be appropriate for Asian populations. Our aim was to identify the usefulness of BMI, waist circumference (WC) and waist-to-hip ratio (WHR) in screening for obesity in an Asian population. METHODS: A cross-sectional sample of 1109 males and 879 females aged 20-45-yr were recruited. Height, weight, WC, hip circumference and percentage body fat (PBF) were measured in all subjects. Then receiver-operating characteristic analyses were used to evaluate the performances of the three anthropometric indices. RESULTS: BMI, WC and WHR showed strong positive correlation with PBF (r=0.47-0.75) in both males and females within both age groups. True-positive rates ranged from 82.4% to 94.1% and 68.8% to 86.3% in males and females, respectively. True-negative rates ranged from 64.1% to 84.7% and from 56.9% to 79.0%, respectively. The areas under the curves (AUCs) for WC and BMI were high (0.76-0.92) in both sexes and divided age groups (20-30-yr and 31-45-yr), and those for WHR were a little lower (0.74-0.88). CONCLUSIONS: BMI and WC are two important predictors for obesity in Chinese, and WHR is an alternative.

Linkage exclusion mapping with bone size in 79 Caucasian pedigrees.

Bone size is an important risk factor of osteoporotic fractures and has strong genetic determination. However, genetic studies on bone size variation are relatively rare. In the present study, we conducted a linkage exclusion mapping for bone size variation on chromosomes 1, 4, 6, and 17 in 79 Caucasian pedigrees. For hip bone size variation, several genomic regions were excluded at effect sizes of 10% or greater, including regions of 61-77cM at 1p35-p34, 43-59cM at 4p15-p13, 1-59cM at 6p25-p21, 82-88cM at 17q23-q24, and 113-114cM at 17q25. For spine bone size, at effect sizes of 10% or greater, we excluded regions of 115-122cM at 1p31-p22, 136-141cM at 1p21, 207-260cM at 1q31-q42, 20-89cM at 4p16-q21, 11-21cM at 6p24-p23, and 1-6cM at 17p13. These results suggested that a number of candidate genes located in the excluded regions, such as the growth hormone (GH) gene, tumor necrosis factor-alpha (TNF-alpha) gene, and bone morphogenetic protein-3 (BMP3) gene, are unlikely to have a substantial effect on bone size variation in this Caucasian population.

Association and haplotype analyses of the COL1A2 and ER-alpha gene polymorphisms with bone size and height in Chinese.

Bone size (BS) is another risk factor of fracture independent of BMD in determining bone strength, and height is highly related with BS. To test the effect of the estrogen receptor-alpha (ER-alpha) and collagen type I alpha 2 (COL1A2) genes on the variation of BS and height, we genotyped the PvuII and XbaI polymorphisms in the intron 1 of the ER-alpha gene and the MspI and (GT)n markers in the intron 47 and intron 1 of the COL1A2 gene in 400 Chinese nuclear families with a total of 1256 individuals. The BS at the hip and spine was measured using a Hologic QDR 2000 dual-energy X-ray absorptiometry (DXA) scanner. Population stratification, total-family association, and within-family association were used to test the relationship of BS (at the spine and hip) and height with the four polymorphisms. We also performed these association analyses with the haplotypes of the MspI and (GT)n polymorphisms in the COL1A2 gene, and with the haplotypes of the PvuII and XbaI markers in the ER-alpha gene. Weak within-family association was found between the COL1A2-MspI (P = 0.05) and the femoral neck BS, between the ER-alpha-PX (P = 0.04) and the intertrochanter BS, and between the COL1A2-(GT)(17) (P = 0.02), COL1A2-m(GT)(17) (P = 0.009) and height. Subsequent permutation tests generally confirmed the suggestive within-family association. For the weak within-family association, the proportions of phenotypic variance accounted by the COL1A2-MspI, ER-alpha-PX, COL1A2-(GT)(17), COL1A2-m(GT)(17) markers were 1.50%, 1.51%, 2.15%, and 2.43% for the corresponding phenotypes. The association results indicate that the (GT)n and MspI markers of COL1A2 gene may have some influence on the variation of both BS and height, and the XbaI and PvuII markers of ER-alpha gene may have some effect on the variation of height in Chinese but not on the variation of BS.

Polymorphisms of four bone mineral density candidate genes in Chinese populations and comparison with other populations of different ethnicity.

Studies on polymorphisms of candidate genes and their association with bone mineral density (BMD) have been reported in many populations, but few have been reported in Chinese populations. We investigated polymorphisms of the following five commonly used markers of four prominent BMD candidate genes with the purpose of identifying useful genetic markers for osteoporosis genetic research in Chinese: the Sp1 and RsaI polymorphisms of the collagen type 1 alpha l (Col1a1) gene, the -174G/C promoter polymorphism of the interleukin 6 (IL-6) gene, the Asn363Ser polymorphism of the glucocorticoid receptor (GR) gene, and the T --> C polymorphism in intron 5 of the transforming growth factor beta(1) (TGF-beta(1)) gene. We evaluated these polymorphisms using PCR-RFLP in samples of at least 124 random individuals. We compared the polymorphisms of these five markers with other populations using the chi(2) test and Fisher's exact two-tailed test. For the RsaI polymorphism, only three heterozygotes but no variant homozygote were identified. For the -174G/C polymorphic site, only one GC heterozygote and no CC homozygote were found. Alleles s, Ser, and A(1) at the Sp1, Asn363Ser, and T --> C marker sites that have been found to be polymorphic in other populations were not found in Chinese. Significant differences of allele and genotype frequency distributions were observed at these polymorphisms ( P < 0.001) after comparing with other populations. Our results suggest that variant alleles of the five markers are absent or too rare to be useful genetic makers in Chinese, despite the fact that they have been commonly used as polymorphic markers in osteoporosis genetic research in other populations.