Genome-wide association study combined with biological context can reveal more disease-related SNPs altering microRNA target seed sites

Background

Emerging studies demonstrate that single nucleotide polymorphisms (SNPs) resided in the microRNA recognition element seed sites (MRESSs) in 3′UTR of mRNAs are putative biomarkers for human diseases and cancers. However, exhaustively experimental validation for the causality of MRESS SNPs is impractical. Therefore bioinformatics have been introduced to predict causal MRESS SNPs. Genome-wide association study (GWAS) provides a way to detect susceptibility of millions of SNPs simultaneously by taking linkage disequilibrium (LD) into account, but the multiple-testing corrections implemented to suppress false positive rate always sacrificed the sensitivity. In our study, we proposed a method to identify candidate causal MRESS SNPs from 12 GWAS datasets without performing multiple-testing corrections. Alternatively, we used biological context to ensure credibility of the selected SNPs.

Results

In 11 out of the 12 GWAS datasets, MRESS SNPs were over-represented in SNPs with p-value ≤ 0.05 (odds ratio (OR) ranged from 1.1 to 2.4). Moreover, host genes of susceptible MRESS SNPs in each of the 11 GWAS dataset shared biological context with reported causal genes. There were 286 MRESS SNPs identified by our method, while only 13 SNPs were identified by multiple-testing corrections with a given threshold of 1 × 10⁻⁵, which is a common cutoff used in GWAS. 27 out of the 286 candidate SNPs have been reported to be deleterious while only 2 out of 13 multiple-testing corrected SNPs were documented in PubMed. MicroRNA-mRNA interactions affected by the 286 candidate SNPs were likely to present negatively correlated expression. These SNPs introduced greater alternation of binding free energy than other MRESS SNPs, especially when grouping by haplotypes (4210 vs. 4105 cal/mol by mean, 9781 vs. 8521 cal/mol by mean, respectively).

Conclusions

MRESS SNPs are promising disease biomarkers in multiple GWAS datasets. The method of integrating GWAS p-value and biological context is stable and effective for selecting candidate causal MRESS SNPs, it reduces the loss of sensitivity compared to multiple-testing corrections. The 286 candidate causal MRESS SNPs provide researchers a credible source to initialize their design of experimental validations in the future.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-669) contains supplementary material, which is available to authorized users.     

Genome-wide association study of body weight in Australian Merino sheep reveals an orthologous region on OAR6 to human and bovine genomic regions affecting height and weight

Background

Body weight (BW) is an important trait for meat production in sheep. Although over the past few years, numerous quantitative trait loci (QTL) have been detected for production traits in cattle, few QTL studies have been reported for sheep, with even fewer on meat production traits. Our objective was to perform a genome-wide association study (GWAS) with the medium-density Illumina Ovine SNP50 BeadChip to identify genomic regions and corresponding haplotypes associated with BW in Australian Merino sheep.

Methods

A total of 1781 Australian Merino sheep were genotyped using the medium-density Illumina Ovine SNP50 BeadChip. Among the 53 862 single nucleotide polymorphisms (SNPs) on this array, 48 640 were used to perform a GWAS using a linear mixed model approach. Genotypes were phased with hsphase; to estimate SNP haplotype effects, linkage disequilibrium blocks were identified in the detected QTL region.

Results

Thirty-nine SNPs were associated with BW at a Bonferroni-corrected genome-wide significance threshold of 1 %. One region on sheep (Ovis aries) chromosome 6 (OAR6) between 36.15 and 38.56 Mb, included 13 significant SNPs that were associated with BW; the most significant SNP was OAR6_41936490.1 (P = 2.37 × 10⁻¹⁶) at 37.69 Mb with an allele substitution effect of 2.12 kg, which corresponds to 0.248 phenotypic standard deviations for BW. The region that surrounds this association signal on OAR6 contains three genes: leucine aminopeptidase 3 (LAP3), which is involved in the processing of the oxytocin precursor; NCAPG non-SMC condensin I complex, subunit G (NCAPG), which is associated with foetal growth and carcass size in cattle; and ligand dependent nuclear receptor corepressor-like (LCORL), which is associated with height in humans and cattle.

Conclusions

