Genome-wide association study in east Asians identifies novel susceptibility loci for breast cancer

Genetic factors play an important role in the etiology of both sporadic and familial breast cancer. We aimed to discover novel genetic susceptibility loci for breast cancer. We conducted a four-stage genome-wide association study (GWAS) in 19,091 cases and 20,606 controls of East-Asian descent including Chinese, Korean, and Japanese women. After analyzing 690,947 SNPs in 2,918 cases and 2,324 controls, we evaluated 5,365 SNPs for replication in 3,972 cases and 3,852 controls. Ninety-four SNPs were further evaluated in 5,203 cases and 5,138 controls, and finally the top 22 SNPs were investigated in up to 17,423 additional subjects (7,489 cases and 9,934 controls). SNP rs9485372, near the TGF-β activated kinase (TAB2) gene in chromosome 6q25.1, showed a consistent association with breast cancer risk across all four stages, with a P-value of 3.8×10⁻¹² in the combined analysis of all samples. Adjusted odds ratios (95% confidence intervals) were 0.89 (0.85–0.94) and 0.80 (0.75–0.86) for the A/G and A/A genotypes, respectively, compared with the genotype G/G. SNP rs9383951 (P = 1.9×10⁻⁶ from the combined analysis of all samples), located in intron 5 of the ESR1 gene, and SNP rs7107217 (P = 4.6×10⁻⁷), located at 11q24.3, also showed a consistent association in each of the four stages. This study provides strong evidence for a novel breast cancer susceptibility locus represented by rs9485372, near the TAB2 gene (6q25.1), and identifies two possible susceptibility loci located in the ESR1 gene and 11q24.3, respectively.     

A family-based association study identified CYP17 as a candidate gene for obesity susceptibility in Caucasians

The cytochrome P450c17α gene (CYP17) encodes a key biosynthesis enzyme of estrogen, which is critical in regulating adipogenesis and adipocyte development in humans. We therefore hypothesized that CYP17 is a candidate gene for predicting obesity. In order to test this hypothesis, we performed a family-based association test to investigate the relationship between the CYP17 gene and obesity phenotypes in a large sample comprising 1873 subjects from 405 Caucasian nuclear families of European origin recruited by the Osteoporosis Research Center of Creighton University, USA. Both single SNPs and haplotypes were tested for associations with obesity-related phenotypes, including body mass index (BMI) and fat mass. We identified three SNPs to be significantly associated with BMI, including rs3740397, rs6163, and rs619824. We further characterized the linkage disequilibrium structure for CYP17 and found that the whole CYP17 gene was located in a single-linkage disequilibrium block. This block was observed to be significantly associated with BMI. A major haplotype in this block was significantly associated with both BMI and fat mass. In conclusion, we suggest that the CYP17 gene has an effect on obesity in the Caucasian population. Further independent studies will be needed to confirm our findings.     

A genome-wide association study on copy-number variation identifies a 11q11 loss as a candidate susceptibility variant for colorectal cancer

Colorectal cancer (CRC) is a complex disease, and therefore its development is determined by the combination of both environmental factors and genetic variants. Although genome-wide association studies (GWAS) of SNP variation have conveniently identified 20 genetic variants so far, a significant proportion of the observed heritability is yet to be explained. Common copy-number variants (CNVs) are one of the most important genomic sources of variability, and hence a potential source to explain part of this missing genetic fraction. Therefore, we have performed a GWAS on CNVs to explore the relationship between common structural variation and CRC development. Phase 1 of the GWAS consisted of 881 cases and 667 controls from a Spanish cohort. Copy-number status was validated by quantitative PCR for each of those common CNVs potentially associated with CRC in phase I. Subsequently, SNPs were chosen as proxies for the validated CNVs for phase II replication (1,342 Spanish cases and 1,874 Spanish controls). Four common CNVs were found to be associated with CRC and were further replicated in Phase II. Finally, we found that SNP rs1944682, tagging a 11q11 CNV, was nominally associated with CRC susceptibility (p value = 0.039; OR = 1.122). This locus has been previously related to extreme obesity phenotypes, which could suggest a relationship between body weight and CRC susceptibility.     

A novel, functional and replicable risk gene region for alcohol dependence identified by genome-wide association study

