Skip to main content
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian
  • Log in
  • My Cart

Main menu

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home

Advanced Search

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses

New Research In

Physical Sciences

Featured Portals

  • Physics
  • Chemistry
  • Sustainability Science

Articles by Topic

  • Applied Mathematics
  • Applied Physical Sciences
  • Astronomy
  • Computer Sciences
  • Earth, Atmospheric, and Planetary Sciences
  • Engineering
  • Environmental Sciences
  • Mathematics
  • Statistics

Social Sciences

Featured Portals

  • Anthropology
  • Sustainability Science

Articles by Topic

  • Economic Sciences
  • Environmental Sciences
  • Political Sciences
  • Psychological and Cognitive Sciences
  • Social Sciences

Biological Sciences

Featured Portals

  • Sustainability Science

Articles by Topic

  • Agricultural Sciences
  • Anthropology
  • Applied Biological Sciences
  • Biochemistry
  • Biophysics and Computational Biology
  • Cell Biology
  • Developmental Biology
  • Ecology
  • Environmental Sciences
  • Evolution
  • Genetics
  • Immunology and Inflammation
  • Medical Sciences
  • Microbiology
  • Neuroscience
  • Pharmacology
  • Physiology
  • Plant Biology
  • Population Biology
  • Psychological and Cognitive Sciences
  • Sustainability Science
  • Systems Biology
Research Article

Susceptibility of brain atrophy to TRIB3 in Alzheimer’s disease, evidence from functional prioritization in imaging genetics

Marco Lorenzi, Andre Altmann, Boris Gutman, Selina Wray, Charles Arber, Derrek P. Hibar, Neda Jahanshad, Jonathan M. Schott, Daniel C. Alexander, Paul M. Thompson, Sebastien Ourselin, and for the Alzheimer’s Disease Neuroimaging Initiative
PNAS March 20, 2018 115 (12) 3162-3167; first published March 6, 2018; https://doi.org/10.1073/pnas.1706100115
Marco Lorenzi
aTranslational Imaging Group, Centre for Medical Image Computing, University College London, London WC1E 7JE, United Kingdom;
bEpione Research Project, Université Côte d’Azur, BP 93 06 902 Inria Sophia Antipolis, France;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: marco.lorenzi@inria.fr
Andre Altmann
aTranslational Imaging Group, Centre for Medical Image Computing, University College London, London WC1E 7JE, United Kingdom;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Boris Gutman
cImaging Genetics Center, University of Southern California, Marina Del Rey, CA 90292;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Selina Wray
dDepartment of Molecular Neuroscience, University College London Institute of Neurology, London WC1N 1PJ, United Kingdom;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Charles Arber
dDepartment of Molecular Neuroscience, University College London Institute of Neurology, London WC1N 1PJ, United Kingdom;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Derrek P. Hibar
cImaging Genetics Center, University of Southern California, Marina Del Rey, CA 90292;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Neda Jahanshad
cImaging Genetics Center, University of Southern California, Marina Del Rey, CA 90292;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jonathan M. Schott
eDepartment of Neurodegeneration, Dementia Research Centre, Institute of Neurology, London WC1N 3BG, United Kingdom;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Daniel C. Alexander
fCentre of Medical Image Computing, University College London, London WC1E 7JE, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Paul M. Thompson
cImaging Genetics Center, University of Southern California, Marina Del Rey, CA 90292;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Sebastien Ourselin
aTranslational Imaging Group, Centre for Medical Image Computing, University College London, London WC1E 7JE, United Kingdom;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  1. Edited by Marcus E. Raichle, Washington University in St. Louis, St. Louis, MO, and approved January 31, 2018 (received for review April 19, 2017)

This article has a Letter. Please see:

  • Disease status affects the association between rs4813620 and the expression of Alzheimer’s disease susceptibility gene TRIB3 - October 24, 2018
  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

Significance

In this study, we use an experimental imaging–genetics approach for investigating the genetic underpinnings of brain atrophy in Alzheimer’s disease. We successfully combined state-of-the-art imaging–genetics methods and experimental gene expression data to uncover biology in brain atrophy. The experimental paradigm highlighted a significant role of tribbles pseudokinase 3 (TRIB3) in modulating the typical pattern of Alzheimer’s brain pathology. This result corroborates through rigorous data-driven statistical methods evidence emerging from previous studies about the role of TRIB3 in modulating known mechanisms of neurodegeneration, such as neuronal death, cellular homeostasis, and interaction with established genes causing autosomal dominant Alzheimer’s disease: APP and PSEN1. The developed integrated statistical–experimental methodology could serve as a roadmap for investigations in other disorders.

Abstract

The joint modeling of brain imaging information and genetic data is a promising research avenue to highlight the functional role of genes in determining the pathophysiological mechanisms of Alzheimer’s disease (AD). However, since genome-wide association (GWA) studies are essentially limited to the exploration of statistical correlations between genetic variants and phenotype, the validation and interpretation of the findings are usually nontrivial and prone to false positives. To address this issue, in this work, we investigate the functional genetic mechanisms underlying brain atrophy in AD by studying the involvement of candidate variants in known genetic regulatory functions. This approach, here termed functional prioritization, aims at testing the sets of gene variants identified by high-dimensional multivariate statistical modeling with respect to known biological processes to introduce a biology-driven validation scheme. When applied to the Alzheimer’s Disease Neuroimaging Initiative (ADNI) cohort, the functional prioritization allowed for identifying a link between tribbles pseudokinase 3 (TRIB3) and the stereotypical pattern of gray matter loss in AD, which was confirmed in an independent validation sample, and that provides evidence about the relation between this gene and known mechanisms of neurodegeneration.

  • imaging–genetics
  • Alzheimer’s disease
  • brain atrophy
  • TRIB3
  • neuroimaging

Alzheimer’s disease (AD) is a devastating neurodegenerative disorder, and its etiology still remains largely concealed. In anticipation of increasing prevalence of AD and other dementias, there is an urgent need for improving the understanding of the disease processes that underlie neurodegeneration. While the knowledge about the genetic and environmental risks underpinning AD is steadily advancing, how these factors interact to lead to the complex pathophysiology that results in dementia is less understood.

