1 Department of Molecular Cell Biology, Leiden University Medical Center, Leiden, the Netherlands
2 Steno Diabetes Center and Hagedorn Research Institute, Gentofte, Denmark
3 Department of Genetic Epidemiology and Biostatistics, Erasmus University Medical Center, Rotterdam, the Netherlands
4 Institute for Research in Extramural Medicine, Free University Medical Center, Amsterdam, the Netherlands
5 Department of Internal Medicine, Erasmus University Medical Center, Rotterdam, the Netherlands
6 Research Center for Prevention and Health, Copenhagen County, Glostrup University Hospital, Glostrup, Denmark
7 Faculty of Health Science, Aarhus University, Aarhus, Denmark
![]() |
ABSTRACT |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
Previously, we have shown that an A3243G mutation in the mitochondrial DNAencoded tRNALeu(UUR)gene is associated with maternally inherited diabetes and deafness (1,2). Carriers of the 3243 mutation have a reduced glucose-stimulated insulin secretion. Biochemical studies indicate that mutant mitochondrial tRNALeucine is less efficiently aminoacylated in comparison with the wild-type tRNA (2). Furthermore, it has been shown that there is an altered balance between mitochondrial and nuclear-encoded proteins in mutant mitochondria, resulting in mitochondrial dysfunction (2).
In this study, we have examined whether variation in the nuclear-encoded mitochondrial leucyl tRNA synthetase gene (LARS2) is associated with type 2 diabetes. The product of the LARS2 gene is the leucyl tRNA synthetase protein (LeuRS, EC 6.1.1.4), which catalyzes the charging of tRNALeu(UUR) with leucine, an essential step in protein synthesis. Subjects (n = 7,836) participating in this study were recruited from four independent population samples in the Netherlands and Denmark to allow independent replication of association results.
![]() |
RESEARCH DESIGN AND METHODS |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
Furthermore, we genotyped two case-control study samples from Denmark, DK1 and DK2. DK1 was a group of unrelated type 2 diabetic patients recruited from the Steno Diabetes Center and a group of unrelated normal glucosetolerant (NGT) subjects sampled at random through public registers at the Steno Diabetes Center and the Research Centre for Prevention and Health (5). In the group of diabetic patients (n = 706, 48% male), the age was 59 ± 10 years and 57 ± 10 years for NGT participants (n = 514, 46% male) (5). DK2 consisted of unrelated type 2 diabetic patients recruited from the Steno Diabetes Center and the Research Centre for Prevention and Health (n = 654, aged 54 ± 10 years, 42% male) and a group of unrelated NGT subjects sampled from the prospective Inter99 study at the Research Centre for Prevention and Health (n = 4,501, aged 45 ± 8 years, 47% male) (6,7). All control subjects participating in the Danish studies underwent a fasting OGTT according to World Health Organization criteria.
Type 2 diabetic subjects were significantly older in all populations (all P 0.01). Questionnaires were used to obtain other relevant information. All participants were Caucasian whites by self-report. Informed written consent was obtained from all subjects before participation. The study was approved by the appropriate medical ethical committees and was in accordance with the principles of the Declaration of Helsinki.
LARS2 sequencing.
The LARS2 gene is mapped to the chromosome 3p21.3 region and spans 160 kb. Its 22 exons encode a 903 aa protein (LocuslinkID 23395). Total RNA was isolated from leukocytes and/or pancreas tissue (n = 7) from 25 Dutch type 2 diabetic subjects using standard procedures. First-strand cDNA was made from these RNA samples using a first-strand cDNA synthesis kit (Invitrogen, Carlsbad, CA). The LARS2 gene was amplified in two overlapping segments (primer sequences and assay conditions available on request). The resulting segments encompassing 2,905 bp of the coding region, including 159 bp of the promoter region, were subsequently sequenced on an automated sequencer (Applied Biosystems, Foster City, CA). Haplotypes of the identified single nucleotide polymorphisms (SNPs) (allele frequency, >2.5%) were calculated using the Phase Program (8).
Genotyping.
Two SNPs in the coding region of the LARS2 gene were studied in further detail for association with type 2 diabetes and related quantitative variables like circulating glucose and insulin levels during an OGTT. The first SNP (g/a) was located at position 109 in the promoter area. A 197-bp fragment surrounding the SNP was amplified using a mismatch primer. The PCR fragment was cleaved with the restriction enzyme MspA1I, resulting in fragments of 178 and 19 bp in the case of the wild-type sequence and 197 bp in the case of the mutant sequence after separation on a 4% agarose gel.
The second SNP is a coding SNP located in exon 10. The C to A variant changes a histidine at position 324 to glutamine (H324Q). Two different assays were developed for genotyping this SNP. The first method was a restriction fragmentlength polymorphismbased method. The genomic region surrounding the SNP was amplified using a mismatch primer introducing an MspA1I site when the C allele is present. Primer sequences for both assays are available on request. The second genotyping method for this SNP was based on Taqman SNP genotyping technology (Applied Biosystems, Foster City, CA). Primer and probe design was performed by the manufacturer (sequences available on request), and reactions were done according to the manufacturers protocol. This method was applied to all samples tested in this study including sufficient replication and control samples. A replication sample of 1,500 Dutch samples using the restriction fragmentlength polymorphismbased method did not identify any mismatches. Furthermore, the accuracy of both measures was further confirmed by direct sequencing of an additional 40 samples.
Overexpression and purification of wild-type and mutant leucyl tRNA synthetase.
The wild-type mature enzyme of 864 amino acid residues (39LeuRS) was overexpressed in the BL21-Codonplus(DE3)RIL strain as described (9). A human cDNA clone containing this His-tagged NH2-terminally truncated LeuRS was obtained from M.P. King (10). The enzyme was purified using nickel affinity chromatography as described (9). Furthermore, we have made a construct with the H324Q polymorphism for production and purification of the mutant enzyme (cloning details available on request).
Aminoacylation and editing reactions.
Aminoacylation properties of the wt and H324Q mutant enzymes were determined by measuring [3H]-leucine incorporation according to a previously described method (9), using either total bovine tRNA and/or human mitochondrial RNA as the tRNA source. Furthermore, we have analyzed the mischarging rate of both enzymes using several noncognate [3H]-labeled amino acids including [3H]-isoleucine (11). A priori power calculations have shown that we should be able to detect differences >1015%.
Statistical analyses.
Hardy-Weinberg equilibrium was tested in all cohorts before further analysis. Differences in allele frequencies were tested by Fishers exact tests. Since testing of the homogeneity of the population-specific odds ratios revealed no significant differences, a common odds ratio was calculated using a Mantel-Haenszel test. Differences in clinical variables were tested by linear regression analysis, with adjustment for age, sex, and BMI. Variables were log transformed before analysis, if necessary. The statistical software packages SPSS 10.0 (SPSS, Chicago, IL) and StatXact 6.0 (Cytel, Cambridge, MA) were used.
![]() |
RESULTS |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
The H324Q variant was examined in all four populations. Allele frequencies for each of the samples are given in Table 1. Significant association between the gene variant and type 2 diabetes was present in the Hoorn cohort (P = 0.04, Table 1). Furthermore, there was a nonsignificant but identical trend toward association in all other cohorts (Table 1). Analysis of the genotype distributions and adjustment for differences in age and sex yielded very similar results (supplementary Table 1 [available at http://diabetes.diabetesjournals.org).
|
|
![]() |
DISCUSSION |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
The risk associated with this gene variant is in the same order of magnitude as shown for the three well-replicated disease-associated genes: KIR6.2, peroxisome proliferatoractivated receptor , and the Calpain10 haplotypes (odds ratio 1.21.4) (12), suggesting again that the contribution of predisposing genes per gene is limited and most likely only a combination of several diabetes risk genes in combination with a detrimental lifestyle results in glucose intolerance. This further implies that only studies with sufficient statistical power will be able to track down and confirm association between such variants and type 2 diabetes.
It is difficult to predict whether this His-to-Gln mutation, present at the surface of the editing domain, affects its function. The editing domain is involved in the repair of mischarged tRNAs, thereby keeping error rates in the synthesis of proteins low. Functional in vitro tests examining aminoacylation and mischarging of the tRNA leucine have not shown major differences. However, additional in vitro and in vivo studies are needed to further elucidate the precise functional impact of this protein variant. Since we have studied only two SNPs in detail, it might be that the H324Q variant is merely reflecting linkage disequilibrium with another unknown causative variant in this gene or a nearby-located gene at the same locus. Further, more detailed studies of this gene locus will be needed to replicate and fully examine the nature of the observed association between the H324Q variant and type 2 diabetes.
In conclusion, we report the association of an amino acid polymorphism in the leucyl tRNA synthetase gene and type 2 diabetes in a combined analysis of four independent study populations in the Netherlands and Denmark. To our knowledge, this is the first report showing association between variation in an aminoacyl tRNA synthetase gene and disease, and it further highlights the importance of correct mitochondrial function in glucose homeostasis (1,1315).
![]() |
ACKNOWLEDGMENTS |
---|
The authors thank all participants for their kind cooperation. M.P. King (Philadelphia, PA) is kindly acknowledged for the gift of the LeuRS expression clone. Annemette Forman, Inge Lise Wantzin, and Marianne Stendal are thanked for dedicated and careful technical assistance and Grete Lademann for secretarial support.
![]() |
FOOTNOTES |
---|
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Address correspondence and reprint requests to L.M. t Hart, PhD, Leiden University Medical Center, Department of Molecular Cell Biology, Wassenaarseweg 72, 2333 AL Leiden, Netherlands. E-mail: l.m.t_hart{at}lumc.nl
Received for publication November 25, 2004 and accepted in revised form March 21, 2005
NGT, normal glucose tolerant; OGTT, oral glucose tolerance test; SNP, single nucleotide polymorphism
![]() |
REFERENCES |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Diabetes | Diabetes Care | Clinical Diabetes | Diabetes Spectrum | DOC News |