The GWAS analysis detected 39 SNPs associated with BW in sheep and a major QTL region was identified on OAR6. In several other mammalian species, regions that are syntenic with this region have been found to be associated with body size traits, which may reflect that the underlying biological mechanisms share a common ancestry. These findings should facilitate the discovery of causative variants for BW and contribute to marker-assisted selection.

Electronic supplementary material

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

A new association test based on disease allele selection for case–control genome-wide association studies

Background

Current robust association tests for case–control genome-wide association study (GWAS) data are mainly based on the assumption of some specific genetic models. Due to the richness of the genetic models, this assumption may not be appropriate. Therefore, robust but powerful association approaches are desirable.

Results

In this paper, we propose a new approach to testing for the association between the genotype and phenotype for case–control GWAS. This method assumes a generalized genetic model and is based on the selected disease allele to obtain a p-value from the more powerful one-sided test. Through a comprehensive simulation study we assess the performance of the new test by comparing it with existing methods. Some real data applications are also used to illustrate the use of the proposed test.

Conclusions

Based on the simulation results and real data application, the proposed test is powerful and robust.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-358) contains supplementary material, which is available to authorized users.     

Association of Genetic Variants with Primary Angle Closure Glaucoma in Two Different Populations

Purpose

A recent large genome-wide association study (GWAS) identified multiple variants associated with primary angle-closure glaucoma (PACG). The present study investigated the role of these variants in two cohorts with PACG recruited from Australia and Nepal.

Method

Patients with PACG and appropriate controls were recruited from eye clinics in Australia (n = 232 cases and n = 288 controls) and Nepal (n = 106 cases and 204 controls). Single nucleotide polymorphisms (SNPs) rs3753841 (COL11A1), rs1015213 (located between PCMTD1 and ST18), rs11024102 (PLEKHA7), and rs3788317 (TXNRD2) were selected and genotyped on the Sequenom. Analyses were conducted using PLINK and METAL.

Results

After adjustment for age and sex, SNP rs3753841 was found to be significantly associated with PACG in the Australian cohort (p = 0.017; OR = 1.34). SNPs rs1015213 (p = 0.014; OR 2.35) and rs11024102 (p = 0.039; OR 1.43) were significantly associated with the disease development in the Nepalese cohort. None of these SNPs survived Bonferroni correction (p = 0.05/4 = 0.013). However, in the combined analysis, of both cohorts, rs3753841 and rs1015213 showed significant association with p-values of 0.009 and 0.004, respectively both surviving Bonferroni correction. SNP rs11024102 showed suggestive association with PACG (p-value 0.035) and no association was found with rs3788317.

Conclusion

The present results support the initial GWAS findings, and confirm the SNP’s contribution to PACG. This is the first study to investigate these loci in both Australian Caucasian and Nepalese populations.     

Association of 25-Hydroxyvitamin D status and genetic variation in the vitamin D metabolic pathway with FEV1 in the Framingham Heart Study

Background

Vitamin D is associated with lung function in cross-sectional studies, and vitamin D inadequacy is hypothesized to play a role in the pathogenesis of chronic obstructive pulmonary disease. Further data are needed to clarify the relation between vitamin D status, genetic variation in vitamin D metabolic genes, and cross-sectional and longitudinal changes in lung function in healthy adults.

Methods

We estimated the association between serum 25-hydroxyvitamin D [25(OH)D] and cross-sectional forced expiratory volume in the first second (FEV₁) in Framingham Heart Study (FHS) Offspring and Third Generation participants and the association between serum 25(OH)D and longitudinal change in FEV₁ in Third Generation participants using linear mixed-effects models. Using a gene-based approach, we investigated the association between 241 SNPs in 6 select vitamin D metabolic genes in relation to longitudinal change in FEV₁ in Offspring participants and pursued replication of these findings in a meta-analyzed set of 4 independent cohorts.

Results

We found a positive cross-sectional association between 25(OH)D and FEV₁ in FHS Offspring and Third Generation participants (P = 0.004). There was little or no association between 25(OH)D and longitudinal change in FEV₁ in Third Generation participants (P = 0.97). In Offspring participants, the CYP2R1 gene, hypothesized to influence usual serum 25(OH)D status, was associated with longitudinal change in FEV₁ (gene-based P < 0.05). The most significantly associated SNP from CYP2R1 had a consistent direction of association with FEV₁ in the meta-analyzed set of replication cohorts, but the association did not reach statistical significance thresholds (P = 0.09).

