Changes in Uric Acid Levels following Bariatric Surgery Are Not Associated with SLC2A9 Variants in the Swedish Obese Subjects Study

Context and Objective

Obesity and SLC2A9 genotype are strong determinants of uric acid levels. However, data on SLC2A9 variants and weight loss induced changes in uric acid levels are missing. We examined whether the changes in uric acid levels two- and ten-years after weight loss induced by bariatric surgery were associated with SLC2A9 single nucleotide polymorphisms (SNPs) in the Swedish Obese Subjects study.

Methods

SNPs (N = 14) identified by genome-wide association studies and exonic SNPs in the SLC2A9 gene locus were genotyped. Cross-sectional associations were tested before (N = 1806), two (N = 1664) and ten years (N = 1201) after bariatric surgery. Changes in uric acid were compared between baseline and Year 2 (N = 1660) and years 2 and 10 (N = 1172). A multiple testing corrected threshold of P = 0.007 was used for statistical significance.

Results

Overall, 11 of the 14 tested SLC2A9 SNPs were significantly associated with cross-sectional uric acid levels at all three time points, with rs13113918 showing the strongest association at each time point (R² = 3.7−5.2%, 3.9×10⁻²²≤p≤7.7×10⁻¹¹). One SNP (rs737267) showed a significant association (R² = 0.60%, P = 0.002) with change in uric acid levels from baseline to Year 2, as common allele homozygotes (C/C, N = 957) showed a larger decrease in uric acid (−61.4 µmol/L) compared to minor allele carriers (A/X: −51.7 µmol/L, N = 702). No SNPs were associated with changes in uric acid from years 2 to 10.

Conclusions

SNPs in the SLC2A9 locus contribute significantly to uric acid levels in obese individuals, and the associations persist even after considerable weight loss due to bariatric surgery. However, we found little evidence for an interaction between genotype and weight change on the response of uric acid to bariatric surgery over ten years. Thus, the fluctuations in uric acid levels among the surgery group appear to be driven by the weight losses and gains, independent of SLC2A9 genotypes.     

Association of Mitochondrial DNA Polymerase γ Gene POLG1 Polymorphisms with Parkinsonism in Chinese Populations

Background

Mitochondrial DNA polymerase gamma (POLG1) mutations were associated with levodopa-responsive Parkinsonism. POLG1 gene contains a number of common nonsynonymous SNPs and intronic regulatory SNPs which may have functional consequences. It is of great interest to discover polymorphisms variants associated with Parkinson''s disease (PD), both in isolation and in combination with specific SNPs.

Materials and Methods

We conducted a case-control study and genotyped twenty SNPs and poly-Q polymorphisms of POLG1 gene including in 344 Chinese sporadic PD patients and 154 healthy controls. All the polymorphisms of POLG1 we found in this study were sequenced by PCR products with dye terminator methods using an ABI-3100 sequencer. Hardy-Weinberg equilibrium and linkage disequilibrium (LD) for association between twenty POLG1 SNPs and PD were calculated using the program Haploview.

Principal Results

We provided evidence for strong association of four intronic SNPs of the POLG1 gene (new report: c.2070-12T>A and rs2307439: c.2070-64G>A in intron 11, P = 0.00011, OR = 1.727; rs2302084: c.3105-11T>C and rs2246900: c.3105-36A>G in intron 19, P = 0.00031, OR = 1.648) with PD. However, we did not identify any significant association between ten exonic SNPs of POLG1 and PD. Linkage disequilibrium analysis indicated that c.2070-12T>A and c.2070-64G>A could be parsed into one block as Haplotype 1 as well as c.3105-11T>C and c.3105-36A>G in Haplotype 2. In addition, case and control study on association of POLG1 CAG repeat (poly-Q) alleles with PD showed a significant association (P = 0.03, OR = 2.16) of the non-10/11Q variants with PD. Although intronic SNPs associated with PD didn''t influence POLG1 mRNA alternative splicing, there was a strong association of c.2070-12T>A and c.2070-64G>A with decreased POLG1 mRNA level and protein levels.

Conclusions

Our findings indicate that POLG1 may play a role in the pathogenesis of PD in Chinese populations.     