Several genome-wide association studies (GWASs) reported tens of risk genes for alcohol dependence, but most of them have not been replicated or confirmed by functional studies. The present study used a GWAS to search for novel, functional and replicable risk gene regions for alcohol dependence. Associations of all top-ranked SNPs identified in a discovery sample of 681 African-American (AA) cases with alcohol dependence and 508 AA controls were retested in a primary replication sample of 1,409 European-American (EA) cases and 1,518 EA controls. The replicable associations were then subjected to secondary replication in a sample of 6,438 Australian family subjects. A functional expression quantitative trait locus (eQTL) analysis of these replicable risk SNPs was followed-up in order to explore their cis-acting regulatory effects on gene expression. We found that within a 90 Mb region around PHF3-PTP4A1 locus in AAs, a linkage disequilibrium (LD) block in PHF3-PTP4A1 formed the only peak associated with alcohol dependence at p<10⁻⁴. Within this block, 30 SNPs associated with alcohol dependence in AAs (1.6×10⁻⁵≤p≤0.050) were replicated in EAs (1.3×10⁻³≤p≤0.038), and 18 of them were also replicated in Australians (1.8×10⁻³≤p≤0.048). Most of these risk SNPs had strong cis-acting regulatory effects on PHF3-PTP4A1 mRNA expression across three HapMap samples. The distributions of −log(p) values for association and functional signals throughout this LD block were highly consistent across AAs, EAs, Australians and three HapMap samples. We conclude that the PHF3-PTP4A1 region appears to harbor a causal locus for alcohol dependence, and proteins encoded by PHF3 and/or PTP4A1 might play a functional role in the disorder.     

A genome-wide association study identifies two novel promising candidate genes affecting Escherichia coli F4ab/F4ac susceptibility in swine

Enterotoxigenic Escherichia coli (ETEC) expressing F4 fimbria is the major pathogenic bacteria causing diarrhoea in neonatal and post-weaning piglets. Previous studies have revealed that the susceptibility to ETEC F4ab/F4ac is an autosomal Mendelian dominant trait and the loci controlling the F4ab/F4ac receptor are located on SSC13q41, between markers SW207 and S0283. To pinpoint these loci and further validate previous findings, we performed a genome-wide association study (GWAS) using a two generation family-based population, consisting of 301 piglets with phenotypes of susceptibility to ETEC F4ab/F4ac by the vitro adhesion test. The DNA of all piglets and their parents was genotyped using the Illumina PorcineSNP60 BeadChip, and 50,972 and 50,483 SNPs were available for F4ab and F4ac susceptibility, respectively, in the association analysis after quality control. In summary, 28 and 18 significant SNPs (p<0.05) were detected associated with F4ab and F4ac susceptibility respectively at genome-wide significance level. From these significant findings, two novel candidate genes, HEG1 and ITGB5, were firstly identified as the most promising genes underlying F4ab/F4ac susceptibility in swine according to their functions and positions. Our findings herein provide a novel evidence for unravelling genetic mechanism of diarrhoea risk in piglets.     

A comparison in association and linkage genome-wide scans for alcoholism susceptibility genes using single-nucleotide polymorphisms

