Meta-analysis of gene-based genome-wide association studies of bone mineral density in Chinese and European subjects

Decoding the genetic relationship between Alzheimer's disease and type 2 diabetes: potential risk variants and future direction for North Africa

Introduction

Alzheimer's disease (AD) and Type 2 diabetes (T2D) are both age-associated diseases. Identification of shared genes could help develop early diagnosis and preventive strategies. Although genetic background plays a crucial role in these diseases, we noticed an underrepresentation tendency of North African populations in omics studies.

Materials and methods

First, we conducted a comprehensive review of genes and pathways shared between T2D and AD through PubMed. Then, the function of the identified genes and variants was investigated using annotation tools including PolyPhen2, RegulomeDB, and miRdSNP. Pathways enrichment analyses were performed with g:Profiler and EnrichmentMap. Next, we analyzed variant distributions in 16 worldwide populations using PLINK2, R, and STRUCTURE software. Finally, we performed an inter-ethnic comparison based on the minor allele frequency of T2D-AD common variants.

Results

A total of 59 eligible papers were included in our study. We found 231 variants and 363 genes shared between T2D and AD. Variant annotation revealed six single nucleotide polymorphisms (SNP) with a high pathogenic score, three SNPs with regulatory effects on the brain, and six SNPs with potential effects on miRNA-binding sites. The miRNAs affected were implicated in T2D, insulin signaling pathways, and AD. Moreover, replicated genes were significantly enriched in pathways related to plasma protein binding, positive regulation of amyloid fibril deposition, microglia activation, and cholesterol metabolism. Multidimensional screening performed based on the 363 shared genes showed that main North African populations are clustered together and are divergent from other worldwide populations. Interestingly, our results showed that 49 SNP associated with T2D and AD were present in North African populations. Among them, 11 variants located in DNM3, CFH, PPARG, ROHA, AGER, CLU, BDNF1, CST9, and PLCG1 genes display significant differences in risk allele frequencies between North African and other populations.

Conclusion

Our study highlighted the complexity and the unique molecular architecture of North African populations regarding T2D-AD shared genes. In conclusion, we emphasize the importance of T2D-AD shared genes and ethnicity-specific investigation studies for a better understanding of the link behind these diseases and to develop accurate diagnoses using personalized genetic biomarkers.

Gene-Set Enrichment Analysis for Identifying Genes and Biological Activities Associated with Growth Traits in Dromedaries

Growth is an important heritable economic trait for dromedaries and necessary for planning a successful breeding program. Until now, genome-wide association studies (GWAS) and QTL-mapping have identified significant single nucleotide polymorphisms (SNPs) associated with growth in domestic animals, but in dromedaries, the number of studies is very low. This project aimed to find biological themes affecting growth in dromedaries. In the first step, 99 candidate SNPs were chosen from a previously established set of SNPs associated with body weight, gain, and birth weight in Iranian dromedaries. Next, 0.5 kb upstream and downstream of each candidate SNP were selected from NCBI (assembly accession: GCA_000803125.3). The annotation of fragments with candidate SNPs regarding the reference genome was retrieved using the Blast2GO tool. Candidate SNPs associated with growth were mapped to 22 genes, and 25 significant biological themes were identified to be related to growth in dromedaries. The main biological functions included calcium ion binding, protein binding, DNA-binding transcription factor activity, protein kinase activity, tropomyosin binding, myosin complex, actin-binding, ATP binding, receptor signaling pathway via JAK-STAT, and cytokine activity. EFCAB5, MTIF2, MYO3A, TBX15, IFNL3, PREX1, and TMOD3 genes are candidates for improving growth in camel breeding programs.

Three functional polymorphisms in CCDC170 were associated with osteoporosis phenotype

MicroRNAs (miRNAs) play essential roles in regulating bone formation and homeostasis. Genomic variations within miRNA target sites may therefore be important sources of genetic differences in osteoporosis risk. The function of CCDC170 in bone biology is still unclear. To verify the function of CCDC170, we knocked down CCDC170 in cells and mice and searched for miRNA recognition sites within CCDC170 using the TargetScan, miRNASNP, and miRBase databases. In this study, our results demonstrated that CCDC170 plays an important role in the positive regulation of bone formation. MiR-153-3p, miR-374b-3p, miR-4274, miR-572 and miR-2964a-5p inhibited CCDC170 expression in an allele-specific manner by binding GWAS lead SNPs rs6932603, rs3757322 and rs3734806. These findings may improve our understanding of the association between CCDC170, miRNAs, GWAS lead SNPs, and osteoporosis pathogenesis and may provide a potential therapeutic target for osteoporosis therapy.

Summary: CCDC170 may provide a potential therapeutic target for osteoporosis therapy.

Genome-wide association study using family-based cohorts identifies the WLS and CCDC170/ESR1 loci as associated with bone mineral density

