Efficiency and power in genetic association studies

Frequency of common polymorphisms in ( ) gene in adults with high serum triglycerides from Colombian Caribbean Coast

Background

Caveolin 1 gene (CAV1) has been associated with insulin resistance, metabolic syndrome and hypertension in humans. Also, it has been related to high serum triglycerides in rodents, however there is little evidence of this relation in humans.

Aim

To describe frequencies of common variations in CAV1 in adults with high serum triglycerides.

Methods

A case-control study was carried out with adults from Colombian Caribbean Coast. A whole blood sample was employed to measure serum concentrations of triglycerides, glucose, total cholesterol and HDLc. Six common Single Nucleotide Polymorphism (SNP) in CAV1 were genotyped (rs926198, rs3779512, rs10270569, rs11773845, rs7804372 and rs1049337). Allelic and genotypic frequencies were determined by direct count and Hardy-Weinberg Equilibrium (HWE) was assessed. Case and control groups were compared with null-hypothesis tests.

Results

A total of 220 cases and 220 controls were included. For rs3779512 an excess in homozygotes frequency was found within case group (40.4% (GG), 41.3% (GT) and 18.1% (TT); F_is=0.13, p=0.03). Another homozygotes excess among case group was found in rs7804372 (59.5% (TT), 32.3% (TA) and 8.2% (AA); F_is= 0.12, p= 0.04). In rs1049337, cases also showed an excess in homozygotes frequency (52.7% (CC), 35.0% (CT) and 12.3% (TT); F_is= 0.16, p= 0.01). Finally, for rs1049337 there were differences in genotype distribution between case and control groups (p <0.05).

Conclusion

An increased frequency of homozygote genotypes was found in subjects with high serum triglycerides. These findings suggest that minor alleles for SNPs rs3779512, rs7804372 and rs1049337 might be associated to higher risk of hypertriglyceridemia.

Effects of genetic variants of ST8SIA2 and NCAM1 genes on seasonal mood changes and circadian preference in the general population

ST8SIA2 and NCAM1 are functionally related genes forming polysialic acid (PSA) - neural cell adhesion molecule (NCAM) complex in suprachiasmatic nucleus (SCN), the regulating site of circadian biological rhythm. In this study, the relationship of ST8SIA2 and NCAM1 with circadian and seasonal rhythms of human behavior was explored. Subjects were 261 healthy Korean adults who were free of any history of clinically significant psychiatric symptoms. The phenotypes were circadian preference and seasonal change of mood and behavior (seasonality) measured by the Composite Scale of Morningness and the Seasonal Pattern Assessment Questionnaire, respectively. Thirty-four single nucleotide polymorphisms (SNPs) across the ST8SIA2 region and 15 SNPs of NCAM1 were analyzed. A nominally significant association with seasonality and circadian preference was observed in 21 variants of both genes. After corrections for multiple testing, associations of 8 SNPs of ST8SIA2 and 2 SNPs of NCAM1 with seasonality remained significant. Some of these SNPs were also associated with psychiatric disorders in previous studies. This study demonstrated a meaningful and/or suggestive evidence of association between behavioral phenotypes reflecting human biological rhythm and two interplaying genes involved in the plasticity of SCN's neuronal network.

Increased expression of IL12B mRNA transcribed from the risk haplotype for Crohn's disease is a risk factor for disease relapse in Japanese patients

BACKGROUND

IL12B is a promising candidate for a susceptibility gene in Crohn's disease (CD). The aim of this study was to perform a candidate gene analysis of IL12B in Japanese CD patients, investigate whether the genotype is associated with disease phenotypes, and determine how the risk allele affects susceptibility to CD.

METHODS

Three hundred seventy-five patients with CD, 265 patients with ulcerative colitis, and 463 healthy controls were examined. Ten single-nucleotide polymorphisms (SNPs) around IL12B were genotyped. Case-control and subphenotype (including disease course) analyses were performed. The allelic expression ratio of IL12B messenger RNA (mRNA) was examined by a SNaPshot analysis in lipopolysaccharide-stimulated monocytes.

RESULTS

Four SNPs located upstream of the IL12B gene were significantly associated with CD. A conditional analysis revealed that these associations included two independent signals tagged by IL12B_1 and IL12B_3 (P = 9.42 × 10^-6 and  1.49 × 10^-4 respectively). IL12B_3 was also associated with earlier relapse in CD (P = 0.0144). The allelic expression ratios of IL12B mRNA transcribed from the risk haplotype to the protective haplotype tagged by IL12B_3 in lipopolysaccharide-stimulated monocytes from ten healthy controls heterozygous for IL12B_3 were significantly higher than that of the respective genomic DNA (P = 0.00923). No SNP was associated with ulcerative colitis.

CONCLUSIONS

We confirmed the association of SNPs located upstream of IL12B with CD in Japanese patients. The demonstrated allelic expression imbalance supports the idea that the IL12B risk haplotype confers susceptibility not only to CD onset but to also relapse through increased IL12B mRNA expression.

Evaluation of the interaction between genetic variants of GAD1 and miRNA in bipolar disorders

BACKGROUND

Glutamic acid dehydrogenase 1 (GAD1) serves as the rate-limiting enzyme for synthesizing GABA, and is reported to be associated with several psychiatric disorders. The present study examined the effects of GAD1 genetic variants on bipolar disorder (BD) and its subtypes. Moreover, we investigated functional interactions between genetic variants and miRNAs via algorithm prediction and experimental validation.

METHODS

A case-control study was conducted with 280 BD patients and 200 healthy controls. Eight tag SNPs in GAD1 were genotyped. For associated markers, we performed in silico prediction for their potential functions through SNP-miRNA interactions by establishing a scoring system to combine information from several miRNA predictive algorithms. We then tested allelic expression differences using Dual-Glo luciferase reporter assays for the selected SNP-miRNA pair. Lastly, we examined the associations of the GAD1 gene and BD in two additional independent datasets with a few thousand samples for replication.

RESULTS

Marker rs3749034 was associated with BD, in particular the BD-II subtype. According to our scoring system, several candidate miRNAs were predicted to interact with rs3749034, and hsa-miR-504 had the highest score. Findings from an in vitro experiment revealed a non-statistically significant trend for lower gene expression level with the A allele of rs3749034 compared with the G allele. The association between rs3749034 and BD was not replicated in either of the independent datasets. Instead, other rarer genetic variants in GAD1 showed suggestive signals (e.g. rs575441409, p-value = 3.8*10^-4, D' = 1 with rs3749034) with BD in the Taiwanese dataset.

LIMITATIONS

The present study considered common genetic variants only. In addition, we only used a 293T cell-line in conducting luciferase reporter assays, as no primary cell-lines from patient samples were available to differentiate the effects between BD subtypes.

CONCLUSIONS

Our results demonstrate a weak effect of the GAD1 gene on the risk of bipolar illness, and the associated marker might represent a proxy for real signals of rare variants.

Genetic variants determining survival and fertility in an adverse African environment: a population-based large-scale candidate gene association study

Human survival probability and fertility decline strongly with age. These life history traits have been shaped by evolution. However, research has failed to uncover a consistent genetic determination of variation in survival and fertility. As an explanation, such genetic determinants have been selected in adverse environments, in which humans have lived during most of their history, but are almost exclusively studied in populations in modern affluent environments. Here, we present a large-scale candidate gene association study in a rural African population living in an adverse environment. In 4387 individuals, we studied 4052 SNPs in 148 genes that have previously been identified as possible determinants of survival or fertility in animals or humans. We studied their associations with survival comparing newborns, middle-age adults, and old individuals. In women, we assessed their associations with reported and observed numbers of children. We found no statistically significant associations of these SNPs with survival between the three age groups nor with women's reported and observed fertility. Population stratification was unlikely to explain these results. Apart from a lack of power, we hypothesise that genetic heterogeneity of complex phenotypes and gene-environment interactions prevent the identification of genetic variants explaining variation in survival and fertility in humans.

Polymorphisms in CYP2C9 are associated with response to indomethacin among neonates with patent ductus arteriosus

BackgroundPatent ductus arteriosus (PDA) is a common complication seen in preterm infants. Indomethacin is routinely used to treat PDA. Evidence suggests that the response of indomethacin is highly heritable. This study investigated the association between single-nucleotide polymorphisms (SNPs) in CYP2C9 and the closure of PDA in response to indomethacin.MethodsSix SNPs in CYP2C9 were analyzed for association with indomethacin response. A case-control analysis was performed among neonates who responded to indomethacin (responders) and among those who required surgical ligation (non-responders). Independent transmission disequilibrium tests were performed among parent-child trios of responders and non-responders.ResultsThe G allele of rs2153628 was associated with increased odds of response to indomethacin in the case-control analysis (odds ratios (OR): 1.918, 95% confidence interval (CI): 1.056, 3.483). Among indomethacin responders, the G allele of rs2153628 and the T allele of rs1799853 were overtransmitted from the parents to their child (OR: 2.667, 95% CI: 1.374, 5.177 and OR: 2.375, 95% CI: 1.040, 5.425, respectively), consistent with the case-control analysis.ConclusionWe identified an association between two SNPs in CYP2C9, rs2153628 and rs1799853, and indomethacin response for the treatment of PDA. These findings suggest that response to indomethacin in the closure of PDA may be influenced by polymorphisms associated with altered indomethacin metabolism.

Association between taste receptor (TAS) genes and the perception of wine characteristics

Several studies have suggested a possible relationship between polymorphic variants of the taste receptors genes and the acceptance, liking and intake of food and beverages. In the last decade investigators have attempted to link the individual ability to taste 6-n-propylthiouracil (PROP) and the sensations, such as astringency and bitterness, elicited by wine or its components, but with contradictory results. We have used the genotype instead of the phenotype (responsiveness to PROP or other tastants), to test the possible relation between genetic variability and the perception of wine characteristic in 528 subjects from Italy and the Czech Republic. We observed several interesting associations, among which the association between several TAS2R38 gene single nucleotide polymorphisms (P = 0.002) and the TAS2R16-rs6466849 polymorphism with wine sourness P = 0.0003). These associations were consistent in both populations, even though the country of origin was an important factor in the two models, thus indicating therefore that genetics alongside cultural factors also play a significant role in the individual liking of wine.

