Nonparametric evaluation of quantitative traits in population-based association studies when the genetic model is unknown

A novel application of data-consistent inversion to overcome spurious inference in genome-wide association studies

The genome-wide association studies (GWAS) typically use linear or logistic regression models to identify associations between phenotypes (traits) and genotypes (genetic variants) of interest. However, the use of regression with the additive assumption has potential limitations. First, the normality assumption of residuals is the one that is rarely seen in practice, and deviation from normality increases the Type-I error rate. Second, building a model based on such an assumption ignores genetic structures, like, dominant, recessive, and protective-risk cases. Ignoring genetic variants may result in spurious conclusions about the associations between a variant and a trait. We propose an assumption-free model built upon data-consistent inversion (DCI), which is a recently developed measure-theoretic framework utilized for uncertainty quantification. This proposed DCI-derived model builds a nonparametric distribution on model inputs that propagates to the distribution of observed data without the required normality assumption of residuals in the regression model. This characteristic enables the proposed DCI-derived model to cover all genetic variants without emphasizing on additivity of the classic-GWAS model. Simulations and a replication GWAS with data from the COPDGene demonstrate the ability of this model to control the Type-I error rate at least as well as the classic-GWAS (additive linear model) approach while having similar or greater power to discover variants in different genetic modes of transmission.

Enterovirus Neutralizing Antibodies, Monocyte Toll Like Receptors Expression and Interleukin Profiles Are Similar Between Non-affected and Affected Siblings From Long-Term Discordant Type 1 Diabetes Multiplex-Sib Families: The Importance of HLA Background

Enteroviruses are main candidates among environmental agents in the development of type 1 diabetes (T1D). However, the relationship between virus and the immune system response during T1D pathogenesis is heterogeneous. This is an interesting paradigm and the search for answers would help to highlight the role of viral infection in the etiology of T1D. The current data is a cross-sectional study of affected and non-affected siblings from T1D multiplex-sib families to analyze associations among T1D, genetic, islet autoantibodies and markers of innate immunity. We evaluated the prevalence of anti-virus antibodies (Coxsackie B and Echo) and its relationships with human leukocyte antigen (HLA) class II alleles, TLR expression (monocytes), serum cytokine profile and islet β cell autoantibodies in 51 individuals (40 T1D and 11 non-affected siblings) from 20 T1D multiplex-sib families and 54 healthy control subjects. The viral antibody profiles were similar among all groups, except for antibodies against CVB2, which were more prevalent in the non-affected siblings. TLR4 expression was higher in the T1D multiplex-sib family's members than in the control subjects. TLR4 expression showed a positive correlation with CBV2 antibody prevalence (rS: 0.45; P = 0.03), CXCL8 (rS: 0.65, P = 0.002) and TNF-α (rS: 0.5, P = 0.01) serum levels in both groups of T1D multiplex-sib family. Furthermore, within these families, there was a positive correlation between HLA class II alleles associated with high risk for T1D and insulinoma-associated protein 2 autoantibody (IA-2A) positivity (odds ratio: 38.8; P = 0.021). However, the HLA protective haplotypes against T1D prevalence was higher in the non-affected than the affected siblings. This study shows that although the prevalence of viral infection is similar among healthy individuals and members from the T1D multiplex-sib families, the innate immune response is higher in the affected and in the non-affected siblings from these families than in the healthy controls. However, autoimmunity against β-islet cells and an absence of protective HLA alleles were only observed in the T1D multiplex-sib members with clinical disease, supporting the importance of the genetic background in the development of T1D and heterogeneity of the interaction between environmental factors and disease pathogenesis despite the high genetic diversity of the Brazilian population.

General regression model: A "model-free" association test for quantitative traits allowing to test for the underlying genetic model

Most genome-wide association studies used genetic-model-based tests assuming an additive mode of inheritance, leading to underpowered association tests in case of departure from additivity. The general regression model (GRM) association test proposed by Fisher and Wilson in 1980 makes no assumption on the genetic model. Interestingly, it also allows formal testing of the underlying genetic model. We conducted a simulation study of quantitative traits to compare the power of the GRM test to the classical linear regression tests, the maximum of the three statistics (MAX), and the allele-based (allelic) tests. Simulations were performed on two samples sizes, using a large panel of genetic models, varying genetic models, minor allele frequencies, and the percentage of explained variance. In case of departure from additivity, the GRM was more powerful than the additive regression tests (power gain reaching 80%) and had similar power when the true model is additive. GRM was also as or more powerful than the MAX or allelic tests. The true simulated model was mostly retained by the GRM test. Application of GRM to HbA1c illustrates its gain in power. To conclude, GRM increases power to detect association for quantitative traits, allows determining the genetic model and is easily applicable.

