Failure to replicate a genetic association may provide important clues about genetic architecture

Analysis of gene-gene interactions using gene-trait similarity regression

OBJECTIVE

Gene-gene interactions (G×G) are important to study because of their extensiveness in biological systems and their potential in explaining missing heritability of complex traits. In this work, we propose a new similarity-based test to assess G×G at the gene level, which permits the study of epistasis at biologically functional units with amplified interaction signals.

METHODS

Under the framework of gene-trait similarity regression (SimReg), we propose a gene-based test for detecting G×G. SimReg uses a regression model to correlate trait similarity with genotypic similarity across a gene. Unlike existing gene-level methods based on leading principal components (PCs), SimReg summarizes all information on genotypic variation within a gene and can be used to assess the joint/interactive effects of two genes as well as the effect of one gene conditional on another.

RESULTS

Using simulations and a real data application to the Warfarin study, we show that the SimReg G×G tests have satisfactory power and robustness under different genetic architecture when compared to existing gene-based interaction tests such as PC analysis or partial least squares. A genome-wide association study with approx. 20,000 genes may be completed on a parallel computing system in 2 weeks.

Why epistasis is important for tackling complex human disease genetics

EDITORIAL SUMMARY

Epistasis has been dismissed by some as having little role in the genetic architecture of complex human disease. The authors argue that this view is the result of a misconception and explain why exploring epistasis is likely to be crucial to understanding and predicting complex disease.

Epiregulin (EREG) and human V-ATPase (TCIRG1): genetic variation, ethnicity and pulmonary tuberculosis susceptibility in Guinea-Bissau and The Gambia

We analyzed two West African samples (Guinea-Bissau: n=289 cases and 322 controls; The Gambia: n=240 cases and 248 controls) to evaluate single-nucleotide polymorphisms (SNPs) in Epiregulin (EREG) and V-ATPase (T-cell immune regulator 1 (TCIRG1)) using single and multilocus analyses to determine whether previously described associations with pulmonary tuberculosis (PTB) in Vietnamese and Italians would replicate in African populations. We did not detect any significant single locus or haplotype associations in either sample. We also performed exploratory pairwise interaction analyses using Visualization of Statistical Epistasis Networks (ViSEN), a novel method to detect only interactions among multiple variables, to elucidate possible interaction effects between SNPs and demographic factors. Although we found no strong evidence of marginal effects, there were several significant pairwise interactions that were identified in either the Guinea-Bissau or the Gambian samples, two of which replicated across populations. Our results indicate that the effects of EREG and TCIRG1 variants on PTB susceptibility, to the extent that they exist, are dependent on gene-gene interactions in West African populations as detected with ViSEN. In addition, epistatic effects are likely to be influenced by inter- and intra-population differences in genetic or environmental context and/or the mycobacterial lineages causing disease.

The association between endothelial nitric oxide synthase gene G894T polymorphism and hypertension in Han Chinese: a case-control study and an updated meta-analysis

Abstract Background: The G894T (rs1799983) polymorphism in endothelial nitric oxide synthase (eNOS/NOS3) gene has been implicated in susceptibility to essential hypertension (EH) in some studies, but no clear consensus has been reached in the Chinese population.

AIMS

This study aimed to investigate the association of the G894T polymorphism and EH in Han Chinese.

SUBJECTS AND METHODS

First, a case-control study was performed involving 1525 subjects in northern Han Chinese to study the association between G894T variants and EH and then a meta-analysis was conducted of all available studies in Han Chinese. A total of 25 studies comprising 13,443 subjects were finally included in this meta-analysis.

RESULTS

The present case-control study failed to show significant association of G894T variant with EH in northern Han Chinese. The subsequent meta-analysis showed that this polymorphism might be associated with EH in Han Chinese (p < 0.001, OR = 1.32), especially in southern Han Chinese (p < 0.001, OR = 1.59), but not in northern Han Chinese (p = 0.12, OR = 1.16). The meta-regression analysis suggested that the geographic difference of subjects was related to heterogeneity (p = 0.029).

CONCLUSIONS

The relationship between the G894T polymorphism and hypertension in Han Chinese may be attributed to the difference in geographic background of subjects. It is necessary to carry out further research with a large sample size and focusing on gene-environment interactions.

Association of Parkinson disease risk loci with mild parkinsonian signs in older persons

IMPORTANCE

Parkinsonian motor signs are common in the aging population and are associated with adverse health outcomes. Compared with Parkinson disease (PD), potential genetic risk factors for mild parkinsonian signs have been largely unexplored.

OBJECTIVE

To determine whether PD susceptibility loci are associated with parkinsonism or substantia nigra pathology in a large community-based cohort of older persons.

DESIGN, SETTING, AND PARTICIPANTS

Eighteen candidate single-nucleotide polymorphisms from PD genome-wide association studies were evaluated in a joint clinicopathologic cohort. Participants included 1698 individuals and a nested autopsy collection of 821 brains from the Religious Orders Study and the Rush Memory and Aging Project, 2 prospective community-based studies.

MAIN OUTCOMES AND MEASURES

The primary outcomes were a quantitative measure of global parkinsonism or component measures of bradykinesia, rigidity, tremor, and gait impairment that were based on the motor Unified Parkinson's Disease Rating Scale. In secondary analyses, we examined associations with additional quantitative motor traits and postmortem indices, including substantia nigra Lewy bodies and neuronal loss.

RESULTS

Parkinson disease risk alleles in the MAPT (rs2942168; P = .0006) and CCDC62 (rs12817488; P = .004) loci were associated with global parkinsonism, and these associations remained after exclusion of patients with a PD diagnosis. Based on motor Unified Parkinson's Disease Rating Scale subscores, MAPT (P = .0002) and CCDC62 (P = .003) were predominantly associated with bradykinesia, and we further discovered associations between SREBF1 (rs11868035; P = .005) and gait impairment, SNCA (rs356220; P = .04) and rigidity, and GAK (rs1564282; P = .03) and tremor. In the autopsy cohort, only NMD3 (rs34016896; P = .03) was related to nigral neuronal loss, and no associations were detected with Lewy bodies.

CONCLUSIONS AND RELEVANCE

In addition to the established link to PD susceptibility, our results support a broader role for several loci in the development of parkinsonian motor signs and nigral pathology in older persons.

Etiology of Obesity Over the Life Span: Ecological and Genetic Highlights from Asian Countries

Obesity is a worldwide pandemic, and the prevalence rate has doubled since the 1980s. Asian countries are also experiencing the global epidemic of obesity with its related health consequences. The prevalence of overweight and obesity are increasing at an alarming rate across all age groups in Asia. These increases are mainly attributed to rapid economic growth, which leads to socio-economic, nutrition and lifestyle transitions, resulting in a positive energy balance. In addition, fat mass and obesity-associated gene variants, copy number variants in chromosomes and epigenetic modifications have shown positive associations with the risk of obesity among Asians. In this review highlights of prevalence and related ecological and genetic factors that could influence the rapid rise in obesity among Asian populations are discussed.

Identifying Potential Regions of Copy Number Variation for Bipolar Disorder

Bipolar disorder is a complex psychiatric disorder with high heritability, but its genetic determinants are still largely unknown. Copy number variation (CNV) is one of the sources to explain part of the heritability. However, it is a challenge to estimate discrete values of the copy numbers using continuous signals calling from a set of markers, and to simultaneously perform association testing between CNVs and phenotypic outcomes. The goal of the present study is to perform a series of data filtering and analysis procedures using a DNA pooling strategy to identify potential CNV regions that are related to bipolar disorder. A total of 200 normal controls and 200 clinically diagnosed bipolar patients were recruited in this study, and were randomly divided into eight control and eight case pools. Genome-wide genotyping was employed using Illumina Human Omni1-Quad array with approximately one million markers for CNV calling. We aimed at setting a series of criteria to filter out the signal noise of marker data and to reduce the chance of false-positive findings for CNV regions. We first defined CNV regions for each pool. Potential CNV regions were reported based on the different patterns of CNV status between cases and controls. Genes that were mapped into the potential CNV regions were examined with association testing, Gene Ontology enrichment analysis, and checked with existing literature for their associations with bipolar disorder. We reported several CNV regions that are related to bipolar disorder. Two CNV regions on chromosome 11 and 22 showed significant signal differences between cases and controls (p < 0.05). Another five CNV regions on chromosome 6, 9, and 19 were overlapped with results in previous CNV studies. Experimental validation of two CNV regions lent some support to our reported findings. Further experimental and replication studies could be designed for these selected regions.

