Mendelian randomization and transcriptome analysis identified immune-related biomarkers for osteoarthritis.

Background: The immune microenvironment assumes a significant role in the pathogenesis of osteoarthritis (OA). However, the current biomarkers for the diagnosis and treatment of OA are not satisfactory. Our study aims to identify new OA immune-related biomarkers to direct the prevention and treatment of OA using multi-omics data. Methods: The discovery dataset integrated the GSE89408 and GSE143514 datasets to identify biomarkers that were significantly associated with the OA immune microenvironment through multiple machine learning methods and weighted gene co-expression network analysis (WGCNA). The identified signature genes were confirmed using two independent validation datasets. We also performed a two-sample mendelian randomization (MR) study to generate causal relationships between biomarkers and OA using OA genome-wide association study (GWAS) summary data (cases n = 24,955, controls n = 378,169). Inverse-variance weighting (IVW) method was used as the main method of causal estimates. Sensitivity analyses were performed to assess the robustness and reliability of the IVW results. Results: Three signature genes (FCER1G, HLA-DMB, and HHLA-DPA1) associated with the OA immune microenvironment were identified as having good diagnostic performances, which can be used as biomarkers. MR results showed increased levels of FCER1G (OR = 1.118, 95% CI 1.031-1.212, P = 0.041), HLA-DMB (OR = 1.057, 95% CI 1.045 -1.069, P = 1.11E-21) and HLA-DPA1 (OR = 1.030, 95% CI 1.005-1.056, P = 0.017) were causally and positively associated with the risk of developing OA. Conclusion: The present study identified the 3 potential immune-related biomarkers for OA, providing new perspectives for the prevention and treatment of OA. The MR study provides genetic support for the causal effects of the 3 biomarkers with OA and may provide new insights into the molecular mechanisms leading to the development of OA.

Investigating the impact of human blood metabolites on the Sepsis development and progression: a study utilizing two-sample Mendelian randomization.

Background: Existing data suggests a potential link between human blood metabolites and sepsis, yet the precise cause-and-effect relationship remains elusive. By using a two-sample Mendelian randomization (MR) analysis, this study aims to establish a causal link between human blood metabolites and sepsis. Methods: A two-sample MR analysis was employed to investigate the relationship between blood metabolites and sepsis. To assess the causal connection between sepsis and human blood metabolites, five different MR methods were employed, A variety of sensitivity analyses were conducted, including Cochrane's Q test, MR-Egger intercept test, MR-PRESSO and leave-one-out (LOO) analysis. In order to ensure the robustness of the causal association between exposure and outcome, the Bonferroni adjustment was employed. Additionally, we conducted analyses of the metabolic pathways of the identified metabolites using the Kyoto Encyclopedia of Genes and Genomes (KEGG) and the Small Molecule Pathway Database (SMPDB) database. Results: The MR analysis revealed a total of 27 metabolites (16 known and 11 unknown) causally linked to the development and progression of sepsis. After applying the Bonferroni correction, 3-carboxy-4-methyl-5-propyl-2-furanpropanoate (CMPF) remained significant in relation to 28-day all-cause mortality in sepsis. By pathway enrichment analysis, we identified four significant metabolic pathways. Notably, the Alpha Linolenic Acid and Linoleic Acid metabolism pathway emerged as a pivotal contributor to the occurrence and progression of sepsis. Conclusion: This study provides preliminary evidence of causal associations between human blood metabolites and sepsis, as ascertained by MR analysis. The findings offer valuable insights into the pathogenesis of sepsis and may provide insight into preventive and therapeutic approaches.

A Mendelian randomization study for drug repurposing reveals bezafibrate and fenofibric acid as potential osteoporosis treatments.

