Identification of protease serine S1 family member 53 as a mitochondrial protein in murine islet beta cells.

The aim of this study was to identify genes that are specifically expressed in pancreatic islet ß-cells (hereafter referred to as ß-cells). Large-scale complementary DNA-sequencing analysis was performed for 3,429 expressed sequence tags derived from murine MIN6 ß-cells, through homology comparisons using the GenBank database. Three individual ESTs were found to code for protease serine S1 family member 53 (Prss53). Prss53 mRNA is processed into both a short and long form, which encode 482 and 552 amino acids, respectively. Transient overexpression of myc-tagged Prss53 in COS-7 cells showed that Prss53 was strongly associated with the luminal surfaces of organellar membranes and that it underwent signal peptide cleavage and N-glycosylation. Immunoelectron microscopy and western blotting revealed that Prss53 localized to mitochondria in MIN6 cells. Short hairpin RNA-mediated Prss53 knockdown resulted in Ppargc1a downregulation and Ucp2 and Glut2 upregulation. JC-1 staining revealed that the mitochondria were depolarized in Prss53-knockdown MIN6 cells; however, no change was observed in glucose-stimulated insulin secretion. Our results suggest that mitochondrial Prss53 expression plays an important role in maintaining the health of ß-cells.

TMEM16E (GDD1) exhibits protein instability and distinct characteristics in chloride channel/pore forming ability.

TMEM16E/GDD1 has been shown to be responsible for the bone-related late-onset disease gnathodiaphyseal dysplasia (GDD), with the dominant allele (TMEM16E(gdd) ) encoding a missense mutation at Cys356. Additionally, several recessive loss-of-function alleles of TMEM16E also cause late-onset limb girdle muscular dystrophy. In this study, we found that TMEM16E was rapidly degraded via the proteasome pathway, which was rescued by inhibition of the PI3K pathway and by the chemical chaperone, sodium butyrate. Moreover, TMEM16E(gdd) exhibited lower stability than TMEM16E, but showed similar propensity to be rescued. TMEM16E did not exhibit cell surface calcium-dependent chloride channel (CaCC) activity, which was originally identified in TMEM16A and TMEM16B, due to their intracellular vesicle distribution. A putative pore-forming domain of TMEM16E, which shared 39.8% similarity in 98 amino acids with TMEM16A, disrupted CaCC activity of TMEM16A via domain swapping. However, the Thr611Cys mutation in the swapped domain, which mimicked conserved cysteine residues between TMEM16A and TMEM16B, reconstituted CaCC activity. In addition, the GDD-causing cysteine mutation made in TMEM16A drastically altered CaCC activity. Based on these findings, TMEM16E possesses distinct function other than CaCC and another protein-stabilizing machinery toward the TMEM16E and TMEM16E(gdd) proteins should be considered for the on-set regulation of their phenotypes in tissues.

Novel compound heterozygous mutations of the growth hormone-releasing hormone receptor gene in a case of isolated growth hormone deficiency.

OBJECTIVE: To elucidate the pathogenesis of isolated growth hormone (GH) deficiency in a Japanese girl without consanguinity. DESIGN: A 2-year-old girl of height 77.2 cm (-3.0 SD for Japanese girls) was found to have an insulin-like growth factor (IGF)-1 level of 7 ng/mL and IGF binding protein-3 (IGFBP-3) level of 0.41 µg/mL. GH responded modestly to a series of pharmacological stimulants, increasing to 2.81 ng/mL with insulin-induced hypoglycemia, 3.78 ng/mL with arginine, and 3.93 with GH-releasing hormone (GHRH). Following direct sequencing of the GHRH receptor (GHRHR) gene, evaluation by the luciferase reporter assay, immunofluorescence study, and in vitro splicing assay with minigene constructs was conducted. RESULTS: Novel compound heterozygous GHRHR gene mutations were identified in the patient. A p.G136V substitution elicited no luciferase activity increment in response to GHRH stimulation, with normal membranous expression. Splicing assay demonstrated that the IVS2+3a>g mutation would lead to aberrant splicing. CONCLUSIONS: A case of isolated GH deficiency due to novel GHRHR gene mutations was identified.

Identification of INS and KCNJ11 gene mutations in type 1B diabetes in Japanese children with onset of diabetes before 5 years of age.