Homology-guided mutational analysis reveals the functional requirements for antinociceptive specificity of collapsin response mediator protein 2-derived peptides

BACKGROUND AND PURPOSE

N-type voltage-gated calcium (Ca_v 2.2) channels are critical determinants of increased neuronal excitability and neurotransmission accompanying persistent neuropathic pain. Although Ca_v 2.2 channel antagonists are recommended as first-line treatment for neuropathic pain, calcium-current blocking gabapentinoids inadequately alleviate chronic pain symptoms and often exhibit numerous side effects. Collapsin response mediator protein 2 (CRMP2) targets Ca_v 2.2 channels to the sensory neuron membrane and allosterically modulates their function. A 15-amino-acid peptide (CBD3), derived from CRMP2, disrupts the functional protein-protein interaction between CRMP2 and Ca_v 2.2 channels to inhibit calcium influx, transmitter release and acute, inflammatory and neuropathic pain. Here, we have mapped the minimal domain of CBD3 necessary for its antinociceptive potential.

EXPERIMENTAL APPROACH

Truncated as well as homology-guided mutant versions of CBD3 were generated and assessed using depolarization-evoked calcium influx in rat dorsal root ganglion neurons, binding between CRMP2 and Ca_v 2.2 channels, whole-cell voltage clamp electrophysiology and behavioural effects in two models of experimental pain: post-surgical pain and HIV-induced sensory neuropathy induced by the viral glycoprotein 120.

KEY RESULTS

The first six amino acids within CBD3 accounted for all in vitro activity and antinociception. Spinal administration of a prototypical peptide (TAT-CBD3-L5M) reversed pain behaviours. Homology-guided mutational analyses of these six amino acids identified at least two residues, Ala1 and Arg4, as being critical for antinociception in two pain models.

CONCLUSIONS AND IMPLICATIONS

These results identify an antinociceptive scaffold core in CBD3 that can be used for development of low MW mimetics of CBD3.

LINKED ARTICLES

This article is part of a themed section on Recent Advances in Targeting Ion Channels to Treat Chronic Pain. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v175.12/issuetoc.

A GWAS in uveal melanoma identifies risk polymorphisms in the CLPTM1L locus

Uveal melanoma, a rare malignant tumor of the eye, is predominantly observed in populations of European ancestry. A genome-wide association study of 259 uveal melanoma patients compared to 401 controls all of European ancestry revealed a candidate locus at chromosome 5p15.33 (region rs421284: OR = 1.7, CI 1.43-2.05). This locus was replicated in an independent set of 276 cases and 184 controls. In addition, risk variants from this region were positively associated with higher expression of CLPTM1L. In conclusion, the CLPTM1L region contains risk alleles for uveal melanoma susceptibility, suggesting that CLPTM1L could play a role in uveal melanoma oncogenesis.

A common risk variant in CACNA1C supports a sex-dependent effect on longitudinal functioning and functional recovery from episodes of schizophrenia-spectrum but not bipolar disorder

Sex is a powerful modulator of disease susceptibility, course and outcome. The gene CACNA1C is among the best replicated vulnerability genes of bipolar disorder and schizophrenia. The aim of the present study was to investigate whether sex and a variant in CACNA1C (rs10774035 as a proxy for the well-acknowledged risk variant rs1006737) influence psychosocial adaptation in a large German patient sample with schizophrenia-spectrum (n=297) and bipolar (n=516) disorders. We analyzed Global Assessment of Functioning (GAF) scores, retrospectively collected for different time points during disease course. We investigated whether CACNA1C sex-dependently modulates longitudinal GAF scores and recovery from episodes of psychiatric disturbance in the above mentioned disorders. Psychosocial recovery was measured as difference score between the current GAF score (assessing the last remission) and the worst GAF score ever during an illness episode. Covariate- adjusted association analyses revealed a sex × rs10774035 genotype interaction on longitudinal GAF and recovery from illness episodes only in schizophrenia-spectrum but not in bipolar disorders. In schizophrenia-spectrum affected males, rs10774035 minor allele (T) carriers had higher GAF scores at three time points (premorbid, worst ever, current). In contrast, females carrying rs10774035 minor alleles had impaired recovery from schizophrenia-spectrum episodes. These results encourage further investigations of gene × sex interactions and longitudinal quantitative phenotypes to unravel the rich variety of behavioral consequences of genetic individuality.