Background: Lipid pathways have been implicated in the pathogenesis of osteoporosis (OP). Lipid-lowering drugs may be used to prevent and treat OP. However, the causal interpretation of results from traditional observational designs is controversial by confounding. We aimed to investigate the causal association between genetically proxied lipid-lowering drugs and OP risk. Methods: We conducted two-step Mendelian randomization (MR) analyses to investigate the causal association of genetically proxied lipid-lowering drugs on the risk of OP. The first step MR was used to estimate the associations of drug target genes expression with low-density lipoprotein cholesterol (LDL-C) levels. The significant SNPs in the first step MR were used as instrumental variables in the second step MR to estimate the associations of LDL-C levels with forearm bone mineral density (FA-BMD), femoral neck BMD (FN-BMD), lumbar spine BMD (LS-BMD) and fracture. The significant lipid-lowering drugs after MR analyses were further evaluated for their effects on bone mineralization using a dexamethasone-induced OP zebrafish model. Results: The first step MR analysis found that the higher expression of four genes (HMGCR, NPC1L1, PCSK9 and PPARG) was significantly associated with a lower LDL-C level. The genetically decreased LDL-C level mediated by the PPARG was significantly associated with increased FN-BMD (BETA = -1.38, p = 0.001) and LS-BMD (BETA = -2.07, p = 3.35 × 10-5) and was marginally significantly associated with FA-BMD (BETA = -2.36, p = 0.008) and reduced fracture risk (OR = 3.47, p = 0.008). Bezafibrate (BZF) and Fenofibric acid (FBA) act as PPARG agonists. Therefore genetically proxied BZF and FBA had significant protective effects on OP. The dexamethasone-induced OP zebrafish treated with BZF and FBA showed increased bone mineralization area and integrated optical density (IOD) with alizarin red staining. Conclusion: The present study provided evidence that BZF and FBA can increase BMD, suggesting their potential effects in preventing and treating OP. These findings potentially pave the way for future studies that may allow personalized selection of lipid-lowering drugs for those at risk of OP.

Functional Enrichment Analysis Identifying Regulatory Information Associated with Human Fracture.

Dozens of loci associated with fracture have been identified by genome-wide association studies (GWASs). However, most of these variants are located in the noncoding regions including introns, long terminal repeats, and intergenic regions. Although combining regulation information helps to identify the causal SNPs and interpret the involvement of these variants in the etiology of human fracture, regulation information which was truly associated with fracture was unknown. A novel functional enrichment method GARFIELD (GWAS Analysis of Regulatory of Functional Information Enrichment with LD correction) was applied to identify fracture-associated regulation information, including transcript factor binding sites, expression quantitative trait loci (eQTLs), chromatin states, enhancer, promoter, dyadic, super enhancer and Epigenome marks. Fracture SNPs were significantly enriched in exon (Bonferroni correction, p value < 7.14 × 10-3) at two GWAS p value thresholds through GARFIELD. High level of fold-enrichment was observed in super enhancer of monocyte and the enhancer of chondrocyte (Bonferroni correction, p value < 4.45 × 10-3). eQTLs of 44 tissues/cells and 10 transcription factors (TFs) were identified to be associated with human fracture. These results provide new insight into the etiology of human fracture, which might increase the identification of the causal SNPs through the fine-mapping study combined with functional annotation, as well as polygenic risk score.

[Construction of fluorescent transgenic zebrafish Tg (ttn.2: EGFP)].

In order to develop a transgenic zebrafish line with green fluorescent protein (enhanced green fluorescent protein, EGFP) expressed specifically in muscle and heart, the recombinant expression vector constructed using the zebrafish ttn.2 gene promoter fragment and EGFP gene coding sequence and the capped mRNA of Tol2 transposase were co-injected into the zebrafish 1-cell stage embryos. The stable genetic Tg (ttn.2: EGFP) transgenic zebrafish line was successfully developed by fluorescence detection, followed by genetic hybridization screening and molecular identification. Fluorescence signals and whole-mount in situ hybridization showed that EGFP expression was located in muscle and heart, the specificity of which was consistent with the expression of ttn.2 mRNA. Inverse PCR showed that EGFP was integrated into chromosomes 4 and 11 of zebrafish in No. 33 transgenic line, while integrated into chromosome 1 in No. 34 transgenic line. The successful construction of this fluorescent transgenic zebrafish line, Tg (ttn.2: EGFP), laid a foundation for the research of muscle and heart development and related diseases. In addition, the transgenic zebrafish lines with strong green fluorescence can also be used as a new ornamental fish.