Vitamin D receptor and calcium-sensing receptor polymorphisms and colorectal cancer survival in the Newfoundland population

Background:

Increased serum levels of vitamin D and calcium have been associated with lower risks of colorectal cancer (CRC) incidence and mortality. These inverse associations may be mediated by the vitamin D receptor (VDR) and the calcium-sensing receptor (CASR). We investigated genetic variants in VDR and CASR for their relevance to CRC prognosis.

Methods:

A population-based cohort of 531 CRC patients diagnosed from 1999 to 2003 in Newfoundland and Labrador, Canada, was followed for mortality and cancer recurrence until April 2010. Germline DNA samples were genotyped with the Illumina Omni-Quad 1 Million chip. Multivariate Cox models assessed 41 tag single-nucleotide polymorphisms and relative haplotypes on VDR and CASR in relation to all-cause mortality (overall survival, OS) and disease-free survival (DFS).

Results:

Gene-level associations were observed between VDR and the DFS of rectal cancer patients (P=0.037) as well as between CASR and the OS of colon cancer patients (P=0.014). Haplotype analysis within linkage blocks of CASR revealed the G-G-G-G-G-A-C haplotype (rs10222633-rs10934578-rs3804592-rs17250717-A986S-R990G-rs1802757) to be associated with a decreased OS of colon cancer (HR, 3.15; 95% CI, 1.66–5.96). Potential interactions were seen among prediagnostic dietary calcium intake with the CASR R990G (P_int=0.040) and the CASR G-T-G-G-G-G-C haplotype for rs10222633-rs10934578-rs3804592-rs17250717-A986S-R990G-rs1802757 (P_int=0.017), with decreased OS time associated with these variants limited to patients consuming dietary calcium below the median, although the stratified results were not statistically significant after correction for multiple testing.

Conclusions:

Polymorphic variations in VDR and CASR may be associated with survival after a diagnosis of CRC.

Variation in Maturity-Onset Diabetes of the Young Genes Influence Response to Interventions for Diabetes Prevention

Context

Variation in genes that cause maturity-onset diabetes of the young (MODY) has been associated with diabetes incidence and glycemic traits.

Objectives

This study aimed to determine whether genetic variation in MODY genes leads to differential responses to insulin-sensitizing interventions.

Design and Setting

This was a secondary analysis of a multicenter, randomized clinical trial, the Diabetes Prevention Program (DPP), involving 27 US academic institutions. We genotyped 22 missense and 221 common variants in the MODY-causing genes in the participants in the DPP.

Participants and Interventions

The study included 2806 genotyped DPP participants randomized to receive intensive lifestyle intervention (n = 935), metformin (n = 927), or placebo (n = 944).

Main Outcome Measures

Association of MODY genetic variants with diabetes incidence at a median of 3 years and measures of 1-year β-cell function, insulinogenic index, and oral disposition index. Analyses were stratified by treatment group for significant single-nucleotide polymorphism × treatment interaction (Pint < 0.05). Sequence kernel association tests examined the association between an aggregate of rare missense variants and insulinogenic traits.

Results

After 1 year, the minor allele of rs3212185 (HNF4A) was associated with improved β-cell function in the metformin and lifestyle groups but not the placebo group; the minor allele of rs6719578 (NEUROD1) was associated with an increase in insulin secretion in the metformin group but not in the placebo and lifestyle groups.

Conclusions

These results provide evidence that genetic variation among MODY genes may influence response to insulin-sensitizing interventions.

Pharmacogenetics of oral antidiabetes drugs: evidence for diverse signals at the IRS1 locus

To investigate the role of IRS1 locus on failure to oral antidiabetes drugs (OADs) we genotyped single-nucleotide polymorphisms (SNPs), rs2943641, rs7578326 (tagging all SNPs genome-wide associated with type 2 diabetes (T2D) and related traits at this locus) and rs1801278 (that is, the loss-of-function IRS1 G972R amino acid substitution) in 2662 patients with T2D. Although no association with OAD failure was observed for rs2943641 and rs7578326 SNPs (odds ratio (OR): 1.04, 95% confidence interval (CI): 0.93-1.16 and OR: 0.97, 95% CI: 0.87-1.09 respectively), a significant association was observed for rs1801278 (OR: 1.34, 95% CI: 1.08-1.66). When meta-analyzed with previous published data, an allelic OR of 1.41 (1.15-1.72; P=0.001) was obtained, so that homozygous R972R individuals have >80% higher risk of failing to OADs as compared with their G972G counterparts. In all, though further studies are needed for confirming this finding, our present data point to IRS1 rs1801278 as a potential biomarker for pursuing the goal of stratified medicine in the field of antihyperglycemic treatment in T2D.

The YWHAE gene confers risk to major depressive disorder in the male group of Chinese Han population

Schizophrenia and major depressive disorder are two major psychiatric illnesses that may share specific genetic risk factors to a certain extent. Increasing evidence suggests that the two disorders might be more closely related than previously considered. To investigate whether YWHAE gene plays a significant role in major depressive disorder in Han Chinese population, we recruited 1135 unrelated major depressive disorder patients (485 males, 650 females) and 989 unrelated controls (296 males, 693 females) of Chinese Han origin. Eleven common SNPs were genotyped using TaqMan® technology. In male-group, the allele and genotype frequencies of rs34041110 differed significantly between patients and control (P_allele=0.036486, OR[95%CI]: 1.249442(1.013988-1.539571); P_genotype=0.045301). Also in this group, allele and genotype frequencies of rs1532976 differed significantly (P_allele=0.013242, OR[95%CI]: 1.302007(1.056501-1.604563); genotype: P=0.039152). Haplotype-analyses showed that, in male-group, positive association with major depressive disorder was found for the A-A-C-G haplotype of rs3752826-rs2131431-rs1873827-rs12452627 (χ²=20.397, P=6.38E-06, OR[95%CI]: 7.442 [2.691-20.583]), its C-A-C-G haplotype (χ²=19.122, P=1.24E-05, OR and 95%CI: 0.402 [0.264-0.612]), its C-C-T-G haplotype (χ²=9.766, P=0.001785, OR[95%CI]: 5.654 [1.664-19.211]). In female-group, positive association was found for the A-A-C-G haplotype of rs3752826-rs2131431-rs1873827-rs12452627 (χ²=78.628, P=7.94E-19, OR[95%CI]: 50.043 [11.087-225.876]), its A-C-T-G haplotype (χ²=38.806, P=4.83E-10, OR[95%CI]: 0.053 [0.015-0.192]), the C-A-C-G haplotype (χ²=18.930, P=1.37E-05, OR[95%CI]: 0.526 [0.392-0.705]), and the C-C-T-G haplotype (χ²=38.668, P=5.18E-10, OR[95%CI]: 6.130 [3.207-11.716]). Our findings support YWHAE being a risk gene for Major Depressive Disorder in the Han Chinese population.

Associations between arsenic (+3 oxidation state) methyltransferase (AS3MT) and N-6 adenine-specific DNA methyltransferase 1 (N6AMT1) polymorphisms, arsenic metabolism, and cancer risk in a chilean population

