Genetic Association of Phosphodiesterases With Human Cognitive Performance

Introduction

Multiple recent, large-scale, genome-wide association studies (GWAS) provide an unbiased view of the genetic fine structure of cognitive function in healthy individuals (Lam et al., 2017; Davies et al., 2018; Lee et al., 2018; Savage et al., 2018). The studies ranged in size from 107,207 to 1.1 million individuals. The collection of phenotypic information across such a large number of subjects required the use of very simple phenotypes of cognitive ability. In the largest GWAS study, which was conducted by the Social Science Genetic Association Consortium (SSGAC) and 23andMe, educational attainment was supplemented by analyses of self-reported math ability, highest math class, and by a composite score for cognitive performance that allowed the pooling of data collected from subjects across multiple studies that utilized different cognitive assessment tools (Lee et al., 2018). The study was limited to subjects of European ancestry and had lower predictive power in African-Americans. Educational attainment more broadly captures cognitive ability, but also personality and motivation, as well as behavior traits such as orderliness, lack of impulsiveness, and perseverance. The conclusions of the study are subject to the concerns that educational attainment is strongly influenced by socioeconomic status, by parental education level, nutrition, early childhood experience, the quality of schooling, and other environmental factors. For example, household income explains 7% of the population variation in educational attainment in the GWAS cohort (Lee et al., 2018). Remarkably, genetic variation was found to account for up to 11% of the population variation in educational attainment.

The genetic contribution to educational attainment comes from many genes of small effect. The GWAS reported by SSGAC and 23andMe identified 1,271 single nucleotide polymorphisms (SNPs) that reached genome-wide significance. The genes identified were overwhelmingly expressed in neurons. Genes expressed principally in glial cells were largely absent (Lee et al., 2018). The GWAS thus provides a rich resource for testing the contribution of candidate genes to human cognition. This allows the interrogation of gene families and pathways for potential targets of therapeutic intervention.

Genetic Association of Phosphodiesterases (PDE) With Educational Attainment and Cognitive Performance

The phosphodiesterase (PDE) enzymes hydrolyze cyclic AMP or cyclic GMP, key intracellular signaling molecules, to the respective monophosphate (AMP or GMP). They comprise a large gene family whose protein products can be modulated by small molecule therapeutics (Bender and Beavo, 2006; Conti and Beavo, 2007). Cilostazol, a small molecule inhibitor of PDE3 is used to treat intermittent claudication in peripheral arterial occlusive disorder; sildenafil, tadalafil, and vardenafil, inhibitors of PDE5, are used to treat male erectile dysfunction; and roflumilast and apremilast, which inhibit all four subtypes of PDE4, are used to treat chronic obstructive pulmonary disorder, psoriasis or psoriatic arthritis (Yuan et al., 2013; Bedenis et al., 2014; Deeks, 2015; Garnock-Jones, 2015). Many of the PDE gene family members are expressed in brain and have been utilized as targets for CNS drug discovery. PDE4, PDE9 and PDE10 inhibitors have been explored in major depression, Alzheimer’s disease and schizophrenia, although clinical data are largely negative (Zhu et al., 2001; Wu et al., 2018).

The large GWAS dataset from the SSGCA and 23andMe study provides a means to test the contribution of genetic variation in PDE gene family members to cognitive performance (Figure 1). SNP in PDE1A, PDE1C, PDE2A, PDE4B and PDE4D reached genome-wide significance (P < 5 × 10⁻⁸) in the primary GWAS which used educational attainment as a phenotype (Table 1). Several informative SNP in PDE1A (rs1835339) and PDE4B (rs556755587, rs11208742) were located immediately upstream of the respective genes, while the remaining SNP were intragenic.

FIGURE 1

Figure 1. (A) Manhattan plot showing the location and P-values for single nucleotide polymorphisms (SNPs) reaching genome-wide significance in the primary genome-wide association genome-wide association study (GWAS) of educational attainment (Lee et al., 2018). The locations of genes encoding phosphodiesterase (PDE) are shown. (B) SNP reaching genome-wide significance on chromosome 5. Informative SNP covering the PDE4D gene are colored blue. (C) Genetic fine structure of the PDE4B gene. SNP achieving genome-wide significance in the GWAS (•), conditional-joint analysis (COJO; +) or multi-trait analysis of GWAS (MTAG; ▾) analyses cluster at the 5’end of the gene. These overlap SNP associating with self-reported math ability (◆) and highest math class taken (•). The 5’ exons of mRNA encoding the major PDE4B isoforms are labeled. A group of 3’ exons encoding the catalytic domain are common to all PDE4B mRNA. (D) Genetic fine structure of the PDE4D gene. SNP achieving genome-wide significance in the GWAS (•), COJO (+) or MTAG (▾) analyses cluster at the 5’end of the gene. These overlap SNP associating with self-reported math ability(◆), cognitive performance (■) and SNP associating with intellect (◆) in an independent study by Savage et al. (2018). Genome position is based on the GRCh37.p13: Annotation Release 105.

TABLE 1

Table 1. Genetic association results for phosphodiesterase (PDE) genes.

