Evaluation of the microsomal glutathione S-transferase 3 (MGST3) locus on 1q23 as a Type 2 diabetes susceptibility gene in Pima Indians

Elevation of plasma glucose concentration may induce generation of oxygen-free radicals, which can play an important role in the progression of diabetes and/or development of its complications. Various glutathione transferases utilize the availability of reduced glutathione for the cellular defense against oxygen-free radicals. One such enzyme is microsomal glutathione S-transferase 3 encoded by MGST3, which maps to chromosome 1q23, a region linked to Type 2 diabetes mellitus (T2DM) in Pima Indians, Caucasian, and Chinese populations. We investigated the MGST3 gene as a potential susceptibility gene for T2DM by screening this locus for single nucleotide polymorphisms (SNPs) in diabetic and non-diabetic Pima Indians. We also measured the skeletal muscle MGST3 mRNA level by Real-Time (RT) PCR and its relationship with insulin action in non-diabetic individuals. We identified 25 diallelic variants, most of which, based on their genotypic concordance, could be divided into three distinct linkage disequilibrium (LD) groups. We genotyped unique representative SNPs in selected diabetic and non-diabetic Pima Indians and found no evidence for association with T2DM. Furthermore, inter-individual variation of skeletal muscle MGST3 mRNA was not correlated with differences in insulin action in non-diabetic subjects. We conclude that alterations of MGST3 are unlikely to contribute to T2DM or differences in insulin sensitivity in the Pima Indians.     

Identification of genetic variation that determines human trehalase activity and its association with type 2 diabetes

A prior linkage scan in Pima Indians identified a putative locus for type two diabetes (T2D) and body mass index (BMI) on chromosome 11q23-25. Association mapping across this region identified single nucleotide polymorphisms (SNPs) in the trehalase gene (TREH) that were associated with T2D. To assess the putative connection between trehalase activity and T2D, we performed a linkage study for trehalase activity in 570 Pima Indians who had measures of trehalase activity. Strong evidence of linkage of plasma trehalase activity (LOD = 7.0) was observed in the TREH locus. Four tag SNPs in TREH were genotyped in these subjects and plasma trehalase activity was highly associated with three SNPs: rs2276064, rs117619140 and rs558907 (p = 2.2 × 10^?11–1.4 × 10^?23), and the fourth SNP, rs10790256, was associated conditionally on these three (p = 2.9 × 10^?7). Together, the four tag SNPs explained 51 % of the variance in plasma trehalase activity and 79 % of the variance attributed to the linked locus. These four tag SNPs were further genotyped in 828 subjects used for association mapping of T2D, and rs558907 was associated with T2D (odds ratio (OR) 1.94, p = 0.002). To assess replication of the T2D association, all four tag SNPs were additionally genotyped in two non-overlapping samples of Native Americans. Rs558907 was reproducibly associated with T2D in 2,942 full-heritage Pima Indians (OR 1.27 p = 0.03) and 3,897 “mixed” heritage Native Americans (OR 1.21, p = 0.03), and the strongest evidence for association came from combining all samples (OR 1.27 p = 1.6 × 10^?4, n = 7,667). However, among 320 longitudinally studied subjects, measures of trehalase activity from a non-diabetic exam did not predict those who would eventually develop diabetes versus those who would remain non-diabetic (hazard ratio 0.94 per SD of trehalase activity, p = 0.29). We conclude that variants in TREH control trehalase activity, and although one of these variants is also reproducibly associated with T2D, it is likely that the effect of the SNP on risk of T2D occurs by a mechanism different than affecting trehalase activity. Alternatively, TREH variants may be tagging a nearby T2D locus.     

Association of single-nucleotide polymorphisms in the suppressor of cytokine signaling 2 (SOCS2) gene with type 2 diabetes in the Japanese