BACKGROUND: The etiology of type 1 diabetes (T1D) is heterogeneous and is according to presence or absence of pancreatic autoantibodies divided into two subtypes: type 1A (autoimmune-mediated) and type 1B (non-autoimmune-mediated). Although several genes have been linked to type 1A diabetes, the genetic cause of type 1B diabetes in Japanese individuals is far from understood. OBJECTIVE: The aim of this study was to test for monogenic forms of diabetes in auto antibody-negative Japanese children with T1D. METHODS: Thirty four (19 males and 15 female) unrelated Japanese children with glutamate decarboxylase (GAD) 65 antibodies and/or IA-2A-negative T1D and diabetes diagnosed at < 5 yr of age were recruited from 17 unrelated hospitals participating in the Japanese Study Group of Insulin Therapy for children and adolescent diabetes (JSGIT). We screened the INS gene and the KCNJ11 gene which encode the ATP-sensitive potassium cannel by direct sequencing in 34 Japanese children with T1D. RESULTS: We identified three novel (C31Y, C96R, and C109F) mutations and one previously reported mutation (R89C) in the INS gene in five children, in addition to one mutation in the KCNJ11 gene (H46R) in one child. These mutations are most likely pathogenic and therefore the cause of diabetes in carriers. CONCLUSION: Our results suggest that monogenic forms of diabetes, particularly INS gene mutations, can be detected in Japanese patients classified with type 1B. Mutation screening, at least of the INS gene, is recommended for Japanese patients diagnosed as autoantibody negative at <5 yr of age.

Bhlhe40, a potential diabetic modifier gene on Dbm1 locus, negatively controls myocyte fatty acid oxidation.

We have previously identified significant quantitative trait loci (QTL) Dbm1 (diabetic modifier QTL 1) on chromosome 6, affecting plasma glucose and insulin concentrations and body weight on F(2) progeny of hypoinsulinemic diabetic Akita mice, with the heterozygous Ins2 gene Cys96Tyr mutation, and non-diabetic A/J mice. To discover diabetic modifier genes on Dbm1, we constructed congenic strain for Dbm1 using the Akita allele as donor in A/J allele genetic background, and compared diabetes-related phenotypes to control mice. The homozygote for Akita allele of Dbm1 was associated with lower plasma glucose concentrations in glucose tolerance test (GTT) in the hypoinsulinemic condition derived from the Ins2 mutation and lower plasma insulin concentrations and body weight in the normoinsulinemic condition without the Ins2 mutation than the homozygote for A/J allele, as we performed QTL analysis on F(2) intercross mice. The Akita allele also decreased the epididymal white adipose tissue (EWAT) weight. According to the analysis of sub-congenic strains, we narrowed down the responsible diabetic modifier region to 9 Mb. As fourteen candidate genes exist in this region, we analyzed genomic variants of these genes and gene expression in the muscle, liver, and EWAT and identified that Bhlhe40 gene expression in muscle is decreased in congenic mice. According to the in vitro functional analyses, Bhlhe40 was shown to negatively control fatty acid oxidation in cultured myocyte. Based on these, we conclude that Bhlhe40 is a possible candidate diabetic modifier gene responsible for Dbm1 locus affecting diabetes and/or obesity through negatively controlling fatty acid oxidation in muscle.

Association Study of Fat-mass and Obesity-associated Gene and Body Mass Index in Japanese Patients with Schizophrenia and Healthy Subjects.

OBJECTIVE: Fat-mass and obesity-associated (FTO) gene is known to be involved in the pathophysiology of obesity and a single-nucleotide polymorphism (SNP) rs9939609 of FTO gene is repeatedly confirmed to be associated with body mass index (BMI) and obesity. The aim of this study is to elucidate effects of FTO gene polymorphism on BMI in Japanese patients with schizophrenia and healthy subjects. METHODS: Three hundred fifty one patients with schizophrenia and 342 age- and sex-matched healthy subjects participated in the study. Information on BMI and antipsychotic medication was also collected from patients and healthy subjects. Genotype of the FTO SNP rs9939609 was determined by TaqMan SNP Genotyping Assays. RESULTS: There was no significant difference in BMI between patients and healthy subjects. No significant difference in BMI was observed among any medications. We observed no significant difference in rs9939609 allele frequencies between patients and healthy subjects. There was a significant difference in BMI between healthy subjects with risk (AA or TA) genotypes and those with TT genotype. We also observed a significant positive correlation between the number of risk allele (A allele) and BMI in healthy subjects. CONCLUSION: Our study suggested that FTO rs9939609 polymorphism might have some impacts on the BMI in healthy subjects, but might not have same impacts on the BMI of patients with schizophrenia.