The SSGCA and 23andMe study also subjected the GWAS dataset to a more stringent, conditional-joint analysis (COJO) (Yang et al., 2012). The genetic association of PDE1A was lost, while SNP covering PDE1C, PDE2A, PDE4B and PDE4D continued to achieve genome-wide significance (Table 1). The SNP data also were analyzed using multi-trait analysis of GWAS (MTAG) which increases the number of SNP identified at genome-wide significance (Turley et al., 2018). The MTAG analysis focused on three complementary cognitive phenotypes: self-reported math ability, highest math class taken, and cognitive performance. The two math phenotypes were self-reported by participants in 23andMe. The measure of cognitive performance was constructed through a meta-analysis of published results from the COGENT consortium (Trampush et al., 2017). MTAG confirms the genetic association with educational years for PDE1A, PDE1C, PDE2A, PDE4B and PDE4D (Table 1). PDE1C also shows genetic association by MTAG across all of the cognitive subanalyses. Multiple SNP over PDE4B and PDE4D reach genome-wide significance for highest math class taken, while only a single SNP upstream of PDE4B continues to reach genome-wide significance for self-reported math ability (rs11208742; P = 9.6 × 10⁻¹⁰) and a single SNP in PDE4D reaches genome-wide significance for cognitive performance (rs7735959; P = 3.89 × 10⁻¹²). As a credibility check, the MTAG for cognitive performance was repeated with the COGENT cohort excluded. SNP in PDE1C and PDE4D continue to reach genome-wide significance. The complete data listings extracted from the SSGCA and 23andMe study are presented in Supplementary Tables S1–S6.

Notably absent from the GWAS are genes previously pursued as targets for cognitive enhancing drugs. Absent are components of the cholinergic neurotransmitter pathway [acetylcholinesterase (ACHE), choline acetyltransferase (CHAT), nicotinic alpha 7 acetylcholine receptor (CHRNA7)], the serotonergic pathway (HT6R, serotonin receptor 6), or the histaminergic (HRH3, histamine receptor 3) and glutamatergic pathways (GRM5 -metabotropic glutamate receptor 5). Nor do components of the amyloid pathway [amyloid protein precursor (APP); apolipoprotein E (APOE)] or components of the complement pathway linked to synaptic pruning by microglial cells (C4A, complement 4A; C4B, complement 4B) show genetic association with cognition in the GWAS data set (Supplementary Table S7).

The genes for PDE4B and PDE4D have dense SNP coverage due to their large size. On average, human genes are 67,000 base pairs (bps) in length (Piovesan et al., 2016). GWAS studies typically are conducted with methods that allow the scanning of 1 million or more SNP per genome. These on average are located every 1,000–3,000 bps across the 3.1 billion bps of the human genome. Thus, a typical gene is covered by 20–30 SNP in a GWAS. Adjacent SNP identify blocks of DNA that are in linkage disequilibrium (LD). The large size of the genes for PDE4B (582,069 bp) and PDE4D (1,519,060 bp) provides proportionally denser coverage in a GWAS. Given that each of these genes utilize multiple promoters and complex splicing to generate multiple protein isoforms, (Gretarsdottir et al., 2003; Fatemi et al., 2008; Guan et al., 2012) the density of SNP coverage provides insight into the genetic fine structure of the PDE4B and PDE4D genes (Figure 1).

PDE4B

PDE4 genes encode multiple protein isoforms that are present as either dimers or monomers within the cell (Figure 2). Informative SNP in the gene encoding PDE4B cluster at the 5’ end of the gene across exons that through differential promoter utilization and splicing encode dimeric forms of the enzyme known as PDE4B4, B1 and B3 with each transcript beginning with a unique, 5’ exon (Bender and Beavo, 2006). The PDE4B gene also encodes transcripts for monomeric forms of PDE4B known as PDE4B2 and PDE4B5 (Campbell et al., 2017). These transcripts utilize a downstream promoter and lack the upstream exons needed to encode the protein domains necessary for dimerization (Richter and Conti, 2002; Cedervall et al., 2015). The absence of genetic association with 3’ SNP indicates that the dimeric forms as compared to the monomeric forms of PDE4B have greater effect on human cognition. Activity of the PDE4B dimer is increased by protein kinase A phosphorylation, (Houslay and Adams, 2003) whereas levels of the monomeric PDE4B isoforms are controlled transcriptionally (Jin and Conti, 2002). Transcriptional up-regulation of PDE4B2 in monocytes and microglia allows production of Tumor Necrosis Factor-α (TNF α) by inflammatory cells in response to Toll receptor-4 pathway activation (e.g., by bacterial lipopolysaccharide or LPS; Jin and Conti, 2002; Zhang et al., 2002; Jin et al., 2005; Wilson et al., 2017).

FIGURE 2

Figure 2. (A) PDE4 protein isoforms. Each of the PDE4 genes produces multiple transcripts that vary in 5’ exon structure (Bender and Beavo, 2006). Dimeric forms of PDE4 enzymes are encoded by long transcripts that vary in 5’ exons but all contain 3’ exons encoding the dimerization, linker and catalytic domains. Short and super-short transcripts encoding monomeric isoforms lack exons encoding the dimerization domain. Basal PDE4 enzymatic activity is increased by PKA phosphorylation of UCR1 (Sette and Conti, 1996; Hoffmann et al., 1998). (B) Protein structure of the PDE4 dimer based on PDB ID:4WZI (Cedervall et al., 2015). The dimerization domain is formed by a four helix bundle comprising the last alpha-helix in UCR1 and the first alpha-helix in UCR2 (Cedervall et al., 2015; Richter and Conti, 2002). A regulatory helix from UCR2 is positioned in trans-across the active site of the opposite monomer and thereby regulates hydrolytic activity by controlling access to cAMP (Burgin et al., 2010).

