Suggestion of GLYAT gene underlying variation of bone size and body lean mass as revealed by a bivariate genome-wide association study

Bone and muscle, two major tissue types of musculoskeletal system, have strong genetic determination. Abnormality in bone and/or muscle may cause musculoskeletal diseases such as osteoporosis and sarcopenia. Bone size phenotypes (BSPs), such as hip bone size (HBS), appendicular bone size (ABS), are genetically correlated with body lean mass (mainly muscle mass). However, the specific genes shared by these phenotypes are largely unknown. In this study, we aimed to identify the specific genes with pleiotropic effects on BSPs and appendicular lean mass (ALM). We performed a bivariate genome-wide association study (GWAS) by analyzing ~690,000 SNPs in 1,627 unrelated Han Chinese adults (802 males and 825 females) followed by a replication study in 2,286 unrelated US Caucasians (558 males and 1,728 females). We identified 14 interesting single nucleotide polymorphisms (SNPs) that may contribute to variation of both BSPs and ALM, with p values <10^?6 in discovery stage. Among them, the association of three SNPs (rs2507838, rs7116722, and rs11826261) in/near GLYAT (glycine-N-acyltransferase) gene was replicated in US Caucasians, with p values ranging from 1.89 × 10^?3 to 3.71 × 10^?4 for ALM–ABS, from 5.14 × 10^?3 to 1.11 × 10^?2 for ALM–HBS, respectively. Meta-analyses yielded stronger association signals for rs2507838, rs7116722, and rs11826261, with pooled p values of 1.68 × 10^?8, 7.94 × 10^?8, 6.80 × 10^?8 for ALB–ABS and 1.22 × 10^?4, 9.85 × 10^?5, 3.96 × 10^?4 for ALM–HBS, respectively. Haplotype allele ATA based on these three SNPs was also associated with ALM–HBS and ALM–ABS in both discovery and replication samples. Interestingly, GLYAT was previously found to be essential to glucose metabolism and energy metabolism, suggesting the gene’s dual role in both bone development and muscle growth. Our findings, together with the prior biological evidence, suggest the importance of GLYAT gene in co-regulation of bone phenotypes and body lean mass.     

Genome-wide association study identifies HMGN3 locus for spine bone size variation in Chinese

Bone size (BS) is one of the major risk factors for osteoporotic fractures. BS variation is genetically determined to a substantial degree with heritability over 50%, but specific genes underlying variation of BS are still largely unknown. To identify specific genes for BS in Chinese, initial genome-wide association scan (GWAS) study and follow-up replication study were performed. In initial GWAS study, a group of 12 contiguous single-nucleotide polymorphism (SNP)s, which span a region of ~25 kb and locate at the upstream of HMGN3 gene (high-mobility group nucleosomal binding domain 3), achieved moderate association signals for spine BS, with P values ranging from 6.2E-05 to 1.8E-06. In the follow-up replication study, eight of the 12 SNPs were detected suggestive replicate associations with BS in 1,728 unrelated female Caucasians, which have well-known differences from Chinese in ethnic genetic background. The SNPs in the region of HMGN3 gene formed a tightly combined haplotype block in both Chinese and Caucasians. The results suggest that the genomic region containing HMGN3 gene may be associated with spine BS in Chinese.     

Pathway analysis of genome-wide association study for bone mineral density