Inter-individual differences in arsenic metabolism have been linked to arsenic-related disease risks. Arsenic (+3) methyltransferase (AS3MT) is the primary enzyme involved in arsenic metabolism, and we previously demonstrated in vitro that N-6 adenine-specific DNA methyltransferase 1 (N6AMT1) also methylates the toxic inorganic arsenic (iAs) metabolite, monomethylarsonous acid (MMA), to the less toxic dimethylarsonic acid (DMA). Here, we evaluated whether AS3MT and N6AMT1 gene polymorphisms alter arsenic methylation and impact iAs-related cancer risks. We assessed AS3MT and N6AMT1 polymorphisms and urinary arsenic metabolites (%iAs, %MMA, %DMA) in 722 subjects from an arsenic-cancer case-control study in a uniquely exposed area in northern Chile. Polymorphisms were genotyped using a custom designed multiplex, ligation-dependent probe amplification (MLPA) assay for 6 AS3MT SNPs and 14 tag SNPs in the N6AMT1 gene. We found several AS3MT polymorphisms associated with both urinary arsenic metabolite profiles and cancer risk. For example, compared to wildtypes, individuals carrying minor alleles in AS3MT rs3740393 had lower %MMA (mean difference = -1.9%, 95% CI: -3.3, -0.4), higher %DMA (mean difference = 4.0%, 95% CI: 1.5, 6.5), and lower odds ratios for bladder (OR = 0.3; 95% CI: 0.1-0.6) and lung cancer (OR = 0.6; 95% CI: 0.2-1.1). Evidence of interaction was also observed for both lung and bladder cancer between these polymorphisms and elevated historical arsenic exposures. Clear associations were not seen for N6AMT1. These results are the first to demonstrate a direct association between AS3MT polymorphisms and arsenic-related internal cancer risk. This research could help identify subpopulations that are particularly vulnerable to arsenic-related disease. Environ. Mol. Mutagen. 58:411-422, 2017. © 2017 Wiley Periodicals, Inc.

Detecting genetic association through shortest paths in a bidirected graph

Genome-wide association studies (GWASs) commonly use marginal association tests for each single-nucleotide polymorphism (SNP). Because these tests treat SNPs as independent, their power will be suboptimal for detecting SNPs hidden by linkage disequilibrium (LD). One way to improve power is to use a multiple regression model. However, the large number of SNPs preclude simultaneous fitting with multiple regression, and subset regression is infeasible because of an exorbitant number of candidate subsets. We therefore propose a new method for detecting hidden SNPs having significant yet weak marginal association in a multiple regression model. Our method begins by constructing a bidirected graph locally around each SNP that demonstrates a moderately sized marginal association signal, the focal SNPs. Vertexes correspond to SNPs, and adjacency between vertexes is defined by an LD measure. Subsequently, the method collects from each graph all shortest paths to the focal SNP. Finally, for each shortest path the method fits a multiple regression model to all the SNPs lying in the path and tests the significance of the regression coefficient corresponding to the terminal SNP in the path. Simulation studies show that the proposed method can detect susceptibility SNPs hidden by LD that go undetected with marginal association testing or with existing multivariate methods. When applied to real GWAS data from the Alzheimer's Disease Neuroimaging Initiative (ADNI), our method detected two groups of SNPs: one in a region containing the apolipoprotein E (APOE) gene, and another in a region close to the semaphorin 5A (SEMA5A) gene.

Genome-Enabled Prediction of Breeding Values for Feedlot Average Daily Weight Gain in Nelore Cattle

Nelore is the most economically important cattle breed in Brazil, and the use of genetically improved animals has contributed to increased beef production efficiency. The Brazilian beef feedlot industry has grown considerably in the last decade, so the selection of animals with higher growth rates on feedlot has become quite important. Genomic selection (GS) could be used to reduce generation intervals and improve the rate of genetic gains. The aim of this study was to evaluate the prediction of genomic-estimated breeding values (GEBV) for average daily weight gain (ADG) in 718 feedlot-finished Nelore steers. Analyses of three Bayesian model specifications [Bayesian GBLUP (BGBLUP), BayesA, and BayesCπ] were performed with four genotype panels [Illumina BovineHD BeadChip, TagSNPs, and GeneSeek High- and Low-density indicus (HDi and LDi, respectively)]. Estimates of Pearson correlations, regression coefficients, and mean squared errors were used to assess accuracy and bias of predictions. Overall, the BayesCπ model resulted in less biased predictions. Accuracies ranged from 0.18 to 0.27, which are reasonable values given the heritability estimates (from 0.40 to 0.44) and sample size (568 animals in the training population). Furthermore, results from Bos taurus indicus panels were as informative as those from Illumina BovineHD, indicating that they could be used to implement GS at lower costs.

Association between the zinc finger protein 804A (ZNF804A) gene and the risk of schizophrenia and bipolar I disorder across diagnostic boundaries

OBJECTIVES

In this study, we aimed to determine the role of genetic variations within the zinc finger protein 804A (ZNF804A) gene, a candidate for a psychosis risk-conferring gene, in the development of schizophrenia (SZ) and bipolar disorder (BP) in the Korean population.

METHODS

A total of 921 patients with SZ, bipolar I (BP-I) and II (BP-II) disorder, and 502 control subjects participated in the study. Twenty-one tag single nucleotide polymorphisms (SNPs) across the genomic region of ZNF804A and seven reference SNPs based on previous reports were genotyped. We applied logistic regression analyses under additive, dominant and recessive models.

RESULTS

Fifteen of the 28 SNPs showed a nominally significant association with at least one diagnostic group. However, none of these associations remained significant after false discovery rate (FDR) correction. As the trend of association was observed mostly in SZ and BP-I with similar patterns, we performed a post hoc analysis for the combined SZ and BP-I group. Five SNPs (rs2369595, rs6755404, rs10931156, rs12476147 and rs1366842) showed a significant association with an FDR-corrected P of <.05.

CONCLUSIONS

This study supports a possible role of ZNF804A in the common susceptibility of major psychoses, and identified additional candidate variants of the gene in the Korean population.

NOTCH4 gene polymorphisms as potential risk factors for brain arteriovenous malformation development and hemorrhagic presentation

OBJECTIVE

Arteriovenous malformations (AVMs) of the brain are a frequent and important cause of intracranial hemorrhage in young adults. Little is known about the molecular-genetic pathomechanisms underlying AVM development. Genes of the NOTCH family control the normal development of vessels and proper arteriovenous specification. Transgenic mice with constitutive expression of active NOTCH4 frequently develop AVMs. Here, the authors report a genetic association study investigating possible associations between NOTCH4 gene polymorphisms and formation and clinical presentation of AVMs.

METHODS

After PCR amplification and direct DNA sequencing or restriction digests, 10 single-nucleotide polymorphisms (SNPs) of the NOTCH4 gene were used for genotyping 153 AVM patients and 192 healthy controls (i.e., blood donors). Pertinent clinical data were available for 129 patients. Uni- and multivariate single-marker and explorative haplotype analyses were performed to identify potential genetic risk factors for AVM development and for hemorrhagic or epileptic presentation.

RESULTS

Eleven calculated haplotypes consisting of 3–4 SNPs (most of which were located in the epidermal growth factor–like domain of the NOTCH4 gene) were observed significantly more often among AVM patients than among controls. Univariate analysis indicated that rs443198_TT and rs915895_AA genotypes both were significantly associated with hemorrhage and that an rs1109771_GG genotype was associated with epilepsy. The association between rs443198_TT and AVM bleeding remained significant in the multivariate regression analysis.

CONCLUSIONS

The authors' results suggest NOTCH4 SNPs as possible genetic risk factors for the development and clinical presentation of AVMs and a role of NOTCH4 in the pathogenesis of this disease.

An adaptive threshold determination method of feature screening for genomic selection

Background

Although the dimension of the entire genome can be extremely large, only a parsimonious set of influential SNPs are correlated with a particular complex trait and are important to the prediction of the trait. Efficiently and accurately selecting these influential SNPs from millions of candidates is in high demand, but poses challenges. We propose a backward elimination iterative distance correlation (BE-IDC) procedure to select the smallest subset of SNPs that guarantees sufficient prediction accuracy, while also solving the unclear threshold issue for traditional feature screening approaches.

Results

Verified through six simulations, the adaptive threshold estimated by the BE-IDC performed uniformly better than fixed threshold methods that have been used in the current literature. We also applied BE-IDC to an Arabidopsis thaliana genome-wide data. Out of 216,130 SNPs, BE-IDC selected four influential SNPs, and confirmed the same FRIGIDA gene that was reported by two other traditional methods.

Conclusions

BE-IDC accommodates both the prediction accuracy and the computational speed that are highly demanded in the genomic selection.