Transcriptomic analysis reveals shared gene signatures and molecular mechanisms between obesity and periodontitis.

Background: Both obesity (OB) and periodontitis (PD) are chronic non-communicable diseases, and numerous epidemiological studies have demonstrated the association between these two diseases. However, the molecular mechanisms that could explain the association between OB and PD are largely unclear. This study aims to investigate the common gene signatures and biological pathways in OB and PD through bioinformatics analysis of publicly available transcriptome datasets. Methods: The RNA expression profile datasets of OB (GSE104815) and PD (GSE106090) were used as training data, and GSE152991 and GSE16134 as validation data. After screening for differentially expressed genes (DEGs) shared by OB and PD, gene enrichment analysis, protein-protein interaction (PPI) network construction, GeneMANIA analysis, immune infiltration analysis and gene set enrichment analysis (GSEA) were performed. In addition, receiver operating characteristic (ROC) curves were used to assess the predictive accuracy of the hub gene. Finally, we constructed the hub gene-associated TF-miRNA-mRNA regulatory network. Results: We identified a total of 147 DEGs shared by OB and PD (38 down-regulated and 109 up-regulated). Functional analysis showed that these genes were mainly enriched in immune-related pathways such as B cell receptor signalling, leukocyte migration and cellular defence responses. 14 hub genes (FGR, MNDA, NCF2, FYB1, EVI2B, LY86, IGSF6, CTSS, CXCR4, LCK, FCN1, CXCL2, P2RY13, MMP7) showed high sensitivity and specificity in the ROC curve analysis. The results of immune infiltration analysis showed that immune cells such as macrophages, activated CD4 T cells and immune B cells were present at high infiltration levels in both OB and PD samples.The results of GeneMANIA analysis and GSEA analysis suggested that five key genes (FGR, LCK, FYB1, LY86 and P2RY13) may be strongly associated with macrophages. Finally, we constructed a TF-miRNA-mRNA regulatory network consisting of 233 transcription factors (TFs), 8 miRNAs and 14 mRNAs based on the validated information obtained from the database. Conclusions: Five key genes (FGR, LCK, FYB1, LY86, P2RY13) may be important biomarkers of OB and PD. These genes may play an important role in the pathogenesis of OB and PD by affecting macrophage activity and participating in immune regulation and inflammatory responses.

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.

Roles of miR-196a and miR-196b in Zebrafish Motor Function.

BACKGROUND: The exertion of motor function depends on various tissues, such as bones and muscles. miR-196 has been widely studied in cancer and other fields, but its effect on bone and skeletal muscle is rarely reported. In order to explore the role of miR-196 family in bone and skeletal muscle, we used the previously successfully constructed miR-196a-1 and miR-196b gene knockout zebrafish animal models for research. METHODS: The behavioral trajectories of zebrafish from 4 days post-fertilization (dpf) to 7 dpf were detected to analyze the effect of miR-196a-1 and miR-196b on motor ability. Hematoxylin-eosin (HE) staining and transmission electron microscopy (TEM) were used to detect the dorsal muscle tissue of zebrafish. The bone tissue of zebrafish was detected by microcomputed tomography (micro-CT). Real-time PCR was used to detect the expression levels of related genes, including vcp, dpm1, acta1b, mylpfb, col1a1a, bmp8a, gdf6a, and fgfr3. RESULTS: The behavioral test showed that the total behavioral trajectory, movement time, and movement speed of zebrafish larvae were decreased in the miR-196a-1 and miR-196b gene knockout lines. Muscle tissue analysis showed that the structure of muscle fibers in the zebrafish lacking miR-196a-1 and miR-196b was abnormal and was characterized by vacuolar degeneration of muscle fibers, intranuclear migration, melanin deposition, and inflammatory cell infiltration. Bone CT examination revealed decreased bone mineral density and trabecular bone number. The real-time PCR results showed that the expression levels of vcp, dpm1, gdf6a, fgfr3, and col1a1a were decreased in the miR-196b gene knockout group. The expression levels of dpm1, acta1b, mylpfb, gdf6a, and col1a1a were decreased, and the expression level of fgfr3 was increased in the miR-196b gene knockout group compared with the wild-type group. CONCLUSIONS: miR-196a-1 and miR-196b play an important role in muscle fiber structure, bone mineral density, and bone trabecular quantity by affecting the expression of vcp, dpm1, acta1b, mylpfb, gdf6a, fgfr3, and col1a1a and then affect the function of the motor system.