Several previous linkage scans in type 2 diabetes (T2D) families indicated a putative susceptibility locus on chromosome 12q15-q22, while the underlying gene for T2D has not yet been identified. We performed a region-wide association analysis on 12q15-q22, using a dense set of >500 single-nucleotide polymorphisms (SNPs), in 1492 unrelated Japanese individuals enrolled in this study. We identified an association between T2D and a haplotype block spanning 13.6 kb of genomic DNA that includes the entire SOCS2 gene. Evolutionary-based haplotype analysis of haplotype-tagging SNPs followed by a "sliding window" haplotypic analysis indicated SNPs that mapped to the 5' region of the SOCS2gene to be associated with T2D with high statistical significance. The SOCS2 gene was expressed ubiquitously in human and murine tissues, including pancreatic beta-cell lines. Adenovirus-mediated expression of the SOCS2 gene in MIN6 cells or isolated rat islets significantly suppressed glucose-stimulated insulin secretion. Our data indicate that SOCS2 may play a role in susceptibility to T2D in the Japanese.     

Association of a promoter variant in the inducible cyclooxygenase-2 gene (PTGS2) with type 2 diabetes mellitus in Pima Indians

Recent studies have suggested that prostaglandin-endoperoxide synthase-2 (PTGS2), also known as cyclo-oxygenase 2, plays an etiological role in the development of type 2 diabetes mellitus (T2DM). PTGS2 generates prostaglandins, which negatively modulate glucose-stimulated insulin secretion, and functions as a mediator of the inflammatory response, which is associated with decreased insulin sensitivity. Moreover, the gene encoding this enzyme, PTGS2, is located on 1q25.2, a region that has been linked with early onset T2DM in Pima Indians. To determine the possible role played by PTGS2 in modulating susceptibility to T2DM, we screened approximately 7.0 kb of the gene, corresponding to the promoter, coding sequence, and flanking exon-intron boundaries, and identified five variants, including three single nucleotide polymorphisms (SNPs) in the promoter, one intronic SNP, and one in the 3' untranslated region. With the exception of one rare promoter SNP (minor allele frequency <0.03), all SNPs were typed in ~1000 Pima Indians. The range of frequencies for the more common alleles was 0.65–0.88, and we found substantial linkage disequilibrium between all PTGS2 SNP pairs (D0.95). Variant alleles at two markers, rs20417 and rs2066826, which are located in the promoter and intron 6, respectively, were in strong linkage disequilibrium with each other (D=0.97) and were associated with a higher prevalence of T2DM. For marker rs20417, individuals with the variant CC genotype had a 30% higher T2DM prevalence compared with subjects with the GG genotype (odds ratio=1.6 per copy of C allele; P=0.01). The variant C allele of rs20417 has been associated with decreased PTGS2 promoter activity, thereby suggesting a possible biological consequence attributable to this polymorphism. These findings indicate that genetic variants in PTGS2 may play a role in mediating susceptibility to T2DM in Pima Indians and are consistent with the hypothesis that chronic inflammation may contribute to the development of T2DM in some individuals.     

Genomics of type I diabetes mellitus and its late complications

In ethnic Russians, MHC (HLA) was shown to be the major locus determining the predisposition to type 1 diabetes mellitus (T1DM). To map the regions linked to T1DM, families with concordant or discordant sib pairs were selected from the Russian population of Moscow. With these families, linkage to T1DM was demonstrated for CTLA4 (IDDM12, 2q32.1-q33), which codes for a T-cell surface antigen, and PDCD2 (IDDM8, 6q25-q27), which is homologous to the mouse programmed cell death activator gene. With polymorphic microsatellites, regions 3q21-q25 (IDDM9) and 10p12.2 (IDDM10) were also linked to T1DM. Case/control and family studies of the polymorphic markers from region 11p13 revealed a new T1DM-associated locus in the vicinity of the catalase gene (CAT); linkage to this locus was not reported earlier for other populations. Diabetic polyneuropathy (DPN) proved to be associated with single-nucleotide polymorphisms Ala(-9)Val (SOD2), Arg213Gly (SOD3), and T(-262)C (CAT) and with a polymorphic microsatellite of the NOS2 promoter. Hence oxidative stress, which results from hyperglycemia, increased mitochondrial production of superoxide radicals, and insufficient activities of antioxidative enzymes, was assumed to play an important part in DPN development in T1DM. Diabetic nephropathy (DN) showed no association with the antioxidative enzyme genes. However, the association was observed for the insertion/deletion (I/D) polymorphism of ACE and the ecNOS34a/4b polymorphism of NOS3, two genes involved in controlling vascular tonicity, and for the I/D polymorphism of APOB and the epsilon 2/epsilon 3/epsilon 4 polymorphism of APOE, two genes involved in lipid transport. In addition, polymorphic microsatellites of chromosome 3q21-q25 proved to be closely associated with DN. The tightest association was established for D3S1550, carriers of allele 12 or genotype 12/14 having high risk of DN (OR = 4.85 and 6.25, respectively). Region 3q21-q25 was assumed to contain a major gene determining DN development, while the other DN-associated genes mostly affect the progression of DN.     