Electronic supplementary material

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

Stress-related psychological symptoms contribute to axial pain persistence after motor vehicle collision: path analysis results from a prospective longitudinal study

Posttraumatic stress disorder (PTSD) symptoms and pain after traumatic events such as motor vehicle collision (MVC) have been proposed to be mutually promoting. We performed a prospective multicenter study that enrolled 948 individuals who presented to the emergency department within 24 hours of MVC and were discharged home after evaluation. Follow-up evaluations were completed 6 weeks, 6 months, and 1 year after MVC. Path analysis results supported the hypothesis that axial pain after MVC consistently promotes the maintenance of hyperarousal and intrusive symptoms, from the early weeks after injury through 1 year. In addition, path analysis results supported the hypothesis that one or more PTSD symptom clusters had an influence on axial pain outcomes throughout the year after MVC, with hyperarousal symptoms most influencing axial pain persistence in the initial months after MVC. The influence of hyperarousal symptoms on pain persistence was only present among individuals with genetic vulnerability to stress-induced pain, suggesting specific mechanisms by which hyperarousal symptoms may lead to hyperalgesia and allodynia. Further studies are needed to better understand the specific mechanisms by which pain and PTSD symptoms enhance one another after trauma, and how such mechanisms vary among specific patient subgroups, to better inform the development of secondary preventive interventions.

Associations Between Genetic Variants of NADPH Oxidase-Related Genes and Blood Pressure Responses to Dietary Sodium Intervention: The GenSalt Study

BACKGROUND

The aim of this study was to comprehensively test the associations of genetic variants of nicotinamide adenine dinucleotide phosphate (NADPH) oxidase-related genes with blood pressure (BP) responses to dietary sodium intervention in a Chinese population.

METHODS

We conducted a 7-day low-sodium intervention followed by a 7-day high-sodium intervention among 1,906 participants in rural China. BP measurements were obtained at baseline and each dietary intervention using a random-zero sphygmomanometer. Linear mixed-effect models were used to assess the additive associations of 63 tag single-nucleotide polymorphisms in 11 NADPH oxidase-related genes with BP responses to dietary sodium intervention. Gene-based analyses were conducted using the truncated product method. The Bonferroni method was used to adjust for multiple testing in all analyses.

RESULTS

Systolic BP (SBP) response to high-sodium intervention significantly decreased with the number of minor T allele of marker rs6967221 in RAC1 (P = 4.51 × 10-4). SBP responses (95% confidence interval) for genotypes CC, CT, and TT were 5.03 (4.71, 5.36), 4.20 (3.54, 4.85), and 0.56 (-1.08, 2.20) mm Hg, respectively, during the high-sodium intervention. Gene-based analyses revealed that RAC1 was significantly associated with SBP response to high-sodium intervention (P = 1.00 × 10-6) and diastolic BP response to low-sodium intervention (P = 9.80 × 10-4).

CONCLUSIONS

These findings suggested that genetic variants of NADPH oxidase-related genes may contribute to the variation of BP responses to sodium intervention in Chinese population. Further replication of these findings is warranted.

Association of Topoisomerase II (TOP2A) and Dual-Specificity Phosphatase 6 (DUSP6) Single Nucleotide Polymorphisms with Radiation Treatment Response and Prognosis of Lung Cancer in Han Chinese

Background

Mutations of DNA topoisomerase II (TOP2A) are associated with chemotherapy resistance, whereas dual-specificity phosphatase 6 (DUSP6) negatively regulates members of the mitogen-activated protein (MAP) kinase superfamily to control cell proliferation. This study assessed TOP2A and DUSP6 single nucleotide polymorphisms (SNPs) in non-small cell lung cancer (NSCLC) patients for association with chemoradiotherapy responses and prognosis.

Material/Methods

A total of 140 Chinese patients with histologically confirmed NSCLC were enrolled and subjected to genotyping of TOP2A rs471692 and DUSP6 rs2279574 using Taqman PCR. An independent sample t test was used to analyze differences in tumor regression after radiotherapy versus SNP risk factors. Kaplan-Meier curves analyzed overall survival, followed by the log-rank test and Cox proportional hazard models.

Results

There were no significant associations of TOP2A rs471692 and DUSP6 rs2279574 polymorphisms or clinicopathological variables with response to chemoradiotherapy (p>0.05). Comparing overall survival of 87 patients with stage I–III NSCLC treated with radiotherapy or chemoradiotherapy to clinicopathological variables, the data showed that tumor regression, weight loss, clinical stage, and cigarette smoking were independent prognostic predictors (p=0.009, 0.043, 0.004, and 0.025, respectively). Tumor regression rate >0.34 was associated with patent survival versus tumor regression rate ≤0.34 (p=0.007).

Conclusions

TOP2A rs471692 and DUSP6 rs2279574 SNPs were not associated with chemoradiotherapy response, whereas tumor regression, weight loss, clinical stage, and cigarette smoking were independent prognostic predictors for these Chinese patients with NSCLC.

Interaction Landscape of Inherited Polymorphisms with Somatic Events in Cancer

Recent studies have characterized the extensive somatic alterations that arise during cancer. However, the somatic evolution of a tumor may be significantly affected by inherited polymorphisms carried in the germline. Here, we analyze genomic data for 5,954 tumors to reveal and systematically validate 412 genetic interactions between germline polymorphisms and major somatic events, including tumor formation in specific tissues and alteration of specific cancer genes. Among germline-somatic interactions, we found germline variants in RBFOX1 that increased incidence of SF3B1 somatic mutation by 8-fold via functional alterations in RNA splicing. Similarly, 19p13.3 variants were associated with a 4-fold increased likelihood of somatic mutations in PTEN. In support of this association, we found that PTEN knockdown sensitizes the MTOR pathway to high expression of the 19p13.3 gene GNA11 Finally, we observed that stratifying patients by germline polymorphisms exposed distinct somatic mutation landscapes, implicating new cancer genes. This study creates a validated resource of inherited variants that govern where and how cancer develops, opening avenues for prevention research.Significance: This study systematically identifies germline variants that directly affect tumor evolution, either by dramatically increasing alteration frequency of specific cancer genes or by influencing the site where a tumor develops. Cancer Discovery; 7(4); 410-23. ©2017 AACR.See related commentary by Geeleher and Huang, p. 354This article is highlighted in the In This Issue feature, p. 339.

THE PPARGC1A - GLY482SER POLYMORPHISM (RS8192678) AND THE METABOLIC SYNDROME IN A CENTRAL ROMANIAN POPULATION

Background

The peroxisome proliferator-activated receptor-γ co-activator 1-α (PPARGC1A), a key transcription factor involved in the control of metabolism and energy homeostasis, is an important biological and positional candidate of the metabolic syndrome. Association studies of its polymorphisms, however, yielded inconsistent sometimes conflicting results, pointing to important ethnic differences, which call for replication in various populations.

Objective

In order to study its most common - potentially functional - polymorphism Gly482Ser (rs8192678), we carried out a case-control study in a central Romanian population.

Material and methods

Two hundred and ninety six patients affected by the metabolic syndrome diagnosed according to the International Diabetes Federation proposed criteria and 166 middle-aged control subjects have been investigated. Genotyping was done by PCR-RFLP, using the restriction enzyme MspI.

Results

While the G(Gly)/A(Ser) allele frequencies (66.89/33.11 vs. 71.68/28.31 %) and GG/GA/AA genotype distribution (45.27-43.24-11.48 vs. 54.21-34.93-10.84 %) differed in the metabolic syndrome and control group, the risk of developing the metabolic syndrome did not reach the limit of statistical significance (OR=1.43; p=0.06, CI 95%: 0.97-2.09). Metabolic parameters in the two study groups did not show significant differences according to the genotype (p>0.05).

Conclusion

rs8192678 could be a functional polymorphism contributing to the development of the metabolic syndrome, but probably its effect is minor, and might depend on gene-gene and gene-environment interactions. Clarification of very small effects would require larger sample sizes.

Genetic variation in leptin and leptin receptor genes is a risk factor for idiopathic recurrent spontaneous abortion

AIM

To determine whether maternal leptin (LEP) and leptin receptor (LEPR) gene polymorphisms are associated with idiopathic recurrent spontaneous abortion (IRSA).

METHODS

This case-control association study conducted from 2010 to 2012 at the Department of Gynecology and Obstetrics, University Hospital Center Osijek and Clinical Institute of Medical Genetics Ljubljana included 178 women with a history of three or more IRSAs before the 22nd week of gestation and 145 women with at least two live births and no history of pathologic pregnancies during reproductive period. Polymorphisms of maternal LEP (rs7799039, rs2122627, rs11761556, rs10244329) and LEPR (rs1137101, rs7516341, rs1186403, rs12062820) were assessed by allele specific real-time polymerase chain reaction. Genotype distribution, allele frequencies, and frequency of haplotypes at LEP and LEPR genetic loci were determined.

