Analysis of TNFAIP3, a feedback inhibitor of nuclear factor-κB and the neighbor intergenic 6q23 region in rheumatoid arthritis susceptibility

Introduction

Genome-wide association studies of rheumatoid arthritis (RA) have identified an association of the disease with a 6q23 region devoid of genes. TNFAIP3, an RA candidate gene, flanks this region, and polymorphisms in both the TNFAIP3 gene and the intergenic region are associated with systemic lupus erythematosus. We hypothesized that there is a similar association with RA, including polymorphisms in TNFAIP3 and the intergenic region.

Methods

To test this hypothesis, we selected tag-single nucleotide polymorphisms (SNPs) in both loci. They were analyzed in 1,651 patients with RA and 1,619 control individuals of Spanish ancestry.

Results

Weak evidence of association was found both in the 6q23 intergenic region and in the TNFAIP3 locus. The rs582757 SNP and a common haplotype in the TNFAIP3 locus exhibited association with RA. In the intergenic region, two SNPs were associated, namely rs609438 and rs13207033. The latter was only associated in patients with anti-citrullinated peptide antibodies. Overall, statistical association was best explained by the interdependent contribution of SNPs from the two loci TNFAIP3 and the 6q23 intergenic region.

Conclusions

Our data are consistent with the hypothesis that several RA genetic factors exist in the 6q23 region, including polymorphisms in the TNFAIP3 gene, like that previously described for systemic lupus erythematosus.     

Association and interaction of the IL4R,IL4, and IL13 loci with type 1 diabetes among Filipinos

In the search for genes involved in type 1 diabetes (T1D), other than the well-established risk alleles at the human leukocyte antigen loci, we have investigated the association and interaction of polymorphisms in genes involved in the IL4/IL13 pathway in a sample of 90 Filipino patients with T1D and 94 controls. Ten single-nucleotide polymorphisms (SNPs), including two promoter SNPs in the IL4R locus on chromosome 16p11, one promoter SNP in the IL4 locus on chromosome 5q31, and four SNPs--including two promoter SNPs--in the IL13 locus on chromosome 5q31 were examined for association, linkage disequilibrium, and interaction. We found that both individual SNPs (IL4R L389L; odds ratio [OR] 0.34; 95% confidence interval [CI] 0.17-0.67; P=.001) and specific haplotypes both in IL4R (OR 0.10; 95% CI 0-0.5; P=.001) and for the five linked IL4 and IL13 SNPs (OR 3.47; P=.004) were strongly associated with susceptibility to T1D. Since IL4 and IL13 both serve as ligands for a receptor composed, in part, of the IL4R alpha chain, we looked for potential epistasis between polymorphisms in the IL4R locus on chromosome 16p11 and the five SNPs in the IL4 and IL13 loci on chromosome 5q31 and found, through use of a logistic-regression model, significant gene-gene interactions (P=.045, corrected for multiple comparisons by permutation analysis). Our data suggest that the risk for T1D is determined, in part, by polymorphisms within the IL4R locus, including promoter and coding-sequence variants, and by specific combinations of genotypes at the IL4R and the IL4 and IL13 loci.     

Identification of Allelic Heterogeneity at Type-2 Diabetes Loci and Impact on Prediction

Although over 60 single nucleotide polymorphisms (SNPs) have been identified by meta-analysis of genome-wide association studies for type-2 diabetes (T2D) among individuals of European descent, much of the genetic variation remains unexplained. There are likely many more SNPs that contribute to variation in T2D risk, some of which may lie in the regions surrounding established SNPs - a phenomenon often referred to as allelic heterogeneity. Here, we use the summary statistics from the DIAGRAM consortium meta-analysis of T2D genome-wide association studies along with linkage disequilibrium patterns inferred from a large reference sample to identify novel SNPs associated with T2D surrounding each of the previously established risk loci. We then examine the extent to which the use of these additional SNPs improves prediction of T2D risk in an independent validation dataset. Our results suggest that multiple SNPs at each of 3 loci contribute to T2D susceptibility (TCF7L2, CDKN2A/B, and KCNQ1; p<5×10⁻⁸). Using a less stringent threshold (p<5×10⁻⁴), we identify 34 additional loci with multiple associated SNPs. The addition of these SNPs slightly improves T2D prediction compared to the use of only the respective lead SNPs, when assessed using an independent validation cohort. Our findings suggest that some currently established T2D risk loci likely harbor multiple polymorphisms which contribute independently and collectively to T2D risk. This opens a promising avenue for improving prediction of T2D, and for a better understanding of the genetic architecture of T2D.     