Conserved physical linkage of GnRH-R and RBM8 in the medaka and human genomes

Candidate genes for human type II gonadotropin-releasing hormone receptor (GnRH-RII) reside on two separate loci, 1q12-q21 and 14q21-23, yet neither locus generates functional GnRH-RII. Instead, their opposite DNA strands encode functional RNA-binding motif protein 8 (RBM8s), which is also encoded by another locus, 5q13-q14. To elucidate the mechanism through which such multiple human GnRH-RII/RBM8 loci arose, here we have defined an RBM8 locus in a comparative model species, the medaka Oryzias latipes. The medaka RBM8, which exists as a single copy gene, is linked to, but does not overlap with, GnRH-R2 on linkage group (LG) 16, demonstrating the ancient origin of the physical linkage between GnRH-R and RBM8. The medaka LG 16 contains orthologous segments to the human chromosome 1 and therefore the 1q12-q21 locus would be an originating human GnRH-RII/RBM8 segment. Furthermore, like the human RBM8s on 1q12-q21 and 5q13-q14 but not that on 14q21-q23, the medaka RBM8 is a multiexon gene, indicating that the 14q21-q23 and 5q13-q14 loci were generated by retrotransposition and segmental genomic duplication, respectively, of the originating 1q12-q21 locus.     

Consistently replicating locus linked to migraine on 10q22-q23

Here, we present the results of two genome-wide scans in two diverse populations in which a consistent use of recently introduced migraine-phenotyping methods detects and replicates a locus on 10q22-q23, with an additional independent replication. No genetic variants have been convincingly established in migraine, and although several loci have been reported, none of them has been consistently replicated. We employed the three known migraine-phenotyping methods (clinical end diagnosis, latent-class analysis, and trait-component analysis) with robust multiple testing correction in a large sample set of 1675 individuals from 210 migraine families from Finland and Australia. Genome-wide multipoint linkage analysis that used the Kong and Cox exponential model in Finns detected a locus on 10q22-q23 with highly significant evidence of linkage (LOD 7.68 at 103 cM in female-specific analysis). The Australian sample showed a LOD score of 3.50 at the same locus (100 cM), as did the independent Finnish replication study (LOD score 2.41, at 102 cM). In addition, four previously reported loci on 8q21, 14q21, 18q12, and Xp21 were also replicated. A shared-segment analysis of 10q22-q23 linked Finnish families identified a 1.6-9.5 cM segment, centered on 101 cM, which shows in-family homology in 95% of affected Finns. This region was further studied with 1323 SNPs. Although no significant association was observed, four regions warranting follow-up studies were identified. These results support the use of symptomology-based phenotyping in migraine and suggest that the 10q22-q23 locus probably contains one or more migraine susceptibility variants.     

Analysis of PBX1 as a candidate gene for type 2 diabetes mellitus in Pima Indians