Advances in imaging technologies have led to noninvasive or minimally invasive imaging biomarkers that capture various aspects of the disease process, including amyloid deposition (1), tau pathology (2), functional decline (3), and neuronal loss (4). Combining such imaging information with genetic measurements—so-called imaging–genetics—provides the means for investigating the effect of genetic variation on underlying biological mechanisms (5).

Genome-wide association studies (GWAS) query millions of SNPs individually for their association with either case–control status (6) or disease-specific quantitative phenotypes [e.g., in the case of AD, regional brain volumes (7) or brain amyloid burden (8)]. Mass univariate analysis of genetic data is still the predominant method in virtue of its ease of use and well-established theoretical framework, albeit that it suffers from significant limitations, including the requirement for multiple testing, redundancies introduced by linkage disequilibrium (LD), and the lack of analysis of epistatic effects (e.g., SNP–SNP interactions), which have to be explicitly modeled and searched for exhaustively (9). Moreover, more than one quantitative phenotype can be derived from the available imaging data (e.g., dozens or hundreds of regional brain volumes or hundreds of thousands of voxel-level metrics) (10). This potentially large number of genotype–phenotypes features of interest generally complicates the problem of reliably detecting statistical associations and thus, hampers the identification of disease-relevant genetic markers by purely statistical means.

Limitations of classical mass univariate statistical methods have, in recent years, been overcome by using multivariate approaches to data analysis in the context of neuroscience studies (11) and GWAS (12). Likewise, in imaging–genetics, meaningful genotype–phenotype interactions (13) are captured by simultaneously modeling sets of genetic variants that are jointly associated with a given imaging phenotype (14⇓⇓–17). Multivariate GWAS have the potential to shed light on the complex genotype–phenotype relationship and may thus highlight novel links between brain physiology and molecular and biological functions. However, although these methods have proven their ability to identify meaningful SNP combinations associated with brain imaging features, the interpretation and validation of the statistical findings remain very challenging tasks. These problems relate directly to the understanding of the functional role of sets of genetic variants and to the difficulty of replicating the statistical results in unseen cohorts.

We approach this technical bottleneck by leveraging multivariate approaches to explore high-dimensional datasets and to generate hypotheses, which are subsequently tested in downstream experiments. High-quality databases of matched genotype and gene expression measurements, such as the Genotype Tissue Expression project (GTEx; gtexportal.org/home/) (18) and BRAINEAC (www.braineac.org/) (19), facilitate the quantification of effects of SNPs on gene expression in numerous tissues, including various brain tissues. Typically, these databases are used to detail the effect of a genetic variant at the very end of an analysis pipeline and to garner evidence for molecular mechanisms of the genetic locus. However, functional information in “convenience” databases can also be used at an earlier stage in the analysis to prioritize a few candidate hypotheses with a clear functional mechanism [e.g., expression quantitative trait loci (eQTL)] for the validation phase and thus, limit the multiple testing burden.

In this work, we apply this investigative approach to study the genetic functional mechanisms underlying brain atrophy in AD. The framework is composed of two steps.

  • i) Statistical discovery. Candidate genetic variants are initially identified through data-driven multivariate statistical analysis of the matched imaging and genetics data (Fig. 1). This is achieved by modeling the joint covariation between 1.1 million SNPs and the cortical and subcortical atrophy represented by 327,684 cortical and 27,120 subcortical thickness values of 639 individuals (either healthy older controls or patients with AD) from the Alzheimer’s Disease Neuroimaging Initiative (ADNI) cohort.

  • ii) Functional prioritization. The candidate genetic variants are subsequently screened for functional relevance by querying high-dimensional gene expression databases, such as GTEx.

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

Cross-validation scheme for the assessment of the genetic loci of maximal genotype–phenotype correlation identified by the PLS model. The whole procedure is repeated 1 million times, and the resulting array is further analyzed. (A) PLS is applied in a split-half setting. For each of the two nonoverlapping randomly sampled groups, the PLS components of joint phenotype and genotype variation are independently estimated. (B, Left) Each chromosome is partitioned in bins of 10-kb size, which are labeled one if they contain an SNP associated with the largest PLS weights (top 10% of absolute values) or zero otherwise. To obtain stable estimates of the loci of maximal weights, the resulting binary arrays independently estimated in the two groups are merged (binwise and operation). The same procedure is applied on the mesh-based PLS weights associated with the phenotype component. (C) Steps A and B are repeated across 1 million folds, and the results are subsequently averaged to obtain the confidence maps associated with genetic and phenotype components (Figs. 2 and 3).

The resulting small set of genetic loci, which are shown to modify gene expression, is then validated in an independent sample of 553 individuals from the ADNI diagnosed with mild cognitive impairment (MCI), a proportion of whom progressed to AD.

Compared with previous approaches, our work (i) analyzes the whole genome and whole brain in a hypothesis-free fashion (i.e., without preselecting SNPs or brain regions) and (ii) uses a functional prioritization step to select genetic loci for validation in an independent cohort.

Starting from the initial ∼1.1 million SNPs, the multivariate statistical analysis allowed for the identification of a relatively small number of genetic loci that are statistically associated with the typical pattern of AD brain pathology. The subsequent functional prioritization step ultimately identified a significant role of tribbles pseudokinase 3 (TRIB3), a gene showing important connections to known mechanisms of neurodegenerative diseases. Indeed, although a role for TRIB3 in dementia has not been extensively explored, there are several aspects of TRIB3 function that have relevance to mechanisms related to neuronal death, cellular homeostasis, and interaction with established AD genes, such as APP and PSEN1.

This study ultimately offers an illustration of the potential of effectively combining multivariate statistical modeling in imaging–genetics with recent instruments available from computational biology to lead to insights on the pathophysiology of neurodegeneration.

Results

Model Training and Estimated Components.

Figs. 2 and 3 show the relevant areas of the identified joint genetic and phenotype variation, respectively, for the first three partial least squares (PLS) components through stability selection. The components were very robust (100% reproducible) during the stability selection procedure (SI Materials and Methods). The fourth and fifth components did not present any relevant locations (i.e., all bins have P < 0.95) after stability selection for both the genetic modality and the imaging modality.

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