Single-nucleotide polymorphisms in PtoCesA7 and their association with growth and wood properties in Populus tomentosa

Cellulose synthase (CesA) genes encode the enzymes that synthesize cellulose; therefore, CesAs play central roles in plant development and affect the yield and quality of wood, essential properties for industrial applications of plant biomass. To effectively manipulate wood biosynthesis in trees and improve wood quality, we thus require a better understanding of the natural variation in CesAs. Association studies have emerged as a powerful tool for identification of variation associated with quantitative traits. Here, we used a candidate gene-based association mapping approach to identify PtoCesA7 allelic variants that associate with growth and wood quality traits in Populus tomentosa. We isolated a full-length PtoCesA7 cDNA and observed high PtoCesA7 expression in xylem, consistent with the xylem-specific expression of CesA7. Nucleotide diversity and linkage disequilibrium (LD) in PtoCesA7, sampled from the P. tomentosa natural distribution, revealed that PtoCesA7 harbors high nucleotide diversity (π(T) = 0.0091) and low LD (r(2) ≥ 0.1, within 800 bp). By association analysis, we identified seven single-nucleotide polymorphisms (SNPs) (false discovery rate Q < 0.10) and 12 haplotypes (Q < 0.10) that associated with growth and wood properties, explaining 3.62-10.59 % of the phenotypic variance. We also validated 9 of the 10 significant marker-trait associations in at least one of three smaller subsets (climatic regions) or in a linkage-mapping population. Thus, our study identified functional PtoCesA7 allelic variants associated with growth and wood quality traits, giving new insights into genes affecting wood quality and quantity. From an applied perspective, the SNPs revealed in this study have potential applications in marker-assisted breeding.

Genetic influence on exercise-induced changes in physical function among mobility-limited older adults

To date, physical exercise is the only intervention consistently demonstrated to attenuate age-related declines in physical function. However, variability exists in seniors' responsiveness to training. One potential source of variability is the insertion (I allele) or deletion (D allele) of a 287 bp fragment in intron 16 of the angiotensin-converting enzyme (ACE) gene. This polymorphism is known to influence a variety of physiological adaptions to exercise. However, evidence is inconclusive regarding the influence of this polymorphism on older adults' functional responses to exercise. This study aimed to evaluate the association of ACE I/D genotypes with changes in physical function among Caucasian older adults (n = 283) following 12 mo of either structured, multimodal physical activity or health education. Measures of physical function included usual-paced gait speed and performance on the Short Physical Performance Battery (SPPB). After checking Hardy-Weinberg equilibrium, we used using linear regression to evaluate the genotype*treatment interaction for each outcome. Covariates included clinic site, body mass index, age, sex, baseline score, comorbidity, and use of angiotensin receptor blockers or ACE inhibitors. Genotype frequencies [II (19.4%), ID (42.4%), DD (38.2%)] were in Hardy-Weinberg equilibrium (P > 0.05). The genotype*treatment interaction was statistically significant for both gait speed (P = 0.002) and SPPB (P = 0.020). Exercise improved gait speed by 0.06 ± 0.01 m/sec and SPPB score by 0.72 ± 0.16 points among those with at least one D allele (ID/DD carriers), but function was not improved among II carriers. Thus, ACE I/D genotype appears to play a role in modulating functional responses to exercise training in seniors.

Replication of results of a genome-wide association study on lung cancer survival in a Korean population

Recently, a genome-wide association study (GWAS) identified single nucleotide polymorphisms (SNPs) that may influence the prognosis of early-stage non-small cell lung cancer (NSCLC) in Caucasians. We attempted to replicate the impact of genetic variants identified in the GWAS on lung cancer survival in a Korean population. A total of 363 patients with surgically resected NSCLCs were enrolled, and 12 SNPs were genotyped using the SEQUENOM MassARRAY iPLEX assay, TaqMan assay, or a polymerase chain reaction-restriction fragment length polymorphism analysis. The association between genotypes and overall survival (OS) was analyzed. Among the 12 SNPs, the rs6034368T>C was associated with OS. Patients with the rs6034368C allele showed a better OS than the patients with the rs6034368T allele (adjusted hazard ratio = 0.72, confidence interval = 0.56-0.93, P = 0.01). The rs12446308A>G had an effect on OS, but it was marginally significant (under a codominant model, adjusted hazard ratio = 1.85, confidence interval = 0.98-3.47, P = 0.06). We identified that the rs6034368T>C was associated with survival in early-stage NSCLC in a Korean population.

Epistasis and quantitative traits: using model organisms to study gene-gene interactions

The role of epistasis in the genetic architecture of quantitative traits is controversial, despite the biological plausibility that nonlinear molecular interactions underpin the genotype-phenotype map. This controversy arises because most genetic variation for quantitative traits is additive. However, additive variance is consistent with pervasive epistasis. In this Review, I discuss experimental designs to detect the contribution of epistasis to quantitative trait phenotypes in model organisms. These studies indicate that epistasis is common, and that additivity can be an emergent property of underlying genetic interaction networks. Epistasis causes hidden quantitative genetic variation in natural populations and could be responsible for the small additive effects, missing heritability and the lack of replication that are typically observed for human complex traits.

Genetic interactions found between calcium channel genes modulate amyloid load measured by positron emission tomography

Late-onset Alzheimer's disease (LOAD) is known to have a complex, oligogenic etiology, with considerable genetic heterogeneity. We investigated the influence of genetic interactions between genes in the Alzheimer's disease (AD) pathway on amyloid-beta (Aβ) deposition as measured by PiB or AV-45 ligand positron emission tomography (PET) to aid in understanding LOAD's genetic etiology. Subsets of the Alzheimer's Disease Neuroimaging Initiative (ADNI) cohorts were used for discovery and for two independent validation analyses. A significant interaction between RYR3 and CACNA1C was confirmed in all three of the independent ADNI datasets. Both genes encode calcium channels expressed in the brain. The results shown here support previous animal studies implicating interactions between these calcium channels in amyloidogenesis and suggest that the pathological cascade of this disease may be modified by interactions in the amyloid-calcium axis. Future work focusing on the mechanisms of such relationships may inform targets for clinical intervention.

Systems analysis of high-throughput data

Modern high-throughput assays yield detailed characterizations of the genomic, transcriptomic, and proteomic states of biological samples, enabling us to probe the molecular mechanisms that regulate hematopoiesis or give rise to hematological disorders. At the same time, the high dimensionality of the data and the complex nature of biological interaction networks present significant analytical challenges in identifying causal variations and modeling the underlying systems biology. In addition to identifying significantly disregulated genes and proteins, integrative analysis approaches that allow the investigation of these single genes within a functional context are required. This chapter presents a survey of current computational approaches for the statistical analysis of high-dimensional data and the development of systems-level models of cellular signaling and regulation. Specifically, we focus on multi-gene analysis methods and the integration of expression data with domain knowledge (such as biological pathways) and other gene-wise information (e.g.,  sequence or methylation data) to identify novel functional modules in the complex cellular interaction network.

Dopaminergic gene polymorphisms and cognitive function in a north Indian schizophrenia cohort

BACKGROUND

Associations of polymorphisms from dopaminergic neurotransmitter pathway genes have mostly been reported in Caucasian ancestry schizophrenia (SZ) samples. As studies investigating single SNPs with SZ have been inconsistent, more detailed analyses utilizing multiple SNPs with the diagnostic phenotype as well as cognitive function may be more informative. Therefore, these analyses were conducted in a north Indian sample.

METHODS

Indian SZ case-parent trios (n = 601 families); unscreened controls (n = 468) and an independent set of 118 trio families were analyzed. Representative SNPs in the Dopamine D3 receptor (DRD3), dopamine transporter (SLC6A3), vesicular monoamine transporter 2 (SLC18A2), catechol-o-methyltransferase (COMT) and dopamine beta-hydroxylase (DBH) were genotyped using SNaPshot/SNPlex assays (n = 59 SNPs). The Trail Making Test (TMT) was administered to a subset of the sample (n = 260 cases and n = 302 parents).

RESULTS