Associations of Variants in CHRNA5/A3/B4 Gene Cluster with Smoking Behaviors in a Korean Population

Multiple genome-wide and targeted association studies reveal a significant association of variants in the CHRNA5-CHRNA3-CHRNB4 (CHRNA5/A3/B4) gene cluster on chromosome 15 with nicotine dependence. The subjects examined in most of these studies had a European origin. However, considering the distinct linkage disequilibrium patterns in European and other ethnic populations, it would be of tremendous interest to determine whether such associations could be replicated in populations of other ethnicities, such as Asians. In this study, we performed comprehensive association and interaction analyses for 32 single-nucleotide polymorphisms (SNPs) in CHRNA5/A3/B4 with smoking initiation (SI), smoking quantity (SQ), and smoking cessation (SC) in a Korean sample (N = 8,842). We found nominally significant associations of 7 SNPs with at least one smoking-related phenotype in the total sample (SI: P = 0.015∼0.023; SQ: P = 0.008∼0.028; SC: P = 0.018∼0.047) and the male sample (SI: P = 0.001∼0.023; SQ: P = 0.001∼0.046; SC: P = 0.01). A spectrum of haplotypes formed by three consecutive SNPs located between rs16969948 in CHRNA5 and rs6495316 in the intergenic region downstream from the 5′ end of CHRNB4 was associated with these three smoking-related phenotypes in both the total and the male sample. Notably, associations of these variants and haplotypes with SC appear to be much weaker than those with SI and SQ. In addition, we performed an interaction analysis of SNPs within the cluster using the generalized multifactor dimensionality reduction method and found a significant interaction of SNPs rs7163730 in LOC123688, rs6495308 in CHRNA3, and rs7166158, rs8043123, and rs11072793 in the intergenic region downstream from the 5′ end of CHRNB4 to be influencing SI in the male sample. Considering that fewer than 5% of the female participants were smokers, we did not perform any analysis on female subjects specifically. Together, our detected associations of variants in the CHRNA5/A3/B4 cluster with SI, SQ, and SC in the Korean smoker samples provide strong evidence for the contribution of this cluster to the etiology of SI, ND, and SC in this Asian population.     

Replication of Association between Schizophrenia and Chromosome 6p21-6p22.1 Polymorphisms in Chinese Han Population

Chromosome 6p21-p22.1, spanning the extended major histocompatibility complex (MHC) region, is a highly polymorphic, gene-dense region. It has been identified as a susceptibility locus of schizophrenia in Europeans, Japanese, and Chinese. In our previous two-stage genome-wide association study (GWAS), polymorphisms of zinc finger with KRAB and SCAN domains 4 (ZKSCAN4), nuclear factor-κB-activating protein-like (NKAPL), and piggyBac transposable element derived 1 (PGBD1), localized to chromosome 6p21-p22.1, were strongly associated with schizophrenia. To further investigate the association between polymorphisms at this locus and schizophrenia in the Chinese Han population, we selected eight other single-nucleotide polymorphisms (SNPs) distributed in or near these genes for a case-control association study in an independent sample of 902 cases and 1,091 healthy controls in an attempt to replicate the GWAS results. Four of these eight SNPs (rs12214383, rs1150724, rs3800324, and rs1997660) displayed a nominal difference in allele frequencies between the case and control groups. The association between two of these SNPs and schizophrenia were significant even after Bonferroni correction (rs12000: allele A>G, P = 2.50E-04, odds ratio [OR] = 1.27, 95% confidence interval [CI] = 1.12–1.45; rs1150722: allele C>T, P = 4.28E-05, OR = 0.55, 95% CI = 0.41–0.73). Haplotype ATTGACGC, comprising these eight SNPs (rs2235359, rs2185955, rs12214383, rs12000, rs1150724, rs1150722, rs3800324, and rs1997660), was significantly associated with schizophrenia (P = 6.60E-05). We also performed a combined study of this replication sample and the first-stage GWAS sample. The combined study revealed that rs12000 and rs1150722 were still strongly associated with schizophrenia (rs12000: allele G>A, P _combined  = 0.0019, OR = 0.81; rs1150722: allele G>A, P _combined  = 3.00E-04, OR = 0.61). These results support our findings that locus 6p21-p22.1 is significantly associated with schizophrenia in the Chinese Han population and encourage further studies of the functions of these genetic factors.     

