Oligogenic basis of premature ovarian insufficiency: an observational study.

BACKGROUND: The etiology of premature ovarian insufficiency, that is, the loss of ovarian activity before 40 years of age, is complex. Studies suggest that genetic factors are involved in 20-25% of cases. The aim of this study was to explore the oligogenic basis of premature ovarian insufficiency. RESULTS: Whole-exome sequencing of 93 patients with POI and whole-genome sequencing of 465 controls were performed. In the gene-burden analysis, multiple genetic variants, including those associated with DNA damage repair and meiosis, were more common in participants with premature ovarian insufficiency than in controls. The ORVAL-platform analysis confirmed the pathogenicity of the RAD52 and MSH6 combination. CONCLUSIONS: The results of this study indicate that oligogenic inheritance is an important cause of premature ovarian insufficiency and provide insights into the biological mechanisms underlying premature ovarian insufficiency.

Gut microbiota impacts bone via Bacteroides vulgatus-valeric acid-related pathways.

Although the gut microbiota has been reported to influence osteoporosis risk, the individual species involved, and underlying mechanisms, remain largely unknown. We performed integrative analyses in a Chinese cohort of peri-/post-menopausal women with metagenomics/targeted metabolomics/whole-genome sequencing to identify novel microbiome-related biomarkers for bone health. Bacteroides vulgatus was found to be negatively associated with bone mineral density (BMD), which was validated in US white people. Serum valeric acid (VA), a microbiota derived metabolite, was positively associated with BMD and causally downregulated by B. vulgatus. Ovariectomized mice fed B. vulgatus demonstrated increased bone resorption and poorer bone micro-structure, while those fed VA demonstrated reduced bone resorption and better bone micro-structure. VA suppressed RELA protein production (pro-inflammatory), and enhanced IL10 mRNA expression (anti-inflammatory), leading to suppressed maturation of osteoclast-like cells and enhanced maturation of osteoblasts in vitro. The findings suggest that B. vulgatus and VA may represent promising targets for osteoporosis prevention/treatment.

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.

Identification of potential genetic causal variants for obesity-related traits using statistical fine mapping.

BACKGROUND: Obesity is highly influenced by heritability and variant effects. While previous genome-wide association studies (GWASs) have successfully identified numerous genetic loci associated with obesity-related traits [body mass index (BMI) and waist-to-hip ratio (WHR)], most causal variants remain unidentified. The high degree of linkage disequilibrium (LD) throughout the genome makes it extremely difficult to distinguish the GWAS-associated SNPs that exert a true biological effect. OBJECTIVE: This study was to identify the potential causal variants having a biological effect on obesity-related traits. METHODS: We used Probabilistic Annotation INTegratOR, a Bayesian fine-mapping method, which incorporated genetic association data (GWAS summary statistics), LD structure, and functional annotations to calculate a posterior probability of causality for SNPs across all loci of interest. Moreover, we performed gene expression analysis using the available public transcriptomic data to validate the corresponding genes of the potential causal SNPs partially. RESULTS: We identified 96 SNPs for BMI and 43 SNPs for WHR with a high posterior probability of causality (> 99%), including 49 BMI SNPs and 24 WHR SNPs which did not reach genome-wide significance in the original GWAS. Finally, we partially validated some genes corresponding to the potential causal SNPs. CONCLUSION: Using a statistical fine-mapping approach, we identified a set of potential causal variants to be prioritized for future functional validation and also detected some novel trait-associated variants. These results provided novel insight into our understanding of the genetics of obesity and also demonstrated that fine mapping may improve upon the results identified by the original GWASs.

Systematic metabolomic studies identified adult adiposity biomarkers with acetylglycine associated with fat loss in vivo.