Imputation of Human Primary Osteoblast Single Cell RNA-Seq Data Identified Three Novel Osteoblastic Subtypes.

BACKGROUND: Recently, single-cell RNA sequencing (scRNA-seq) technology was increasingly used to study transcriptomics at a single-cell resolution, scRNA-seq analysis was complicated by the "dropout", where the data only captures a small fraction of the transcriptome. This phenomenon can lead to the fact that the actual expressed transcript may not be detected. We previously performed osteoblast subtypes classification and dissection on freshly isolated human osteoblasts. MATERIALS AND METHODS: Here, we used the scImpute method to impute the missing values of dropout genes from a scRNA-seq dataset generated on freshly isolated human osteoblasts. RESULTS: Based on the imputed gene expression patterns, we discovered three new osteoblast subtypes. Specifically, these newfound osteoblast subtypes are osteoblast progenitors, and two undetermined osteoblasts. Osteoblast progenitors showed significantly high expression of proliferation related genes (FOS, JUN, JUNB and JUND). Analysis of each subtype showed that in addition to bone formation, these undetermined osteoblasts may involve osteoclast and adipocyte differentiation and have the potential function of regulate immune activation. CONCLUSIONS: Our findings provided a new perspective for studying the osteoblast heterogeneity and potential biological functions of these freshly isolated human osteoblasts at the single-cell level, which provides further insight into osteoblasts subtypes under various (pathological) physiological conditions.

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.

Loss of Wnt16 Leads to Skeletal Deformities and Downregulation of Bone Developmental Pathway in Zebrafish.

Wingless-type MMTV integration site family, member 16 (wnt16), is a wnt ligand that participates in the regulation of vertebrate skeletal development. Studies have shown that wnt16 can regulate bone metabolism, but its molecular mechanism remains largely undefined. We obtained the wnt16-/- zebrafish model using the CRISPR-Cas9-mediated gene knockout screen with 11 bp deletion in wnt16, which led to the premature termination of amino acid translation and significantly reduced wnt16 expression, thus obtaining the wnt16-/- zebrafish model. The expression of wnt16 in bone-related parts was detected via in situ hybridization. The head, spine, and tail exhibited significant deformities, and the bone mineral density and trabecular bone decreased in wnt16-/- using light microscopy and micro-CT analysis. RNA sequencing was performed to explore the differentially expressed genes (DEGs). Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis found that the down-regulated DEGs are mainly concentrated in mTOR, FoxO, and VEGF pathways. Protein-protein interaction (PPI) network analysis was performed with the detected DEGs. Eight down-regulated DEGs including akt1, bnip4, ptena, vegfaa, twsg1b, prkab1a, prkab1b, and pla2g4f.2 were validated by qRT-PCR and the results were consistent with the RNA-seq data. Overall, our work provides key insights into the influence of wnt16 gene on skeletal development.

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.

Genistein inhibits lung cancer cell stem-like characteristics by modulating MnSOD and FoxM1 expression.