Eight SNPs were nominally associated with SZ in either case-control or family based analyses (p < 0.05, rs7631540 and rs2046496 in DRD3; rs363399 and rs10082463 in SLC18A2; rs4680, rs4646315 and rs9332377 in COMT). rs6271 at DBH was associated in both analyses. Haplotypes of DRD3 SNPs incorporating rs7631540-rs2134655-rs3773678-rs324030-rs6280-rs905568 showed suggestive associations in both case-parent and trio samples. At SLC18A2, rs10082463 was nominally associated with psychomotor performance and rs363285 with executive functions using the TMT but did not withstand multiple corrections.

CONCLUSIONS

Suggestive associations with dopaminergic genes were detected in this study, but convincing links between dopaminergic polymorphisms and SZ or cognitive function were not observed.

Genome-wide association study of autistic-like traits in a general population study of young adults

Lay abstract: It has been proposed that autistic-like traits in the general population lie on a continuum, with clinical Autism Spectrum Disorder (ASD), representing the extreme end of this distribution. The current study undertook a genome-wide association (GWA) scan of 965 young Western Australian adults to identify novel risk variants associated with autistic-like traits. No associations reached genome-wide significance; however, a review of nominally associated single nucleotide polymorphisms (SNPs) indicated two positional candidate loci that have been previously implicated in autistic-like trait etiology.

Scientific abstract: Research has proposed that autistic-like traits in the general population lie on a continuum, with clinical ASD representing the extreme end of this distribution. Inherent in this proposal is that biological mechanisms associated with clinical ASD may also underpin variation in autistic-like traits within the general population. A GWA study using 2,462,046 SNPs was undertaken for ASD in 965 individuals from the Western Australian Pregnancy Cohort (Raine) Study. No SNP associations reached genome-wide significance (p < 5.0 × 10⁻⁸). However, investigations into nominal observed SNP associations (p < 1.0 × 10⁻⁵) add support to two positional candidate genes previously implicated in ASD etiology, PRKCB1, and CBLN1. The rs198198 SNP (p = 9.587 × 10⁻⁶), is located within an intron of the protein kinase C, beta 1 (PRKCB1) gene on chromosome 16p11. The PRKCB1 gene has been previously reported in linkage and association studies for ASD, and its mRNA expression has been shown to be significantly down regulated in ASD cases compared with controls. The rs16946931 SNP (p = 1.78 × 10⁻⁶) is located in a region flanking the Cerebellin 1 (CBLN1) gene on chromosome 16q12.1. The CBLN1 gene is involved with synaptogenesis and is part of a gene family previously implicated in ASD. This GWA study is only the second to examine SNPs associated with autistic-like traits in the general population, and provides evidence to support roles for the PRKCB1 and CBLN1 genes in risk of clinical ASD.

Meta-analysis of genome-wide association data identifies novel susceptibility loci for obesity

Obesity is a major public health problem with strong genetic determination. Multiple genetic variants have been implicated for obesity by conducting genome-wide association (GWA) studies, primarily focused on body mass index (BMI). Fat body mass (FBM) is phenotypically more homogeneous than BMI and is more appropriate for obesity research; however, relatively few studies have been conducted on FBM. Aiming to identify variants associated with obesity, we carried out meta-analyses of seven GWA studies for BMI-related traits including FBM, and followed these analyses by de novo replication. The discovery cohorts consisted of 21 969 individuals from diverse ethnic populations and a total of over 4 million genotyped or imputed SNPs. The de novo replication cohorts consisted of 6663 subjects from two independent samples. To complement individual SNP-based association analyses, we also carried out gene-based GWA analyses in which all variations within a gene were considered jointly. Individual SNP-based association analyses identified a novel locus 1q21 [rs2230061, CTSS (Cathepsin S)] that was associated with FBM after the adjustment of lean body mass (LBM) (P = 3.57 × 10(-8)) at the genome-wide significance level. Gene-based association analyses identified a novel gene NLK (nemo-like kinase) in 17q11 that was significantly associated with FBM adjusted by LBM. In addition, we confirmed three previously reported obesity susceptibility loci: 16q12 [rs62033400, P = 1.97 × 10(-14), FTO (fat mass and obesity associated)], 18q22 [rs6567160, P = 8.09 × 10(-19), MC4R (melanocortin 4 receptor)] and 2p25 [rs939583, P = 1.07 × 10(-7), TMEM18 (transmembrane protein 18)]. We also found that rs6567160 may exert pleiotropic effects to both FBM and LBM. Our results provide additional insights into the molecular genetic basis of obesity and may provide future targets for effective prevention and therapeutic intervention.

Importance of different types of prior knowledge in selecting genome-wide findings for follow-up

Biological plausibility and other prior information could help select genome-wide association (GWA) findings for further follow-up, but there is no consensus on which types of knowledge should be considered or how to weight them. We used experts' opinions and empirical evidence to estimate the relative importance of 15 types of information at the single-nucleotide polymorphism (SNP) and gene levels. Opinions were elicited from 10 experts using a two-round Delphi survey. Empirical evidence was obtained by comparing the frequency of each type of characteristic in SNPs established as being associated with seven disease traits through GWA meta-analysis and independent replication, with the corresponding frequency in a randomly selected set of SNPs. SNP and gene characteristics were retrieved using a specially developed bioinformatics tool. Both the expert and the empirical evidence rated previous association in a meta-analysis or more than one study as conferring the highest relative probability of true association, whereas previous association in a single study ranked much lower. High relative probabilities were also observed for location in a functional protein domain, although location in a region evolutionarily conserved in vertebrates was ranked high by the data but not by the experts. Our empirical evidence did not support the importance attributed by the experts to whether the gene encodes a protein in a pathway or shows interactions relevant to the trait. Our findings provide insight into the selection and weighting of different types of knowledge in SNP or gene prioritization, and point to areas requiring further research.

Role of genetic heterogeneity and epistasis in bladder cancer susceptibility and outcome: a learning classifier system approach

BACKGROUND AND OBJECTIVE

Detecting complex patterns of association between genetic or environmental risk factors and disease risk has become an important target for epidemiological research. In particular, strategies that provide multifactor interactions or heterogeneous patterns of association can offer new insights into association studies for which traditional analytic tools have had limited success.

MATERIALS AND METHODS

To concurrently examine these phenomena, previous work has successfully considered the application of learning classifier systems (LCSs), a flexible class of evolutionary algorithms that distributes learned associations over a population of rules. Subsequent work dealt with the inherent problems of knowledge discovery and interpretation within these algorithms, allowing for the characterization of heterogeneous patterns of association. Whereas these previous advancements were evaluated using complex simulation studies, this study applied these collective works to a 'real-world' genetic epidemiology study of bladder cancer susceptibility.

RESULTS AND DISCUSSION

We replicated the identification of previously characterized factors that modify bladder cancer risk--namely, single nucleotide polymorphisms from a DNA repair gene, and smoking. Furthermore, we identified potentially heterogeneous groups of subjects characterized by distinct patterns of association. Cox proportional hazard models comparing clinical outcome variables between the cases of the two largest groups yielded a significant, meaningful difference in survival time in years (survivorship). A marginally significant difference in recurrence time was also noted. These results support the hypothesis that an LCS approach can offer greater insight into complex patterns of association.

CONCLUSIONS

This methodology appears to be well suited to the dissection of disease heterogeneity, a key component in the advancement of personalized medicine.

Genomic risk models improve prediction of longitudinal lipid levels in children and young adults

In clinical medicine, lipids are commonly measured biomarkers used to assess an individual's risk for cardiovascular disease, heart attack, and stroke. Accurately predicting longitudinal lipid levels based on genomic information can inform therapeutic practices and decrease cardiovascular risk by identifying high-risk patients prior to onset. Using genotyped and imputed genetic data from 523 unrelated Caucasian Americans from the Bogalusa Heart Study, surveyed on 4,026 occasions from 4 to 48 years of age, we generated various lipid genomic risk models based on previously reported markers. We observed a significant improvement in prediction over non-genetic risk models in high density lipoprotein cholesterol (increase in the squared correlation between observed and predicted values, ΔR (2) = 0.032), low density lipoprotein cholesterol (ΔR (2) = 0.053), total cholesterol (ΔR (2) = 0.043), and triglycerides (ΔR (2) = 0.031). Many of our approaches are based on an n-fold cross-validation procedure that are, by design, adaptable to a clinical environment.

Toll-like receptor-associated sequence variants and prostate cancer risk among men of African descent