PLS genotype component: the circular plots show the probability of a given genetic locus being associated with the phenotype components shown in Fig. 3. Outer to inner circles represents components 1–3, respectively. The plots show the probability of a given genetic bin of size 10 kb being relevant in the PLS model (i.e., to contain an SNP that is ranked in the top 10% of the absolute weights of the genotype component). The genes with eQTL close to the important loci (P > 0.95) are listed in the innermost circle depending on their genomic position. The red radial lines are located in correspondence of known AD genes: ABCA7, APOE, APP, BIN1, CASS4, CD2AP, CD33, CELF1, CLU, CR1, DSG2, EPHA1, FERMT2, HLA-DRB5, INPP5D, MAPT, MEF2C, MS4, NME8, PICALM, PSEN1, PSEN2, PTK2B, SLC24A4, SORL1, and ZCWPW1. High-resolution circular plots for each component are provided in Figs. S2–S4.

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

PLS phenotype component: Upper shows the topographical distribution of the PLS weights associated with the cortical and subcortical brain areas. The absolute value of the weights is proportional to the importance of the underlying brain areas. The relevance of the brain areas is quantified in Lower. The colors (red to white) indicate the probability of a brain area being associated with the genotype component shown in Fig. 2 and quantify the probability of each cortical mesh point being relevant in the PLS model (i.e., to be ranked among the top 10% of the absolute weights of the phenotype component). L, left; R, right.

Genetic Components.

The circular Manhattan plot [Circos v0.96 (20)] of Fig. 2 shows the selection frequency for the PLS genotype components and describes the importance of the genetic loci associated with cortical thickness variation for components 1–3. The plot shows the probability of a given genetic bin of size 10 kb being relevant in the PLS model (i.e., to contain an SNP that is ranked in the top 10% of the absolute weights of the genotype component). Spatially contiguous loci generally show similar importance values, which are caused by LD of these regions.

In the genetic components 1–3, a total of 118 bins exceeded the selection frequency threshold (61, 50, and 7 for components 1–3, respectively). From these bins, 402 (196, 181, and 25) influential SNPs were extracted and annotated with 98 genes through the Ensembl Variant Effect Predictor (VEP) for GRCh37 (date accessed: October 17, 2016) (21). The extended Apolipoprotein E (APOE) locus comprising APOE and TOMM40 was selected as the highest scoring region in component 1. A total of 3,956 candidate SNP–gene pairs were considered for the GTEx-based eQTL analysis in six tissues. However, a few genes did not show sufficient expression levels in some tissues, and these combinations were excluded from the analysis, resulting in 1,598 unique SNP–gene tissue tests, of which 104 were significant at the Bonferroni-corrected P value threshold (P = 3.1e−5) (Dataset S1) linking to 14 genes (Dataset S2 and Fig. S5): CAPN9, CRYL1, FAM135B, IL-10 receptor subunit alpha (IL10RA), IP6K3, ITGA1, KIN, LAMC1, LINC00941, LYSMD4, RBPMS2, RP11-181K3.4, TM2 domain-containing 1 (TM2D1), and TRIB3. These genes are listed in the innermost circle of Fig. 2 depending on their genomic position.

The independent validation of those 14 genes in the MCI cohort confirmed TRIB3 (P = 0.0034) (Table 1). Three additional genes were close to nominal significance: TM2D1 (P = 0.053), LAMC1 (P = 0.062), and RP11-181K3.4 (P = 0.053) (Table 1). Of note, the top eQTL SNP for TRIB3 rs4813620 received a P = 0.06175 in stage 1 of a large AD GWAS (6). However, rs62191440, an SNP in strong LD with rs4813620 (D′ = 0.8469; r2 = 0.6559) in the European population (22), received a P value of 0.00601 (Fig. S6) and also constitutes an eQTL for TRIB3 in various tissues in GTEx, including brain tissues cortex and caudate ganglia (Fig. S7). Interestingly, when estimating the PLS components on the subcohort of 279 training individuals positive to amyloid in the cerebrospinal fluid (CSF) (Table 2), we achieved compatible validation results in the independent MCI group. Within this setting, TRIB3 still leads to marginally significant differences (P = 0.0134) between progressing and stable MCI (Table S1).

View this table:
  • View inline
  • View popup
Table 1.

Statistical comparison of the genes scores in training and testing groups (Kruskal–Wallis nonparametric test)

View this table:
  • View inline
  • View popup
Table 2.

Sociodemographic, clinical, and genetic characteristics of the study cohort [mean (SD)]

Morphometric Components.

Fig. 3 shows the PLS phenotype components 1–3 along with the associated selection frequency describing the loci of brain atrophy associated with genetic variation. The first component is mainly associated with the thinning of the cortical mantle and is localized in temporal and posterior cingulate cortices. The relevant areas at the subcortical level are primarily associated with amygdalae and thalami. The second component is mostly associated with the thinning of the subcortical areas (hippocampi and amygdalae) and with the cortical thinning of the temporal areas at the cortical level. The third component is similar to component 2 and describes a subcortical thickness pattern prevalent in hippocampi, amygdalae, and thalami. At the cortical level, the component is associated with the thinning of frontal cortice and with isolated spots located in the parahippocampal gyrus.

Discussion

In this work, we modeled high-dimensional genome-wide SNP data and brain-wide cortical thickness data via joint multivariate statistical modeling and functional prioritization of genes through bioinformatics annotation and a large eQTL database.

Our study ultimately identified a link between TRIB3 and the stereotypical pattern of gray matter loss in AD (cortical thinning in temporal and posterior cingulate regions and subcortical atrophy). TRIB3 is a pseudokinase that acts as a regulator of several signaling pathways. For example, it can interact directly with Akt and inhibit the prosurvival Akt pathway (23). TRIB3 expression is induced during neuronal cell death (24), and recently increased levels of the TRIB3 protein were found in dopaminergic neurons of the substantia nigra pars compacta in patients with Parkinson’s disease (25). TRIB3 expression is stress induced and increases in response to nerve growth factor deprivation, endoplasmatic reticulum stress, and amino acid deprivation (24). Although a role for TRIB3 in dementia has not been extensively explored, there are several aspects of TRIB3 function that have relevance to known mechanisms of neurodegenerative disease. TRIB3 can interact directly with P62 to modulate autophagic flux (26), an important process in maintaining cellular homeostasis that is known to be disrupted in neurodegeneration (27). Knockdown of TRIB3 modulates PSEN1 stability (26), and a yeast two-hybrid screen identified progranulin as a direct interaction partner of TRIB3 (28). Intriguingly, it has recently been shown that TRIB3 induces both apoptosis and autophagy in Aβ-induced neuronal death, and silencing of TRIB3 was strongly neuroprotective (29). These links warrant additional investigation for a functional role of TRIB3 in neuronal death in dementia.