Family-Based Association Study of Pulmonary Function in a Population in Northeast Asia

The spirometric measurement of pulmonary function by measuring the forced expiratory volume in one second (FEV₁) is a heritable trait that reflects the physiological condition of the lung and airways. Genome-wide linkage and association studies have identified a number of genes and genetic loci associated with pulmonary function. However, limited numbers of studies have been reported for Asian populations. In this study, we aimed to investigate genetic evidence of pulmonary function in a population in northeast Asia. We conducted a family-based association test with 706 GENDISCAN study participants from 72 Mongolian families to determine candidate genetic determinants of pulmonary function. For the replication, we chose seven candidate single nucleotide polymorphisms (SNPs) from the 5 loci, and tested 1062 SNPs for association with FEV₁ from 2,729 subjects of the Korea Healthy Twin study. We identified TMEM132C as a potential candidate gene at 12q24.3, which is a previously reported locus of asthma and spirometric indices. We also found two adjacent candidate genes (UNC93A and TTLL2) in the 6q27 region, which has been previously identified as a pulmonary function locus in the Framingham cohort study. Our findings suggest that novel candidate genes (TMEM132C, UNC93A and TTLL2) in two different regions are associated with pulmonary function in a population in northeast Asia.     

Immunochip Analysis Identifies Multiple Susceptibility Loci for Systemic Sclerosis

In this study, 1,833 systemic sclerosis (SSc) cases and 3,466 controls were genotyped with the Immunochip array. Classical alleles, amino acid residues, and SNPs across the human leukocyte antigen (HLA) region were imputed and tested. These analyses resulted in a model composed of six polymorphic amino acid positions and seven SNPs that explained the observed significant associations in the region. In addition, a replication step comprising 4,017 SSc cases and 5,935 controls was carried out for several selected non-HLA variants, reaching a total of 5,850 cases and 9,401 controls of European ancestry. Following this strategy, we identified and validated three SSc risk loci, including DNASE1L3 at 3p14, the SCHIP1-IL12A locus at 3q25, and ATG5 at 6q21, as well as a suggested association of the TREH-DDX6 locus at 11q23. The associations of several previously reported SSc risk loci were validated and further refined, and the observed peak of association in PXK was related to DNASE1L3. Our study has increased the number of known genetic associations with SSc, provided further insight into the pleiotropic effects of shared autoimmune risk factors, and highlighted the power of dense mapping for detecting previously overlooked susceptibility loci.     

Overlapping genetic susceptibility variants between three autoimmune disorders: rheumatoid arthritis, type 1 diabetes and coeliac disease

Introduction

Genome wide association studies, replicated by numerous well powered validation studies, have revealed a large number of loci likely to play a role in susceptibility to many multifactorial diseases. It is now well established that some of these loci are shared between diseases with similar aetiology. For example, a number of autoimmune diseases have been associated with variants in the PTPN22, TNFAIP3 and CTLA4 genes. Here we have attempted to define overlapping genetic variants between rheumatoid arthritis (RA), type 1 diabetes (T1D) and coeliac disease (CeD).

Methods

We selected eight SNPs previously identified as being associated with CeD and six T1D-associated SNPs for validation in a sample of 3,962 RA patients and 3,531 controls. Genotyping was performed using the Sequenom MassArray platform and comparison of genotype and allele frequencies between cases and controls was undertaken. A trend test P-value < 0.004 was regarded as significant.

Results