Obesity is associated with various adverse health outcomes. Body fat (BF) distribution is recognized as an important factor of negative health consequences of obesity. Although metabolomics studies, mainly focused on body mass index (BMI) and waist circumference, have explored the biological mechanisms involved in the development of obesity, these proxy composite measures are not accurate and cannot reflect BF distribution, and thus may hinder accurate assessment of metabolic alterations and differential risk of metabolic disorders among individuals presenting adiposity differently throughout the body. Thus, the exact relations between metabolites and BF remain to be elucidated. Here, we aim to examine the associations of metabolites and metabolic pathways with BF traits which reflect BF distribution. We performed systematic untargeted serum metabolite profiling and dual-energy X-ray absorptiometry (DXA) whole body fat scan for 517 Chinese women. We jointly analyzed DXA-derived four BF phenotypes to detect cross-phenotype metabolite associations and to prioritize important metabolomic factors. Topology-based pathway analysis was used to identify important BF-related biological processes. Finally, we explored the relationships of the identified BF-related candidate metabolites with BF traits in different sex and ethnicity through two independent cohorts. Acetylglycine, the top distinguished finding, was validated for its obesity resistance effect through in vivo studies of various diet-induced obese (DIO) mice. Eighteen metabolites and fourteen pathways were discovered to be associated with BF phenotypes. Six of the metabolites were validated in varying sex and ethnicity. The obesity-resistant effects of acetylglycine were observed to be highly robust and generalizable in both human and DIO mice. These findings demonstrate the importance of metabolites associated with BF distribution patterns and several biological pathways that may contribute to obesity and obesity-related disease etiology, prevention, and intervention. Acetylglycine is highlighted as a potential therapeutic candidate for preventing excessive adiposity in future studies.

Single-cell RNA sequencing reveals in vivo osteoimmunology interactions between the immune and skeletal systems.

Background: While osteoimmunology interactions between the immune and skeletal systems are known to play an important role in osteoblast development, differentiation and bone metabolism related disease like osteoporosis, such interactions in either bone microenvironment or peripheral circulation in vivo at the single-cell resolution have not yet been characterized. Methods: We explored the osteoimmunology communications between immune cells and osteoblastic lineage cells (OBCs) by performing CellphoneDB and CellChat analyses with single-cell RNA sequencing (scRNA-seq) data from human femoral head. We also explored the osteoimmunology effects of immune cells in peripheral circulation on skeletal phenotypes. We used a scRNA-seq dataset of peripheral blood monocytes (PBMs) to perform deconvolution analysis. Then weighted gene co-expression network analysis (WGCNA) was used to identify monocyte subtype-specific subnetworks. We next used cell-specific network (CSN) and the least absolute shrinkage and selection operator (LASSO) to analyze the correlation of a gene subnetwork identified by WGCNA with bone mineral density (BMD). Results: We constructed immune cell and OBC communication networks and further identified L-R genes, such as JAG1 and NOTCH1/2, with ossification related functions. We also found a Mono4 related subnetwork that may relate to BMD variation in both older males and postmenopausal female subjects. Conclusions: This is the first study to identify numerous ligand-receptor pairs that likely mediate signals between immune cells and osteoblastic lineage cells. This establishes a foundation to reveal advanced and in-depth osteoimmunology interactions to better understand the relationship between local bone microenvironment and immune cells in peripheral blood and the impact on bone phenotypes.

Repurposing Approved Drugs for Sarcopenia Based on Transcriptomics Data in Humans.

Sarcopenia, characterized by age-related loss of muscle mass, strength, and decreased physical performance, is a growing public health challenge amid the rapidly ageing population. As there are no approved drugs that target sarcopenia, it has become increasingly urgent to identify promising pharmacological interventions. In this study, we conducted an integrative drug repurposing analysis utilizing three distinct approaches. Firstly, we analyzed skeletal muscle transcriptomic sequencing data in humans and mice using gene differential expression analysis, weighted gene co-expression analysis, and gene set enrichment analysis. Subsequently, we employed gene expression profile similarity assessment, hub gene expression reversal, and disease-related pathway enrichment to identify and repurpose candidate drugs, followed by the integration of findings with rank aggregation algorithms. Vorinostat, the top-ranking drug, was also validated in an in vitro study, which demonstrated its efficacy in promoting muscle fiber formation. Although still requiring further validation in animal models and human clinical trials, these results suggest a promising drug repurposing prospect in the treatment and prevention of sarcopenia.