Diabetic modifier QTL, Dbm4, affecting elevated fasting blood glucose concentrations in congenic mice.

We have previously identified four significant quantitative trait loci (QTL) affecting plasma glucose concentrations on F(2) progeny of hypoinsulinemic diabetic Akita mice, heterozygous for the Ins2 gene Cys96Tyr mutation, and non-diabetic A/J mice, one of which on chromosome 15 named Dbm4 (diabetic modifier QTL 4) was shown to affect fasting plasma glucose concentrations with a maximum LOD score of 6.17. To estimate the influence of Dbm4 itself to the diabetes-related phenotypes, we constructed congenic strain with heterozygous Ins2 mutation using the Akita allele as donor of Dbm4 locus in the A/J genetic background, and measured quantitative traits including plasma glucose concentrations in glucose tolerance test (GTT). In this study, we found the Akita allele in Dbm4 was associated with higher fasting plasma glucose concentrations as in previous QTL analysis. According to gene expression assay, key enzymes of hepatic gluconeogenesis were expressed to the more increased degree in the liver of congenic mice compared to the A/J allele based control mice. Based on these results, we concluded that diabetic modifier gene(s) exist on Dbm4 locus affecting fasting plasma glucose concentrations via regulation of gluconeogenic gene expression in the hypoinsulinemic diabetic condition. Identification of the modifier gene responsible for Dbm4 would provide new drug development targets for human type 2 diabetes with hepatic insulin resistance.

Identification of a novel mutation in the exon 2 splice donor site of the POU1F1/PIT-1 gene in Japanese identical twins with mild combined pituitary hormone deficiency.

CONTEXT: To date, approximately 35 different POU1F1 mutations have been described in patients with familial and sporadic combined pituitary hormone deficiency (CPHD) from different ethnic backgrounds. The majority are missense mutations clustered within the conserved POU-specific and POU-homeo domains, encoded by exons 4 and 6, respectively. OBJECTIVES: This study aimed to identify the molecular basis and clinical characteristics of a Japanese CPHD family with a novel POU1F1 mutation. DESIGN: The POU1F1 gene was sequenced in identical twin brothers with mild CPHD. The mutation identified was also evaluated in family members as well as 188 Japanese controls and then examined in functional studies. RESULTS: A novel heterozygous splice site mutation (Ex2 + 1G>T; c.214 + 1G>T) was detected. This mutation was also present in their undiagnosed mother, but not in any of the controls. In vitro splicing studies suggested this mutation to result in an in-frame skipping of exon 2, thus producing an internally deleted protein lacking most of the R2 transactivation subdomain (TAD-R2). Heterologous expression studies of the mutated POU1F1 protein showed only modest reductions in its transactivation activities in HEK293T cells, while acting as a dominant-negative inhibitor of the endogenous activities of POU1F1 in pituitary GH3 cells. CONCLUSIONS: This is the first report of a mutation at the exon 2 donor splice site of POU1F1, affecting TAD-R2. The addition of this mutation to the growing list of pathological POU1F1 mutations may provide deeper insights into clinical heterogeneity in the expressions of individual mutations and a better understanding of the structure-function relationships of POU1F1.

Analysis of expression and structure of the rat GH-secretagogue/ghrelin receptor (Ghsr) gene: roles of epigenetic modifications in transcriptional regulation.

In the current study, to elucidate the molecular basis of cell type-specific expression of the GH-secretagogue/ghrelin receptor type 1A (GHSR1A), we characterized the structure and putative promoter region of the rat Ghsr gene. We identified an alternative 5'-untranslated first exon that contains multiple transcription start sites, and confirmed a 200-bp sequence proximal to this exon to be sufficient for basal promoter activity. A promoter-associated CpG island conserved across different species was found to be hypomethylated in Ghsr1a-expressing cell lines, while being heavily methylated in non-expressing cells. In cells with low or absent Ghsr1a expression, treatment with demethylating agents activated Ghsr1a transcription. Chromatin immunoprecipitation assays demonstrated Ghsr1a-expressing cells to display active histone modifications, whereas repressive modifications were present exclusively in other cell types. These results suggest epigenetic modifications at GHSR to play important roles in determining GHSR1A expression and abundance, and therefore the consequent sensitivity of cells to ghrelin.

