Effect of training-sample size and classification difficulty on the accuracy of genomic predictors.

Introduction: As part of the MicroArray Quality Control (MAQC)-II project, this analysis examines how the choice of univariate feature-selection methods and classification algorithms may influence the performance of genomic predictors under varying degrees of prediction difficulty represented by three clinically relevant endpoints. Methods: We used gene-expression data from 230 breast cancers (grouped into training and independent validation sets), and we examined 40 predictors (five univariate feature-selection methods combined with eight different classifiers) for each of the three endpoints. Their classification performance was estimated on the training set by using two different resampling methods and compared with the accuracy observed in the independent validation set. Results: A ranking of the three classification problems was obtained, and the performance of 120 models was estimated and assessed on an independent validation set. The bootstrapping estimates were closer to the validation performance than were the cross-validation estimates. The required sample size for each endpoint was estimated, and both gene-level and pathway-level analyses were performed on the obtained models. Conclusions: We showed that genomic predictor accuracy is determined largely by an interplay between sample size and classification difficulty. Variations on univariate feature-selection methods and choice of classification algorithm have only a modest impact on predictor performance, and several statistically equally good predictors can be developed for any given classification problem.

A mega-analysis of genome-wide association studies for major depressive disorder.

Prior genome-wide association studies (GWAS) of major depressive disorder (MDD) have met with limited success. We sought to increase statistical power to detect disease loci by conducting a GWAS mega-analysis for MDD. In the MDD discovery phase, we analyzed more than 1.2 million autosomal and X chromosome single-nucleotide polymorphisms (SNPs) in 18 759 independent and unrelated subjects of recent European ancestry (9240 MDD cases and 9519 controls). In the MDD replication phase, we evaluated 554 SNPs in independent samples (6783 MDD cases and 50 695 controls). We also conducted a cross-disorder meta-analysis using 819 autosomal SNPs with P<0.0001 for either MDD or the Psychiatric GWAS Consortium bipolar disorder (BIP) mega-analysis (9238 MDD cases/8039 controls and 6998 BIP cases/7775 controls). No SNPs achieved genome-wide significance in the MDD discovery phase, the MDD replication phase or in pre-planned secondary analyses (by sex, recurrent MDD, recurrent early-onset MDD, age of onset, pre-pubertal onset MDD or typical-like MDD from a latent class analyses of the MDD criteria). In the MDD-bipolar cross-disorder analysis, 15 SNPs exceeded genome-wide significance (P<5 × 10(-8)), and all were in a 248 kb interval of high LD on 3p21.1 (chr3:52 425 083-53 822 102, minimum P=5.9 × 10(-9) at rs2535629). Although this is the largest genome-wide analysis of MDD yet conducted, its high prevalence means that the sample is still underpowered to detect genetic effects typical for complex traits. Therefore, we were unable to identify robust and replicable findings. We discuss what this means for genetic research for MDD. The 3p21.1 MDD-BIP finding should be interpreted with caution as the most significant SNP did not replicate in MDD samples, and genotyping in independent samples will be needed to resolve its status.

Genetic variants in novel pathways influence blood pressure and cardiovascular disease risk.

Blood pressure is a heritable trait influenced by several biological pathways and responsive to environmental stimuli. Over one billion people worldwide have hypertension (≥140 mm Hg systolic blood pressure or  ≥90 mm Hg diastolic blood pressure). Even small increments in blood pressure are associated with an increased risk of cardiovascular events. This genome-wide association study of systolic and diastolic blood pressure, which used a multi-stage design in 200,000 individuals of European descent, identified sixteen novel loci: six of these loci contain genes previously known or suspected to regulate blood pressure (GUCY1A3-GUCY1B3, NPR3-C5orf23, ADM, FURIN-FES, GOSR2, GNAS-EDN3); the other ten provide new clues to blood pressure physiology. A genetic risk score based on 29 genome-wide significant variants was associated with hypertension, left ventricular wall thickness, stroke and coronary artery disease, but not kidney disease or kidney function. We also observed associations with blood pressure in East Asian, South Asian and African ancestry individuals. Our findings provide new insights into the genetics and biology of blood pressure, and suggest potential novel therapeutic pathways for cardiovascular disease prevention.