Investigation of Genetic Variants,Birthweight and Hypothalamic-Pituitary-Adrenal Axis Function Suggests a Genetic Variant in the SERPINA6 Gene Is Associated with Corticosteroid Binding Globulin in the Western Australia Pregnancy Cohort (Raine) Study

Background

The hypothalamic-pituitary-adrenal (HPA) axis regulates stress responses and HPA dysfunction has been associated with several chronic diseases. Low birthweight may be associated with HPA dysfunction in later life, yet human studies are inconclusive. The primary study aim was to identify genetic variants associated with HPA axis function. A secondary aim was to evaluate if these variants modify the association between birthweight and HPA axis function in adolescents.

Methods

Morning fasted blood samples were collected from children of the Western Australia Pregnancy Cohort (Raine) at age 17 (n = 1077). Basal HPA axis function was assessed by total cortisol, corticosteroid binding globulin (CBG), and adrenocorticotropic hormone (ACTH). The associations between 124 tag single nucleotide polymorphisms (SNPs) within 16 HPA pathway candidate genes and each hormone were evaluated using multivariate linear regression and penalized linear regression analysis using the HyperLasso method.

Results

The penalized regression analysis revealed one candidate gene SNP, rs11621961 in the CBG encoding gene (SERPINA6), significantly associated with total cortisol and CBG. No other candidate gene SNPs were significant after applying the penalty or adjusting for multiple comparisons; however, several SNPs approached significance. For example, rs907621 (p = 0.002) and rs3846326 (p = 0.003) in the mineralocorticoid receptor gene (NR3C2) were associated with ACTH and SERPINA6 SNPs rs941601 (p = 0.004) and rs11622665 (p = 0.008), were associated with CBG. To further investigate our findings for SERPINA6, rare and common SNPs in the gene were imputed from the 1,000 genomes data and 8 SNPs across the gene were significantly associated with CBG levels after adjustment for multiple comparisons. Birthweight was not associated with any HPA outcome, and none of the gene-birthweight interactions were significant after adjustment for multiple comparisons.

Conclusions

Our study suggests that genetic variation in the SERPINA6 gene may be associated with altered CBG levels during adolescence. Replication of these findings is required.     

Replication and Predictive Value of SNPs Associated with Melanoma and Pigmentation Traits in a Southern European Case-Control Study

Background

Genetic association studies have revealed numerous polymorphisms conferring susceptibility to melanoma. We aimed to replicate previously discovered melanoma-associated single-nucleotide polymorphisms (SNPs) in a Greek case-control population, and examine their predictive value.

Methods

Based on a field synopsis of genetic variants of melanoma (MelGene), we genotyped 284 patients and 284 controls at 34 melanoma-associated SNPs of which 19 derived from GWAS. We tested each one of the 33 SNPs passing quality control for association with melanoma both with and without accounting for the presence of well-established phenotypic risk factors. We compared the risk allele frequencies between the Greek population and the HapMap CEU sample. Finally, we evaluated the predictive ability of the replicated SNPs.

Results

Risk allele frequencies were significantly lower compared to the HapMap CEU for eight SNPs (rs16891982 – SLC45A2, rs12203592 – IRF4, rs258322 – CDK10, rs1805007 – MC1R, rs1805008 - MC1R, rs910873 - PIGU, rs17305573- PIGU, and rs1885120 - MTAP) and higher for one SNP (rs6001027 – PLA2G6) indicating a different profile of genetic susceptibility in the studied population. Previously identified effect estimates modestly correlated with those found in our population (r = 0.72, P<0.0001). The strongest associations were observed for rs401681-T in CLPTM1L (odds ratio [OR] 1.60, 95% CI 1.22–2.10; P = 0.001), rs16891982-C in SCL45A2 (OR 0.51, 95% CI 0.34–0.76; P = 0.001), and rs1805007-T in MC1R (OR 4.38, 95% CI 2.03–9.43; P = 2×10⁻⁵). Nominally statistically significant associations were seen also for another 5 variants (rs258322-T in CDK10, rs1805005-T in MC1R, rs1885120-C in MYH7B, rs2218220-T in MTAP and rs4911442-G in the ASIP region). The addition of all SNPs with nominal significance to a clinical non-genetic model did not substantially improve melanoma risk prediction (AUC for clinical model 83.3% versus 83.9%, p = 0.66).