RESULTS

We observed more frequent genotype for rs7516341 (nominal P=0.034, odds ratio [OR] 0.61, 95% confidence interval [CI] 0.38-0.97) and rs1137101 (nominal P=0.048, OR 1.66, 95% CI 1.00-2.80) in the LEPR gene in patients than in controls, but these results did not remain significant after correction for multiple testing according to Bonferroni (adjusted P value threshold was set at 0.05). We did not observe differential distribution of genotype frequencies in the LEP gene between cases and controls. In patients with IRSA, GTCC haplotype in the LEPR gene locus was significantly less frequent than in controls (PP=0.00865, OR 0.45), contrary to ACTC haplotype (PP=0.0087, OR 1.98).

CONCLUSIONS

We showed that genetic variability in the LEPR gene was associated with IRSA, warranting confirmation on a greater number of patients and pathogenesis investigation.

Quercetin intake, MATE1 polymorphism, and metabolic syndrome in Korean population: Hallym aging study

Multidrug and toxic compound extrusion transporter-1 (MATE1) is a quercetin transporter. We examined the associations of quercetin intake and polymorphism of MATE1 in relation to metabolic syndrome (MetS) in Hallym Aging Study. Quercetin intake and the measurements for MetS were assessed in 2004. Six tagging single nucleotide polymorphisms (SNPs) at MATE1 gene were genotyped in 428 Korean adults in 2012. We found a lower prevalence of MetS with quercetin intake; compared to the lowest quartile, odds ratios (ORs, 95% confidence intervals; CIs) were 0.44 (0.24-0.84) for the 3^rd quartile. Individuals with the minor allele of MATE1, rs2453589, tended to have a lower prevalence of MetS compared to those with the major allele (OR=0.69; CI=0.36-1.34). However, interactions between quercetin intake and six MATE1 polymorphisms in relation to MetS were not significant (p for interaction ≥0.37). In conclusion, intake of quercetin was associated with MetS in Korean populations.

Modification of the association between early adversity and obsessive-compulsive disorder by polymorphisms in the MAOA, MAOB and COMT genes

The monoamine oxidases (MAOA/B) and catechol-O-methyltransferase (COMT) enzymes break down regulatory components within serotonin and dopamine pathways, and polymorphisms within these genes are candidates for OCD susceptibility. Childhood trauma has been linked OCD psychopathology, but little attention has been paid to the interactions between genes and environment in OCD aetiology. This pilot study investigated gene-by-environment interactions between childhood trauma and polymorphisms in the MAOA, MAOB and COMT genes in OCD. Ten polymorphisms (MAOA: 3 variants, MAOB: 4 variants, COMT: 3 variants) were genotyped in a cohort of OCD patients and controls. Early-life trauma was assessed using the Childhood Trauma Questionnaire (CTQ). Gene-by-gene (GxG) and gene-by-environment interactions (GxE) of the variants and childhood trauma were assessed using logistic regression models. Significant GxG interactions were found between rs362204 (COMT) and two independent polymorphisms in the MAOB gene (rs1799836 and rs6651806). Haplotype associations for OCD susceptibility were found for MAOB. Investigation of GxE interactions indicated that the sexual abuse sub-category was significantly associated with all three genes in haplotype x environment interaction analyses. Preliminary findings indicate that polymorphisms within the MAOB and COMT genes interact resulting in risk for OCD. Childhood trauma interacts with haplotypes in COMT, MAOA and MAOB, increasing risk for OCD.

Discovering Genome-Wide Tag SNPs Based on the Mutual Information of the Variants

Exploring linkage disequilibrium (LD) patterns among the single nucleotide polymorphism (SNP) sites can improve the accuracy and cost-effectiveness of genomic association studies, whereby representative (tag) SNPs are identified to sufficiently represent the genomic diversity in populations. There has been considerable amount of effort in developing efficient algorithms to select tag SNPs from the growing large-scale data sets. Methods using the classical pairwise-LD and multi-locus LD measures have been proposed that aim to reduce the computational complexity and to increase the accuracy, respectively. The present work solves the tag SNP selection problem by efficiently balancing the computational complexity and accuracy, and improves the coverage in genomic diversity in a cost-effective manner. The employed algorithm makes use of mutual information to explore the multi-locus association between SNPs and can handle different data types and conditions. Experiments with benchmark HapMap data sets show comparable or better performance against the state-of-the-art algorithms. In particular, as a novel application, the genome-wide SNP tagging is performed in the 1000 Genomes Project data sets, and produced a well-annotated database of tagging variants that capture the common genotype diversity in 2,504 samples from 26 human populations. Compared to conventional methods, the algorithm requires as input only the genotype (or haplotype) sequences, can scale up to genome-wide analyses, and produces accurate solutions with more information-rich output, providing an improved platform for researchers towards the subsequent association studies.

The role of GRIP1 and ephrin B3 in blood pressure control and vascular smooth muscle cell contractility

Several erythropoietin-producing hepatocellular receptor B family (EPHB) and their ligands, ephrinBs (EFNBs), are involved in blood pressure regulation in animal models. We selected 528 single nucleotide polymorphisms (SNPs) within the genes of EPHB6, EFNB2, EFNB3 and GRIP1 in the EPH/EFN signalling system to query the International Blood Pressure Consortium dataset. A SNP within the glutamate receptor interacting protein 1 (GRIP1) gene presented a p-value of 0.000389, approaching the critical p-value of 0.000302, for association with diastolic blood pressure of 60,396 individuals. According to echocardiography, we found that Efnb3 gene knockout mice showed enhanced constriction in the carotid arteries. In vitro studies revealed that in mouse vascular smooth muscle cells, siRNA knockdown of GRIP1, which is in the EFNB3 reverse signalling pathway, resulted in increased contractility of these cells. These data suggest that molecules in the EPHB/EFNB signalling pathways, specifically EFNB3 and GRIP1, are involved blood pressure regulation.

Evaluation of association of common variants in HTR1A and HTR5A with schizophrenia and executive function

The 5-HT1A receptor (HTR1A) and the 5-HT5A receptor (HTR5A) are key 5-HT receptors with distinct inhibitory functions. Studies have been conducted to investigate the association of a few HTR1A polymorphisms with schizophrenia, producing conflicting results, and the relationship between HTR5A and schizophrenia has not yet been well investigated. We aimed to examine the association of HTR1A and HTR5A with schizophrenia and executive function. The study included a discovery stage with 1,115 patients and 2,289 controls and a replication stage with 2,128 patients and 3,865 controls. A total of 30 common SNPs in HTR1A and HTR5A were genotyped in the discovery stage, and significantly associated SNPs were genotyped in the replication stage. We identified that two SNPs (rs878567 in HTR1A and rs1800883 in HTR5A) were significantly associated with schizophrenia in both datasets, and similar results were observed in imputation and haplotype association analyses. Moreover, we found that SNP rs1800883 significantly interacted with executive function when processing the perseverative error of Wisconsin Card Sorting Test in patients. Our results provide further supportive evidence of the effect of HTR1A and HTR5A on the etiology of schizophrenia and suggest that the selected genetic variations in HTR5A may be involved in impaired executive function.

Inherited Variants in Wnt Pathway Genes Influence Outcomes of Prostate Cancer Patients Receiving Androgen Deprivation Therapy

Aberrant Wnt signaling has been associated with many types of cancer. However, the association of inherited Wnt pathway variants with clinical outcomes in prostate cancer patients receiving androgen deprivation therapy (ADT) has not been determined. Here, we comprehensively studied the contribution of common single nucleotide polymorphisms (SNPs) in Wnt pathway genes to the clinical outcomes of 465 advanced prostate cancer patients treated with ADT. Two SNPs, adenomatous polyposis coli (APC) rs2707765 and rs497844, were significantly (p ≤ 0.009 and q ≤ 0.043) associated with both prostate cancer progression and all-cause mortality, even after multivariate analyses and multiple testing correction. Patients with a greater number of favorable alleles had a longer time to disease progression and better overall survival during ADT (p for trend ≤ 0.003). Additional, cDNA array and in silico analyses of prostate cancer tissue suggested that rs2707765 affects APC expression, which in turn is correlated with tumor aggressiveness and patient prognosis. This study identifies the influence of inherited variants in the Wnt pathway on the efficacy of ADT and highlights a preclinical rationale for using APC as a prognostic marker in advanced prostate cancer.

Genetic variation in leptin and leptin receptor genes as a risk factor for idiopathic male infertility

