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.

Identification of Novel Pleiotropic SNPs Associated with Osteoporosis and Rheumatoid Arthritis.

Genome-wide association studies (GWASs) have identified hundreds of genetic loci for osteoporosis (OP) and rheumatoid arthritis (RA), individually, however, a large proportion of the total trait heritability remains unexplained. Previous studies demonstrated that these two diseases may share some common genetic determination and risk factors, but they were generally focused on individual trait and failed to identify the common variants that play key functional roles in the etiology of these two diseases. Here, we performed a conditional false discovery rate (cFDR) analysis to identify novel pleiotropic variants shared between them by integrating two independent GWASs with summary statistics for total body bone mineral density (TB-BMD, a major risk factor for osteoporosis) (n = 66,628) and RA (n = 58,284). A fine-mapping approach was also applied to identify the most probable causal variants with biological effects on both TB-BMD and RA. As a result, we found 47 independent pleiotropic SNPs shared between TB-BMD and RA, and 40 of them were validated in heel ultrasound estimated BMD (eBMD), femoral neck BMD (FN-BMD) or lumbar spine (LS-BMD). We detected one SNP (rs13299616) was novel and not identified by previous BMD or RA-related studies. Combined with fine-mapping and GWAS-eQTL colocalization analyses, our results suggested that locus 1p13.2 (including PTPN22, MAGI3, PHTF1, and RSBN1) was an important region to regulate TB-BMD and RA simultaneously. These findings provide new insights into the shared biological mechanisms and functional genetic determinants between OP and RA, and novel potential targets for treatment development.

Identifying Pleiotropic SNPs Associated With Femoral Neck and Heel Bone Mineral Density.

BACKGROUND: Genome-wide association studies (GWASs) routinely identify loci associated with risk factors for osteoporosis. However, GWASs with relatively small sample sizes still lack sufficient power to ascertain the majority of genetic variants with small to modest effect size, which may together truly influence the phenotype. The loci identified only account for a small percentage of the heritability of osteoporosis. This study aims to identify novel genetic loci associated with DXA-derived femoral neck (FNK) bone mineral density (BMD) and quantitative ultrasound of the heel calcaneus estimated BMD (eBMD), and to detect shared/causal variants for the two traits, to assess whether the SNPs or putative causal SNPs associated with eBMD were also associated with FNK-BMD. METHODS: Novel loci associated with eBMD and FNK-BMD were identified by the genetic pleiotropic conditional false discovery rate (cFDR) method. Shared putative causal variants between eBMD and FNK-BMD and putative causal SNPs for each trait were identified by the colocalization method. Mendelian randomization analysis addresses the causal relationship between eBMD/FNK-BMD and fracture. RESULTS: We identified 9,500 (cFDR < 9.8E-6), 137 (cFDR < 8.9E-4) and 124 SNPs associated with eBMD, FNK-BMD, and both eBMD and FNK-BMD, respectively, with 37 genomic regions where there was a SNP that influences both eBMD and FNK-BMD. Most genomic regions only contained putative causal SNPs associated with eBMD and 3 regions contained two distinct putative causal SNPs influenced both traits, respectively. We demonstrated a causal effect of FNK-BMD/eBMD on fracture. CONCLUSION: Most of SNPs or putative causal SNPs associated with FNK-BMD were also associated with eBMD. However, most of SNPs or putative causal SNPs associated with eBMD were not associated with FNK-BMD. The novel variants we identified may help to account for the additional proportion of variance of each trait and advance our understanding of the genetic mechanisms underlying osteoporotic fracture.

Identification of novel functional CpG-SNPs associated with type 2 diabetes and coronary artery disease.

Genome-wide association studies (GWASs) have identified hundreds of single nucleotide polymorphisms (SNPs) associated with type 2 diabetes (T2D) and coronary artery disease (CAD), respectively. Nevertheless, these studies were generally performed for single-trait/disease and failed to assess the pleiotropic role of the identified variants. To identify novel functional loci and the pleiotropic relationship between CAD and T2D, the targeted cFDR analysis on CpG-SNPs was performed by integrating two independent large and multi-centered GWASs with summary statistics of T2D (26,676 cases and 132,532 controls) and CAD (60,801 cases and 123,504 controls). Applying the cFDR significance threshold of 0.05, we observed a pleiotropic enrichment between T2D and CAD by incorporating pleiotropic effects into a conditional analysis framework. We identified 79 novel CpG-SNPs for T2D, 61 novel CpG-SNPs for CAD, and 18 novel pleiotropic loci for both traits. Among these novel CpG-SNPs, 33 of them were annotated as methylation quantitative trait locus (meQTL) in whole blood, and ten of them showed expression QTL (eQTL), meQTL, and metabolic QTL (metaQTL) effects simultaneously. To the best of our knowledge, we performed the first targeted cFDR analysis on CpG-SNPs, and our findings provided novel insights into the shared biological mechanisms and overlapped genetic heritability between T2D and CAD.

Additional common loci associated with stroke and obesity identified using pleiotropic analytical approach.