Conclusion

Overall, our study has validated genetic variants that are likely to contribute to melanoma susceptibility in the Greek population.     

A Large Scale Gene-Centric Association Study of Lung Function in Newly-Hired Female Cotton Textile Workers with Endotoxin Exposure

Background

Occupational exposure to endotoxin is associated with decrements in pulmonary function, but how much variation in this association is explained by genetic variants is not well understood.

Objective

We aimed to identify single nucleotide polymorphisms (SNPs) that are associated with the rate of forced expiratory volume in one second (FEV₁) decline by a large scale genetic association study in newly-hired healthy young female cotton textile workers.

Methods

DNA samples were genotyped using the Illumina Human CVD BeadChip. Change rate in FEV₁ was modeled as a function of each SNP genotype in linear regression model with covariate adjustment. We controlled the type 1 error in study-wide level by permutation method. The false discovery rate (FDR) and the family-wise error rate (FWER) were set to be 0.10 and 0.15 respectively.

Results

Two SNPs were found to be significant (P<6.29×10⁻⁵), including rs1910047 (P = 3.07×10⁻⁵, FDR = 0.0778) and rs9469089 (P = 6.19×10⁻⁵, FDR = 0.0967), as well as other eight suggestive (P<5×10⁻⁴) associated SNPs. Gene-gene and gene-environment interactions were also observed, such as rs1910047 and rs1049970 (P = 0.0418, FDR = 0.0895); rs9469089 and age (P = 0.0161, FDR = 0.0264). Genetic risk score analysis showed that the more risk loci the subjects carried, the larger the rate of FEV₁ decline occurred (P _trend = 3.01×10⁻¹⁸). However, the association was different among age subgroups (P = 7.11×10⁻⁶) and endotoxin subgroups (P = 1.08×10⁻²). Functional network analysis illustrates potential biological connections of all interacted genes.

Conclusions

Genetic variants together with environmental factors interact to affect the rate of FEV₁ decline in cotton textile workers.     

A Genome-Wide Association Study Reveals Variants in ARL15 that Influence Adiponectin Levels

The adipocyte-derived protein adiponectin is highly heritable and inversely associated with risk of type 2 diabetes mellitus (T2D) and coronary heart disease (CHD). We meta-analyzed 3 genome-wide association studies for circulating adiponectin levels (n = 8,531) and sought validation of the lead single nucleotide polymorphisms (SNPs) in 5 additional cohorts (n = 6,202). Five SNPs were genome-wide significant in their relationship with adiponectin (P≤5×10⁻⁸). We then tested whether these 5 SNPs were associated with risk of T2D and CHD using a Bonferroni-corrected threshold of P≤0.011 to declare statistical significance for these disease associations. SNPs at the adiponectin-encoding ADIPOQ locus demonstrated the strongest associations with adiponectin levels (P-combined = 9.2×10⁻¹⁹ for lead SNP, rs266717, n = 14,733). A novel variant in the ARL15 (ADP-ribosylation factor-like 15) gene was associated with lower circulating levels of adiponectin (rs4311394-G, P-combined = 2.9×10⁻⁸, n = 14,733). This same risk allele at ARL15 was also associated with a higher risk of CHD (odds ratio [OR] = 1.12, P = 8.5×10⁻⁶, n = 22,421) more nominally, an increased risk of T2D (OR = 1.11, P = 3.2×10⁻³, n = 10,128), and several metabolic traits. Expression studies in humans indicated that ARL15 is well-expressed in skeletal muscle. These findings identify a novel protein, ARL15, which influences circulating adiponectin levels and may impact upon CHD risk.     

Possible Single-Nucleotide Polymorphism Loci Associated with Systemic Sclerosis Susceptibility: A Genetic Association Study in a Chinese Han Population

Objective

The aim of this study was to confirm the association of RHOB and FAM167A-BLK gene polymorphisms with susceptibility to systemic sclerosis (SSc) in a Chinese Han population.

Methods