These earlier findings align with our eQTL analysis, where carriers of the minor allele show increased TRIB3 expression (Fig. S5), which potentially lowers the threshold to TRIB3-mediated neuronal cell death. TRIB3 expression was modulated by the identified SNP in various other tissues, including the caudate (Fig. S7), a region affected in Parkinson’s disease and Huntington’s disease. A recent study of Trib3 expression in mice concluded that “Trib3 has a pathophysiological role in diabetes” (30); diabetes itself is a known risk factor for dementia (31), perhaps through shared metabolic processes with AD (32). Interestingly, one of three SNPs (rs1555318) selected in the PLS model and attributed to TRIB3 showed a strong association with type 2 diabetes in stage 1 of a large GWAS (P = 4.4e−4) (Fig. S8) (33). Other GWAS showed links between TRIB3 and information processing speed (P = 1.7e−7) (34) and AD (P = 0.006) (6). An earlier genetic study on AD in Swedish men found an association in TRIB3 as well (P = 0.044) (35), which was replicated in a Canadian cohort (P < 0.001) (36). Lastly, TRIB3 was reported to physically interact with APP (37), and it shares numerous functional annotations for biological processes regarding lipid metabolism with APOE.

The functional prioritization component of the analysis successfully reduced the set of candidate genetic variants for the independent validation; however, this prioritization has a shortcoming: it hypothesizes that identified SNPs alter the expression of a nearby gene. Although this scheme led to the identification of TRIB3 in the cortical thickness phenotype, it did miss a long-established AD risk gene: APOE. SNPs belonging to APOE (rs429358 and rs7412) were selected as the highest-scoring SNPs in component 1. However, none of them were detected as an eQTL, and thus, APOE was excluded from the downstream analysis. Other types of functional prioritizations based on exonic function prediction may have retained APOE and other genes in the pipeline. However, SNP data typically feature only a few nonsynonymous exonic variants, and their high frequency (minor allele frequency > 5%) renders them unlikely to receive significant “damaging” scores in these predictions. Thus, for this scenario, the use of these function predictions would be limited.

The list of genes that we identified contains other interesting candidates. For instance, IL10RA is a receptor for IL-10, a cytokine that controls inflammatory response (38). Carriers of the minor allele show increased IL10RA expression (Fig. S5), and Il10ra expression is increased in affected brain regions with increasing age and presence of AD pathology in transgenic mouse models of AD [MOUSEAC (39)] (Fig. S9). Moreover, a link between down-regulation of IL10RA and TRIB3 in TRIB3-silenced HepG2 cells was reported in ref. 26 along with increased abundance of Presenilin 1, ApoE3, and Clusterin. Finally, blocking IL-10 response was recently suggested as a therapeutic mechanism in AD (40). A gene that showed a statistical trend in the validation sample was TM2D1, which is a beta-amyloid binding protein and may be involved in beta-amyloid–induced apoptosis (41). Furthermore, Myocyte Enhancer Factor 2A (MEF2A), like APOE, was filtered out by the functional prioritization. However, MEF2A is a paralog of MEF2C, which is an established AD gene (6). Noteworthy, bins covering MEF2C only barely missed the selection threshold in component 2 for additional analysis (maximum P = 0.926) (Fig. 2).

This study illustrates the potential of effectively combining multivariate statistical modeling in imaging–genetics with recent instruments available from computational biology to lead to insights on the disease pathophysiology. Thanks to the ever-growing data-driven knowledge based on the vast quantities of information now available to the research community, the paradigm proposed in this study may represent a promising avenue for linking imaging–genetics findings to the current knowledge on functional genetics mechanisms involved in neurodegeneration.

Materials and Methods

This section describes the study data, the statistical setting, and the methodology used in this study. Additional details and discussion about the methodological aspects can be found in SI Materials and Methods.

Study Participants.

Data used in the preparation of this article were obtained from the ADNI database (adni.loni.usc.edu). The ADNI was launched in 2003 as a public–private partnership led by Principal Investigator Michael W. Weiner. Up-to-date information is available at www.adni-info.org. This research mainly involves further processing of previously collected personal data according to the US ethics regulations. Each subject provided signed informed consent before participation. We have explicit authorization for the use of the ADNI dataset, and we have signed the relevant papers guaranteeing that we abide by the ethics standards.

We selected genotype and phenotype data available in the ADNI-1/GO/2 datasets for 1,192 subjects. Summary sociodemographic, clinical, and genetic information is available in Table 2. At the time of study entry, subjects were diagnosed as healthy individuals (n = 401), MCI (n = 553), or AD (n = 238). A total of 212 (38.3%) MCI patients subsequently converted to AD over the course of the study (6 y). All participants were non-Hispanic Caucasian. AD and MCI groups show significant cognitive decline measured by the Mini Mental State Examination (MMSE) and the Alzheimer’s Disease Assessment Scale (ADAS) Cognitive Subscale (COG) compared with the healthy individuals (P < 1e−2, two-sample t test for groupwise comparison). There was also a significant increase in individuals with pathological levels of Aβ1–42 in the CSF (Aβ1–42 < 192 pg/mL) across the clinical groups, with proportions ranging from 43% in healthy individuals to 93% in AD patients (P < 1e−2). Similarly, we observed a higher prevalence of APOE4 carriers in AD and progressing MCI individuals compared with healthy and MCI stable groups. For this analysis, the 639 healthy and AD subjects form the discovery set, while the MCI converters and nonconverters form the independent validation set.

Statistical Discovery.

The joint relationship between the genetic and imaging modalities was investigated through PLS modeling (42⇓–44). Among the several PLS versions proposed in the literature, we focus on the symmetric formulation of PLS computed through the singular value decomposition of the cross-covariance matrix (Fig. S1) (43, 44, 45). Within this setting, the aim of PLS is to estimate the latent components that maximize the global covariance between the two input modalities. Each input feature receives a weight in the latent component that represents its relative importance for describing the global joint multimodal relationship. Analyzing these weights helps identify SNPs that are linked to the patterns of cortical thinning in the brain.