Background

Osteoporosis is a common and debilitating bone disease that is characterised by a low bone mineral density (BMD), a highly heritable trait. Genome-wide association studies (GWAS) have proven to be very successful in identifying common genetic variants associated with BMD adjusted for age, gender and weight, however a large portion of the genetic variance for this trait remains unexplained. There is evidence to suggest significant genetic correlation between body size traits and BMD. It has also recently been suggested that unintended bias can be introduced as a result of adjusting a phenotype for a correlated trait. We performed a GWAS meta-analysis in two populations (total n = 6,696) using BMD data adjusted for only age and gender, in an attempt to identify genetic variants associated with BMD including those that may have potential pleiotropic effects on BMD and body size traits.

Results

We observed a single variant, rs2566752, associated with spine BMD at the genome-wide significance level in the meta-analysis (P = 3.36 × 10⁻⁰⁹). Logistic regression analysis also revealed an association between rs2566752 and fracture rate in one of our study cohorts (P = 0.017, n = 5,654). This is an intronic variant located in the wntless Wnt ligand secretion mediator (WLS) gene (1p31.3), a known BMD locus which encodes an integral component of the Wnt ligand secretion pathway. Bioinformatics analyses of variants in moderate LD with rs2566752 produced strong evidence for a regulatory role for the variants rs72670452, rs17130567 and rs1430738. Expression quantitative trait locus (eQTL) analysis suggested that the variants rs12568456 and rs17130567 are associated with expression of the WLS gene in whole blood, cerebellum and temporal cortex brain tissue (P = 0.034–1.19 × 10⁻²³). Gene-wide association testing using the VErsatile Gene-based Association Study 2 (VEGAS2) software revealed associations between the coiled-coil domain containing 170 (CCDC170) gene, located adjacent to the oestrogen receptor 1 (ESR1) gene, and BMD at the spine, femoral neck and total hip sites (P = 1.0 × 10⁻⁰⁶, 2.0 × 10⁻⁰⁶ and 2.0 × 10⁻⁰⁶ respectively).

Conclusions

Genetic variation at the WLS and CCDC170/ESR1 loci were found to be significantly associated with BMD adjusted for only age and gender at the genome-wide level in this meta-analysis.

Gene-based association analysis identified novel genes associated with bone mineral density

Genetic factors contribute to the variation of bone mineral density (BMD), which is a major risk factor of osteoporosis. The aim of this study was to identify more “novel” genes for BMD. Based on the publicly available SNP-based P values, we performed an initial gene-based analysis in a total of 32,961 individuals. Furthermore, we performed differential expression, pathway and protein-protein interaction analyses to find supplementary evidence to support the significance of the identified genes. About 21,695 genes for femoral neck (FN)-BMD and 21,683 genes for lumbar spine (LS)-BMD were analyzed using gene-based association analysis. A total of 35 FN-BMD associated genes and 53 LS-BMD associated genes were identified (P < 2.3×10^-6) after Bonferroni correction. Among them, 64 genes have not been reported in previous SNP-based genome-wide association studies. Differential expression analysis further supported the significant associations of 14 genes with FN-BMD and 19 genes with LS-BMD. Especially, WNT3 and WNT9B in the Wnt signaling pathway for FN-BMD were further supported by pathway analysis and protein-protein interaction analysis. The present study took the advantage of gene-based association method to perform a supplementary analysis of the GWAS dataset and found some BMD-associated genes. The evidence taken together supported the importance of Wnt signaling pathway genes in determining osteoporosis. Our findings provided more insights into the genetic basis of osteoporosis.

Genome-wide and gene-based association studies of anxiety disorders in European and African American samples

Anxiety disorders (ADs) are common mental disorders caused by a combination of genetic and environmental factors. Since ADs are highly comorbid with each other, partially due to shared genetic basis, studying AD phenotypes in a coordinated manner may be a powerful strategy for identifying potential genetic loci for ADs. To detect these loci, we performed genome-wide association studies (GWAS) of ADs. In addition, as a complementary approach to single-locus analysis, we also conducted gene- and pathway-based analyses. GWAS data were derived from the control sample of the Molecular Genetics of Schizophrenia (MGS) project (2,540 European American and 849 African American subjects) genotyped on the Affymetrix GeneChip 6.0 array. We applied two phenotypic approaches: (1) categorical case-control comparisons (CC) based upon psychiatric diagnoses, and (2) quantitative phenotypic factor scores (FS) derived from a multivariate analysis combining information across the clinical phenotypes. Linear and logistic models were used to analyse the association with ADs using FS and CC traits, respectively. At the single locus level, no genome-wide significant association was found. A trans-population gene-based meta-analysis across both ethnic subsamples using FS identified three genes (MFAP3L on 4q32.3, NDUFAB1 and PALB2 on 16p12) with genome-wide significance (false discovery rate (FDR] <5%). At the pathway level, several terms such as transcription regulation, cytokine binding, and developmental process were significantly enriched in ADs (FDR <5%). Our approaches studying ADs as quantitative traits and utilizing the full GWAS data may be useful in identifying susceptibility genes and pathways for ADs.