A total of 248 SSc patients and 251 healthy controls of Chinese Han ethnicity, which visited the department of dermatology of Peking Union Medical College Hospital, were included in the study. Six selected single nucleotide polymorphisms (SNPs) in the RHOB and FAM167A-BLK regions were selected as markers and were genotyped using a MassARRAY system, which is based on the matrix-assisted laser desorption/ionization time of flight mass spectrometry technique.

Results

Three SNPs in the coding regions of the RHOB and FAM167A-BLK genes displayed an association with SSc: (1) rs1062292T, which is a newly discovered SNP in the RHOB gene (P = 0.03, odds ratio [OR] = 1.62, 95% confidence interval (CI) = 1.05–2.50), (2) rs2736340T (P = 0.03, OR = 1.39, 95%CI = 1.03–1.85), and (3) rs13277113A (P = 0.04, OR = 1.34, 95%CI = 1.01–1.76), both in the FAM167A-BLK gene. Our results support previous findings that vaiants in the RHOB and FAM167A-BLK genes may be associated with susceptibility to SSc. However, the loci of the SNPs in RHOB region that displayed an association with SSc are quite different from the loci which were identified in studies of Caucasian populations.

Conclusion

Our results confirm that RHOB and FAM167A-BLK polymorphisms exist in Chinese Han SSc patients. Therefore, variants of the RHOB and FAM167A-BLK genes are promising genetic markers for SSc.     

Multiple Independent Genetic Factors at NOS1AP Modulate the QT Interval in a Multi-Ethnic Population

Extremes of electrocardiographic QT interval are associated with increased risk for sudden cardiac death (SCD); thus, identification and characterization of genetic variants that modulate QT interval may elucidate the underlying etiology of SCD. Previous studies have revealed an association between a common genetic variant in NOS1AP and QT interval in populations of European ancestry, but this finding has not been extended to other ethnic populations. We sought to characterize the effects of NOS1AP genetic variants on QT interval in the multi-ethnic population-based Dallas Heart Study (DHS, n = 3,072). The SNP most strongly associated with QT interval in previous samples of European ancestry, rs16847548, was the most strongly associated in White (P = 0.005) and Black (P = 3.6×10⁻⁵) participants, with the same direction of effect in Hispanics (P = 0.17), and further showed a significant SNP × sex-interaction (P = 0.03). A second SNP, rs16856785, uncorrelated with rs16847548, was also associated with QT interval in Blacks (P = 0.01), with qualitatively similar results in Whites and Hispanics. In a previously genotyped cohort of 14,107 White individuals drawn from the combined Atherosclerotic Risk in Communities (ARIC) and Cardiovascular Health Study (CHS) cohorts, we validated both the second locus at rs16856785 (P = 7.63×10⁻⁸), as well as the sex-interaction with rs16847548 (P = 8.68×10⁻⁶). These data extend the association of genetic variants in NOS1AP with QT interval to a Black population, with similar trends, though not statistically significant at P<0.05, in Hispanics. In addition, we identify a strong sex-interaction and the presence of a second independent site within NOS1AP associated with the QT interval. These results highlight the consistent and complex role of NOS1AP genetic variants in modulating QT interval.     

Association Analysis of Dyslipidemia-Related Genes in Diabetic Nephropathy

Type 1 diabetes (T1D) increases risk of the development of microvascular complications and cardiovascular disease (CVD). Dyslipidemia is a common risk factor in the pathogenesis of both CVD and diabetic nephropathy (DN), with CVD identified as the primary cause of death in patients with DN. In light of this commonality, we assessed single nucleotide polymorphisms (SNPs) in thirty-seven key genetic loci previously associated with dyslipidemia in a T1D cohort using a case-control design. SNPs (n = 53) were genotyped using Sequenom in 1467 individuals with T1D (718 cases with proteinuric nephropathy and 749 controls without nephropathy i.e. normal albumin excretion). Cases and controls were white and recruited from the UK and Ireland. Association analyses were performed using PLINK to compare allele frequencies in cases and controls. In a sensitivity analysis, samples from control individuals with reduced renal function (estimated glomerular filtration rate<60 ml/min/1.73 m²) were excluded. Correction for multiple testing was performed by permutation testing. A total of 1394 samples passed quality control filters. Following regression analysis adjusted by collection center, gender, duration of diabetes, and average HbA1c, two SNPs were significantly associated with DN. rs4420638 in the APOC1 region (odds ratio [OR]  = 1.51; confidence intervals [CI]: 1.19–1.91; P = 0.001) and rs1532624 in CETP (OR = 0.82; CI: 0.69–0.99; P = 0.034); rs4420638 was also significantly associated in a sensitivity analysis (P = 0.016) together with rs7679 (P = 0.027). However, no association was significant following correction for multiple testing. Subgroup analysis of end-stage renal disease status failed to reveal any association. Our results suggest common variants associated with dyslipidemia are not strongly associated with DN in T1D among white individuals. Our findings, cannot entirely exclude these key genes which are central to the process of dyslipidemia, from involvement in DN pathogenesis as our study had limited power to detect variants of small effect size. Analysis in larger independent cohorts is required.     