Defining the role of common variation in the genomic and biological architecture of adult human height.

Using genome-wide data from 253,288 individuals, we identified 697 variants at genome-wide significance that together explained one-fifth of the heritability for adult height. By testing different numbers of variants in independent studies, we show that the most strongly associated ∼2,000, ∼3,700 and ∼9,500 SNPs explained ∼21%, ∼24% and ∼29% of phenotypic variance. Furthermore, all common variants together captured 60% of heritability. The 697 variants clustered in 423 loci were enriched for genes, pathways and tissue types known to be involved in growth and together implicated genes and pathways not highlighted in earlier efforts, such as signaling by fibroblast growth factors, WNT/β-catenin and chondroitin sulfate-related genes. We identified several genes and pathways not previously connected with human skeletal growth, including mTOR, osteoglycin and binding of hyaluronic acid. Our results indicate a genetic architecture for human height that is characterized by a very large but finite number (thousands) of causal variants.

The MicroArray Quality Control (MAQC)-II study of common practices for the development and validation of microarray-based predictive models.

Gene expression data from microarrays are being applied to predict preclinical and clinical endpoints, but the reliability of these predictions has not been established. In the MAQC-II project, 36 independent teams analyzed six microarray data sets to generate predictive models for classifying a sample with respect to one of 13 endpoints indicative of lung or liver toxicity in rodents, or of breast cancer, multiple myeloma or neuroblastoma in humans. In total, >30,000 models were built using many combinations of analytical methods. The teams generated predictive models without knowing the biological meaning of some of the endpoints and, to mimic clinical reality, tested the models on data that had not been used for training. We found that model performance depended largely on the endpoint and team proficiency and that different approaches generated models of similar performance. The conclusions and recommendations from MAQC-II should be useful for regulatory agencies, study committees and independent investigators that evaluate methods for global gene expression analysis.

A role of peripheral myelin protein 2 in lipid homeostasis of myelinating Schwann cells.

Peripheral myelin protein 2 (Pmp2, P2 or Fabp8), a member of the fatty acid binding protein family, was originally described together with myelin basic protein (Mbp or P1) and myelin protein zero (Mpz or P0) as one of the most abundant myelin proteins in the peripheral nervous system (PNS). Although Pmp2 is predominantly expressed in myelinated Schwann cells, its role in glia is currently unknown. To study its function in PNS biology, we have generated a complete Pmp2 knockout mouse (Pmp2(-/-) ). Comprehensive characterization of Pmp2(-/-) mice revealed a temporary reduction in their motor nerve conduction velocity (MNCV). While this change was not accompanied by any defects in general myelin structure, we detected transitory alterations in the myelin lipid profile of Pmp2(-/-) mice. It was previously proposed that Pmp2 and Mbp have comparable functions in the PNS suggesting that the presence of Mbp can partially mask the Pmp2(-/-) phenotype. Indeed, we found that Mbp lacking Shi(-/-) mice, similar to Pmp2(-/-) animals, have preserved myelin structure and reduced MNCV, but this phenotype was not aggravated in Pmp2(-/-) /Shi(-/-) mutants indicating that Pmp2 and Mbp do not substitute each other's functions in the PNS. These data, together with our observation that Pmp2 binds and transports fatty acids to membranes, uncover a role for Pmp2 in lipid homeostasis of myelinating Schwann cells.

Guidelines for the use and interpretation of assays for monitoring autophagy.