The human proto-oncogene PBX1 codes for a homeodomain containing protein that modulates expression of several genes, including those contributing to regulation of insulin action and glucose metabolism. PBX1 is located on chromosome 1q22, a region linked with type 2 diabetes in Pima Indians, Caucasians, and an Old Order Amish population. We have investigated the PBX1 genomic sequence to identify polymorphisms that may contribute to diabetes susceptibility in the Pimas. PBX1 is composed of nine exons spanning approx. 117 kb and is located within 300 kb of microsatellite D1S1677, which marks the peak of linkage to diabetes susceptibility in the Pima Indians. We detected 16 single nucleotide polymorphisms in PBX1 including one causing a glycine to serine substitution at residue 21. Comparison of the frequencies of the polymorphisms between affected and unaffected Pima Indians did not detect any significant differences, indicating that mutations in PBX1 do not explain the linkage of 1q with type 2 diabetes in this population. The genomic structure of PBX1 provides a basis for similar systematic examinations of this candidate locus in other populations in relation to both type 2 diabetes and other metabolic disorders.     

Heterogeneity in gene loci associated with type 2 diabetes on human chromosome 20q13.1

Human chromosome 20q12-q13.1 has been linked to type 2 diabetes mellitus (T2DM) in multiple studies. We screened a 5.795-Mb region for diabetes-related susceptibility genes in a Caucasian cohort of 310 controls and 300 cases with T2DM and end-stage renal disease (ESRD), testing 390 SNPs for association with T2DM-ESRD. The most significant SNPs were found in the perigenic regions: HNF4A (hepatocyte nuclear factor 4alpha), SLC12A5 (potassium-chloride cotransporter member 5), CDH22 (cadherin-like 22), ELMO2 (engulfment and cell motility 2), SLC13A3 (sodium-dependent dicarboxylate transporter member 3), and PREX1 (phosphatidylinositol 3,4,5-triphosphate-dependent RAC exchanger 1). Haplotype analysis found six haplotype blocks globally associated with disease (p<0.05). We replicated the PREX1 SNP association in an independent case-control T2DM population and inferred replication of CDH22, ELMO2, SLC13A3, SLC12A5, and PREX1 using in silico perigenic analysis of two T2DM Genome-Wide Association Study data sets. We found substantial heterogeneity between study results.     

Analysis of the SAG5 locus reveals a distinct genomic organisation in virulent and avirulent strains of Toxoplasma gondii

We have recently characterised, in the virulent strain RH of Toxoplasma gondii, three glycosylphosphatidylinositol-anchored surface antigens related to SAG1 (p30) and encoded by highly homologous, tandemly arrayed genes named SAG5A, SAG5B and SAG5C. In the present study, we compared the genomic organisation of the SAG5 locus in strains belonging to the three major genotypes of T. gondii. Southern blot analysis using a SAG5-specific probe produced two related but distinct hybridisation patterns, one exclusive of genotype I virulent strains, the other shared by avirulent strains of either genotype II or genotype III. To understand the molecular bases of this intergenotypic heterogeneity, we cloned and sequenced the SAG5 locus in the genotype II strain Me49. We found that in this isolate the SAG5B gene is missing, with SAG5A and SAG5C laying contiguously. This genomic arrangement explains the hybridisation profiles observed for all the avirulent strains examined and indicates that the presence of SAG5B is a distinctive trait of genotype I. Furthermore, we identified two novel SAG1-related genes, SAG5D and SAG5E, mapping respectively 1.8 and 4.0 kb upstream of SAG5A. SAG5D is transcribed in tachyzoites and encodes a polypeptide of 362 amino acids sharing 50% identity with SAG5A-C, whereas SAG5E is a transcribed pseudogene. We also evaluated polymorphisms at the SAG5 locus by comparing the coding regions of SAG5A-E from strains representative of the three archetypal genotypes. In agreement with the strict allelic dimorphism of T. gondii, we identified two alleles for SAG5D, whereas SAG5A, SAG5C and SAG5E were found to be three distinct nucleotide variants. The higher intergenotypic polymorphism of SAG5A, SAG5C and SAG5E suggests that these genes underwent a more rapid genetic drift than the other members of the SAG1 family. Finally, we developed a new PCR–restriction fragment length polymorphism method based on the SAG5C gene that is able to discriminate between strains of genotype I, II and III by a single endonuclease digestion.     

