Diabetes 54:268-274, 2005 © 2005 by the American Diabetes Association, Inc. Linkage of Plasma Adiponectin Levels to 3q27 Explained by Association With Variation in the APM1 Gene
1 University of Maryland School of Medicine, Baltimore, Maryland
We performed a genome-wide linkage scan of plasma adiponectin levels in 569 nondiabetic participants in the Amish Family Diabetes Study. The highest logarithm of odds (LOD) score (2.13; P = 0.0009) occurred on chromosome 3q27 between markers D3S1602 and D3S1580, which flank APM1/ACDC, the adiponectin gene. The APM1 +2019 A/ insertion/deletion polymorphism in the 3' untranslated region (single nucleotide polymorphism [SNP] +2019; deletion allele frequency 0.30 in Amish) showed strong association with adiponectin levels in a dosage-dependent manner in a direction consistent with that reported in previous studies, with deletion heterozygosity increasing adiponectin levels by 1.3 ± 0.5 µg/ml and deletion homozygosity increasing levels by 3.0 ± 0.8 µg/ml (P < 0.0001). Two other SNPs, rs2241766 and rs1501299, showed moderate association. In a subset of 523 subjects genotyped for both SNP +2019 and rs2241766, including the APM1 SNP +2019 genotype as a covariate reduced the linkage signal at 3q27 by 1.26 LOD units (from 2.22 to 0.96) and including both SNPs reduced the signal by 1.51 LOD units (to 0.71). These findings, combined with a two-point LOD score of 2.35 for SNP +2019, provide evidence that variation in APM1 is responsible for linkage of adiponectin levels to 3q27 in the Old Order Amish.
Adiponectin is a fat-secreted hormone involved in insulin sensitivity, free fatty acid oxidation, and inhibition of inflammation (1,2). Circulating adiponectin levels are decreased in obesity (3), diabetes (4), hypertension (5), and dyslipidemia (6), which are components of the metabolic syndrome, or insulin resistance syndrome. Adiponectin is one of a number of adipose tissuesecreted proteins collectively called adipocytokines, the discovery of which in the past decade has radically changed the view of adipose tissue from that of a simple energy storage depot to that of an active endocrine organ. Also known as ACRP30, AdipoQ, and GBP28, adiponectin was discovered in 1995 and 1996 as the most abundantly expressed sequence tag in an adipose mRNA library (79). The chromosomal region 3q27, which contains APM1 (also known as ACDC), the gene encoding adiponectin, has been linked to several features of the metabolic syndrome (10) and diabetes (11). Adiponectin levels have been significantly linked to chromosomes 3 (12,13), 5 and 14 (14), and 9 (15) and suggestively linked to others (14,15). Several investigators have characterized single nucleotide polymorphisms (SNPs) in the adiponectin gene and found one or more of these SNPs to be associated with diabetes (16) and other metabolic syndrome features, as well as with adiponectin levels themselves (1618). We measured adiponectin levels in participants in the Amish Family Diabetes Study (AFDS) to 1) confirm and further define the relation between adiponectin levels and other metabolic syndrome components (19) and 2) identify loci harboring genetic variants that may influence adiponectin levels (and thus insulin resistance and type 2 diabetes). After finding that our strongest linkage signal for adiponectin levels in our genome scan occurred in the chromosomal region containing APM1, the adiponectin structural gene, we performed association and linkage analysis between APM1 variation and adiponectin levels. Our findings suggest that sequence variation in APM1 influences adiponectin levels.
Recruitment for the AFDS began in early 1995 with the goal of identifying genes influencing the risk of type 2 diabetes and related traits. The study protocol was approved by the Institutional Review Board at the University of Maryland School of Medicine, and informed consent was obtained from each study participant. With the help of liaisons from the Old Order Amish community, we identified individuals with type 2 diabetes. These probands and their first- and second-degree relatives and spouses age 18 years or older were invited to participate. All subjects could be connected into a single 14-generation pedigree (20). Among the 691 individuals recruited and included in a 10-cM genomewide scan, there were 70 subjects with type 2 diabetes (25 prevalent cases and 45 incident cases) contained in 43 sibships (including 40 affected subjects in 13 multiplex sibships). Stored fasting plasma from 627 subjects (including 58 with diabetes) of the 691 scanned subjects was available for measuring adiponectin levels. Because of the potential influence of diabetes and its treatment on adiponectin levels, we excluded type 2 diabetic subjects from the analysis, leaving a total of 569 subjects for the linkage analysis. Fasting plasma adiponectin levels were measured in an additional 23 nondiabetic subjects not included in the original genome-wide scan/linkage analysis; these subjects were included in the association analysis (for a total of 592 subjects).
Phenotypes.
Short tandem repeat genotyping.
Candidate gene SNP genotyping. The five SNPs were genotyped using the SNPstream Ultra High Throughput Genotyping System (Beckman Coulter, Fullerton, CA) according to the manufacturers instructions. Briefly, the target genomic sequences containing the SNPs were amplified as part of 12-plex PCRs, which were subjected to an extension reaction using 5'-tagged extension primers and fluorescent dyelabeled terminators. In a thermal cycled extension step, the primers hybridized to the specific amplicons one base adjacent to the SNP site and were extended by one base at the 3' end with a fluorescently labeled nucleotide. These reaction products were transferred to an array plate, where each of the 12 extension products in the multiplex reaction were sorted by hybridization of the unique 5'-tag sequence to its complementary probe immobilized in a mini-array within each well. The plates were imaged by the SNPscope (Beckman Coulter). SNPs were called by comparing the two fluorescent signals. The error rates for these SNPs based on 5% blind replicates were 0%, except for SNPs rs2241766 and rs1501299, which had error rates of 1.2 and 4.8%, respectively. SNP +2019 was genotyped using the PSQ HS 96A Pyrosequencer according to the manufacturers methods (Pyrosequencing, Uppsala, Sweden; www.pyrosequencing.com). Briefly, a PCR product was generated from a primer pair that included one primer covalently coupled to biotin, the biotinylated template was bound to streptavidin-coated Sepharose HP beads, and this mixture was then annealed to a sequencing primer. Stepwise elongation of a sequencing primer strand upon sequential addition of a specified sequence of deoxynucleotide triphosphates and the degradation of nucleotides by apyrase were carried out simultaneously. As the sequencing reaction progressed, the DNA strand was extended and the sequence was determined from the measured signal output of light upon nucleotide incorporation. The resulting peaks in the pyrogram were analyzed using Pyrosequencing software. The error rate for SNP +2019 based on 5% blind replicates was 0%.
Sequencing.
Statistical analysis.
Heritability, linkage, and association analyses were carried out using a variance components methodology. For heritability and linkage analyses, we partitioned variation in the adiponectin trait into components attributable to environmental covariates, the additive effects of genes (i.e., residual heritability), and a specific quantitative trait locus (QTL, the linkage component). These analyses were conducted using maximum likelihood procedures as implemented in the Sequential Oligogenic Linkage Analysis Routines (SOLAR) program (21). The residual heritability was modeled as a function of the expected genetic covariances between relatives, and the QTL effect was modeled as a function of the identity-by-descent probabilities at the marker locus. Allele frequencies were estimated from the data using maximum likelihood methods (21). The hypothesis of linkage was evaluated by the likelihood ratio test, which determines whether the locus-specific effect is significantly >0 (Ho [null hypothesis]: Before association analysis, genotypes were checked for Mendelian consistency using the Pedcheck software package (23) and inconsistencies (<0.5% of genotypes) were resolved or removed. Using variance components, we estimated the effects of the SNP genotype on adiponectin levels, adjusting for sex and sex-specific age and age squared). Genotype was parameterized according to the number of copies of the minor allele the subject carried (genotype = 0 if major allele homozygote, 1 if heterozygote, and 2 if minor allele homozygote). Parameter estimates were obtained conditional on the pedigree structure. For each SNP, the likelihood ratio test was used to compare a model that included the APM1 SNP as a covariate with a model that did not include the SNP. Allele frequencies were estimated by the maximum likelihood estimates conditioned on the pedigree structure as implemented in SOLAR. Linkage disequilibrium (LD) was estimated using an algorithm implemented in the software package Zaplo, which uses both pedigree and LD information to find all possible haplotype configurations of the data under a zero-recombinant assumption (24). To test whether the APM1 association explained the linkage to 3q27, variance components linkage analysis was repeated with the associated APM1 polymorphism genotype as a covariate. The gene-dropping option of Mendel 5.0 (25) was used to simulate the segregation 10,000 times of an unlinked SNP with the same frequency as the APM1 polymorphism. These simulated SNP genotypes were then used as covariates in 10,000 repetitions of the linkage analysis to attempt to assess the significance of the LOD change that occurred with the inclusion of the APM1 variant as a covariate.
Power estimation.
Heritability and genome scan. The heritability of the plasma adiponectin trait, adjusted for sex and sex-specific age and age squared, was estimated at 0.55 ± 0.08. The maximum LOD score in the entire genome scan was 2.13 (P = 0.0009), occurring on chromosome 3q between markers D3S1602 and D3S1580 (Figs. 1 and Table 1). We observed four other regions of suggestive linkage (Table 1): 7q1121, near marker D7S669 (LOD = 1.87, P = 0.0017); 9q34, between markers D9S290 and D9S164 (LOD = 1.54, P = 0.0039); 10p1214, near marker D10S548 (LOD = 1.73, P = 0.0024); and 16p13, near marker D16S407 (LOD = 1.69, P = 0.0026).
Power estimation. The results of the power simulation indicated that the sample provided reasonable power (78%) to detect linkage (at LOD scores >3.0) for a QTL accounting for 30% of the total trait variance. For a QTL accounting for 25% of the variance, power would be 49% to detect a LOD score >3.0 and 80% for a LOD score >2.0. The power was much lower for a QTL accounting for only 20% of the variance (27% for a LOD score >3.0 and 57% for a LOD score >2.0).
Candidate gene analysis.
Including the additive effect of SNP +2019 in the variance components model decreased the variance in adiponectin levels by 2.2%, and including both SNP +2019 and rs2241766 decreased the variance by 3.2% from baseline. The combination of all six SNPs resulted in a decrease of 4.4% from baseline. As shown in Table 3, the three associated SNPs are in strong LD with one another, as estimated by pairwise D' coefficients (0.810.96), but weakly predictive of one another (because of differing allele frequencies), as indicated by low to moderate pairwise R2 estimates (0.020.61). The rs2241766 G allele, which is much less common than the SNP +2019-delA allele (0.08 vs. 0.30), is found almost exclusively in one of the haplotypes containing the SNP +2019-delA allele (Table 4). The rs1501299-T allele is found on the majority of rs2241766-SNP2019 T-delA haplotypes and is absent from the G-delA haplotype. Taken with the association findings discussed in the preceding paragraph, the association of rs2241766 with adiponectin levels may be an independent effect or may be an enhancement of the SNP +2019 effect. Unfortunately, these two possibilities cannot be distinguished by the association data. In contrast, the modest rs1501299 association is most likely driven by its LD with its flanking SNPs.
To assess whether the association with APM1 polymorphisms explained the linkage at 3q27, we reran the genomewide scan using the SNP +2019 and/or rs2241766 genotypes as covariates. When SNP +2019 was included in the model, the LOD score decreased 1.26 units from 2.22 to 0.96, corresponding to an 18.2-fold (i.e., 101.26-fold) decrease in the evidence for linkage (LOD scores computed on the sample containing 523 subjects genotyped for SNP +2019 and rs2241766) (Fig. 2). When SNP rs2241766 was added to the base linkage model, there was a 1.03-unit drop in the LOD score, and when both SNP +2019 and rs2241766 were added simultaneously to the model, there was a 1.51-unit drop in the LOD score, for a total 32.4-fold decrease in the evidence for linkage. The two-SNP model provided a significantly improved fit over the model including SNP +2019 only (P = 0.024).
Assessing whether SNP +2019 accounted for a "significant" proportion of the linkage was not straightforward. To gain some insights into this issue, we performed a simulation in which we dropped an unlinked SNP of the same frequency (0.30) 10,000 times through a 14-generation pedigree connecting virtually all subjects. We then repeated the linkage analysis using the simulated SNP as a covariate. None of the replicates resulted in a reduction as great as 1.26 LOD units, suggesting that association with APM1 SNP +2019, and possibly rs2241766 and/or nearby APM1 SNPs, is quite likely the cause of the linkage to 3q27. Furthermore, the two-point LOD score for APM1 SNP +2019 was 2.35, whereas those for the other polymorphisms ranged from 0.01 (rs1501299) to 1.73 (rs2241766, which was also the next most significantly associated polymorphism). Because we did not initially sequence APM1 in an Amish sample, but rather selected polymorphisms for genotyping based on published information, we followed up the positive association results with direct sequencing analysis of the coding region (the 3' end of exon 2 and the 5' end of exon 3) to determine whether our positive associations might be tracking a potentially functional coding variant. The sample for sequencing consisted of a set of 12 relatively unrelated married couples (24 individuals, 48 alleles) in upper generations in which at least one member was heterozygous for the G allele at rs2241766. We did not detect any variants in the coding region that would explain the association.
Although we were unable to detect chromosomal loci linked to adiponectin levels significant at the LOD 3.0 level, five regions (3q27, 7q1121, 9q34, 10p12-p14, and 16p13) provided evidence of suggestive linkage. The highest LOD score, 2.13, was observed in the region containing APM1, the adiponectin structural gene. Although APM1 might seem an obvious candidate gene for control of adiponectin levels, a genome-wide scan for adiponectin levels in Pima Indians showed no evidence for linkage to this region (although suggestive linkage was observed between markers D3S4018 and D3S1769, which map to 3q13) (15). In contrast, strong linkage of adiponectin levels to 3q27 has been recently reported in the IRAS (Insulin Resistance Atherosclerosis Study) Family Study (LOD = 3.84) (13). One interesting finding of that study was that evidence for linkage was provided entirely by the Hispanic-American sample (LOD = 5.25); no evidence for linkage was seen in the African-American families (LOD = 0.0) (13). Weak linkage to 3q27 (LOD = 1.3) was found in another genomewide scan of adiponectin levels in Caucasians, where much higher signals on chromosomes 5 and 14 were seen (14). Linkage in the 3q2224 region (LOD = 3.2), about 40 cM centromeric to 3q27, has been observed in Mexican Americans (12). Based on these findings, it is possible that variation in APM1 is an important determinant of circulating adiponectin levels in some populations, but not others. Our linkages for adiponectin on chromosomes 7q1121, 9q34, 10p1214, and 16p13 are close to linkages we observed in the Amish for other obesity-related traits, such as total serum cholesterol, BMI-adjusted leptin, and BMI (Table 1). Although some of these traits are somewhat correlated with adiponectin levels (19), they provide additional support for the idea that these regions are true linkages rather than false-positives. The 10p region has also been reported by others to be close to linkages for adiponectin (27), obesity (28), and BMI (29), thus providing additional evidence that a gene influencing obesity and related traits resides in this region. To confirm our suspicion that the linkage signal on 3q27 resulted from variation in the adiponectin gene itself, we performed association analysis with six previously reported SNPs (Table 2). SNP +2019 showed strong association with adiponectin levels and appeared to explain most of the linkage signal at 3q27, as accounting for this SNP led to a 57% drop in the LOD score, corresponding to an 18.2-fold drop in the evidence for linkage. This insertion/deletion polymorphism has been previously reported to be associated with adiponectin levels (16). Another group has found the corresponding insertion to be in LD with the 45T/276G haplotype previously associated with type 2 diabetes (17) and has suggested that SNP +2019 may be a causal variant via mechanisms involving regulatory elements in the 3' region of the APM1 gene. Based on searches in UTRSite and UTRScan, the polymorphism does not reside within a known functional sequence; the closest known functional element identified is a K-box binding region 371 bp downstream. However, functional elements of the 3' untranslated region (UTR) continue to be elucidated, and the possibility that the variant alters a protein or micro-RNA binding site and/or RNA secondary structure cannot be not ruled out (30). Studies of the effect of the polymorphism on RNA production and stability along with addi-tional characterization of the 3' UTR sequence will be needed to confirm the causative nature of this polymorphism. In addition to association with SNP +2019, there appeared to be a further source of association in rs2241766 independent from the SNP +2019 association, as evidenced by 1) an apparent additive effect between the two polymorphisms that reached statistical significance in the linkage analysis and 2) the persistence of the SNP +2019 association when the sample was restricted to those with the rs2241766 TT genotype. No common variants were detected in a set of 24 relatively unrelated individuals enriched for the presence of the rs2241766 and SNP +2019 minor alleles, increasing the likelihood that these variants directly influence adiponectin levels, although the possibility remains that variants outside the coding region are influencing the phenotype. As in the Amish, the rs1501299-T allele has been associated with higher adiponectin levels in a Japanese population (18). The Amish data also replicate associations of the rs2241766-G, rs1501299-T, and SNP +2019-del alleles with increased adiponectin levels in the French population (16). Unlike the French (16), the Amish did not show association of the rs266729-G allele with decreased adiponectin levels; however, a follow-up analysis of the Amish in which subgroups in the lower and upper 50th BMI percentile were analyzed separately did reveal a recessive association between the rs266729-G allele and lower adiponectin levels in the high BMI (>26.4) group only (P = 0.014). The two-point linkage analysis of the relation of the SNPs themselves to adiponectin levels provides additional evidence that the sequence variation within APM1 explains the 3q27 linkage, as the highest LOD score among the SNPs, 2.35, which was almost identical to the multipoint peak score, was observed at the polymorphism (SNP +2019) with the strongest association result. In summary, we have provided evidence for linkage of adiponectin levels to several distinct regions in the genome in the Amish. Among them, 3q27 contains the structural gene for adiponectin itself, APM1, in which sequence variation appears to be an important determinant of adiponectin levels. Because adiponectin is a key regulator of insulin sensitivity and other metabolic syndromerelated traits, we would expect this gene, and likely specific SNPs therein (e.g., SNP +2019, rs2241766) to be important contributors to these traits. Replication of our findings along with further investigation of sequence variation within APM1 and functional analysis of these variants will be required to further discern the role of APM1 as a susceptibility gene for obesity, type 2 diabetes, and the metabolic syndrome.
Electronic database information for Internet-based resources used in this study: 1) NCBI UniSTS: http://www.ncbi.nih.gov/entrez/query.fcgi?db=unists; 2) UCSC Human Genome Browser Gateway: http://genome.ucsc.edu/cgi-bin/hgGateway; and 3) UTRScan/UTRSite: http://bighost.area.ba.cnr.it/BIG/UTRHome/.
This study was supported in part by a research grant from GlaxoWellcome and Axys Pharmaceuticals; National Institutes of Health Grants DK-54361 and DK-02673; an American Diabetes Association Research Award to A.R.S.; an American Heart Association Mid-Atlantic Predoctoral Fellowship to T.I.P.; University of Maryland General Clinical Research Center Grant M01-RR-165001; and the General Clinical Research Centers Program, National Center for Research Resources, National Institutes of Health. The authors thank Wendy Warren, Mary Ann Drolet, Denis Massey, Mary Morrissey, Janet Reedy, and our Amish liaisons for their energetic efforts in study subject recruitment and characterization; Drs. Alejandro Schäffer and Richa Agarwala for assistance in pedigree construction; and Jack Shelton for sequencing. This work would not have been possible without the outstanding cooperation of the Amish community. Address correspondence and reprint requests to Toni I. Pollin, PhD, Division of Endocrinology, Diabetes, and Nutrition, University of Maryland School of Medicine, 660 W. Redwood St., Rm. 478, Baltimore, MD 21201 E-mail: tpollin{at}medicine.umaryland.edu Received for publication July 1, 2004 and accepted in revised form September 15, 2004
Abbreviations: AFDS, Amish Family Diabetes Study; LD, linkage disequilibrium; QTL, quantitative trait locus; SNP, single nucleotide polymorphism; SOLAR, Sequential Oligogenic Linkage Analysis Routines; UTR, untranslated region
This article has been cited by other articles:
|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||