miR-485-5p Binding Site SNP rs8752 in HPGD Gene Is Associated with Breast Cancer Risk

Background

Single nucleotide polymorphisms (SNPs) that reside in microRNA target sites may play an important role in breast cancer development and progression. To reveal the association between microRNA target site SNPs and breast cancer risk, we performed a large case-control study in China.

Methods

We performed a two-stage case-control study including 2744 breast cancer cases and 3125 controls. In Stage I, we genotyped 192 SNPs within microRNA binding sites identified from the “Patrocles” database using custom Illumina GoldenGate VeraCode assays on the Illumina BeadXpress platform. In Stage II, genotyping was performed on SNPs potentially associated with breast cancer risk using the TaqMan platform in an independent replication set.

Results

In stage I, 15 SNPs were identified to be significantly associated with breast cancer risk (P<0.05). In stage II, one SNP rs8752 was replicated at P<0.05. This SNP is located in the 3’ untranslated region (UTR) of the 15-hydroxyprostaglandin dehydrogenase (HPGD) gene at 4q34-35, a miR-485-5p binding site. Compared with the GG genotype, the combined GA+AA genotypes has a significantly higher risk of breast cancer (OR = 1.18; 95% CI: 1.06-1.31, P = 0.002). Specifically, this SNP was associated with estrogen receptor (ER) positive breast cancer (P = 0.0007), but not with ER negative breast cancer (P = 0.23), though p for heterogeneity not significant.

Conclusion

Through a systematic case-control study of microRNA binding site SNPs, we identified a new breast cancer risk variant rs8752 in HPGD in Chinese women. Further studies are warranted to investigate the underling mechanism for this association.     

Genome Wide Meta-analysis Highlights the Role of Genetic Variation in RARRES2 in the Regulation of Circulating Serum Chemerin

Chemerin is an adipokine proposed to link obesity and chronic inflammation of adipose tissue. Genetic factors determining chemerin release from adipose tissue are yet unknown. We conducted a meta-analysis of genome-wide association studies (GWAS) for serum chemerin in three independent cohorts from Europe: Sorbs and KORA from Germany and PPP-Botnia from Finland (total N = 2,791). In addition, we measured mRNA expression of genes within the associated loci in peripheral mononuclear cells by micro-arrays, and within adipose tissue by quantitative RT-PCR and performed mRNA expression quantitative trait and expression-chemerin association studies to functionally substantiate our loci. Heritability estimate of circulating chemerin levels was 16.2% in the Sorbs cohort. Thirty single nucleotide polymorphisms (SNPs) at chromosome 7 within the retinoic acid receptor responder 2 (RARRES2)/Leucine Rich Repeat Containing (LRRC61) locus reached genome-wide significance (p<5.0×10⁻⁸) in the meta-analysis (the strongest evidence for association at rs7806429 with p = 7.8×10⁻¹⁴, beta = −0.067, explained variance 2.0%). All other SNPs within the cluster were in linkage disequilibrium with rs7806429 (minimum r² = 0.43 in the Sorbs cohort). The results of the subgroup analyses of males and females were consistent with the results found in the total cohort. No significant SNP-sex interaction was observed. rs7806429 was associated with mRNA expression of RARRES2 in visceral adipose tissue in women (p<0.05 after adjusting for age and body mass index). In conclusion, the present meta-GWAS combined with mRNA expression studies highlights the role of genetic variation in the RARRES2 locus in the regulation of circulating chemerin concentrations.     