Identification and functional analysis of novel human growth hormone secretagogue receptor (GHSR) gene mutations in Japanese subjects with short stature.

CONTEXT: Short stature (SS) is a multifactorial developmental condition with a significant genetic component. Recent studies have revealed that rare deleterious mutations in the GH-secretagogue receptor type 1A (GHSR1A) gene could be a cause of familial SS or GH deficiency. OBJECTIVE: The aim of this study was to evaluate the contribution of GHSR1A mutations to the molecular mechanism underlying SS in Japanese subjects. METHODS: We performed mutational screening of the GHSR1A gene in 127 unrelated Japanese SS patients diagnosed with either isolated GH deficiency or idiopathic SS. Identified mutations were analyzed in 188 control subjects, and their functional properties were examined in a heterologous expression system. RESULTS: Four novel heterozygous GHSR1A mutations were identified (ΔQ36, P108L, C173R, and D246A). Expression studies demonstrated that these mutations had varying functional consequences: 1) all mutations showed a loss-of-function effect on the constitutive signaling activity of GHSR1A, but the degree of loss varied widely; 2) C173R caused intracellular retention of the mutated protein, resulting in total loss of receptor function; 3) P108L resulted in a large decrease in binding affinity to ghrelin, without affecting its surface expression; 4) D246A uniquely impaired agonist- and inverse agonist-stimulated receptor signaling; and 5) ΔQ36 showed only a subtle reduction in constitutive activity. The cumulative frequency of these putative functional mutations was significantly higher in the patient group than in controls (4.72 vs. 0.53%; P = 0.019; odds ratio = 9.28; 95% confidence interval, 1.10-78.0). CONCLUSIONS: Our results suggest that GHSR1A mutations contribute to the genetic etiology of SS in the Japanese population.

Identification and functional analysis of novel human growth hormone-releasing hormone receptor (GHRHR) gene mutations in Japanese subjects with short stature.

CONTEXT: Growth hormone-releasing hormone receptor (GHRHR) gene mutations have been identified in patients of different ethnic origins with isolated GH deficiency (IGHD) type IB. However, the prevalence of these mutations in the Japanese population has yet to be fully determined. OBJECTIVES: This study aimed to evaluate the contributions of GHRHR mutations to the molecular mechanism underlying short stature in Japanese subjects. DESIGN: The GHRHR gene was sequenced in 127 unrelated Japanese patients with either IGHD (n = 14) or idiopathic short stature (ISS; n = 113). Sequence variants were evaluated in family members and 188 controls, and then examined in functional studies. RESULTS: A novel homozygous E382E (c.1146G>A) synonymous variant, at the last base of exon 12, was identified in an IGHD family with two affected sisters. In vitro splicing studies showed this mutation to result in skipping of exon 12. In one ISS patient, a heterozygous ATG-166T>C variant was found in the distal Pit-1 P2 binding element of the GHRHR promoter. In two control subjects, a close but distinct variant, ATG-164T>C, was detected. Functional studies showed that both promoter variants diminish promoter activity by altering Pit-1 binding ability. Four missense variants were also found in both patient and control groups but had no detectable functional consequences. CONCLUSIONS: The homozygous GHRHR mutation was rare, being detected in only one Japanese IGHD family. Future research is needed to clarify the genetic contributions of heterozygous functional promoter variants to GHD, ISS and normal-stature variations.

Recessive mutations in the putative calcium-activated chloride channel Anoctamin 5 cause proximal LGMD2L and distal MMD3 muscular dystrophies.

