|
|
||||||||
CLINICAL STUDY |
1 Clinical Research Group, Department of Child and Adolescent Psychiatry, 2 Department of Child and Adolescent Psychiatry, Rheinische Kliniken Essen, University of Duisburg-Essen, 45147 Essen, Germany, 3 Genome Analysis, Institute of Molecular Biotechnology, 07745 Jena, Germany, 4 Institute of Pharmacology and Toxicology and 5 Institute of Medical Biometry and Epidemiology, University of Marburg, 35039 Marburg, Germany, 6 Institute of Anatomy II, Friedrich-Schiller-University, 07743 Jena, Germany, 7 Aventis Pharma AG, 65812 Bad Soden a Ts, Germany, 8 Treatment Center Insula, 83489 Berchtesgaden, Germany, 9 Spessartklinik Bad Orb, 63619 Bad Orb, Germany, 10 Center for Neurobiology and Behavior, University of Pennsylvania, Philadelphia, Pennsylvania 19104-6140, USA, 11 Department of Paediatric Endocrinology, Charité Childrens Hospital, Humboldt University, 13353 Berlin, Germany, 12 GSF National Research Center for Environment and Health, Institute of Epidemiology, 85764 Neuherberg, Germany, 13 Danish Epidemiology Science Centre, Institute of Preventive Medicine, Copenhagen University Hospital, DK-1399 Copenhagen, Denmark, 14 Institute of Human Nutrition, Royal Veterinary and Agricultural University, DK-1958 Frederiksberg Copenhagen, Denmark, 15 Steno Diabetes Center, DK-2820 Gentofte, Denmark, 16 Hotel Dieu, Laboratoire de Nutrition, EA, 3502 INSERM Avenir, 75004 Paris, France, 17 PsychoBiology Research Unit, University of Leeds, Leeds LS2 9JT, UK and 18 Department of Pediatrics and Adolescent Medicine, University of Ulm, 89075 Ulm, Germany
(Correspondence should be addressed to J Hebebrand, Rheinische Kliniken Essen, Department of Child and Adolescent Psychiatry, University of Duisburg-Essen, Virchowstr. 174, 45147 Essen, Germany; Email: Johannes.Hebebrand{at}lvr.de)
| Abstract |
|---|
|
|
|---|
Design: This consisted of genomic screening of 13.4 kb encompassing the MCHR1 in extremely obese German children and adolescents and association analyses for two coding single nucleotide polymorphisms (SNPs). To confirm initial positive association results, additional association studies and transmission disequilibrium tests in further German, Danish, French and American samples were conducted. Selected SNPs were investigated using functional in vitro studies and reporter gene assays.
Methods: Single-stranded conformation polymorphism analysis, re-sequencing, PCR-restriction fragment length polymorphism analyses, tetra-primer amplification refractory mutation systems, matrix-assisted laser desorption/ionization time of flight mass spectrometry and reporter gene assays were carried out as well as measuring inositol phosphate formation, inhibition of cAMP formation and activation of p42/44 MAP kinase.
Results: We identified 11 infrequent variations and two SNPs in the MCHR1 coding sequence and 18 SNPs (eight novel) in the flanking sequence. Association and transmission disequilibrium with obesity were detected for several SNPs in independent study groups of German obese children and adolescents and controls. In two German samples, encompassing 4056 and 295 individuals, trends towards association with obesity were detected. Findings in a second epidemiological German sample and in Danish, French and American samples were negative. Functional in vitro studies as well as reporter gene assays revealed no significant results.
Conclusion: Our initial association of MCHR1 alleles/haplotype detected might be related to juvenile-onset obesity, conditional on a particular genetic and/or environmental background. Alternatively, we could not exclude the possibility that the initially detected association represented a false positive finding.
| Introduction |
|---|
|
|
|---|
Both Pmch / and Mchr1 / mice have a reduced fat mass and an increased metabolic rate (1214). Furthermore, Mchr1 / mice were described to either have a similar (13) or a reduced (14) body weight compared with wild-type mice and to be less susceptible (13) or resistant to diet-induced obesity (14). Pmch /; ob/ob mice (double null) had a significantly reduced body weight and fat content compared with ob/ob mice caused by an increased metabolic rate and increased physical activity but not by a reduced food intake (15).
MCHR1 (official gene symbol: GPR24) is expressed peripherally and centrally in hypothalamic nuclei and hippocampal regions involved in olfaction and regulation of feeding behaviour and body weight (79).
Recently, two infrequent missense variations (Y181H and R248Q) were identified in a mutation screen of MCHR1 in 106 individuals with severe, early onset obesity (16); however, in vitro studies showed no evidence for functional implications. Furthermore, no association for single nucleotide polymorphisms (SNPs) rs133072 and rs133073 with obesity-related phenotypes in 541 white adult probands with type 2 diabetes and related disorders was detected (16).
The human reference genome at the University of California Santa Cruz (UCSC) geneme browser (http://genome.ucsc.edu/cgi-bin/hgGateway) (hg17) currently includes seven human MCHR1 mRNA entries (AF490537
[GenBank]
, AB063174
[GenBank]
, AY562945
[GenBank]
, BC001736
[GenBank]
, BC021146
[GenBank]
, BT006725
[GenBank]
and CR456497
[GenBank]
). Each harbours three potential translational start codons (Met1, Met6 and Met70) corresponding to proteins of 422 amino acids (aa) (17), 417 aa and 353 aa respectively (Fig. 1a
). In these entries as well as in the reference protein (Q99705
[GenBank]
), an open reading frame corresponding to a 422 aa protein is stated. The translational start codon of orthologous Mchr1 genes in mice (AF498248
[GenBank]
; Fig. 1b
) (18) and rats (AF008650
[GenBank]
) (19), however, corresponds to human Met70. Therefore, MCHR1 proteins longer than 353 aa seem to be specific to humans.
|
| Materials and methods |
|---|
|
|
|---|
The association and linkage disequilibrium studies were based on six study groups in Marburg (study group 1), Berlin, Munich, Copenhagen, Paris and Philadelphia (study groups 26) respectively:
Study group 1
90th percentile (20; mean BMI 33.4 ± 6.6 kg/m2; mean BMI percentile 99.3; mean age 14.0 ± 2.7 years). For the initial mutation screen we chose a subgroup of 215 extremely obese children and adolescents (127 females: mean BMI 39.8 ± 5.3 kg/m2; mean age 15.3 ± 2.4 years). For transmission disequilibrium tests (TDT; 21) both parents of 525 of the children and adolescents (304 females: mean BMI 32.2 ± 6.0 kg/m2; mean age 13.7 ± 2.8 years) were genotyped for four SNPs (sample described in 22); TDTs for these and 12 other SNPs were also based on a smaller subgroup of 61 families who are comprised of one obese index patient and both biological parents (trios).
15th percentile; mean BMI 18.3 ± 1.1 kg/m2; mean age 25.2 ± 3.7 years).
40th and
60th percentile; mean BMI 21.9 ± 1.1 kg/m2; mean age 24.7 ± 2.6 years).
90th percentile; mean BMI 29.1 ± 3.4 kg/m2; mean age 25.3 ± 3.7 years).
Study group 2
196 (98 females) unrelated German obese children (BMI
90th percentile) and 99 (57 females) unrelated normal weight infants with inborn endocrine disorders were collected in the Department of Paediatric Endocrinology of the Charité Childrens Hospital, Berlin.
Study group 3
50th percentile.
50th percentile.
Study group 4
231 men (BMI
31 kg/m2) with juvenile-onset obesity and 318 male controls (BMI < 31 kg/m2) from the draft board examinations carried out at 1826 years of age (mean age 19.9 ± 1.9 years) in the eastern part of Denmark between 1953 and 1977 (24, 25). Within the draft board population, the criterion for the obese group corresponded approximately to the 99.5th percentile. DNA was sampled at a follow-up examination in 19982000.
Study group 5 468 obese (353 females: mean BMI 47.8 ± 7.8 kg/m2; mean age 43.3 ± 12.2 years) adults and 322 controls (166 females: mean BMI 21.5 ± 1.8 kg/m2; mean age 51.9 ± 10.8 years) were collected in the Laboratoire de Nutrition, Hotel Dieu in Paris, France.
Study group 6 187 trios (153 female index cases) based on 160 extremely obese cases of European American origin and 27 extremely obese cases of African American origin (mean BMI 50.5 ± 9.4 kg/m2; mean age 37.2 ± 9.3 years) and their parents associated in the Center for Neurobiology and Behavior of the University of Pennsylvania, Philadelphia, PA, USA.
Written informed consent was given by all participants and, in the case of minors, by their parents. The ethics committee of the University of Marburg and the local ethics committees of the participating investigators approved genotyping of obese subjects, their parents and controls.
Mutation screen
Single-stranded conformation polymorphism analysis (SSCP) was performed as described previously (26). PCR products of primers corresponding to amplicon 1 of exon 2 (see supplementary Table 1
) were digested by both AluI and MspI (Fermentas, St. Leon-Rot, Germany). For analysis of the second amplicon of exon 2 Crf13I (Fermentas) was used. All amplicons were electrophoresed on 21% acrylamide gels (37.5:1; Q Biogene, Heidelberg, Germany). Gels were run at room temperature for 17 h at 4007 V and additionally at 4 °C for 18 h at 500 V. Gels were silver stained. PCR products of individuals showing aberrant SSCP patterns were sequenced bi-directionally on a LiCor 4200-2 automatic sequencer using the Base ImagIR 4.0 software (MWG Biotech, Ebersberg, Germany) as described previously (26) and by Sequence Laboratories (Seq Lab, Göttingen, Germany) respectively.
|
For initial and confirmatory comparison of genotype and allele frequencies we used the CochranArmitage trend test and the asymptotic Pearson Chi-square test respectively.
TDTs
The initial, confirmatory and post hoc TDTs in the German and American families were carried out using the program Genehunter, version 2.0 beta (http://www.fhcrc.org/labs/kruglyak/Downloads/index.html (27)). Initially, association studies based on the index patients of the German trios had been performed, so that the TDT results are not independent. We performed haplotype TDTs with up to four markers for all 525 trios (study group 1a), where phase can be determined unambiguously.
Sequencing
A genomic region of 13 379 bp encompassing MCHR1 was covered by eight overlapping PCR products named regions AH (see supplementary Table 2
). For each region four to eight nested/semi-nested PCRs were performed. PCR products were sequenced using PCR primers and BigDye Terminator Cycle Sequencing v2.0 kit (Applied Biosystems, Weiterstadt, Germany). Sequencing reactions were electrophoresed on ABI 377 automated sequencers. Base calling was performed using phred (28, 29). Sequence assembly was done using phrap (http://www.phrap.org/phrap.docs/phrap.html). Trace files were inspected visually in gap4 (30).
|
Thirteen SNPs (see supplementary Tables 3 and 4) were genotyped by restriction fragment length polymorphism analyses (RFLP) or tetra-primer amplification refractory mutation system (ARMS; 31). SNP rs133072 was genotyped in study groups 1, 2 and 46 by PCR-RFLP (see supplementary Table 3). SNP rs133072 was genotyped in studygroup 3 by PCR followed by a MassEXTEND reaction and subsequent matrix-assisted laser desorption/ionization time of flight mass spectrometry (MALDITOF MS) analysis according to the manufacturers protocol (hME; Sequenom, San Diego, CA, USA). Four SNPs (SNP001745614, rs133063, rs133064 and SNP001745615) were genotyped by genomic sequencing as described above (see supplementary Table 5).
Haplotype frequencies
These were estimated for four SNPs (rs133068, rs133072, rs133073 and rs133074) in 525 trios (study group 1a) and linkage disequilibria between these SNPs were tested with the program EH, version 1.11 (http://linkage.rockefeller.edu/ott/eh.htm) (32).
Receptor constructs
A nested PCR in a pool containing all tissues of human MTC panel I (BD Clontech, Heidelberg, Germany) was performed (see supplementary Table 6). The PCR product was cloned into pCR2.1-TOPO (Invitrogen, Karlsruhe, Germany). All recombinant plasmids were re-sequenced as described above. Subsequently, the coding sequence of MCHR1 starting with Met1 was cloned into the eukaryotic expression vector pcDNA3.1 (Invitrogen). To allow for immunological detection, MCHR1 clones were tagged with an N-terminal 9 aa residue epitope (YPYDVPDYA) derived from influenza virus haemagglutinin protein (HA-tag) after Met1. We subsequently refer to the numbering of aa according to the 422 aa receptor protein which is in agreement with the annotation in the UCSC genome browser (hg17), the Human Genome Variation database and dbSNP.
Heterologous expression
The MCHR1 variants D32- and N32-MCHR1 were expressed in transiently transfected COS-7 and HEK293 cells. COS-7 cells were grown in Dulbeccos modified Eagles medium (DMEM) supplemented with 10% fetal calf serum (FCS), 100 U/ml penicillin and 100 µg/ml streptomycin, HEK293 cells were maintained in Iscoves Modified Dulbeccos Medium (IMDM) with the same additives. Transient transfection was performed using metafectene (Biontex, Munich, Germany). Cell surface expression of the receptor variants was estimated using an indirect cellular ELISA as previously described (33). HA-tagged receptor variants were detected with a horseradish peroxidase (HRP)-labelled anti-HA monoclonal antibody (12CA5; Roche Diagnostics, Mannheim, Germany). Visualization was carried out using 3,3',5,5' tetramethylbenzidine as a substrate for HRP (1-Step Turbo TMB-ELISA; Pierce Biotechnology, Rockford, IL, USA).
Inositol phosphate (IP) formation
Transfected COS-7 cells were incubated with 2 µCi/ml myo-[2-3H]inositol (15.0 Ci/mmol; Amersham Phar-macia Biotech, Amersham, Bucks, UK) for 18 h. Thereafter, cells were washed once with serum-free DMEM containing 10 mM LiCl and stimulated with MCH (1011 to 105 M) for 1 h at 37 °C. Agonist-induced increases in intracellular IP levels were determined as described (34).
Inhibition of cAMP accumulation
In HEK293 cells, D32- and N32-MCHR1 were transiently coexpressed with the human ß2-adrenoceptor (ß2-AR). Cells were treated with 10 µM isoproterenol to stimulate cAMP formation and 1 µM MCH. As a positive control, ß2-AR was coexpressed with the Gi-coupled M2 muscarinic acetylcholine receptor which was stimulated with 100 µM carbachol. HEK293 cells were labelled with 2 µCi/ml [2,8-3H]-adenine and intracellular [3H]cAMP was determined by anion-exchange chromatography (34, 35).
Western blot analysis
Activation of p42/44 MAP kinase (MAPK) was determined by Western blot analysis as previously described (36). Phosphorylated p42/44 MAPK was detected using a phosphospecific anti-p42/44 MAPK rabbit monoclonal antibody (Cell Signalling, Beverly, MA, USA), reblots as loading control were performed using an anti-p42 MAPK mouse monoclonal antibody (Santa Cruz Biotechnology, Santa Cruz, CA, USA) binding to phosphorylated and unphosphorylated MAPK. Proteins were revealed by chemiluminescence (ECL; Amersham) using peroxidase-coupled secondary antibodies.
Promoter studies
Genomic fragments of 1181 bp (position 98 943-100122 in GenBank entry Z86090 [GenBank] , i.e. upstream of Met1) were PCR amplified with primers introducing NcoI and XhoI restriction sites respectively (see supplementary Table 6) and directionally cloned into pGL3-Basic Luciferase reporter vector (Promega, Mannheim, Germany). PC12 rat pheochromocytoma cells were purchased from Deutsche Sammlung von Mikroorganismen und Zellkulturen (Braunschweig, Germany) and cultured in RPMI 1640 with 10% horse serum and 5% FCS (Biochrom, Germany) in collagen-coated flasks. Transfections and luciferase assays were performed and calculated as described (37), except for electroporation which was carried out with 3 x 106 cells in 0.15 ml and induction for which 1 x 106 cells/ml per well were cultured in the absence or presence of 500 µM dibutyryl-cAMP (Sigma) for 24 h.
| Results |
|---|
|
|
|---|
Mutation screen
SSCP screen of the MCHR1 coding sequence (CDS) in 215 extremely obese children and adolescents (subgroup of study group 1a) and 230 underweight students (study group 1b) revealed 11 infrequent (allele frequencies < 1%; Table 1
) sequence variations leading to seven missense aa exchanges, six non-conservative and four silent sequence variations. In addition, two known SNPs were found (Table 1
): rs133072 (G > A) introduces the nonconservative aa exchange D32N in the N-terminus of the 422 aa receptor protein whereas SNP rs133073 (T > C) is silent (N82). These SNPs are in tight linkage disequilibrium (D' = 0.989, P < 0.001).
Infrequent sequence variations
All infrequent variations were detected in the heterozygous state. Five of the six non-conservative substitutions affect aa that are conserved between humans, mice and rats (D28V, T305M, R317Q, P377S and T411M). Three obese adolescents each harboured two missense variations (D28V and T411M). Since two of them inherited both variations from their mothers (DNA from the parents of the third patient was not available), these two variations form a haplotype [D28V; T411M]. Apart from the common missense variation D32N (rs133072) only one other missense variation (D28V) was detected in both extremely obese (one female: age 16.6 years; BMI 35.83 kg/m2; two males: age 18.83 years; BMI 44.56 kg/m2 and age 15.33 years; BMI 41.33 kg/m) and underweight individuals (age 22.80 years; BMI 19.15 kg/m2 and age 23.16 years; BMI 16.80 kg/m2). Four infrequent missense variations occurred only in extremely obese individuals (BMI
99th percentile): T25M (one male: age 17.13 years; BMI 39.26 kg/m2), T305M (one female: age 13.11 years; BMI 43.32 kg/m2), R317Q (one male: age 15.33 years; BMI 43.24 kg/m2), T411M (one female: age 16.6 years; BMI 35.83 kg/m2; two males: age 18.83 years; BMI 44.56 kg/m2 and age 15.33 years; BMI 41.33 kg/m2). Two infrequent missense variations (R210H and P377S; both affect conserved aa) were identified in underweight students only (R210H: one male: age 21.87 years; BMI 19.15 kg/m2; P377S: one female: age 23.40 years; BMI 16.54 kg/m2). For three of the obese carriers of non-conservative missense variations at conserved aa positions, DNA of one or both parents was also available and genotyped. Onlyone of the parents with such a variation had a BMI in the normal range (age 34.68 years; BMI 22.58 kg/m2; one of the mothers harbouring the [D28V; T411M] haplotype). The other mother with the [D28V; T411M] haplotype was obese (age 46.51 years; BMI 35.21 kg/m2). The father with the non-conservative missense variation T305M was extremely obese (age 41.86 years; BMI 36.11 kg/m2); a sib of the respective index patient with the same variation was also extremely obese (age 10.69 years; BMI 25.10 kg/m2). The father who transmitted the R317Q variation was also extremely obese (age 42.83 years; BMI 44.66 kg/m2).
Common polymorphisms; initial and confirmatory association studies
Comparison of allele and genotype frequencies of both SNPs, rs133072 and rs133073, revealed initial association with obesity in 215 obese probands (subsample of study group 1a) versus 230 underweight probands (study group 1b; all two-sided P values < 0.04), which was independently confirmed upon comparison of an independent obesity sample comprising 504 individuals (the remaining 405 obese probands of study group 1a and the 99 subjects of study group 1d) with 96 normal weight probands (study group 1c; all one-sided P values < 0.03). The final comparison based on all 719 obese (study groups 1a and 1d) and 326 non-obese subjects (study groups 1b and 1c) substantiated association of the A-allele of rs133072 (39.1% in obese and 32.4% in controls) and of the C-allele of rs133073 (41.0% in obese and 33.6% in controls) with obesity (two-sided P values < 0.003; Table 2
).
Common polymorphisms; TDTs
In light of the confirmed association, a TDT based on 525 obese children (study group 1a), all of whom had been included in either the initial or confirmatory association study, and their 1050 parents was performed. Preferential transmissions of the alleles that were shown to be associated with obesity (A-allele of rs133072 and C-allele of rs133073) were observed (two-sided P values < 0.0003; Table 2
). The transmission rates for the A-allele and the C-allele in the trio sample were 58.5% and 58.3% respectively.
Attempts to re-confirm association/transmission disequilibrium in additional study groups
Pursuing the hypothesis that the A-allele of SNP rs133072 is associated with obesity, we genotyped this SNP in additional German and non-German study groups.
50th percentile were negative; the corresponding one-sided P value was 0.26 (A-allele 37.2% in obese and 36.0% in controls). More detailed analyses based on subdivisions for age and BMI decile revealed trends for different allele frequencies between those 124 subjects with BMI >90th percentile and the 633 subjects with a BMI
50th percentile in probands
40 years (A-allele 40.3% in obese and 36.3% in controls; one-sided P value 0.12; Table 2
50th percentile) were detected. In contrast to our hypothesis and the finding in study group 3a, the frequency of the A-allele was higher in controls than in obese individuals (34.0% in obese and 35.2% in controls; one-sided P value = 1; Table 2Identification of additional SNPs, definition of haplotypes and haplotype-based TDTs
To identify additional SNPs which are in linkage disequilibrium with rs133072 we re-sequenced the MCHR1 transcribed region according to AF490537
[GenBank]
as well as 9807 bp of its 5' region and the intron of 1214 bp in nine obese children and adolescents (subgroup of study group 1a) who were homozygous for A at SNP rs133072 and who had received at least one A-allele from a heterozygous parent, as well as in ten obese children and adolescents homozygous for G. This selection was based on the assumption that these two genotype-based groups should be most divergent for functionally relevant SNPs. Twenty SNPs were identified, 12 of which had been reported previously in dbSNP (http://www.ncbi.nlm.nih.gov/SNP/; rs86708, rs133062, rs2032512, rs133063, rs133064, rs133068, rs133069, rs133070, rs133072, rs133073, rs133074 and rs3087592). The eight newly detected SNPs were submitted to the Human Genome Variation database (http://hgvbase.cgb.ki.se/; 001745613, 001745614, 001745615, 001745616, 001745617, 001745618, 001745619 and 001745620; Fig. 1a
).
All nine patients who were selected for homozygosity of the A-allele of rs133072 were also homozygous at 15 of the 20 SNPs; the ten obese individuals selected for homozygosity of the G-allele were homozygous for all SNPs except for 001745613, 001745616 and 001745620. This indicated the existence of two ancestral haplotypes differing in at least 13 sites spanning the re-sequenced region of 13 378 bp.
TDTs for 16 SNPs (including SNPs rs133072 and rs133073; Table 3
) were performed in 61 trios (including the 19 individuals selected for homozygosity at rs133072) and/or in the total of 525 trios (subgroup of study group 1a). All but one of these SNPs revealed transmission rates greater than 60%.
|
Functional in vitro studies
Upon heterologous expression in COS-7 cells, D32- as well as N32-MCHR1 were correctly inserted into the cell membrane and were detectable in equal amounts by cell surface ELISA with 92.5 ± 18.0% relative expression of N32- vs D32-MCHR1 (means ± S.E.M., n
3, optical density for D32-MCHR1 = 100%). Following stimulation with MCH, both receptor proteins were functionally active as indicated by dose-dependent IP accumulation (Fig. 2
), inhibition of isoproterenol-induced cAMP formation (see supplementary Fig. 1
), and phosphorylation of p42/44 MAPK (see supplementary Fig. 2
). Yet, with regard to IP formation no significant functional differences were detectable (Fig. 2
), since both receptor variants showed equal basal activity and ligand potency (EC50 values 1.74 ± 0.45 nM and 3.37 ± 0.97 nM for N32- and D32-MCHR1 respectively). Similar results were obtained using HEK293 cells (data not shown).
|
| Discussion |
|---|
|
|
|---|
A crucial requirement for positive association studies is confirmation in independent samples (38, 39). To substantiate our initially identified and independently confirmed positive association results we determined SNP rs133072 genotypes in five additional study groups (study groups 26) comprised of obese children and adults and controls from within and outside of Germany.
The comparison of allele and genotype frequencies in the Berlin sample comprised of children (study group 2; Table 9) revealed trends similar to those observed in the Marburg sample, but the results were not significant. There was no association in an epidemiological German sample (study group 3a, KORA S4) comprising 393 subjects with a BMI in the highest decile and 2052 subjects with a BMI
50th percentile (Table 2
). However, the initial positive finding was supported upon inclusion of younger probands only (
40 years). In a second German epidemiological sample (study group 3b, KORA S3) there were no significant differences between allele and genotype frequencies of SNP rs133072 in obese and control individuals (Table 2
).
In an adult French sample (study group 5), the A-allele of rs133072 was slightly more frequent in obese individuals but there was no significant association (Table 2
). In contrast, the A-allele of rs133072 was less frequent in obese Danish probands than in controls (study group 4; Table 2
). The TDT was negative in American trios (study group 6; Table 2
).
These findings agree with Gibson et al. (16) who did not detect association of SNPs rs133072 or rs133073 with any obesity-related trait in adult white individuals (age 4065 years) with type 2 diabetes mellitus and related metabolic disorders.
As our initial results were positive, we analysed the MCHR1 in more depth, both genetically and functionally. We re-sequenced ~13.4 kb encompassing the MCHR1. Within the re-sequenced segment we identified two ancestral haplotypes which were deduced from 13 SNPs (see supplementary Table 7). Each of the ancestral haplotypes harbours seven additional SNPs.
In addition to SNPs rs133072 and rs133073 we detected 11 infrequent sequence variations in the CDS of which five non-conservative missense variations affect aa conserved between humans, mice and rats (Table 1
). There was only one variation (D28V) detected in both lean and obese subjects. Because all other variations each occurred in either cases or controls, a role in weight regulation cannot be excluded but remains to be elucidated. We have not characterized the detected infrequent missense variations functionally, so that their role in weight regulation is currently unclear.
Two of the detected non-conservative missense variations (R317Q and T305M) co-segregated with obesity. R317Q was identified in an extremely obese female adolescent and her extremely obese father but not in her extremely obese mother (BMI 40.82 kg/m2;
99th percentile). Gibson et al. (16) identified this variation in one severely obese individual, and it co-segregated with obesity across two generations but there was no evidence for a functional relevance. The non-conservative missense variation T305M was detected in one extremely obese female adolescent, her extremely obese sister and her extremely obese father but not in her slightly obese mother (BMI 25.56 kg/m2; 80th percentile). Recently, the functional relevance of a number of aa variations of MCHR1 was determined; however, none of the studied variations or regions from these studies overlap with the variations identified in this study (4043).
Altered function of N32-MCHR1, i.e. altered IP formation, inhibition of cAMP formation and p42/44 MAPK activation compared with wild-type D32-MCHR1 was not detectable. Nonetheless, since MCHR1 couples to multiple G proteins and thus instigates a plethora of signalling pathways upon stimulation (44, 45), differences in receptor signalling apart from investigated parameters may exist. Also, variations in receptor function might be subtle and not easily discernible with current experimental setups.
In reporter gene assays, we found that MCHR1 promoter fragments of ~1.2 kb comprising either the haplotype associated with obesity or the non-relevant haplotype do not cause differences in either basal or cAMP-induced luciferase expression.
We have concluded that, in accordance with criteria for a solid association study in common disease (38, 39, 46), we have initially detected association for a highly plausible candidate gene (allele) with obesity, which we were able to independently confirm in a second case-control study. These results are supported by family-based studies. However, we cannot exclude the possibility that our original and confirmed association findings represent false positives, since attempts to re-confirm the association and functional studies were negative. It is conceivable that the association applies only to juvenile-onset obesity, conditional on a particular genetic and/or environmental background. In the light of the different results pertaining to MCHR1 SNP rs133072 in six different study groups we consider epigenetic phenomena (47) as one possible explanation. Several genes are involved in the development of obesity and formal genetic studies indicate that both additive and non-additive gene effects contribute to the disease (48, 49). We suppose that a minor effect of the A-allele of rs133072 in MHCR1 can, together with supposedly many other minor effects of sequence variations in different genes and/or their regulatory regions, lead to an increased susceptibility for obesity.
| Acknowledgements |
|---|
| Footnotes |
|---|
|
|
|---|
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
P. Pissios, R. L. Bradley, and E. Maratos-Flier Expanding the Scales: The Multiple Roles of MCH in Regulating Energy Balance and Other Biological Functions Endocr. Rev., October 1, 2006; 27(6): 606 - 620. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||