The aim of this study was to identify the candidate causal single nucleotide polymorphisms (SNPs) and candidate causal mechanisms that contribute to bone mineral density (BMD) and to generate a SNP to gene to pathway hypothesis using an analytical pathway-based approach. We used hip BMD GWAS data of the genotypes of 301,019 SNPs in 5,715 Europeans. ICSNPathway (identify candidate causal SNPs and pathways) analysis was applied to the BMD GWAS dataset. The first stage involved the pre-selection of candidate causal SNPs by linkage disequilibrium analysis and the functional SNP annotation of the most significant SNPs found. The second stage involved the annotation of biological mechanisms for the pre-selected candidate causal SNPs using improved-gene set enrichment analysis. ICSNPathway analysis identified seven candidate SNPs, eight candidate pathways, and seven hypothetical biological mechanisms. Eight pathways are as follows; gamma-hexachlorocyclohexane degradation (nominal p-value < 0.001, false discovery rate (FDR) <0.001), regulation of the smoothened signaling pathway (nominal p-value < 0.001, FDR = 0.016), TACI and BCMA stimulation of B cell immune response (nominal p-value < 0.001, FDR = 0.021), endonuclease activity (nominal p-value = 0.001, FDR = 0,026), regulation of defense response to virus (nominal p-value = 0.001, FDR = 0.028), serine_type_endopeptidase_inhibitor_activity (nominal p-value = 0.001, FDR = 0.044), endoribonuclease activity (nominal p-value = 0.002, FDR = 0.045), and myeloid leukocyte differentiation (nominal p-value = 0.001, FDR = 0.050). The most significant causal pathway was gamma-hexachlorocyclohexane degradation. CYP3A5, PON2, PON3, CMBL, PON1, ALPL, CYP3A43, CYP3A7, ACP6, ACPP, and ALPI (p < 0.05) are involved in the pathway of gamma-hexachlorocyclohexane degradation. Further examination of the gene contents revealed that DBR1, DICER1, EXO1, FEN1, POP1, POP4, RPP30, and RPP38 were involved in 2 of the 8 pathways (p < 0.05). By applying ICSNPathway analysis to BMD GWAS data, we identified seven candidate SNPs and eight pathways involving gamma-hexachlorocyclohexane degradation, which may contribute to low BMD.     

A genome-wide association analysis implicates SOX6 as a candidate gene for wrist bone mass

Osteoporosis is a highly heritable common bone disease leading to fractures that severely impair the life quality of patients. Wrist fractures caused by osteoporosis are largely due to the scarcity of wrist bone mass. Here we report the results of a genome-wide association study (GWAS) of wrist bone mineral density (BMD). We examined ～500000 SNP markers in 1000 unrelated homogeneous Caucasian subjects and found a novel allelic association with wrist BMD at rs11023787 in the SOX6 (SRY (sex determining region Y)-box 6) gene (P=9.00×10(-5)). Subjects carrying the C allele of rs11023787 in SOX6 had significantly higher mean wrist BMD values than those with the T allele (0.485:0.462 g cm(-2) for C allele vs. T allele carriers). For validation, we performed SOX6 association for BMD in an independent Chinese sample and found that SNP rs11023787 was significantly associated with wrist BMD in the Chinese sample (P=6.41×10(-3)). Meta-analyses of the GWAS scan and the replication studies yielded P-values of 5.20×10(-6) for rs11023787. Results of this study, together with the functional relevance of SOX6 in cartilage formation, support the SOX6 gene as an important gene for BMD variation.     

Identification of novel susceptibility Loci for kawasaki disease in a Han chinese population by a genome-wide association study

Kawasaki disease (KD) is an acute systemic vasculitis syndrome that primarily affects infants and young children. Its etiology is unknown; however, epidemiological findings suggest that genetic predisposition underlies disease susceptibility. Taiwan has the third-highest incidence of KD in the world, after Japan and Korea. To investigate novel mechanisms that might predispose individuals to KD, we conducted a genome-wide association study (GWAS) in 250 KD patients and 446 controls in a Han Chinese population residing in Taiwan, and further validated our findings in an independent Han Chinese cohort of 208 cases and 366 controls. The most strongly associated single-nucleotide polymorphisms (SNPs) detected in the joint analysis corresponded to three novel loci. Among these KD-associated SNPs three were close to the COPB2 (coatomer protein complex beta-2 subunit) gene: rs1873668 (p = 9.52×10⁻⁵), rs4243399 (p = 9.93×10⁻⁵), and rs16849083 (p = 9.93×10⁻⁵). We also identified a SNP in the intronic region of the ERAP1 (endoplasmic reticulum amino peptidase 1) gene (rs149481, p_best = 4.61×10⁻⁵). Six SNPs (rs17113284, rs8005468, rs10129255, rs2007467, rs10150241, and rs12590667) clustered in an area containing immunoglobulin heavy chain variable regions genes, with p_best-values between 2.08×10⁻⁵ and 8.93×10⁻⁶, were also identified. This is the first KD GWAS performed in a Han Chinese population. The novel KD candidates we identified have been implicated in T cell receptor signaling, regulation of proinflammatory cytokines, as well as antibody-mediated immune responses. These findings may lead to a better understanding of the underlying molecular pathogenesis of KD.     