Conclusions

Serum 25(OH)D status was associated with cross-sectional FEV₁, but not longitudinal change in FEV₁. The inconsistent associations may be driven by differences in the groups studied. CYP2R1 demonstrated a gene-based association with longitudinal change in FEV₁ and is a promising candidate gene for further studies.

Electronic supplementary material

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

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 of production traits in a commercial population of Large White pigs: evidence of haplotypes affecting meat quality

Background

Numerous quantitative trait loci (QTL) have been detected in pigs over the past 20 years using microsatellite markers. However, due to the low density of these markers, the accuracy of QTL location has generally been poor. Since 2009, the dense genome coverage provided by the Illumina PorcineSNP60 BeadChip has made it possible to more accurately map QTL using genome-wide association studies (GWAS). Our objective was to perform high-density GWAS in order to identify genomic regions and corresponding haplotypes associated with production traits in a French Large White population of pigs.

Methods

Animals (385 Large White pigs from 106 sires) were genotyped using the PorcineSNP60 BeadChip and evaluated for 19 traits related to feed intake, growth, carcass composition and meat quality. Of the 64 432 SNPs on the chip, 44 412 were used for GWAS with an animal mixed model that included a regression coefficient for the tested SNPs and a genomic kinship matrix. SNP haplotype effects in QTL regions were then tested for association with phenotypes following phase reconstruction based on the Sscrofa10.2 pig genome assembly.

Results

Twenty-three QTL regions were identified on autosomes and their effects ranged from 0.25 to 0.75 phenotypic standard deviation units for feed intake and feed efficiency (four QTL), carcass (12 QTL) and meat quality traits (seven QTL). The 10 most significant QTL regions had effects on carcass (chromosomes 7, 10, 16, 17 and 18) and meat quality traits (two regions on chromosome 1 and one region on chromosomes 8, 9 and 13). Thirteen of the 23 QTL regions had not been previously described. A haplotype block of 183 kb on chromosome 1 (six SNPs) was identified and displayed three distinct haplotypes with significant (0.0001 < P < 0.03) associations with all evaluated meat quality traits.

Conclusions

GWAS analyses with the PorcineSNP60 BeadChip enabled the detection of 23 QTL regions that affect feed consumption, carcass and meat quality traits in a LW population, of which 13 were novel QTL. The proportionally larger number of QTL found for meat quality traits suggests a specific opportunity for improving these traits in the pig by genomic selection.     

Educational attainment: a genome wide association study in 9538 Australians

Background

Correlations between Educational Attainment (EA) and measures of cognitive performance are as high as 0.8. This makes EA an attractive alternative phenotype for studies wishing to map genes affecting cognition due to the ease of collecting EA data compared to other cognitive phenotypes such as IQ.

Methodology

In an Australian family sample of 9538 individuals we performed a genome-wide association scan (GWAS) using the imputed genotypes of ∼2.4 million single nucleotide polymorphisms (SNP) for a 6-point scale measure of EA. Top hits were checked for replication in an independent sample of 968 individuals. A gene-based test of association was then applied to the GWAS results. Additionally we performed prediction analyses using the GWAS results from our discovery sample to assess the percentage of EA and full scale IQ variance explained by the predicted scores.

Results

The best SNP fell short of having a genome-wide significant p-value (p = 9.77×10⁻⁷). In our independent replication sample six SNPs among the top 50 hits pruned for linkage disequilibrium (r²<0.8) had a p-value<0.05 but only one of these SNPs survived correction for multiple testing - rs7106258 (p = 9.7*10⁻⁴) located in an intergenic region of chromosome 11q14.1. The gene based test results were non-significant and our prediction analyses show that the predicted scores explained little variance in EA in our replication sample.

Conclusion

While we have identified a polymorphism chromosome 11q14.1 associated with EA, further replication is warranted. Overall, the absence of genome-wide significant p-values in our large discovery sample confirmed the high polygenic architecture of EA. Only the assembly of large samples or meta-analytic efforts will be able to assess the implication of common DNA polymorphisms in the etiology of EA.     