Regulon active landscape reveals cell development and functional state changes of human primary osteoblasts in vivo.

BACKGROUND: While transcription factor (TF) regulation is known to play an important role in osteoblast development, differentiation, and bone metabolism, the molecular features of TFs in human osteoblasts at the single-cell resolution level have not yet been characterized. Here, we identified modules (regulons) of co-regulated genes by applying single-cell regulatory network inference and clustering to the single-cell RNA sequencing profiles of human osteoblasts. We also performed cell-specific network (CSN) analysis, reconstructed regulon activity-based osteoblast development trajectories, and validated the functions of important regulons both in vivo and in vitro. RESULTS: We identified four cell clusters: preosteoblast-S1, preosteoblast-S2, intermediate osteoblasts, and mature osteoblasts. CSN analysis results and regulon activity-based osteoblast development trajectories revealed cell development and functional state changes of osteoblasts. CREM and FOSL2 regulons were mainly active in preosteoblast-S1, FOXC2 regulons were mainly active in intermediate osteoblast, and RUNX2 and CREB3L1 regulons were most active in mature osteoblasts. CONCLUSIONS: This is the first study to describe the unique features of human osteoblasts in vivo based on cellular regulon active landscapes. Functional state changes of CREM, FOSL2, FOXC2, RUNX2, and CREB3L1 regulons regarding immunity, cell proliferation, and differentiation identified the important cell stages or subtypes that may be predominantly affected by bone metabolism disorders. These findings may lead to a deeper understanding of the mechanisms underlying bone metabolism and associated diseases.

An autoencoder-based deep learning method for genotype imputation.

Genotype imputation has a wide range of applications in genome-wide association study (GWAS), including increasing the statistical power of association tests, discovering trait-associated loci in meta-analyses, and prioritizing causal variants with fine-mapping. In recent years, deep learning (DL) based methods, such as sparse convolutional denoising autoencoder (SCDA), have been developed for genotype imputation. However, it remains a challenging task to optimize the learning process in DL-based methods to achieve high imputation accuracy. To address this challenge, we have developed a convolutional autoencoder (AE) model for genotype imputation and implemented a customized training loop by modifying the training process with a single batch loss rather than the average loss over batches. This modified AE imputation model was evaluated using a yeast dataset, the human leukocyte antigen (HLA) data from the 1,000 Genomes Project (1KGP), and our in-house genotype data from the Louisiana Osteoporosis Study (LOS). Our modified AE imputation model has achieved comparable or better performance than the existing SCDA model in terms of evaluation metrics such as the concordance rate (CR), the Hellinger score, the scaled Euclidean norm (SEN) score, and the imputation quality score (IQS) in all three datasets. Taking the imputation results from the HLA data as an example, the AE model achieved an average CR of 0.9468 and 0.9459, Hellinger score of 0.9765 and 0.9518, SEN score of 0.9977 and 0.9953, and IQS of 0.9515 and 0.9044 at missing ratios of 10% and 20%, respectively. As for the results of LOS data, it achieved an average CR of 0.9005, Hellinger score of 0.9384, SEN score of 0.9940, and IQS of 0.8681 at the missing ratio of 20%. In summary, our proposed method for genotype imputation has a great potential to increase the statistical power of GWAS and improve downstream post-GWAS analyses.

Identification of genetic loci shared between Alzheimer's disease and hypertension.