PDE4B gene variants previously were linked to risk for schizophrenia in a study of 878 schizophrenic or schizoaffective patients and 604 controls that included both Caucasian and African American subjects. In contrast to SNP associated with cognition, SNP associated with schizophrenia clustered at the 3’ end of the PDE4B gene in a block of LD extending from 66,750,560 to 66,829,764 bp on Chr 1. Genetic association of the 3’ end of the PDE4B gene with schizophrenia was replicated in Japanese and Han Chinese (Numata et al., 2008; Guan et al., 2012). This LD blocks contains exons encoding transcripts for the PDE4B2 and PDE4B5 monomeric forms of PDE4B. If variant monomeric forms of PDE4B are causal in schizophrenia, their expression in brain microglia could be relevant to the etiology of the disease though an effect on synaptic pruning (Sekar et al., 2016).

Dimeric forms of PDE4B interact with a binding partner known as disrupted in schizophrenia 1 (DISC1), which holds the enzyme in a closed, inactive conformation (Millar et al., 2005). DISC1 was identified in a Scottish kindred in association with psychiatric disorders that also affect cognition (e.g., schizophrenia, bipolar disorder and major depression; Millar et al., 2000; Blackwood and Muir, 2004). However, no SNP covering DISC1 reach genome-wide significance in the GWAS data set.

PDE4D

The much larger PDE4D gene is covered by 16 informative SNP that are all intragenic (Figures 1, 2). These cluster over the 5’ end of the gene and four of the seven SNP cover the first three exons of the transcript for a dimeric form of PDE4D known as PDE4D7 (Gretarsdottir et al., 2003; Bender and Beavo, 2006). Two additional SNP cover the 5’ end exons for transcripts encoding the PDE4D4 and PDE4D5 dimeric forms of the enzyme. PDE4D7, D4 and D5 are expressed in brain and their enzymatic activity is increased by PKA phosphorylation (Richter et al., 2005).

An independent GWAS in 269,867 individuals of European ancestry using a meta-analysis of scores on neurocognitive tests independently confirmed association of gene variants in PDE4D with cognitive performance (Savage et al., 2018). The analysis identified independent significant SNP and then grouped these under a lead SNP to identify regions in approximate linkage equilibrium with each other. For PDE4D, the lead SNP was identified as rs34426618. This identified an LD block on chromosome 5 (bps 59,572,804–59,666,806) at the 5’ end of the PDE4D gene between the first two exons of the PDE4D7 transcript (D7A1 and D7A2) with genome-wide significance for cognition (P = 2.16 × 10⁻⁸). SNP over other PDE were not associated with neurocognitive performance. Davies et al. (2018) in a third GWAS of 300,486 also identified a 5’ SNP over PDE4D (rs36004779) at genome-wide significance. Lam et al. (2017) previously reported genome-wide significance for PDE4D in a GWAS incorporating 107,207 individuals of European ancestry. Xu et al. (2017) further reported suggestive association of PDE4D with cognition in a Chinese study of adult twins. Thus, the 5’end of the PDE4D gene that encodes dimeric isoforms of the enzyme is a consistent hit across multiple GWAS using varying measures of neurocognitive performance.

Genetic association of PDE4D with stroke was shown previously in an Icelandic GWAS study (Gretarsdottir et al., 2003). As in the SSGCA and 23andMe study, informative SNP were located over 5’ exons encoding dimeric forms of PDE4D including three LD blocks over the 5’ exons encoding the dimeric PDE4D7 isoform. The Icelandic study identified both at-risk and protective haplotypes and showed that the at-risk haplotype altered expression of particular PDE4D transcripts such as PDE4D7. Thus, genetic association of 5’ exons of PDE4D with educational attainment and cognitive performance may also be due to subtle changes in expression of particular PDE4D isoforms.

Exome sequencing in children with a rare neurodevelopmental disorder identified PDE4D mutations as causative for acrodysostosis type 2 with or without hormone resistance (ACRDYS2; Lee et al., 2012; Linglart et al., 2012; Michot et al., 2012). ACRDYS2 causes severe intellectual disability accompanied by brachydactyly (short fingers and toes) and nasal hypoplasia. These ultra-rare mutations in PDE4D arise de novo during gametogenesis. The ACRDYS2 mutations are autosomal dominant, and thus of large effect, and are exclusively missense mutations. These mutations contrast with the PDE4D gene variants in the GWAS as these are common variants of small effect. As do the SNP showing genetic association in the GWAS scan, the ACRDYS2 mutations implicate the PDE4D dimer in normal brain function. Mutations alter amino acid residues needed for PDE4D dimerization and affect regulation of enzymatic activity by PKA phosphorylation (Cedervall et al., 2015; Gurney et al., 2015; Motte et al., 2017).