Manganese superoxide dismutase (MnSOD) promotes invasive and migratory activities by upregulating Forkhead box protein M1 (FoxM1) expression. The present study investigated whether modulation of MnSOD and FoxM1 expression was responsible for the antitumor effects of genistein on cancer stem-like cells (CSLCs) derived from non-small cell lung cancer cells (NSCLCs). Spheroids prepared from H460 or A549 cells were defined as lung cancer stem-like cells (LCSLCs) and were treated with genistein. The Cell Counting Kit-8 assay was performed to assess human lung fibroblast IMR-90 cell proliferation, as well as NSCLC H460 and A549 cell proliferation following treatment with genistein. MnSOD, FoxM1, cluster of differentiation (CD)133, CD44, BMI1 proto-oncogene, polycomb ring finger (Bmi1) and Nanog homeobox (Nanog) protein expression levels were examined via western blotting. The sphere formation assay was conducted to evaluate LCSLC self-renewal potential, and LSCLC migratory and invasive activities were analyzed using the wound healing and Transwell invasion assays, respectively. Knockdown and overexpression of MnSOD and FOXM1 via short hairpin-RNA or cDNA transfection were performed. The results indicated that genistein (80 and 100 µM) suppressed H460 and A549 cell viability compared with IMR-90 cells. Sub-cytotoxic concentrations of genistein (20 and 40 µM) inhibited sphere formation activity and decreased the protein expression levels of CD133, CD44, Bmi1 and Nanog in LCSLCs compared with the control group. Genistein also suppressed the migratory and invasive activities of LCSLCs compared with the control group. MnSOD and FoxM1 overexpression antagonized the effects of genistein (40 µM), whereas MnSOD and FoxM1 knockdown enhanced the inhibitory effects of genistein (20 µM) on CSLC characteristics of LCSLCs. Overall, the results suggested that genistein suppressed lung cancer cell CSLC characteristics by modulating MnSOD and FoxM1 expression levels.

The DNMT1/miR-34a/FOXM1 Axis Contributes to Stemness of Liver Cancer Cells.

BACKGROUND: Whether DNA methyltransferase 1 (DNMT1)/miR-34a/FoxM1 signaling promotes the stemness of liver cancer stem cells (LCSCs) remains unclear. This study aimed to assess whether methylation-based silencing of miR-34a by DNMT1 contributes to stemness features via FoxM1 upregulation in LCSCs. METHODS: The CD133+ subgroup of MHCC97H cells sorted by MACS was used as LCSCs. DNMT1, BMI1, SOX2, and OCT4 mRNA levels, and miR-34a amounts were determined by qRT-PCR. DNMT1, CD44, and FoxM1 proteins were analyzed by immunoblot. Sphere and colony formation abilities were detected by respective assays. CD133+ cell percentages were assessed by flow cytometry. In vivo oncogenicity was evaluated using a tumor xenograft model in mice. The effects of DNMT1/miR-34a signaling on the stemness of LCSCs were examined by knockdown or overexpression of DNMT1 and/or transfection of miR-34a mimic or inhibitor using lentivirus-delivery systems. FoxM1 association with miR-34a was detected by a reporter assay. RESULTS: We here showed that LCSCs exhibited elevated DNMT1 activity and expression, lower miR-34a expression with higher promoter methylation, and stronger stemness, compared with the parental liver cancer cells. DNMT1 knockdown repressed DNMT1, increased miR-34a amounts by promoter demethylation, and reduced stemness in LCSCs, whereas DNMT1 overexpression had the opposite effects in liver cancer cells. Transfection with miR-34a mimic repressed the stemness of LCSCs, while miR-34a inhibitor significantly downregulated miR-34a and enhanced stemness, without affecting DNMT1 in liver cancer cells. MiR-34a mimic rescued the effects of DNMT1 overexpression on the stemness of LCSCs, without affecting DNMT1 expression. Finally, FOXM1 was identified as a direct target by miR-34a in LCSCs. CONCLUSIONS: We revealed that aberrant activation of DNMT1 causes miR-34a promoter methylation and suppression, leading to FoxM1 upregulation by disinhibition and promotion of LCSC stemness. These findings suggest that blockage of DNMT1/miR-34a-mediated FOXM1 upregulation might suppress liver cancer by targeting LCSCs.

A trans-ethnic two-stage polygenetic scoring analysis detects genetic correlation between osteoporosis and schizophrenia.