Genome-wide association analysis of seedling root development in maize (Zea mays L.)

Background

Plants rely on the root system for anchorage to the ground and the acquisition and absorption of nutrients critical to sustaining productivity. A genome wide association analysis enables one to analyze allelic diversity of complex traits and identify superior alleles. 384 inbred lines from the Ames panel were genotyped with 681,257 single nucleotide polymorphism markers using Genotyping-by-Sequencing technology and 22 seedling root architecture traits were phenotyped.

Results

Utilizing both a general linear model and mixed linear model, a GWAS study was conducted identifying 268 marker trait associations (p ≤ 5.3×10^-7). Analysis of significant SNP markers for multiple traits showed that several were located within gene models with some SNP markers localized within regions of previously identified root quantitative trait loci. Gene model GRMZM2G153722 located on chromosome 4 contained nine significant markers. This predicted gene is expressed in roots and shoots.

Conclusion

This study identifies putatively associated SNP markers associated with root traits at the seedling stage. Some SNPs were located within or near (<1 kb) gene models. These gene models identify possible candidate genes involved in root development at the seedling stage. These and respective linked or functional markers could be targets for breeders for marker assisted selection of seedling root traits.

Electronic supplementary material

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

Genome-wide association and biological pathway analysis for milk-fat composition in Danish Holstein and Danish Jersey cattle

Background

The milk fat profile of the Danish Holstein (DH) and Danish Jersey (DJ) show clear differences. Identification of the genomic regions, genes and biological pathways underlying the milk fat biosynthesis will improve the understanding of the biology underlying bovine milk fat production and may provide new possibilities to change the milk fat composition by selective breeding. In this study a genome wide association scan (GWAS) in the DH and DJ was performed for a detailed milk fatty acid (FA) profile using the HD bovine SNP array and subsequently a biological pathway analysis based on the SNP data was performed.

Results

The GWAS identified in total 1,233 SNPs (FDR < 0.10) spread over 18 chromosomes for nine different FA traits for the DH breed and 1,122 SNPs (FDR < 0.10) spread over 26 chromosomes for 13 different FA traits were detected for the DJ breed. Of these significant SNPs, 108 SNP markers were significant in both DH and DJ (C14-index, BTA26; C16, BTA14; fat percentage (FP), BTA14). This was supported by an enrichment test. The QTL on BTA14 and BTA26 represented the known candidate genes DGAT and SCD. In addition we suggest ACSS3 to be a good candidate gene for the QTL on BTA5 for C10:0 and C15:0. In addition, genetic correlations between the FA traits within breed showed large similarity across breeds. Furthermore, the biological pathway analysis revealed that fat digestion and absorption (KEGG04975) plays a role for the traits FP, C14:1, C16 index and C16:1.

Conclusion

There was a clear similarity between the underlying genetics of FA in the milk between DH and DJ. This was supported by the fact that there was substantial overlap between SNPs for FP, C14 index, C14:1, C16 index and C16:1. In addition genetic correlations between FA showed a similar pattern across DH and DJ. Furthermore the biological pathway analysis suggested that fat digestion and absorption KEGG04975 is important for the traits FP, C14:1, C16 index and C16:1.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1112) contains supplementary material, which is available to authorized users.     

A bi-filtering method for processing single nucleotide polymorphism array data improves the quality of genetic map and accuracy of quantitative trait locus mapping in doubled haploid populations of polyploid Brassica napus

Background

Single nucleotide polymorphism (SNP) markers have a wide range of applications in crop genetics and genomics. Due to their polyploidy nature, many important crops, such as wheat, cotton and rapeseed contain a large amount of repeat and homoeologous sequences in their genomes, which imposes a huge challenge in high-throughput genotyping with sequencing and/or array technologies. Allotetraploid Brassica napus (AACC, 2n = 4x = 38) comprises of two highly homoeologous sub-genomes derived from its progenitor species B. rapa (AA, 2n = 2x = 20) and B. oleracea (CC, 2n = 2x = 18), and is an ideal species to exploit methods for reducing the interference of extensive inter-homoeologue polymorphisms (mHemi-SNPs and Pseudo-simple SNPs) between closely related sub-genomes.

Results