The high sequence conservation of PDE4D is noteworthy. PDE4D7, a representative dimeric form of PDE4D, is highly conserved across species with >99% amino acid sequence conservation across the 748 amino acid length of the polypeptide. The core amino acid sequence of PDE4D7 from UCR1 through the end of the catalytic domain is absolutely conserved across all species except for a phenylalanine in UCR2 that is unique to primate (phenylalanine271) which in non-primate species is a tyrosine (Burgin et al., 2010). That single amino acid difference may be an important biological adaption in primates as phenylalanine271 projects from UCR2 into the catalytic site. Modeling studies suggest that in non-primate species, the tyrosine hydroxyl forms a hydrogen bond with the cAMP ribose (Burgin et al., 2010). This affects K_M as shown be reciprocal exchange of the residues between PDE4D and PDE4B with the phenylalanine lowering K_M from 6 μM to 1 μM. Thus, primate PDE4D may function within cells at lower cAMP concentration as compared to the non-primate enzyme. This could be relevant to adaptions in primates that enhance dendritic compartmentalization in cortical neurons (Beaulieu-Laroche et al., 2018).

PDE4D-PKA-CREB Pathway Analysis

The role of PDE4 in learning and memory has been studied extensively in the Drosophila fruit fly and in mice. Mutational studies in fruit flies identified PDE4 as a key modulator of a signaling pathway critical for associative learning, courtship behavior and neurodevelopment (Byers et al., 1981). Deletions of the single PDE4 gene in Drosophila, known as dnc -dunce, disrupt learning and memory by preventing the hydrolysis of cAMP and thereby altering the normal spatial and temporal patterning of cAMP signaling (Byers et al., 1981; Gervasi et al., 2010). This is a key aspect of cAMP signaling across multiple cells and tissues (Willoughby et al., 2006; Castro et al., 2010; Blackman et al., 2011). Mutations in a second gene, known as rut-rutabaga, identified calcium-calmodulin dependent adenylate cyclase (ADCY) as the upstream source of cAMP important for learning and memory (Levin et al., 1992). In mammals, two calcium-calmodulin dependent ADCYs (ADCY1 and ADCY8) link Ca²⁺ influx through the N-methyl D-aspartate (NMDA) receptor to early stages of memory formation (Wong et al., 1999). Further genetic analysis in mice and flies identified CREB -cAMP response element binding protein, as the key downstream effector (Milner et al., 1998). CREB acts with a second protein, CREBBP—CREB binding protein, to regulate gene transcription (Chrivia et al., 1993). Haplo-insufficiency in humans causes Rubinstein-Taybi Syndrome, a neurodevelopmental disorder that causes intellectual disability, while the homozygous gene deletions are lethal (Petrij et al., 1995; Tanaka et al., 2000). However, SNP over ADCY1, ADCY8, CREB1, and CREBBP do not reach genome-wide significance in the GWAS. This lack of genetic association with upstream and downstream effector molecules highlights the tight biochemical regulation of the key signaling molecule, cAMP, by a critical hydrolytic enzyme, PDE4D.

The gene for brain derived neurotrophic factor (BDNF), contains promoter elements known as cAMP response elements that are required for CREB-mediated gene expression (Tabuchi et al., 2002). Studies in mice indicate that PDE4D modulates expression of BDNF through phosphorylation of CREB (Zhang et al., 2018). SNP in BDNF reach genome-wide significance in the GWAS, COJO and MTAG for educational attainment and in the MTAG for highest math class (Table 2). BDNF is processed to the mature form by proteolysis of a pro-peptide (Leibrock et al., 1989). A single amino acid polymorphism in the pro-peptide, Val66Met is associated with major depression and response to anti-depressant treatment (Tsai, 2018). Two of the SNP identified in the GWAS map to the 5’ end of the BDNF (rs962369 and rs12273363) and therefore potentially mark an LD block containing the Val66Met variant. The BDNF Val66Met variant is associated with increased decline in cognitive performance in Alzheimer’s disease (Gomar et al., 2016; Boots et al., 2017). The GWAS implies that genetic variation in BDNF may impact more broadly on cognition in healthy individuals.

TABLE 2

Table 2. Association Results for brain derived neurotrophic factor (BDNF).

Conclusions

The SSGCA and 23andMe GWAS study of educational attainment and cognitive performance identified 1,271 independent SNP at genome-wide significance which cover 10% of the 20,000 or so genes encoded by the human genome. The vast majority of these SNP occur in intragenic regions of uncertain function. Recombination events during individual meiosis to create a sperm or an egg are relatively rare. Thus, adjacent regions of DNA tend to be transmitted as blocks from parents to child. This causes adjacent genetic markers on a single chromosome to be in LD due to their lower probability of being separated by recombination during meiosis. Across a population, meiotic events are pooled in a GWAS study, such that the blocks of DNA in LD become smaller. This has two effects in a GWAS study, genome-wide significance may be “inflated” by genetic markers that are in LD, and more importantly, the informative SNP may not be the cause of the biological variation, but instead mark a block of DNA that contains the key variant.

This methodology was good enough to reflect known biology as shown by an analysis of the PDE gene family. PDE1A, PDE1C, PDE2, PDE4B and PDE4D have plausible links to cognitive function. Selective inhibitors improve cognitive performance in diverse animal models although human clinical data are not yet available (Heckman et al., 2015). Further genetic data link PDE4B and PDE4D to human cognition. Genetic variation in PDE4B has been linked to schizophrenia, while ACRYDS2 mutations in PDE4D severely affect cognitive function. (Fatemi et al., 2008; Lee et al., 2012; Linglart et al., 2012; Michot et al., 2012). As shown in the Icelandic study, SNP in the PDE4D gene thought to be associated with stroke had subtle effects on PDE4D expression and even on the expression of transcripts encoding different PDE4D isoforms (Gretarsdottir et al., 2003). Thus, the PDE enzymes provide a rich source of targets for drug discovery and development. The recent large scale GWAS studies of human cognitive ability suggest the potential of PDE1C, PDE4B and PDE4D inhibitors to address cognitive impairment across multiple CNS disorders.