Alzheimer's disease (AD) and high blood pressure (BP) are prevalent age-related diseases with significant unexplained heritability. A thorough analysis of genetic pleiotropy between AD and BP will lay a foundation for the study of the associated molecular mechanisms, leading to a better understanding of the development of each phenotype. We used the conditional false discovery rate (cFDR) method to identify novel genetic loci associated with both AD and BP. The cFDR approach improves the effective sample size for association testing by combining GWAS summary statistics for correlated phenotypes. We identified 50 pleiotropic SNPs for AD and BP, 7 of which are novel and have not previously been reported to be associated with either AD or BP. The novel SNPs located at STK3 are particularly noteworthy, as this gene may influence AD risk via the Hippo signaling network, which regulates cell death. Bayesian colocalization analysis demonstrated that although AD and BP are associated, they do not appear to share the same causal variants. We further performed two sample Mendelian randomization analysis, but could not detect a causal effect of BP on AD. Despite the inability to establish a causal link between AD and BP, our findings report some potential novel pleiotropic loci that may influence disease susceptibility. In summary, we identified 7 SNPs that annotate to 4 novel genes which have not previously been reported to be associated with AD nor with BP and discuss the possible role of one of these genes, STK3 in the Hippo signaling network.

Gut microbiota accelerates obesity in peri-/post-menopausal women via Bacteroides fragilis and acetic acid.

OBJECTIVE: Many animal experiments and epidemiological studies have shown that the gut microbiota (GM) plays an important role in the development of obesity, but the specific biological mechanism involved in the pathogenesis of disease remain unknown. We aimed to examine the relationships and functional mechanisms of GM on obesity in peri- and post-menopausal women. METHODS: We recruited 499 Chinese peri- and post-menopausal women and performed comprehensive analyses of the gut microbiome, targeted metabolomics for short-chain fatty acids in serum, and host whole-genome sequencing by various association analysis methods. RESULTS: Through constrained linear regression analysis, we found that an elevated abundance of Bacteroides fragilis (B. fragilis) was associated with obesity. We also found that serum levels of acetic acid were negatively associated with obesity, and that B. fragilis was negatively associated with serum acetic acid levels by partial Spearman correlation analysis. Mendelian randomization analysis indicated that B. fragilis increases the risk of obesity and may causally down-regulate acetic acid levels. CONCLUSIONS: We found the gut with B. fragilis may accelerate obesity, in part, by suppressing acetic acid levels. Therefore, B. fragilis and acetic acid may represent important therapeutic targets for obesity intervention in peri- and post-menopausal women.

Integration of the Human Gut Microbiome and Serum Metabolome Reveals Novel Biological Factors Involved in the Regulation of Bone Mineral Density.

While the gut microbiome has been reported to play a role in bone metabolism, the individual species and underlying functional mechanisms have not yet been characterized. We conducted a systematic multi-omics analysis using paired metagenomic and untargeted serum metabolomic profiles from a large sample of 499 peri- and early post-menopausal women to identify the potential crosstalk between these biological factors which may be involved in the regulation of bone mineral density (BMD). Single omics association analyses identified 22 bacteria species and 17 serum metabolites for putative association with BMD. Among the identified bacteria, Bacteroidetes and Fusobacteria were negatively associated, while Firmicutes were positively associated. Several of the identified serum metabolites including 3-phenylpropanoic acid, mainly derived from dietary polyphenols, and glycolithocholic acid, a secondary bile acid, are metabolic byproducts of the microbiota. We further conducted a supervised integrative feature selection with respect to BMD and constructed the inter-omics partial correlation network. Although still requiring replication and validation in future studies, the findings from this exploratory analysis provide novel insights into the interrelationships between the gut microbiome and serum metabolome that may potentially play a role in skeletal remodeling processes.

Identification of PDXDC1 as a novel pleiotropic susceptibility locus shared between lumbar spine bone mineral density and birth weight.