We found statistically significant evidence for association of the TAGAP locus with RA (P = 5.0 × 10^-4). A marker at one other locus, C1QTNF6, previously associated with T1D, showed nominal association with RA in the current study but did not remain statistically significant at the corrected threshold.

Conclusions

In exploring the overlap between T1D, CeD and RA, there is strong evidence that variation within the TAGAP gene is associated with all three autoimmune diseases. Interestingly a number of loci appear to be specific to one of the three diseases currently studied suggesting that they may play a role in determining the particular autoimmune phenotype at presentation.     

Novel Loci for Non-Syndromic Coarctation of the Aorta in Sporadic and Familial Cases

BackroundCoarctation of the aorta (CoA) accounts for 5-8% of all congenital heart defects. CoA can be detected in up to 20% of patients with Ullrich-Turner syndrome (UTS), in which a part or all of one of the X chromosomes is absent. The etiology of non-syndromic CoA is poorly understood. In the present work, we test the hypothesis that rare copy number variation (CNV) especially on the gonosomes, contribute to the etiology of non-syndromic CoA.MethodsWe performed high-resolution genome-wide CNV analysis using the Affymetrix SNP 6.0 microarray platform for 70 individuals with sporadic CoA, 3 families with inherited CoA (n=13) and 605 controls. Our analysis comprised genome wide association, CNV burden and linkage. CNV was validated by multiplex ligation-dependent probe amplification.ResultsWe identified a significant abundance of large (>100 kb) CNVs on the X chromosome in males with CoA (p=0.005). 11 out of 51 (~ 22%) male cases had these large CNVs. Association analysis in the sporadic cohort revealed 14 novel loci for CoA. The locus on 21q22.3 in the sporadic CoA cohort overlapped with a gene locus identified in all familial cases of CoA (candidate gene TRPM2). We identified one CNV locus within a locus with high multipoint LOD score from a linkage analysis of the familial cases (SEPT9); another locus overlapped with a region implicated in Kabuki syndrome. In the familial cases, we identified a total of 7 CNV loci that were exclusively present in cases but not in unaffected family members.ConclusionOf all candidate loci identified, the TRPM2 locus was the most frequently implicated autosomal locus in sporadic and familial cases. However, the abundance of large CNVs on the X chromosome of affected males suggests that gonosomal aberrations are not only responsible for syndromic CoA but also involved in the development of sporadic and non-syndromic CoA and their male dominance.     

Effects of GWAS-Associated Genetic Variants on lncRNAs within IBD and T1D Candidate Loci

Long non-coding RNAs are a new class of non-coding RNAs that are at the crosshairs in many human diseases such as cancers, cardiovascular disorders, inflammatory and autoimmune disease like Inflammatory Bowel Disease (IBD) and Type 1 Diabetes (T1D). Nearly 90% of the phenotype-associated single-nucleotide polymorphisms (SNPs) identified by genome-wide association studies (GWAS) lie outside of the protein coding regions, and map to the non-coding intervals. However, the relationship between phenotype-associated loci and the non-coding regions including the long non-coding RNAs (lncRNAs) is poorly understood. Here, we systemically identified all annotated IBD and T1D loci-associated lncRNAs, and mapped nominally significant GWAS/ImmunoChip SNPs for IBD and T1D within these lncRNAs. Additionally, we identified tissue-specific cis-eQTLs, and strong linkage disequilibrium (LD) signals associated with these SNPs. We explored sequence and structure based attributes of these lncRNAs, and also predicted the structural effects of mapped SNPs within them. We also identified lncRNAs in IBD and T1D that are under recent positive selection. Our analysis identified putative lncRNA secondary structure-disruptive SNPs within and in close proximity (+/−5 kb flanking regions) of IBD and T1D loci-associated candidate genes, suggesting that these RNA conformation-altering polymorphisms might be associated with diseased-phenotype. Disruption of lncRNA secondary structure due to presence of GWAS SNPs provides valuable information that could be potentially useful for future structure-function studies on lncRNAs.     