Integrative analysis of GWASs, human protein interaction, and gene expression identified gene modules associated with BMDs

CONTEXT

To date, few systems genetics studies in the bone field have been performed. We designed our study from a systems-level perspective by integrating genome-wide association studies (GWASs), human protein-protein interaction (PPI) network, and gene expression to identify gene modules contributing to osteoporosis risk.

METHODS

First we searched for modules significantly enriched with bone mineral density (BMD)-associated genes in human PPI network by using 2 large meta-analysis GWAS datasets through a dense module search algorithm. One included 7 individual GWAS samples (Meta7). The other was from the Genetic Factors for Osteoporosis Consortium (GEFOS2). One was assigned as a discovery dataset and the other as an evaluation dataset, and vice versa.

RESULTS

In total, 42 modules and 129 modules were identified significantly in both Meta7 and GEFOS2 datasets for femoral neck and spine BMD, respectively. There were 3340 modules identified for hip BMD only in Meta7. As candidate modules, they were assessed for the biological relevance to BMD by gene set enrichment analysis in 2 expression profiles generated from circulating monocytes in subjects with low versus high BMD values. Interestingly, there were 2 modules significantly enriched in monocytes from the low BMD group in both gene expression datasets (nominal P value <.05). Two modules had 16 nonredundant genes. Functional enrichment analysis revealed that both modules were enriched for genes involved in Wnt receptor signaling and osteoblast differentiation.

CONCLUSION

We highlighted 2 modules and novel genes playing important roles in the regulation of bone mass, providing important clues for therapeutic approaches for osteoporosis.

Genome-wide approaches for identifying genetic risk factors for osteoporosis

Osteoporosis, the most common type of bone disease worldwide, is clinically characterized by low bone mineral density (BMD) and increased susceptibility to fracture. Multiple genetic and environmental factors and gene-environment interactions have been implicated in its pathogenesis. Osteoporosis has strong genetic determination, with the heritability of BMD estimated to be as high as 60%. More than 80 genes or genetic variants have been implicated in risk of osteoporosis by hypothesis-free genome-wide studies. However, these genes or genetic variants can only explain a small portion of BMD variation, suggesting that many other genes or genetic variants underlying osteoporosis risk await discovery. Here, we review recent progress in genome-wide studies of osteoporosis and discuss their implications for medicine and the major challenges in the field.

The transcriptional profile of mesenchymal stem cell populations in primary osteoporosis is distinct and shows overexpression of osteogenic inhibitors

Primary osteoporosis is an age-related disease characterized by an imbalance in bone homeostasis. While the resorptive aspect of the disease has been studied intensely, less is known about the anabolic part of the syndrome or presumptive deficiencies in bone regeneration. Multipotent mesenchymal stem cells (MSC) are the primary source of osteogenic regeneration. In the present study we aimed to unravel whether MSC biology is directly involved in the pathophysiology of the disease and therefore performed microarray analyses of hMSC of elderly patients (79–94 years old) suffering from osteoporosis (hMSC-OP). In comparison to age-matched controls we detected profound changes in the transcriptome in hMSC-OP, e.g. enhanced mRNA expression of known osteoporosis-associated genes (LRP5, RUNX2, COL1A1) and of genes involved in osteoclastogenesis (CSF1, PTH1R), but most notably of genes coding for inhibitors of WNT and BMP signaling, such as Sclerostin and MAB21L2. These candidate genes indicate intrinsic deficiencies in self-renewal and differentiation potential in osteoporotic stem cells. We also compared both hMSC-OP and non-osteoporotic hMSC-old of elderly donors to hMSC of ∼30 years younger donors and found that the transcriptional changes acquired between the sixth and the ninth decade of life differed widely between osteoporotic and non-osteoporotic stem cells. In addition, we compared the osteoporotic transcriptome to long term-cultivated, senescent hMSC and detected some signs for pre-senescence in hMSC-OP.

Our results suggest that in primary osteoporosis the transcriptomes of hMSC populations show distinct signatures and little overlap with non-osteoporotic aging, although we detected some hints for senescence-associated changes. While there are remarkable inter-individual variations as expected for polygenetic diseases, we could identify many susceptibility genes for osteoporosis known from genetic studies. We also found new candidates, e.g. MAB21L2, a novel repressor of BMP-induced transcription. Such transcriptional changes may reflect epigenetic changes, which are part of a specific osteoporosis-associated aging process.