BACKGROUND

Recent advances demonstrate a relationship between chronic/recurrent inflammation and prostate cancer (PCA). Among inflammatory regulators, toll-like receptors (TLRs) play a critical role in innate immune responses. However, it remains unclear whether variant TLR genes influence PCA risk among men of African descent. Therefore, we evaluated the impact of 32 TLR-associated single nucleotide polymorphisms (SNPs) on PCA risk among African-Americans and Jamaicans.

METHODS

SNP profiles of 814 subjects were evaluated using Illumina’s Veracode genotyping platform. Single and combined effects of SNPs in relation to PCA risk were assessed using age-adjusted logistic regression and entropy-based multifactor dimensionality reduction (MDR) models.

RESULTS

Seven sequence variants detected in TLR6, TOLLIP, IRAK4, IRF3 were marginally related to PCA. However, none of these effects remained significant after adjusting for multiple hypothesis testing. Nevertheless, MDR modeling revealed a complex interaction between IRAK4 rs4251545 and TLR2 rs1898830 as a significant predictor of PCA risk among U.S. men (permutation testing p-value = 0.001).

CONCLUSIONS

MDR identified an interaction between IRAK4 and TLR2 as the best two factor model for predicting PCA risk among men of African descent. However, these findings require further assessment and validation.

Association between putative functional variants in the PSMB9 gene and risk of melanoma--re-analysis of published melanoma genome-wide association studies

To mine possibly hidden causal single-nucleotide polymorphisms (SNPs) of melanoma, we investigated the association of SNPs in 76 M/G1 transition genes with melanoma risk using our published genome-wide association study (GWAS) data set with 1804 melanoma cases and 1026 cancer-free controls. We found multiple SNPs with P < 0.01 and performed validation studies for 18 putative functional SNPs in PSMB9 in two other GWAS data sets. Two SNPs (rs1351383 and rs2127675) were associated with melanoma risk in the GenoMEL data set (P = 0.013 and 0.004, respectively), but failed in validation using the Australian data set. Genotype-phenotype analysis revealed these two SNPs were significantly correlated with mRNA expression level of PSMB9. Further experiments revealed that SNP rs2071480, which is in high LD with rs1351383 and rs2127675, may have a weak effect on the promoter activity of PSMB9. Taken together, our data suggested that functional variants in PSMB9 may contribute to melanoma susceptibility.

The UDP-glucuronate decarboxylase gene family in Populus: structure, expression, and association genetics

In woody crop plants, the oligosaccharide components of the cell wall are essential for important traits such as bioenergy content, growth, and structural wood properties. UDP-glucuronate decarboxylase (UXS) is a key enzyme in the synthesis of UDP-xylose for the formation of xylans during cell wall biosynthesis. Here, we isolated a multigene family of seven members (PtUXS1-7) encoding UXS from Populus tomentosa, the first investigation of UXSs in a tree species. Analysis of gene structure and phylogeny showed that the PtUXS family could be divided into three groups (PtUXS1/4, PtUXS2/5, and PtUXS3/6/7), consistent with the tissue-specific expression patterns of each PtUXS. We further evaluated the functional consequences of nucleotide polymorphisms in PtUXS1. In total, 243 single-nucleotide polymorphisms (SNPs) were identified, with a high frequency of SNPs (1/18 bp) and nucleotide diversity (π_T = 0.01033, θ_w =0.01280). Linkage disequilibrium (LD) analysis showed that LD did not extend over the entire gene (r²<0.1, P<0.001, within 700 bp). SNP- and haplotype-based association analysis showed that nine SNPs (Q <0.10) and 12 haplotypes (P<0.05) were significantly associated with growth and wood property traits in the association population (426 individuals), with 2.70% to 12.37% of the phenotypic variation explained. Four significant single-marker associations (Q <0.10) were validated in a linkage mapping population of 1200 individuals. Also, RNA transcript accumulation varies among genotypic classes of SNP10 was further confirmed in the association population. This is the first comprehensive study of the UXS gene family in woody plants, and lays the foundation for genetic improvements of wood properties and growth in trees using genetic engineering or marker-assisted breeding.

A survey of putative anxiety-associated genes in panic disorder patients with and without bladder symptoms

BACKGROUND

We have previously described a subtype of panic disorder (PD) that we termed 'bladder syndrome', characterized by urological and bladder symptoms (and possibly interstitial cystitis) in the patients and/or their family members and confirmed the validity of this subset in family linkage and association analysis. In this study, we determine (a) whether 20 single-nucleotide polymorphisms (SNPs) reported in the literature can be replicated in a new PD dataset and (b) whether dividing the sample into those with and without the 'bladder syndrome' can help to resolve the genetic heterogeneity within this new sample.

METHODS

We selected 20 putative associated SNPs from the literature, taken from studies published since 2004. We tested these SNPs for association in a sample of 351 PD patients and 552 controls, and then divided them into subgroups of 92 patients from bladder families and 259 from nonbladder families.

RESULTS

(a) When analyzed in all PD patients, none of the 20 SNPs appeared to be replicated (except for SLC6A4 from our previous study, but in a sample that overlaps substantially with that in our previous report). (b) However, some intriguing findings emerged when we separated bladder from nonbladder families: SLC6A4, reported by us previously, yielded stronger evidence than before (P=0.0018) when examined only in nonbladder families, and in contrast, is not statistically significant in bladder families. Two other markers yielded nominally significant results in bladder families - rs5751876 in ADORA2A (P=0.046) and rs12579350 in TMEM16B (P=0.035) - but were not significant in nonbladder families. (c) Two markers had noticeably lower P-values when we differentiated the women and analyzed them separately - rs12579350 in TMEM16B (P-value decreased from 0.035, as above, to 0.00055) and a different SNP in ADORA2A, rs4822492 (P-value decreases from 0.07 to 0.028).

SIGNIFICANCE

Our results indicate that most of the 20 reported associations do not hold up when PD is analyzed as one group. However, our findings provide further evidence that PD with bladder symptoms may be genetically different from PD without bladder. We suggest that it is worth pursuing SLC6A4 in nonbladder PD, and ADORA2A and TMEM16B in bladder PD. Also, the possibility of a male-female difference in PD is worth pursuing. We also briefly discuss issues of replication and multiple tests.

Association between Interferon-Inducible Protein 6 (IFI6) Polymorphisms and Hepatitis B Virus Clearance

CD8+ T cells are key factors mediating hepatitis B virus (HBV) clearance. However, these cells are killed through HBV-induced apoptosis during the antigen-presenting period in HBV-induced chronic liver disease (CLD) patients. Interferon-inducible protein 6 (IFI6) delays type I interferon-induced apoptosis in cells. We hypothesized that single nucleotide polymorphisms (SNPs) in the IFI6 could affect the chronicity of CLD. The present study included a discovery stage, in which 195 CLD patients, including chronic hepatitis B (HEP) and cirrhosis patients and 107 spontaneous recovery (SR) controls, were analyzed. The genotype distributions of rs2808426 (C > T) and rs10902662 (C > T) were significantly different between the SR and HEP groups (odds ratio [OR], 6.60; 95% confidence interval [CI], 1.64 to 26.52, p = 0.008 for both SNPs) and between the SR and CLD groups (OR, 4.38; 95% CI, 1.25 to 15.26; p = 0.021 and OR, 4.12; 95% CI, 1.18 to 14.44; p = 0.027, respectively). The distribution of diplotypes that contained these SNPs was significantly different between the SR and HEP groups (OR, 6.58; 95% CI, 1.63 to 25.59; p = 0.008 and OR, 0.15; 95% CI, 0.04 to 0.61; p = 0.008, respectively) and between the SR and CLD groups (OR, 4.38; 95% CI, 1.25 to 15.26; p = 0.021 and OR, 4.12; 95% CI, 1.18 to 14.44; p = 0.027, respectively). We were unable to replicate the association shown by secondary enrolled samples. A large-scale validation study should be performed to confirm the association between IFI6 and HBV clearance.

HLA distribution in COPD patients

BACKGROUND

Auto-immunity may contribute to the pathogenesis of chronic obstructive pulmonary disease (COPD), particularly to the presence of emphysema. Auto-immune diseases are characterized by an abnormal distribution of HLA class II alleles (DR and DQ). The distribution of DRB1 and DQB1 alleles has not been investigated in COPD.

METHODS