Author Contributions

MG wrote the manuscript.

Funding

This work was supported by a research grant from the National Institute of Mental Health (grant number MH091791) and by Tetra Discovery Partners, Inc.

Conflict of Interest Statement

MG is an employee of Tetra Discovery Partners, Inc. which has a financial interest in the discovery and development of PDE4B and PDE4D allosteric inhibitors for the treatment of Fragile X Syndrome, Alzheimer’s disease and other CNS disorders.

Acknowledgments

I wish to thank Dr. Alex Burgin, Dr. Lance Stewart, Dr. Lawrence Wennogle, Prof. Kathryn Paige Harden, and Dr. Aysu Okbay for commenting on the manuscript.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnmol.2019.00022/full#supplementary-material

References

Beaulieu-Laroche, L., Toloza, E. H. S., van der Goes, M. S., Lafourcade, M., Barnagian, D., Williams, Z. M., et al. (2018). Enhanced dendritic compartmentalization in human cortical neurons. Cell 175, 643.e14–651.e14. doi: 10.1016/j.cell.2018.08.045

PubMed Abstract | CrossRef Full Text | Google Scholar

Bedenis, R., Stewart, M., Cleanthis, M., Robless, P., Mikhailidis, D. P., and Stansby, G. (2014). Cilostazol for intermittent claudication. Cochrane Database Syst. Rev. 10:CD003748. doi: 10.1002/14651858.CD003748.pub4

PubMed Abstract | CrossRef Full Text | Google Scholar

Bender, A. T., and Beavo, J. A. (2006). Cyclic nucleotide phosphodiesterases: molecular regulation to clinical use. Pharmacol. Rev. 58, 488–520. doi: 10.1124/pr.58.3.5

PubMed Abstract | CrossRef Full Text | Google Scholar

Blackman, B. E., Horner, K., Heidmann, J., Wang, D., Richter, W., Rich, T. C., et al. (2011). PDE4D and PDE4B function in distinct subcellular compartments in mouse embryonic fibroblasts. J. Biol. Chem. 286, 12590–12601. doi: 10.1074/jbc.m110.203604

PubMed Abstract | CrossRef Full Text | Google Scholar

Blackwood, D. H., and Muir, W. J. (2004). Clinical phenotypes associated with DISC1, a candidate gene for schizophrenia. Neurotox Res. 6, 35–41. doi: 10.1007/bf03033294

PubMed Abstract | CrossRef Full Text | Google Scholar

Boots, E. A., Schultz, S. A., Clark, L. R., Racine, A. M., Darst, B. F., Koscik, R. L., et al. (2017). BDNF Val66Met predicts cognitive decline in the Wisconsin Registry for Alzheimer’s Prevention. Neurology 88, 2098–2106. doi: 10.1212/WNL.0000000000003980

PubMed Abstract | CrossRef Full Text | Google Scholar

Burgin, A. B., Magnusson, O. T., Singh, J., Witte, P., Staker, B. L., Bjornsson, J. M., et al. (2010). Design of phosphodiesterase type 4D (PDE4D) allosteric modulators for cognition with improved safety. Nat. Biotechnol. 28, 63–70. doi: 10.1038/nbt.1598

PubMed Abstract | CrossRef Full Text | Google Scholar

Byers, D., Davis, R. L., and Kiger, J. A. Jr. (1981). Defect in cyclic AMP phosphodiesterase due to the dunce mutation of learning in Drosophila melanogaster. Nature 289, 79–81. doi: 10.1038/289079a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Campbell, S. L., van Groen, T., Kadish, I., Smoot, L. H. M., and Bolger, G. B. (2017). Altered phosphorylation, electrophysiology, and behavior on attenuation of PDE4B action in hippocampus. BMC Neurosci. 18:77. doi: 10.1186/s12868-017-0396-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Castro, L. R., Gervasi, N., Guiot, E., Cavellini, L., Nikolaev, V. O., Paupardin-Tritsch, D., et al. (2010). Type 4 phosphodiesterase plays different integrating roles in different cellular domains in pyramidal cortical neurons. J. Neurosci. 30, 6143–6151. doi: 10.1523/JNEUROSCI.5851-09.2010

PubMed Abstract | CrossRef Full Text | Google Scholar

Cedervall, P., Aulabaugh, A., Geoghegan, K. F., McLellan, T. J., and Pandit, J. (2015). Engineered stabilization and structural analysis of the autoinhibited conformation of PDE4. Proc. Natl. Acad. Sci. U S A 112, E1414–E1422. doi: 10.1073/pnas.1419906112

PubMed Abstract | CrossRef Full Text | Google Scholar

Chrivia, J. C., Kwok, R. P., Lamb, N., Hagiwara, M., Montminy, M. R., and Goodman, R. H. (1993). Phosphorylated CREB binds specifically to the nuclear protein CBP. Nature 365, 855–859. doi: 10.1038/365855a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Conti, M., and Beavo, J. (2007). Biochemistry and physiology of cyclic nucleotide phosphodiesterases: essential components in cyclic nucleotide signaling. Annu. Rev. Biochem. 76, 481–511. doi: 10.1146/annurev.biochem.76.060305.150444