Polymorphisms in the trace amine receptor 4 (TRAR4) gene on chromosome 6q23.2 are associated with susceptibility to schizophrenia

Several linkage studies across multiple population groups provide convergent support for a susceptibility locus for schizophrenia--and, more recently, for bipolar disorder--on chromosome 6q13-q26. We genotyped 192 European-ancestry and African American (AA) pedigrees with schizophrenia from samples that previously showed linkage evidence to 6q13-q26, focusing on the MOXD1-STX7-TRARs gene cluster at 6q23.2, which contains a number of prime candidate genes for schizophrenia. Thirty-one screening single-nucleotide polymorphisms (SNPs) were selected, providing a minimum coverage of at least 1 SNP/20 kb. The association observed with rs4305745 (P=.0014) within the TRAR4 (trace amine receptor 4) gene remained significant after correction for multiple testing. Evidence for association was proportionally stronger in the smaller AA sample. We performed database searches and sequenced genomic DNA in a 30-proband subsample to obtain a high-density map of 23 SNPs spanning 21.6 kb of this gene. Single-SNP analyses and also haplotype analyses revealed that rs4305745 and/or two other polymorphisms in perfect linkage disequilibrium (LD) with rs4305745 appear to be the most likely variants underlying the association of the TRAR4 region with schizophrenia. Comparative genomic analyses further revealed that rs4305745 and/or the associated polymorphisms in complete LD with rs4305745 could potentially affect gene expression. Moreover, RT-PCR studies of various human tissues, including brain, confirm that TRAR4 is preferentially expressed in those brain regions that have been implicated in the pathophysiology of schizophrenia. These data provide strong preliminary evidence that TRAR4 is a candidate gene for schizophrenia; replication is currently being attempted in additional clinical samples.     

Genome-wide scan and fine-mapping linkage study of androgenetic alopecia reveals a locus on chromosome 3q26

Androgenetic alopecia (AGA, male pattern baldness) is the most common form of hair loss. The origin of AGA is genetic, with the X chromosome located androgen receptor gene (AR) being the only risk gene identified to date. We present the results of a genome-wide linkage study of 95 families and linkage fine mapping of the 3q21-q29, 11q14-q25, 18p11-q23, and 19p13-q13 regions in an extended sample of 125 families of German descent. The locus with strongest evidence for linkage was mapped to 3q26 with a nonparametric linkage (NPL) score of 3.97 (empirical p value = 0.00055). This is the first step toward the identification of new susceptibility genes in AGA, a process which will provide important insights into the molecular and cellular basis of scalp hair loss.     

Genomic structure of karyopherin α2 (KPNA2) within a low-copy repeat on chromosome 17q23-q24 and mutation analysis in patients with Russell-Silver syndrome

Human karyopherin alpha2 (KPNA2), a member of the karyopherin alpha family, plays a key role in the nuclear import of proteins with a classical nuclear localization signal (NLS). KPNA2, as part of a karyopherin alpha-beta heterodimer, directly binds to the NLS of proteins and functions as an adaptor that binds NLS-containing proteins via karyopherin beta to the nuclear pore complex. The NLS protein-receptor complex is translocated through the pore by an energy-dependent mechanism. Recently, we have identified and mapped the gene for KPNA2 in close proximity to a translocation breakpoint within 17q23-q24 associated with Russell-Silver syndrome (RSS). Therefore, we considered KPNA2 as a positional candidate gene for this heterogeneous disorder. RSS is mainly characterized by pre- and postnatal growth retardation, lateral asymmetry, and other dysmorphic features. Here, we present the genomic organization of the human KPNA2 gene with 11 exons spanning approximately 10 kb on chromosome 17q23-q24. Screening for mutations within all exons and adjacent intronic sequences from 31 unrelated RSS patients revealed three single nucleotide polymorphisms (SNPs) in exons 1, 5, and 7, and five SNPs in introns 1, 4 (2 SNPs), 8, and 9, respectively. No disease-related mutation was identified by comparing the sequence data of the RSS patients with their clinically normal parents and controls.     