An increasing number of epidemiological studies have suggested that birth weight (BW) may be a determinant of bone health later in life, although the underlying genetic mechanism remains unclear. Here, we applied a pleiotropic conditional false discovery rate (cFDR) approach to the genome-wide association study (GWAS) summary statistics for lumbar spine bone mineral density (LS BMD) and BW, aiming to identify novel susceptibility variants shared between these two traits. We detected 5 novel potential pleiotropic loci which are located at or near 7 different genes (NTAN1, PDXDC1, CACNA1G, JAG1, FAT1P1, CCDC170, ESR1), among which PDXDC1 and FAT1P1 have not previously been linked to these phenotypes. To partially validate the findings, we demonstrated that the expression of PDXDC1 was dramatically reduced in ovariectomized (OVX) mice in comparison with sham-operated (SHAM) mice in both the growth plate and trabecula bone. Furthermore, immunohistochemistry assay with serial sections showed that both osteoclasts and osteoblasts express PDXDC1, supporting its potential role in bone metabolism. In conclusion, our study provides insights into some shared genetic mechanisms for BMD and BW as well as a novel potential therapeutic target for the prevention of OP in the early stages of the disease development. KEY MESSAGES : We investigated pleiotropy-informed enrichment between LS BMD and BW. We identified genetic variants related to both LS BMD and BW by utilizing a cFDR approach. PDXDC1 is a novel pleiotropic gene which may be related to both LS BMD and BW. Elevated expression of PDXDC1 is related to higher BMD and lower ratio n-6/n-3 PUFA indicating a bone protective effect of PDXDC1.

A multiethnic whole genome sequencing study to identify novel loci for bone mineral density.

At present, there have only been a few DNA sequencing-based studies to explore the genetic determinants of bone mineral density (BMD). We carried out the largest whole genome sequencing analysis to date for femoral neck and spine BMD (n = 4981), with one of the highest average sequencing depths implemented thus far at 22×, in a multiethnic sample (58% Caucasian and 42% African American) from the Louisiana Osteoporosis Study (LOS). The LOS samples were combined with summary statistics from the GEFOS consortium and several independent samples of various ethnicities to perform GWAS meta-analysis (n = 44 506). We identified 31 and 30 genomic risk loci for femoral neck and spine BMD, respectively. The findings substantiate many previously reported susceptibility loci (e.g. WNT16 and ESR1) and reveal several others that are either novel or have not been widely replicated in GWAS for BMD, including two for femoral neck (IGF2 and ZNF423) and one for spine (SIPA1). Although we were not able to uncover ethnicity specific differences in the genetic determinants of BMD, we did identify several loci which demonstrated sex-specific associations, including two for women (PDE4D and PIGN) and three for men (TRAF3IP2, NFIB and LYSMD4). Gene-based rare variant association testing detected MAML2, a regulator of the Notch signaling pathway, which has not previously been suggested, for association with spine BMD. The findings provide novel insights into the pathophysiological mechanisms of osteoporosis.

Single-cell RNA sequencing deconvolutes the in vivo heterogeneity of human bone marrow-derived mesenchymal stem cells.

Bone marrow-derived mesenchymal stem cells (BM-MSCs) are multipotent stromal cells that have a critical role in the maintenance of skeletal tissues such as bone, cartilage, and the fat in bone marrow. In addition to providing microenvironmental support for hematopoietic processes, BM-MSCs can differentiate into various mesodermal lineages including osteoblast/osteocyte, chondrocyte, and adipocyte that are crucial for bone metabolism. While BM-MSCs have high cell-to-cell heterogeneity in gene expression, the cell subtypes that contribute to this heterogeneity in vivo in humans have not been characterized. To investigate the transcriptional diversity of BM-MSCs, we applied single-cell RNA sequencing (scRNA-seq) on freshly isolated CD271+ BM-derived mononuclear cells (BM-MNCs) from two human subjects. We successfully identified LEPRhiCD45low BM-MSCs within the CD271+ BM-MNC population, and further codified the BM-MSCs into distinct subpopulations corresponding to the osteogenic, chondrogenic, and adipogenic differentiation trajectories, as well as terminal-stage quiescent cells. Biological functional annotations of the transcriptomes suggest that osteoblast precursors induce angiogenesis coupled with osteogenesis, and chondrocyte precursors have the potential to differentiate into myocytes. We also discovered transcripts for several clusters of differentiation (CD) markers that were either highly expressed (e.g., CD167b, CD91, CD130 and CD118) or absent (e.g., CD74, CD217, CD148 and CD68) in BM-MSCs, representing potential novel markers for human BM-MSC purification. This study is the first systematic in vivo dissection of human BM-MSCs cell subtypes at the single-cell resolution, revealing an insight into the extent of their cellular heterogeneity and roles in maintaining bone homeostasis.