PubMed Abstract | CrossRef Full Text | Google Scholar

Davies, G., Lam, M., Harris, S. E., Trampush, J. W., Luciano, M., Hill, W. D., et al. (2018). Study of 300,486 individuals identifies 148 independent genetic loci influencing general cognitive function. Nat. Commun. 9:2098. doi: 10.1038/s41467-018-04362-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Deeks, E. D. (2015). Apremilast: a review in psoriasis and psoriatic arthritis. Drugs 75, 1393–1403. doi: 10.1007/s40265-015-0439-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Fatemi, S. H., King, D. P., Reutiman, T. J., Folsom, T. D., Laurence, J. A., Lee, S., et al. (2008). PDE4B polymorphisms and decreased PDE4B expression are associated with schizophrenia. Schizophr. Res. 101, 36–49. doi: 10.1016/j.schres.2008.01.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Garnock-Jones, K. P. (2015). Roflumilast: a review in COPD. Drugs 75, 1645–1656. doi: 10.1007/s40265-015-0463-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Gervasi, N., Tchenio, P., and Preat, T. (2010). PKA dynamics in a Drosophila learning center: coincidence detection by rutabaga adenylyl cyclase and spatial regulation by dunce phosphodiesterase. Neuron 65, 516–529. doi: 10.1016/j.neuron.2010.01.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Gomar, J. J., Conejero-Goldberg, C., Huey, E. D., Davies, P., Goldberg, T. E., and Initiative Alzheimer’s Disease Neuroimaging. (2016). Lack of neural compensatory mechanisms of BDNF val66met met carriers and APOE E4 carriers in healthy aging, mild cognitive impairment, and Alzheimer’s disease. Neurobiol. Aging 39, 165–173. doi: 10.1016/j.neurobiolaging.2015.12.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Gretarsdottir, S., Thorleifsson, G., Reynisdottir, S. T., Manolescu, A., Jonsdottir, S., Jonsdottir, T., et al. (2003). The gene encoding phosphodiesterase 4D confers risk of ischemic stroke. Nat. Genet. 35, 131–138. doi: 10.1038/ng1245

PubMed Abstract | CrossRef Full Text | Google Scholar

Guan, F., Zhang, C., Wei, S., Zhang, H., Gong, X., Feng, J., et al. (2012). Association of PDE4B polymorphisms and schizophrenia in Northwestern Han Chinese. Hum. Genet. 131, 1047–1056. doi: 10.1007/s00439-011-1120-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Gurney, M. E., D’Amato, E. C., and Burgin, A. B. (2015). Phosphodiesterase-4 (PDE4) molecular pharmacology and Alzheimer’s disease. Neurotherapeutics 12, 49–56. doi: 10.1007/s13311-014-0309-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Heckman, P. R., Wouters, C., and Prickaerts, J. (2015). Phosphodiesterase inhibitors as a target for cognition enhancement in aging and Alzheimer’s disease: a translational overview. Curr. Pharm. Des. 21, 317–331. doi: 10.2174/1381612820666140826114601

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoffmann, R., Wilkinson, I. R., McCallum, J. F., Engels, P., and Houslay, M. D. (1998). cAMP-specific phosphodiesterase HSPDE4D3 mutants which mimic activation and changes in rolipram inhibition triggered by protein kinase A phosphorylation of Ser-54: generation of a molecular model. Biochem. J. 333, 139–149. doi: 10.1042/bj3330139

PubMed Abstract | CrossRef Full Text | Google Scholar

Houslay, M. D., and Adams, D. R. (2003). PDE4 cAMP phosphodiesterases: modular enzymes that orchestrate signalling cross-talk, desensitization and compartmentalization. Biochem. J. 370, 1–18. doi: 10.1042/bj20021698

PubMed Abstract | CrossRef Full Text | Google Scholar

Jin, S. L., and Conti, M. (2002). Induction of the cyclic nucleotide phosphodiesterase PDE4B is essential for LPS-activated TNF-α responses. Proc. Natl. Acad. Sci. U S A 99, 7628–7633. doi: 10.1073/pnas.122041599

PubMed Abstract | CrossRef Full Text | Google Scholar

Jin, S. L., Lan, L., Zoudilova, M., and Conti, M. (2005). Specific role of phosphodiesterase 4B in lipopolysaccharide-induced signaling in mouse macrophages. J. Immunol. 175, 1523–1531. doi: 10.4049/jimmunol.175.3.1523

PubMed Abstract | CrossRef Full Text | Google Scholar