Applying novel genome-wide linkage strategies to search for loci influencing type 2 diabetes and adult height in American Samoa

Type 2 diabetes mellitus (T2DM) is a common complex phenotype that by the year 2010 is predicted to affect 221 million people globally. In the present study we performed a genome-wide linkage scan using the allele-sharing statistic Sall implemented in Allegro and a novel two-dimensional genome-wide strategy implemented in Merloc that searches for pairwise interaction between genetic markers located on different chromosomes linked to T2DM. In addition, we used a robust score statistic from the newly developed QTL-ALL software to search for linkage to variation in adult height. The strategies were applied to a study sample consisting of 238 sib-pairs affected with T2DM from American Samoa. We did not detect any genome-wide significant susceptibility loci for T2DM. However, our two-dimensional linkage investigation detected several loci pairs of interest, including 11q22 and 21q21, 9q21 and 11q22, 1p22-p21 and 4p15, and 4p15 and 15q11-q14, with a two-loci maximum LOD score (MLS) greater than 2.00. Most detected individual loci have previously been identified as susceptibility loci for diabetes-related traits. Our two-dimensional linkage results may facilitate the selection of potential candidate genes and molecular pathways for further diabetes studies because these results, besides providing candidate loci, also demonstrate that polygenic effects may play an important role in T2DM. Linkage was detected (p value of 0.005) for variation in adult height on chromosome 9q31, which was reported previously in other populations. Our finding suggests that the 9q31 region may be a strong quantitative trait locus for adult height, which is likely to be of importance across populations.     

Use of a High-Density Protein Microarray to Identify Autoantibodies in Subjects with Type 2 Diabetes Mellitus and an HLA Background Associated with Reduced Insulin Secretion

New biomarkers for type 2 diabetes mellitus (T2DM) may aid diagnosis, drug development or clinical treatment. Evidence is increasing for the adaptive immune system’s role in T2DM and suggests the presence of unidentified autoantibodies. While high-density protein microarrays have emerged as a useful technology to identify possible novel autoantigens in autoimmune diseases, its application in T2DM has lagged. In Pima Indians, the HLA haplotype (HLA-DRB1*02) is protective against T2DM and, when studied when they have normal glucose tolerance, subjects with this HLA haplotype have higher insulin secretion compared to those without the protective haplotype. Possible autoantibody biomarkers were identified using microarrays containing 9480 proteins in plasma from Pima Indians with T2DM without the protective haplotype (n = 7) compared with those with normal glucose regulation (NGR) with the protective haplotype (n = 11). A subsequent validation phase involving 45 cases and 45 controls, matched by age, sex and specimen storage time, evaluated 77 proteins. Eleven autoantigens had higher antibody signals among T2DM subjects with the lower insulin-secretion HLA background compared with NGR subjects with the higher insulin-secretion HLA background (p<0.05, adjusted for multiple comparisons). PPARG2 and UBE2M had lowest p-values (adjusted p = 0.023) while PPARG2 and RGS17 had highest case-to-control antibody signal ratios (1.7). A multi-protein classifier involving the 11 autoantigens had sensitivity, specificity, and area under the receiver operating characteristics curve of 0.73, 0.80, and 0.83 (95% CI 0.74–0.91, p = 3.4x10^-8), respectively. This study identified 11 novel autoantigens which were associated with T2DM and an HLA background associated with reduced insulin secretion. While further studies are needed to distinguish whether these antibodies are associated with insulin secretion via the HLA background, T2DM more broadly, or a combination of the two, this study may aid the search for autoantibody biomarkers by narrowing the list of protein targets.