The aim of this study was to examine whether there is an association among genetic variability in leptin (LEP) and leptin receptor (LEPR) genes and male infertility. We performed a case-control study and were searching for an association between polymorphisms of LEP and LEPR genes and male infertility. The study group consisted of 317 patients with idiopathic infertility and a control group of 241 fertile men from Slovenia. Four single nucleotide polymorphisms (SNPs) in LEP gene and four single nucleotide polymorphisms (SNPs) in LEPR gene were chosen and genotyped. Statistically significant SNP was further validated in additional 255 infertile patients and 168 controls from Serbia and Macedonia. In the Slovenian population, we found a statistically significant difference in genotype distribution for rs10244329 polymorphism in LEP gene (recessive genotype model, p value = 0.048). The trend toward statistically significant difference in genotype distribution for rs10244329 polymorphism was confirmed in the Serbian and Macedonian populations (p value = 0.07). Our data suggest that genetic variability in the LEP gene might be associated with male infertility warranting further confirmation and mechanistic investigations.

Integrating haplotype-specific linkage maps in tetraploid species using SNP markers

Linkage mapping can help unravel the complexities of polyploid genomes. Here, we integrate haplotype-specific linkage maps in autotetraploid potato and explore the possibilities for mapping in other polyploid species. High-density linkage mapping in autopolyploid species has become possible in recent years given the increasing number of molecular markers now available through modern genotyping platforms. Such maps along with larger experimental populations are needed before we can obtain sufficient accuracy to make marker-trait association studies useful in practice. Here, we describe a method to create genetic linkage maps for an autotetraploid species with large numbers of markers and apply it to an F1 population of tetraploid potato (Solanum tuberosum L.) of 235 individuals genotyped using a 20K SNP array. SNP intensity values were converted to allele dosages after which we calculated pairwise maximum likelihood estimates of recombination frequencies between all marker segregation types under the assumption of random bivalent pairing. These estimates were used in the clustering of markers into linkage groups and their subsequent ordering into 96 homologue maps. The homologue maps were integrated per chromosome, resulting in a total map length of 1061 cM from 6910 markers covering all 12 potato chromosomes. We examined the questions of marker phasing and binning and propose optimal strategies for both. We also investigated the effect of quadrivalent formation and preferential pairing on recombination frequency estimation and marker phasing, which is of great relevance not only for potato but also for genetic studies in other tetraploid species for which the meiotic pairing behaviour is less well understood.

Interactive Influence of N6AMT1 and As3MT Genetic Variations on Arsenic Metabolism in the Population of Inner Mongolia, China

Chronic arsenic exposure via drinking water has become a worldwide public health concern. In humans, inorganic arsenic (iAs) is metabolized to monomethylarsonic acid (MMA) and dimethylarsinic acid (DMA) mainly mediated by arsenic (+3 oxidation state) methyltransferase (As3MT). We reported recently that N-6 adenine-specific DNA methyltransferase 1 (N6AMT1) was involved in arsenic metabolism, and examined its interactive effect with As3MT on arsenic metabolism in vitro To further evaluate the interactive effect of N6AMT1 and As3MT on arsenic biomethylation in humans, we conducted a human population-based study including 289 subjects living in rural villages in Inner Mongolia, China, and assessed their urinary arsenic metabolites profiles in relation to genetic polymorphisms and haplotypes of N6AMT1 and As3MT Five N6AMT1 single nucleotide polymorphisms (SNPs; rs1003671, rs7282257, rs2065266, rs2738966, rs2248501) and the N6AMT1 haplotype 2_GGCCAT were significantly associated with the percentage of iAs (% iAs) in urine (e.g., for rs7282257, mean was 9.62% for TT, 6.73% for AA). Rs1003671 was also in a significant relationship with urinary MMA and DMA (the mean of %MMA was 24.95% for GA, 31.69% for GG; the mean of % DMA was 69.21% for GA, 59.82% for GG). The combined effect of N6AMT1 haplotype 2_GGCCAT and As3MT haplotype 2_GCAC showed consistence with the additive significance of each haplotype on % iAs: the mean was 5.47% and 9.36% for carriers with both and null haplotypes, respectively. Overall, we showed that N6AMT1 genetic polymorphisms were associated with arsenic biomethylation in the Chinese population, and its interaction with As3MT was observed in specific haplotype combinations.

Determinants of acquired activated protein C resistance and D-dimer in breast cancer

BACKGROUND

We have previously reported acquired activated protein C (APC) resistance and elevated plasma D-dimer levels in breast cancer patients. Here, we aimed to identify phenotypic and genetic determinants that contribute to the acquired APC resistance and increased D-dimer levels in breast cancer. Healthy controls served as reference. We also addressed whether higher APC resistance or D-dimer levels could be potential markers of clinicopathological breast cancer characteristics.

MATERIALS AND METHODS

358 breast cancer patients and 273 healthy controls were enrolled and hemostatic plasma parameters were determined; factor (F) V, FVIII, FIX, FX, fibrinogen, von Willebrand factor (VWF), normalized APC sensitivity ratio (n-APC-sr), protein C, protein S, antithrombin, tissue factor pathway inhibitor (TFPI), and D-dimer. Common single nucleotide polymorphisms were genotyped in coagulation-related genes in the breast cancer patients.

RESULTS

The phenotypic hemostatic factors explained 25% and 31% of the variability in acquired APC resistance and D-dimer levels, respectively, in the breast cancer patients. Fibrinogen (β=-0.35, P<0.001), protein C (β=0.28, P<0.001), and FIX (β=0.22, P=0.026) were identified as determinants of n-APC-sr (in FV Leiden non-carriers), whereas TFPI (β=0.28, P<0.001), antithrombin (β=-0.25, P<0.001), and FX (β=0.15, P=0.040) were the major determinants of D-dimer. Moreover, borderline higher APC resistance (>75th percentile) was found in patients with triple negative tumors (odds ratio (OR) 1.97, 95% CI 0.99-3.90).

CONCLUSIONS

This study reports phenotypic hemostatic parameters that determine acquired APC resistance and D-dimer levels in breast cancer patients. The explanatory power was modest, however, our findings are hypothesis generating and may contribute to further understand the background for cancer associated-coagulopathy and thrombosis.

An enzyme in the kynurenine pathway that governs vulnerability to suicidal behavior by regulating excitotoxicity and neuroinflammation

Emerging evidence suggests that inflammation has a key role in depression and suicidal behavior. The kynurenine pathway is involved in neuroinflammation and regulates glutamate neurotransmission. In the cerebrospinal fluid (CSF) of suicidal patients, levels of inflammatory cytokines and the kynurenine metabolite quinolinic acid (QUIN), an N-methyl-d-aspartate receptor agonist, are increased. The enzyme amino-β-carboxymuconate-semialdehyde-decarboxylase (ACMSD) limits QUIN formation by competitive production of the neuroprotective metabolite picolinic acid (PIC). Therefore, decreased ACMSD activity can lead to excess QUIN. We tested the hypothesis that deficient ACMSD activity underlies suicidal behavior. We measured PIC and QUIN in CSF and plasma samples from 137 patients exhibiting suicidal behavior and 71 healthy controls. We used DSM-IV and the Montgomery-Åsberg Depression Rating Scale and Suicide Assessment Scale to assess behavioral changes. Finally, we genotyped ACMSD tag single-nucleotide polymorphisms (SNPs) in 77 of the patients and 150 population-based controls. Suicide attempters had reduced PIC and a decreased PIC/QUIN ratio in both CSF (P<0.001) and blood (P=0.001 and P<0.01, respectively). The reductions of PIC in CSF were sustained over 2 years after the suicide attempt based on repeated measures. The minor C allele of the ACMSD SNP rs2121337 was more prevalent in suicide attempters and associated with increased CSF QUIN. Taken together, our data suggest that increased QUIN levels may result from reduced activity of ACMSD in suicidal subjects. We conclude that measures of kynurenine metabolites can be explored as biomarkers of suicide risk, and that ACMSD is a potential therapeutic target in suicidal behavior.

Association between AUTS2 haplotypes and alcohol dependence in a Japanese population

OBJECTIVE

Recent genome-wide analysis has indicated that the autism susceptibility candidate 2 (AUTS2) gene is involved in the regulation of alcohol consumption. We hypothesised that AUTS2 might be associated with the development of alcohol dependence. Therefore, in this exploratory study, we compared the genotype and allele frequencies of the polymorphisms rs6943555 and rs9886351 in the AUTS2 gene between patients with alcohol dependence and healthy control subjects living in a Japanese provincial prefecture. We also examined whether or not the haplotypes consisting of these polymorphisms are related to alcohol dependence.

METHODS