To this end, HLA medium-low resolution typing was performed following standardized protocols in 320 clinically stable COPD patients included in the PAC-COPD study. Results were compared with controls of the same geographical and ethnic origin, and potential relationships with the severity of airflow limitation and lung diffusing capacity impairment were explored in patients with COPD.

RESULTS

The distribution of DRB1 and DQB1 alleles in COPD was similar to that of controls except for a significantly higher prevalence of DRB1*14 in patients with severe airflow limitation and low diffusing capacity.

CONCLUSIONS

By and large, HLA distribution was similar in COPD patients and controls, but the HLA class II allele DRB1*14 may contribute to the pathogenesis of severe COPD with emphysema.

Genetic variants and their interactions in disease risk prediction - machine learning and network perspectives

A central challenge in systems biology and medical genetics is to understand how interactions among genetic loci contribute to complex phenotypic traits and human diseases. While most studies have so far relied on statistical modeling and association testing procedures, machine learning and predictive modeling approaches are increasingly being applied to mining genotype-phenotype relationships, also among those associations that do not necessarily meet statistical significance at the level of individual variants, yet still contributing to the combined predictive power at the level of variant panels. Network-based analysis of genetic variants and their interaction partners is another emerging trend by which to explore how sub-network level features contribute to complex disease processes and related phenotypes. In this review, we describe the basic concepts and algorithms behind machine learning-based genetic feature selection approaches, their potential benefits and limitations in genome-wide setting, and how physical or genetic interaction networks could be used as a priori information for providing improved predictive power and mechanistic insights into the disease networks. These developments are geared toward explaining a part of the missing heritability, and when combined with individual genomic profiling, such systems medicine approaches may also provide a principled means for tailoring personalized treatment strategies in the future.

An evolutionary perspective on epistasis and the missing heritability

The relative importance between additive and non-additive genetic variance has been widely argued in quantitative genetics. By approaching this question from an evolutionary perspective we show that, while additive variance can be maintained under selection at a low level for some patterns of epistasis, the majority of the genetic variance that will persist is actually non-additive. We propose that one reason that the problem of the “missing heritability” arises is because the additive genetic variation that is estimated to be contributing to the variance of a trait will most likely be an artefact of the non-additive variance that can be maintained over evolutionary time. In addition, it can be shown that even a small reduction in linkage disequilibrium between causal variants and observed SNPs rapidly erodes estimates of epistatic variance, leading to an inflation in the perceived importance of additive effects. We demonstrate that the perception of independent additive effects comprising the majority of the genetic architecture of complex traits is biased upwards and that the search for causal variants in complex traits under selection is potentially underpowered by parameterising for additive effects alone. Given dense SNP panels the detection of causal variants through genome-wide association studies may be improved by searching for epistatic effects explicitly.

In this study we have shown that two independent problems may have a common cause. Why do traits under selection exhibit additive genetic variance, and why is the proportion of the heritability explained by additive effects much smaller than the total heritability estimated to exist? Our results indicate that epistatic interactions can allow deleterious mutations to persist under selection and that these interactions can abate the depletion of additive genetic variation. Furthermore, a much larger element of non-additive genetic variance is maintained, which supports the notion that the heritability estimated from family studies could be a mixture of both additive and non-additive components. We show that searching directly for epistatic effects greatly improves the discovery of variants under selection, despite the multiple testing penalty being much larger. Finally, we demonstrate that common practices in genome-wide association studies could lead to both an ascertainment bias in detecting additive effects and a confirmation bias in perceiving that most of the genetic variance is additive.

Tamoxifen-induced venothromboembolic events: exploring validation of putative genetic association

OBJECTIVE

A pilot study to examine accrual rates, efficiency of data capture approaches, study design and genotyping capacity for a future genetic validation study was undertaken.

DESIGN

The process pilot evaluated feasibility of applying a matched case-control design to validate association of two candidate estrogen receptor (ER) single nucleotide polymorphisms (SNPs) with incidence of venothromboembolic events (VTE) in breast cancer patients treated with tamoxifen where criteria included frequency matching by age, number of years diagnosed with breast cancer within 4-year intervals, and geographic residency.

SETTING

The study was conducted at Marshfield Clinic, in central Wisconsin.

PARTICIPANTS

Study-eligible cases with a breast cancer diagnosis between 1994 and 2006 who experienced a VTE within 5 years of last tamoxifen exposure were matched at a ratio of 1:4 to controls with a breast cancer diagnosed between 1994 and 2006 with no VTE history following tamoxifen exposure for ≥2 years.

METHODS

Feasibility of enrolling, phenotyping, and genotyping 20% of the total number of validated eligible cases and controls was tested in order to project enrollment rates and assess probability of enrolling the projected sample size for the prospective validation study and adequacy of planned data capture. Conditional logistic regression analysis was conducted for the matched case-control study design.

RESULTS

Enrollment accruals included 19 of 24 targeted cases (79%), and 74 of 96 (77%) targeted controls. Electronic data capture for most variables was nearly 100%. No unexpected statistically significant differences were observed between cases and controls. Capacity to conduct in-house screening for rs2234689 (ER1 PvuII), rs9340799 (ER1 XbaI), rs13146272 (CYP4V2), rs2227589 (SERPINC 1) and rs1613662 (GP6) was successfully established. Association of GP6 with VTE was further validated (P=0.0403; OR, 0.19).

CONCLUSION

Accrual rates to the larger prospective study will require a multi-center design to ensure enrollment of adequate numbers of cases and controls for achieving the projected sample size required to validate association of the ER SNPs. To prevent study failure due to poor accrual, the importance of conducting feasibility studies before launching large scale validation studies of genetic association and adverse drug events, is discussed.

Polymorphic simple sequence repeat (SSR) loci within cellulose synthase (PtoCesA) genes are associated with growth and wood properties in Populus tomentosa

Chinese white poplar (Populus tomentosa), an important commercial tree species for timber and pulp production in northern China, has been used to examine the individual genes and allelic diversity responsible for complex traits controlling growth and lignocellulosic biosynthesis. Taking advantage of the low degree of linkage disequilibrium (LD) within P. tomentosa association populations, we examined associations between 15 cellulose synthase (PtoCesA) genes and traits including growth and wood properties. Thirty-six novel simple sequence repeat (SSR) markers within PtoCesA genes were detected by re-sequencing and genotyped in an association population (460 individuals). Single-marker and haplotype-based LD approaches were used to identify significant marker-trait associations. Family-based linkage studies and real-time PCR testing were conducted to validate the functional significance of SSR variation. Fifteen single-marker associations from seven PtoCesA genes and nine haplotype-based associations within six genes were identified in the association population (false discovery rate Q < 0.05). Next, five SSR marker-trait associations (Q < 0.05) from four PtoCesA genes were successfully validated in a linkage mapping population (1200 individuals). The results imply a functional role for these genes in mediating wood properties, demonstrating the potential of combining single-marker and haplotype-based LD approaches to detect functional allelic variation underlying quantitative traits in a low-LD population.

Analysis of natural variation reveals neurogenetic networks for Drosophila olfactory behavior

Understanding the relationship between genetic variation and phenotypic variation for quantitative traits is necessary for predicting responses to natural and artificial selection and disease risk in human populations, but is challenging because of large sample sizes required to detect and validate loci with small effects. Here, we used the inbred, sequenced, wild-derived lines of the Drosophila melanogaster Genetic Reference Panel (DGRP) to perform three complementary genome-wide association (GWA) studies for natural variation in olfactory behavior. The first GWA focused on single nucleotide polymorphisms (SNPs) associated with mean differences in olfactory behavior in the DGRP, the second was an extreme quantitative trait locus GWA on an outbred advanced intercross population derived from extreme DGRP lines, and the third was for SNPs affecting the variance among DGRP lines. No individual SNP in any analysis was associated with variation in olfactory behavior by using a strict threshold accounting for multiple tests, and no SNP overlapped among the analyses. However, combining the top SNPs from all three analyses revealed a statistically enriched network of genes involved in cellular signaling and neural development. We used mutational and gene expression analyses to validate both candidate genes and network connectivity at a high rate. The lack of replication between the GWA analyses, small marginal SNP effects, and convergence on common cellular networks were likely attributable to epistasis. These results suggest that fully understanding the genotype-phenotype relationship requires a paradigm shift from a focus on single SNPs to pathway associations.