Adiponectin concentrations: a genome-wide association study

Adiponectin is associated with obesity and insulin resistance. To date, there has been no genome-wide association study (GWAS) of adiponectin levels in Asians. Here we present a GWAS of a cohort of Korean volunteers. A total of 4,001 subjects were genotyped by using a genome-wide marker panel in a two-stage design (979 subjects initially and 3,022 in a second stage). Another 2,304 subjects were used for follow-up replication studies with selected markers. In the discovery phase, the top SNP associated with mean log adiponectin was rs3865188 in CDH13 on chromosome 16 (p = 1.69 × 10(-15) in the initial sample, p = 6.58 × 10(-39) in the second genome-wide sample, and p = 2.12 × 10(-32) in the replication sample). The meta-analysis p value for rs3865188 in all 6,305 individuals was 2.82 × 10(-83). The association of rs3865188 with high-molecular-weight adiponectin (p = 7.36 × 10(-58)) was even stronger in the third sample. A reporter assay that evaluated the effects of a CDH13 promoter SNP in complete linkage disequilibrium with rs3865188 revealed that the major allele increased expression 2.2-fold. This study clearly shows that genetic variants in CDH13 influence adiponectin levels in Korean adults.     

Population genetic variation in genome-wide gene expression

Evolutionary biologists seek to understand which traits display variation, are heritable, and influence differential reproduction, because such traits respond to natural selection and underlie organic evolution. Selection acts upon individual differences within a population. Whether individual differences within a natural population include variation in gene expression levels has not yet been addressed on a genome-wide scale. Here we use DNA microarray technology for measuring comparative gene expression and a refined statistical analysis for the purpose of comparing gene expression levels in natural isolates of the wine yeast Saccharomyces cerevisiae. A method for the Bayesian analysis of gene expression levels is used to compare four natural isolates of S. cerevisiae from Montalcino, Italy. Widespread variation in amino acid metabolism, sulfur assimilation and processing, and protein degradation-primarily consisting of differences in expression level smaller than a factor of 2-is demonstrated. Genetic variation in gene expression among isolates from a natural population is present on a genomic scale. It remains to be determined what role differential gene expression may play in adaptation to new or changing environments.     

Pathway analysis of a genome-wide association study in schizophrenia

Objective

The aim of this study was to identify the candidate single nucleotide polymorphisms (SNPs) and candidate mechanisms that contribute to schizophrenia susceptibility and to generate a SNP to gene to pathway hypothesis using an analytical pathway-based approach.

Methods

We used schizophrenia GWAS data of the genotypes of 660,259 SNPs in 1378 controls and 1351 cases of European descent after quality control filtering. ICSNPathway (Identify candidate Causal SNPs and Pathways) analysis was applied to the schizophrenia GWAS dataset. The first stage involved the pre-selection of candidate SNPs by linkage disequilibrium analysis and the functional SNP annotation of the most significant SNPs found. The second stage involved the annotation of biological mechanisms for the pre-selected candidate SNPs using improved-gene set enrichment analysis.

Results

ICSNPathway analysis identified fifteen candidate SNPs, ten candidate pathways, and nine hypothetical biological mechanisms. The most strongly associated potential pathways were as follows. First, rs1644731 and rs1644730 to RDH8 to estrogen biosynthetic process (p < 0.001, FDR < 0.001). The genes involved in this pathway are RDH8 and HSD3B1 (p < 0.05). All-trans-retinol dehydrogenase (RDH8) is a visual cycle enzyme that reduces all-trans-retinal to all-trans-retinol in the presence of NADPH. The chemical reactions and pathways involved result in the formation of estrogens, which are C18 steroid hormones that can stimulate the development of female sexual characteristics. Second, rs1146031 to ACVR1 to mesoderm formation and activin binding (p < 0.001, FDR = 0.032, 0.034). Two of 15 candidate genes are known genes associated with schizophrenia: KCNQ2 and APOL2. One of the 10 candidate pathways, estrogen biosynthetic process, is known to be associated with schizophrenia (p < 0.001, FDR < 0.001). However, 13 of candidate genes (RDH8, ACVR1, PSMD9, KCNAB1, SLC17A3, ARCN1, COG7, STAB2, LRPAP1, STAB1, CXCL16, COL4A4, EXOSC3) and 9 of candidate pathways were novel.

