- Research Paper
- Open access
- Published:
Olfactory receptor genes cooperate with protocadherin genes in human extreme obesity
Genes & Nutrition volume 10, Article number: 16 (2015)
Abstract
Worldwide, the incidence of obesity has increased dramatically over the past decades. More knowledge about the complex etiology of obesity is needed in order to find additional approaches for treatment and prevention. Investigating the exome sequencing data of 30 extremely obese subjects (BMI 45–65 kg/m2) shows that predicted damaging missense variants in olfactory receptor genes on chromosome 1q and rare predicted damaging variants in the protocadherin (PCDH) beta-cluster genes on chromosome 5q31, reported in our previous work, co-localize in subjects with extreme obesity. This implies a synergistic effect between genetic variation in these gene clusters in the predisposition to extreme obesity. Evidence for a general involvement of the olfactory transduction pathway on itself could not be found. Bioinformatic analysis indicates a specific involvement of the PCDH beta-cluster genes in controlling tissue development. Further mechanistic insight needs to await the identification of the ligands of the 1q olfactory receptors. Eventually, this may provide the possibility to manipulate food flavor in a way to reduce the risk of overeating and of extreme obesity in genetically predisposed subjects.
Introduction
Worldwide, the incidence of obesity has increased dramatically over the past decades and all kinds of measures ranging from dietary and lifestyle interventions to bariatric surgery are being taken to revert this situation, albeit with limited success. Apparently, more knowledge about the complex etiology of obesity is needed in order to find additional approaches for treatment and prevention. Obesity and parameters of adiposity are to a considerable extent determined by genetic factors. Heritability estimates for obesity in different studies amount to at least 0.4 (Walley et al. 2006; Farooqi and O’Rahilly 2007; Berndt et al. 2013). Knowledge of those factors and the processes that they influence provides possibilities for diagnostic testing and counseling in addition to development of novel intervention methods. Unfortunately, more than 90 % of the genetic background is still unexplained (Speliotes et al. 2010; Loos 2012; Wheeler et al. 2013); thus, a further search for obesity-related genes is warranted. Recently, using functional clustering analysis of genes with rare variants, we have shown that in people with extreme obesity (BMI 45–65 kg/m2), the protocadherin (PCDH) genes on chromosome 5q31 harbor a surplus of rare variation with a moderate-to-high predicted biological effect (Mariman et al. 2014).
Studies in the mouse have provided evidence for a prominent role of the PCDH gene cluster in the development of the olfactory system (Hasegawa et al. 2012; Ledderose et al. 2013; Hirano et al. 2012). We therefore wondered whether the previously identified rare PCDH gene variation (Mariman et al. 2014) and damaging variation in the genes of the olfactory transduction pathway might work synergistically in the predisposition to extreme obesity. In that case, damaging missense variants in the genes of the olfactory transduction pathway might occur more often in the extreme obese subjects carrying the rare potentially damaging variation in the PCDH genes than in the non-carrier subjects. In the present study, we analyzed genetic variants of the olfactory system in relation to the earlier reported variants in the PCDH genes in 30 extreme obese subjects.
Materials and methods
Subjects
The selection of the subjects has been described before (Mariman et al. 2014). In short, from 561 obese individuals who attended the private obesity clinic CO-EUR (http://www.co-eur.eu/), 30 (19 females, 11 males) relatively young (average age 29.7 years, range 19–40.4 years) extremely obese (average BMI 51.1 kg/m2, range 45.3–65.1 kg/m2) Caucasian subjects were selected. Written informed consent for genetic studies was obtained from all participants, and permission was granted by the Medical Research Ethics Committee of the Maastricht University Medical Center. For these subjects, only weight, length, age and sex were registered, and no clinical biochemical data were available.
DNA isolation and sequencing
From peripheral blood leukocytes, genomic DNA was isolated using the QIAamp DNA blood kit (Qiagen, Amsterdam, The Netherlands) and was then outsourced for whole exome sequencing in a CLIA-certified laboratory (EdgeBio, http://www.edgebio.com/). The Nimblegen capture kit was used for the selection of genomic coding regions, followed by sequencing on the Illumina HiSeq 2000.
For our exome sequencing, the target 64 Mb consisted of 34-Mb exons. The average median whole targeted sequencing read depth was 50. For all samples, 94.6 % of bases had a read depth >10, and 87.2 % of bases had a read depth >20. 84–86 % of the reads mapped to the target regions, and 14–16 % was off target.
Sequence results were compared with GRCh37 (b37) as the reference genome (http://www.ncbi.nlm.nih.gov/projects/genome/assembly/grc/human/). For every subject, data were returned in excel-readable format as a SNP file and an indel file (insertion/deletion). Only the SNP files were used for the present study. Information was provided concerning genetic data (chromosome, cytoband, reference position, gene, zygosity, rs number, alleles, etc.), quality of the sequencing (Phred quality score, depth) and predicted biological effect by various methods (SNPeffect (Reumers et al. 2005), SIFT score (http://sift.jcvi.org/), PolyPhen2 class and score (http://genetics.bwh.harvard.edu/pph2/), etc.). In addition, the corresponding KEGG pathway identifier(s) (http://www.genome.jp/kegg/) were given for the genes.
Identification of PCDH-variant carriers
Previously, we analyzed the sequencing data for the presence of rare variants with a moderate-to-high predicted biological damaging effect in functional gene clusters (Mariman et al. 2014). PCDH gene clusters on chromosome 5q31 displayed a high frequency of such variants in the extreme obese cohort. In the present study, we compared genetic data between carriers of these PCDH variants and the individuals not carrying such variants in our cohort (non-carriers). In Table 1, PCDH-variant carriers are listed with information of the PCDH genes in which their rare variants occurred.
Selection of variants
For the present study, we focused on missense variants with a predicted damaging impact on protein function present in genes that had already been functionally annotated in KEGG pathways. The selection route for these variants is depicted in Fig. 1. In short, all missense variants detected in the 30 extremely obese subjects with a Phred quality score >30 (equivalent to a base call accuracy >99.9 %) and a read depth of at least 20, with KEGG pathway ID(s) and with scores on biological impact of the variant allele were selected. From the resulting list, a subset of predicted damaging variants was further selected based on the presence of both a SIFT score of 0.0 and a PolyPhen2 class ‘D’. The SIFT score is based on the premise that the amino acid order of a protein is related to evolution. A SIFT score below 0.05 means that the amino acid is highly conserved and substitution may result in damaging of the function of the protein. In this regard, a SIFT score of 0.0 indicates a high probability of damaging (Ng and Henikoff 2003). PolyPhen2 takes structural features of the protein into account to classify amino acid substitutions like occurrence in a functional module and/or impact on the 3D structure (Adzhubei et al. 2010). Substitutions are quantitatively scored, but also qualitatively classified as ‘benign’ (B), ‘possibly damaging’ (P) and ‘probably damaging’ (D). We took the last category into account.
Thus, two SNP data files were composed: file A: all annotated missense variants and file B: predicted damaging missense variants. The amount of damaging missense variants were evaluated by two parameters: (1) the ratio of damaging missense variants to total missense variants, which equals the number of missense variants per pathway in file B/number of missense variants of the same pathway in file A and (2) the burden, which equals the sum of damaging missense variants per pathway per person.
Computational assessment of the function of PCDH gene clusters
The genes of the PCDH alpha-, beta- and gamma-clusters were used to infer coexpression data. Publically available gene expression data (GEO http://www.ncbi.nlm.nih.gov/geo/) for human and mouse were used to find genes that coexpress with PCDH clusters. The approach is an extension of the one published by Baughman et al. (2009). Shortly, using an in-house computational pipeline, all microarray experiment consisting of multiple expression measurements was evaluated separately. Evaluation consists of calculating the average Pearson’s coexpression of each gene of the three PCDH clusters. If the coexpression of the system (e.g., alpha-cluster) is congruent (on average, genes are up- and downregulated together), a high weight is assigned to the experiment. Analogously, if genes within the system do not appear to coexpress with each other in the experiment, the experiment itself receives a low weight. After assigning weights to experiments, for each experiment, the coexpression of each human gene is calculated with each of the clusters as average Pearson’s correlation of expression. The results of coexpression with the clusters are then integrated for all experiments (in a weight-sensitive manner, as previously calculated) and for both human and mouse. This gives the final list of genes coexpressing with each cluster for which the top 100 was selected as genes most likely functionally related to PCDH gene sets. As can be expected, those PCDH cluster-linked gene sets contained several PCDH genes. For further analysis, those genes were removed resulting in a set of 83 genes for the alpha-cluster, 70 genes for the beta-cluster and 65 genes for the gamma-cluster.
Next, these gene sets were analyzed by the program DAVID (http://david.abcc.ncifcrf.gov/) to find functional aspects that are specific for the different gene sets, which would serve as an indicator for the function of the PCDH cluster genes. Particular analyses involved the beta-linked gene set against all genes in three sets as the background, beta-linked gene set against all human genes, and alpha- and beta-linked gene sets against all human genes.
The Genome of the Netherlands as control
We used the general Dutch population data as the control. The Genome of the Netherlands database is represented by the general Dutch population (www.nlgenome.nl) (Consortium 2014). That project covered an area of 41,543Â km2 with more than 17Â million inhabitants and included the province Limburg from where the present cohort was assembled. The number of subjects is 498 in the extractable genome database. A set of 353 non-pseudo-genes of olfactory receptors (OR) was used to extract genome information. When no information of the variant was available in the Genome of the Netherlands, we used the data from large populations in the dbSNP database (www.ncbi.nlm.nih.gov/snp).
Statistics
From the 30 extreme obese subjects, 12 are PCDH-variant carriers, whereas 18 are non-carriers (Table 1). Chi-square analysis was used to test whether the total number of damaging missense variant hits in OR genes or in subsets of those genes was different between the two groups. A hit is defined as a hetero- or homozygous carrier of the altered allele. A similar analysis was done to investigate the hits of damaging missense variants in OR genes between carriers of the previously identified PCDH alpha-, beta- or gamma-cluster variants and the rest of the subjects. Mann–Whitney test was used as an alternative method to compare PCDH-variant carriers with non-carriers. A P = 0.05 was regarded as the threshold for significant deviation from the null hypothesis.
Genetic association was assessed by Chi-square analysis comparing the observed genotype frequencies in this extreme obese cohort with the genotype frequencies reported in the Genome of The Netherlands database or reported in the dbSNP database.
The comparison on the burden of damaging missense alleles was analyzed by Chi-square with Yates’ correction. In total, 2266 missense polymorphisms and 319 damaging missense polymorphisms were detected either in the present cohort or in the controls. Those were used to extract the number of total missense and damaging missense alleles as the input for the data shown in Table 3.
Results
Relation between the PCDH cluster and the olfactory transduction pathway
In the olfactory transduction pathway, we scored in total 883 hits (hetero- or homozygous carrier of the altered allele) of predicted damaging missense variants in 110 polymorphisms. All genes appeared to belong to the superfamily of OR genes (Malnic et al. 2004). Almost half of these hits, 432, were in genes on chromosome 11 and 146 in genes on chromosome 1, whereas 307 were in genes distributed over 15 other chromosomes. For the total set of hits of predicted damaging missense variants in the olfactory transduction pathway, and in the three subgroups according to the chromosome location, we tested by χ 2 analysis whether there might be a higher number of OR-variant hits in PCDH-variant carriers than in non-carriers. The results showed that significance was only reached for the chromosome 1 subgroup (P = 0.023). Mann–Whitney test gave similar results with significance only for the chromosome 1 subgroup (P = 0.031, Fig. 2b).
To check the specificity of such co-damaging missense variation between two pathways, we tested the tight junction and ABC receptor pathways in a similar way. These two pathways were selected due to their high ratio of damaging variants as mentioned later. For both pathways, no significant relation between their variants and PCDH variation was detected by χ 2 analysis (P = 0.78 and P = 0.16, respectively) nor by Mann–Whitney test (P = 0.79 and P = 0.06, respectively). Because of the relatively low number of variants in these two pathways, no subdivision according to chromosomes or genetic loci was made.
Detailed analysis of olfactory transduction genes on chromosome 1q
The 146 hits of predicted damaging missense variants in OR genes on chromosome 1 were all located on the long arm and divided over two loci at 1q23 (n = 26) and 1q44 (n = 120). In order to assess whether the observed relation with PCDH variation could be ascribed to one of those loci or to both, χ 2 analysis and the Mann–Whitney test were used to compare the hits of variation between PCDH-variant carriers and non-carriers. χ 2 analysis resulted in a trend for difference for 1q23 (P = 0.062) and 1q44 (P = 0.092), whereas the Mann–Whitney test resulted in significance for 1q23 (P = 0.037) and a trend for 1q44 (P = 0.060). Apparently, both loci contribute to the effect suggesting that genes from both locations cooperate with PCDH genes in the genetic predisposition to extreme obesity.
Table 2 gives a detailed overview of the involved OR genes together with the polymorphisms and the IDs of the subjects who are carriers of the risk alleles. It shows that relatively common polymorphisms as well as rare variants are involved. For all except one variant, the observed allele frequencies were in line with that of the control general populations. A comparison of the observed genotype distribution of the G/T variant in the OR14C36 gene (rs28545014; D231Y) with that reported for the Genome of the Netherlands indicated a genetic association between this missense variant and extreme obesity (P < 0.0001). Using frequencies of the CS-Agilent population or of the ESP-cohort populations from dbSNP did not alter this finding (P = 0.001).
Detailed analysis of 1q OR genes and the PCDH clusters
To determine whether the variation in the 1q OR genes was linked with one of the PCDH clusters, the occurrence of hits of 1q OR predicted damaging missense variants was compared per cluster among the subjects with and without a PCDH variant. Similar analyses were performed also for the 1q23 and 1q44 loci separately. Subjects (ID # 14 and #19) with variants in genes of two different PCDH clusters were counted in both of the clusters. Significance for χ 2 analysis was only observed for subjects with a PCDH beta-cluster variant regarding all 1q OR variants (P < 0.0001, Fig. 3b), variants on 1q23 (P = 0.014, Fig. 3c) and on 1q44 (P = 0.004, Fig. 3d). Mann–Whitney test gave similar P values. It indicates that the genetic cooperation between damaging variation in the OR genes on chromosome 1q and variation in the PCDH genes is confined to the PCDH beta-cluster genes and does not occur with the alpha- and gamma-clusters.
Extended functional assessment of PCDH gene clusters
The functions of PCDH gene clusters were indirectly assessed by the functions of their linked genes. Sets of coexpressing genes were generated for each of the PCDH clusters (Table S1) and functionally annotated to obtain an indication of the specific function of the different PCDH gene clusters. The outcome (Table S2) indicates that the PCDH beta-cluster genes are linked to the control of tissue development (‘DNA binding,’ ‘transcription regulation,’ ‘cell cycle,’ ‘cell cycle phase’), whereas the PCDH-alpha and gamma genes seem to be more engaged in neuron organization and function (‘neuron projection,’ ‘cytoplasm,’ ‘synapse,’ ‘axon’).
The relevance of high number of damaging missense variants in the olfactory transduction pathway in extreme obesity
Olfactory transduction pathways showed high number of damaging missense variants in our extreme obesity cohort. We evaluated the ratio of predicted damaging missense variant hits to total missense variant hits for various pathways. The olfactory transduction pathway together with the tight junction pathway and ABC transporters pathway formed the top three (Table S3). However, comparing our cohort with the control on the number of damaging missense alleles in the olfactory transduction pathway per person showed that extreme obese subjects even had a lower burden than the general population (Table 3).
Discussion
The main olfactory system consists of the olfactory bulb, where at the glomeruli neuronal cells are connected with OR cells. When volatile chemicals are recognized by ORs, a signal is transduced into perception and generates the sensation of smell (Cleland 2014). Obviously, smell is an important determinant of food selection and consumption (Jaeger et al. 2013; Blundell and Finlayson 2004; Shepherd 2006; Reed and Knaapila 2010; Sorensen et al. 2003). Close to 400 human OR genes distributed over 21 human chromosomes constitute one of the largest gene families (Hasin-Brumshtein et al. 2009; Malnic et al. 2004). The high redundancy of ORs may make the system more ‘tolerant’ to damaging mutations. The embedded high genetic variability of OR genes is well recognized (Mainland et al. 2014).
In the present study, we also confirmed that a relatively large amount of missense variants with a predicted damaging impact on protein function occurred in extreme obesity. However, compared to the general population, it even seems that extreme obese subjects carried less damaging alleles. Nevertheless, different number of subjects per sequenced cohort and different methods used in exome capture and sequencing may eventually lead to the difference. In the general Dutch population, only 233 out of total 319 damaging missense polymorphisms were detected. Because low-frequent damaging variants in the general population are easy to be missed in our small cohort, this might lead to an under-estimated burden in our extreme obesity cohort in comparison with that of the control. A rough calculation taking this under-estimation into account would lead to no significant difference on the burden. For the common 86 damaging missense polymorphisms, there was also no difference on the burden between the extreme obese and general Dutch population (29.1 vs. 29.5). These results suggest at least that the high damaging genetic variation of the olfactory transduction pathway, in either ratio or absolute burden, is not specific to extreme obesity. However, as indicated in the present and in few other studies, there can still be genetic association between particular ORs and obesity (Choquette et al. 2012; Jaeger et al. 2013; Lunde et al. 2012; Jarick et al. 2011).
Here, we have shown that common and rare damaging missense variants in OR genes on chromosome 1q are carried more often by extreme obese subjects who also carry variants with a moderate-to-low predicted biological effect in the genes of the PCDH cluster on chromosome 5q31, in particular in the genes of the beta-cluster. Our observations suggest a synergistic effect of genetic variation in genes of the PCDH beta-cluster and of the olfactory transduction pathway on chromosome 1. The fact that the rare variants in the PCDH genes were previously shown to be implicated in extreme obesity [8] suggests that variation in the OR genes and PCDH genes exerts an additive effect on the risk of extreme obesity.
The link between the 1q OR and PCDH beta-cluster genes suggests that they are involved in the same biological process. The PCDH genes play a role in establishing neuronal connectivity and synaptic specificity and maturation in various regions of the central nervous system (Li et al. 2010). Although no clear specific functions have been annotated on the PCDH clusters, some experimental evidence has shown non-exclusive roles of PCDH clusters in the development of the olfactory system. During development of the mouse, the maturation of the olfactory system requires the coalescence of homotypic olfactory sensory neurons into glomeruli of the olfactory bulb, which are specified by ORs. Hasegawa et al. (2012) have proposed that this is accomplished by two counteracting processes, of which the coalescence is led by axon guidance molecules like ephrin, semaphorin, plexin, BIG-2 and kirrel2/3, and the repulsion of homotypic axons is mediated by the alpha-protocadherins. However, the development of the olfactory system is more complex and is also accompanied by the expression of the PCDH beta- and gamma-cluster genes (Hirano et al. 2012; Ledderose et al. 2013). The beta-cluster has been proposed to play an important role in neuron cell individuality [10]. It is in line with our bioinformatics analysis that suggests a specific role of those genes in the control of gene expression and cell division, and as such in tissue development. Our results also suggest that the PCDH beta-cluster genes are involved in the development of glomeruli of the olfactory system that in humans are specified by the OR genes on chromosome 1q.
While studies on gene function in the mouse indicate that the PCDH and OR genes mediate the development and fine-tuning of the olfactory sensory system, our results suggest that in humans OR genes on chromosome 1q and PCDH beta-cluster genes may form two parts of one axis. The co-impairment of these two parts can contribute to the risk of extreme obesity at a relatively young age. Indeed, the human olfactory system and its connection with the brain are complex, and its importance for food selection and food intake is obvious but has been largely overlooked (Shepherd 2006; Palouzier-Paulignan et al. 2012). As a consequence, it is unknown how the co-impairment of OR genes and PCDH genes influences human olfaction. Theoretically, a poor performance of this axis may inhibit the maturation of certain olfactory glomeruli and induce less strong synaptic connections with certain sensory neurons leading to impairment of particular sense signals (Reed and Knaapila 2010), which may reduce the food-rewarding circuits and stimulate the body to eat more for compensation. On the other hand, the inhibition of certain types of olfactory glomeruli may lead to oversensitivity of others triggering stronger hedonic signals to the brain (Blundell and Finlayson 2004). In this sense, altered olfaction may also stimulate the food attraction. Unfortunately, in humans, only a few ORs have been annotated for their ligands. The gene OR10J5 located at 1q23 is the human homolog of the murine gene for OR23, which recognizes the odorant lyral with a perception of lemony odor (Malnic et al. 2004). The OR2M7 gene on 1q44 is likely to be involved in the ability to smell methanethiol, which is set free in the urine after eating asparagus (Eriksson et al. 2010). Obviously, knowledge about the ligands of the specific ORs on chromosome 1q seems crucial for solving the puzzle of these genes and their role in food intake.
The number of genetic studies looking into inter-individual variation in smell perception is now rapidly increasing (Keller et al. 2007; Reed and Knaapila 2010). However, studies making the link with food intake are still very limited. Lunde et al. (2012) showed that genetic variation in the OR7D4 gene determines the level of aversion to consume pork from male pigs containing androstenone. Moreover, genetic variation in this gene and other OR7 genes was reported to display genetic association with eating behavior and adiposity (Choquette et al. 2012). In another study, Jaeger et al. (2013) were able to link preference for foods and products with added beta-ionone to genetic variation in the OR5A1 gene. Interestingly, a common copy number variant region was identified on chromosome 11q11 associated with obesity and exclusively covering three OR genes, OR4p4, OR4S2 and OR4C6 (Jarick et al. 2011). Our present results may add to those studies by suggesting genetic association between variation in the OR14C36 gene on 1q44 and extreme obesity. Altogether, it is surprising how little genetic evidence has been reported for the involvement of the OR genes in weight regulation. This may be due to the fact that complex gene–gene interactions play an important role instead of single variations in the OR genes. In this respect, a polygenic model will be more suitable to analyze OR gene effects for obesity just like we have shown here.
It has been shown that morbidly obese subjects have an olfactory dysfunction in comparison with moderately obese individuals (Richardson et al. 2004). This has led to the debate whether weight gain leads to impairment of the olfactory function or whether it is the other way round. The fact that weight loss after a gastric bypass does not improve the olfactory sensation suggests that it is the olfactory dysfunction that contributes to extreme obesity (Richardson et al. 2012). Our results may support the concept that weight gain can be a consequence of an impaired olfactory system, since a genetic influence from the olfactory system together with PCDH genes underlying the risk of extreme obesity is suggested here.
Recently, it was shown that ORs are expressed in other tissues than the olfactory bulb. After feeding obesity-prone and obesity-resistant rats a high-fat diet for 14Â days, microarray analysis showed up-regulation of four OR genes in duodenal enterocytes of the obesity-prone animals (Primeaux et al. 2013). It suggests that these receptors play a role in sensing nutrients or their metabolites as they pass through the lumen of the gut, and assist in regulating their uptake. In lung epithelium of the mouse, it was found that beta-carotene simultaneously up-regulates 26 OR genes and 9 genes of the PCDH beta-cluster. It was proposed that those genes might be active in the highly innervated pulmonary neuroendocrine cell clusters (van Helden et al. 2010), which may have various functions like lung tissue repair, sensing volatile chemicals and monitoring oxygen levels in the inhaled air (Song et al. 2012; Gu et al. 2014). A link with chronic obstructive pulmonary disease (COPD) has been suggested (Gu et al. 2014), but not with obesity. Obviously, our knowledge about the functional cooperation of PCDH and OR genes in the human body is still very limited, and this cooperation may not only occur in the extreme obesity context. Until we obtain more knowledge, our present findings cannot be definitely ascribed to aspects of smell sensation.
The present findings are based on the analysis of whole exome sequencing data and partly comply with the suggested guidelines for quality (MacArthur et al. 2014). For instance, we predicted missense variants to be damaging by using a combination of two functional annotation programs, SIFT and PolyPhen2. Also, we did not assume that the implicated variants or combinations thereof were fully penetrant. Nevertheless, our work is based on a limited number of subjects and our findings need to be confirmed in a similar, independent study.
Conclusions
In summary, we found that predicted damaging missense variants in OR genes on chromosome 1q and rare damaging variants in PCDHÂ beta-cluster genes on 5q31 co-localize in subjects with extreme obesity, suggesting a synergistic effect in the predisposition to this condition. Further, mechanistic insight may come from the identification of the ligands of the 1q ORs linked to the PCDH beta-cluster. Eventually, this may provide the possibility to manipulate food flavor in a way to reduce the risk of overeating and extreme obesity in some genetically predisposed subjects.
References
Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, Kondrashov AS, Sunyaev SR (2010) A method and server for predicting damaging missense mutations. Nat Methods 7(4):248–249. doi:10.1038/nmeth0410-248
Baughman JM, Nilsson R, Gohil VM, Arlow DH, Gauhar Z, Mootha VK (2009) A computational screen for regulators of oxidative phosphorylation implicates SLIRP in mitochondrial RNA homeostasis. PLoS Genet 5(8):e1000590. doi:10.1371/journal.pgen.1000590
Berndt SI, Gustafsson S, Magi R, Ganna A, Wheeler E, Feitosa MF, Justice AE, Monda KL, Croteau-Chonka DC, Day FR, Esko T, Fall T, Ferreira T, Gentilini D, Jackson AU, Luan J, Randall JC, Vedantam S, Willer CJ, Winkler TW, Wood AR, Workalemahu T, Hu YJ, Lee SH, Liang L, Lin DY, Min JL, Neale BM, Thorleifsson G, Yang J, Albrecht E, Amin N, Bragg-Gresham JL, Cadby G, den Heijer M, Eklund N, Fischer K, Goel A, Hottenga JJ, Huffman JE, Jarick I, Johansson A, Johnson T, Kanoni S, Kleber ME, Konig IR, Kristiansson K, Kutalik Z, Lamina C, Lecoeur C, Li G, Mangino M, McArdle WL, Medina-Gomez C, Muller-Nurasyid M, Ngwa JS, Nolte IM, Paternoster L, Pechlivanis S, Perola M, Peters MJ, Preuss M, Rose LM, Shi J, Shungin D, Smith AV, Strawbridge RJ, Surakka I, Teumer A, Trip MD, Tyrer J, Van Vliet-Ostaptchouk JV, Vandenput L, Waite LL, Zhao JH, Absher D, Asselbergs FW, Atalay M, Attwood AP, Balmforth AJ, Basart H, Beilby J, Bonnycastle LL, Brambilla P, Bruinenberg M, Campbell H, Chasman DI, Chines PS, Collins FS, Connell JM, Cookson WO, de Faire U, de Vegt F, Dei M, Dimitriou M, Edkins S, Estrada K, Evans DM, Farrall M, Ferrario MM, Ferrieres J, Franke L, Frau F, Gejman PV, Grallert H, Gronberg H, Gudnason V, Hall AS, Hall P, Hartikainen AL, Hayward C, Heard-Costa NL, Heath AC, Hebebrand J, Homuth G, Hu FB, Hunt SE, Hypponen E, Iribarren C, Jacobs KB, Jansson JO, Jula A, Kahonen M, Kathiresan S, Kee F, Khaw KT, Kivimaki M, Koenig W, Kraja AT, Kumari M, Kuulasmaa K, Kuusisto J, Laitinen JH, Lakka TA, Langenberg C, Launer LJ, Lind L, Lindstrom J, Liu J, Liuzzi A, Lokki ML, Lorentzon M, Madden PA, Magnusson PK, Manunta P, Marek D, Marz W, Mateo Leach I, McKnight B, Medland SE, Mihailov E, Milani L, Montgomery GW, Mooser V, Muhleisen TW, Munroe PB, Musk AW, Narisu N, Navis G, Nicholson G, Nohr EA, Ong KK, Oostra BA, Palmer CN, Palotie A, Peden JF, Pedersen N, Peters A, Polasek O, Pouta A, Pramstaller PP, Prokopenko I, Putter C, Radhakrishnan A, Raitakari O, Rendon A, Rivadeneira F, Rudan I, Saaristo TE, Sambrook JG, Sanders AR, Sanna S, Saramies J, Schipf S, Schreiber S, Schunkert H, Shin SY, Signorini S, Sinisalo J, Skrobek B, Soranzo N, Stancakova A, Stark K, Stephens JC, Stirrups K, Stolk RP, Stumvoll M, Swift AJ, Theodoraki EV, Thorand B, Tregouet DA, Tremoli E, Van der Klauw MM, van Meurs JB, Vermeulen SH, Viikari J, Virtamo J, Vitart V, Waeber G, Wang Z, Widen E, Wild SH, Willemsen G, Winkelmann BR, Witteman JC, Wolffenbuttel BH, Wong A, Wright AF, Zillikens MC, Amouyel P, Boehm BO, Boerwinkle E, Boomsma DI, Caulfield MJ, Chanock SJ, Cupples LA, Cusi D, Dedoussis GV, Erdmann J, Eriksson JG, Franks PW, Froguel P, Gieger C, Gyllensten U, Hamsten A, Harris TB, Hengstenberg C, Hicks AA, Hingorani A, Hinney A, Hofman A, Hovingh KG, Hveem K, Illig T, Jarvelin MR, Jockel KH, Keinanen-Kiukaanniemi SM, Kiemeney LA, Kuh D, Laakso M, Lehtimaki T, Levinson DF, Martin NG, Metspalu A, Morris AD, Nieminen MS, Njolstad I, Ohlsson C, Oldehinkel AJ, Ouwehand WH, Palmer LJ, Penninx B, Power C, Province MA, Psaty BM, Qi L, Rauramaa R, Ridker PM, Ripatti S, Salomaa V, Samani NJ, Snieder H, Sorensen TI, Spector TD, Stefansson K, Tonjes A, Tuomilehto J, Uitterlinden AG, Uusitupa M, van der Harst P, Vollenweider P, Wallaschofski H, Wareham NJ, Watkins H, Wichmann HE, Wilson JF, Abecasis GR, Assimes TL, Barroso I, Boehnke M, Borecki IB, Deloukas P, Fox CS, Frayling T, Groop LC, Haritunian T, Heid IM, Hunter D, Kaplan RC, Karpe F, Moffatt MF, Mohlke KL, O’Connell JR, Pawitan Y, Schadt EE, Schlessinger D, Steinthorsdottir V, Strachan DP, Thorsteinsdottir U, van Duijn CM, Visscher PM, Di Blasio AM, Hirschhorn JN, Lindgren CM, Morris AP, Meyre D, Scherag A, McCarthy MI, Speliotes EK, North KE, Loos RJ, Ingelsson E (2013) Genome-wide meta-analysis identifies 11 new loci for anthropometric traits and provides insights into genetic architecture. Nat Genet 45(5):501–512. doi:10.1038/ng.2606
Blundell JE, Finlayson G (2004) Is susceptibility to weight gain characterized by homeostatic or hedonic risk factors for overconsumption? Physiol Behav 82(1):21–25. doi:10.1016/j.physbeh.2004.04.021
Choquette AC, Bouchard L, Drapeau V, Lemieux S, Tremblay A, Bouchard C, Vohl MC, Perusse L (2012) Association between olfactory receptor genes, eating behavior traits and adiposity: results from the Quebec Family Study. Physiol Behav 105(3):772–776. doi:10.1016/j.physbeh.2011.10.015
Cleland TA (2014) Construction of odor representations by olfactory bulb microcircuits. Prog Brain Res 208:177–203. doi:10.1016/b978-0-444-63350-7.00007-3
Consortium GotN (2014) Whole-genome sequence variation, population structure and demographic history of the Dutch population. Nat Genet 46(8):818–825. doi:10.1038/ng.3021
Eriksson N, Macpherson JM, Tung JY, Hon LS, Naughton B, Saxonov S, Avey L, Wojcicki A, Pe’er I, Mountain J (2010) Web-based, participant-driven studies yield novel genetic associations for common traits. PLoS Genet 6(6):e1000993. doi:10.1371/journal.pgen.1000993
Farooqi IS, O’Rahilly S (2007) Genetic factors in human obesity. Obes Rev 8(Suppl 1):37–40. doi:10.1111/j.1467-789X.2007.00315.x
Gu X, Karp PH, Brody SL, Pierce RA, Welsh MJ, Holtzman MJ, Ben-Shahar Y (2014) Chemosensory functions for pulmonary neuroendocrine cells. Am J Respir Cell Mol Biol 50(3):637–646. doi:10.1165/rcmb.2013-0199OC
Hasegawa S, Hirabayashi T, Kondo T, Inoue K, Esumi S, Okayama A, Hamada S, Yagi T (2012) Constitutively expressed Protocadherin-alpha regulates the coalescence and elimination of homotypic olfactory axons through its cytoplasmic region. Front Mol Neurosci 5:97. doi:10.3389/fnmol.2012.00097
Hasin-Brumshtein Y, Lancet D, Olender T (2009) Human olfaction: from genomic variation to phenotypic diversity. Trends Genet TIG 25(4):178–184. doi:10.1016/j.tig.2009.02.002
Hirano K, Kaneko R, Izawa T, Kawaguchi M, Kitsukawa T, Yagi T (2012) Single-neuron diversity generated by Protocadherin-beta cluster in mouse central and peripheral nervous systems. Front Mol Neurosci 5:90. doi:10.3389/fnmol.2012.00090
Jaeger SR, McRae JF, Bava CM, Beresford MK, Hunter D, Jia Y, Chheang SL, Jin D, Peng M, Gamble JC, Atkinson KR, Axten LG, Paisley AG, Tooman L, Pineau B, Rouse SA, Newcomb RD (2013) A Mendelian trait for olfactory sensitivity affects odor experience and food selection. Curr Biol CB 23(16):1601–1605. doi:10.1016/j.cub.2013.07.030
Jarick I, Vogel CI, Scherag S, Schafer H, Hebebrand J, Hinney A, Scherag A (2011) Novel common copy number variation for early onset extreme obesity on chromosome 11q11 identified by a genome-wide analysis. Hum Mol Genet 20(4):840–852. doi:10.1093/hmg/ddq518
Keller A, Zhuang H, Chi Q, Vosshall LB, Matsunami H (2007) Genetic variation in a human odorant receptor alters odour perception. Nature 449(7161):468–472. doi:10.1038/nature06162
Ledderose J, Dieter S, Schwarz MK (2013) Maturation of postnatally generated olfactory bulb granule cells depends on functional gamma-protocadherin expression. Sci Rep 3:1514. doi:10.1038/srep01514
Li Y, Serwanski DR, Miralles CP, Fiondella CG, Loturco JJ, Rubio ME, De Blas AL (2010) Synaptic and nonsynaptic localization of protocadherin-gammaC5 in the rat brain. J Comp Neurol 518(17):3439–3463. doi:10.1002/cne.22390
Loos RJF (2012) Genetic determinants of common obesity and their value in prediction. Best Pract Res Clin Endocrinol Metab 26(2):211–226. doi:10.1016/j.beem.2011.11.003
Lunde K, Egelandsdal B, Skuterud E, Mainland JD, Lea T, Hersleth M, Matsunami H (2012) Genetic variation of an odorant receptor OR7D4 and sensory perception of cooked meat containing androstenone. PLoS ONE 7(5):e35259. doi:10.1371/journal.pone.0035259
MacArthur DG, Manolio TA, Dimmock DP, Rehm HL, Shendure J, Abecasis GR, Adams DR, Altman RB, Antonarakis SE, Ashley EA, Barrett JC, Biesecker LG, Conrad DF, Cooper GM, Cox NJ, Daly MJ, Gerstein MB, Goldstein DB, Hirschhorn JN, Leal SM, Pennacchio LA, Stamatoyannopoulos JA, Sunyaev SR, Valle D, Voight BF, Winckler W, Gunter C (2014) Guidelines for investigating causality of sequence variants in human disease. Nature 508(7497):469–476. doi:10.1038/nature13127
Mainland JD, Keller A, Li YR, Zhou T, Trimmer C, Snyder LL, Moberly AH, Adipietro KA, Liu WL, Zhuang H, Zhan S, Lee SS, Lin A, Matsunami H (2014) The missense of smell: functional variability in the human odorant receptor repertoire. Nat Neurosci 17(1):114–120. doi:10.1038/nn.3598
Malnic B, Godfrey PA, Buck LB (2004) The human olfactory receptor gene family. Proc Natl Acad Sci USA 101(8):2584–2589
Mariman EC, Bouwman FG, Aller EE, van Baak MA, Wang P (2014) High frequency of rare variants with a moderate-to-high predicted biological effect in protocadherin genes of extremely obese. Genes Nutr 9(3):399. doi:10.1007/s12263-014-0399-1
Ng PC, Henikoff S (2003) SIFT: predicting amino acid changes that affect protein function. Nucleic Acids Res 31(13):3812–3814
Palouzier-Paulignan B, Lacroix MC, Aime P, Baly C, Caillol M, Congar P, Julliard AK, Tucker K, Fadool DA (2012) Olfaction under metabolic influences. Chem Senses 37(9):769–797. doi:10.1093/chemse/bjs059
Primeaux SD, Braymer HD, Bray GA (2013) High fat diet differentially regulates the expression of olfactory receptors in the duodenum of obesity-prone and obesity-resistant rats. Dig Dis Sci 58(1):72–76. doi:10.1007/s10620-012-2421-z
Reed DR, Knaapila A (2010) Genetics of taste and smell: poisons and pleasures. Prog Mol Biol Transl Sci 94:213–240. doi:10.1016/b978-0-12-375003-7.00008-x
Reumers J, Schymkowitz J, Ferkinghoff-Borg J, Stricher F, Serrano L, Rousseau F (2005) SNPeffect: a database mapping molecular phenotypic effects of human non-synonymous coding SNPs. Nucleic acids research 33(Database issue):D527–D532. doi:10.1093/nar/gki086
Richardson BE, Vander Woude EA, Sudan R, Thompson JS, Leopold DA (2004) Altered olfactory acuity in the morbidly obese. Obes Surg 14(7):967–969. doi:10.1381/0960892041719617
Richardson BE, Vanderwoude EA, Sudan R, Leopold DA, Thompson JS (2012) Gastric bypass does not influence olfactory function in obese patients. Obes Surg 22(2):283–286. doi:10.1007/s11695-011-0487-x
Shepherd GM (2006) Smell images and the flavour system in the human brain. Nature 444(7117):316–321. doi:10.1038/nature05405
Song H, Yao E, Lin C, Gacayan R, Chen MH, Chuang PT (2012) Functional characterization of pulmonary neuroendocrine cells in lung development, injury, and tumorigenesis. Proc Natl Acad Sci USA 109(43):17531–17536. doi:10.1073/pnas.1207238109
Sorensen LB, Moller P, Flint A, Martens M, Raben A (2003) Effect of sensory perception of foods on appetite and food intake: a review of studies on humans. Int J Obes Relat Metab Disord 27(10):1152–1166. doi:10.1038/sj.ijo.0802391
Speliotes EK, Willer CJ, Berndt SI, Monda KL, Thorleifsson G, Jackson AU, Lango Allen H, Lindgren CM, Luan J, Magi R, Randall JC, Vedantam S, Winkler TW, Qi L, Workalemahu T, Heid IM, Steinthorsdottir V, Stringham HM, Weedon MN, Wheeler E, Wood AR, Ferreira T, Weyant RJ, Segre AV, Estrada K, Liang L, Nemesh J, Park JH, Gustafsson S, Kilpelainen TO, Yang J, Bouatia-Naji N, Esko T, Feitosa MF, Kutalik Z, Mangino M, Raychaudhuri S, Scherag A, Smith AV, Welch R, Zhao JH, Aben KK, Absher DM, Amin N, Dixon AL, Fisher E, Glazer NL, Goddard ME, Heard-Costa NL, Hoesel V, Hottenga JJ, Johansson A, Johnson T, Ketkar S, Lamina C, Li S, Moffatt MF, Myers RH, Narisu N, Perry JR, Peters MJ, Preuss M, Ripatti S, Rivadeneira F, Sandholt C, Scott LJ, Timpson NJ, Tyrer JP, van Wingerden S, Watanabe RM, White CC, Wiklund F, Barlassina C, Chasman DI, Cooper MN, Jansson JO, Lawrence RW, Pellikka N, Prokopenko I, Shi J, Thiering E, Alavere H, Alibrandi MT, Almgren P, Arnold AM, Aspelund T, Atwood LD, Balkau B, Balmforth AJ, Bennett AJ, Ben-Shlomo Y, Bergman RN, Bergmann S, Biebermann H, Blakemore AI, Boes T, Bonnycastle LL, Bornstein SR, Brown MJ, Buchanan TA, Busonero F, Campbell H, Cappuccio FP, Cavalcanti-Proenca C, Chen YD, Chen CM, Chines PS, Clarke R, Coin L, Connell J, Day IN, den Heijer M, Duan J, Ebrahim S, Elliott P, Elosua R, Eiriksdottir G, Erdos MR, Eriksson JG, Facheris MF, Felix SB, Fischer-Posovszky P, Folsom AR, Friedrich N, Freimer NB, Fu M, Gaget S, Gejman PV, Geus EJ, Gieger C, Gjesing AP, Goel A, Goyette P, Grallert H, Grassler J, Greenawalt DM, Groves CJ, Gudnason V, Guiducci C, Hartikainen AL, Hassanali N, Hall AS, Havulinna AS, Hayward C, Heath AC, Hengstenberg C, Hicks AA, Hinney A, Hofman A, Homuth G, Hui J, Igl W, Iribarren C, Isomaa B, Jacobs KB, Jarick I, Jewell E, John U, Jorgensen T, Jousilahti P, Jula A, Kaakinen M, Kajantie E, Kaplan LM, Kathiresan S, Kettunen J, Kinnunen L, Knowles JW, Kolcic I, Konig IR, Koskinen S, Kovacs P, Kuusisto J, Kraft P, Kvaloy K, Laitinen J, Lantieri O, Lanzani C, Launer LJ, Lecoeur C, Lehtimaki T, Lettre G, Liu J, Lokki ML, Lorentzon M, Luben RN, Ludwig B, Magic Manunta P, Marek D, Marre M, Martin NG, McArdle WL, McCarthy A, McKnight B, Meitinger T, Melander O, Meyre D, Midthjell K, Montgomery GW, Morken MA, Morris AP, Mulic R, Ngwa JS, Nelis M, Neville MJ, Nyholt DR, O’Donnell CJ, O’Rahilly S, Ong KK, Oostra B, Pare G, Parker AN, Perola M, Pichler I, Pietilainen KH, Platou CG, Polasek O, Pouta A, Rafelt S, Raitakari O, Rayner NW, Ridderstrale M, Rief W, Ruokonen A, Robertson NR, Rzehak P, Salomaa V, Sanders AR, Sandhu MS, Sanna S, Saramies J, Savolainen MJ, Scherag S, Schipf S, Schreiber S, Schunkert H, Silander K, Sinisalo J, Siscovick DS, Smit JH, Soranzo N, Sovio U, Stephens J, Surakka I, Swift AJ, Tammesoo ML, Tardif JC, Teder-Laving M, Teslovich TM, Thompson JR, Thomson B, Tonjes A, Tuomi T, van Meurs JB, van Ommen GJ, Vatin V, Viikari J, Visvikis-Siest S, Vitart V, Vogel CI, Voight BF, Waite LL, Wallaschofski H, Walters GB, Widen E, Wiegand S, Wild SH, Willemsen G, Witte DR, Witteman JC, Xu J, Zhang Q, Zgaga L, Ziegler A, Zitting P, Beilby JP, Farooqi IS, Hebebrand J, Huikuri HV, James AL, Kahonen M, Levinson DF, Macciardi F, Nieminen MS, Ohlsson C, Palmer LJ, Ridker PM, Stumvoll M, Beckmann JS, Boeing H, Boerwinkle E, Boomsma DI, Caulfield MJ, Chanock SJ, Collins FS, Cupples LA, Smith GD, Erdmann J, Froguel P, Gronberg H, Gyllensten U, Hall P, Hansen T, Harris TB, Hattersley AT, Hayes RB, Heinrich J, Hu FB, Hveem K, Illig T, Jarvelin MR, Kaprio J, Karpe F, Khaw KT, Kiemeney LA, Krude H, Laakso M, Lawlor DA, Metspalu A, Munroe PB, Ouwehand WH, Pedersen O, Penninx BW, Peters A, Pramstaller PP, Quertermous T, Reinehr T, Rissanen A, Rudan I, Samani NJ, Schwarz PE, Shuldiner AR, Spector TD, Tuomilehto J, Uda M, Uitterlinden A, Valle TT, Wabitsch M, Waeber G, Wareham NJ, Watkins H, Procardis C, Wilson JF, Wright AF, Zillikens MC, Chatterjee N, McCarroll SA, Purcell S, Schadt EE, Visscher PM, Assimes TL, Borecki IB, Deloukas P, Fox CS, Groop LC, Haritunians T, Hunter DJ, Kaplan RC, Mohlke KL, O’Connell JR, Peltonen L, Schlessinger D, Strachan DP, van Duijn CM, Wichmann HE, Frayling TM, Thorsteinsdottir U, Abecasis GR, Barroso I, Boehnke M, Stefansson K, North KE, McCarthy MI, Hirschhorn JN, Ingelsson E, Loos RJ (2010) Association analyses of 249,796 individuals reveal 18 new loci associated with body mass index. Nat Genet 42(11):937–948. doi:10.1038/ng.686
van Helden YG, Godschalk RW, Heil SG, Bunschoten A, Hessel S, Amengual J, Bonet ML, von Lintig J, van Schooten FJ, Keijer J (2010) Downregulation of Fzd6 and Cthrc1 and upregulation of olfactory receptors and protocadherins by dietary beta-carotene in lungs of Bcmo1-/-mice. Carcinogenesis 31(8):1329–1337. doi:10.1093/carcin/bgq083
Walley AJ, Blakemore AI, Froguel P (2006) Genetics of obesity and the prediction of risk for health. Human molecular genetics 15(Spec No 2):R124–R130. doi:10.1093/hmg/ddl215
Wheeler E, Huang N, Bochukova EG, Keogh JM, Lindsay S, Garg S, Henning E, Blackburn H, Loos RJ, Wareham NJ, O’Rahilly S, Hurles ME, Barroso I, Farooqi IS (2013) Genome-wide SNP and CNV analysis identifies common and low-frequency variants associated with severe early-onset obesity. Nat Genet 45(5):513–517. doi:10.1038/ng.2607
Conflict of interest
All authors state that there is no conflict of interest. All authors commented on the results and the manuscript.
Author information
Authors and Affiliations
Corresponding author
Electronic supplementary material
Below is the link to the electronic supplementary material.
Table S1
List of genes coexpressed with PCDH gene clusters. Microarray analysis results stored in public databases were analyzed for genes showing coexpression with the different PCDH gene clusters (XLSX 13Â kb)
Table S2
Functional annotation analysis of genes coexpressed with PCDH gene clusters. The genes that coexpress with PCDH gene clusters were analyzed with regard to annotated function using DAVID (XLSX 21Â kb)
Table S3
Total and damaging missense variation in pathways. The pathways that had at least 500 hits for total missense variants and at least a 1Â % ratio of damaging missense variant hits to total missense variant hits are listed (DOCX 16Â kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
Mariman, E.C.M., Szklarczyk, R., Bouwman, F.G. et al. Olfactory receptor genes cooperate with protocadherin genes in human extreme obesity. Genes Nutr 10, 16 (2015). https://doi.org/10.1007/s12263-015-0465-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s12263-015-0465-3