Yongyong Shi

Identification of loci associated with schizophrenia by genome-wide association and follow-up

We carried out a genome-wide association study of schizophrenia (479 cases, 2,937 controls) and tested loci with P < 10−5 in up to 16,726 additional subjects. Of 12 loci followed up, 3 had strong independent support (P < 5 × 10−4), and the overall pattern of replication was unlikely to occur by chance (P = 9 × 10−8). Meta-analysis provided strongest evidence for association around ZNF804A (P = 1.61 × 10−7) and this strengthened when the affected phenotype included bipolar disorder (P = 9.96 × 10−9).

A partition-ligation-combination-subdivision EM algorithm for haplotype inference with multiallelic markers: update of the SHEsis (http://analysis.bio-x.cn)

Haplotypic information in diploid organisms provides valuable information on human evolutionary history and plays an important role in identifying a candidate gene in the etiology of complex genetic diseases. However, haplotypes of diploid individuals cannot be acquired easily. Molecular haplotyping methods are very costly and have low throughput, and current genotyping and sequencing methods do not provide information on the linkage phase in diploid organisms. The application of statistical methods to infer the haplotype phase in samples of diploid sequences is a very cost-effective approach. Several computational and statistical methods have been developed for haplotype inference, including Clark’s algorithm [1], the Expectation Maximization (EM) algorithm [2], and Gibbs sampler [3]. Because of its interpretability and stability, the EM algorithm has become one of the most widely used statistical algorithms. However, the standard EM algorithm has several weaknesses, including the inability to handle a large number of markers and convergence to the local optimum. To overcome these problems, various derivative methods have been developed, such as the Partition-Ligation EM (PLEM) algorithm to handle many more linked loci [4], the Optimal Step Length EM (OSLEM) algorithm to accelerate the calculations [5], and the Stochastic EM (SEM) algorithm to deal with missing genotypic data and to avoid convergence to local maxima [6]. However, most packages are intended for use with single-nucleotide polymorphism (SNP) data in a biallelic manner. More and more researchers are analyzing both multiallelic and biallelic markers in the linkage and/or association studies of complex diseases. The analysis of linkage disequilibrium (LD) between multiallelic loci and haplotype inference of many loci (including biand multiallelic markers) present a number of common problems. The major difficulty for the haplotype inference problem npg Cell Research (2009) 19:519-523.

Genome-wide association study identifies susceptibility loci for polycystic ovary syndrome on chromosome 2p16.3, 2p21 and 9q33.3

Polycystic ovary syndrome (PCOS) is a common metabolic disorder in women. To identify causative genes, we conducted a genome-wide association study (GWAS) of PCOS in Han Chinese. The discovery set included 744 PCOS cases and 895 controls; subsequent replications involved two independent cohorts (2,840 PCOS cases and 5,012 controls from northern Han Chinese; 498 cases and 780 controls from southern and central Han Chinese). We identified strong evidence of associations between PCOS and three loci: 2p16.3 (rs13405728; combined P-value by meta-analysis Pmeta = 7.55 × 10−21, odds ratio (OR) 0.71); 2p21 (rs13429458, Pmeta = 1.73 × 10−23, OR 0.67); and 9q33.3 (rs2479106, Pmeta = 8.12 × 10−19, OR 1.34). These findings provide new insight into the pathogenesis of PCOS. Follow-up studies of the candidate genes in these regions are recommended.

A genome-wide association study identifies two new lung cancer susceptibility loci at 13q12.12 and 22q12.2 in Han Chinese

Lung cancer is the leading cause of cancer-related deaths worldwide. To identify genetic factors that modify the risk of lung cancer in individuals of Chinese ancestry, we performed a genome-wide association scan in 5,408 subjects (2,331 individuals with lung cancer (cases) and 3,077 controls) followed by a two-stage validation among 12,722 subjects (6,313 cases and 6,409 controls). The combined analyses identified six well-replicated SNPs with independent effects and significant lung cancer associations (P < 5.0 × 10−8) located in TP63 (rs4488809 at 3q28, P = 7.2 × 10−26), TERT-CLPTM1L (rs465498 and rs2736100 at 5p15.33, P = 1.2 × 10−20 and P = 1.0 × 10−27, respectively), MIPEP-TNFRSF19 (rs753955 at 13q12.12, P = 1.5 × 10−12) and MTMR3-HORMAD2-LIF (rs17728461 and rs36600 at 22q12.2, P = 1.1 × 10−11 and P = 6.2 × 10−13, respectively). Two of these loci (13q12.12 and 22q12.2) were newly identified in the Chinese population. These results suggest that genetic variants in 3q28, 5p15.33, 13q12.12 and 22q12.2 may contribute to the susceptibility of lung cancer in Han Chinese.

Genome-wide association study in Han Chinese identifies four new susceptibility loci for coronary artery disease.

We performed a meta-analysis of 2 genome-wide association studies of coronary artery disease comprising 1,515 cases and 5,019 controls followed by replication studies in 15,460 cases and 11,472 controls, all of Chinese Han ancestry. We identify four new loci for coronary artery disease that reached the threshold of genome-wide significance (P < 5 × 10−8). These loci mapped in or near TTC32-WDR35, GUCY1A3, C6orf10-BTNL2 and ATP2B1. We also replicated four loci previously identified in European populations (in or near PHACTR1, TCF21, CDKN2A-CDKN2B and C12orf51). These findings provide new insights into pathways contributing to the susceptibility for coronary artery disease in the Chinese Han population.

Common variants on 8p12 and 1q24.2 confer risk of schizophrenia

Schizophrenia is a severe mental disorder affecting ∼1% of the world population, with heritability of up to 80%. To identify new common genetic risk factors, we performed a genome-wide association study (GWAS) in the Han Chinese population. The discovery sample set consisted of 3,750 individuals with schizophrenia and 6,468 healthy controls (1,578 cases and 1,592 controls from northern Han Chinese, 1,238 cases and 2,856 controls from central Han Chinese, and 934 cases and 2,020 controls from the southern Han Chinese). We further analyzed the strongest association signals in an additional independent cohort of 4,383 cases and 4,539 controls from the Han Chinese population. Meta-analysis identified common SNPs that associated with schizophrenia with genome-wide significance on 8p12 (rs16887244, P = 1.27 × 10−10) and 1q24.2 (rs10489202, P = 9.50 × 10−9). Our findings provide new insights into the pathogenesis of schizophrenia.

A genome-wide association study identifies new susceptibility loci for non-cardia gastric cancer at 3q13.31 and 5p13.1

Gastric cancer, including the cardia and non-cardia types, is the second leading cause of cancer-related deaths worldwide. To identify genetic risk variants for non-cardia gastric cancer, we performed a genome-wide association study in 3,279 individuals (1,006 with non-cardia gastric cancer and 2,273 controls) of Chinese descent. We replicated significant associations in an additional 6,897 subjects (3,288 with non-cardia gastric cancer and 3,609 controls). We identified two new susceptibility loci for non-cardia gastric cancer at 5p13.1 (rs13361707 in the region including PTGER4 and PRKAA1; odds ratio (OR) = 1.41; P = 7.6 × 10−29) and 3q13.31 (rs9841504 in ZBTB20; OR = 0.76; P = 1.7 × 10−9). Imputation analyses also confirmed previously reported associations of rs2294008 and rs2976392 on 8q24, rs4072037 on 1q22 and rs13042395 on 20p13 with non-cardia gastric cancer susceptibility in the Han Chinese population.

GWAS Identifies Novel Susceptibility Loci on 6p21.32 and 21q21.3 for Hepatocellular Carcinoma in Chronic Hepatitis B Virus Carriers

Genome-wide association studies (GWAS) have recently identified KIF1B as susceptibility locus for hepatitis B virus (HBV)–related hepatocellular carcinoma (HCC). To further identify novel susceptibility loci associated with HBV–related HCC and replicate the previously reported association, we performed a large three-stage GWAS in the Han Chinese population. 523,663 autosomal SNPs in 1,538 HBV–positive HCC patients and 1,465 chronic HBV carriers were genotyped for the discovery stage. Top candidate SNPs were genotyped in the initial validation samples of 2,112 HBV–positive HCC cases and 2,208 HBV carriers and then in the second validation samples of 1,021 cases and 1,491 HBV carriers. We discovered two novel associations at rs9272105 (HLA-DQA1/DRB1) on 6p21.32 (OR = 1.30, P = 1.13×10−19) and rs455804 (GRIK1) on 21q21.3 (OR = 0.84, P = 1.86×10−8), which were further replicated in the fourth independent sample of 1,298 cases and 1,026 controls (rs9272105: OR = 1.25, P = 1.71×10−4; rs455804: OR = 0.84, P = 6.92×10−3). We also revealed the associations of HLA-DRB1*0405 and 0901*0602, which could partially account for the association at rs9272105. The association at rs455804 implicates GRIK1 as a novel susceptibility gene for HBV–related HCC, suggesting the involvement of glutamate signaling in the development of HBV–related HCC.

A case control and family based association study of the neuregulin1 gene and schizophrenia

Data from twin, family, and adoption studies provide strong evidence that genetic factors play a major aetiological role in schizophrenia. By a series of linkage studies, chromosome 8p has been implicated as a region harbouring a schizophrenia susceptibility gene.1–4 Recently, Stefansson and colleagues reported that neuregulin 1 (NRG1 ), located in 8p21-12, may be involved in the aetiology of schizophrenia.4,5 In their linkage and association studies, a 290 kb core at risk haplotype at the 5′ end of NRG1 was found to be strongly associated with schizophrenia in Icelandic and Scottish populations. This haplotype contains the first exon of NRG1 , which encodes a part of glial growth factor 2 ( GGF2 ). Deficiency of glial growth factors has been presumed to be implicated in the pathogenesis of schizophrenia.6 Futhermore, NRG1 mutant mice have fewer functional N-methyl D-aspartate(NMDA) receptors than wild type mice, and display stereotypic behavioural abnormalities similar to those of normal mice treated with the psychogenic drug phenylcyclidine.4 This core at risk haplotype was defined by five single nucleotide polymorphisms (SNP8NRG221132, SNP8NRG221533, SNP8NRG241930, SNP8NRG243177, SNP8NRG433E1006) and two microsatellites (478B14-848, 420M91395). The frequency of this haplotype in schizophrenic individuals was higher than in controls; in Icelandic samples the frequency was 15.4 (7.5%; p  =  0.000087).4 The first replication using Scottish samples revealed a similar result at 10.2 (5.9%; p  =  0.00031).5 Another replication performed by Williams et al with British or Irish samples used one SNP and the two microsatellites of the core at risk haplotype. However, the association was much weaker at 9.5 (7.5%; p  =  0.04).7 Yang et al reported other markers located in the middle of NRG1 and associated with schizophrenia, in a Chinese population.8 Another independent analysis using 13 microsatellites found two groups of haplotypes, which were significantly …

Association analyses identify multiple new lung cancer susceptibility loci and their interactions with smoking in the Chinese population

To find additional susceptibility loci for lung cancer, we tested promising associations from our previous genome-wide association study (GWAS) of lung cancer in the Chinese population in an extended validation sample size of 7,436 individuals with lung cancer (cases) and 7,483 controls. We found genome-wide significant (P < 5.0 × 10−8) evidence for three additional lung cancer susceptibility loci at 10p14 (rs1663689, close to GATA3, P = 2.84 × 10−10), 5q32 (rs2895680 in PPP2R2B-STK32A-DPYSL3, P = 6.60 × 10−9) and 20q13.2 (rs4809957 in CYP24A1, P = 1.20 × 10−8). We also found consistent associations for rs247008 at 5q31.1 (IL3-CSF2-P4HA2, P = 7.68 × 10−8) and rs9439519 at 1p36.32 (AJAP1-NPHP4, P = 3.65 × 10−6). Four of these loci showed evidence for interactions with smoking dose (P = 1.72 × 10−10, P = 5.07 × 10−3, P = 6.77 × 10−3 and P = 4.49 × 10−2 for rs2895680, rs4809957, rs247008 and rs9439519, respectively). These results advance our understanding of lung cancer susceptibility and highlight potential pathways that integrate genetic variants and smoking in the development of lung cancer.