The recently described human anion channel Anoctamin (ANO) protein family comprises at least ten members, many of which have been shown to correspond to calcium-activated chloride channels. To date, the only reported human mutations in this family of genes are dominant mutations in ANO5 (TMEM16E, GDD1) in the rare skeletal disorder gnathodiaphyseal dysplasia. We have identified recessive mutations in ANO5 that result in a proximal limb-girdle muscular dystrophy (LGMD2L) in three French Canadian families and in a distal non-dysferlin Miyoshi myopathy (MMD3) in Dutch and Finnish families. These mutations consist of a splice site, one base pair duplication shared by French Canadian and Dutch cases, and two missense mutations. The splice site and the duplication mutations introduce premature-termination codons and consequently trigger nonsense-mediated mRNA decay, suggesting an underlining loss-of-function mechanism. The LGMD2L phenotype is characterized by proximal weakness, with prominent asymmetrical quadriceps femoris and biceps brachii atrophy. The MMD3 phenotype is associated with distal weakness, of calf muscles in particular. With the use of electron microscopy, multifocal sarcolemmal lesions were observed in both phenotypes. The phenotypic heterogeneity associated with ANO5 mutations is reminiscent of that observed with Dysferlin (DYSF) mutations that can cause both LGMD2B and Miyoshi myopathy (MMD1). In one MMD3-affected individual, defective membrane repair was documented on fibroblasts by membrane-resealing ability assays, as observed in dysferlinopathies. Though the function of the ANO5 protein is still unknown, its putative calcium-activated chloride channel function may lead to important insights into the role of deficient skeletal muscle membrane repair in muscular dystrophies.

Diabetes mellitus in a Japanese girl with HDR syndrome and GATA3 mutation.

We report on a Japanese girl with HDR (hypoparathyroidism, sensorineural deafness, and renal dysplasia) syndrome who developed diabetes mellitus (DM) at three years of age (blood glucose 713 mg/dL, HbA(1c) 8.0%) in the absence of anti-glutamic acid decarboxylase autoantibodies. Mutation analysis revealed a de novo heterozygous two base pair deletion at exon 6 of the GATA3 gene (c.1200_1201delCA; p.H400fsX506). GATA3 expression was identified by PCR amplification for human pancreas cDNA, and mouse Gata3 was weekly but unequivocally expressed in pancreatic beta cells. The results, in conjunction with the previous findings indicating the critical role of GATA3 in lymphocyte function, GATA3 haploinsufficiency may affect the function of beta cells and/or lymphocytes, leading to the development of DM in relatively exceptional patients with high susceptibility to DM.

Association study between the pericentrin (PCNT) gene and schizophrenia.

Disrupted-in-schizophrenia 1 (DISC1), a known genetic risk factor for schizophrenia (SZ) and major depressive disorder (MDD), interacts with several proteins and some of them are reported to be genetically associated with SZ. Pericentrin (PCNT) also interacts with DISC1 and recently single-nucleotide polymorphisms (SNPs) within the PCNT gene have been found to show significant associations with SZ and MDD. In this study, case-controlled association analysis was performed to determine if the PCNT gene is implicated in SZ. Nine SNPs were analyzed in 1,477 individuals (726 patients with SZ and 751 healthy controls). No significant difference was observed between the controls and the patients in allelic frequencies or genotypic distributions of eight SNPs. Although allelic distribution of rs11702684 was different between the two groups (P = 0.042), the difference did not reach statistical significance after permutation correction for multiple comparisons. In the haplotypic analysis, we could not find any significant association in our subjects, either. This gene may not play a major role independently in the etiology of SZ in the Japanese population.

YKL-40 secreted from adipose tissue inhibits degradation of type I collagen.

Obesity is considered a chronic low-grade inflammatory status and the stromal vascular fraction (SVF) cells of adipose tissue (AT) are considered a source of inflammation-related molecules. We identified YKL-40 as a major protein secreted from SVF cells in human visceral AT. YKL-40 expression levels in SVF cells from visceral AT were higher than in those from subcutaneous AT. Immunofluorescence staining revealed that YKL-40 was exclusively expressed in macrophages among SVF cells. YKL-40 purified from SVF cells inhibited the degradation of type I collagen, a major extracellular matrix of AT, by matrix metalloproteinase (MMP)-1 and increased rate of fibril formation of type I collagen. The expression of MMP-1 in preadipocytes and macrophages was enhanced by interaction between these cells. These results suggest that macrophage/preadipocyte interaction enhances degradation of type I collagen in AT, meanwhile, YKL-40 secreted from macrophages infiltrating into AT inhibits the type I collagen degradation.

Positive association of the pericentrin (PCNT) gene with major depressive disorder in the Japanese population.