In 2008 we published the first set of guidelines for standardizing research in autophagy. Since then, research on this topic has continued to accelerate, and many new scientists have entered the field. Our knowledge base and relevant new technologies have also been expanding. Accordingly, it is important to update these guidelines for monitoring autophagy in different organisms. Various reviews have described the range of assays that have been used for this purpose. Nevertheless, there continues to be confusion regarding acceptable methods to measure autophagy, especially in multicellular eukaryotes. A key point that needs to be emphasized is that there is a difference between measurements that monitor the numbers or volume of autophagic elements (e.g., autophagosomes or autolysosomes) at any stage of the autophagic process vs. those that measure flux through the autophagy pathway (i.e., the complete process); thus, a block in macroautophagy that results in autophagosome accumulation needs to be differentiated from stimuli that result in increased autophagic activity, defined as increased autophagy induction coupled with increased delivery to, and degradation within, lysosomes (in most higher eukaryotes and some protists such as Dictyostelium) or the vacuole (in plants and fungi). In other words, it is especially important that investigators new to the field understand that the appearance of more autophagosomes does not necessarily equate with more autophagy. In fact, in many cases, autophagosomes accumulate because of a block in trafficking to lysosomes without a concomitant change in autophagosome biogenesis, whereas an increase in autolysosomes may reflect a reduction in degradative activity. Here, we present a set of guidelines for the selection and interpretation of methods for use by investigators who aim to examine macroautophagy and related processes, as well as for reviewers who need to provide realistic and reasonable critiques of papers that are focused on these processes. These guidelines are not meant to be a formulaic set of rules, because the appropriate assays depend in part on the question being asked and the system being used. In addition, we emphasize that no individual assay is guaranteed to be the most appropriate one in every situation, and we strongly recommend the use of multiple assays to monitor autophagy. In these guidelines, we consider these various methods of assessing autophagy and what information can, or cannot, be obtained from them. Finally, by discussing the merits and limits of particular autophagy assays, we hope to encourage technical innovation in the field.

Alum interaction with dendritic cell membrane lipids is essential for its adjuvanticity.

As an approved vaccine adjuvant for use in humans, alum has vast health implications, but, as it is a crystal, questions remain regarding its mechanism. Furthermore, little is known about the target cells, receptors, and signaling pathways engaged by alum. Here we report that, independent of inflammasome and membrane proteins, alum binds dendritic cell (DC) plasma membrane lipids with substantial force. Subsequent lipid sorting activates an abortive phagocytic response that leads to antigen uptake. Such activated DCs, without further association with alum, show high affinity and stable binding with CD4(+) T cells via the adhesion molecules intercellular adhesion molecule-1 (ICAM-1) and lymphocyte function-associated antigen-1 (LFA-1). We propose that alum triggers DC responses by altering membrane lipid structures. This study therefore suggests an unexpected mechanism for how this crystalline structure interacts with the immune system and how the DC plasma membrane may behave as a general sensor for solid structures.

Genome-wide association analyses identify 18 new loci associated with serum urate concentrations.

Elevated serum urate concentrations can cause gout, a prevalent and painful inflammatory arthritis. By combining data from >140,000 individuals of European ancestry within the Global Urate Genetics Consortium (GUGC), we identified and replicated 28 genome-wide significant loci in association with serum urate concentrations (18 new regions in or near TRIM46, INHBB, SFMBT1, TMEM171, VEGFA, BAZ1B, PRKAG2, STC1, HNF4G, A1CF, ATXN2, UBE2Q2, IGF1R, NFAT5, MAF, HLF, ACVR1B-ACVRL1 and B3GNT4). Associations for many of the loci were of similar magnitude in individuals of non-European ancestry. We further characterized these loci for associations with gout, transcript expression and the fractional excretion of urate. Network analyses implicate the inhibins-activins signaling pathways and glucose metabolism in systemic urate control. New candidate genes for serum urate concentration highlight the importance of metabolic control of urate production and excretion, which may have implications for the treatment and prevention of gout.

Molecular determinants of vesicular storage and release of the gliotransmitter D-serine from cortical astrocytes