Based on a recent B. napus 6K SNP array, we developed a bi-filtering procedure to identify unauthentic lines in a DH population, and mHemi-SNPs and Pseudo-simple SNPs in an array data matrix. The procedure utilized both monomorphic and polymorphic SNPs in the DH population and could effectively distinguish the mHemi-SNPs and Pseudo-simple SNPs that resulted from superposition of the signals from multiple SNPs. Compared with conventional procedure for array data processing, the bi-filtering method could minimize the pseudo linkage relationship caused by the mHemi-SNPs and Pseudo-simple SNPs, thus improving the quality of SNP genetic map. Furthermore, the improved genetic map could increase the accuracies of mapping of QTLs as demonstrated by the ability to eliminate non-real QTLs in the mapping population.

Conclusions

The bi-filtering analysis of the SNP array data represents a novel approach to effectively assigning the multi-loci SNP genotypes in polyploid B. napus and may find wide applications to SNP analyses in polyploid crops.

Electronic supplementary material

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

Genome-Wide Association Study for Incident Myocardial Infarction and Coronary Heart Disease in Prospective Cohort Studies: The CHARGE Consortium

Background

Data are limited on genome-wide association studies (GWAS) for incident coronary heart disease (CHD). Moreover, it is not known whether genetic variants identified to date also associate with risk of CHD in a prospective setting.

Methods

We performed a two-stage GWAS analysis of incident myocardial infarction (MI) and CHD in a total of 64,297 individuals (including 3898 MI cases, 5465 CHD cases). SNPs that passed an arbitrary threshold of 5×10⁻⁶ in Stage I were taken to Stage II for further discovery. Furthermore, in an analysis of prognosis, we studied whether known SNPs from former GWAS were associated with total mortality in individuals who experienced MI during follow-up.

Results

In Stage I 15 loci passed the threshold of 5×10⁻⁶; 8 loci for MI and 8 loci for CHD, for which one locus overlapped and none were reported in previous GWAS meta-analyses. We took 60 SNPs representing these 15 loci to Stage II of discovery. Four SNPs near QKI showed nominally significant association with MI (p-value<8.8×10⁻³) and three exceeded the genome-wide significance threshold when Stage I and Stage II results were combined (top SNP rs6941513: p = 6.2×10⁻⁹). Despite excellent power, the 9p21 locus SNP (rs1333049) was only modestly associated with MI (HR = 1.09, p-value = 0.02) and marginally with CHD (HR = 1.06, p-value = 0.08). Among an inception cohort of those who experienced MI during follow-up, the risk allele of rs1333049 was associated with a decreased risk of subsequent mortality (HR = 0.90, p-value = 3.2×10⁻³).

Conclusions

QKI represents a novel locus that may serve as a predictor of incident CHD in prospective studies. The association of the 9p21 locus both with increased risk of first myocardial infarction and longer survival after MI highlights the importance of study design in investigating genetic determinants of complex disorders.     

Investigation of a genome wide association signal for obesity: synthetic association and haplotype analyses at the melanocortin 4 receptor gene locus

Background

Independent genome-wide association studies (GWAS) showed an obesogenic effect of two single nucleotide polymorphisms (SNP; rs12970134 and rs17782313) more than 150 kb downstream of the melanocortin 4 receptor gene (MC4R). It is unclear if the SNPs directly influence MC4R function or expression, or if the SNPs are on a haplotype that predisposes to obesity or includes functionally relevant genetic variation (synthetic association). As both exist, functionally relevant mutations and polymorphisms in the MC4R coding region and a robust association downstream of the gene, MC4R is an ideal model to explore synthetic association.

Methodology/Principal Findings

We analyzed a genomic region (364.9 kb) encompassing the MC4R in GWAS data of 424 obesity trios (extremely obese child/adolescent and both parents). SNP rs12970134 showed the lowest p-value (p = 0.004; relative risk for the obesity effect allele: 1.37); conditional analyses on this SNP revealed that 7 of 78 analyzed SNPs provided independent signals (p≤0.05). These 8 SNPs were used to derive two-marker haplotypes. The three best (according to p-value) haplotype combinations were chosen for confirmation in 363 independent obesity trios. The confirmed obesity effect haplotype includes SNPs 3′ and 5′ of the MC4R. Including MC4R coding variants in a joint model had almost no impact on the effect size estimators expected under synthetic association.