A systematic dissection of human primary osteoblasts in vivo at single-cell resolution.

Human osteoblasts are multifunctional bone cells, which play essential roles in bone formation, angiogenesis regulation, as well as maintenance of hematopoiesis. However, the categorization of primary osteoblast subtypes in vivo in humans has not yet been achieved. Here, we used single-cell RNA sequencing (scRNA-seq) to perform a systematic cellular taxonomy dissection of freshly isolated human osteoblasts from one 31-year-old male with osteoarthritis and osteopenia after hip replacement. Based on the gene expression patterns and cell lineage reconstruction, we identified three distinct cell clusters including preosteoblasts, mature osteoblasts, and an undetermined rare osteoblast subpopulation. This novel subtype was found to be the major source of the nuclear receptor subfamily 4 group A member 1 and 2 (NR4A1 and NR4A2) in primary osteoblasts, and the expression of NR4A1 was confirmed by immunofluorescence staining on mouse osteoblasts in vivo. Trajectory inference analysis suggested that the undetermined cluster, together with the preosteoblasts, are involved in the regulation of osteoblastogenesis and also give rise to mature osteoblasts. Investigation of the biological processes and signaling pathways enriched in each subpopulation revealed that in addition to bone formation, preosteoblasts and undetermined osteoblasts may also regulate both angiogenesis and hemopoiesis. Finally, we demonstrated that there are systematic differences between the transcriptional profiles of human and mouse osteoblasts, highlighting the necessity for studying bone physiological processes in humans rather than solely relying on mouse models. Our findings provide novel insights into the cellular heterogeneity and potential biological functions of human primary osteoblasts at the single-cell level.

A transcriptome-wide association study to detect novel genes for volumetric bone mineral density.

Transcriptome-wide association studies (TWAS) systematically investigate the association of genetically predicted gene expression with disease risk, providing an effective approach to identify novel susceptibility genes. Osteoporosis is the most common metabolic bone disease, associated with reduced bone mineral density (BMD) and increased risk of osteoporotic fractures, whereas genetic factors explain approximately 70% of the variance in phenotypes associated with bone. BMD is commonly assessed using dual-energy X-ray absorptiometry (DXA) to obtain measurements (g/cm2) of areal BMD. However, quantitative computed tomography (QCT) measured 3D volumetric BMD (vBMD) (g/cm3) has important advantages compared with DXA since it can evaluate cortical and trabecular microstructural features of bone quality, which can be used to directly predict fracture risk. Here, we performed the first TWAS for volumetric BMD (vBMD) by integrating genome-wide association studies (GWAS) data from two independent cohorts, namely the Framingham Heart Study (FHS, n = 3298) and the Osteoporotic Fractures in Men (MrOS, n = 4641), with tissue-specific gene expression data from the Genotype-Tissue Expression (GTEx) project. We first used stratified linkage disequilibrium (LD) score regression approach to identify 12 vBMD-relevant tissues, for which vBMD heritability is enriched in tissue-specific genes of the given tissue. Focusing on these tissues, we subsequently leveraged GTEx expression reference panels to predict tissue-specific gene expression levels based on the genotype data from FHS and MrOS. The associations between predicted gene expression levels and vBMD variation were then tested by MultiXcan, an innovative TWAS method that integrates information available across multiple tissues. We identified 70 significant genes associated with vBMD, including some previously identified osteoporosis-related genes such as LYRM2 and NME8, as well as some novel loci such as DNAAF2 and SPAG16. Our findings provide novel insights into the pathophysiological mechanisms of osteoporosis and highlight several novel vBMD-associated genes that warrant further investigation.