Neuron-astrocyte reciprocal communication at synapses has emerged as a novel signalling pathway in brain function. Astrocytes sense the level of synaptic activity and, in turn, influence its efficacy through the regulated release of 'gliotransmitters' such as glutamate, ATP or D-serine. A calcium-dependent exocytosis is proposed to drive the release of gliotransmitters but its existence is still debated. Over the last years, we have been studying the molecular determinants governing D-serine release from glia using different approaches. Using a novel bioassay for D-serine, we have been able to show that D-serine release occurs mainly through a calcium- and SNARE proteindependent mechanism just supporting the idea that this amino acid is released by exocytosis from glia. We next have pursued our exploration by confocal imaging and tracking of the exocytotic routes for Dserine- mediated gliotransmission and have shown that D-serine releasable pools are confined to synaptobrevin2/cellubrevin-bearing vesicles. To shed light onto the mechanisms controlling the storage and the release of gliotransmitters and namely D-serine, we have developed a new method for the immunoisolation of synaptobrevin 2- positive vesicles from rat cortical astrocytes in culture while preserving their content in gliotransmitters. The purified organelles are clear round shape vesicles of excellent purity with homogeneous size (40 nm) as judged by electron microscopy. Immunoblotting analysis revealed that isolated vesicles contain most of the major proteins already described for neuron-derived vesicles like synaptic vesicle protein 2 (SV2) and the proton pump H?-ATPase. In addition, we have analyzed the content for various amino acids of these vesicles by means of chiral capillary electrophoresis coupled to laser-induced fluorescence detection. The purified vesicles contain large amount of D-serine. We also detect peaks corresponding to unidentified compounds that may correspond to others amino acids. Postembedding immunogold labelling of the rat neocortex further revealed the expression of D-serine in astrocytes processes contacting excitatory synapses. Finally, we have examined the uptake properties for Dserine and glutamate inside the isolated glial vesicles. Our results provide significant support for the existence of an uptake system for D-serine in secretory glial vesicles and for the storage of chemical substances like D-serine and glutamate. 11th International Congress on Amino Acids, Peptides and Proteins 763 123

Genome-wide association study of major recurrent depression in the U.K. population.

OBJECTIVE: Studies of major depression in twins and families have shown moderate to high heritability, but extensive molecular studies have failed to identify susceptibility genes convincingly. To detect genetic variants contributing to major depression, the authors performed a genome-wide association study using 1,636 cases of depression ascertained in the U.K. and 1,594 comparison subjects screened negative for psychiatric disorders. METHOD: Cases were collected from 1) a case-control study of recurrent depression (the Depression Case Control [DeCC] study; N=1346), 2) an affected sibling pair linkage study of recurrent depression (probands from the Depression Network [DeNT] study; N=332), and 3) a pharmacogenetic study (the Genome-Based Therapeutic Drugs for Depression [GENDEP] study; N=88). Depression cases and comparison subjects were genotyped at Centre National de Génotypage on the Illumina Human610-Quad BeadChip. After applying stringent quality control criteria for missing genotypes, departure from Hardy-Weinberg equilibrium, and low minor allele frequency, the authors tested for association to depression using logistic regression, correcting for population ancestry. RESULTS: Single nucleotide polymorphisms (SNPs) in BICC1 achieved suggestive evidence for association, which strengthened after imputation of ungenotyped markers, and in analysis of female depression cases. A meta-analysis of U.K. data with previously published results from studies in Munich and Lausanne showed some evidence for association near neuroligin 1 (NLGN1) on chromosome 3, but did not support findings at BICC1. CONCLUSIONS: This study identifies several signals for association worthy of further investigation but, as in previous genome-wide studies, suggests that individual gene contributions to depression are likely to have only minor effects, and very large pooled analyses will be required to identify them.

Gene-Age Interactions in Blood Pressure Regulation: A Large-Scale Investigation with the CHARGE, Global BPgen, and ICBP Consortia.

Although age-dependent effects on blood pressure (BP) have been reported, they have not been systematically investigated in large-scale genome-wide association studies (GWASs). We leveraged the infrastructure of three well-established consortia (CHARGE, GBPgen, and ICBP) and a nonstandard approach (age stratification and metaregression) to conduct a genome-wide search of common variants with age-dependent effects on systolic (SBP), diastolic (DBP), mean arterial (MAP), and pulse (PP) pressure. In a two-staged design using 99,241 individuals of European ancestry, we identified 20 genome-wide significant (p ≤ 5 × 10(-8)) loci by using joint tests of the SNP main effect and SNP-age interaction. Nine of the significant loci demonstrated nominal evidence of age-dependent effects on BP by tests of the interactions alone. Index SNPs in the EHBP1L1 (DBP and MAP), CASZ1 (SBP and MAP), and GOSR2 (PP) loci exhibited the largest age interactions, with opposite directions of effect in the young versus the old. The changes in the genetic effects over time were small but nonnegligible (up to 1.58 mm Hg over 60 years). The EHBP1L1 locus was discovered through gene-age interactions only in whites but had DBP main effects replicated (p = 8.3 × 10(-4)) in 8,682 Asians from Singapore, indicating potential interethnic heterogeneity. A secondary analysis revealed 22 loci with evidence of age-specific effects (e.g., only in 20 to 29-year-olds). Age can be used to select samples with larger genetic effect sizes and more homogenous phenotypes, which may increase statistical power. Age-dependent effects identified through novel statistical approaches can provide insight into the biology and temporal regulation underlying BP associations.

