- Research article
- Open Access
- Open Peer Review
The obesity gene, TMEM18, is of ancient origin, found in majority of neuronal cells in all major brain regions and associated with obesity in severely obese children
BMC Medical Geneticsvolume 11, Article number: 58 (2010)
TMEM18 is a hypothalamic gene that has recently been linked to obesity and BMI in genome wide association studies. However, the functional properties of TMEM18 are obscure.
The evolutionary history of TMEM18 was inferred using phylogenetic and bioinformatic methods. The gene's expression profile was investigated with real-time PCR in a panel of rat and mouse tissues and with immunohistochemistry in the mouse brain. Also, gene expression changes were analyzed in three feeding-related mouse models: food deprivation, reward and diet-induced increase in body weight. Finally, we genotyped 502 severely obese and 527 healthy Swedish children for two SNPs near TMEM18 (rs6548238 and rs756131).
TMEM18 was found to be remarkably conserved and present in species that diverged from the human lineage over 1500 million years ago. The TMEM18 gene was widely expressed and detected in the majority of cells in all major brain regions, but was more abundant in neurons than other cell types. We found no significant changes in the hypothalamic and brainstem expression in the feeding-related mouse models. There was a strong association for two SNPs (rs6548238 and rs756131) of the TMEM18 locus with an increased risk for obesity (p = 0.001 and p = 0.002).
We conclude that TMEM18 is involved in both adult and childhood obesity. It is one of the most conserved human obesity genes and it is found in the majority of all brain sites, including the hypothalamus and the brain stem, but it is not regulated in these regions in classical energy homeostatic models.
Obesity is an increasing global health problem accompanied by several significant health impairments, such as type-2 diabetes and cardiovascular diseases. Although excessive energy intake together with low levels of physical activity contribute to the ongoing obesity epidemic, genetic research has shown that heritage might account for as much as 70% of the population variation in BMI . Genome-wide association (GWA) studies have linked several genomic loci with a high BMI, and several new candidate genes have been identified . Most recently, the FTO gene has been given special attention because of its strong association with obesity in multiple cohorts differing in age . Basic research showed that FTO-deficient mice have less adipose tissue and a lower body mass compared with wild type mice, which supports the biological role proposed for the gene by the GWA analyses .
Aside from FTO, GWA studies associated with obesity yet another gene, named TMEM18. TMEM18 was first proposed as an important obesity-related locus by the GIANT consortium . That study strongly associated the nearby SNP, rs6548238, to increased BMI and body weight. Subsequently, these results were replicated in the Icelandic GWA study for three other SNPs located in the close proximity of TMEM18 (rs2867125, rs4854344 and rs7561317), which all showed a strong association to both BMI and body weight . The FTO locus' association to BMI was confirmed therein. Furthermore, the TMEM18 locus' association to BMI has been replicated in two additional studies [7, 8] and it has also been associated to the age of menarche . It should be noted that while the risk contribution of the TMEM18 locus to obesity in adults has been examined (6), no attempt to determine such association has been made in obese children. However, one previous study has associated the locus with BMI in children .
The TMEM18 gene codes for a poorly characterized transmembrane protein. One study indicated that this protein is located in the nuclear envelope in neural stem cells . TMEM18 may be involved in cell migration as overexpression of the protein increases the migration of neural stem cells towards glioma in the rat brain. A preliminary expression profiling of TMEM18, that accompanied the first GWA study, suggested that it is ubiquitously expressed, but with certain differences between tissues . In line with the proposed role in the regulation of body weight, TMEM18 was shown to be expressed in the brain, including the hypothalamus, the region responsible for the control of energy homeostasis. Only few other brain areas were investigated in that project.
The current project provides the first detailed characterization of the human TMEM18 protein's sequence features, evolutionary history, gene expression profile in the rat and mouse, distribution in the mouse brain, and gene expression regulation in several murine feeding/body weight models. We also investigated the risk of obesity for variants upstream of the TMEM18 locus in the cohort of severely obese children and healthy controls, and performed the first examination of the locus association to several obesity-related traits.
Evolutionary and sequence analysis
Retrieval of protein sequences
The proteomes of eleven different organisms representing different lineages of eukaryotes, with the focus on the metazoans and vertebrates, were gathered from different sources as follows. All protein sequences for Homo sapiens, Mus musculus, Gallus gallus, Xenopus tropicalis, Caenorhabditis elegans and Drosophila melanogaster were downloaded from Ensembl 53 ; Branciostoma floridae ,Trichoplax adhaerens  and Thalassiosira pseudonana  were downloaded from the Joint Genome Institute's website, Saccharomyces cerevisae was downloaded from the Saccharomyces genome database , Arabidopsis thaliana from The Arabidopsis Information Resource (TAIR)  and all bacterial proteomes were downloaded from NCBI's ftp-site .
Mining for TMEM18
The human TMEM18 sequence was collected from UniProt  and used to query the OMA database  (release September 2008). The OMA is a database which automatically groups protein sequences predicted to be orthologous. To allow identification of members of the TMEM18 family that are not orthologous to the human gene, we used the following procedure. Sequences that belong to the orthologous group of human TMEM18 were downloaded from the database and aligned using mafft-einsi  with default parameters. The multiple sequence alignment (MSA) was viewed and edited in Jalview . A hidden-markov model (HMM) was created and calibrated from the MSA using the HMMER  package with default parameters. The HMM was applied on the listed proteomes using the hmmpfam program of the HMMER package with the default parameters. The results were evaluated for false-positives and hits with an E-value below 10-6 were considered as significant. The protein and nucleotide sequences for the hits were collected and redundancy was manually removed so that only one protein/transcript was kept for each gene.
Multiple sequence alignments and Phylogenetic analysis
The protein sequences identified in the mining procedure were aligned using mafft-einsi with default parameters. The resulting MSA was viewed, analyzed and edited in Jalview. The protein sequences of the MSA were reversely translated to the corresponding mRNA transcripts by using RevTrans  with the default parameters and the retrieved nucleotide MSA was used for the phylogenetic analysis. The phylogenetic inference was performed with the maximum likelihood method using RAxML (v. 7.0.4)  according to the "Hard & Slow Way" described in the program's manual, the EasyRAx script found at the program's website was utilized to execute the analysis. In brief, the best known maximum likelihood tree was found in 100 maximum parsimony trees, 1000 non-parametric bootstraps were carried out and frequencies were written on the best known maximum likelihood tree. The final tree was drawn using Mesquite (v 2.6) and edited in Inkscape (v. 0.46), using the T. pseudonana sequence as the outgroup. Pairwise sequence alignments were made using EMBOSS'  implementation of the Needleman-Wunsch  algorithm with the default parameters.
Sequence analysis of the human TMEM18 protein
Poly Phobius  (v. 1.04) was used to predict transmembrane topology for the human TMEM18 sequence. The protein MSA made for the evolutionary analysis was used as input to increase accuracy of the prediction. Glycosylation sites of type O-β-GlcNAc, which occurs in the nucleus, was predicted using the YinOYang  (v. 1.2) server with the default parameters. Phosphorylation sites were predicted with the NetPhos  (v. 2.0) server using the default parameters. The prediction of nuclear localization signals (NLS) was performed with PredictNLS and the NLSdb database . Finally, the UniProt and Ensembl websites were consulted for additional annotation.
TMEM18 tissue expression profiling with quantitative real-time PCR
The TMEM18 mRNA expression was detected in tissue panels from both the mouse and the rat. Animal handling and tissue panel preparation was performed as previously described in Haitina et al. together with details of the real-time PCR analysis . The primers used can be found in Additional file 1. All animal experiments has been approved by "Uppsala's Ethical Committee on Animal Experiments" (reference number: C285/5 and C84/8) and follows international guidelines.
TMEM18-like immunoreactivity in the murine brain
C57Bl6/J adult male mice were housed in the controlled environment (21°C, 12:12 LD cycle), with ad libitum access to chow and water. An intraperitonial injection of mixture of Dormitor (70 μg/g b. wt., Orion, Finland) and Ketalar (7 μg/g b. wt., Pfizer, Sweden) was used to anesthetize the mice. Phosphate buffered saline (PBS) followed by 4% paraformaldehyde were used to perform transcardial perfusion through the left ventricle. The fixed brain was excised and incubated overnight in 4% paraformaldehyde; subsequently the tissue was dehydrated and infiltrated with paraffin (Tissue Tek vacuum infiltration processor; Miles Scientific, Elkhart, IN). Sectioning (7 μm) of the paraffin-embedded brain was carried out on the Microm microtome onto superfrost slides (Menzel-Gläser, Braunschweig, Germany).
Sections were deparaffinized in X-tra solve (Medite Histotechnic, Burgdorf, Germany) and rehydrated with a series of washes in descending concentrations of ethanol followed by autoclaved MilliQ water. Antigen retrieval was performed by boiling the sections in 0.01 M citric acid (pH 6.0) for 10 minutes. Slides were cooled and rinsed with 1× PBS pH 7.4. They were transferred to the humidity (PBS) chamber. Primary antibodies: rabbit anti-TMEM18, (ProSci, Poway, CA) and mouse anti-NeuN (Neuronal Nuclei, Millipore, Temecula CA, USA) diluted 1:400 in PBS+0.3% Triton X100 (Sigma-Aldrich, USA) were added onto the slides and incubation was carried out overnight at 4°C. The slides were then washed with PBS and incubated in the secondary antibodies: donkey-anti-rabbit 594 (1:200; Alexa Fluor, Invitrogen, USA) and goat-anti-mouse 488 (1:400; Alexa Fluor, Invitrogen, USA), respectively in PBS+0.3% Triton X100 in the dark chamber for 4 hours. The slides were washed with PBS and stained with DAPI. They were coverslipped using DTG (2.5% DABCO (Sigma), 50 mM Tris-HCl pH 8.0, 90% glycerol). Sections were analyzed under the Zeiss XBO75 microscope.
Analysis of results
Pictures of the cortex, amygdala, hypothalamus, thalamus and hippocampus were taken with the Zeiss AxioCam HRm under 20× magnification. Photomicrographs show a random portion of each region of interest. The images were analyzed and cells counted through appropriate filters using ImageJ 1.41o software (National Institute of Health, USA). A total of 5 pictures from the cortex, 4 from the amygdala, and 3 each from the hypothalamus, thalamus and hippocampus were analyzed. The total number of all cells, NeuN cells and TMEM18-positive cells as well as the number of TMEM18 expressing neurons and non-neurons were counted separately. Percentages of neurons expressing TMEM18, non-neurons expressing TMEM18, and all cells combined expressing TMEM18 were calculated.
Male C57BL/6J mice (Scanbur, Sweden), housed individually or in groups of two (Exp. 1, 2 and 3) in macrolon cages, under LD 12:12 (lights on at 0700), weighed ca. 28 g at the beginning of the experiment. Tap water and chow (Lactamin, Sweden) were available ad libitum unless specified otherwise. Procedures described herein were approved by the Uppsala Animal Welfare Committee and followed the EU and Swedish guidelines.
Experiment 1. Hypothalamic and brainstem TMEM18 mRNA levels following 16-h and 24-h of food deprivation
Chow was removed before the onset of the dark phase and mice (n = 8) were decapitated between 1000 and 1100 on the next day. Control mice (n = 8) had ad libitum access to food.
Experiment 2. Hypothalamic and brainstem TMEM18 mRNA levels following 48-h consumption of palatable sucrose or Intralipid
Mice gained access to a bottle containing 10% sucrose or 4.1% Intralipid (Fresenius, Sweden) for 48 h; control animals had access to chow only (n = 8/group). Intralipid, a palatable lipid emulsion, has been used in experiments utilizing liquid diets [32, 33]. Sucrose and Intralipid were isocaloric (0.4 kcal/g, 1 kcal = 4.2 kJ); energy density of chow was 3.6 kcal/g. The solutions were similar in palatability: each mouse ingested on average 8.1 kcal of Intralipid and 7.4 kcal of sucrose per day. Total energy intake per animal was 10.3 kcal in the chow group, 14.1 kcal in the sucrose group and 12.9 kcal in the Intralipid group. Mice were decapitated after 48 h (1100 - 1200).
Experiment 3. Hypothalamic and brainstem TMEM18 mRNA levels upon increased body weight
Mice received 10% sucrose in addition to chow for 3 weeks; controls had chow only (n = 8/group). The initial body weights did not differ significantly between the two groups (controls: 27.6 ± 0.3 g; sucrose: 27.8 ± 0.5 g). At endpoint, sucrose-fed mice weighed 32.1 ± 0.4 g, whereas the control animals weighed 29.5 ± 0.6 g (P < 0.05; t-test). Energy intake was similar to that described for the corresponding groups in Experiment 1. Mice were decapitated between 1100 and 1200.
Analysis of genetic association to obesity and related traits
Birth weights and lengths were obtained from growth charts. Actual body weights and lengths were measured to the nearest 0.1 kg and 1 cm, respectively. BMI standard deviation score (BMI SDS) was calculated from weight, height, standardized for age and gender . For the obese subjects, levels of plasma glucose, serum insulin, serum triglycerides and cholesterol were analysed from blood samples drawn after 16-h overnight fasting as previously described .
We genotyped 1027 children and adolescents comprising two study groups as described earlier . Briefly, a case group of 502 obese children (262 girls and 240 boys) was enrolled at the National Childhood Obesity Centre at Karolinska University Hospital, Huddinge, Sweden. The obese subjects were between 6 and 20 years old and mean age was 12.6 ± 3.3 years. BMI SDS ranged from 2.7 to 11.1 and mean BMI SDS was 6.2 ± 1.4 2) The control group with 525 healthy Swedish adolescents (268 girls and 257 boys) was recruited from 17 upper secondary schools in the Stockholm area. The subjects in the control group were between 15 and 20 years old and mean age was 17.1 ± 0.8 years. BMI SDS ranged from -1.9 to 2.4 and mean BMI SDS was 0.2 ± 0.8. Subjects with overweight/obesity or chronic diseases were excluded from the control group and seven children with type-2 diabetes were excluded from the obese group. The study was approved by the Regional Committee of Ethics, Stockholm, (reference number KI702/03) and follows the Helsinki declaration and international guidelines. All participants or their legal guardians gave their written informed consent of inclusion in the study and publication of the results.
In order to test for deviation from Hardy-Weinberg equilibrium, the Pearson's χ2-test (1 d.f) was applied. Genotype and allele frequencies were calculated and logistic regression was used to calculate the odds ratio (OR) with a 95% confidence interval (CI) assuming an additive model. Association with obesity was determined by comparing subjects with normal weight and obesity. Associations between genotypes and phenotypes were analyzed with linear regression, assuming an additive model. Quantitative skewed variables were log-transformed before analysis. Covariates such as birth weight, birth length, weight, length, gender, BMI SDS and age were tested for dependence on the response variables and included in the model if significant. All the analysis was performed using PLINK .
Genomic DNA was extracted from peripheral blood using QiaGen Maxiprep kit (Qiagen, Hilden, Germany). Genotyping was performed with pre-designed Taqman single-nucleotide polymorphism genotyping assay (Applied Biosystems, Foster City, USA) and an ABI7900 genetic analyser with SDS 2.2 software.
Sequence analysis of the human TMEM18 protein
A comprehensive picture of the human TMEM18 protein can be found in Figure 1. The human TMEM18 gene is located on chromosome 2. It is 9466 bp-long, including introns, and it has seven protein coding exons. The gene has three transcripts of which the longest one is translated into a 140 amino acid-long protein, that were predicted to have three transmembrane α-helices. The N-terminus of the protein, which probably faces the inside of the nuclear envelope, was predicted to have five O-β-GlcNAc glycosylation sites, and the second loop to have one phosphorylation site. The protein has a nuclear localization signal at the C-terminal, which according to NLSdb, is identical to, among others, human homeobox proteins and it targets the protein for transport to the nucleus. A coiled-coil domain was also found in the same part of the protein.
We searched for homologues in a wide range of genomes in all kingdoms of life and we present herein the specific mining of eleven genomes for genes homologous to TMEM18. We found that all genomes contained one single positive hit for the TMEM18 family, except for S. cerevisae and C. elegans where no homologous genes were identified. The non-redundant proteins sequences used that were received in the mining is available in Additional file 2. We investigated this further by mining the NCBI's protein database nr using protein PSI-BLAST with the human TMEM18 protein sequence as query and with the default parameters. We found no sequences from the fungi or pseudoocoelomata (where C. elegans belongs) lineages in these searches, suggesting that these taxonomic groups have lost the TMEM18 gene. The multiple sequence alignment (Figure 2) showed that the TMEM18 sequence is well conserved, eleven residues were identical in all the investigated organisms and another 30 were identical with 80% of the sequences. The average pair wise sequence identity (similarity) was 37% (82%) for all species and 41% (88%) within vertebrates. Thus, the TMEM18 family has been highly conserved during the evolutionary history of eukaryotes.
We created a phylogenetic tree for the TMEM18 family (Figure 3). T. pseudonana was used as the outgroup and, as expected; the plant A. thaliana was found to be basal of the metazoan group. Within the metazoan group, the order of invertebrates was not as anticipated, since T. adhaerens should be most basal, followed by D. melanogaster and B. floridae. However, the supporting bootstrap values were not strong for this part of the topology, which may reflect high conservation of the sequences. The vertebrates formed a clear cluster with an expected topology, apart from the reversal of the X. tropicalis and D. rerio positions, but here the bootstrap values were low.
We investigated the TMEM18's mRNA expression pattern in a wide range of rat and mouse tissues using real-time PCR. We found that TMEM18 was widely expressed, but varied in expression level between different tissues and was abundant in feeding and body weight regulatory brain regions such as the hypothalamus and brainstem, although higher in other tissues (Figure 4). As a comparison, we downloaded microarray expression data for the D. melanogaster TMEM18 from FlyAtlas , which confirmed the wide expression pattern found in the rat and mouse (Figure 4)
Immunohistochemical detection of TMEM18 in the mouse brain
We stained a large number of sections from the mouse brain. No region was devoid of TMEM18 immunoreactivity and the TMEM18 protein showed a widespread distribution throughout the mouse brain. We also performed co-labeling of TMEM18 with the cellular marker, DAPI, and a neuronal marker, NeuN (Figure 5) and in total 7276 cells from the cortex (1904 cells), amygdala (1505 cells), hypothalamus (1314 cells), thalamus (1033 cells) and hippocampus (1520 cells) were counted. The results revealed that TMEM18 was present in 71% of all cell types, 73% of NeuN-positive cells and 63% of NeuN-negative cells (Figure 5).
Central expression of TMEM18 in mice differing in feeding or body weight profiles
Hypothalamic and brainstem expression of TMEM18 was unchanged in all of the feeding paradigms utilized herein, including 16- and 24-h energy deprivation, as well as short- (48 h) and long-term (3 weeks long; leading to increased body weight) exposure to palatable tastants that differ in macronutrient composition (see Figure 6).
Analysis of genetic association to obesity and related traits
None of the two SNPs analyzed caused deviation from Hardy-Weinberg equilibrium (Table 1) and the variants were almost in complete linkage disequilibrium, r2 = 0.95. Table 1 shows obesity risk estimates (odds ratios; OR) for the two genotyped variants located 23 kbp upstream TMEM18, rs6548238 and rs7561317. Both SNPs were significantly (p = 0.002 and p = 0.001) associated with obesity in directions consistent with prior reports of association with BMI (OR = 1.450 and OR = 1.507) [5–7]. Table 2 summarizes the analyzed phenotypes in the normal weight and obese children and adolescents for the two SNPs. No significant association was found after adjustment for multiple testing. There was, however, a clear trend for BMI SDS among the obese children: the homozygote carriers of the major allele had 0.5 SD higher BMI SDS levels compared to the homozygote carriers of the minor allele (p = 0.073 and p = 0.038). This trend was also observed among the obese children for a shorter birth length for carriers of the major allele for both SNPs (p = 0.048 and p = 0.032).
The TMEM18 gene has a remarkably long evolutionary history as it is found in most of the eukaryotic genomes we have looked at. It has therefore been present for at least 1500 MYA since the divergence of animals and plants/T. pseudonana . Thus, it is one of the most ancient genes implicated in human obesity, such as FTO  and MC4R , that we are aware of. Most of the main genes involved in body weight regulation appeared after the radiation of vertebrates and very few are present beyond the animal kingdom. The amino acid sequences are highly conserved (Figures 2 and 3) with relatively small differences between vertebrates and other animals. Another very interesting feature (of curiosity, shared by the FTO gene), is that the gene is found in a single copy with no close relative in the genomes where it is present, despite that some genomes have undergone recent whole genome duplications such as A. thaliana (38 MYA) . Surprisingly, we found two separate lineages where the TMEM18 gene is absent: fungi (yeasts, molds and mushrooms) and pseudoocoelomata (C. elegans). Both contain well analyzed genomes of high quality, suggesting that this gene has been lost at least two times on separate occasions in evolution and that the gene is thus not essential for life in all eukaryotes. The TMEM18 protein is predicted to have three transmembrane helices, a relatively uncommon property among human membrane proteins, as fewer than 5% of them share this feature . The membrane topology is mainly found among proteins with poorly characterized functions that do not have close relatives and are, therefore, sole representatives of their protein family in the human proteome. The TMEM18 protein has a nuclear localization signal that targets the protein for transportation to the nucleus, which also has been confirmed experimentally . Moreover, it contains a coiled-coil domain (Figure 1) that is common among DNA or RNA binding proteins, such a transcription factors, e.g. c-Fos . However, transcription factors are generally not membrane-bound proteins, but at least one exception exists in the gene named Myelin gene regulatory factor (MRF) that is predicted to have a transmembrane helix and plays a central regulatory role in myelination of the CNS . We found that there are only seven proteins, unrelated to each other, in the human proteome that contain three transmembrane helices and one coiled-coil domain, according to our queries in UniProt. Thus, the molecular structure of the TMEM18 protein is very uncommon in terms of other proteins in the human genome and surely very much different compared to other proteins involved in obesity, such as the enzyme FTO and the G protein-coupled melanocortin receptor, MC4R.
Other interesting features of the TMEM18 protein are the evolutionary conserved potential O-β-GlcNAc glycosylation and the phosphorylation sites, which are targets for the important posttranslational regulation of many proteins' functions and activities. Noteworthy, these sites are found in parts of the proteins that are supposed to face the inside of the nuclear envelope, i.e., the N-terminal and second loop (Figure 1) which, together with the beginning of the first and second transmembrane helices, are the most evolutionary conserved parts of the TMEM18 protein. This suggests that these parts of the protein constitute a potential active site, binding site or another important functional structure of the protein.
The TMEM18 gene is expressed in all the tissues we analyzed in the rat and mouse. Moreover, microarray data on TMEM18 from the fly (D. melanogaster), found in the FlyAtlas database, (Figure 4) correspond well to the mammalian expression data and confirm the widespread expression with the abundant presence in the brain, although higher in several other tissues. Thus, the expression pattern is in agreement with the evolutionary conservation of the protein sequence. We performed a comprehensive immunohistochemical analysis as well as in situ hybridization (Additional file 3). The immunohistochemical results showed that the TMEM18 protein is expressed throughout the mouse brain, which agrees with the real-time PCR expression profiles (Figure 4) and the in situ hybridization (See Additional file 3). We performed co-localization with specific neuronal marker and assessed the extent of the expression in five major brain areas. We found that TMEM18 is expressed in most cells (71%), with a somewhat higher, but significant (p = 0.011; t-test), presence in neurons (73%) than in non-neurons (64%). The outcome is consistent for different brain areas studied (Figure 5), while the hypothalamus is the only region that has a relatively lower percentage of expression in neurons. The TMEM18 gene was presented as a hypothalamic gene by the earlier study reporting the GWA to obesity . We did not however see any evidence of significant enrichment in the TMEM18 expression level in the hypothalamus compared to other central regions, which is clear for, e.g., the aforementioned FTO gene and most neuroendocrine genes involved in body weight regulation. Thus, our results show that TMEM18 is likely to have a function in multiple tissues rather than being a solely hypothalamic gene. Nonetheless, the real-time PCR and immunohistochemical analyses confirmed that the TMEM18 is expressed abundantly in the hypothalamus and brainstem, regions that play a crucial role in the regulation of energy homeostasis: the hypothalamus, as the area that controls primarily feeding for calories and contains also components of the feeding reward system, and the brainstem that serves as the "relay station" for the gut-CNS signal cascade. Changes in energy needs of the organism, body weight as well as the rewarding value of food affect expression of genes known to participate in the regulation of consumption and body weight. Therefore, we studied the potential regulation of the TMEM18 gene expression in three different animal models that explore classical feeding paradigms. Feeding-related gene expression profile has previously been shown to be significantly changed in these tissues. For example, palatability affected hypothalamic expression of melanocortin receptor-4 (MC4R) and kappa (KOR) and mu (MOR) opioid receptors; an increased body weight influenced proopiomelanocortin (POMC), melanin concentrating hormone (MCH) and dynorphin (DYN)  while restriction of energy intake through deprivation upregulated expression of the obesity-associated gene FTO. Somewhat surprisingly, the TMEM18's expression level in the hypothalamus and brainstem remained the same in all models utilized herein (Figure 6). In fact, it shows a more stable expression (GeNorm expression stability value; M = 0.493) than the commonly used housekeeping genes: Tubulin beta (M = 0.575) in the brainstem samples from the deprivation and increased body weight models (Experiment 1 and 3) and Cyclophilin (M = 0.505, TMEM18 M = 0.458) in the hypothalamus samples from the Sucrose vs Intralipid model (Experiment 2). The average expression stability value is 0.469, which is remarkable as it parallels housekeeping genes in these models. Still, TMEM18 would not be considered a true housekeeping gene as it is not found in all cells. It is however clear that TMEM18 is differently regulated compared to the classical neuropeptides and receptors, whose mRNA levels are altered in these models. It is interesting in this context that TMEM18 contains several potential phosphorylation and glycosylation sites (Figure 1) as both these processes are common regulatory pathways for nuclear proteins. It is also possible that this gene is so prevalent in central neurons (73% of them express it) that feeding-related changes in TMEM18 expression might have affected only a small subpopulation of TMEM18 cells and they were not detectable with the real-time PCR analysis of a large region, such as the hypothalamus or brainstem.
Our human genetic study showed a strong significance for the association of two SNPs (Table 1) located upstream of TMEM18 with childhood obesity in a Swedish cohort of obese children and a clear trend for an increase in BMI SDS for homozygote carriers of the major allele. Further, the estimations of the odds-ratios suggest that there are significantly higher risks for carriers of the major allele of both SNPs to develop obesity (Table 1). The odds-ratio for rs6548238 (OR = 1.450) is in the same magnitude as the calculations made by Renström and colleagues (OR = 1.100) for a cohort of 4923 Swedish adults . Moreover, the odds-ratios for the locus near TMEM18 are also as high as reported for the FTO locus (OR = 1.593) for the same cohort . This strengthens previous observations that this region is as strongly associated with obesity as FTO. We do not find however any significant association for the variations with any of the investigated obesity-related traits, many relevant for the detection of a diabetic phenotype, while there is a trend for a decrease in birth weight and length for homozygote carriers of the major alleles (Table 2). In conclusion, the results show that a region near TMEM18 is significantly contributing to obesity in severely obese children, a group that arguably has the highest genetic component and greatest need for healthcare among different population subgroups of obese individuals.
In conclusion, we show that common variations nearby the TMEM18 gene are associated not only to adult obesity but also obesity in severely obese children in a similar magnitude as FTO, the first gene associated with the common form of obesity. TMEM18 has a remarkably long evolutionary history that spans at least 1500 MY, longer than most other genes implicated in body weight regulation that we are aware of. The gene is widely distributed in the brain, found in the majority of all cells in the brain with slight, yet significant, enrichment in neurons compared to non-neurons. Expression of the gene does not seem to be affected in classical homeostatic models where several known regulators of appetite and reward are regulated. The protein has highly evolutionary conserved structural properties that induce the targeting to the nucleus of the cell, which may clarify why this newly identified gene has different expression and regulation patterns compared to the known neuroendocrine factors thus far associated with human obesity.
Maes HH, Neale MC, Eaves LJ: Genetic and environmental factors in relative body weight and human adiposity. Behavior genetics. 1997, 27: 325-351. 10.1023/A:1025635913927.
Hofker M, Wijmenga C: A supersized list of obesity genes. Nature genetics. 2009, 41: 139-140. 10.1038/ng0209-139.
Frayling TM, Timpson NJ, Weedon MN, Zeggini E, Freathy RM, Lindgren CM, Perry JR, Elliott KS, Lango H, Rayner NW, Shields B, Harries LW, Barrett JC, Ellard S, Groves CJ, Knight B, Patch AM, Ness AR, Ebrahim S, Lawlor DA, Ring SM, Ben-Shlomo Y, Jarvelin MR, Sovio U, Bennett AJ, Melzer D, Ferrucci L, Loos RJ, Barroso I, Wareham NJ, et al: A common variant in the FTO gene is associated with body mass index and predisposes to childhood and adult obesity. Science (New York, NY). 2007, 316: 889-894.
Fischer J, Koch L, Emmerling C, Vierkotten J, Peters T, Bruning JC, Ruther U: Inactivation of the Fto gene protects from obesity. Nature. 2009, 458: 894-898. 10.1038/nature07848.
Willer CJ, Speliotes EK, Loos RJ, Li S, Lindgren CM, Heid IM, Berndt SI, Elliott AL, Jackson AU, Lamina C, Lettre G, Lim N, Lyon HN, McCarroll SA, Papadakis K, Qi L, Randall JC, Roccasecca RM, Sanna S, Scheet P, Weedon MN, Wheeler E, Zhao JH, Jacobs LC, Prokopenko I, Soranzo N, Tanaka T, Timpson NJ, Almgren P, Bennett A, et al: Six new loci associated with body mass index highlight a neuronal influence on body weight regulation. Nature genetics. 2009, 41: 25-34. 10.1038/ng.287.
Thorleifsson G, Walters GB, Gudbjartsson DF, Steinthorsdottir V, Sulem P, Helgadottir A, Styrkarsdottir U, Gretarsdottir S, Thorlacius S, Jonsdottir I, Jonsdottir T, Olafsdottir EJ, Olafsdottir GH, Jonsson T, Jonsson F, Borch-Johnsen K, Hansen T, Andersen G, Jorgensen T, Lauritzen T, Aben KK, Verbeek AL, Roeleveld N, Kampman E, Yanek LR, Becker LC, Tryggvadottir L, Rafnar T, Becker DM, Gulcher J, et al: Genome-wide association yields new sequence variants at seven loci that associate with measures of obesity. Nature genetics. 2009, 41: 18-24. 10.1038/ng.274.
Renstrom F, Payne F, Nordstrom A, Brito EC, Rolandsson O, Hallmans G, Barroso I, Nordstrom P, Franks PW: Replication and extension of genome-wide association study results for obesity in 4923 adults from northern Sweden. Human molecular genetics. 2009, 18: 1489-1496. 10.1093/hmg/ddp041.
Zhao J, Bradfield JP, Li M, Wang K, Zhang H, Kim CE, Annaiah K, Glessner JT, Thomas K, Garris M, Frackelton EC, Otieno FG, Shaner JL, Smith RM, Chiavacci RM, Berkowitz RI, Hakonarson H, Grant SF: The Role of Obesity-associated Loci Identified in Genome-wide Association Studies in the Determination of Pediatric BMI. Obesity Silver Spring, Md. 2009, 17 (12): 2254-7. 10.1038/oby.2009.159.
Perry JR, Stolk L, Franceschini N, Lunetta KL, Zhai G, McArdle PF, Smith AV, Aspelund T, Bandinelli S, Boerwinkle E, Cherkas L, Eiriksdottir G, Estrada K, Ferrucci L, Folsom AR, Garcia M, Gudnason V, Hofman A, Karasik D, Kiel DP, Launer LJ, van Meurs J, Nalls MA, Rivadeneira F, Shuldiner AR, Singleton A, Soranzo N, Tanaka T, Visser JA, Weedon MN, et al: Meta-analysis of genome-wide association data identifies two loci influencing age at menarche. Nature genetics. 2009,
Jurvansuu J, Zhao Y, Leung DS, Boulaire J, Yu YH, Ahmed S, Wang S: Transmembrane protein 18 enhances the tropism of neural stem cells for glioma cells. Cancer research. 2008, 68: 4614-4622. 10.1158/0008-5472.CAN-07-5291.
Hubbard TJ, Aken BL, Ayling S, Ballester B, Beal K, Bragin E, Brent S, Chen Y, Clapham P, Clarke L, Coates G, Fairley S, Fitzgerald S, Fernandez-Banet J, Gordon L, Graf S, Haider S, Hammond M, Holland R, Howe K, Jenkinson A, Johnson N, Kahari A, Keefe D, Keenan S, Kinsella R, Kokocinski F, Kulesha E, Lawson D, Longden I, et al: Ensembl 2009. Nucleic acids research. 2009, 37: D690-697. 10.1093/nar/gkn828.
Putnam NH, Butts T, Ferrier DE, Furlong RF, Hellsten U, Kawashima T, Robinson-Rechavi M, Shoguchi E, Terry A, Yu JK, Benito-Gutierrez EL, Dubchak I, Garcia-Fernandez J, Gibson-Brown JJ, Grigoriev IV, Horton AC, de Jong PJ, Jurka J, Kapitonov VV, Kohara Y, Kuroki Y, Lindquist E, Lucas S, Osoegawa K, Pennacchio LA, Salamov AA, Satou Y, Sauka-Spengler T, Schmutz J, Shin IT, et al: The amphioxus genome and the evolution of the chordate karyotype. Nature. 2008, 453: 1064-1071. 10.1038/nature06967.
Srivastava M, Begovic E, Chapman J, Putnam NH, Hellsten U, Kawashima T, Kuo A, Mitros T, Salamov A, Carpenter ML, Signorovitch AY, Moreno MA, Kamm K, Grimwood J, Schmutz J, Shapiro H, Grigoriev IV, Buss LW, Schierwater B, Dellaporta SL, Rokhsar DS: The Trichoplax genome and the nature of placozoans. Nature. 2008, 454: 955-960. 10.1038/nature07191.
Armbrust EV, Berges JA, Bowler C, Green BR, Martinez D, Putnam NH, Zhou S, Allen AE, Apt KE, Bechner M, Brzezinski MA, Chaal BK, Chiovitti A, Davis AK, Demarest MS, Detter JC, Glavina T, Goodstein D, Hadi MZ, Hellsten U, Hildebrand M, Jenkins BD, Jurka J, Kapitonov VV, Kroger N, Lau WW, Lane TW, Larimer FW, Lippmeier JC, Lucas S, et al: The genome of the diatom Thalassiosira pseudonana: ecology, evolution, and metabolism. Science (New York, NY). 2004, 306: 79-86.
Saccharomyces genome database. [http://www.yeastgenome.org]
Huala E, Dickerman AW, Garcia-Hernandez M, Weems D, Reiser L, LaFond F, Hanley D, Kiphart D, Zhuang M, Huang W, Mueller LA, Bhattacharyya D, Bhaya D, Sobral BW, Beavis W, Meinke DW, Town CD, Somerville C, Rhee SY: The Arabidopsis Information Resource (TAIR): a comprehensive database and web-based information retrieval, analysis, and visualization system for a model plant. Nucleic acids research. 2001, 29: 102-105. 10.1093/nar/29.1.102.
NCBI ftp-page. [ftp://ftp.ncbi.nlm.nih.gov]
Bairoch A, Bougueleret L, Altairac S, Amendolia V, Auchincloss A, Argoud-Puy G, Axelsen K, Baratin D, Blatter MC, Boeckmann B, Bolleman J, Bollondi L, Boutet E, Quintaje SB, Breuza L, Bridge A, deCastro E, Ciapina L, Coral D, Coudert E, Cusin I, Delbard G, Dornevil D, Roggli PD, Duvaud S, Estreicher A, Famiglietti L, Feuermann M, Gehant S, Farriol-Mathis N, et al: The Universal Protein Resource (UniProt) 2009. Nucleic acids research. 2009, 37: D169-174. 10.1093/nar/gkn664.
Roth AC, Gonnet GH, Dessimoz C: Algorithm of OMA for large-scale orthology inference. BMC bioinformatics. 2008, 9: 518-10.1186/1471-2105-9-518.
Katoh K, Toh H: Recent developments in the MAFFT multiple sequence alignment program. Briefings in bioinformatics. 2008, 9: 286-298. 10.1093/bib/bbn013.
Waterhouse AM, Procter JB, Martin DM, Clamp M, Barton GJ: Jalview Version 2--a multiple sequence alignment editor and analysis workbench. Bioinformatics (Oxford, England). 2009, 25: 1189-1191. 10.1093/bioinformatics/btp033.
Eddy SR: Profile hidden Markov models. Bioinformatics (Oxford, England). 1998, 14: 755-763. 10.1093/bioinformatics/14.9.755.
Wernersson R, Pedersen AG: RevTrans: Multiple alignment of coding DNA from aligned amino acid sequences. Nucleic acids research. 2003, 31: 3537-3539. 10.1093/nar/gkg609.
Stamatakis A: RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics (Oxford, England). 2006, 22: 2688-2690. 10.1093/bioinformatics/btl446.
Rice P, Longden I, Bleasby A: EMBOSS: the European Molecular Biology Open Software Suite. Trends Genet. 2000, 16: 276-277. 10.1016/S0168-9525(00)02024-2.
Needleman SB, Wunsch CD: A general method applicable to the search for similarities in the amino acid sequence of two proteins. Journal of molecular biology. 1970, 48: 443-453. 10.1016/0022-2836(70)90057-4.
Kall L, Krogh A, Sonnhammer EL: An HMM posterior decoder for sequence feature prediction that includes homology information. Bioinformatics (Oxford, England). 2005, 21 (Suppl 1): i251-257. 10.1093/bioinformatics/bti1014.
Gupta R, Brunak S: Prediction of glycosylation across the human proteome and the correlation to protein function. Pacific Symposium on Biocomputing. 2002, 310-322.
Blom N, Gammeltoft S, Brunak S: Sequence and structure-based prediction of eukaryotic protein phosphorylation sites. Journal of molecular biology. 1999, 294: 1351-1362. 10.1006/jmbi.1999.3310.
Nair R, Carter P, Rost B: NLSdb: database of nuclear localization signals. Nucleic acids research. 2003, 31: 397-399. 10.1093/nar/gkg001.
Haitina T, Olsson F, Stephansson O, Alsio J, Roman E, Ebendal T, Schioth HB, Fredriksson R: Expression profile of the entire family of Adhesion G protein-coupled receptors in mouse and rat. BMC neuroscience. 2008, 9: 43-10.1186/1471-2202-9-43.
Miedlar JA, Rinaman L, Vollmer RR, Amico JA: Oxytocin gene deletion mice overconsume palatable sucrose solution but not palatable lipid emulsions. Am J Physiol Regul Integr Comp Physiol. 2007, 293: R1063-1068.
Sclafani A, Rinaman L, Vollmer RR, Amico JA: Oxytocin knockout mice demonstrate enhanced intake of sweet and nonsweet carbohydrate solutions. Am J Physiol Regul Integr Comp Physiol. 2007, 292: R1828-1833.
Rolland-Cachera MF, Sempe M, Guilloud-Bataille M, Patois E, Pequignot-Guggenbuhl F, Fautrad V: Adiposity indices in children. The American journal of clinical nutrition. 1982, 36: 178-184.
Jacobsson JA, Danielsson P, Svensson V, Klovins J, Gyllensten U, Marcus C, Schioth HB, Fredriksson R: Major gender difference in association of FTO gene variant among severely obese children with obesity and obesity related phenotypes. Biochemical and biophysical research communications. 2008, 368: 476-482. 10.1016/j.bbrc.2008.01.087.
Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, Sham PC: PLINK: a tool set for whole-genome association and population-based linkage analyses. American journal of human genetics. 2007, 81: 559-575. 10.1086/519795.
Chintapalli VR, Wang J, Dow JA: Using FlyAtlas to identify better Drosophila melanogaster models of human disease. Nature genetics. 2007, 39: 715-720. 10.1038/ng2049.
Hedges SB, Dudley J, Kumar S: TimeTree: a public knowledge-base of divergence times among organisms. Bioinformatics (Oxford, England). 2006, 22: 2971-2972. 10.1093/bioinformatics/btl505.
Fredriksson R, Hagglund M, Olszewski PK, Stephansson O, Jacobsson JA, Olszewska AM, Levine AS, Lindblom J, Schioth HB: The obesity gene, FTO, is of ancient origin, up-regulated during food deprivation and expressed in neurons of feeding-related nuclei of the brain. Endocrinology. 2008, 149: 2062-2071. 10.1210/en.2007-1457.
Haitina T, Klovins J, Takahashi A, Lowgren M, Ringholm A, Enberg J, Kawauchi H, Larson ET, Fredriksson R, Schioth HB: Functional characterization of two melanocortin (MC) receptors in lamprey showing orthology to the MC1 and MC4 receptor subtypes. BMC evolutionary biology. 2007, 7: 101-10.1186/1471-2148-7-101.
Ermolaeva MD, Wu M, Eisen JA, Salzberg SL: The age of the Arabidopsis thaliana genome duplication. Plant molecular biology. 2003, 51: 859-866. 10.1023/A:1023001130337.
Almen MS, Nordstrom KJ, Fredriksson R, Schioth HB: Mapping the human membrane proteome: a majority of the human membrane proteins can be classified according to function and evolutionary origin. BMC biology. 2009, 7: 50-10.1186/1741-7007-7-50.
Glover JN, Harrison SC: Crystal structure of the heterodimeric bZIP transcription factor c-Fos-c-Jun bound to DNA. Nature. 1995, 373: 257-261. 10.1038/373257a0.
Emery B, Agalliu D, Cahoy JD, Watkins TA, Dugas JC, Mulinyawe SB, Ibrahim A, Ligon KL, Rowitch DH, Barres BA: Myelin gene regulatory factor is a critical transcriptional regulator required for CNS myelination. Cell. 2009, 138: 172-185. 10.1016/j.cell.2009.04.031.
Olszewski PK, Fredriksson R, Olszewska AM, Stephansson O, Alsio J, Radomska KJ, Levine AS, Schioth HB: Hypthalamic FTO is associated with the regulation of energy intake. BMC neuroscience. 2009, 10: 129-10.1186/1471-2202-10-129.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/11/58/prepub
The studies were supported by the Swedish Research Council (VR, medicine), AFA insurance, Swedish Brain Research Foundation, Svenska Läkaresällskapet, Åhlens Foundation, Novo Nordisk Foundation, National Institute of Drug Abuse (R01DA021280), and National Institute of Diabetes and Ingestive and Kidney Diseases (P30DK50456).
The authors declare that they have no competing interests.
MSA designed, conceived and performed the bioinformatical and evolutionary analysis and drafted the manuscript. JAJ carried out the genetic studies and drafted the manuscript. JHAS did the immunohistochemical and in situ hybdridization experiments. PKO designed, performed and analyzed the feeding experiments and participated in the writing of the draft. JC and JA designed, performed and analyzed the expression studies and contributed to the discussion of the manuscript. ASL designed and conceived the feeding experiments. RF, MC and HBS conceived the study, and participated in its design and the writing of the draft. All authors read and approved the final manuscript.
Markus Sällman Almén, Josefin A Jacobsson contributed equally to this work.