The subjects of this study consisted of 64 patients with alcohol dependence and 75 unrelated healthy people. The AUTS2 genotypes were determined by the polymerase chain reaction (PCR)-restriction fragment length polymorphism (RFLP) method.

RESULTS

No significant differences in the genotype and allele frequencies of the polymorphisms AUTS2 rs6943555 and rs9886351 were found between alcohol dependence and control subjects. On the other hand, the frequencies of the AUTS2 haplotypes were significantly different between them, and the rs6943555 and rs9886351 A-A haplotype was associated with alcohol dependence (p=0.0187).

CONCLUSION

This suggests that the rs6943555 and rs9886351 A-A haplotype might affect the vulnerability to alcohol dependence pathogenesis. Further studies are needed to confirm the reproducibility of the results of this study with increased numbers of subjects.

Preeclampsia does not share common risk alleles in 9p21 with coronary artery disease and type 2 diabetes

INTRODUCTION

Preeclampsia is a common and partially genetic pregnancy complication characterized by hypertension and proteinuria. Association with cardiovascular disease and type 2 diabetes has been reported in 9p21 by several genome-wide association studies. It has been hypothesized that cardiometabolic diseases may share common etiology with preeclampsia.

MATERIALS AND METHODS

We tested association with the 9p21 region to preeclampsia in the Finnish population by genotyping 23 tagging single nucleotide polymorphisms (SNPs) in 15 extended preeclampsia families and in a nationwide cohort consisting of 281 cases and 349 matched controls. Replication was conducted in additional datasets.

RESULTS

Four SNPs (rs7044859, rs496892, rs564398 and rs7865618) showed nominal association (p ≤ 0.024 uncorrected) with preeclampsia in the case-control cohort. To increase power, we genotyped two SNPs in additional 388 cases and 341 controls from the Finnish Genetics of Preeclampsia Consortium (FINNPEC) cohort. Partial replication was also attempted in a UK cohort (237 cases and 199 controls) and in 74 preeclamptic families from Australia/New Zealand. We were unable to replicate the initial association in the extended Finnish dataset or in the two international cohorts.

CONCLUSIONS

Our study did not find evidence for the involvement of the 9p21 region in the risk of preeclampsia. Key Message Chromosome 9p21 is not associated with preeclampsia.

CYP4F2 affects phenotypic outcome in adrenoleukodystrophy by modulating the clearance of very long-chain fatty acids

X-linked adrenoleukodystrophy (ALD) is a severe neurodegenerative disorder caused by the accumulation of very long-chain fatty acids (VLCFA) due to mutations in the ABCD1 gene. The phenotypic spectrum ranges from a fatal cerebral demyelinating disease in childhood (cerebral ALD) to a progressive myelopathy without cerebral involvement in adulthood (adrenomyeloneuropathy). Because ABCD1 mutations have no predictive value with respect to clinical outcome a role for modifier genes was postulated. We report that the CYP4F2 polymorphism rs2108622 increases the risk of developing cerebral ALD in Caucasian patients. The rs2108622 polymorphism (c.1297G>A) results in an amino acid substitution valine for methionine at position 433 (p.V433M). Using cellular models of VLCFA accumulation, we show that p.V433M decreases the conversion of VLCFA into very long-chain dicarboxylic acids by ω-oxidation, a potential escape route for the deficient peroxisomal β-oxidation of VLCFA in ALD. Although p.V433M does not affect the catalytic activity of CYP4F2 it reduces CYP4F2 protein levels markedly. These findings open perspectives for therapeutic interventions in a disease with currently limited treatment options.

Common Genetic Variation in CYP17A1 and Response to Abiraterone Acetate in Patients with Metastatic Castration-Resistant Prostate Cancer

Treatment with abiraterone acetate and prednisone (AA/P) prolongs survival in metastatic castration-resistant prostate cancer (mCRPC) patients. We evaluated the genetic variation in CYP17A1 as predictive of response to AA/P. A prospective collection of germline DNA prior to AA/P initiation and follow-up of a mCRPC cohort was performed. Five common single-nucleotide polymorphisms (SNPs) in CYP17A1 identified using a haplotype-based tagging algorithm were genotyped. Clinical outcomes included biochemical response and time to biochemical progression on AA/P. Logistic regression was used to assess the association between tag SNPs and biochemical response. Proportional hazards regression was used to assess the association between tag SNPs and time to biochemical progression. Odds or hazard ratio per minor allele were estimated and p-values below 0.05 were considered statistically significant. Germline DNA was successfully genotyped for four tag SNPs in 87 patients. The median age was 73 years (54-90); the median prostate-specific antigen was 66 ng/dL (0.1-99.9). A single SNP, rs2486758, was associated with lower odds of experiencing a biochemical response (Odds ratio 0.22, 95% confidence interval 0.07-0.63, p = 0.005) and a shorter time to biochemical progression (Hazard ratio 2.23, 95% confidence interval 1.39-3.56, p < 0.001). This tag SNP located in the promoter region of CYP17A1 will need further validation as a predictive biomarker for AA/P therapy.

Enabling Privacy-Preserving GWASs in Heterogeneous Human Populations

The proliferation of large genomic databases offers the potential to perform increasingly larger-scale genome-wide association studies (GWASs). Due to privacy concerns, however, access to these data is limited, greatly reducing their usefulness for research. Here, we introduce a computational framework for performing GWASs that adapts principles of differential privacy-a cryptographic theory that facilitates secure analysis of sensitive data-to both protect private phenotype information (e.g., disease status) and correct for population stratification. This framework enables us to produce privacy-preserving GWAS results based on EIGENSTRAT and linear mixed model (LMM)-based statistics, both of which correct for population stratification. We test our differentially private statistics, PrivSTRAT and PrivLMM, on simulated and real GWAS datasets and find they are able to protect privacy while returning meaningful results. Our framework can be used to securely query private genomic datasets to discover which specific genomic alterations may be associated with a disease, thus increasing the availability of these valuable datasets.

Polymorphisms of Insulin-Like Growth Factor 1 Pathway Genes and Breast Cancer Risk

Genetic variants of insulin-like growth factor 1 (IGF1) pathway genes have been shown to be associated with breast density and IGF1 levels and, therefore, may also influence breast cancer risk via pro-survival signaling cascades. The aim of this study was to investigate associations between IGF1 pathway single nucleotide polymorphisms (SNPs) and breast cancer risk among European and East Asian women, and potential interactions with menopausal status and breast tumor subtype. Stratified analyses of 1,037 cases and 1,050 controls from a population-based case–control study were conducted to assess associations with breast cancer for 22 SNPs across 5 IGF1 pathway genes in European and East Asian women. Odds ratios were calculated using logistic regression in additive genetic models. Polytomous logistic regression was used to assess heterogeneity by breast tumor subtype. Two SNPs of the IGF1 gene (rs1019731 and rs12821878) were associated with breast cancer risk among European women. Four highly linked IGF1 SNPs (rs2288378, rs17727841, rs7136446, and rs7956547) were modified by menopausal status among East Asian women only and associated with postmenopausal breast cancers. The association between rs2288378 and breast cancer risk was also modified by breast tumor subtype among East Asian women. Several IGF1 polymorphisms were found to be associated with breast cancer risk and some of these associations were modified by menopausal status or breast tumor subtype. Such interactions should be considered when assessing the role of these variants in breast cancer etiology.

Evaluating the Role of Genetic Variants on first-line antiepileptic drug response in North India: Significance of SCN1A and GABRA1 Gene Variants in Phenytoin Monotherapy and its Serum Drug Levels

AIM

The present study aimed to evaluate association of genetic variants on drug response and therapy optimization parameters in patients treated with first-line antiepileptic drugs (AEDs). Genetic variants from ion channels, their functionally related genes, and synaptic vesicle cycle (SVC) genes with a potential role in epilepsy pathophysiology were thus prioritized.

METHODS

A total of 12 genes from ion channels and related gene set and seven genes from SVC comprising 155 SNPs were genotyped and evaluated with drug response, dose levels, and drug levels in 408 patients with epilepsy.

RESULTS

Both GABRA1 and SCN1A variants showed haplotypic and diplotypic associations in response to phenytoin (PHT). Diplotype analysis of GABRA1 variants revealed association of rs12658835|rs7735530 (AG/AG) (P-valuecorrected  = 0.034, OR = 3.75, 95% CI = 1.36-11.05) and rs12658835|rs7735530|rs7732641|rs2279020 (AGCA/AGCA) (P-valuecorrected  = 0.035, OR = 2.48, 95% CI = 0.96-6.41) with recurrent seizures. SCN1A haplotype rs6432860|rs3812718 (AC: P-valuecorrected  = 0.022, OR = 2.72, 95% CI = 1.39-5.35) and diplotype (AC/AC: P-valuecorrected  = 0.034, OR = 6.42, 95% CI = 1.10-65.76) were further observed to be associated with recurrent seizures. With respect to therapy optimization parameters, we observed significantly lower dose-adjusted drug levels at maximum dose of PHT in patients carrying AC/AC diplotype (P-value = 0.021).