Conclusions/Significance

A haplotype reaching from a region 5′ of the MC4R to a region at least 150 kb from the 3′ end of the gene showed a stronger association to obesity than single SNPs. Synthetic association analyses revealed that MC4R coding variants had almost no impact on the association signal. Carriers of the haplotype should be enriched for relevant mutations outside the MC4R coding region and could thus be used for re-sequencing approaches. Our data also underscore the problems underlying the identification of relevant mutations depicted by GWAS derived SNPs.     

SNP rs10248565 in HDAC9 as a novel genomic aberration biomarker of lung adenocarcinoma in non-smoking women

Background

Numerous efforts have been made to elucidate the etiology and improve the treatment of lung cancer, but the overall five-year survival rate is still only 15%. Although cigarette smoking is the primary risk factor for lung cancer, only 7% of female lung cancer patients in Taiwan have a history of smoking. Since cancer results from progressive accumulation of genetic aberrations, genomic rearrangements may be early events in carcinogenesis.

Results

In order to identify biomarkers of early-stage adenocarcinoma, the genome-wide DNA aberrations of 60 pairs of lung adenocarcinoma and adjacent normal lung tissue in non-smoking women were examined using Affymetrix Genome-Wide Human SNP 6.0 arrays. Common copy number variation (CNV) regions were identified by ≥30% of patients with copy number beyond 2 ± 0.5 of copy numbers for each single nucleotide polymorphism (SNP) and at least 100 continuous SNP variant loci. SNPs associated with lung adenocarcinoma were identified by McNemar’s test. Loss of heterozygosity (LOH) SNPs were identified in ≥18% of patients with LOH in the locus. Aberration of SNP rs10248565 at HDAC9 in chromosome 7p21.1 was identified from concurrent analyses of CNVs, SNPs, and LOH.

Conclusion

The results elucidate the genetic etiology of lung adenocarcinoma by demonstrating that SNP rs10248565 may be a potential biomarker of cancer susceptibility.     

Association of Genetic Markers with CSF Oligoclonal Bands in Multiple Sclerosis Patients

Objective

to explore the association between genetic markers and Oligoclonal Bands (OCB) in the Cerebro Spinal Fluid (CSF) of Italian Multiple Sclerosis patients.

Methods

We genotyped 1115 Italian patients for HLA-DRB1*15 and HLA-A*02. In a subset of 925 patients we tested association with 52 non-HLA SNPs associated with MS susceptibility and we calculated a weighted Genetic Risk Score. Finally, we performed a Genome Wide Association Study (GWAS) with OCB status on a subset of 562 patients. The best associated SNPs of the Italian GWAS were replicated in silico in Scandinavian and Belgian populations, and meta-analyzed.

Results

HLA-DRB1*15 is associated with OCB+: p = 0.03, Odds Ratio (OR) = 1.6, 95% Confidence Limits (CL) = 1.1–2.4. None of the 52 non-HLA MS susceptibility loci was associated with OCB, except one SNP (rs2546890) near IL12B gene (OR: 1.45; 1.09–1.92). The weighted Genetic Risk Score mean was significantly (p = 0.0008) higher in OCB+ (7.668) than in OCB− (7.412) patients. After meta-analysis on the three datasets (Italian, Scandinavian and Belgian) for the best associated signals resulted from the Italian GWAS, the strongest signal was a SNP (rs9320598) on chromosome 6q (p = 9.4×10⁻⁷) outside the HLA region (65 Mb).

Discussion

genetic factors predispose to the development of OCB.     

Are C-Reactive Protein Associated Genetic Variants Associated with Serum Levels and Retinal Markers of Microvascular Pathology in Asian Populations from Singapore?

Introduction

C-reactive protein (CRP) levels are associated with cardiovascular disease and systemic inflammation. We assessed whether CRP-associated loci were associated with serum CRP and retinal markers of microvascular disease, in Asian populations.

Methods