Common Single Nucleotide Polymorphisms in Genes Related to Immune Function and Risk of Papillary Thyroid Cancer

Accumulating evidence suggests that alterations in immune function may be important in the etiology of papillary thyroid cancer (PTC). To identify genetic markers in immune-related pathways, we evaluated 3,985 tag single nucleotide polymorphisms (SNPs) in 230 candidate gene regions (adhesion-extravasation-migration, arachidonic acid metabolism/eicosanoid signaling, complement and coagulation cascade, cytokine signaling, innate pathogen detection and antimicrobials, leukocyte signaling, TNF/NF-kB pathway or other) in a case-control study of 344 PTC cases and 452 controls. We used logistic regression models to estimate odds ratios (OR) and calculate one degree of freedom P values of linear trend (P_SNP-trend) for the association between genotype (common homozygous, heterozygous, variant homozygous) and risk of PTC. To correct for multiple comparisons, we applied the false discovery rate method (FDR). Gene region- and pathway-level associations (P_Region and P_Pathway) were assessed by combining individual P_SNP-trend values using the adaptive rank truncated product method. Two SNPs (rs6115, rs6112) in the SERPINA5 gene were significantly associated with risk of PTC (P_SNP-FDR/P_SNP-trend = 0.02/6×10⁻⁶ and P_SNP-FDR/P_SNP-trend = 0.04/2×10⁻⁵, respectively). These associations were independent of a history of autoimmune thyroiditis (OR = 6.4; 95% confidence interval: 3.0–13.4). At the gene region level, SERPINA5 was suggestively associated with risk of PTC (P_Region-FDR/P_Region =0.07/0.0003). Overall, the complement and coagulation cascade pathway was the most significant pathway (P_Pathway = 0.02) associated with PTC risk largely due to the strong effect of SERPINA5. Our results require replication but suggest that the SERPINA5 gene, which codes for the protein C inhibitor involved in many biological processes including inflammation, may be a new susceptibility locus for PTC.     

Endometriosis Is Associated with Rare Copy Number Variants

Endometriosis is a complex gynecological condition that affects 6–10% of women in their reproductive years and is defined by the presence of endometrial glands and stroma outside the uterus. Twin, family, and genome-wide association (GWA) studies have confirmed a genetic role, yet only a small part of the genetic risk can be explained by SNP variation. Copy number variants (CNVs) account for a greater portion of human genetic variation than SNPs and include more recent mutations of large effect. CNVs, likely to be prominent in conditions with decreased reproductive fitness, have not previously been examined as a genetic contributor to endometriosis. Here we employ a high-density genotyping microarray in a genome-wide survey of CNVs in a case-control population that includes 2,126 surgically confirmed endometriosis cases and 17,974 population controls of European ancestry. We apply stringent quality filters to reduce the false positive rate common to many CNV-detection algorithms from 77.7% to 7.3% without noticeable reduction in the true positive rate. We detected no differences in the CNV landscape between cases and controls on the global level which showed an average of 1.92 CNVs per individual with an average size of 142.3 kb. On the local level we identify 22 CNV-regions at the nominal significance threshold (P<0.05), which is greater than the 8.15 CNV-regions expected based on permutation analysis (P<0.001). Three CNV''s passed a genome-wide P-value threshold of 9.3×10⁻⁴; a deletion at SGCZ on 8p22 (P = 7.3×10⁻⁴, OR = 8.5, Cl = 2.3–31.7), a deletion in MALRD1 on 10p12.31 (P = 5.6×10⁻⁴, OR = 14.1, Cl = 2.7–90.9), and a deletion at 11q14.1 (P = 5.7×10⁻⁴, OR = 33.8, Cl = 3.3–1651). Two SNPs within the 22 CNVRs show significant genotypic association with endometriosis after adjusting for multiple testing; rs758316 in DPP6 on 7q36.2 (P = 0.0045) and rs4837864 in ASTN2 on 9q33.1 (P = 0.0002). Together, the CNV-loci are detected in 6.9% of affected women compared to 2.1% in the general population.     