Identification of New Genetic Risk Variants for Type 2 Diabetes

Although more than 20 genetic susceptibility loci have been reported for type 2 diabetes (T2D), most reported variants have small to moderate effects and account for only a small proportion of the heritability of T2D, suggesting that the majority of inter-person genetic variation in this disease remains to be determined. We conducted a multistage, genome-wide association study (GWAS) within the Asian Consortium of Diabetes to search for T2D susceptibility markers. From 590,887 SNPs genotyped in 1,019 T2D cases and 1,710 controls selected from Chinese women in Shanghai, we selected the top 2,100 SNPs that were not in linkage disequilibrium (r²<0.2) with known T2D loci for in silico replication in three T2D GWAS conducted among European Americans, Koreans, and Singapore Chinese. The 5 most promising SNPs were genotyped in an independent set of 1,645 cases and 1,649 controls from Shanghai, and 4 of them were further genotyped in 1,487 cases and 3,316 controls from 2 additional Chinese studies. Consistent associations across all studies were found for rs1359790 (13q31.1), rs10906115 (10p13), and rs1436955 (15q22.2) with P-values (per allele OR, 95%CI) of 6.49×10⁻⁹ (1.15, 1.10–1.20), 1.45×10⁻⁸ (1.13, 1.08–1.18), and 7.14×10⁻⁷ (1.13, 1.08–1.19), respectively, in combined analyses of 9,794 cases and 14,615 controls. Our study provides strong evidence for a novel T2D susceptibility locus at 13q31.1 and the presence of new independent risk variants near regions (10p13 and 15q22.2) reported by previous GWAS.     

Molecular genetics of type 1 diabetes mellitus: Achievements and future trends

Type 1 diabetes mellitus (T1DM) is a widespread severe disease that results from autoimmune destruction of β cells in Langerhans islets of the pancreas. To date, several loci involved in T1DM have been reliably identified using various approaches: the MHC locus, VNTR within the 5′-nontranscibed region of the insulin gene (INS), CTLA4 (T-cell surface receptor), PTPN22, PTPN2 (T-cell tyrosine phosphatases), IL2 (interleukin 2, IL-2), IL2RA (IL-2 receptor α chain), KIAA0350 (unknown function), and IFIH1 (receptor for double-stranded DNA generated in virus infections). Functional analysis of their protein products confirmed the hypothesis that T1DM is underlain by deregulation of the mechanisms of immune tolerance and, on the other hand, a destructive immune response against the body’s own proteins after virus infection or some other immune stress. Thus the protein products of MHC, INS, PTPN22, and PTPN2 are involved in the intrathymic formation of the T-cell repertoire, responsible for immune defense of the body. On the other hand, nonspecific activation of T cells, which starts autoimmune destruction of pancreatic β cells, is most likely associated with the protein products of CTLA4, IL2, IL2RA, and, possibly, PTPN22 and PTPN2. Apart from the genes with unknown functions, the only exception is IFIH1, but its association with T1DM confirms that certain virus infections can activate autoreactive T cells and lead to T1DM.     

Genetic variants influencing lipid levels and risk of dyslipidemia in Chinese population

Recently, several human genetic and genomewide association studies (GWAS) have discovered many genetic loci that are associated with the concentration of the blood lipids. To confirm the reported loci in Chinese population, we conducted a cross-section study to analyse the association of 25 reported SNPs, genotyped by the ABI SNaPshot method, with the blood levels of total cholesterol (TC), low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C) and triglycerides (TG) in 1900 individuals by multivariate analysis. Logistic regression was applied to assess the association of the genetic loci with the risk of different types of dyslipidemia. Our study has convincingly identified that 12 of 25 studied SNPs were strongly associated with one or more blood lipid parameters (TC, LDL, HDL and TG). Among the 12 associated SNPs, 10 significantly influence the risk of one or more types of dyslipidemia. We firstly found four SNPs (rs12654264 in HMGCR; rs2479409 in PCSK9; rs16996148 in CILP2, PBX4; rs4420638 in APOE-C1-C4-C2) robustly and independently associate with four types of dyslipidemia (MHL, mixed hyperlipidemia; IHTC, isolated hypercholesterolemia; ILH, isolated low HDL-C; IHTG, isolated hypertriglyceridemia). Our results suggest that genetic susceptibility is different on the same candidate locus for the different populations. Meanwhile, most of the reported genetic variants strongly influence one or more plasma lipid levels and the risk of dyslipidemia in Chinese population.     