Stroke is a complex disease with multiple etiologies. Numerous studies suggest an established association between obesity and stroke, which may partly arise from the shared genetic components between the two phenotypes. Despite genome-wide association studies (GWASs) have identified some loci associated with stroke and obesity individually, the estimated genetic variability explained by these loci is limited (especially for stroke) and the pleiotropic loci between them are largely unknown. In this study, we jointly applied the pleiotropy-informed conditional false discovery rate (cFDR) method and the genetic analysis incorporating pleiotropy and annotation (GPA) method on summary statistics of two large GWASs to detect the genetic overlap between stroke (n = 446,696) and obesity (n = 681,275). Stratified Q-Q and fold-enrichment plots showed strong pleiotropic enrichment between the two phenotypes. With cFDR < 0.05 and fdr.GPA < 0.2, we identified 24 (16 novel) stroke-associated SNPs and 12 (10 novel) of them to be potentially pleiotropic SNPs for both phenotypes. The corresponding genes were enriched in trait-associated gene ontology (GO) terms "brain development" and "negative regulation of transport". In conclusion, our study demonstrated the feasibility and effectivity of the two pleiotropic methods which successfully improved the genetic discovery by incorporating related GWAS datasets and validated the genetic intercommunity between stroke and obesity. The identification of pleiotropic loci may provide us any new insights into potential genetic and etiology mechanism between them for the further studies.

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.

Identifying potentially common genes between dyslipidemia and osteoporosis using novel analytical approaches.

Dyslipidemia (DL) is closely related to osteoporosis (OP), while the exact common genetic mechanisms are still largely unknown. We proposed to use novel genetic analysis methods with pleiotropic information to identify potentially novel and/or common genes for the potential shared pathogenesis associated with OP and/or DL. We assessed the pleiotropy between plasma lipid (PL) and femoral neck bone mineral density (FNK BMD). We jointly applied the conditional false discovery rate (cFDR) method and the genetic analysis incorporating pleiotropy and annotation (GPA) method to the summary statistics provided by genome-wide association studies (GWASs) of FNK BMD (n = 49,988) and PL (n = 188,577) to identify potentially novel and/or common genes for BMD/PL. We found strong pleiotropic enrichment between PL and FNK BMD. Two hundred and forty-five PL SNPs were identified as potentially novel SNPs by cFDR and GPA. The corresponding genes were enriched in gene ontology (GO) terms "phospholipid homeostasis" and "chylomicron remnant clearance". Three SNPs (rs2178950, rs9939318, and rs9368716) might be the pleiotropic ones and the corresponding genes NLRC5 (rs2178950) and TRPS1 (rs9939318) were involved in NF-κB signaling pathway and Wnt signaling pathway as well as inflammation and innate immune processes. Our study validated the pleiotropy between PL and FNK BMD, and corroborated the reliability and high-efficiency of cFDR and GPA methods in further analyses of existing GWASs with summary statistics. We identified potentially common and/or novel genes for PL and/or FNK BMD, which may provide new insight and direction for further research.

Improved detection of genetic loci in estimated glomerular filtration rate and type 2 diabetes using a pleiotropic cFDR method.

Genome-wide association studies (GWAS) have been shown to have the potential of explaining more of the "missing heritability" of complex human phenotypes by improving statistical approaches. Here, we applied a genetic-pleiotropy-informed conditional false discovery rate (cFDR) to capture additional polygenic effects associated with estimated glomerular filtration rate (creatinine) (eGFRcrea) and type 2 diabetes (T2D). The cFDR analysis improves the identification of pleiotropic variants by incorporating potentially shared genetic mechanisms between two related traits. The Q-Q and fold-enrichment plots were used to assess the enrichment of SNPs associated with eGFRcrea or T2D, and Manhattan plots were used for showing chromosomal locations of the significant loci detected. By applying the cFDR method, we newly identified 74 loci for eGFRcrea and 3 loci for T2D with the cFDR criterion of 0.05 compared with previous related GWAS studies. Four shared SNPs were detected to be associated with both eGFRcrea and T2D at the significant conjunction cFDR level of 0.05, and these shared SNPs had not been reported in previous studies. In addition, we used DAVID analysis to perform functional analysis of the shared SNPs' annotated genes and found their potential hidden associations with eGFRcrea and T2D. In this study, the cFDR method shows the feasibility to detect more genetic variants underlying the heritability of eGFRcrea and T2D, and the overlapping SNPs identified could be regarded as candidate loci that provide a thread of genetic mechanisms between eGFRcrea and T2D in future research.

Linking Alzheimer's disease and type 2 diabetes: Novel shared susceptibility genes detected by cFDR approach.

BACKGROUND: Both type 2 diabetes (T2D) and Alzheimer's disease (AD) occur commonly in the aging populations and T2D has been considered as an important risk factor for AD. The heritability of both diseases is estimated to be over 50%. However, common pleiotropic single-nucleotide polymorphisms (SNPs)/loci have not been well-defined. The aim of this study is to analyze two large public accessible GWAS datasets to identify novel common genetic loci for T2D and/or AD. METHODS AND MATERIALS: The recently developed novel conditional false discovery rate (cFDR) approach was used to analyze the summary GWAS datasets from International Genomics of Alzheimer's Project (IGAP) and Diabetes Genetics Replication And Meta-analysis (DIAGRAM) to identify novel susceptibility genes for AD and T2D. RESULTS: We identified 78 SNPs (including 58 novel SNPs) that were associated with AD in Europeans conditional on T2D (cFDR<0.05). 66 T2D SNPs (including 40 novel SNPs) were identified by conditioning on SNPs association with AD (cFDR<0.05). A conjunction-cFDR (ccFDR) analysis detected 8 pleiotropic SNPs with a significance threshold of ccFDR<0.05 for both AD and T2D, of which 5 SNPs (rs6982393, rs4734295, rs7812465, rs10510109, rs2421016) were novel findings. Furthermore, among the 8 SNPs annotated at 6 different genes, 3 corresponding genes TP53INP1, TOMM40 and C8orf38 were related to mitochondrial dysfunction, critically involved in oxidative stress, which potentially contribute to the etiology of both AD and T2D. CONCLUSION: 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.