Lam, M., Trampush, J. W., Yu, J., Knowles, E., Davies, G., Liewald, D. C., et al. (2017). Large-scale cognitive GWAS meta-analysis reveals tissue-specific neural expression and potential nootropic drug targets. Cell Rep. 21, 2597–2613. doi: 10.1016/j.celrep.2017.11.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, H., Graham, J. M. Jr., Rimoin, D. L., Lachman, R. S., Krejci, P., Tompson, S. W., et al. (2012). Exome sequencing identifies PDE4D mutations in acrodysostosis. Am. J. Hum. Genet. 90, 746–751. doi: 10.1016/j.ajhg.2012.03.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, J. J., Wedow, R., Okbay, A., Kong, E., Maghzian, O., Zacher, M., et al. (2018). Gene discovery and polygenic prediction from a genome-wide association study of educational attainment in 1.1 million individuals. Nat. Genet. 50, 1112–1121. doi: 10.1038/s41588-018-0147-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Leibrock, J., Lottspeich, F., Hohn, A., Hofer, M., Hengerer, B., Masiakowski, P., et al. (1989). Molecular cloning and expression of brain-derived neurotrophic factor. Nature 341, 149–152. doi: 10.1038/341149a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Levin, L. R., Han, P. L., Hwang, P. M., Feinstein, P. G., Davis, R. L., and Reed, R. R. (1992). The Drosophila learning and memory gene rutabaga encodes a Ca2+/Calmodulin-responsive adenylyl cyclase. Cell 68, 479–489. doi: 10.1016/0092-8674(92)90185-f

PubMed Abstract | CrossRef Full Text | Google Scholar

Linglart, A., Fryssira, H., Hiort, O., Holterhus, P. M., Perez de Nanclares, G., Argente, J., et al. (2012). PRKAR1A and PDE4D mutations cause acrodysostosis but two distinct syndromes with or without GPCR-signaling hormone resistance. J. Clin. Endocrinol. Metab. 97, E2328–E2338. doi: 10.1210/jc.2012-2326

PubMed Abstract | CrossRef Full Text | Google Scholar

Michot, C., Le Goff, C., Goldenberg, A., Abhyankar, A., Klein, C., Kinning, E., et al. (2012). Exome sequencing identifies PDE4D mutations as another cause of acrodysostosis. Am. J. Hum. Genet. 90, 740–745. doi: 10.1016/j.ajhg.2012.03.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Millar, J. K., Pickard, B. S., Mackie, S., James, R., Christie, S., Buchanan, S. R., et al. (2005). DISC1 and PDE4B are interacting genetic factors in schizophrenia that regulate cAMP signaling. Science 310, 1187–1191. doi: 10.1126/science.1112915

PubMed Abstract | CrossRef Full Text | Google Scholar

Millar, J. K., Wilson-Annan, J. C., Anderson, S., Christie, S., Taylor, M. S., Semple, C. A., et al. (2000). Disruption of two novel genes by a translocation co-segregating with schizophrenia. Hum. Mol. Genet. 9, 1415–1423. doi: 10.1093/hmg/9.9.1415

PubMed Abstract | CrossRef Full Text | Google Scholar

Milner, B., Squire, L. R., and Kandel, E. R. (1998). Cognitive neuroscience and the study of memory. Neuron 20, 445–468. doi: 10.1016/S0896-6273(00)80987-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Motte, E., Le Stunff, C., Briet, C., Dumaz, N., and Silve, C. (2017). Modulation of signaling through GPCR-cAMP-PKA pathways by PDE4 depends on stimulus intensity: possible implications for the pathogenesis of acrodysostosis without hormone resistance. Mol. Cell. Endocrinol. 442, 1–11. doi: 10.1016/j.mce.2016.11.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Numata, S., Ueno, S., Iga, J., Song, H., Nakataki, M., Tayoshi, S., et al. (2008). Positive association of the PDE4B (phosphodiesterase 4B) gene with schizophrenia in the Japanese population. J. Psychiatr. Res. 43, 7–12. doi: 10.1016/j.jpsychires.2008.01.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Petrij, F., Giles, R. H., Dauwerse, H. G., Saris, J. J., Hennekam, R. C., Masuno, M., et al. (1995). Rubinstein-Taybi syndrome caused by mutations in the transcriptional co-activator CBP. Nature 376, 348–351. doi: 10.1038/376348a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Piovesan, A., Caracausi, M., Antonaros, F., Pelleri, M. C., and Vitale, L. (2016). GeneBase 1.1: a tool to summarize data from NCBI gene datasets and its application to an update of human gene statistics. Database 2016:28025344. doi: 10.1093/database/baw153

PubMed Abstract | CrossRef Full Text | Google Scholar

Richter, W., and Conti, M. (2002). Dimerization of the type 4 cAMP-specific phosphodiesterases is mediated by the upstream conserved regions (UCRs). J. Biol. Chem. 277, 40212–40221. doi: 10.1074/jbc.M203585200

PubMed Abstract | CrossRef Full Text | Google Scholar

Richter, W., Jin, S. L., and Conti, M. (2005). Splice variants of the cyclic nucleotide phosphodiesterase PDE4D are differentially expressed and regulated in rat tissue. Biochem. J. 388, 803–811. doi: 10.1042/bj20050030

PubMed Abstract | CrossRef Full Text | Google Scholar

Savage, J. E., Jansen, P. R., Stringer, S., Watanabe, K., Bryois, J., de Leeuw, C. A., et al. (2018). Genome-wide association meta-analysis in 269,867 individuals identifies new genetic and functional links to intelligence. Nat. Genet. 50, 912–919. doi: 10.1038/s41588-018-0152-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Sekar, A., Bialas, A. R., de Rivera, H., Davis, A., Hammond, T. R., Kamitaki, N., et al. (2016). Schizophrenia risk from complex variation of complement component 4. Nature 530, 177–183. doi: 10.1038/nature16549

PubMed Abstract | CrossRef Full Text | Google Scholar