A reassessment of the possible effects of the serotonin transporter gene linked polymorphism 5-HTTLPR on premature ejaculation

A population-based sample of 1673 (valid phenotypic and genotypic data were available from 1412 individuals) Finnish male twins and siblings of twins aged 18-45 years provided questionnaire data regarding ejaculatory function as well as saliva samples for genotype analyses. Genotypic analyses were conducted controlling for between-subjects dependence. No significant association was found between the 5-HTTLPR polymorphism and a composite variable measuring premature ejaculation or between this polymorphism and a self-report measure of ejaculation latency time. Previously conducted studies have found contradicting results regarding the possible role of 5-HTTLPR in premature ejaculation. Methodological inconsistencies have been pointed out in these studies, which have all been conducted on rather small samples. While differences in terms of measurement of ejaculatory function could partly explain why positive findings from some earlier studies could not be replicated, the present study, given the large sample size and multifactorial measures used, indicated that the 5-HTTLPR polymorphism has a limited, if any, impact on ejaculatory function.

Allelic variation in PtGA20Ox associates with growth and wood properties in Populus spp

Populus tomentosa is an economically important tree crop that produces wood for lumber, pulp, paper, and biofuels. Wood quality traits are likely to be strongly affected by the plant hormone gibberellic acid (GA), which regulates growth. GA20Ox encodes one of the major regulatory enzymes of GA biosynthesis and may therefore play a large role in growth and wood quality. Here, linkage disequilibrium (LD) studies were used to identify significant associations between single nucleotide polymorphisms (SNPs) within PtGA20Ox and growth and wood-quality traits of P. tomentosa. We isolated a full-length GA20Ox cDNA from Populus tomentosa by reverse transcription (RT)-PCR; this 1401 bp cDNA clone had an open reading frame of 1158 bp and encoded a protein of 385 amino acids. PtGA20Ox transcripts were maximally expressed in the mature xylem of vascular tissues, suggesting that PtGA20Ox is highly expressed and specifically associated with secondary xylem formation. Resequencing the PtGA20Ox locus of 36 individuals identified 55 SNPs, and the frequency of SNPs was 1/31 bp. The 29 most common SNPs (frequency>0.1) were genotyped in an association population (426 individuals) that was also phenotyped for key growth and wood quality traits. LD did not extend over the entire gene (r(2)<0.1, within 500 bp), demonstrating that a candidate-gene-based LD approach may the best way to understand the molecular basis underlying quantitative variation in this species. SNP- and haplotype-based association analyses indicated that four SNPs (false discovery rate Q<0.05) and 14 haplotypes (P<0.05) were significantly associated with growth and wood properties. The phenotypic variance explained by each SNP ranged from 3.44% to 14.47%. The SNP markers identified in this study can be applied to breeding programs for the improvement of growth and wood-property traits by marker-assisted selection.

Identification and confirmation of an exonic splicing enhancer variation in exon 5 of the Alzheimer disease associated PICALM gene

Alzheimer's disease (AD) is a progressive neurodegenerative disorder characterized by memory and cognitive impairment and is the leading cause of dementia in the elderly. A number of genome wide association studies and subsequent replication studies have been published recently on late onset AD (LOAD). These studies identified several new susceptibility genes including phosphatidylinositol-binding clathrin assembly protein (PICALM) on chromosome 11. The aim of our study was to examine the entire coding sequence of PICALM to determine if the association could be explained by any previously undetected sequence variation. Therefore, we sequenced 48 cases and 48 controls homozygous for the risk allele in the signal SNP rs3851179. We did not find any new variants; however, rs592297, a known coding synonymous SNP that is part of an exonic splice enhancer region in exon 5, is in strong linkage disequilibrium with rs3851179 and should be examined for functional significance in Alzheimer pathophysiology.

Chromosome 1p13 genetic variants antagonize the risk of myocardial infarction associated with high ApoB serum levels

Background

Genetic variation at 1p13 modulates serum lipid levels and the risk of coronary heart disease through the regulation of serum lipid levels. Here we investigate if the interaction between genetic variants at 1p13 and serum lipid levels affects the risk of non-fatal myocardial infarction (MI) in the Stockholm Heart Epidemiology Program (SHEEP), a large population based case control study.

Methods

In the present study only non fatal MI cases (n = 1213, men/women: 852/361) and controls (n = 1516, men/women =1054/507) matched by age, sex and residential area, were included. Three SNPs 12740374 G/T, rs599839A/G and rs646776T/C mapping at 1p13 were analysed for association with serum lipid levels and the risk of MI by a weighted least square regression and logistic regression analyses, respectively. To analyse the effect of the interaction between genetic variants and serum lipid levels on the risk of MI, we applied the biological model of interaction that estimates the difference in risk, expressed as OR (95%CI), observed in the presence and in the absence of both exposures. One derived measure is the Synergy index (S) and 95%CI, where S > 1 indicates synergy and S < 1 antagonism between the two interaction terms.

Results

Rs12740374G/T and rs646776T/C were in strong linkage disequilibrium (LD) (r² = 0.99), therefore only rs599839A/G and rs646776 were included in the analysis. Consistently with published data, presence of the rare genotypes was associated with reduced total-, LDL-cholesterol and ApoB serum levels (all p < 0.05) as compared to the reference genotype, but was not associated with the risk of MI.

However, the increased risk of MI observed in individual exposed to high (≥75^th percentile) serum lipid levels was offset in subjects carrying the rare alleles G and C. In particular, the risk of MI associated with high ApoB serum levels OR (95%CI) 2.27 (1.86-2.77) was reduced to 1.76 (1.33-2.34) in the presence of the G allele at rs599839 with an S of 0.47 (0.20-0.90).

Conclusions

These results indicate that an antagonism between ApoB serum levels and genetic variants at 1p13 contributes to reduce the risk of non-fatal MI in the presence of high ApoB serum levels.

Association genetics in Corymbia citriodora subsp. variegata identifies single nucleotide polymorphisms affecting wood growth and cellulosic pulp yield

Wood is an important biological resource which contributes to nutrient and hydrology cycles through ecosystems, and provides structural support at the plant level. Thousands of genes are involved in wood development, yet their effects on phenotype are not well understood. We have exploited the low genomic linkage disequilibrium (LD) and abundant phenotypic variation of forest trees to explore allelic diversity underlying wood traits in an association study. Candidate gene allelic diversity was modelled against quantitative variation to identify SNPs influencing wood properties, growth and disease resistance across three populations of Corymbia citriodora subsp. variegata, a forest tree of eastern Australia. Nine single nucleotide polymorphism (SNP) associations from six genes were identified in a discovery population (833 individuals). Associations were subsequently tested in two smaller populations (130-160 individuals), 'validating' our findings in three cases for actin 7 (ACT7) and COP1 interacting protein 7 (CIP7). The results imply a functional role for these genes in mediating wood chemical composition and growth, respectively. A flip in the effect of ACT7 on pulp yield between populations suggests gene by environment interactions are at play. Existing evidence of gene function lends strength to the observed associations, and in the case of CIP7 supports a role in cortical photosynthesis.

A replication study confirms the association of dendritic cell immunoreceptor (DCIR) polymorphisms with ACPA - negative RA in a large Asian cohort

Objectives

Dendritic cell immunoreceptor (DCIR) has been implicated in development of autoimmune disorders in rodent and DCIR polymorphisms were associated with anti-citrullinated proteins antibodies (ACPA)-negative rheumatoid arthritis (RA) in Swedish Caucasians. This study was undertaken to further investigate whether DCIR polymorphisms are also risk factors for the development of RA in four Asian populations originated from China and Malaysia.

Methods

We genotyped two DCIR SNPs rs2377422 and rs10840759 in Han Chinese population (1,193 cases, 1,278 controls), to assess their association with RA. Subsequently, rs2377422 was further genotyped in three independent cohorts of Malaysian-Chinese subjects (MY_Chinese, 254 cases, 206 controls), Malay subjects (MY_ Malay, 515 cases, 986 controls), and Malaysian-Indian subjects (MY_Indian, 378 cases, 285 controls), to seek confirmation of association in various ethnic groups. Meta-analysis was preformed to evaluate the contribution of rs2377422 polymorphisms to the development of ACPA-negative RA in distinct ethnic groups. Finally, we carried out association analysis of rs2377422 polymorphisms with DCIR mRNA expression levels.

Results