Drug Repurposing for COVID-19 Treatment by Integrating Network Pharmacology and Transcriptomics.

Since coronavirus disease 2019 (COVID-19) is a serious new worldwide public health crisis with significant morbidity and mortality, effective therapeutic treatments are urgently needed. Drug repurposing is an efficient and cost-effective strategy with minimum risk for identifying novel potential treatment options by repositioning therapies that were previously approved for other clinical outcomes. Here, we used an integrated network-based pharmacologic and transcriptomic approach to screen drug candidates novel for COVID-19 treatment. Network-based proximity scores were calculated to identify the drug-disease pharmacological effect between drug-target relationship modules and COVID-19 related genes. Gene set enrichment analysis (GSEA) was then performed to determine whether drug candidates influence the expression of COVID-19 related genes and examine the sensitivity of the repurposing drug treatment to peripheral immune cell types. Moreover, we used the complementary exposure model to recommend potential synergistic drug combinations. We identified 18 individual drug candidates including nicardipine, orantinib, tipifarnib and promethazine which have not previously been proposed as possible treatments for COVID-19. Additionally, 30 synergistic drug pairs were ultimately recommended including fostamatinib plus tretinoin and orantinib plus valproic acid. Differential expression genes of most repurposing drugs were enriched significantly in B cells. The findings may potentially accelerate the discovery and establishment of an effective therapeutic treatment plan for COVID-19 patients.

Identification of novel functional CpG-SNPs associated with Type 2 diabetes and birth weight.

Genome-wide association studies (GWASs) have identified hundreds of genetic loci for type 2 diabetes (T2D) and birth weight (BW); however, a large proportion of the total trait heritability remains unexplained. The previous studies were generally focused on individual traits and largely failed to identify the majority of the variants that play key functional roles in the etiology of the disease. Here, we aim to identify novel functional loci for T2D, BW and the pleiotropic variants shared between them by performing a targeted conditional false discovery rate (cFDR) analysis that integrates two independent GWASs with summary statistics for T2D (n = 26,676 cases and 132,532 controls) and BW (n = 153,781) which entails greater statistical power than individual trait analyses. In this analysis, we considered CpG-SNPs, which are SNPs that may influence DNA methylation status, and are therefore considered to be functionally important. We identified 103 novel CpG-SNPs for T2D, 182 novel CpG-SNPs for BW (cFDR < 0.05), and 52 novel pleiotropic loci for both (conjunction cFDR [ccFDR] < 0.05). Among the identified novel CpG-SNPs, 33 were annotated as methylation quantitative trait loci (meQTLs) in whole blood, and 145 displayed at least some effects on meQTL, metabolic QTL (metaQTL), and/or expression QTL (eQTL). These findings may provide further insights into the shared biological mechanisms and functional genetic determinants that overlap between T2D and BW, thereby providing novel potential targets for treatment/intervention development.

Cell-specific network analysis of human folliculogenesis reveals network rewiring in antral stage oocytes.

Although previous studies have explored the gene expression profiles of human oocytes and granulosa cells by single-cell RNA sequencing (scRNA-seq), the dynamic regulatory network at a single-cell resolution during folliculogenesis remains largely unknown. We identified 10 functional modules by WGCNA, four of which were significantly correlated with primary/antral oocyte and antral/pre-ovulatory granulosa cells. Functional enrichment analysis showed that the brown module, which was correlated with antral oocyte, was enriched in oocyte differentiation, and two core subnetworks identified by MCODE were involved in cell cycle (blue subnetwork) and oogenesis (red subnetwork). The cell-specific network (CSN) analysis demonstrated a distinct gene network structure associated with the antral follicular stage, which was notably different from other developmental stages. To our knowledge, this is the first study to explore gene functions during folliculogenesis at single-cell network level. We uncovered two potential gene subnetworks, which may play an important role in oocyte function beginning at the antral stage, and further established their rewiring process at intra-network/whole transcriptome level. The findings provide crucial insights from a novel network perspective to be further explored in functional mechanistic studies.