BACKGROUNDS: To explore the genetic correlation between schizophrenia (SCZ) and osteoporosis (OP). DESIGN, SETTING, PARTICIPANTS, MEASUREMENTS: We conducted a trans-ethnic two-stage genetic correlation analysis of OP and SCZ, totally invoking 2286 Caucasia subjects in discovery stage and 4124 Chinese subjects in replication stage. The bone mineral density (BMD) and bone area values of ulna & radius, hip and spine were measured using Hologic 4500W dual energy X-ray absorptiometry machine. SCZ was diagnosed according to DSM-IV criteria. For the genome-wide association study (GWAS) of Caucasian OP, Chinese OP and Chinese SCZ, SNP genotyping was performed using Affymetrix SNP 6.0 array. For the GWAS of Caucasian SCZ, SNP genotyping was conducted using the Affymetrix 5.0 array, Affymetrix 6.0 array and Illumina 550 K array. Polygenetic risk scoring (PRS) analysis was conducted by PRSice software. Also, Linkage disequilibrium score regression (LD Score regression) analysis was performed to evaluate the genetic correlation between OP and SCZ. Multi-trait analysis of GWAS (MTAG) was performed to detect novel candidate genes for osteoporosis and SCZ. RESULTS: In the Caucasia discovery samples, significant genetic correlations were observed for ulna & radius BMD vs. SCZ (P value = 0.010), ulna & radius area vs. SCZ (P value = 0.031). In the Chinese replication samples, we observed significant correlation for ulna & radius area vs. SCZ (P value = 0.019). In addition, LD Score regression also identified significant genetic correlations between SCZ and bone phenotypes in Caucasian and Chinese sample respectively. MTAG analysis identified several novel candidate genes, such as CTNNA2 (MTAG P value = 2.24 × 10-6) for SCZ and FADS2 (MTAG P value = 2.66 × 10-7) for osteoporosis. CONCLUSIONS: Our study results support the overlapped genetic basis for osteoporosis and SCZ, and provide novel clues for elucidating the biological mechanism of increased osteoporosis risk in SCZ patients.

Isovitexin reduces carcinogenicity and stemness in hepatic carcinoma stem-like cells by modulating MnSOD and FoxM1.

BACKGROUND: Manganese superoxide dismutase (MnSOD) upregulating FoxM1 have previously been demonstrated promoting lung cancer stemness. Isovitexin exhibits antitumor activities in various cancers. This study aimed to assess whether isovitexin inhibits hepatic carcinoma stem-like cells (HCSLCs) features via regulating MnSOD and FoxM1 expression. METHODS: Second-generation spheres from the hepatic carcinoma cell lines, respectively, were used as HCSLCs. Protein amounts of MnSOD, FoxM1 and stemness-associated markers (CD133, CD44, ALDH1, Bmi1, Nanog and Oct4) were determined by immunoblotting. In vitro carcinogenicity was evaluated by sphere- and colony-formation assays. The effects of isovitexin on HCSLC carcinogenicity and stemness were examined in vitro and in xenograft models. An adenoviral delivery system was employed to manipulate MnSOD and/or FoxM1. Luciferase reporter assay was performed to verify isovitexin downregulated FoxM1 by inhibiting MnSOD-mediated effects of E2F1 and/or Sp1 on activation of FoxM1 promoter. RESULTS: FoxM1 upregulation by MnSOD contributed to carcinogenicity and stemness, with increased sphere- and colony-formation capabilities, upregulated stemness-associated markers and CD133+ subpopulation as well as elevated oncogenicity in vivo in HCSLCs compared with hepatic carcinoma cells. Isovitexin substantially decreased sphere and colony formation rates, and stemness-associated markers in cultured HCSLCs by suppressing MnSOD and FoxM1 expression. Importantly, isovitexin significantly inhibited tumor growth of in nude mice bearing HCSLCs and reduced CD133 protein expression of xenograft in nude mice. MnSOD or FoxM1 knockdown enhanced the effects of isovitexin suppression on carcinogenicity and stemness in HCSLC. MnSOD or FoxM1 overexpression attenuated the effects of isovitexin. Additionally, isovitexin and MnSOD knockdown could inhibit FoxM1 reporter activity via a decreased binding of E2F1 and/or Sp1 onto FoxM1 promoter. FoxM1 overexpression reversed the effects of isovitexin combined with MnSOD knockdown, without affecting MnSOD expression. Moreover, MnSOD knockdown plus thiostrepton, a FoxM1 specific inhibitor, cooperated with isovitexin to repress xenograft tumor growth and downregulate MnSOD and FoxM1 in nude mice bearing HCSLCs from MHCC97H cells. CONCLUSIONS: Isovitexin inhibits carcinogenicity and stemness in HCSLCs by downregulating FoxM1via inhibition of MnSOD.