Nonparametric Risk and Nonparametric Odds in Quantitative Genetic Association Studies

The coefficient in a linear regression model is commonly employed to evaluate the genetic effect of a single nucleotide polymorphism associated with a quantitative trait under the assumption that the trait value follows a normal distribution or is appropriately normally distributed after a certain transformation. When this assumption is violated, the distribution-free tests are preferred. In this work, we propose the nonparametric risk (NR) and nonparametric odds (NO), obtain the asymptotic normal distribution of estimated NR and then construct the confidence intervals. We also define the genetic models using NR, construct the test statistic under a given genetic model and a robust test, which are free of the genetic uncertainty. Simulation studies show that the proposed confidence intervals have satisfactory cover probabilities and the proposed test can control the type I error rates and is more powerful than the exiting ones under most of the considered scenarios. Application to gene of PTPN22 and genomic region of 6p21.33 from the Genetic Analysis Workshop 16 for association with the anticyclic citrullinated protein antibody further show their performances.

Suppression of Dwarf and irregular xylem Phenotypes Generates Low-Acetylated Biomass Lines in Arabidopsis

eskimo1-5 (esk1-5) is a dwarf Arabidopsis (Arabidopsis thaliana) mutant that has a constitutive drought syndrome and collapsed xylem vessels, along with low acetylation levels in xylan and mannan. ESK1 has xylan O-acetyltransferase activity in vitro. We used a suppressor strategy on esk1-5 to screen for variants with wild-type growth and low acetylation levels, a favorable combination for ethanol production. We found a recessive mutation in the KAKTUS (KAK) gene that suppressed dwarfism and the collapsed xylem character, the cause of decreased hydraulic conductivity in the esk1-5 mutant. Backcrosses between esk1-5 and two independent knockout kak mutants confirmed suppression of the esk1-5 effect. kak single mutants showed larger stem diameters than the wild type. The KAK promoter fused with a reporter gene showed activity in the vascular cambium, phloem, and primary xylem in the stem and hypocotyl. However, suppression of the collapsed xylem phenotype in esk1 kak double mutants was not associated with the recovery of cell wall O-acetylation or any major cell wall modifications. Therefore, our results indicate that, in addition to its described activity as a repressor of endoreduplication, KAK may play a role in vascular development. Furthermore, orthologous esk1 kak double mutants may hold promise for ethanol production in crop plants.

(S)-Lacosamide Binding to Collapsin Response Mediator Protein 2 (CRMP2) Regulates CaV2.2 Activity by Subverting Its Phosphorylation by Cdk5

The neuronal circuit remodels during development as well as in human neuropathologies such as epilepsy. Neurite outgrowth is an obligatory step in these events. We recently reported that alterations in the phosphorylation state of an axon specification/guidance protein, the collapsin response mediator protein 2 (CRMP2), play a major role in the activity-dependent regulation of neurite outgrowth. We also identified (S)-LCM, an inactive stereoisomer of the clinically used antiepileptic drug (R)-LCM (Vimpat®), as a novel tool for preferentially targeting CRMP2-mediated neurite outgrowth. Here, we investigated the mechanism by which (S)-LCM affects CRMP2 phosphorylation by two key kinases, cyclin-dependent kinase 5 (Cdk5) and glycogen synthase kinase 3β (GSK-3β). (S)-LCM application to embryonic cortical neurons resulted in reduced levels of Cdk5- and GSK-3β-phosphorylated CRMP2. Mechanistically, (S)-LCM increased CRMP2 binding to both Cdk5- and GSK-3β without affecting binding of CRMP2 to its canonical partner tubulin. Saturation transfer difference nuclear magnetic resonance (STD NMR) and differential scanning fluorimetry (DSF) experiments demonstrated direct binding of (S)-LCM to CRMP2. Using an in vitro luminescent kinase assay, we observed that (S)-LCM specifically inhibited Cdk5-mediated phosphorylation of CRMP2. Cross-linking experiments and analytical ultracentrifugation showed no effect of (S)-LCM on the oligomerization state of CRMP2. The increased association between Cdk5-phosphorylated CRMP2 and CaV2.2 was reduced by (S)-LCM in vitro and in vivo. This reduction translated into a decrease of calcium influx via CaV2.2 in (S)-LCM-treated neurons compared to controls. (S)-LCM, to our knowledge, is the first molecule described to directly inhibit CRMP2 phosphorylation and may be useful for delineating CRMP2-facilitated functions.