Conclusion

By applying ICSNPathway analysis to schizophrenia GWAS data, we identified candidate SNPs, genes like KCNQ2 and APOL2 and pathways involving the estrogen biosynthetic process may contribute to schizophrenia susceptibility. Further analyses are needed to validate the results of this analysis.     

A genome-wide association study for racing performances in Thoroughbreds clarifies a candidate region near the MSTN gene

Using 1400 microsatellites, a genome-wide association study (GWAS) was performed to identify genomic regions associated with lifetime earnings and performance ranks, as determined by the Japan Racing Association (JRA). The minimum heritability (h(2) ) was estimated at 7-8% based on the quantitative trait model, suggesting that the racing performance is heritable. Following GWAS with microsatellites, fine mapping led to identification of three SNPs on ECA18, namely, g.65809482T>C (P=1.05E-18), g.65868604G>T (P=6.47E-17), and g.66539967A>G (P=3.35E-14) associated with these performance measures. The haplotype of these SNPs, together with a recently published nearby SNP, g.66493737C>T (P=9.06E-16) in strong linkage disequilibrium, also showed a very clear association with the performance (P<1E-05). The candidate genomic region contained eight genes annotated by ENSEMBL, including the myostatin gene (MSTN). These findings suggest the presence of a gene affecting the racing performance in Thoroughbred racehorses in this region on ECA18.     

Prioritization of SNPs for genome-wide association studies using an interaction model of genetic variation,gene expression,and trait variation

The identification of true causal loci to unravel the statistical evidence of genotype-phenotype correlations and the biological relevance of selected single-nucleotide polymorphisms (SNPs) is a challenging issue in genome-wide association studies (GWAS). Here, we introduced a novel method for the prioritization of SNPs based on p-values from GWAS. The method uses functional evidence from populations, including phenotype-associated gene expressions. Based on the concept of genetic interactions, such as perturbation of gene expression by genetic variation, phenotype and gene expression related SNPs were prioritized by adjusting the p-values of SNPs. We applied our method to GWAS data related to drug-induced cytotoxicity. Then, we prioritized loci that potentially play a role in druginduced cytotoxicity. By generating an interaction model, our approach allowed us not only to identify causal loci, but also to find intermediate nodes that regulate the flow of information among causal loci, perturbed gene expression, and resulting phenotypic variation.     

肺结核全基因组关联研究进展

肺结核是由结核分枝杆菌感染引起的一类古老但仍对人类造成巨大影响的传染性疾病。到目前为止, 肺结核依然是由单一病原菌导致死亡人数最多的疾病, 并且随着耐药菌株的出现而呈现死灰复燃之势。近几年, 肺结核全基因组关联研究在世界范围内取得了阶段性成果, 发现了与肺结核相关联的遗传易感位点和区域, 使肺结核的遗传学研究进入了一个崭新的阶段, 为后续肺结核的早期和综合防治提供了重要线索。然而, 由于人群遗传结构差异和宿主/病原体相互作用, 与其他复杂疾病相比, 肺结核全基因组关联研究依旧面临重重困难, 进展缓慢。文章对不同人群肺结核全基因组关联研究及其验证进行综述, 并系统阐述了目前研究中存在的困难及可能的应对策略。     

Tools for efficient epistasis detection in genome-wide association study

Background  

Genome-wide association study (GWAS) aims to find genetic factors underlying complex phenotypic traits, for which epistasis or gene-gene interaction detection is often preferred over single-locus approach. However, the computational burden has been a major hurdle to apply epistasis test in the genome-wide scale due to a large number of single nucleotide polymorphism (SNP) pairs to be tested.     

Functional inactivation of the genome-wide association study obesity gene neuronal growth regulator 1 in mice causes a body mass phenotype