DCIR rs2377422 was found to be significantly associated with ACPA -negative RA in Han Chinese (OR 1.92, 95% CI 1.27–2.90, P = 0.0020). Meta-analysis confirms DCIR rs2377422 as a risk factor for ACPA-negative RA across distinct ethnic groups (OR_overall = 1.17, 95% CI 1.06–1.30, P = 0.003). The SNP rs2377422 polymorphism showed significant association with DCIR mRNA expression level, i.e. RA-risk CC genotype exhibit a significant increase in the expression of DCIR (P = 0.0023, Kruskal–Wallis).

Conclusions

Our data provide evidence for association between DCIR rs2377422 and RA in non-Caucasian populations and confirm the influence of DCIR polymorphisms on RA susceptibility, especially on ACPA-negative RA.

Interaction among apoptosis-associated sequence variants and joint effects on aggressive prostate cancer

Background

Molecular and epidemiological evidence demonstrate that altered gene expression and single nucleotide polymorphisms in the apoptotic pathway are linked to many cancers. Yet, few studies emphasize the interaction of variant apoptotic genes and their joint modifying effects on prostate cancer (PCA) outcomes. An exhaustive assessment of all the possible two-, three- and four-way gene-gene interactions is computationally burdensome. This statistical conundrum stems from the prohibitive amount of data needed to account for multiple hypothesis testing.

Methods

To address this issue, we systematically prioritized and evaluated individual effects and complex interactions among 172 apoptotic SNPs in relation to PCA risk and aggressive disease (i.e., Gleason score ≥ 7 and tumor stages III/IV). Single and joint modifying effects on PCA outcomes among European-American men were analyzed using statistical epistasis networks coupled with multi-factor dimensionality reduction (SEN-guided MDR). The case-control study design included 1,175 incident PCA cases and 1,111 controls from the prostate, lung, colo-rectal, and ovarian (PLCO) cancer screening trial. Moreover, a subset analysis of PCA cases consisted of 688 aggressive and 488 non-aggressive PCA cases. SNP profiles were obtained using the NCI Cancer Genetic Markers of Susceptibility (CGEMS) data portal. Main effects were assessed using logistic regression (LR) models. Prior to modeling interactions, SEN was used to pre-process our genetic data. SEN used network science to reduce our analysis from > 36 million to < 13,000 SNP interactions. Interactions were visualized, evaluated, and validated using entropy-based MDR. All parametric and non-parametric models were adjusted for age, family history of PCA, and multiple hypothesis testing.

Results

Following LR modeling, eleven and thirteen sequence variants were associated with PCA risk and aggressive disease, respectively. However, none of these markers remained significant after we adjusted for multiple comparisons. Nevertheless, we detected a modest synergistic interaction between AKT3 rs2125230-PRKCQ rs571715 and disease aggressiveness using SEN-guided MDR (p = 0.011).

Conclusions

In summary, entropy-based SEN-guided MDR facilitated the logical prioritization and evaluation of apoptotic SNPs in relation to aggressive PCA. The suggestive interaction between AKT3-PRKCQ and aggressive PCA requires further validation using independent observational studies.

C-reactive protein (CRP) promoter polymorphisms influence circulating CRP levels in a genome-wide association study of African Americans

C-reactive protein (CRP) is an acute phase reactant protein produced primarily by the liver. Circulating CRP levels are influenced by genetic and non-genetic factors, including infection and obesity. Genome-wide association studies (GWAS) provide an unbiased approach towards identifying loci influencing CRP levels. None of the six GWAS for CRP levels has been conducted in an African ancestry population. The present study aims to: (i) identify genetic variants that influence serum CRP in African Americans (AA) using a genome-wide association approach and replicate these findings in West Africans (WA), (ii) assess transferability of major signals for CRP reported in European ancestry populations (EA) to AA and (iii) use the weak linkage disequilibrium (LD) structure characteristic of African ancestry populations to fine-map the previously reported CRP locus. The discovery cohort comprised 837 unrelated AA, with the replication of significant single-nucleotide polymorphisms (SNPs) assessed in 486 WA. The association analysis was conducted with 2 366 856 genotyped and imputed SNPs under an additive genetic model with adjustment for appropriate covariates. Genome-wide and replication significances were set at P < 5 × 10(-8) and P < 0.05, respectively. Ten SNPs in (CRP pseudogene-1) CRPP1 and CRP genes were associated with serum CRP (P = 2.4 × 10(-09) to 4.3 × 10(-11)). All but one of the top-scoring SNPs associated with CRP in AA were successfully replicated in WA. CRP signals previously identified in EA samples were transferable to AAs, and we were able to fine-map this signal, reducing the region of interest from the 25 kb of LD around the locus in the HapMap CEU sample to only 8 kb in our AA sample.

Variation in MSRA modifies risk of neonatal intestinal obstruction in cystic fibrosis

Meconium ileus (MI), a life-threatening intestinal obstruction due to meconium with abnormal protein content, occurs in approximately 15 percent of neonates with cystic fibrosis (CF). Analysis of twins with CF demonstrates that MI is a highly heritable trait, indicating that genetic modifiers are largely responsible for this complication. Here, we performed regional family-based association analysis of a locus that had previously been linked to MI and found that SNP haplotypes 5' to and within the MSRA gene were associated with MI (P = 1.99 × 10(-5) to 1.08 × 10(-6); Bonferroni P = 0.057 to 3.1 × 10(-3)). The haplotype with the lowest P value showed association with MI in an independent sample of 1,335 unrelated CF patients (OR = 0.72, 95% CI [0.53-0.98], P = 0.04). Intestinal obstruction at the time of weaning was decreased in CF mice with Msra null alleles compared to those with wild-type Msra resulting in significant improvement in survival (P = 1.2 × 10(-4)). Similar levels of goblet cell hyperplasia were observed in the ilea of the Cftr(-/-) and Cftr(-/-)Msra(-/-) mice. Modulation of MSRA, an antioxidant shown to preserve the activity of enzymes, may influence proteolysis in the developing intestine of the CF fetus, thereby altering the incidence of obstruction in the newborn period. Identification of MSRA as a modifier of MI provides new insight into the biologic mechanism of neonatal intestinal obstruction caused by loss of CFTR function.

Genome-wide association scan identifies a risk locus for preeclampsia on 2q14, near the inhibin, beta B gene

Elucidating the genetic architecture of preeclampsia is a major goal in obstetric medicine. We have performed a genome-wide association study (GWAS) for preeclampsia in unrelated Australian individuals of Caucasian ancestry using the Illumina OmniExpress-12 BeadChip to successfully genotype 648,175 SNPs in 538 preeclampsia cases and 540 normal pregnancy controls. Two SNP associations (rs7579169, p = 3.58×10⁻⁷, OR = 1.57; rs12711941, p = 4.26×10⁻⁷, OR = 1.56) satisfied our genome-wide significance threshold (modified Bonferroni p<5.11×10⁻⁷). These SNPs reside in an intergenic region less than 15 kb downstream from the 3′ terminus of the Inhibin, beta B (INHBB) gene on 2q14.2. They are in linkage disequilibrium (LD) with each other (r² = 0.92), but not (r²<0.80) with any other genotyped SNP ±250 kb. DNA re-sequencing in and around the INHBB structural gene identified an additional 25 variants. Of the 21 variants that we successfully genotyped back in the case-control cohort the most significant association observed was for a third intergenic SNP (rs7576192, p = 1.48×10⁻⁷, OR = 1.59) in strong LD with the two significant GWAS SNPs (r²>0.92). We attempted to provide evidence of a putative regulatory role for these SNPs using bioinformatic analyses and found that they all reside within regions of low sequence conservation and/or low complexity, suggesting functional importance is low. We also explored the mRNA expression in decidua of genes ±500 kb of INHBB and found a nominally significant correlation between a transcript encoded by the EPB41L5 gene, ∼250 kb centromeric to INHBB, and preeclampsia (p = 0.03). We were unable to replicate the associations shown by the significant GWAS SNPs in case-control cohorts from Norway and Finland, leading us to conclude that it is more likely that these SNPs are in LD with as yet unidentified causal variant(s).

Appetite regulation genes are associated with body mass index in black South African adolescents: a genetic association study

BACKGROUND