Analysis of genetic association and epistasis interactions between circadian clock genes and symptom dimensions of bipolar affective disorder

Bipolar affective disorder (BD) is a severe psychiatric disorder characterized by periodic changes in mood from depression to mania. Disruptions of biological rhythms increase risk of mood disorders. Because clinical representation of disease is heterogeneous, homogenous sets of patients are suggested to use in the association analyses. In our study, we aimed to apply previously computed structure of bipolar disorder symptom dimension for analyses of genetic association. We based quantitative trait on: main depression, sleep disturbances, appetite disturbances, excitement and psychotic dimensions consisted of OPCRIT checklist items. We genotyped 42 polymorphisms from circadian clock genes: PER3, ARNTL, CLOCK and TIMELSSS from 511 patients BD (n = 292 women and n = 219 men). As quantitative trait we used clinical dimensions, described above. Genetic associations between alleles and quantitative trait were performed using applied regression models applied in PLINK. In addition, we used the Kruskal-Wallis test to look for associations between genotypes and quantitative trait. During second stage of our analyses, we used multidimensional scaling (multifactor dimensionality reduction) for quantitative trait to compute pairwise epistatic interactions between circadian gene variants. We found association between ARNTL variant rs11022778 main depression (p = 0.00047) and appetite disturbances (p = 0.004). In epistatic interaction analyses, we observed two locus interactions between sleep disturbances (p = 0.007; rs11824092 of ARNTL and rs11932595 of CLOCK) as well as interactions of subdimension in main depression and ARNTL variants (p = 0.0011; rs3789327, rs10766075) and appetite disturbances in depression and ARNTL polymorphism (p = 7 × 10(-4); rs11022778, rs156243).

System for stable β-estradiol-inducible gene expression in the moss Physcomitrella patens

Inducible transgene expression provides a useful tool to analyze gene function. The moss Physcomitrellapatens is a model basal land plant with well-developed research tools, including a high efficiency of gene targeting and substantial genomics resources. However, current systems for controlled transgene expression remain limited. Here we report the development of an estrogen receptor mediated inducible gene expression system, based on the system used in flowering plants. After identifying the appropriate promoters to drive the chimeric transducer, we succeeded in inducing transcription over 1,000-fold after 24 h incubation with β-estradiol. The P. patens system was also effective for high-level long-term induction of gene expression; transcript levels of the activated gene were maintained for at least seven days on medium containing β-estradiol. We also established two potentially neutral targeting sites and a set of vectors for reproducible expression of two transgenes. This β-estradiol-dependent system will be useful to test genes individually or in combination, allowing stable, inducible transgenic expression in P. patens.

The psychiatric vulnerability gene CACNA1C and its sex-specific relationship with personality traits, resilience factors and depressive symptoms in the general population

Genome-wide association studies have reported an association between the A-allele of rs1006737 within CACNA1C and affective disorders and schizophrenia. The aim of the present study was to investigate the relationship between rs1006737 and established and potential endophenotypes for these disorders in a population-based cohort of 3793 subjects, using an analytical method designed to assess a previously reported sex-specific effect of CACNA1C. The investigated endophenotypes included personality traits and resilience factors. At 10-year follow-up, subjects were screened for depressive symptoms. All subjects were genotyped for rs1006737. The direction of the effect and mode of inheritance of rs1006737 differed between the sexes. In men, the A-allele was associated with higher emotional lability and lower resilience, that is, lower sense of coherence (P=0.021), lower perceived social support (P=0.018), lower dispositional optimism (P=0.032) and more depressive symptoms at follow-up (P=0.007). In women, the A-allele was associated with lower emotional lability and stronger resilience, that is, higher sense of coherence (P=0.00028), higher perceived social support (P=0.010), lower neuroticism (P=0.022) and fewer depressive symptoms at follow-up (P=0.035). After conservative Bonferroni correction for 32 tests, results only remained significant for sense of coherence in women (P=0.009). These results suggest that CACNA1C is involved in the genetic architecture of endophenotypes for affective disorders and schizophrenia, and that it shows a distinct sex-specific effect. Comprehensive phenotype characterization in case-control samples and the general population, as well as an adequate modeling of sex-specific genetic effects, may be warranted to elucidate the pathogenetic mechanisms conferred by robustly identified susceptibility genes.