Sette, C., and Conti, M. (1996). Phosphorylation and activation of a cAMP-specific phosphodiesterase by the cAMP-dependent protein kinase. Involvement of serine 54 in the enzyme activation. J. Biol. Chem. 271, 16526–16534. doi: 10.1074/jbc.271.28.16526

PubMed Abstract | CrossRef Full Text | Google Scholar

Tabuchi, A., Sakaya, H., Kisukeda, T., Fushiki, H., and Tsuda, M. (2002). Involvement of an upstream stimulatory factor as well as cAMP-responsive element-binding protein in the activation of brain-derived neurotrophic factor gene promoter I. J. Biol. Chem. 277, 35920–35931. doi: 10.1074/jbc.m204784200

PubMed Abstract | CrossRef Full Text | Google Scholar

Tanaka, Y., Naruse, I., Hongo, T., Xu, M., Nakahata, T., Maekawa, T., et al. (2000). Extensive brain hemorrhage and embryonic lethality in a mouse null mutant of CREB-binding protein. Mech. Dev. 95, 133–145. doi: 10.1016/s0925-4773(00)00360-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Trampush, J. W., Yang, M. L. Z., Yu, J., Knowles, E., Davies, G., Liewald, D. C., et al. (2017). GWAS meta-analysis reveals novel loci and genetic correlates for general cognitive function: a report from the COGENT consortium. Mol. Psychiatry 22, 1651–1652. doi: 10.1016/j.euroneuro.2016.09.587

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsai, S. J. (2018). Critical issues in BDNF Val66Met genetic studies of neuropsychiatric disorders. Front. Mol. Neurosci. 11:156. doi: 10.3389/fnmol.2018.00156

PubMed Abstract | CrossRef Full Text | Google Scholar

Turley, P., Walters, R. K., Maghzian, O., Okbay, A., Lee, J. J., Fontana, M. A., et al. (2018). Multi-trait analysis of genome-wide association summary statistics using MTAG. Nat. Genet. 50, 229–237. doi: 10.1038/s41588-017-0009-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Willoughby, D., Wong, W., Schaack, J., Scott, J. D., and Cooper, D. M. (2006). An anchored PKA and PDE4 complex regulates subplasmalemmal cAMP dynamics. EMBO J. 25, 2051–2061. doi: 10.1038/sj.emboj.7601113

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilson, N. M., Gurney, M. E., Dietrich, W. D., and Atkins, C. M. (2017). Therapeutic benefits of phosphodiesterase 4B inhibition after traumatic brain injury. PLoS One 12:e0178013. doi: 10.1371/journal.pone.0178013

PubMed Abstract | CrossRef Full Text | Google Scholar

Wong, S. T., Athos, J., Figueroa, X. A., Pineda, V. V., Schaefer, M. L., Chavkin, C. C., et al. (1999). Calcium-stimulated adenylyl cyclase activity is critical for hippocampus-dependent long-term memory and late phase LTP. Neuron 23, 787–798. doi: 10.1016/s0896-6273(01)80036-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, Y., Li, Z., Huang, Y. Y., Wu, D., and Luo, H. B. (2018). Novel phosphodiesterase inhibitors for cognitive improvement in Alzheimer’s disease. J. Med. Chem. 61, 5467–5483. doi: 10.1021/acs.jmedchem.7b01370

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, C., Zhang, D., Wu, Y., Tian, X., Pang, Z., Li, S., et al. (2017). A genome-wide association study of cognitive function in Chinese adult twins. Biogerontology 18, 811–819. doi: 10.1007/s10522-017-9725-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, J., Ferreira, T., Morris, A. P., Medland, S. E., Genetic Investigation of ANthropometric Traits (GIANT) Consortium, DIAbetes Genetics Replication and Meta-analysis (DIAGRAM) Consortium, Madden, et al. (2012). Conditional and joint multiple-SNP analysis of GWAS summary statistics identifies additional variants influencing complex traits. Nat. Genet. 44, 369–375, S1–S3. doi: 10.1038/ng.2213

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, J., Zhang, R., Yang, Z., Lee, J., Liu, Y., Tian, J., et al. (2013). Comparative effectiveness and safety of oral phosphodiesterase type 5 inhibitors for erectile dysfunction: a systematic review and network meta-analysis. Eur. Urol. 63, 902–912. doi: 10.1016/j.eururo.2013.01.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, C., Xu, Y., Chowdhary, A., Fox, D. III., Gurney, M. E., Zhang, H. T., et al. (2018). Memory enhancing effects of BPN14770, an allosteric inhibitor of phosphodiesterase-4D, in wild-type and humanized mice. Neuropsychopharmacology 43, 2299–2309. doi: 10.1038/s41386-018-0178-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, B., Yang, L., Konishi, Y., Maeda, N., Sakanaka, M., and Tanaka, J. (2002). Suppressive effects of phosphodiesterase type IV inhibitors on rat cultured microglial cells: comparison with other types of cAMP-elevating agents. Neuropharmacology 42, 262–269. doi: 10.1016/s0028-3908(01)00174-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, J., Mix, E., and Winblad, B. (2001). The antidepressant and antiinflammatory effects of rolipram in the central nervous system. CNS Drug Rev. 7, 387–398. doi: 10.1111/j.1527-3458.2001.tb00206.x

PubMed Abstract | CrossRef Full Text | Google Scholar