In this study, we applied a robust approach for the stable estimation and interpretation of PLS weights in genome-wide genotyping data aimed at promoting sparsity (i.e., selecting only a few features for simplified interpretation) and regularity (by aggregating SNPs within the same genetic neighborhood). This is achieved through a stability selection procedure, in which the reproducibility and robustness of the PLS parameters are assessed through a split-half cross-validation–based scheme on 1 million repetitions of the models on randomly sampled subgroups (Fig. 1 and SI Materials and Methods). By considering a predefined partition of each chromosome into contiguous loci of size 10 kb, the procedure leads to the estimation of a confidence measure taking values ranging between 0.0 and 1.0, indicating the probability of each genetic loci containing highly reproducible PLS weights and therefore, serving as a measure of importance of the genomic location (Fig. 2). A similar procedure was used to assess the importance of the phenotype component (Fig. 1). However, no regional binning was used (Fig. 3). The procedure was applied to assess the parameter reproducibility of the first five PLS modes; subsequent analyses were performed only on components with relevant genetic and brain regions (i.e., reproducible PLS weights with selection frequency >95%).

Gene Identification.

We analyzed the 10-kb bins (genetic loci) with the selection frequency exceeding 0.95 (i.e., bins selected in 95% or more of the 1 million replications). Within these bins, we then identified the influential SNPs: an SNP was declared influential if it was associated with the weights of greatest magnitude in the PLS components estimated on the full data sample (i.e., SNPs with absolute weights exceeding the 99th quantile of all weights in the component). These weights are the ones contributing to the high selection frequency in the split-half procedure and are representative of the significant variation modeled in the data.

To link SNPs to corresponding genes, we used the computational VEP for GRCh37 with the GENCODE gene annotation. SNPs tagged as “regulatory” were manually investigated and annotated with the nearby genes.

Functional Prioritization.

All SNPs successfully annotated with a gene were subjected to functional prioritization through eQTL analysis based on the GTEx data. The sample size in GTEx for relevant brain tissues in AD was rather small (e.g., n = 81 for hippocampus). Therefore, we added five more tissues with large sample sizes that were more distantly relevant to AD. Nerve tibial (n = 256) was added as a proxy for nervous tissue, whole blood (n = 338) and artery tibial (n = 285) were included to cover blood-based changes and effects on blood vessels (46), and adipose s.c. (n = 298) was selected due to links between AD and obesity, type 2 diabetes, and metabolic disease (47, 48). Finally, transformed fibroblasts (n = 272) were included as a general purpose cell line. P values were corrected for multiple testing using the Bonferroni method.

Model Validation in Independent MCI Subjects.

The genes that were found to be under expression control by the identified SNPs were validated for their capacity to predict clinical conversion in MCI subjects. To this end, for each identified gene, we applied the PLS weights estimated on the discovery set on the validation set, with the genetic component restricted to SNPs ±20 kb from the gene borders. The identified latent projections (i.e., a weighted sum of SNPs) result in one score per subject per gene. For each gene, the association of the projection score with conversion status was assessed by statistically comparing the scores distribution between healthy individuals and AD patients and between MCI converters and nonconverters (Kruskal–Wallis nonparametric test for two-sample comparison, Bonferroni correction for multiple comparisons).

Acknowledgments