Genome-wide association analysis (GWAS) for serum CRP was performed in East-Asian Chinese (N = 2,434) and Malays (N = 2,542) and South-Asian Indians (N = 2,538) from Singapore. Leveraging on GWAS data, we assessed, in silico, association levels among the Singaporean datasets for 22 recently identified CRP-associated loci. At loci where directional inconsistencies were observed, quantification of inter-ethnic linkage disequilibrium (LD) difference was determined. Next, we assessed association for a variant at CRP and retinal vessel traits [central retinal artery equivalent (CRAE) and central retinal vein equivalent (CRVE)] in a total of 24,132 subjects of East-Asian, South-Asian and European ancestry.

Results

Serum CRP was associated with SNPs in/near APOE, CRP, HNF1A and LEPR (p-values ≤4.7×10⁻⁸) after meta-analysis of Singaporean populations. Using a candidate-SNP approach, we further replicated SNPs at 4 additional loci that had been recently identified to be associated with serum CRP (IL6R, GCKR, IL6 and IL1F10) (p-values ≤0.009), in the Singaporean datasets. SNPs from these 8 loci explained 4.05% of variance in serum CRP. Two SNPs (rs2847281 and rs6901250) were detected to be significant (p-value ≤0.036) but with opposite effect directions in the Singaporean populations as compared to original European studies. At these loci we did not detect significant inter-population LD differences. We further did not observe a significant association between CRP variant and CRVE or CRAE levels after meta-analysis of all Singaporean and European datasets (p-value >0.058).

Conclusions

Common variants associated with serum CRP, first detected in primarily European studies, are also associated with CRP levels in East-Asian and South-Asian populations. We did not find a causal link between CRP and retinal measures of microvascular disease.     

A Single Nucleotide Polymorphism within DUSP9 Is Associated with Susceptibility to Type 2 Diabetes in a Japanese Population

Aims

The DUSP9 locus on chromosome X was identified as a susceptibility locus for type 2 diabetes in a meta-analysis of European genome-wide association studies (GWAS), and GWAS in South Asian populations identified 6 additional single nucleotide polymorphism (SNP) loci for type 2 diabetes. However, the association of these loci with type 2 diabetes have not been examined in the Japanese. We performed a replication study to investigate the association of these 7 susceptibility loci with type 2 diabetes in the Japanese population.

Methods

We genotyped 11,319 Japanese participants (8,318 with type 2 diabetes and 3,001 controls) for each of the 7 SNPs–rs5945326 near DUSP9, rs3923113 near GRB14, rs16861329 in ST6GAL1, rs1802295 in VPS26A, rs7178572 in HMG20A, rs2028299 near AP3S2, and rs4812829 in HNF4A–and examined the association of each of these 7 SNPs with type 2 diabetes by using logistic regression analysis.

Results

All SNPs had the same direction of effect (odds ratio [OR]>1.0) as in the original reports. One SNP, rs5945326 near DUSP9, was significantly associated with type 2 diabetes at a genome-wide significance level (p = 2.21×10⁻⁸; OR 1.39, 95% confidence interval [CI]: 1.24−1.56). The 6 SNPs derived from South Asian GWAS were not significantly associated with type 2 diabetes in the Japanese population by themselves (p≥0.007). However, a genetic risk score constructed from 6 South Asian GWAS derived SNPs was significantly associated with Japanese type 2 diabetes (p = 8.69×10⁻⁴, OR  = 1.06. 95% CI; 1.03−1.10).

Conclusions/interpretation

These results indicate that the DUSP9 locus is a common susceptibility locus for type 2 diabetes across different ethnicities, and 6 loci identified in South Asian GWAS also have significant effect on susceptibility to Japanese type 2 diabetes.     

Genome-wide association analyses for meat quality traits in Chinese Erhualian pigs and a Western Duroc?×?(Landrace × Yorkshire) commercial population

Background

Understanding the genetic mechanisms that underlie meat quality traits is essential to improve pork quality. To date, most quantitative trait loci (QTL) analyses have been performed on F₂ crosses between outbred pig strains and have led to the identification of numerous QTL. However, because linkage disequilibrium is high in such crosses, QTL mapping precision is unsatisfactory and only a few QTL have been found to segregate within outbred strains, which limits their use to improve animal performance. To detect QTL in outbred pig populations of Chinese and Western origins, we performed genome-wide association studies (GWAS) for meat quality traits in Chinese purebred Erhualian pigs and a Western Duroc × (Landrace × Yorkshire) (DLY) commercial population.