Obesity is a complex trait with both environmental and genetic contributors. Genome-wide association studies have identified several variants that are robustly associated with obesity and body mass index (BMI), many of which are found within genes involved in appetite regulation. Currently, genetic association data for obesity are lacking in Africans-a single genome-wide association study and a few replication studies have been published in West Africa, but none have been performed in a South African population.

OBJECTIVE

To assess the association of candidate loci with BMI in black South Africans. The authors focused on single nucleotide polymorphisms (SNPs) in the FTO, LEP, LEPR, MC4R, NPY2R and POMC genes.

DESIGN

A genetic association study.

PARTICIPANTS

990 randomly selected individuals from the larger Birth to Twenty cohort (a longitudinal birth cohort study of health and development in Africans).

MEASURES

The authors genotyped 44 SNPs within the six candidate genes that included known BMI-associated SNPs and tagSNPs based on linkage disequilibrium in an African population for FTO, LEP and NPY2R. To assess population substructure, the authors included 18 ancestry informative markers. Weight, height, sex, sex-specific pubertal stage and exact age collected during adolescence (13 years) were used to identify loci that predispose to obesity early in life.

RESULTS

Sex, sex-specific pubertal stage and exact age together explain 14.3% of the variation in log(BMI) at age 13. After adjustment for these factors, four SNPs were individually significantly associated with BMI: FTO rs17817449 (p=0.022), LEP rs10954174 (p=0.0004), LEP rs6966536 (p=0.012) and MC4R rs17782313 (p=0.045). Together the four SNPs account for 2.1% of the variation in log(BMI). Each risk allele was associated with an estimated average increase of 2.5% in BMI.

CONCLUSIONS

The study highlighted SNPs in FTO and MC4R as potential genetic markers of obesity risk in South Africans. The association with two SNPs in the 3' untranslated region of the LEP gene is novel.

Performance analysis of novel methods for detecting epistasis

Background

Epistasis is recognized fundamentally important for understanding the mechanism of disease-causing genetic variation. Though many novel methods for detecting epistasis have been proposed, few studies focus on their comparison. Undertaking a comprehensive comparison study is an urgent task and a pathway of the methods to real applications.

Results

This paper aims at a comparison study of epistasis detection methods through applying related software packages on datasets. For this purpose, we categorize methods according to their search strategies, and select five representative methods (TEAM, BOOST, SNPRuler, AntEpiSeeker and epiMODE) originating from different underlying techniques for comparison. The methods are tested on simulated datasets with different size, various epistasis models, and with/without noise. The types of noise include missing data, genotyping error and phenocopy. Performance is evaluated by detection power (three forms are introduced), robustness, sensitivity and computational complexity.

Conclusions

None of selected methods is perfect in all scenarios and each has its own merits and limitations. In terms of detection power, AntEpiSeeker performs best on detecting epistasis displaying marginal effects (eME) and BOOST performs best on identifying epistasis displaying no marginal effects (eNME). In terms of robustness, AntEpiSeeker is robust to all types of noise on eME models, BOOST is robust to genotyping error and phenocopy on eNME models, and SNPRuler is robust to phenocopy on eME models and missing data on eNME models. In terms of sensitivity, AntEpiSeeker is the winner on eME models and both SNPRuler and BOOST perform well on eNME models. In terms of computational complexity, BOOST is the fastest among the methods. In terms of overall performance, AntEpiSeeker and BOOST are recommended as the efficient and effective methods. This comparison study may provide guidelines for applying the methods and further clues for epistasis detection.

A common missense single nucleotide polymorphism in the E-selectin gene is significantly associated with essential hypertension in the Han population but only weakly associated in the Uygur population

Experimental and clinical observations suggest that E-selectin could have an important role in essential hypertension (EH), but the relationship between common E-selectin variants and EH has not been extensively studied in the Chinese population. In this study, we explored the association between two common variants in the E-selectin gene (rs5361A/C and rs5355C/T) and EH in the Uygur, Kazakh and Han populations in the Xinjiang area. A case-control study was conducted to explore the association between these two single-nucleotide polymorphisms and EH in a large sample size, including 941 EH subjects (309 Uygur, 264 Kazakh and 368 Han individuals) and 924 control subjects (300 Uygur, 275 Kazakh and 349 Han individuals). Univariate analysis showed that the rs5361 A/C polymorphism was significantly associated with EH in the Uygur (P=0.002) and Han (P=3.6 × 10(-7)) populations. The CC genotype of this SNP was present only in patients with EH in all of the three nationalities studied. Han individuals who carry the CC genotype of rs5361 were more susceptible to EH, according to the dominant models (P=1.13 × 10(-4), odds ratio=3.812, 95% confidence interval: 1.685-7.792), but there was no association of genotype with EH for the other ethnicities. No significant difference in rs5355 C/T polymorphism rate was found between the EH and control groups. Our results indicate that the common variant rs5361 is strongly associated with EH in Han individuals and weakly associated in Uygur individuals. The CC genotype of rs5361 might be an independent risk factor for EH among Uygur, Kazakh and Han individuals in the Xinjiang area.

What should the genome-wide significance threshold be? Empirical replication of borderline genetic associations

BACKGROUND

Robust replication is a sine qua non for the rigorous documentation of proposed associations in the genome-wide association (GWA) setting. Currently, associations of common variants reaching P ≤ 5 × 10(-8) are considered replicated. However, there is some ambiguity about the most suitable threshold for claiming genome-wide significance.

METHODS

We defined as 'borderline' associations those with P > 5 × 10(-8) and P ≤ 1 × 10(-7). The eligible associations were retrieved using the 'Catalog of Published Genome-Wide Association Studies'. For each association we assessed whether it reached P ≤ 5 × 10(-8) with inclusion of additional data from subsequent GWA studies.

RESULTS

Thirty-four eligible genotype-phenotype associations were evaluated with data and clarifications contributed from diverse investigators. Replication data from subsequent GWA studies could be obtained for 26 of them. Of those, 19 associations (73%) reached P ≤ 5 × 10(-8) for the same or a related trait implicating either the exact same allele or one in very high linkage disequilibrium and 17 reached P < 10(-8). If the seven associations that did not reach P ≤ 5 × 10(-8) when additional data were considered are assumed to have been false-positives, the false-discovery rate for borderline associations is estimated to be 27% [95% confidence interval (CI) 12-48%]. For five associations, the current P-value is > 10(-6) [corresponding false-discovery rate 19% (95% CI 7-39%)].

CONCLUSION

A substantial proportion, but not all, of the associations with borderline genome-wide significance represent replicable, possibly genuine associations. Our empirical evaluation suggests a possible relaxation in the current GWS threshold.

Meta-analysis of the NAD(P)H: quinine oxidoreductase 1 gene 609 C>T polymorphism with esophageal cancer risk

Association between the NAD(P)H: quinine oxidoreductase 1 (NQO1) gene 609 C>T polymorphism and esophageal cancer (EC) has been widely evaluated; however, the results are often irreproducible. We thus aimed to comprehensively evaluate this association through a meta-analysis. Data were extracted from 10 study populations involving 1390 EC patients and 1812 controls, and were analyzed using STATA software. Random-effects model was applied irrespective of between-study heterogeneity, which was assessed by the inconsistency index (I(2)) statistic. Publication bias was weighted by the funnel plot and Egger's test. Genotype distributions of the NQO1 gene 609 C>T polymorphism met Hardy-Weinberg equilibrium in controls for all studies. Allelic comparison indicated that NQO1 609 T allele conferred an increased risk (odds ratio [OR]=1.23; 95% confidence interval [CI]: 1.02-1.49; p=0.035), accompanying significant heterogeneity (I(2)=63.4%, p=0.003) and no publication bias (p(Egger)=0.391). This association was potentially enhanced in homozygous comparison (OR=1.58; 95% CI: 1.03-2.41; p=0.035; I(2)= 55.4%, p(heterogeneity)=0.017 and p(Egger)=0.461). Under dominant and recessive models, similar associations were obtained with an increased, although marginally significant risk. Subgroup analysis by ethnicity supported the risk profiles of the NQO1 gene 609 T allele and 609 TT genotype with EC in Eastern Asians, not in Europeans. Meta-regression analysis indicated that association between the NQO1 gene 609 C>T polymorphism and EC risk was significantly decreased with aging in case-patients (R(2)=-0.57; p=0.042). We expand previous studies by showing that the NQO1 gene 609 C>T polymorphism might contribute to EC occurrence, especially in Eastern Asians.