CONCLUSION

The results further substantiate the role of GABRA1 in PHT mode of action and contribution of SCN1A in response and therapy optimization with PHT monotherapy.

Association study of PNPLA2 gene with histological parameters of NAFLD in an obese population

INTRODUCTION

The prevalence of non-alcoholic fatty liver disease (NAFLD) and the closely associated metabolic syndrome is high and is related to risk factors such as obesity and type 2 diabetes. A genetic basis for NAFLD has been suggested, but only few causal genes have been identified. The most significant association reported to date is the robust association of the PNPLA3 I148M variant with susceptibility to NAFLD. We therefore hypothesized that the PNPLA2 gene might also be involved in NAFLD pathogenesis, because of its close sequence similarity with PNPLA3 and its possible involvement in ectopic fat accumulation.

METHODS

In this study, we investigated the association of PNPLA2 polymorphisms with the development of non-alcoholic fatty liver disease in a prospectively recruited Belgian obese population comprising 633 individuals with varying degrees of fatty liver disease. We selected 3 PNPLA2 SNPs for genotyping, including 2 tagSNPs that cover most information on common genetic variation in the selected region.

RESULTS

After performing linear regression analysis, we found that 2 of the analyzed PNPLA2 SNPs were associated with anthropometric and metabolic parameters. In our subcohort of patients that underwent liver biopsy (n=372/633 or 58.7%), we assessed the influence of the PNPLA2 variants on the severity of histologically determined liver damage, but we did not find convincing evidence for association.

CONCLUSION

Although we found evidence for moderate association between PNPLA2 tagSNPs and anthropometric and metabolic parameters in our cohort, no evidence for association between polymorphisms in the PNPLA2 gene and the presence and severity of NAFLD was identified.

Variants in Adjacent Oxytocin/Vasopressin Gene Region and Associations with ASD Diagnosis and Other Autism Related Endophenotypes

BACKGROUND

There has been increasing interest in oxytocin (peptide: OT, gene: OXT) as a treatment pathway for neurodevelopmental disorders such as Autism Spectrum Disorder (ASD). Neurodevelopmental disorders affect functional, social, and intellectual abilities. With advances in molecular biology, research has connected multiple gene regions to the clinical presentation of ASD. Studies have also shown that the neuropeptide hormones OT and arginine vasopressin (AVP) influence mammalian social and territorial behaviors and may have treatment potential for neurodevelopmental disorders. Published data examining molecular and phenotypic variation in ASD, such as cognitive abilities, are limited. Since most studies have focused on the receptors in the OT-AVP system, we investigated genetic variation within peptide genes for association with phenotypic ASD features that help identify subgroups within the spectrum.

METHODS

In this study, TDT analysis was carried out utilizing FBAT in 207 probands (156 trios) and a European Ancestry (EA) subsample (108 trios).The evolutionarily related and adjacent genes of OXT and AVP were studied for associations between the tagged single nucleotide polymorphisms and ASD diagnosis, social abilities, restrictive and repetitive behaviors, and IQ for cognitive abilities. Additionally, relationships with whole blood serotonin (WB5HT) were explored because of the developmental relationships connecting plasma levels of OT and WB5HT within ASD.

RESULTS

RESULTS indicate significant association between OXT rs6084258 (p = 0.001) and ASD. Associations with several endophenotypes were also noted: OXT rs6133010 was associated with IQ (full scale IQ, p = 0.008; nonverbal IQ, p = 0.010, verbal IQ, p = 0.006); and OXT rs4813625 and OXT rs877172 were associated with WB5HT levels (EA, p = 0.027 and p = 0.033, respectively). Additionally, we measured plasma OT (pOT) levels in a subsample (N = 54). RESULTS show the three polymorphisms, OXT rs6084258, OXT rs11697250, and OXT rs877172, have significant association with pOT (EA, p = 0.011, p = 0.010, and p = 0.002, respectively).

CONCLUSIONS

These findings suggest that SNPs near OXT and AVP are associated with diagnosis of ASD, social behaviors, restricted and repetitive behaviors, IQ, pOT, and WB5HT. Future studies need to replicate these findings and examine gene-interactions in other neurodevelopmental disorders. Mechanisms of action may influence early social and cognitive development that may or may not be limited to ASD diagnosis.

Investigation of common and rare genetic variation in the BAMBI genomic region in light of human obesity

The aim of this study was to confirm the previously identified link between BAMBI and human obesity by means of a genetic and functional analysis. We performed both a mutation analysis, using high-resolution melting curve analysis, and a genetic association study, including 8 common tagSNPs in the BAMBI gene region. Three of the identified genetic variants (R151W, H201R, and C229R) were evaluated for their Wnt signaling enhancing capacity in a Wnt luciferase reporter assay. Mutation screening of the BAMBI coding region and exon-intron boundaries on our population of 677 obese children and adolescents and 529 lean control subjects resulted in the identification of 18 variants, 10 of which were not previously reported and 12 of which were exclusively found in obese individuals. The difference in variant frequency, not taking into account common polymorphisms, between obese (3.1 %) and lean (0.9 %) subjects was statistically significant (p = 0.004). Our Wnt luciferase assay, using WT and mutant BAMBI constructs, showed a significantly reduced activity for all of the investigated variants. Logistic and linear regression analysis on our Caucasian population of 1022 obese individuals and 606 lean controls, did not identify associations with obesity parameters (p values >0.05). We found several rare genetic variations, which represent the first naturally occurring missense variants of BAMBI in obese patients. Three variants (R151W, H201R, and C229R) were shown to reduce Wnt signaling enhancing capacity of BAMBI and we believe this result should encourage further study of this gene in other obese populations. In addition, we did not find evidence for the involvement of BAMBI common variation in human obesity in our population.

Role played by the SP4 gene in schizophrenia and major depressive disorder in the Han Chinese population

BACKGROUND

Psychiatric disorders such as schizophrenia and major depressive disorder (MDD) are likely to be caused by multiple susceptibility genes, each with small effects in increasing the risk of illness. Identifying DNA variants associated with schizophrenia and MDD is a crucial step in understanding the pathophysiology of these disorders.

AIMS

To investigate whether the SP4 gene plays a significant role in schizophrenia or MDD in the Han Chinese population.

METHOD

We focused on nine single nucleotide polymorphisms (SNPs) harbouring the SP4 gene and carried out case-control studies in 1235 patients with schizophrenia, 1045 patients with MDD and 1235 healthy controls recruited from the Han Chinese population.

RESULTS

We found that rs40245 was significantly associated with schizophrenia in both allele and genotype distributions (Pallele = 0.0005, Pallele = 0.004 after Bonferroni correction; Pgenotype = 0.0023, Pgenotype = 0.0184 after Bonferroni correction). The rs6461563 SNP was significantly associated with schizophrenia in the allele distributions (Pallele = 0.0033, Pallele = 0.0264 after Bonferroni correction).

CONCLUSIONS

Our results suggest that common risk factors in the SP4 gene are associated with schizophrenia, although not with MDD, in the Han Chinese population.

DNA repair variants and breast cancer risk

A functional DNA repair system has been identified as important in the prevention of tumour development. Previous studies have hypothesized that common polymorphisms in DNA repair genes could play a role in breast cancer risk and also identified the potential for interactions between these polymorphisms and established breast cancer risk factors such as physical activity. Associations with breast cancer risk for 99 single nucleotide polymorphisms (SNPs) from genes in ten DNA repair pathways were examined in a case-control study including both Europeans (644 cases, 809 controls) and East Asians (299 cases, 160 controls). Odds ratios in both additive and dominant genetic models were calculated separately for participants of European and East Asian ancestry using multivariate logistic regression. The impact of multiple comparisons was assessed by correcting for the false discovery rate within each DNA repair pathway. Interactions between several breast cancer risk factors and DNA repair SNPs were also evaluated. One SNP (rs3213282) in the gene XRCC1 was associated with an increased risk of breast cancer in the dominant model of inheritance following adjustment for the false discovery rate (P < 0.05), although no associations were observed for other DNA repair SNPs. Interactions of six SNPs in multiple DNA repair pathways with physical activity were evident prior to correction for FDR, following which there was support for only one of the interaction terms (P < 0.05). No consistent associations between variants in DNA repair genes and breast cancer risk or their modification by breast cancer risk factors were observed.