Methods

Three hundred and thirty six Chinese Erhualian and 610 DLY pigs were genotyped using the Illumina PorcineSNP60K Beadchip and evaluated for 20 meat quality traits. After quality control, 35 985 and 56 216 single nucleotide polymorphisms (SNPs) were available for the Chinese Erhualian and DLY datasets, respectively, and were used to perform two separate GWAS. We also performed a meta-analysis that combined P-values and effects of 29 516 SNPs that were common to Erhualian, DLY, F₂ and Sutai pig populations.

Results

We detected 28 and nine suggestive SNPs that surpassed the significance level for meat quality in Erhualian and DLY pigs, respectively. Among these SNPs, ss131261254 on pig chromosome 4 (SSC4) was the most significant (P = 7.97E-09) and was associated with drip loss in Erhualian pigs. Our results suggested that at least two QTL on SSC12 and on SSC15 may have pleiotropic effects on several related traits. All the QTL that were detected by GWAS were population-specific, including 12 novel regions. However, the meta-analysis revealed seven novel QTL for meat characteristics, which suggests the existence of common underlying variants that may differ in frequency across populations. These QTL regions contain several relevant candidate genes.

Conclusions

These findings provide valuable insights into the molecular basis of convergent evolution of meat quality traits in Chinese and Western breeds that show divergent phenotypes. They may contribute to genetic improvement of purebreds for crossbred performance.

Electronic supplementary material

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

A Common Variant in the SIAH2 Locus Is Associated with Estrogen Receptor-Positive Breast Cancer in the Chinese Han Population

Background

Genome-wide association studies (GWAS) have identified many loci associated with breast cancer risk. These studies have primarily been conducted in populations of European descent.

Objective

To determine whether previously reported susceptibility loci in other ethnic groups are also risk factors for breast cancer in a Chinese population.

Method

We genotyped 21 previously reported single nucleotide polymorphisms (SNPs) within a female Chinese cohort of 1203 breast cancer cases and 2525 healthy controls using the Sequenom iPlex platform. Fourteen SNPs passed the quality control test. These SNPs were subjected to statistical analysis for the entire cohort and were further analyzed for estrogen receptor (ER) status. The associations of the SNPs with disease susceptibility were assessed using logistic regression, adjusting for age. The Bonferroni correction was used to conservatively account for multiple testing, and the threshold for statistical significance was P<3.57×10⁻³ (0.05/14).

Result

Although none of the SNPs showed an overall association with breast cancer, an analysis of the ER status of the breast cancer patients revealed that the SIAH2 locus (rs6788895; P = 5.73×10⁻⁴, odds ratio [OR] = 0.81) is associated with ER-positive breast cancer.

Conclusion

A common variant in the SIAH2 locus is associated with ER-positive breast cancer in the Chinese Han population. The replication of published GWAS results in other ethnic groups provides important information regarding the genetic etiology of breast cancer.     

Genome-Wide Association Study Identifies Novel Pharmacogenomic Loci For Therapeutic Response to Montelukast in Asthma

Background

Genome-wide association study (GWAS) is a powerful tool to identify novel pharmacogenetic single nucleotide polymorphisms (SNPs). Leukotriene receptor antagonists (LTRAs) are a major class of asthma medications, and genetic factors contribute to variable responses to these drugs. We used GWAS to identify novel SNPs associated with the response to the LTRA, montelukast, in asthmatics.

Methods

Using genome-wide genotype and phenotypic data available from American Lung Association - Asthma Clinical Research Center (ALA-ACRC) cohorts, we evaluated 8-week change in FEV₁ related to montelukast administration in a discovery population of 133 asthmatics. The top 200 SNPs from the discovery GWAS were then tested in 184 additional samples from two independent cohorts.

Results

Twenty-eight SNP associations from the discovery GWAS were replicated. Of these, rs6475448 achieved genome-wide significance (combined P = 1.97 x 10^-09), and subjects from all four studies who were homozygous for rs6475448 showed increased ΔFEV₁ from baseline in response to montelukast.

Conclusions

Through GWAS, we identified a novel pharmacogenomic locus related to improved montelukast response in asthmatics.