BACKGROUND: Pericentrin (PCNT) interacts with disruption-in-schizophrenia 1 (DISC1), a known genetic risk factor for schizophrenia, bipolar disorder and major depressive disorder (MDD). We sought to determine whether the PCNT gene is implicated in MDD. METHODS: We performed case-control association analyses in the Japanese population. We analyzed 9 single nucleotide polymorphisms (SNPs) in 173 patients with MDD and 348 healthy controls. RESULTS: We found a significant allelic association between 3 SNPs (rs3788265, rs2073376 and rs2073380) of the PCNT gene and MDD (p = 0.006, 0.005 and 0.021, respectively). After correction for multiple testing, 2 SNPs (rs3788265 and rs2073376) retained significant allelic associations with MDD. In addition, we found a significant association between the 2 marker haplotypes (r3788265 and rs2073376) and MDD (permutation p = 0.011). LIMITATIONS: Our sample was small and comprised only Japanese participants. In addition, owing to the late onset of MDD, it is possible that the disorder will develop in at least some participants in our control group. Finally, we did not show how SNPs of the PCNT gene alter its function. CONCLUSION: Our results suggest that genetic variations in the PCNT gene may play a significant role in the etiology of MDD in the Japanese population.

Overexpression of HMGA2 relates to reduction of the let-7 and its relationship to clinicopathological features in pituitary adenomas.

High-mobility group A2 is highly expressed during embryogenesis and in various benign and malignant tumors. Recent studies report that high-mobility group A2 is negatively regulated by the let-7 microRNAs (miRNAs) family in vitro. The development of pituitary adenomas in high-mobility group A2 transgenic mice showed that high-mobility group A2 may be involved in pituitary tumorigenesis. However, no studies have investigated the clinical significance of high-mobility group A2 and its relationship to the let-7 miRNA family in human pituitary adenomas. Using immunohistochemistry, we analyzed high-mobility group A2 expression with respect to various clinicopathologic factors in 98 pituitary adenomas. Overexpression of high-mobility group A2 was observed in 39% (38/98) of pituitary adenomas compared with normal adenohypophysial tissue and was frequently found in adenomas including prolactin (PRL), adrenocorticotrophic hormone, or follicle-stimulating hormone/luteinizing hormone and in null cell adenomas, but relatively rare in growth hormone (GH) and mixed GH/PRL adenomas. High-mobility group A2 expression was significantly associated with tumor invasion (P<0.05) and was significantly higher in grade IV than in grades I, II, and III adenomas (P<0.05). High levels of high-mobility group A2 expression were more frequently observed in macroadenomas than in microadenomas (P<0.05). High levels of high-mobility group A2 expression also significantly correlated with the proliferation marker Ki-67 (P<0.0001). Real-time quantitative RT-PCR analysis was carried out to evaluate the expression of let-7 in 55 pituitary adenomas. Subsequently, decreased expression of let-7 was confirmed in 23 of 55 (42%) adenomas and was correlated with high-grade tumors (P<0.05). An inverse correlation between let-7 and high-mobility group A2 expression was evident (R=-0.33, P<0.05). These findings support a causal link between let-7 and high-mobility group A2 whereby loss of let-7 expression induces high-mobility group A2 upregulation that represents an important mechanism in pituitary tumorigenesis and progression.

Gene expression and association analyses of the phosphodiesterase 4B (PDE4B) gene in major depressive disorder in the Japanese population.

The phosphodiesterase 4B (PDE4B) interacts with disrupted-in-schizophrenia 1 (DISC1), which is a known genetic risk factor for schizophrenia, bipolar disorder and major depressive disorder (MDD). PDE4B is also important in the regulation of cAMP signaling, a second messenger implicated in learning, memory, and mood. In this study, we determined mRNA expression levels of the PDE4B gene in the peripheral blood leukocytes of patients with MDD and control subjects (n = 33, each). Next we performed two-stage case-controlled association analyses (first set; case = 174, controls = 348; second set; case = 481, controls = 812) in the Japanese population to determine if the PDE4B gene is implicated in MDD. In the leukocytes, a significantly higher expression of the PDE4B mRNA was observed in the drug-naïve MDD patients compared with control subjects (P < 0.0001) and the expression of the MDD patients significantly decreased after antidepressant treatment (P = 0.030). In the association analysis, we observed significant allelic associations of four SNPs (the most significant, rs472952; P = 0.002) and a significant haplotypic association (permutation P = 0.019) between the PDE4B gene and MDD in the first-set samples. However, we could not confirm these significant associations in the following independent second-set of samples. Our results suggest that the PDE4B gene itself does not link to MDD but the elevated mRNA levels of PDE4B might be implicated in the pathophysiology of MDD.