Association of kidney structure-related gene variants with type 2 diabetes-attributed end-stage kidney disease in African Americans

African Americans (AAs) are at higher risk for developing end-stage kidney disease (ESKD) compared to European Americans. Genome-wide association studies have identified variants associated with diabetic and non-diabetic kidney diseases. Nephropathy loci, including SLC7A9, UMOD, and SHROOM3, have been implicated in the maintenance of normal glomerular and renal tubular structure and function. Herein, 47 genes important in podocyte, glomerular basement membrane, mesangial cell, mesangial matrix, renal tubular cell, and renal interstitium structure were examined for association with type 2 diabetes (T2D)-attributed ESKD in AAs. Single-variant association analysis was performed in the discovery stage, including 2041 T2D-ESKD cases and 1140 controls (non-diabetic, non-nephropathy). Discrimination analyses in 667 T2D cases-lacking nephropathy excluded T2D-associated SNPs. Nominal associations were tested in an additional 483 T2D-ESKD cases and 554 controls in the replication stage. Meta-analysis of 4218 discovery and replication samples revealed three significant associations with T2D-ESKD at CD2AP and MMP2 (P _corr < 0.05 corrected for effective number of SNPs in each locus). Removal of APOL1 renal-risk genotype carriers revealed additional association at five loci, TTC21B, COL4A3, NPHP3-ACAD11, CLDN8, and ARHGAP24 (P _corr < 0.05). Genetic variants at COL4A3, CLDN8, and ARHGAP24 were potentially pathogenic. Gene-based associations revealed suggestive significant aggregate effects of coding variants at four genes. Our findings suggest that genetic variation in kidney structure-related genes may contribute to T2D-attributed ESKD in the AA population.     

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.     

Meta-Analysis of Genome-Wide Association Studies in African Americans Provides Insights into the Genetic Architecture of Type 2 Diabetes

Type 2 diabetes (T2D) is more prevalent in African Americans than in Europeans. However, little is known about the genetic risk in African Americans despite the recent identification of more than 70 T2D loci primarily by genome-wide association studies (GWAS) in individuals of European ancestry. In order to investigate the genetic architecture of T2D in African Americans, the MEta-analysis of type 2 DIabetes in African Americans (MEDIA) Consortium examined 17 GWAS on T2D comprising 8,284 cases and 15,543 controls in African Americans in stage 1 analysis. Single nucleotide polymorphisms (SNPs) association analysis was conducted in each study under the additive model after adjustment for age, sex, study site, and principal components. Meta-analysis of approximately 2.6 million genotyped and imputed SNPs in all studies was conducted using an inverse variance-weighted fixed effect model. Replications were performed to follow up 21 loci in up to 6,061 cases and 5,483 controls in African Americans, and 8,130 cases and 38,987 controls of European ancestry. We identified three known loci (TCF7L2, HMGA2 and KCNQ1) and two novel loci (HLA-B and INS-IGF2) at genome-wide significance (4.15×10⁻⁹⁴<P<5×10⁻⁸, odds ratio (OR) = 1.09 to 1.36). Fine-mapping revealed that 88 of 158 previously identified T2D or glucose homeostasis loci demonstrated nominal to highly significant association (2.2×10⁻²³ < locus-wide P<0.05). These novel and previously identified loci yielded a sibling relative risk of 1.19, explaining 17.5% of the phenotypic variance of T2D on the liability scale in African Americans. Overall, this study identified two novel susceptibility loci for T2D in African Americans. A substantial number of previously reported loci are transferable to African Americans after accounting for linkage disequilibrium, enabling fine mapping of causal variants in trans-ethnic meta-analysis studies.     

