Our study provided evidence for shared genetic loci between T2D and AD in European subjects by using cFDR and ccFDR analyses. These results may provide novel insight into the etiology and potential therapeutic targets of T2D and/or AD.

Previous studies have demonstrated the genetic correlations between type 2 diabetes, obesity and dyslipidemia, and indicated that many genes have pleiotropic effects on them. However, these pleiotropic genes have not been well-defined. It is essential to identify pleiotropic genes using systematic approaches because systematically analyzing correlated traits is an effective way to enhance their statistical power. To identify potential pleiotropic genes for these three disorders, we performed a systematic analysis by incorporating GWAS (genome-wide associated study) datasets of six correlated traits related to type 2 diabetes, obesity and dyslipidemia using Meta-CCA (meta-analysis using canonical correlation analysis). Meta-CCA is an emerging method to systematically identify potential pleiotropic genes using GWAS summary statistics of multiple correlated traits. 2,720 genes were identified as significant genes after multiple testing (Bonferroni corrected p value < 0.05). Further, to refine the identified genes, we tested their relationship to the six correlated traits using VEGAS-2 (versatile gene-based association study-2). Only the genes significantly associated (Bonferroni corrected p value < 0.05) with more than one trait were kept. Finally, 25 genes (including two confirmed pleiotropic genes and eleven novel pleiotropic genes) were identified as potential pleiotropic genes. They were enriched in 5 pathways including the statin pathway and the PPAR (peroxisome proliferator-activated receptor) Alpha pathway. In summary, our study identified potential pleiotropic genes and pathways of type 2 diabetes, obesity and dyslipidemia, which may shed light on the common biological etiology and pathogenesis of these three disorders and provide promising insights for new therapies.

Context Although metabolic profiles appear to play an important role in menopausal bone loss, the functional mechanisms by which metabolites influence bone mineral density (BMD) during menopause are largely unknown. Objective We aimed to systematically identify metabolites associated with BMD variation and their potential functional mechanisms in peri-/post-menopausal women. Design and Methods We performed serum metabolomic profiling and whole-genome sequencing for 517 perimenopausal (16%) and early postmenopausal (84%) women aged 41 to 64 years in this cross-sectional study. Partial least squares (PLS) regression and general linear regression analysis were applied to identify BMD-associated metabolites, and weighted gene co-expression network analysis was performed to construct co-functional metabolite modules. Furthermore, we performed Mendelian randomization analysis to identify causal relationships between BMD-associated metabolites and BMD variation. Finally, we explored the effects of a novel prominent BMD-associated metabolite on bone metabolism through both in vivo/in vitro experiments. Results Twenty metabolites and a co-functional metabolite module (consisting of fatty acids) were significantly associated with BMD variation. We found dodecanoic acid (DA), within the identified module, causally decreased total hip BMD. Subsequently, the in vivo experiments might support that dietary supplementation with DA could promote bone loss, as well as increase the osteoblast and osteoclast numbers in normal/ovariectomized mice. DA treatment differentially promoted osteoblast and osteoclast differentiation, especially for osteoclast differentiation at higher concentrations in vitro (e.g.,10, 100μM). Conclusions This study sheds light on metabolomic profiles associated with postmenopausal osteoporosis risk, highlighting the potential importance of fatty acids, as exemplified by DA, in regulating BMD.

Interactions between multiple genes are involved in the development of complex diseases. However, there are few analyses of gene interactions associated with papillary thyroid cancer (PTC). Weighted gene co-expression network analysis (WGCNA) is a novel and powerful method that detects gene interactions according to their co-expression similarities. In the present study, WGCNA was performed in order to identify functional genes associated with PTC using R package. First, differential gene expression analysis was conducted in order to identify the differentially expressed genes (DEGs) between PTC and normal samples. Subsequently, co-expression networks of the DEGs were constructed for the two sample groups, respectively. The two networks were compared in order to identify a poorly preserved module. Concentrating on the significant module, validation analysis was performed to confirm the identified genes and combined functional enrichment analysis was conducted in order to identify more functional associations of these genes with PTC. As a result, 1062 DEGs were identified for network construction. A brown module containing 118 highly related genes was selected as it exhibited the lowest module preservation. After validation analysis, 61 genes in the module were confirmed to be associated with PTC. Following the enrichment analysis, two PTC-related pathways were identified: Wnt signal pathway and transcriptional misregulation in cancer. LRP4, KLK7, PRICKLE1, ETV4 and ETV5 were predicted to be candidate genes regulating the pathogenesis of PTC. These results provide novel insights into the etiology of PTC and the identification of potential functional genes.

Genome-wide association studies (GWASs) have discovered >50 risk loci for type 1 diabetes (T1D). However, those variations only have modest effects on the genetic risk of T1D. In recent years, accumulated studies have suggested that gene–gene interactions might explain part of the missing heritability. The purpose of our research was to identify potential and novel risk genes for T1D by systematically considering the gene–gene interactions through network analyses. We carried out a novel system network analysis of summary GWAS statistics jointly with transcriptomic gene expression data to identify some of the missing heritability for T1D using weighted gene co-expression network analysis (WGCNA). Using WGCNA, seven modules for 1852 nominally significant (P ≤ 0.05) GWAS genes were identified by analyzing microarray data for gene expression profile. One module (tagged as green module) showed significant association (P ≤ 0.05) between the module eigengenes and the trait. This module also displayed a high correlation (r = 0.45, P ≤ 0.05) between module membership (MM) and gene significant (GS), which indicated that the green module of co-expressed genes is of significant biological importance for T1D status. By further describing the module content and topology, the green module revealed a significant enrichment in the “regulation of immune response” (GO:0050776), which is a crucially important pathway in T1D development. Our findings demonstrated a module and several core genes that act as essential components in the etiology of T1D possibly via the regulation of immune response, which may enhance our fundamental knowledge of the underlying molecular mechanisms for T1D.