We conducted genome-wide linkage scans using both microsatellite and single-nucleotide polymorphism (SNP) markers. Regions showing the strongest evidence of linkage to alcoholism susceptibility genes were identified. Haplotype analyses using a sliding-window approach for SNPs in these regions were performed. In addition, we performed a genome-wide association scan using SNP data. SNPs in these regions with evidence of association (P 相似文献   

Genome-wide association study using diversity outcross mice identified candidate genes of pancreatic cancer

To understand the genetic causes of pancreatic cancer (PC), we conducted a genome-wide association study (GWAS) using the diversity outbred (DO) mice population to identify susceptibility genes underlying 7,12-dimethylbenzanthraene (DMBA) induced PC. The phenotype studied was the percent PC lesion area in the DO mice population. We genotyped 7851 SNP markers specifically designed for DO mice across the whole mouse genome. Four susceptibility genes with P values exceeding the genome-wide threshold for percent PC lesion area (P < 2.37 × 10⁻⁶) were identified, i.e., Epha4, Gpc5, Kcnj6, Arid1b. The most significant SNP of Gpc5 (UNC140360310) that is associated with PC lesion area in mice also significantly influences the Gpc5 expression, suggesting that this Gpc5 SNP exerts its role in PC through cis-regulating the gene expression of Gpc5. Together, our data supported that Gpc5 as a tumor suppressor gene involved in the etiology of PC.     

DNA microarray analysis of rheumatoid arthritis susceptibility genes identified by genome-wide association studies

In a recent interesting review, Alex Clarke and Timothy Vyse described the genetics of rheumatic disease [1]. In the past several years, genome-wide association studies (GWAS) have led to the identification of six high-risk rheumatoid arthritis (RA) susceptibility genes - namely, CD244, PADI4, SLC22A2, PTPN22, CTLA4, and STAT4 (summarized in [2]). In vitro studies using mutant alleles and cultured cells have revealed the individual upregulation of CD244, PADI4, SLC22A2, and PTPN22 [2-6]; however, studies on the expression of RA susceptibility genes in RA patients are rare. We therefore investigated the expression of the above-mentioned six RA susceptibility genes in 112 RA patients using DNA microarray analysis. This study aims to clarify whether DNA microarray analysis and GWAS produce comparable results with respect to RA susceptibility genes.Total RNA extracted from total peripheral blood cells obtained from 112 RA patients and 45 healthy individuals was used to prepare aminoallyl RNA. As a reference, mixed RNA from 45 healthy individuals was used. The aminoallyl RNA of each individual and the reference was subjected to Cy3 and Cy5 labeling, respectively, and was hybridized with an oligonucleotide-based DNA microarray. The data obtained were analyzed by nonparametric statistical group comparison. The intensities of the noprobe spots were used as the background. The median and standard deviation of the background intensity were calculated. The genes with an intensity value that was less than the median plus 2 standard deviation of the background intensity were identified as null. The Cy3/Cy5 ratios of all spots on the DNA microarray were normalized using the global ratio median. Only gene expression data that were collected from at least 80% of samples from each group were selected for further analysis. The unpaired Mann-Whitney test was used to determine statistically significant differences in the mRNA expression levels between the RA and healthy groups. Statistical significance was set at P < 0.05.The results of our DNA microarray analysis showed that the expressions of four out of the six RA susceptibility genes were significantly higher in RA patients than in healthy individuals (1.0 × 10^-16 to 2.32 × 10^-5) (Table ​(Table1).1). As described above, the upregulation of these four genes (CD244, PADI4, SLC22A2, and PTPN22) has been previously confirmed in in vitro studies. We found, however, that CTLA4 expression levels were similar between the RA and control groups, whereas STAT4 expression was significantly downregulated in the RA group (1.38 × 10^-8). We investigated the expression of other RA susceptibility genes - namely, TRF1/C5 [7], CD40 [8], and CCL21 [8] - and found that their expressions were similar in both groups. The genetic risk factors for RA were recently reported to differ between Caucasian and Asian (Korean) populations [9]. The samples used in our microarray analysis were derived from the same Asian (Japanese) cohort. The expression profiles for these three genes may therefore not be consistent with the profiles determined by GWAS.

Table 1

Candidate genes identified from rheumatoid arthritis genome-wide association studies

A genome-wide association study uncovers novel genomic regions and candidate genes of yield-related traits in upland cotton

Key message

A total of 62 SNPs associated with yield-related traits were identified by a GWAS. Based on significant SNPs, two candidate genes pleiotropically increase lint yield.

Abstract

Improved fibre yield is considered a constant goal of upland cotton (Gossypium hirsutum) breeding worldwide, but the understanding of the genetic basis controlling yield-related traits remains limited. To better decipher the molecular mechanism underlying these traits, we conducted a genome-wide association study to determine candidate loci associated with six yield-related traits in a population of 719 upland cotton germplasm accessions; to accomplish this, we used 10,511 single-nucleotide polymorphisms (SNPs) genotyped by an Illumina CottonSNP63K array. Six traits, including the boll number, boll weight, lint percentage, fruit branch number, seed index and lint index, were assessed in multiple environments; large variation in all phenotypes was detected across accessions. We identified 62 SNP loci that were significantly associated with different traits on chromosomes A07, D03, D05, D09, D10 and D12. A total of 689 candidate genes were screened, and 27 of them contained at least one significant SNP. Furthermore, two genes (Gh_D03G1064 and Gh_D12G2354) that pleiotropically increase lint yield were identified. These identified SNPs and candidate genes provide important insights into the genetic control underlying high yields in G. hirsutum, ultimately facilitating breeding programmes of high-yielding cotton.

     

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.     

A next generation sequencing combined genome-wide association study identifies novel tuberculosis susceptibility loci in Chinese population

The genetic factors of tuberculosis (TB) susceptibility have been widely recognized. Here we performed a two-stage study in 616 TB patients and 709 healthy controls to systematically identify the genetic markers of TB susceptibility. In the discovery stage, we identified 93 single nucleotide polymorphisms (SNPs) and 3 human leucocyte antigen (HLA) class II alleles that had potential associations with TB susceptibility. In the validation stage, we confirmed that 6 nominally significant SNPs, including 2 novel missense variants at RAB17 and DCTN4 (odds ratio (OR) = 1.40, P = 4.98 × 10⁻³ and OR = 2.30, P = 3.17 × 10⁻² respectively), were associated with the predisposition to TB. Moreover, our study found that HLA-II allele DQA1*05:05 (P = 0.0011, OR = 1.44, 95%CI = 1.15–1.77) was a TB susceptibility locus for the first time. This study comprehensively investigated the genetic variants that were associated with TB susceptibility and provided insight into the tuberculosis pathogenesis.     

An approach based on a genome-wide association study reveals candidate loci for narcolepsy

Narcolepsy is a sleep disorder characterized by excessive daytime sleepiness, cataplexy, and a pathological manifestation of rapid eye movement during sleep. Narcoleptic pathogenesis is triggered by both genetic and environmental factors. Recently, development of genome-wide association studies (GWAS) has identified new genetic factors, with many more susceptibility genes yet to be elucidated. Using a new approach that consists of a combination of GWAS and an extensive database search for candidate genes, we picked up 202 candidate genes and performed a replication study in 222 narcoleptic patients and 380 controls. Statistical analysis indicated that six genes, NFATC2, SCP2, CACNA1C, TCRA, POLE, and FAM3D, were associated with narcolepsy (P < 0.001). Some of these associations were further supported by gene expression analyses and an association study in essential hypersomnia (EHS), CNS hypersonia similar to narcolepsy. This novel approach will be applicable to other GWAS in the search of disease-related susceptibility genes.     

Designing candidate gene and genome-wide case-control association studies

This protocol describes how to appropriately design a genetic association case-control study, either focusing on a candidate gene (CG) or region or implementing a genome-wide approach. The steps described involve: (i) defining the case phenotype in adequate detail; (ii) checking the heritability of the disease in question; (iii) considering whether a population-based study is the appropriate design for the research question; (iv) the appropriate selection of controls; (v) sample size calculations and (vi) giving due consideration to whether it is a de novo or replication study. General guidelines are given, as well as specific examples of a CG and a genome-wide association study into type 2 diabetes. Software and websites used in this protocol include the International HapMap Consortium website, Genetic Power Calculator, CaT, and SNPSpD. Running each of the programs takes only a few seconds; the rate-limiting steps involve thinking through the designs and parameters in the disease models.     

Genetic polymorphisms and cancer susceptibility of breast cancer in Korean women

Breast cancer is the most prevalent cancer among women in Western countries, and its prevalence is also increasing in Asia. The major risk factor for breast cancer can be traced to reproductive events that influence the lifetime levels of hormones. However, a large percentage of breast cancer cases cannot, be explained by these risk factors. The identification of susceptibility factors that predispose individuals to breast cancer (for instance, if they are exposed to particular environmental agents) could possibly give further insight into the etiology of this malignancy and provide targets for the future development of therapeutics. The most interesting candidate genes include those that mediate a range of functions. These include carcinogen metabolism, DNA repair, steroid hormone metabolism, signal transduction, and cell cycle control. we conducted a hospital-based case-control study on South Korea to evaluate the potential modifying role of the genetic pollymprphisms of selected low penetrance gens that are involved carcinogen metabolisms (i.e., CYP1A1, CYP2E1, GSTM1/T1/P1, NAT1/2, etc.), estrogen synthesis and metabolism (i.e., CYP19, CYP17, CYP1B1, COMT, ER-alpha, etc.), DNA repair (i.e., XRCC1/3, ERCC2/4, ATM, AGT, etc.), and signal transduction as well as others (i.e., TGF- beta, IGF-1, TNF- beta, IL-1B, IL-1RN, etc.). We also took into account the potential interaction between these and the known risk factors of breast cancer. The results of selected genes will be presented in this mini-review.     

A genome-wide association study follow-up suggests a possible role for PPARG in systemic sclerosis susceptibility

Introduction

A recent genome-wide association study (GWAS) comprising a French cohort of systemic sclerosis (SSc) reported several non-HLA single-nucleotide polymorphisms (SNPs) showing a nominal association in the discovery phase. We aimed to identify previously overlooked susceptibility variants by using a follow-up strategy.

Methods

Sixty-six non-HLA SNPs showing a P value <10^-4 in the discovery phase of the French SSc GWAS were analyzed in the first step of this study, performing a meta-analysis that combined data from the two published SSc GWASs. A total of 2,921 SSc patients and 6,963 healthy controls were included in this first phase. Two SNPs, PPARG rs310746 and CHRNA9 rs6832151, were selected for genotyping in the replication cohort (1,068 SSc patients and 6,762 healthy controls) based on the results of the first step. Genotyping was performed by using TaqMan SNP genotyping assays.

Results

We observed nominal associations for both PPARG rs310746 (P_MH = 1.90 × 10^-6, OR, 1.28) and CHRNA9 rs6832151 (P_MH = 4.30 × 10^-6, OR, 1.17) genetic variants with SSc in the first step of our study. In the replication phase, we observed a trend of association for PPARG rs310746 (P value = 0.066; OR, 1.17). The combined overall Mantel-Haenszel meta-analysis of all the cohorts included in the present study revealed that PPARG rs310746 remained associated with SSc with a nominal non-genome-wide significant P value (P_MH = 5.00 × 10^-7; OR, 1.25). No evidence of association was observed for CHRNA9 rs6832151 either in the replication phase or in the overall pooled analysis.

Conclusion

Our results suggest a role of PPARG gene in the development of SSc.     

A genome-wide association study identified new variants associated with mathematical abilities in Chinese children

Mathematical ability is moderately heritable, and it is a complex trait which can be evaluated in several different categories. A few genetic studies have been published on general mathematical ability. However, no genetic study focused on specific mathematical ability categories. In this study, we separately performed genome-wide association studies on 11 mathematical ability categories in 1146 students from Chinese elementary schools. We identified seven genome-wide significant single nucleotide polymorphisms (SNPs) with strong linkage disequilibrium among each other (all r² > 0.8) associated with mathematical reasoning ability (top SNP: rs34034296, p = 2.01 × 10⁻⁸, nearest gene: CUB and Sushi multiple domains 3, CSMD3). We replicated one SNP (rs133885) from 585 SNPs previously reported to be associated with general mathematical ability associated with division ability in our data (p = 1.053 × 10⁻⁵). In the gene- and gene-set enrichment analysis by MAGMA, we found three significant enrichments of associations with three mathematical ability categories for three genes (LINGO2, OAS1 and HECTD1). We also observed four significant enrichments of associations with four mathematical ability categories for three gene sets. Our results suggest new candidate genetic loci for the genetics of mathematical ability.     

A novel nonsense mutation in the DMP1 gene identified by a genome-wide association study is responsible for inherited rickets in Corriedale sheep

Inherited rickets of Corriedale sheep is characterized by decreased growth rate, thoracic lordosis and angular limb deformities. Previous outcross and backcross studies implicate inheritance as a simple autosomal recessive disorder. A genome wide association study was conducted using the Illumina OvineSNP50 BeadChip on 20 related sheep comprising 17 affected and 3 carriers. A homozygous region of 125 consecutive single-nucleotide polymorphism (SNP) loci was identified in all affected sheep, covering a region of 6 Mb on ovine chromosome 6. Among 35 candidate genes in this region, the dentin matrix protein 1 gene (DMP1) was sequenced to reveal a nonsense mutation 250C/T on exon 6. This mutation introduced a stop codon (R145X) and could truncate C-terminal amino acids. Genotyping by PCR-RFLP for this mutation showed all 17 affected sheep were "T T" genotypes; the 3 carriers were "C T"; 24 phenotypically normal related sheep were either "C T" or "C C"; and 46 unrelated normal control sheep from other breeds were all "C C". The other SNPs in DMP1 were not concordant with the disease and can all be ruled out as candidates. Previous research has shown that mutations in the DMP1 gene are responsible for autosomal recessive hypophosphatemic rickets in humans. Dmp1_knockout mice exhibit rickets phenotypes. We believe the R145X mutation to be responsible for the inherited rickets found in Corriedale sheep. A simple diagnostic test can be designed to identify carriers with the defective "T" allele. Affected sheep could be used as animal models for this form of human rickets, and for further investigation of the role of DMP1 in phosphate homeostasis.