MAVS of triploid hybrid of red crucian carp and allotetraploid possesses the improved antiviral activity compared with the counterparts of its parents.

Triploid hybrid (3n = 150) of red crucian carp (♀, 2n = 100) and allotetraploid (♂, 4n = 200) presents the obviously stronger disease resistance than its parents. To elucidate the innate immunity of triploid hybrid, the MAVS homologues of triploid hybrid (3nMAVS), red crucian carp (2nMAVS) and allotetraploid (4nMAVS) have been identified and characterized separately in this study. 2nMAVS and 4nMAVS were evolutionarily conserved; however, 3nMAVS showed lower amino acid similarity and differently predicted structure to 2nMAVS or 4nMAVS. 3nMAVS transcription increase rate in host cells were obviously higher than 2nMAVS or 4nMAVS in response to different stimuli, which included spring viraemia of carp virus (SVCV), grass carp reovirus (GCRV) and poly (I:C). The reporter assay in EPC cells showed that 3nMAVS owned much stronger ability to induce the production of DrIFNφ1 and eIFN than either 2nMAVS or 4nMAVS. Accordingly, EPC cells transfected with 3nMAVS presented obviously stronger antiviral activity against both GCRV and SVCV than the cells expressing 2nMAVS or 4nMAVS. All the data support the conclusion that 3nMAVS-mediated antiviral signaling during innate immune activation was stronger than those of 2nMAVS and 4nMAVS, which provided us the new insight on the innate immune system of triploid hybrid.

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.

Integration of summary data from GWAS and eQTL studies identified novel causal BMD genes with functional predictions.

PURPOSE: Osteoporosis is a common global health problem characterized by low bone mineral density (BMD) and increased risk of fracture. Genome-wide association studies (GWAS) have identified >100 genetic loci associated with BMD. However, the functional genes responsible for most associations remain largely unknown. We conducted an innovative summary statistic data-based Mendelian randomization (SMR) analysis to identify novel causal genes associated with BMD and explored their potential functional significance. METHODS: After quality control of the largest GWAS meta-analysis data of BMD and the largest expression quantitative trait loci (eQTL) meta-analysis data from peripheral blood samples, 5967 genes were tested using the SMR method. Another eQTL data was used to verify the results. Next we performed a fine-mapping association analysis to investigate the functional SNP in the identified loci. Weighted gene co-expression network analysis (WGCNA) was used to explore functional relationships for the identified novel genes with known putative osteoporosis genes. Further, we assessed functions of the identified genes through in vitro cellular study or previous functional studies. RESULTS: We identified two potentially causal genes (ASB16-AS1 and SYN2) associated with BMD. SYN2 was a novel osteoporosis candidate gene and ASB16-AS1 locus was known to be associated with BMD but was not the nearest gene to the top GWAS SNP. Fine-mapping association analysis showed that rs184478 and rs795000 was predicted to be possible causal SNPs in ASB16-AS1 and SYN2, respectively. ASB16-AS1 co-expressed with several known putative osteoporosis risk genes. In vitro cellular study showed that over-expressed ASB16-AS1 increased the expression of osteoblastogenesis related genes (BMP2 and ALPL), indicating its functional significance. CONCLUSION: Our findings support that ASB16-AS1 and SYN2 may represent two novel functional genes underlying BMD variation. The findings provide a basis for further functional mechanistic studies.