To date, genome-wide association studies (GWAS) have identified at least 32 novel loci for obesity and body mass-related traits. However, the causal genetic variant and molecular mechanisms of specific susceptibility genes in relation to obesity are yet to be fully confirmed and characterised. Here, we examined whether the candidate gene NEGR1 encoding the neuronal growth regulator 1, also termed neurotractin or Kilon, accounts for the obesity association. To characterise the function of NEGR1 for body weight control in vivo, we generated two novel mutant mouse lines, including a constitutive NEGR1-deficient mouse line as well as an ENU-mutagenised line carrying a loss-of-function mutation (Negr1-I87N) and performed metabolic phenotypic analyses. Ablation of NEGR1 results in a small but steady reduction of body mass in both mutant lines, accompanied with a small reduction in body length in the Negr1-I87N mutants. Magnetic resonance scanning reveals that the reduction of body mass in Negr1-I87N mice is due to a reduced proportion of lean mass. Negr1-I87N mutants display reduced food intake and physical activity while normalised energy expenditure remains unchanged. Expression analyses confirmed the brain-specific distribution of NEGR1 including strong expression in the hypothalamus. In vitro assays show that NEGR1 promotes cell-cell adhesion and neurite growth of hypothalamic neurons. Our results indicate a role of NEGR1 in the control of body weight and food intake. This study provides evidence that supports the link of the GWAS candidate gene NEGR1 with body weight control.     

全基因组关联研究现状

在过去的5年中, 全基因组关联研究(Genome-wide association study, GWAS)方法已被证明是研究复杂疾病和性状遗传易感变异的一种有效手段。目前, 各国科学家在多种复杂疾病和性状中开展了大量的GWAS, 对肿瘤、糖尿病、心脏病、神经精神疾病、自身免疫及免疫相关疾病等复杂疾病以及一些常见性状(如身高、体重、血脂、色素等)的遗传易感基因研究取得了重大成果。截止到2010年9月11日, 运用GWAS开展了对近200种复杂疾病/性状的研究, 发现了3 000多个疾病相关的遗传变异。文章就GWAS的发展及其在复杂疾病/性状中的应用做一综述。     

CONAN: copy number variation analysis software for genome-wide association studies

Background  

Genome-wide association studies (GWAS) based on single nucleotide polymorphisms (SNPs) revolutionized our perception of the genetic regulation of complex traits and diseases. Copy number variations (CNVs) promise to shed additional light on the genetic basis of monogenic as well as complex diseases and phenotypes. Indeed, the number of detected associations between CNVs and certain phenotypes are constantly increasing. However, while several software packages support the determination of CNVs from SNP chip data, the downstream statistical inference of CNV-phenotype associations is still subject to complicated and inefficient in-house solutions, thus strongly limiting the performance of GWAS based on CNVs.     

Identification of potential mechanism and hub genes for neuropathic pain by expression-based genome-wide association study

Neuropathic pain (NP) is a common pathological pain state with limited effective treatments. This study was designed to identify potential mechanisms and candidate genes using gene expression–based genome-wide association study (eGWAS). All NP-related microarray experiments were obtained from Gene Expression Omnibus and ArrayExpress. Significantly dysregulated genes were identified between experimental and untreated groups, and the number of microarray experiments in which each gene was dysregulated was calculated. Significantly dysregulated genes were ranked according to P values of the chi-square test. Using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes database, we performed functional and pathway enrichment analysis. Protein-protein interaction (PPI) network and module analysis was performed using Cytoscape software. A total of 115 candidate genes were identified from 19 independent microarray experiments by eGWAS based on the Bonferroni threshold ( P < 2.97 × 10 ⁻⁶). Immune and inflammatory responses, and complement and coagulation cascades, were respectively the most enriched biological process and pathways for candidate genes. The hub genes with highest connectivity in PPI network and two modules Ccl2 and Jun, and Ctss application of the eGWAS methodology can identify mechanisms and candidate genes associated with NP. Our results support the validity and prevalence of inflammatory and immune mechanisms across different NP models, and Ccl2, Jun, and Ctss may be the hub genes for NP.