Edoxaban versus warfarin for the treatment of symptomatic venous thromboembolism.

BACKGROUND: Whether the oral factor Xa inhibitor edoxaban can be an alternative to warfarin in patients with venous thromboembolism is unclear. METHODS: In a randomized, double-blind, noninferiority study, we randomly assigned patients with acute venous thromboembolism, who had initially received heparin, to receive edoxaban at a dose of 60 mg once daily, or 30 mg once daily (e.g., in the case of patients with creatinine clearance of 30 to 50 ml per minute or a body weight below 60 kg), or to receive warfarin. Patients received the study drug for 3 to 12 months. The primary efficacy outcome was recurrent symptomatic venous thromboembolism. The principal safety outcome was major or clinically relevant nonmajor bleeding. RESULTS: A total of 4921 patients presented with deep-vein thrombosis, and 3319 with a pulmonary embolism. Among patients receiving warfarin, the time in the therapeutic range was 63.5%. Edoxaban was noninferior to warfarin with respect to the primary efficacy outcome, which occurred in 130 patients in the edoxaban group (3.2%) and 146 patients in the warfarin group (3.5%) (hazard ratio, 0.89; 95% confidence interval [CI], 0.70 to 1.13; P<0.001 for noninferiority). The safety outcome occurred in 349 patients (8.5%) in the edoxaban group and 423 patients (10.3%) in the warfarin group (hazard ratio, 0.81; 95% CI, 0.71 to 0.94; P=0.004 for superiority). The rates of other adverse events were similar in the two groups. A total of 938 patients with pulmonary embolism had right ventricular dysfunction, as assessed by measurement of N-terminal pro-brain natriuretic peptide levels; the rate of recurrent venous thromboembolism in this subgroup was 3.3% in the edoxaban group and 6.2% in the warfarin group (hazard ratio, 0.52; 95% CI, 0.28 to 0.98). CONCLUSIONS: Edoxaban administered once daily after initial treatment with heparin was noninferior to high-quality standard therapy and caused significantly less bleeding in a broad spectrum of patients with venous thromboembolism, including those with severe pulmonary embolism. (Funded by Daiichi-Sankyo; Hokusai-VTE ClinicalTrials.gov number, NCT00986154.).

Salt-sensitive hypertension and cardiac hypertrophy in mice deficient in the ubiquitin ligase Nedd4-2.

Nedd4-2 has been proposed to play a critical role in regulating epithelial Na+ channel (ENaC) activity. Biochemical and overexpression experiments suggest that Nedd4-2 binds to the PY motifs of ENaC subunits via its WW domains, ubiquitinates them, and decreases their expression on the apical membrane. Phosphorylation of Nedd4-2 (for example by Sgk1) may regulate its binding to ENaC, and thus ENaC ubiquitination. These results suggest that the interaction between Nedd4-2 and ENaC may play a crucial role in Na+ homeostasis and blood pressure (BP) regulation. To test these predictions in vivo, we generated Nedd4-2 null mice. The knockout mice had higher BP on a normal diet and a further increase in BP when on a high-salt diet. The hypertension was probably mediated by ENaC overactivity because 1) Nedd4-2 null mice had higher expression levels of all three ENaC subunits in kidney, but not of other Na+ transporters; 2) the downregulation of ENaC function in colon was impaired; and 3) NaCl-sensitive hypertension was substantially reduced in the presence of amiloride, a specific inhibitor of ENaC. Nedd4-2 null mice on a chronic high-salt diet showed cardiac hypertrophy and markedly depressed cardiac function. Overall, our results demonstrate that in vivo Nedd4-2 is a critical regulator of ENaC activity and BP. The absence of this gene is sufficient to produce salt-sensitive hypertension. This model provides an opportunity to further investigate mechanisms and consequences of this common disorder.