Genome wide CNV analysis reveals additional variants associated with milk production traits in Holsteins

Background

Milk production is an economically important sector of global agriculture. Much attention has been paid to the identification of quantitative trait loci (QTL) associated with milk, fat, and protein yield and the genetic and molecular mechanisms underlying them. Copy number variation (CNV) is an emerging class of variants which may be associated with complex traits.

Results

In this study, we performed a genome-wide association between CNVs and milk production traits in 26,362 Holstein bulls and cows. A total of 99 candidate CNVs were identified using Illumina BovineSNP50 array data, and association tests for each production trait were performed using a linear regression analysis with PCA correlation. A total of 34 CNVs on 22 chromosomes were significantly associated with at least one milk production trait after false discovery rate (FDR) correction. Some of those CNVs were located within or near known QTL for milk production traits. We further investigated the relationship between associated CNVs with neighboring SNPs. For all 82 combinations of traits and CNVs (less than 400 kb in length), we found 17 cases where CNVs directly overlapped with tag SNPs and 40 cases where CNVs were adjacent to tag SNPs. In 5 cases, CNVs located were in strong linkage disequilibrium with tag SNPs, either within or adjacent to the same haplotype block. There were an additional 20 cases where CNVs did not have a significant association with SNPs, suggesting that the effects of those CNVs were probably not captured by tag SNPs.

Conclusion

We conclude that combining CNV with SNP analyses reveals more genetic variations underlying milk production traits than those revealed by SNPs alone.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-683) contains supplementary material, which is available to authorized users.     

Family-Based Association Analysis Confirms the Role of the Chromosome 9q21.32 Locus in the Susceptibility of Diabetic Nephropathy

A genome-wide association scan of type 1 diabetic patients from the GoKinD collections previously identified four novel diabetic nephropathy susceptibility loci that have subsequently been shown to be associated with diabetic nephropathy in unrelated patients with type 2 diabetes. To expand these findings, we examined whether single nucleotide polymorphisms (SNPs) at these susceptibility loci were associated with diabetic nephropathy in patients from the Joslin Study of Genetics of Nephropathy in Type 2 Diabetes Family Collection. Six SNPs across the four loci identified in the GoKinD collections and 7 haplotype tagging SNPs, were genotyped in 66 extended families of European ancestry. Pedigrees from this collection contained an average of 18.5 members, including 2 to 14 members with type 2 diabetes. Among diabetic family members, the 9q21.32 locus approached statistical significance with advanced diabetic nephropathy (P = 0.037 [adjusted P = 0.222]). When we expanded our definition of diabetic nephropathy to include individuals with high microalbuminuria, the strength of this association improved significantly (P = 1.42×10⁻³ [adjusted P = 0.009]). This same locus also trended toward statistical significance with variation in urinary albumin excretion in family members with type 2 diabetes (P = 0.032 [adjusted P = 0.192]) and in analyses expanded to include all relatives (P = 0.019 [adjusted P = 0.114]). These data increase support that SNPs identified in the GoKinD collections on chromosome 9q21.32 are true diabetic nephropathy susceptibility loci.     

A study of CNVs as trait-associated polymorphisms and as expression quantitative trait loci

We conducted a comprehensive study of copy number variants (CNVs) well-tagged by SNPs (r(2)≥ 0.8) by analyzing their effect on gene expression and their association with disease susceptibility and other complex human traits. We tested whether these CNVs were more likely to be functional than frequency-matched SNPs as trait-associated loci or as expression quantitative trait loci (eQTLs) influencing phenotype by altering gene regulation. Our study found that CNV-tagging SNPs are significantly enriched for cis eQTLs; furthermore, we observed that trait associations from the NHGRI catalog show an overrepresentation of SNPs tagging CNVs relative to frequency-matched SNPs. We found that these SNPs tagging CNVs are more likely to affect multiple expression traits than frequency-matched variants. Given these findings on the functional relevance of CNVs, we created an online resource of expression-associated CNVs (eCNVs) using the most comprehensive population-based map of CNVs to inform future studies of complex traits. Although previous studies of common CNVs that can be typed on existing platforms and/or interrogated by SNPs in genome-wide association studies concluded that such CNVs appear unlikely to have a major role in the genetic basis of several complex diseases examined, our findings indicate that it would be premature to dismiss the possibility that even common CNVs may contribute to complex phenotypes and at least some common diseases.     

