Genome-wide association scan for stature in Chinese: evidence for ethnic specific loci

In Caucasian, several studies have identified some common variants associated with human stature variation. However, no such study was performed in Chinese, which is the largest population in the world and evidently differs from Caucasian in genetic background. To identify common or ethnic specific genes for stature in Chinese, an initial GWAS and follow-up replication study were performed. Our initial GWAS study found that a group of 13 contiguous SNPs, which span a region of ~150 kb containing two neighboring genes, zinc finger protein (ZNP) 510 and ZNP782, achieved strong signals for association with stature, with P values ranging from 9.71 × 10⁻⁵ to 3.11 × 10⁻⁶. After false discovery rate correction for multiple testing, 9 of the 13 SNPs remain significant (FDR q = 0.036−0.046). The follow-up replication study in an independent 2,953 unrelated southern Chinese confirmed the association of rs10816533 with stature (P = 0.029). All the13 SNPs were in consistently strong linkage disequilibrium (D′ > 0.99) and formed a single perfect haplotype block. The minor allele frequencies for the 13 contiguous SNPs have evidently ethnic difference, which range from 0.21 to 0.33 in Chinese but have as low as ~0.017 reported in dbSNP database in Caucasian. The present results suggest that the genomic region containing the ZNP510 and ZNP782 genes is an ethnic specific locus associated with stature variation in Chinese.     

Identification of PLCL1 gene for hip bone size variation in females in a genome-wide association study

Osteoporosis, the most prevalent metabolic bone disease among older people, increases risk for low trauma hip fractures (HF) that are associated with high morbidity and mortality. Hip bone size (BS) has been identified as one of the key measurable risk factors for HF. Although hip BS is highly genetically determined, genetic factors underlying the trait are still poorly defined. Here, we performed the first genome-wide association study (GWAS) of hip BS interrogating approximately 380,000 SNPs on the Affymetrix platform in 1,000 homogeneous unrelated Caucasian subjects, including 501 females and 499 males. We identified a gene, PLCL1 (phospholipase c-like 1), that had four SNPs associated with hip BS at, or approaching, a genome-wide significance level in our female subjects; the most significant SNP, rs7595412, achieved a p value of 3.72x10(-7). The gene's importance to hip BS was replicated using the Illumina genotyping platform in an independent UK cohort containing 1,216 Caucasian females. Two SNPs of the PLCL1 gene, rs892515 and rs9789480, surrounded by the four SNPs identified in our GWAS, achieved p values of 8.62x10(-3) and 2.44x10(-3), respectively, for association with hip BS. Imputation analyses on our GWAS and the UK samples further confirmed the replication signals; eight SNPs of the gene achieved combined imputed p values<10(-5) in the two samples. The PLCL1 gene's relevance to HF was also observed in a Chinese sample containing 403 females, including 266 with HF and 177 control subjects. A SNP of the PLCL1 gene, rs3771362 that is only approximately 0.6 kb apart from the most significant SNP detected in our GWAS (rs7595412), achieved a p value of 7.66x10(-3) (odds ratio = 0.26) for association with HF. Additional biological support for the role of PLCL1 in BS comes from previous demonstrations that the PLCL1 protein inhibits IP3 (inositol 1,4,5-trisphosphate)-mediated calcium signaling, an important pathway regulating mechanical sensing of bone cells. Our findings suggest that PLCL1 is a novel gene associated with variation in hip BS, and provide new insights into the pathogenesis of HF.     

Copy number variation in CNP267 region may be associated with hip bone size