Polymorphisms within the Novel Type 2 Diabetes Risk Locus MTNR1B Determine β-Cell Function

Background

Very recently, a novel type 2 diabetes risk gene, i.e., MTNR1B, was identified and reported to affect fasting glycemia. Using our thoroughly phenotyped cohort of subjects at an increased risk for type 2 diabetes, we assessed the association of common genetic variation within the MTNR1B locus with obesity and prediabetes traits, namely impaired insulin secretion and insulin resistance.

Methodology/Principal Findings

We genotyped 1,578 non-diabetic subjects, metabolically characterized by oral glucose tolerance test, for five tagging single nucleotide polymorphisms (SNPs) covering 100% of common genetic variation (minor allele frequency >0.05) within the MTNR1B locus (rs10830962, rs4753426, rs12804291, rs10830963, rs3781638). In a subgroup (N = 513), insulin sensitivity was assessed by hyperinsulinemic-euglycemic clamp, and in a further subgroup (N = 301), glucose-stimulated insulin secretion was determined by intravenous glucose tolerance test. After appropriate adjustment for confounding variables and Bonferroni correction for multiple comparisons, none of the tagging SNPs was reliably associated with measures of adiposity. SNPs rs10830962, rs4753426, and rs10830963 were significantly associated with higher fasting plasma glucose concentrations (p<0.0001) and reduced OGTT- and IVGTT-induced insulin release (p≤0.0007 and p≤0.01, respectively). By contrast, SNP rs3781638 displayed significant association with lower fasting plasma glucose levels and increased OGTT-induced insulin release (p<0.0001 and p≤0.0002, respectively). Moreover, SNP rs3781638 revealed significant association with elevated fasting- and OGTT-derived insulin sensitivity (p≤0.0021). None of the MTNR1B tagging SNPs altered proinsulin-to-insulin conversion.

Conclusions/Significance

In conclusion, common genetic variation within MTNR1B determines glucose-stimulated insulin secretion and plasma glucose concentrations. Their impact on β-cell function might represent the prevailing pathomechanism how MTNR1B variants increase the type 2 diabetes risk.     

Significant Impact of the MTHFR Polymorphisms and Haplotypes on Male Infertility Risk

Background

Methylenetetrahydrofolate reductase (MTHFR) converts 5,10-methylene tetrahydrofolate to 5-methyl tetrahydrofolate and affects the activity of cellular cycles participating in nucleotide synthesis, DNA repair, genome stability, maintenance of methyl pool, and gene regulation. Genetically compromised MTHFR activity has been suggested to affect male fertility. The objective of the present study was to find the impact on infertility risk of c.203G>A, c.1298A>C, and c.1793G>A polymorphisms in the MTHFR gene.

Methods

PCR-RFLP and DNA sequencing were used to genotype the common SNPs in the MTHFR gene in 630 infertile and 250 fertile males. Chi-square test was applied for statistical comparison of genotype data. Linkage disequilibrium between the SNPs and the frequency of common haplotypes were assessed using Haploview software. Biochemical levels of total homocysteine (tHcy) and folic acid were measured. Meta-analysis on c.1298A>C polymorphism was performed using data from ten studies, comprising 2734 cases and 2737 controls.

Results

c.203G>A and c.1298A>C were found to be unrelated to infertility risk. c.1793G>A was protective against infertility (P = 0.0008). c.677C>T and c.1793G>A were in significant LD (D’ = 0.9). Folic acid and tHcy level did not correlate with male infertility. Pooled estimate on c.1298A>C data from all published studies including our data showed no association of this polymorphism with male infertility (Odds ratio = 1.035, P = 0.56), azoospermia (Odds ratio = 0.97, P = 0.74), or oligoasthenoteratozoospermia (Odds ratio = 0.92, p = 0.29). Eight haplotypes with more than 1% frequency were detected, of which CCGA was protective against infertility (p = 0.02), but the significance of the latter was not seen after applying Bonferroni correction.

Conclusion

Among MTHFR polymorphisms, c.203G>A and c.1298A>C do not affect infertility risk and c.1793G>A is protective against infertility. Haplotype analysis suggested that risk factors on the MTHFR locus do not extend too long on the DNA string.