Variation in CHI3LI in Relation to Type 2 Diabetes and Related Quantitative Traits

Background

CHI3LI encoding the inflammatory glycoprotein YKL-40 is located on chromosome 1q32.1. YKL-40 is involved in inflammatory processes and patients with Type 2 Diabetes (T2D) have elevated circulating YKL-40 levels which correlate with their level of insulin resistance. Interestingly, it has been reported that rs10399931 (−329 G/A) of CHI3LI contributes to the inter-individual plasma YKL-40 levels in patients with sarcoidosis, and that rs4950928 (−131 C/G) is a susceptibility polymorphism for asthma and a decline in lung function. We hypothesized that single nucleotide polymorphisms (SNPs) or haplotypes thereof the CHI3LI locus might influence risk of T2D. The aim of the present study was to investigate the putative association between SNPs and haplotype blocks of CHI3LI and T2D and T2D related quantitative traits.

Methods/Principal Findings

Eleven SNPs of CHI3LI were genotyped in 6514 individuals from the Inter99 cohort and 2924 individuals from the outpatient clinic at Steno Diabetes Center. In cas-control studies a total of 2345 T2D patients and 5302 individuals with a normal glucose tolerance test were examined.We found no association between rs10399931 (OR, 0.98 (CI, 0.88–1.10), p = 0.76), rs4950928 (0.98 (0.87–1.10), p = 0.68) or any of the other SNPs with T2D. Similarly, we found no significant association between any of the 11 tgSNPs and T2D related quantitative traits, all p>0.14. None of the identified haplotype blocks of CHI3LI showed any association with T2D, all p>0.16.

Conclusions/Significance

None of the examined SNPs or haplotype blocks of CHI3LI showed any association with T2D or T2D related quantitative traits. Estimates of insulin resistance and dysregulated glucose homeostasis in T2D do not seem to be accounted for by the examined variations of CHI3LI.     

A comprehensive resequence analysis of the KLK15–KLK3–KLK2 locus on chromosome 19q13.33

Single nucleotide polymorphisms (SNPs) in the KLK3 gene on chromosome 19q13.33 are associated with serum prostate-specific antigen (PSA) levels. Recent genome wide association studies of prostate cancer have yielded conflicting results for association of the same SNPs with prostate cancer risk. Since the KLK3 gene encodes the PSA protein that forms the basis for a widely used screening test for prostate cancer, it is critical to fully characterize genetic variation in this region and assess its relationship with the risk of prostate cancer. We have conducted a next-generation sequence analysis in 78 individuals of European ancestry to characterize common (minor allele frequency, MAF >1%) genetic variation in a 56 kb region on chromosome 19q13.33 centered on the KLK3 gene (chr19:56,019,829–56,076,043 bps). We identified 555 polymorphic loci in the process including 116 novel SNPs and 182 novel insertion/deletion polymorphisms (indels). Based on tagging analysis, 144 loci are necessary to tag the region at an r ² threshold of 0.8 and MAF of 1% or higher, while 86 loci are required to tag the region at an r ² threshold of 0.8 and MAF >5%. Our sequence data augments coverage by 35 and 78% as compared to variants in dbSNP and HapMap, respectively. We observed six non-synonymous amino acid or frame shift changes in the KLK3 gene and three changes in each of the neighboring genes, KLK15 and KLK2. Our study has generated a detailed map of common genetic variation in the genomic region surrounding the KLK3 gene, which should be useful for fine-mapping the association signal as well as determining the contribution of this locus to prostate cancer risk and/or regulation of PSA expression.