Osteoporotic hip fracture (HF) is a serious global public health problem associated with high morbidity and mortality. Hip bone size (BS) has been identified as one of key measurable risk factors for HF, independent of bone mineral density (BMD). Hip BS is highly genetically determined, but genetic factors underlying BS variation are still poorly defined. Here, we performed an initial genome-wide copy number variation (CNV) association analysis for hip BS in 1,627 Chinese Han subjects using Affymetrix GeneChip Human Mapping SNP 6.0 Array and a follow-up replicate study in 2,286 unrelated US Caucasians sample. We found that a copy number polymorphism (CNP267) located at chromosome 2q12.2 was significantly associated with hip BS in both initial Chinese and replicate Caucasian samples with p values of 4.73E-03 and 5.66E-03, respectively. An important candidate gene, four and a half LIM domains 2 (FHL2), was detected at the downstream of CNP267, which plays important roles in bone metabolism by binding to several bone formation regulator, such as insulin-like growth factor-binding protein 5 (IGFBP-5) and androgen receptor (AR). Our findings suggest that CNP267 region may be associated with hip BS which might influence the FHL2 gene downstream.     

Genome-Wide Association Analyses Identify SPOCK as a Key Novel Gene Underlying Age at Menarche

For females, menarche is a most significant physiological event. Age at menarche (AAM) is a trait with high genetic determination and is associated with major complex diseases in women. However, specific genes for AAM variation are largely unknown. To identify genetic factors underlying AAM variation, a genome-wide association study (GWAS) examining about 380,000 SNPs was conducted in 477 Caucasian women. A follow-up replication study was performed to validate our major GWAS findings using two independent Caucasian cohorts with 854 siblings and 762 unrelated subjects, respectively, and one Chinese cohort of 1,387 unrelated subjects—all females. Our GWAS identified a novel gene, SPOCK (Sparc/Osteonectin, CWCV, and Kazal-like domains proteoglycan), which had seven SNPs associated with AAM with genome-wide false discovery rate (FDR) q<0.05. Six most significant SNPs of the gene were selected for validation in three independent replication cohorts. All of the six SNPs were replicated in at least one cohort. In particular, SNPs rs13357391 and rs1859345 were replicated both within and across different ethnic groups in all three cohorts, with p values of 5.09×10⁻³ and 4.37×10⁻³, respectively, in the Chinese cohort and combined p values (obtained by Fisher's method) of 5.19×10⁻⁵ and 1.02×10⁻⁴, respectively, in all three replication cohorts. Interestingly, SPOCK can inhibit activation of MMP-2 (matrix metalloproteinase-2), a key factor promoting endometrial menstrual breakdown and onset of menstrual bleeding. Our findings, together with the functional relevance, strongly supported that the SPOCK gene underlies variation of AAM.     

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.     

Meta-analysis of two Chinese populations identifies an autoimmune disease risk allele in 22q11.21 as associated with systemic lupus erythematosus

IntroductionSystemic lupus erythematosus (SLE) is a heterogeneous disease with a diverse spectrum of clinical symptoms, ranging from skin rash to end-organ damage. 22q11.21 has been identified as a susceptibility region for several autoimmune diseases, including SLE. However, detailed information for SLE association and the underlying functional mechanism(s) is still lacking.MethodsThrough meta-analysis of two genome-wide association studies (GWAS) on Han Chinese populations, comprising a total of 1,659 cases and 3,398 controls matched geographically, we closely examined the 22q11.21 region, especially on the reported single-nucleotide polymorphisms (SNPs) associated with different autoimmune diseases and their relationships. We further replicated the most significant associations of SNPs with SLE using 2,612 cases and 2,323 controls of Asian ancestry.ResultsAll reported SNPs in the 22q11.21 region with different autoimmune diseases were examined using the two GWAS data and meta-analysis results, and supportive evidence of association with SLE was found (meta-analysis: P_meta ≤ 7.27E-05), which might require further investigation. SNP rs2298428 was identified as the most significant SNP associated with SLE in this region (P_meta =2.70E-09). It showed independent effects through both stepwise and conditional logistic regression, and there is no evidence of other independent association signals for SLE in this region. The association of rs2298428 was further replicated in three cohorts from Hong Kong, Anhui and Thailand comprising a total of 2,612 cases and 2,323 controls (joint analysis of GWAS and replication result: P_all =1.31E-11, odds ratio =1.23). SNP rs2298428 was shown to be an expression quantitative locus for UBE2L3 gene in different cell types, with the risk allele (T) being correlated with higher expression of UBE2L3. This is consistent with earlier reports on higher expression of UBE2L3 in patients with SLE.ConclusionsAssociation with distinct autoimmune diseases highlights the significance of this region in autoreactive responses and potentially shared functional mechanisms in these diseases.

Electronic supplementary material

The online version of this article (doi:10.1186/s13075-015-0577-6) contains supplementary material, which is available to authorized users.     

PHACTR3基因与中国汉族和美国高加索人群骨密度的关联研究

研究组前期的全基因组关联研究发现PHACTR3基因与骨折关联,为了检测该基因与骨密度的关联关系,采用精细定位关联研究来检测PHACTR3基因内及其附近的SNPs与骨密度的关系。首先在中国样本（1627个不相关的汉族样本）和美国样本（2286个不相关高加索样本）中对PHACTR3基因的140个SNPs进行基因分型,然后采用Plink软件检测PHACTR3基因与腰椎和髋部骨密度的关联关系。发现研究组以前报道的与骨折关联的SNPs rs1555364和rs6064822与腰椎和髋部骨密度关联（P=4．89×10^-2-1．26×10^-2）。另外还发现位于PHACTR3基因内含子中3个SNPs位点（rs6027138,rs1182531和rs1182532）与中国人群和白人腰椎骨密度均显著关联,将中国人与白人样本合并起来进行荟萃分析（Meta—analysis）,得到合并P值为1．40×10^-3到4．00×10^-4,另外发现rs6064820与髋部BMD相关联,合并P值为6．70×10^-3。本研究进一步证实了PHACTR3基因在骨密度变异中的作用,对骨质疏松发病机制的认识提供了新的理论依据。     

A Post-GWAS Replication Study Confirming the PTK2 Gene Associated with Milk Production Traits in Chinese Holstein

Our initial genome-wide association study (GWAS) demonstrated that two SNPs (ARS-BFGL-NGS-33248, UA-IFASA-9288) within the protein tyrosine kinase 2 (PTK2) gene were significantly associated with milk production traits in Chinese Holstein dairy cattle. To further validate if the statistical evidence provided in GWAS were true-positive findings, a replication study was performed herein through genotype-phenotype associations. The two tested SNPs were found to show significant associations with milk production traits, which confirmed the associations observed in the original study. Specifically, SNPs lying in the PTK2 gene were also detected by sequencing 14 unrelated sires in Chinese Holsteins and a total of thirty-three novel SNPs were identified. Thirteen out of these identified SNPs were genotyped and tested for association with milk production traits in an independent resource population. After Bonferroni correction for multiple testing, twelve SNPs were statistically significant for more than two milk production traits. Analyses of pairwise D’ measures of linkage disequilibrium (LD) between all SNPs were also explored. Two haplotype blocks were inferred and the association study at haplotype level revealed similar effects on milk production traits. In addition, the RNA expression analyses revealed that a non-synonymous coding SNP (g.4061098T>G) was involved in the regulation of gene expression. Thus the findings presented here provide strong evidence for associations of PTK2 variants with dairy production traits and may be applied in Chinese Holstein breeding program.     

Identifying consensus disease pathways in Parkinson's disease using an integrative systems biology approach

Parkinson's disease (PD) has had six genome-wide association studies (GWAS) conducted as well as several gene expression studies. However, only variants in MAPT and SNCA have been consistently replicated. To improve the utility of these approaches, we applied pathway analyses integrating both GWAS and gene expression. The top 5000 SNPs (p<0.01) from a joint analysis of three existing PD GWAS were identified and each assigned to a gene. For gene expression, rather than the traditional comparison of one anatomical region between sets of patients and controls, we identified differentially expressed genes between adjacent Braak regions in each individual and adjusted using average control expression profiles. Over-represented pathways were calculated using a hyper-geometric statistical comparison. An integrated, systems meta-analysis of the over-represented pathways combined the expression and GWAS results using a Fisher's combined probability test. Four of the top seven pathways from each approach were identical. The top three pathways in the meta-analysis, with their corrected p-values, were axonal guidance (p = 2.8E-07), focal adhesion (p = 7.7E-06) and calcium signaling (p = 2.9E-05). These results support that a systems biology (pathway) approach will provide additional insight into the genetic etiology of PD and that these pathways have both biological and statistical support to be important in PD.     

A genome-wide association study for quantitative traits in schizophrenia in China

Few genome-wide association studies (GWAS) of schizophrenia have included Chinese populations, and verification of positive genetic findings from other ethnic groups is rare in Chinese groups. We used fluid intelligence as the quantitative trait reflecting schizophrenia dysfunction in Chinese populations, and determined the impact of genetic variation on fluid intelligence phenotypic patterns to identify genetic influences in schizophrenia. The study sample comprised 98 patients with schizophrenia and 60 healthy controls. The general fluid intelligence of participants was assessed with Cattell's Culture-Free Intelligence Test (CCFIT). Subjects were genotyped using the Illumina HumanHap 660 beadchip. We identified the methionine sulfoxide reductase A (MSRA) gene on chromosome 8 as having an association with fluid intelligence. However, only CCFIT subtest 1 (series score) demonstrated a significant result for the interaction term using the criteria of the quantitative trait (QT) analysis of 10(-5) for at least three SNPs. There were 15 haplotype blocks of MSRA gene SNPs identified using Haploview 4.2 with solid spine D' > 0.80. The strongest QT interaction was noted in Block 3, with the most common haplotypes being AAACAGCAG and CGCAGAAGA. In conclusion, we report data from a GWAS with quantitative traits design from Chinese first-episode schizophrenia patients and matched controls. Although the gene identified requires confirmation in an independent sample, the MSRA gene located on chromosome 8 was found to be associated with the phenotype of schizophrenia.     

Genome‐wide association study for 13 agronomic traits reveals distribution of superior alleles in bread wheat from the Yellow and Huai Valley of China

Bread wheat is a leading cereal crop worldwide. Limited amount of superior allele loci restricted the progress of molecular improvement in wheat breeding. Here, we revealed new allelic variation distribution for 13 yield‐related traits in series of genome‐wide association studies (GWAS) using the wheat 90K genotyping assay, characterized in 163 bread wheat cultivars. Agronomic traits were investigated in 14 environments at three locations over 3 years. After filtering SNP data sets, GWAS using 20 689 high‐quality SNPs associated 1769 significant loci that explained, on average, ~20% of the phenotypic variation, both detected already reported loci and new promising genomic regions. Of these, repetitive and pleiotropic SNPs on chromosomes 6AS, 6AL, 6BS, 5BL and 7AS were significantly linked to thousand kernel weight, for example BS00021705_51 on 6BS and wsnp_Ex_c32624_41252144 on 6AS, with phenotypic variation explained (PVE) of ~24%, consistently identified in 12 and 13 of the 14 environments, respectively. Kernel length‐related SNPs were mainly identified on chromosomes 7BS, 6AS, 5AL and 5BL. Plant height‐related SNPs on chromosomes 4DS, 6DL, 2DS and 1BL were, respectively, identified in more than 11 environments, with averaged PVE of ~55%. Four SNPs were confirmed to be important genetic loci in two RIL populations. Based on repetivity and PVE, a total of 41 SNP loci possibly played the key role in modulating yield‐related traits of the cultivars surveyed. Distribution of superior alleles at the 41 SNP loci indicated that superior alleles were getting popular with time and modern cultivars had integrated many superior alleles, especially for peduncle length‐ and plant height‐related superior alleles. However, there were still 19 SNP loci showing less than percentages of 50% in modern cultivars, suggesting they should be paid more attention to improve yield‐related traits of cultivars in the Yellow and Huai wheat region. This study could provide useful information for dissection of yield‐related traits and valuable genetic loci for marker‐assisted selection in Chinese wheat breeding programme.     

Replication of Association between Schizophrenia and Chromosome 6p21-6p22.1 Polymorphisms in Chinese Han Population

Chromosome 6p21-p22.1, spanning the extended major histocompatibility complex (MHC) region, is a highly polymorphic, gene-dense region. It has been identified as a susceptibility locus of schizophrenia in Europeans, Japanese, and Chinese. In our previous two-stage genome-wide association study (GWAS), polymorphisms of zinc finger with KRAB and SCAN domains 4 (ZKSCAN4), nuclear factor-κB-activating protein-like (NKAPL), and piggyBac transposable element derived 1 (PGBD1), localized to chromosome 6p21-p22.1, were strongly associated with schizophrenia. To further investigate the association between polymorphisms at this locus and schizophrenia in the Chinese Han population, we selected eight other single-nucleotide polymorphisms (SNPs) distributed in or near these genes for a case-control association study in an independent sample of 902 cases and 1,091 healthy controls in an attempt to replicate the GWAS results. Four of these eight SNPs (rs12214383, rs1150724, rs3800324, and rs1997660) displayed a nominal difference in allele frequencies between the case and control groups. The association between two of these SNPs and schizophrenia were significant even after Bonferroni correction (rs12000: allele A>G, P = 2.50E-04, odds ratio [OR] = 1.27, 95% confidence interval [CI] = 1.12–1.45; rs1150722: allele C>T, P = 4.28E-05, OR = 0.55, 95% CI = 0.41–0.73). Haplotype ATTGACGC, comprising these eight SNPs (rs2235359, rs2185955, rs12214383, rs12000, rs1150724, rs1150722, rs3800324, and rs1997660), was significantly associated with schizophrenia (P = 6.60E-05). We also performed a combined study of this replication sample and the first-stage GWAS sample. The combined study revealed that rs12000 and rs1150722 were still strongly associated with schizophrenia (rs12000: allele G>A, P _combined  = 0.0019, OR = 0.81; rs1150722: allele G>A, P _combined  = 3.00E-04, OR = 0.61). These results support our findings that locus 6p21-p22.1 is significantly associated with schizophrenia in the Chinese Han population and encourage further studies of the functions of these genetic factors.     

Genome-Wide Association Study Implicates Chromosome 9q21.31 as a Susceptibility Locus for Asthma in Mexican Children

Many candidate genes have been studied for asthma, but replication has varied. Novel candidate genes have been identified for various complex diseases using genome-wide association studies (GWASs). We conducted a GWAS in 492 Mexican children with asthma, predominantly atopic by skin prick test, and their parents using the Illumina HumanHap 550 K BeadChip to identify novel genetic variation for childhood asthma. The 520,767 autosomal single nucleotide polymorphisms (SNPs) passing quality control were tested for association with childhood asthma using log-linear regression with a log-additive risk model. Eleven of the most significantly associated GWAS SNPs were tested for replication in an independent study of 177 Mexican case–parent trios with childhood-onset asthma and atopy using log-linear analysis. The chromosome 9q21.31 SNP rs2378383 (p = 7.10×10⁻⁶ in the GWAS), located upstream of transducin-like enhancer of split 4 (TLE4), gave a p-value of 0.03 and the same direction and magnitude of association in the replication study (combined p = 6.79×10⁻⁷). Ancestry analysis on chromosome 9q supported an inverse association between the rs2378383 minor allele (G) and childhood asthma. This work identifies chromosome 9q21.31 as a novel susceptibility locus for childhood asthma in Mexicans. Further, analysis of genome-wide expression data in 51 human tissues from the Novartis Research Foundation showed that median GWAS significance levels for SNPs in genes expressed in the lung differed most significantly from genes not expressed in the lung when compared to 50 other tissues, supporting the biological plausibility of our overall GWAS findings and the multigenic etiology of childhood asthma.     

A knowledge-based weighting framework to boost the power of genome-wide association studies

Background

We are moving to second-wave analysis of genome-wide association studies (GWAS), characterized by comprehensive bioinformatical and statistical evaluation of genetic associations. Existing biological knowledge is very valuable for GWAS, which may help improve their detection power particularly for disease susceptibility loci of moderate effect size. However, a challenging question is how to utilize available resources that are very heterogeneous to quantitatively evaluate the statistic significances.

Methodology/Principal Findings

We present a novel knowledge-based weighting framework to boost power of the GWAS and insightfully strengthen their explorative performance for follow-up replication and deep sequencing. Built upon diverse integrated biological knowledge, this framework directly models both the prior functional information and the association significances emerging from GWAS to optimally highlight single nucleotide polymorphisms (SNPs) for subsequent replication. In the theoretical calculation and computer simulation, it shows great potential to achieve extra over 15% power to identify an association signal of moderate strength or to use hundreds of whole-genome subjects fewer to approach similar power. In a case study on late-onset Alzheimer disease (LOAD) for a proof of principle, it highlighted some genes, which showed positive association with LOAD in previous independent studies, and two important LOAD related pathways. These genes and pathways could be originally ignored due to involved SNPs only having moderate association significance.

Conclusions/Significance

With user-friendly implementation in an open-source Java package, this powerful framework will provide an important complementary solution to identify more true susceptibility loci with modest or even small effect size in current GWAS for complex diseases.     

Meta-analysis of new genome-wide association studies of colorectal cancer risk

Colorectal cancer is the second leading cause of cancer death in developed countries. Genome-wide association studies (GWAS) have successfully identified novel susceptibility loci for colorectal cancer. To follow up on these findings, and try to identify novel colorectal cancer susceptibility loci, we present results for GWAS of colorectal cancer (2,906 cases, 3,416 controls) that have not previously published main associations. Specifically, we calculated odds ratios and 95% confidence intervals using log-additive models for each study. In order to improve our power to detect novel colorectal cancer susceptibility loci, we performed a meta-analysis combining the results across studies. We selected the most statistically significant single nucleotide polymorphisms (SNPs) for replication using ten independent studies (8,161 cases and 9,101 controls). We again used a meta-analysis to summarize results for the replication studies alone, and for a combined analysis of GWAS and replication studies. We measured ten SNPs previously identified in colorectal cancer susceptibility loci and found eight to be associated with colorectal cancer (p value range 0.02 to 1.8?×?10(-8)). When we excluded studies that have previously published on these SNPs, five SNPs remained significant at p?相似文献   

An analysis of single nucleotide polymorphisms of 125 DNA repair genes in the Texas genome-wide association study of lung cancer with a replication for the XRCC4 SNPs

DNA repair genes are important for maintaining genomic stability and limiting carcinogenesis. We analyzed all single nucleotide polymorphisms (SNPs) of 125 DNA repair genes covered by the Illumina HumanHap300 (v1.1) BeadChips in a previously conducted genome-wide association study (GWAS) of 1154 lung cancer cases and 1137 controls and replicated the top-hits of XRCC4 SNPs in an independent set of 597 cases and 611 controls in Texas populations. We found that six of 20 XRCC4 SNPs were associated with a decreased risk of lung cancer with a P-value of 0.01 or lower in the discovery dataset, of which the most significant SNP was rs10040363 (P for allelic test=4.89 x 10??). Moreover, the data in this region allowed us to impute a potentially functional SNP rs2075685 (imputed P for allelic test=1.3 x 10?3). A luciferase reporter assay demonstrated that the rs2075685G>T change in the XRCC4 promoter increased expression of the gene. In the replication study of rs10040363, rs1478486, rs9293329, and rs2075685, however, only rs10040363 achieved a borderline association with a decreased risk of lung cancer in a dominant model (adjusted OR=0.80, 95% CI=0.62-1.03 and P=0.079). In the final combined analysis of both the Texas GWAS discovery and replication datasets, the strength of the association was increased for rs10040363 (adjusted OR=0.77, 95% CI=0.66-0.89, P(dominant)=5 x 10?? and P for trend=5 x 10??) and rs1478486 (adjusted OR=0.82, 95% CI=0.71-0.94, P(dominant)=6 x 10?3 and P for trend=3.5 x 10?3). Finally, we conducted a meta-analysis of these XRCC4 SNPs with available data from published GWA studies of lung cancer with a total of 12,312 cases and 47,921 controls, in which none of these XRCC4 SNPs was associated with lung cancer risk. It appeared that rs2075685, although associated with increased expression of a reporter gene and lung cancer risk in the Texas populations, did not have an effect on lung cancer risk in other populations. This study underscores the importance of replication using published data in larger populations.     

Extensive natural variation for cellular hydrogen peroxide release is genetically controlled

Natural variation in DNA sequence contributes to individual differences in quantitative traits. While multiple studies have shown genetic control over gene expression variation, few additional cellular traits have been investigated. Here, we investigated the natural variation of NADPH oxidase-dependent hydrogen peroxide (H₂O₂ release), which is the joint effect of reactive oxygen species (ROS) production, superoxide metabolism and degradation, and is related to a number of human disorders. We assessed the normal variation of H₂O₂ release in lymphoblastoid cell lines (LCL) in a family-based 3-generation cohort (CEPH-HapMap), and in 3 population-based cohorts (KORA, GenCord, HapMap). Substantial individual variation was observed, 45% of which were associated with heritability in the CEPH-HapMap cohort. We identified 2 genome-wide significant loci of Hsa12 and Hsa15 in genome-wide linkage analysis. Next, we performed genome-wide association study (GWAS) for the combined KORA-GenCord cohorts (n = 279) using enhanced marker resolution by imputation (>1.4 million SNPs). We found 5 significant associations (p<5.00×10−8) and 54 suggestive associations (p<1.00×10−5), one of which confirmed the linked region on Hsa15. To replicate our findings, we performed GWAS using 58 HapMap individuals and ∼2.1 million SNPs. We identified 40 genome-wide significant and 302 suggestive SNPs, and confirmed genome signals on Hsa1, Hsa12, and Hsa15. Genetic loci within 900 kb from the known candidate gene p67phox on Hsa1 were identified in GWAS in both cohorts. We did not find replication of SNPs across all cohorts, but replication within the same genomic region. Finally, a highly significant decrease in H₂O₂ release was observed in Down Syndrome (DS) individuals (p<2.88×10−12). Taken together, our results show strong evidence of genetic control of H₂O₂ in LCL of healthy and DS cohorts and suggest that cellular phenotypes, which themselves are also complex, may be used as proxies for dissection of complex disorders.     

Common Variant rs9939609 in Gene FTO Confers Risk to Polycystic Ovary Syndrome

Background

Fat mass and obesity-associated gene (FTO) has been associated with obesity, especially the common variant rs9939609. Polycystic ovary syndrome (PCOS) is a complex endocrine-metabolic disorder and over 50% of patients are overweight/obese. Thus FTO is a potential candidate gene for PCOS but their relationship is confusing and remains to be clarified in different population with a large sample size.

Method

This study was performed adopting a two-stage design by genotyping SNP rs9939609. The first set comprise of 741 PCOS and 704 control subjects, with data from our previous GWAS. The second phase of replication study was performed among another independent group of 2858 PCOS and 2358 control subjects using TaqMan-MGB probe assay. All subjects are from Han Chinese.

Results

The less meaningful association of FTO rs9939609 and PCOS discovered in GWAS (P = 2.47E-03), was further confirmed in the replication study (P = 1.86E-09). Using meta-analysis, the P-meta value has reached 6.89E-12, over-exceeding the genome-wide association level of 5.00E-8. By combination, the P value was 1.26E-11 and after BMI adjustment it remained significant(P = 1.82E-06). To further elucidate whether this association is resulted from obesity or PCOS per se, the samples were divided into two groups–obese and non-obese PCOS, and the results were still positive in obese group (P obese = 5.81E-05, OR = 1.55), as well as in non-obese PCOS group (P non-obese = 7.06E-04, OR = 1.28).

Conclusion

Variant rs9939609 in FTO is associated with PCOS in Chinese women, not only in obese PCOS subjects, but also in non-obese cases.     

A fine-mapping study of 7 top scoring genes from a GWAS for major depressive disorder

Major depressive disorder (MDD) is a psychiatric disorder that is characterized--amongst others--by persistent depressed mood, loss of interest and pleasure and psychomotor retardation. Environmental circumstances have proven to influence the aetiology of the disease, but MDD also has an estimated 40% heritability, probably with a polygenic background. In 2009, a genome wide association study (GWAS) was performed on the Dutch GAIN-MDD cohort. A non-synonymous coding single nucleotide polymorphism (SNP) rs2522833 in the PCLO gene became only nominally significant after post-hoc analysis with an Australian cohort which used similar ascertainment. The absence of genome-wide significance may be caused by low SNP coverage of genes. To increase SNP coverage to 100% for common variants (m.a.f.>0.1, r(2)>0.8), we selected seven genes from the GAIN-MDD GWAS: PCLO, GZMK, ANPEP, AFAP1L1, ST3GAL6, FGF14 and PTK2B. We genotyped 349 SNPs and obtained the lowest P-value for rs2715147 in PCLO at P?=?6.8E-7. We imputed, filling in missing genotypes, after which rs2715147 and rs2715148 showed the lowest P-value at P?=?1.2E-6. When we created a haplotype of these SNPs together with the non-synonymous coding SNP rs2522833, the P-value decreased to P?=?9.9E-7 but was not genome wide significant. Although our study did not identify a more strongly associated variant, the results for PCLO suggest that the causal variant is in high LD with rs2715147, rs2715148 and rs2522833.     

A genome-wide SNP scan reveals novel loci for egg production and quality traits in white leghorn and brown-egg dwarf layers

Availability of the complete genome sequence as well as high-density SNP genotyping platforms allows genome-wide association studies (GWAS) in chickens. A high-density SNP array containing 57,636 markers was employed herein to identify associated variants underlying egg production and quality traits within two lines of chickens, i.e., White Leghorn and brown-egg dwarf layers. For each individual, age at first egg (AFE), first egg weight (FEW), and number of eggs (EN) from 21 to 56 weeks of age were recorded, and egg quality traits including egg weight (EW), eggshell weight (ESW), yolk weight (YW), eggshell thickness (EST), eggshell strength (ESS), albumen height(AH) and Haugh unit(HU) were measured at 40 and 60 weeks of age. A total of 385 White Leghorn females and 361 brown-egg dwarf dams were selected to be genotyped. The genome-wide scan revealed 8 SNPs showing genome-wise significant (P<1.51E-06, Bonferroni correction) association with egg production and quality traits under the Fisher's combined probability method. Some significant SNPs are located in known genes including GRB14 and GALNT1 that can impact development and function of ovary, but more are located in genes with unclear functions in layers, and need to be studied further. Many chromosome-wise significant SNPs were also detected in this study and some of them are located in previously reported QTL regions. Most of loci detected in this study are novel and the follow-up replication studies may be needed to further confirm the functional significance for these newly identified SNPs.