M.L., J.M.S., D.C.A., and S.O. received support from European Union’s Horizon 2020 Research and Innovation Programme Grant 666992 (EuroPOND) for this work. A.A. holds an MRC eMedLab Medical Bioinformatics Career Development Fellowship. This work was supported by Medical Research Council Grant MR/L016311/1. The contribution to this work by B.G. and P.M.T. was funded by NIH “Big Data to Knowledge” Grant U54 EB020403 (principal investigator: P.M.T.). S.W. and C.A. are supported by the National Institute for Health Research (NIHR) Queen Square Biomedical Research Unit in Dementia and Alzheimer’s Research UK and the NIHR University College London (UCL) Hospitals Biomedical Research Centre. J.M.S. acknowledges the support of the NIHR UCL Hospitals Biomedical Research Centre; the Wolfson Foundation; Engineering and Physical Sciences Research Council (EPSRC) Grant EP/J020990/1; Medical Research Council (MRC) Grant MR/L023784/1; Alzheimer Research UK (ARUK) Grants ARUK-Network 2012-6-ICE, ARUK-PG2017-1946, and ARUK-PG2017-1946; Brain Research Trust Grant UCC14191; and European Union’s Horizon 2020 Research and Innovation Programme Grant 666992. EPSRC Grants EP/J020990/01 and EP/M020533/1 support the work of D.C.A. and S.O. on this topic. S.O. receives funding from EPSRC Grants EP/H046410/1 and EP/K005278, MRC Grant MR/J01107X/1, EU-FP7 Project VPH-DARE@IT Grant FP7-ICT-2011-9-601055, the NIHR Biomedical Research Unit (Dementia) at the UCL, and NIHR University College London Hospitals Biomedical Research Centre (NIHR BRC UCLH/UCL High Impact Initiative-BW.mn.BRC10269). Data collection and sharing for this project were funded by the ADNI (NIH Grant U01 AG024904) and the Department of Defense (DOD) ADNI (DOD Award W81XWH-12-2-0012). The ADNI is funded by the National Institute on Aging and the National Institute of Biomedical Imaging and Bioengineering, and it is also funded through generous contributions from the following: AbbVie, Alzheimer’s Association; Alzheimer’s Drug Discovery Foundation; Araclon Biotech; BioClinica, Inc.; Biogen; Bristol-Myers Squibb Company; CereSpir, Inc.; Eisai Inc.; Elan Pharmaceuticals, Inc.; Eli Lilly and Company; EuroImmun; F. Hoffmann-La Roche Ltd. and its affiliated company Genentech, Inc.; Fujirebio; GE Healthcare; IXICO Ltd.; Janssen Alzheimer Immunotherapy Research & Development, LLC; Johnson & Johnson Pharmaceutical Research & Development LLC; Lumosity; Lundbeck; Merck & Co., Inc.; Meso Scale Diagnostics, LLC; NeuroRx Research; Neurotrack Technologies; Novartis Pharmaceuticals Corporation; Pfizer Inc.; Piramal Imaging; Servier; Takeda Pharmaceutical Company; and Transition Therapeutics. The Canadian Institutes of Health Research is providing funds to support ADNI clinical sites in Canada. Private sector contributions are facilitated by the Foundation for the NIH (https://fnih.org/). The grantee organization is the Northern California Institute for Research and Education, and the study is coordinated by the Alzheimer’s Disease Cooperative Study at the University of California, San Diego. The ADNI data are disseminated by the Laboratory for NeuroImaging at the University of Southern California.

Footnotes

  • ↵1M.L. and A.A. contributed equally to this work.

  • ↵2To whom correspondence should be addressed. Email: marco.lorenzi{at}inria.fr.
  • ↵3Data used in preparation of this article were obtained from the Alzheimer’s Disease Neuroimaging Initiative (ADNI) database (adni.loni.usc.edu). As such, the investigators within the ADNI contributed to the design and implementation of ADNI and/or provided data but did not participate in analysis or writing of this report. A complete listing of ADNI investigators can be found at adni.loni.usc.edu/wp-content/uploads/how_to_apply/ADNI_Acknowledgement_List.pdf.

  • Author contributions: M.L., A.A., S.W., D.P.H., N.J., J.M.S., D.C.A., P.M.T., and S.O. designed research; M.L., A.A., and C.A. performed research; M.L., A.A., and B.G. analyzed data; and M.L., A.A., S.W., C.A., and J.M.S. wrote the paper.

  • The authors declare no conflict of interest.

  • This article is a PNAS Direct Submission.

  • Data deposition: All the study data are available through the Alzheimer’s Disease Neuroimaging Initiative (ADNI, adni.loni.usc.edu/). The experimental material, including PLS components and probability measures, has been deposited in Figshare (https://doi.org/10.6084/m9.figshare.5914375.v1).

  • This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1706100115/-/DCSupplemental.

Published under the PNAS license.

View Abstract

References

  1. ↵
    1. Rabinovici GD,
    2. Jagust WJ
    (2009) Amyloid imaging in aging and dementia: Testing the amyloid hypothesis in vivo. Behav Neurol 21:117–128.
    OpenUrlCrossRefPubMed
  2. ↵
    1. Villemagne VL,
    2. Okamura N
    (2014) In vivo tau imaging: Obstacles and progress. Alzheimers Dement 10(3 Suppl):S254–S264.
    OpenUrlCrossRefPubMed
  3. ↵
    1. Mosconi L, et al.
    (2010) Pre-clinical detection of Alzheimer’s disease using FDG-PET, with or without amyloid imaging. J Alzheimers Dis 20:843–854.
    OpenUrlCrossRefPubMed
  4. ↵
    1. Frisoni GB,
    2. Fox NC,
    3. Jack CR Jr,
    4. Scheltens P,
    5. Thompson PM
    (2010) The clinical use of structural MRI in Alzheimer disease. Nat Rev Neurol 6:67–77.
    OpenUrlCrossRefPubMed
  5. ↵
    1. Bigos KL,
    2. Hariri AR,
    3. Weinberger DR
    (2016) Neuroimaging Genetics: Principles and Practices (Oxford Univ Press, Oxford).
  6. ↵
    1. Lambert JC, et al., European Alzheimer’s Disease Initiative (EADI); Genetic and Environmental Risk in Alzheimer’s Disease; Alzheimer’s Disease Genetic Consortium; Cohorts for Heart and Aging Research in Genomic Epidemiology
    (2013) Meta-analysis of 74,046 individuals identifies 11 new susceptibility loci for Alzheimer’s disease. Nat Genet 45:1452–1458.
    OpenUrlCrossRefPubMed
  7. ↵
    1. Potkin SG, et al.
    (2009) Hippocampal atrophy as a quantitative trait in a genome-wide association study identifying novel susceptibility genes for Alzheimer’s disease. PLoS ONE 4:e6501.
    OpenUrlCrossRefPubMed
  8. ↵
    1. Ramanan VK, et al., Alzheimer’s Disease Neuroimaging Initiative
    (2014) APOE and BCHE as modulators of cerebral amyloid deposition: A florbetapir PET genome-wide association study. Mol Psychiatry 19:351–357.
    OpenUrlCrossRefPubMed
  9. ↵
    1. Kam-Thong T, et al.
    (2012) GLIDE: GPU-based linear regression for detection of epistasis. Hum Hered 73:220–236.
    OpenUrlCrossRefPubMed
  10. ↵
    1. Stein JL, et al., Alzheimer’s Disease Neuroimaging Initiative
    (2010) Voxelwise genome-wide association study (vGWAS). Neuroimage 53:1160–1174.
    OpenUrlCrossRefPubMed
  11. ↵
    1. Schrouff J, et al.
    (2013) PRoNTo: Pattern recognition for neuroimaging toolbox. Neuroinformatics 11:319–337.
    OpenUrlCrossRefPubMed
  12. ↵
    1. Szymczak S, et al.
    (2009) Machine learning in genome-wide association studies. Genet Epidemiol 33:S51–S57.
    OpenUrlCrossRefPubMed
  13. ↵
    1. Liu J,
    2. Calhoun VD
    (2014) A review of multivariate analyses in imaging genetics. Front Neuroinform 8:29.
    OpenUrl
  14. ↵
    1. Le Floch E, et al.
    (2012) Significant correlation between a set of genetic polymorphisms and a functional brain network revealed by feature selection and sparse partial least squares. Neuroimage 63:11–24.
    OpenUrlCrossRefPubMed
  15. ↵
    1. Vounou M,
    2. Nichols TE,
    3. Montana G, Alzheimer’s Disease Neuroimaging Initiative
    (2010) Discovering genetic associations with high-dimensional neuroimaging phenotypes: A sparse reduced-rank regression approach. Neuroimage 53:1147–1159.
    OpenUrlCrossRefPubMed
  16. ↵
    1. Silver M,
    2. Janousova E,
    3. Hua X,
    4. Thompson PM,
    5. Montana G, Alzheimer’s Disease Neuroimaging Initiative
    (2012) Identification of gene pathways implicated in Alzheimer’s disease using longitudinal imaging phenotypes with sparse regression. Neuroimage 63:1681–1694.
    OpenUrlCrossRefPubMed
  17. ↵
    1. Liu J, et al.
    (2009) Combining fMRI and SNP data to investigate connections between brain function and genetics using parallel ICA. Hum Brain Mapp 30:241–255.
    OpenUrlCrossRefPubMed
  18. ↵
    1. Carithers LJ, et al., GTEx Consortium
    (2015) A novel approach to high-quality postmortem tissue procurement: The GTEx project. Biopreserv Biobank 13:311–319.
    OpenUrlCrossRefPubMed
  19. ↵
    1. Trabzuni D, et al.
    (2011) Quality control parameters on a large dataset of regionally dissected human control brains for whole genome expression studies. J Neurochem 119:275–282.
    OpenUrlCrossRefPubMed
  20. ↵
    1. Krzywinski M, et al.
    (2009) Circos: An information aesthetic for comparative genomics. Genome Res 19:1639–1645.
    OpenUrlAbstract/FREE Full Text
  21. ↵
    1. Aken BL, et al.
    (2016) The Ensembl gene annotation system. Database (Oxford) 2016:baw093.
    OpenUrlCrossRefPubMed
  22. ↵
    1. Machiela MJ,
    2. Chanock SJ
    (2015) LDlink: A web-based application for exploring population-specific haplotype structure and linking correlated alleles of possible functional variants. Bioinformatics 31:3555–3557.
    OpenUrlCrossRefPubMed
  23. ↵
    1. Du K,
    2. Herzig S,
    3. Kulkarni RN,
    4. Montminy M
    (2003) TRB3: A tribbles homolog that inhibits Akt/PKB activation by insulin in liver. Science 300:1574–1577.
    OpenUrlAbstract/FREE Full Text
  24. ↵
    1. Zareen N,
    2. Biswas SC,
    3. Greene LA
    (2013) A feed-forward loop involving Trib3, Akt and FoxO mediates death of NGF-deprived neurons. Cell Death Differ 20:1719–1730.
    OpenUrlCrossRefPubMed
  25. ↵
    1. Aimé P, et al.
    (2015) Trib3 Is elevated in Parkinson’s disease and mediates death in Parkinson’s disease models. J Neurosci 35:10731–10749.
    OpenUrlAbstract/FREE Full Text
  26. ↵
    1. Hua F, et al.
    (2015) TRB3 links insulin/IGF to tumour promotion by interacting with p62 and impeding autophagic/proteasomal degradations. Nat Commun 13:7951.
    OpenUrl
  27. ↵
    1. Menzies FM,
    2. Fleming A,
    3. Rubinsztein DC
    (2015) Compromised autophagy and neurodegenerative diseases. Nat Rev Neurosci 16:345–357.
    OpenUrlCrossRefPubMed
  28. ↵
    1. Zhou Y, et al.
    (2008) E3 ubiquitin ligase SIAH1 mediates ubiquitination and degradation of TRB3. Cell Signal 20:942–948.
    OpenUrlCrossRefPubMed
  29. ↵
    1. Saleem S,
    2. Biswas SC
    (2017) Tribbles Pseudokinase 3 induces both apoptosis and autophagy in amyloid-β induced neuronal death. J Biol Chem 292:2571–2585.
    OpenUrlAbstract/FREE Full Text
  30. ↵
    1. Zhang W, et al.
    (2016) Skeletal muscle TRIB3 mediates glucose toxicity in diabetes and high- fat diet-induced insulin resistance. Diabetes 65:2380–2391.
    OpenUrlAbstract/FREE Full Text
  31. ↵
    1. Sims-Robinson C,
    2. Kim B,
    3. Rosko A,
    4. Feldman EL
    (2010) How does diabetes accelerate Alzheimer disease pathology? Nat Rev Neurol 6:551–559.
    OpenUrlCrossRefPubMed
  32. ↵
    1. Ribe EM,
    2. Lovestone S
    (2016) Insulin signalling in Alzheimer’s disease and diabetes: From epidemiology to molecular links. J Intern Med 280:430–442.
    OpenUrl
  33. ↵
    1. Morris AP, et al., Wellcome Trust Case Control Consortium; Meta-Analyses of Glucose and Insulin-related traits Consortium (MAGIC) Investigators; Genetic Investigation of ANthropometric Traits (GIANT) Consortium; Asian Genetic Epidemiology Network–Type 2 Diabetes (AGEN-T2D) Consortium; South Asian Type 2 Diabetes (SAT2D) Consortium; DIAbetes Genetics Replication And Meta-analysis (DIAGRAM) Consortium
    (2012) Large-scale association analysis provides insights into the genetic architecture and pathophysiology of type 2 diabetes. Nat Genet 44:981–990.
    OpenUrlCrossRefPubMed
  34. ↵
    1. Luciano M, et al.
    (2011) Whole genome association scan for genetic polymorphisms influencing information processing speed. Biol Psychol 86:193–202.
    OpenUrlCrossRefPubMed
  35. ↵
    1. Giedraitis V, et al.
    (2009) Genetic analysis of Alzheimer’s disease in the uppsala longitudinal study of adult men. Dement Geriatr Cogn Disord 27:59–68.
    OpenUrlCrossRefPubMed
  36. ↵
    1. Li H, et al.
    (2008) Candidate single-nucleotide polymorphisms from a genomewide association study of Alzheimer disease. Arch Neurol 65:45–53.
    OpenUrlCrossRefPubMed
  37. ↵
    1. Oláh J, et al.
    (2011) Interactions of pathological hallmark proteins: Tubulin polymerization promoting protein/p25, β-amyloid, and α-synuclein. J Biol Chem 286:34088–34100.
    OpenUrlAbstract/FREE Full Text
  38. ↵
    1. Li MO,
    2. Flavell RA
    (2008) Contextual regulation of inflammation: A duet by transforming growth factor-β and interleukin-10. Immunity 28:468–476.
    OpenUrlCrossRefPubMed
  39. ↵
    1. Matarin M, et al.
    (2015) A genome-wide gene-expression analysis and database in transgenic mice during development of amyloid or tau pathology. Cell Rep 10:633–644.
    OpenUrlPubMed
  40. ↵
    1. Guillot-Sestier MV, et al.
    (2015) Il10 deficiency rebalances innate immunity to mitigate Alzheimer-like pathology. Neuron 85:534–548.
    OpenUrlCrossRefPubMed
  41. ↵
    1. Kajkowski EM, et al.
    (2001) β -Amyloid peptide-induced apoptosis regulated by a novel protein containing a g protein activation module. J Biol Chem 276:18748–18756.
    OpenUrlAbstract/FREE Full Text
  42. ↵
    1. Wold H
    (1966) Estimation of principal components and related models by iterative least squares. Multivariate Analysis (Academic, New York), pp 391–420.
  43. ↵
    1. McIntosh AR,
    2. Bookstein FL,
    3. Haxby JV,
    4. Grady CL
    (1996) Spatial pattern analysis of functional brain images using partial least squares. Neuroimage 3:143–157.
    OpenUrlCrossRefPubMed
  44. ↵
    1. Worsley KJ
    (1997) An overview and some new developments in the statistical analysis of PET and fMRI data. Hum Brain Mapp 5:254–258.
    OpenUrlPubMed
  45. ↵
    1. Friston KJ,
    2. Frith CD,
    3. Liddle PF,
    4. Frackowiak RSJ
    (1993) Functional connectivity: The principal-component analysis of large (PET) data sets. J Cereb Blood Flow Metab 13:5–14.
    OpenUrlCrossRefPubMed
  46. ↵
    1. Kimbrough IF,
    2. Robel S,
    3. Roberson ED,
    4. Sontheimer H
    (2015) Vascular amyloidosis impairs the gliovascular unit in a mouse model of Alzheimer’s disease. Brain 138:3716–3733.
    OpenUrlCrossRefPubMed
  47. ↵
    1. Luchsinger JA,
    2. Gustafson DR
    (2009) Adiposity, type 2 diabetes, and Alzheimer’s disease. J Alzheimers Dis 4:693–704.
    OpenUrl
  48. ↵
    1. Ferreira ST,
    2. Clarke JR,
    3. Bomfim TR,
    4. De Felice FG
    (2014) Inflammation, defective insulin signaling, and neuronal dysfunction in Alzheimer’s disease. Alzheimers Dement 10(1 Suppl):S76–S83.
    OpenUrlCrossRefPubMed
PreviousNext
Back to top
Article Alerts
Email Article

Thank you for your interest in spreading the word on PNAS.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Susceptibility of brain atrophy to TRIB3 in Alzheimer’s disease, evidence from functional prioritization in imaging genetics
(Your Name) has sent you a message from PNAS
(Your Name) thought you would like to see the PNAS web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Susceptibility of brain atrophy to TRIB3 in Alzheimer’s disease, evidence from functional prioritization in imaging genetics
Marco Lorenzi, Andre Altmann, Boris Gutman, Selina Wray, Charles Arber, Derrek P. Hibar, Neda Jahanshad, Jonathan M. Schott, Daniel C. Alexander, Paul M. Thompson, Sebastien Ourselin, for the Alzheimer’s Disease Neuroimaging Initiative
Proceedings of the National Academy of Sciences Mar 2018, 115 (12) 3162-3167; DOI: 10.1073/pnas.1706100115

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
Susceptibility of brain atrophy to TRIB3 in Alzheimer’s disease, evidence from functional prioritization in imaging genetics
Marco Lorenzi, Andre Altmann, Boris Gutman, Selina Wray, Charles Arber, Derrek P. Hibar, Neda Jahanshad, Jonathan M. Schott, Daniel C. Alexander, Paul M. Thompson, Sebastien Ourselin, for the Alzheimer’s Disease Neuroimaging Initiative
Proceedings of the National Academy of Sciences Mar 2018, 115 (12) 3162-3167; DOI: 10.1073/pnas.1706100115
Digg logo Reddit logo Twitter logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Mendeley logo Mendeley
Proceedings of the National Academy of Sciences: 115 (12)
Table of Contents

Submit

Sign up for Article Alerts

Article Classifications

  • Biological Sciences
  • Neuroscience
  • Physical Sciences
  • Computer Sciences

Jump to section

  • Article
    • Abstract
    • Results
    • Discussion
    • Materials and Methods
    • Acknowledgments
    • Footnotes
    • References
  • Figures & SI
  • Info & Metrics
  • PDF

You May Also be Interested in

Abstract depiction of a guitar and musical note
Science & Culture: At the nexus of music and medicine, some see disease treatments
Although the evidence is still limited, a growing body of research suggests music may have beneficial effects for diseases such as Parkinson’s.
Image credit: Shutterstock/agsandrew.
Large piece of gold
News Feature: Tracing gold's cosmic origins
Astronomers thought they’d finally figured out where gold and other heavy elements in the universe came from. In light of recent results, they’re not so sure.
Image credit: Science Source/Tom McHugh.
Dancers in red dresses
Journal Club: Friends appear to share patterns of brain activity
Researchers are still trying to understand what causes this strong correlation between neural and social networks.
Image credit: Shutterstock/Yeongsik Im.
Yellow emoticons
Learning the language of facial expressions
Aleix Martinez explains why facial expressions often are not accurate indicators of emotion.
Listen
Past PodcastsSubscribe
Goats standing in a pin
Transplantation of sperm-producing stem cells
CRISPR-Cas9 gene editing can improve the effectiveness of spermatogonial stem cell transplantation in mice and livestock, a study finds.
Image credit: Jon M. Oatley.

Similar Articles

Site Logo
Powered by HighWire
  • Submit Manuscript
  • Twitter
  • Facebook
  • RSS Feeds
  • Email Alerts

Articles

  • Current Issue
  • Special Feature Articles – Most Recent
  • List of Issues

PNAS Portals

  • Anthropology
  • Chemistry
  • Classics
  • Front Matter
  • Physics
  • Sustainability Science
  • Teaching Resources

Information

  • Authors
  • Editorial Board
  • Reviewers
  • Librarians
  • Press
  • Site Map
  • PNAS Updates

Feedback    Privacy/Legal

Copyright © 2021 National Academy of Sciences. Online ISSN 1091-6490