Bacteria and Genes Involved in Arsenic Speciation in Sediment Impacted by Long-Term Gold Mining

The bacterial community and genes involved in geobiocycling of arsenic (As) from sediment impacted by long-term gold mining were characterized through culture-based analysis of As-transforming bacteria and metagenomic studies of the arsC, arrA, and aioA genes. Sediment was collected from the historically gold mining impacted Mina stream, located in one of the world’s largest mining regions known as the “Iron Quadrangle”. A total of 123 As-resistant bacteria were recovered from the enrichment cultures, which were phenotypically and genotypically characterized for As-transformation. A diverse As-resistant bacteria community was found through phylogenetic analyses of the 16S rRNA gene. Bacterial isolates were affiliated with Proteobacteria, Firmicutes, and Actinobacteria and were represented by 20 genera. Most were AsV-reducing (72%), whereas AsIII-oxidizing accounted for 20%. Bacteria harboring the arsC gene predominated (85%), followed by aioA (20%) and arrA (7%). Additionally, we identified two novel As-transforming genera, Thermomonas and Pannonibacter. Metagenomic analysis of arsC, aioA, and arrA sequences confirmed the presence of these genes, with arrA sequences being more closely related to uncultured organisms. Evolutionary analyses revealed high genetic similarity between some arsC and aioA sequences obtained from isolates and clone libraries, suggesting that those isolates may represent environmentally important bacteria acting in As speciation. In addition, our findings show that the diversity of arrA genes is wider than earlier described, once none arrA-OTUs were affiliated with known reference strains. Therefore, the molecular diversity of arrA genes is far from being fully explored deserving further attention.     

Distribution of Microbial Arsenic Reduction,Oxidation and Extrusion Genes along a Wide Range of Environmental Arsenic Concentrations

The presence of the arsenic oxidation, reduction, and extrusion genes arsC, arrA, aioA, and acr3 was explored in a range of natural environments in northern Chile, with arsenic concentrations spanning six orders of magnitude. A combination of primers from the literature and newly designed primers were used to explore the presence of the arsC gene, coding for the reduction of As (V) to As (III) in one of the most common detoxification mechanisms. Enterobacterial related arsC genes appeared only in the environments with the lowest As concentration, while Firmicutes-like genes were present throughout the range of As concentrations. The arrA gene, involved in anaerobic respiration using As (V) as electron acceptor, was found in all the systems studied. The As (III) oxidation gene aioA and the As (III) transport gene acr3 were tracked with two primer sets each and they were also found to be spread through the As concentration gradient. Sediment samples had a higher number of arsenic related genes than water samples. Considering the results of the bacterial community composition available for these samples, the higher microbial phylogenetic diversity of microbes inhabiting the sediments may explain the increased number of genetic resources found to cope with arsenic. Overall, the environmental distribution of arsenic related genes suggests that the occurrence of different ArsC families provides different degrees of protection against arsenic as previously described in laboratory strains, and that the glutaredoxin (Grx)-linked arsenate reductases related to Enterobacteria do not confer enough arsenic resistance to live above certain levels of As concentrations.     

Microbiological Oxidation of Antimony(III) with Oxygen or Nitrate by Bacteria Isolated from Contaminated Mine Sediments

Bacterial oxidation of arsenite [As(III)] is a well-studied and important biogeochemical pathway that directly influences the mobility and toxicity of arsenic in the environment. In contrast, little is known about microbiological oxidation of the chemically similar anion antimonite [Sb(III)]. In this study, two bacterial strains, designated IDSBO-1 and IDSBO-4, which grow on tartrate compounds and oxidize Sb(III) using either oxygen or nitrate, respectively, as a terminal electron acceptor, were isolated from contaminated mine sediments. Both isolates belonged to the Comamonadaceae family and were 99% similar to previously described species. We identify these novel strains as Hydrogenophaga taeniospiralis strain IDSBO-1 and Variovorax paradoxus strain IDSBO-4. Both strains possess a gene with homology to the aioA gene, which encodes an As(III)-oxidase, and both oxidize As(III) aerobically, but only IDSBO-4 oxidized Sb(III) in the presence of air, while strain IDSBO-1 could achieve this via nitrate respiration. Our results suggest that expression of aioA is not induced by Sb(III) but may be involved in Sb(III) oxidation along with an Sb(III)-specific pathway. Phylogenetic analysis of proteins encoded by the aioA genes revealed a close sequence similarity (90%) among the two isolates and other known As(III)-oxidizing bacteria, particularly Acidovorax sp. strain NO1. Both isolates were capable of chemolithoautotrophic growth using As(III) as a primary electron donor, and strain IDSBO-4 exhibited incorporation of radiolabeled [¹⁴C]bicarbonate while oxidizing Sb(III) from Sb(III)-tartrate, suggesting possible Sb(III)-dependent autotrophy. Enrichment cultures produced the Sb(V) oxide mineral mopungite and lesser amounts of Sb(III)-bearing senarmontite as precipitates.     

Association of Cytokines in Individuals Sensitive and Insensitive to Dust Mites in a Brazilian Population

Introduction

Allergic reaction to dust mites is a relatively common condition among children, triggering cutaneous and respiratory responses that have a great impact on the health of this population. Anaphylactic hypersensitivity is characterized by an exacerbated response involving the production of regulatory cytokines responsible for stimulating the production of IgE antibodies.

Objective

To investigate an association of variants in cytokine genes (IL1A ⁻⁸⁸⁹, IL1B ^{−511, +3962}, IL1R ¹⁹⁷⁰, IL1RA ¹¹¹⁰⁰, IL4RA ⁺¹⁹⁰², IL12 ⁻¹¹⁸⁸, IFNG ⁺⁸⁷⁴, TGFB1 ^{codon 10, codon 25}, TNFA ^{−308, −238}, IL2 ^{−330, +166}, IL4 ^{−1098, −590, −33}, IL6 ^{−174, nt565}, and IL10 ^{−1082, −819, −592}) between patients sensitive to dust mites and a control group.

Methods

A total of 254 patients were grouped as atopic and non-atopic according to sensitivity as evaluated by the Prick Test and to cytokine genotyping by the polymerase chain reaction-sequence specific primers (PCR-SSP) method using the Cytokine Genotyping Kit.

Results

A comparison between individuals allergic to Dermatophagoides farinae, Dermatophagoides pteronyssinus, and Blomia tropicalis and a non-atopic control group showed significant differences between allele and genotype frequencies in the regulatory regions of cytokine genes, with important evidence for IL4 ⁻⁵⁹⁰ in T/C (10.2% vs. 43.1%, odd ratio [OR] = 0.15, p = 5.2 10⁻⁸, pc = 0.0000011, and 95% confidence interval [95%CI] = 0.07–0.32) and T/T genotypes (42.9% vs. 13.8%, OR = 4.69, p = 2.5 10⁻⁶, pc = 0.000055, and 95%CI = 2.42–9.09). Other associations were observed in the pro-inflammatory cytokines IL1A ⁻⁸⁸⁹ (T/T, C, and T) and IL2 ⁻³³⁰ (G/T and T/T) and the anti-inflammatory cytokines IL4RA ⁺¹⁹⁰² (A and G), IL4 ⁻⁵⁹⁰ (T/C, T/T, C, and T), and IL10 ⁻⁵⁹² (A/A, C/A, A, and C).

Conclusion

Our results suggest a possible association between single nucleotide polymorphisms (SNPs) in cytokine genes and hypersensitivity to dust mites.     

A Divergent Artiodactyl MYADM-like Repeat Is Associated with Erythrocyte Traits and Weight of Lamb Weaned in Domestic Sheep

A genome-wide association study (GWAS) was performed to investigate seven red blood cell (RBC) phenotypes in over 500 domestic sheep (Ovis aries) from three breeds (Columbia, Polypay, and Rambouillet). A single nucleotide polymorphism (SNP) showed genome-wide significant association with increased mean corpuscular hemoglobin concentration (MCHC, P = 6.2×10⁻¹⁴) and genome-wide suggestive association with decreased mean corpuscular volume (MCV, P = 2.5×10⁻⁶). The ovine HapMap project found the same genomic region and the same peak SNP has been under extreme historical selective pressure, demonstrating the importance of this region for survival, reproduction, and/or artificially selected traits. We observed a large (>50 kb) variant haplotype sequence containing a full-length divergent artiodactyl MYADM-like repeat in strong linkage disequilibrium with the associated SNP. MYADM gene family members play roles in membrane organization and formation in myeloid cells. However, to our knowledge, no member of the MYADM gene family has been identified in development of morphologically variant RBCs. The specific RBC differences may be indicative of alterations in morphology. Additionally, erythrocytes with altered morphological structure often exhibit increased structural fragility, leading to increased RBC turnover and energy expenditure. The divergent artiodactyl MYADM-like repeat was also associated with increased ewe lifetime kilograms of lamb weaned (P = 2×10⁻⁴). This suggests selection for normal RBCs might increase lamb weights, although further validation is required before implementation in marker-assisted selection. These results provide clues to explain the strong selection on the artiodactyl MYADM-like repeat locus in sheep, and suggest MYADM family members may be important for RBC morphology in other mammals.     

Genetic Architecture of Skin and Eye Color in an African-European Admixed Population

Variation in human skin and eye color is substantial and especially apparent in admixed populations, yet the underlying genetic architecture is poorly understood because most genome-wide studies are based on individuals of European ancestry. We study pigmentary variation in 699 individuals from Cape Verde, where extensive West African/European admixture has given rise to a broad range in trait values and genomic ancestry proportions. We develop and apply a new approach for measuring eye color, and identify two major loci (HERC2[OCA2] P = 2.3×10⁻⁶², SLC24A5 P = 9.6×10⁻⁹) that account for both blue versus brown eye color and varying intensities of brown eye color. We identify four major loci (SLC24A5 P = 5.4×10⁻²⁷, TYR P = 1.1×10⁻⁹, APBA2[OCA2] P = 1.5×10⁻⁸, SLC45A2 P = 6×10⁻⁹) for skin color that together account for 35% of the total variance, but the genetic component with the largest effect (∼44%) is average genomic ancestry. Our results suggest that adjacent cis-acting regulatory loci for OCA2 explain the relationship between skin and eye color, and point to an underlying genetic architecture in which several genes of moderate effect act together with many genes of small effect to explain ∼70% of the estimated heritability.     

Modeling of Environmental Effects in Genome-Wide Association Studies Identifies SLC2A2 and HP as Novel Loci Influencing Serum Cholesterol Levels

Genome-wide association studies (GWAS) have identified 38 larger genetic regions affecting classical blood lipid levels without adjusting for important environmental influences. We modeled diet and physical activity in a GWAS in order to identify novel loci affecting total cholesterol, LDL cholesterol, HDL cholesterol, and triglyceride levels. The Swedish (SE) EUROSPAN cohort (N _SE = 656) was screened for candidate genes and the non-Swedish (NS) EUROSPAN cohorts (N _NS = 3,282) were used for replication. In total, 3 SNPs were associated in the Swedish sample and were replicated in the non-Swedish cohorts. While SNP rs1532624 was a replication of the previously published association between CETP and HDL cholesterol, the other two were novel findings. For the latter SNPs, the p-value for association was substantially improved by inclusion of environmental covariates: SNP rs5400 (p _{SE,unadjusted} = 3.6×10⁻⁵, p _SE,adjusted = 2.2×10⁻⁶, p _{NS,unadjusted} = 0.047) in the SLC2A2 (Glucose transporter type 2) and rs2000999 (p _{SE,unadjusted} = 1.1×10⁻³, p _SE,adjusted = 3.8×10⁻⁴, p _{NS,unadjusted} = 0.035) in the HP gene (Haptoglobin-related protein precursor). Both showed evidence of association with total cholesterol. These results demonstrate that inclusion of important environmental factors in the analysis model can reveal new genetic susceptibility loci.     

The Associations between Immunity-Related Genes and Breast Cancer Prognosis in Korean Women

We investigated the role of common genetic variation in immune-related genes on breast cancer disease-free survival (DFS) in Korean women. 107 breast cancer patients of the Seoul Breast Cancer Study (SEBCS) were selected for this study. A total of 2,432 tag single nucleotide polymorphisms (SNPs) in 283 immune-related genes were genotyped with the GoldenGate Oligonucleotide pool assay (OPA). A multivariate Cox-proportional hazard model and polygenic risk score model were used to estimate the effects of SNPs on breast cancer prognosis. Harrell’s C index was calculated to estimate the predictive accuracy of polygenic risk score model. Subsequently, an extended gene set enrichment analysis (GSEA-SNP) was conducted to approximate the biological pathway. In addition, to confirm our results with current evidence, previous studies were systematically reviewed. Sixty-two SNPs were statistically significant at p-value less than 0.05. The most significant SNPs were rs1952438 in SOCS4 gene (hazard ratio (HR) = 11.99, 95% CI = 3.62–39.72, P = 4.84E-05), rs2289278 in TSLP gene (HR = 4.25, 95% CI = 2.10–8.62, P = 5.99E-05) and rs2074724 in HGF gene (HR = 4.63, 95% CI = 2.18–9.87, P = 7.04E-05). In the polygenic risk score model, the HR of women in the 3^rd tertile was 6.78 (95% CI = 1.48–31.06) compared to patients in the 1^st tertile of polygenic risk score. Harrell’s C index was 0.813 with total patients and 0.924 in 4-fold cross validation. In the pathway analysis, 18 pathways were significantly associated with breast cancer prognosis (P<0.1). The IL-6R, IL-8, IL-10RB, IL-12A, and IL-12B was associated with the prognosis of cancer in data of both our study and a previous study. Therefore, our results suggest that genetic polymorphisms in immune-related genes have relevance to breast cancer prognosis among Korean women.     

Sequence Analysis of Six Blood Pressure Candidate Regions in 4,178 Individuals: The Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Targeted Sequencing Study

Background

Genome-wide association studies (GWAS) identified multiple loci for blood pressure (BP) and hypertension. Six genes – ATP2B1, CACNB2, CYP17A1, JAG1, PLEKHA7, and SH2B3 – were evaluated for sequence variation with large effects on systolic blood pressure (SBP), diastolic blood pressure (DBP), pulse pressure (PP), and mean arterial pressure (MAP).

Methods and Results

Targeted genomic sequence was determined in 4,178 European ancestry participants from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium. Common variants (≥50 minor allele copies) were evaluated individually and rare variants (minor allele frequency, MAF≤1%) were aggregated by locus. 464 common variants were identified across the 6 genes. An upstream CYP17A1 variant, rs11191416 (MAF = 0.09), was the most significant association for SBP (P = 0.0005); however the association was attenuated (P = 0.0469) after conditioning on the GWAS index single nucleotide polymorphism (SNP). A PLEKHA7 intronic variant was the strongest DBP association (rs12806040, MAF = 0.007, P = 0.0006) and was not in LD (r² = 0.01) with the GWAS SNP. A CACNB2 intronic SNP, rs1571787, was the most significant association with PP (MAF = 0.27, P = 0.0003), but was not independent from the GWAS SNP (r² = 0.34). Three variants (rs6163 and rs743572 in the CYP17A1 region and rs112467382 in PLEKHA7) were associated with BP traits (P<0.001). Rare variation, aggregately assessed in the 6 regions, was not significantly associated with BP measures.

Conclusion

Six targeted gene regions, previously identified by GWAS, did not harbor novel variation with large effects on BP in this sample.     

Genetic Association Study of TNFAIP3, IFIH1, IRF5 Polymorphisms with Polymyositis/Dermatomyositis in Chinese Han Population

Background

Single-nucleotide polymorphisms (SNPs) in the TNFAIP3, IFIH1, and IRF5 genes have been associated with several auto-inflammation diseases, while the susceptibility between these genes and idiopathic inflammatory myopathies (IIMs) were not reported. This study aimed to investigate whether TNFAIP3, IFIH1, and IRF5 gene polymorphisms confer susceptibility for the IIMs in Chinese Han population.

Methods

A large case–control study of Chinese subjects with polymyositis (PM) (n = 298) and dermatomyositis (DM) (n = 530) was accomplished. 968 healthy and ethnically matched controls were available for comparison. Six SNPs in the TNFAIP3 region (rs2230926 and rs5029939), the IFIH1 gene (rs1990760 and rs3747517) and the IRF5 region (rs4728142 and rs729302) were assessed and genotyped using the Sequenom MassArray iPLEX platform.

Results

Our study indicated a strong allele association was observed in PM/DM and PM patients for rs2230926 (OR: 1.61, 95%CI: 1.20–2.16, P_c = 7.5×10⁻³; OR: 1.88, 95%CI: 1.30–2.74, P_c = 4.0×10⁻³, respectively) and rs5029939 (OR: 1.64, 95%CI: 1.21–2.21, P_c = 6.0×10⁻³; OR: 1.88, 95%CI: 1.28–2.76, P_c = 5.5×10⁻³,respectively). And rs2230926 and rs5029939 were significantly associated with interstitial lung disease (ILD) in PM/DM and PM patients (P_c = 0.04 and P_c = 0.016; P_c = 0.02 and P_c = 0.03, respectively). In addition, rs4728142 allele and genotype had significant association with PM/DM patients (P_c = 0.026 and P_c = 0.048, respectively). Further analysis with three logistic regression genetic models revealed statistically significant difference in the genotypic distribution in the PM/DM, PM or DM patients when the additive and dominant models were used.

Conclusions

This was the first study to reveal TNFAIP3 and IRF5 polymorphisms were associated with PM/DM patients or these patients with ILD, indicating that TNFAIP3 and IRF5 might be the susceptibility gene for PM/DM patients in Chinese Han population.     

Quantitative Assessment of the Association between rs2046210 at 6q25.1 and Breast Cancer Risk

Genome-wide association studies (GWAS) have identified several genetic susceptibility loci for breast cancer (BC). One of them, conducted among Chinese women, found an association of rs2046210 at 6q25.1 with the risk of BC recently. Since then, numerous association studies have been carried out to investigate the relationship between this polymorphism and BC risk in various populations. However, these have yielded contradictory results. We therefore performed a meta-analysis to clarify this inconsistency. Overall, a total of 235003 subjects based on 13 studies were included in our study. Significantly increased BC risk was detected in the pooled analysis [allele contrast: OR = 1.13, 95%CI = 1.10–1.17, P(Z) <10⁻⁵, P(Q) <10⁻⁴; dominant model: OR = 1.21, 95%CI = 1.14–1.27, P(Z) <10⁻⁵, P(Q) <10⁻⁴; recessive model: OR = 1.18, 95%CI = 1.12–1.24, P(Z) <10⁻⁵, P(Q) = 0.04]. In addition, our data revealed that rs2046210 conferred greater risk in estrogen receptor (ER)-negative tumors [OR = 1.27, 95%CI = 1.15–1.40, P(Z) <10⁻⁵, P(Q) <10⁻⁴] than in ER-positive ones [OR = 1.18, 95%CI = 1.09–1.28, P(Z) <10⁻⁴, P(Q) = 0.0003]. When stratified by ethnicity, significant associations were found in Caucasian and Asian populations, but not detected among Africans. There was evidence of heterogeneity (P<0.05), however, the heterogeneity largely disappeared after stratification by ethnicity. The present meta-analysis demonstrated that the rs2046210 polymorphism may be associated with increased BC susceptibility, but this association varies in different ethnicities.     

Different Genetic Associations of the IgE Production among Fetus,Infancy and Childhood

Elevation of serum IgE levels has long been associated with allergic diseases. Many genes have been linked to IgE production, but few have been linked to the developmental aspects of genetic association with IgE production. To clarify developmental genetic association, we investigated what genes and gene-gene interactions affect IgE levels among fetus, infancy and childhood in Taiwan individuals. A birth cohort of 571 children with completion of IgE measurements from newborn to 1.5, 3, and 6 years of age was subject to genetic association analysis on the 384-customized SNPs of 159 allergy candidate genes. Fifty-three SNPs in 37 genes on innate and adaptive immunity, and stress and response were associated with IgE production. Polymorphisms of the IL13, and the HLA-DPA1 and HLA-DQA1 were, respectively, the most significantly associated with the IgE production at newborn and 6 years of age. Analyses of gene-gene interactions indentified that the combination of NPSR1, rs324981 TT with FGF1, rs2282797 CC had the highest risk (85.7%) of IgE elevation at 1.5 years of age (P = 1.46×10⁻⁴). The combination of IL13, CYFIP2 and PDE2A was significantly associated with IgE elevation at 3 years of age (P = 5.98×10⁻⁷), and the combination of CLEC2D, COLEC11 and CCL2 was significantly associated with IgE elevation at 6 years of age (P = 6.65×10⁻⁷). Our study showed that the genetic association profiles of the IgE production among fetus, infancy and childhood are different. Genetic markers for early prediction and prevention of allergic sensitization may rely on age-based genetic association profiles.     

Assessing the Microbial Community and Functional Genes in a Vertical Soil Profile with Long-Term Arsenic Contamination

Arsenic (As) contamination in soil and groundwater has become a serious problem to public health. To examine how microbial communities and functional genes respond to long-term arsenic contamination in vertical soil profile, soil samples were collected from the surface to the depth of 4 m (with an interval of 1 m) after 16-year arsenic downward infiltration. Integrating BioLog and functional gene microarray (GeoChip 3.0) technologies, we showed that microbial metabolic potential and diversity substantially decreased, and community structure was markedly distinct along the depth. Variations in microbial community functional genes, including genes responsible for As resistance, carbon and nitrogen cycling, phosphorus utilization and cytochrome c oxidases were detected. In particular, changes in community structures and activities were correlated with the biogeochemical features along the vertical soil profile when using the rbcL and nifH genes as biomarkers, evident for a gradual transition from aerobic to anaerobic lifestyles. The C/N showed marginally significant correlations with arsenic resistance (p = 0.069) and carbon cycling genes (p = 0.073), and significant correlation with nitrogen fixation genes (p = 0.024). The combination of C/N, NO₃ ⁻ and P showed the highest correlation (r = 0.779, p = 0.062) with the microbial community structure. Contradict to our hypotheses, a long-term arsenic downward infiltration was not the primary factor, while the spatial isolation and nutrient availability were the key forces in shaping the community structure. This study provides new insights about the heterogeneity of microbial community metabolic potential and future biodiversity preservation for arsenic bioremediation management.     

Investigation of Type 2 Diabetes Risk Alleles Support CDKN2A/B,CDKAL1, and TCF7L2 As Susceptibility Genes in a Han Chinese Cohort

Background

Recent genome-wide association studies (GWASs) have reported several genetic variants to be reproducibly associated with type 2 diabetes. Additional variants have also been detected from a metaanalysis of three GWASs, performed in populations of European ancestry. In the present study, we evaluated the influence of 17 genetic variants from 15 candidate loci, identified in type 2 diabetes GWASs and the metaanalysis, in a Han Chinese cohort.

Methodology/Principal Findings

Selected type 2 diabetes–associated genetic variants were genotyped in 1,165 type 2 diabetic patients and 1,136 normoglycemic control individuals of Southern Han Chinese ancestry. The OR for risk of developing type 2 diabetes was calculated using a logistic regression model adjusted for age, sex, and BMI. Genotype-phenotype associations were tested using a multivariate linear regression model. Genetic variants in CDKN2A/B, CDKAL1, TCF7L2, TCF2, MC4R, and PPARG showed a nominal association with type 2 diabetes (P≤0.05), of whom the three first would stand correction for multiple testing: CDKN2A/B rs10811661, OR: 1.26 (1.12–1.43) P = 1.8*10⁻⁴; CDKAL1 rs10946398, OR: 1.23 (1.09–1.39); P = 7.1*10⁻⁴, and TCF7L2 rs7903146, OR: 1.61 (1.19–2.18) P = 2.3 * 10⁻³. Only nominal phenotype associations were observed, notably for rs8050136 in FTO and fasting plasma glucose (P = 0.002), postprandial plasma glucose (P = 0.002), and fasting C-peptide levels (P = 0.006) in the diabetic patients, and with BMI in controls (P = 0.033).

Conclusions/Significance

We have identified significant association between variants in CDKN2A/B, CDKAL1 and TCF7L2, and type 2 diabetes in a Han Chinese cohort, indicating these genes as strong candidates conferring susceptibility to type 2 diabetes across different ethnicities.     

Polymorphisms of Serotonin Receptor 2A and 2C Genes and COMT in Relation to Obesity and Type 2 Diabetes

Background

Candidate genes of psychological importance include 5HT2A, 5HT2C, and COMT, implicated in the serotonin, noradrenaline and dopamine pathways, which also may be involved in regulation of energy balance. We investigated the associations of single nucleotide polymorphisms (SNPs) of these genes with obesity and metabolic traits.

Methodology/Principal Findings

In a population of 166 200 young men examined at the draft boards, obese men (n = 726, BMI≥31.0 kg/m²) and a randomly selected group (n = 831) were re-examined at two surveys at mean ages 46 and 49 years (S-46, S-49). Anthropometric, physiological and biochemical measures were available. Logistic regression analyses were used to assess age-adjusted odds ratios. No significant associations were observed of 5HT2A rs6311, 5HT2C rs3813929 and COMT rs4680 with obesity, except that COMT rs4680 GG-genotype was associated with fat-BMI (OR = 1.08, CI = 1.01–1.16). The SNPs were associated with a number of physiological variables; most importantly 5HT2C rs3813929 T-allele was associated with glucose (OR = 4.56, CI = 1.13–18.4) and acute insulin response (OR = 0.65, CI = 0.44–0.94) in S-49. COMT rs4680 GG-genotype was associated with glucose (OR = 1.04, CI = 1.00–1.09). Except for an association between 5HT2A rs6311 and total-cholesterol at both surveys, significant in S-46 (OR = 2.66, CI = 1.11–6.40), no significant associations were observed for the other phenotypes. Significant associations were obtained when combined genotype of 5HT2C rs3813929 and COMT rs4680 were examined in relation to BMI (OR = 1.12, CI = 1.03–1.21), fat-BMI (OR = 1.22, CI = 1.08–1.38), waist (OR = 1.13, CI = 1.04–1.22), and cholesterol (OR = 5.60, CI = 0.99–31.4). Analyses of impaired glucose tolerance (IGT) and type 2 diabetes (T2D) revealed, a 12.3% increased frequency of 5HT2C rs3813929 T-allele and an 11.6% increased frequency of COMT rs4680 GG-genotype in individuals with IGT or T2D (χ², p = 0.05 and p = 0.06, respectively). Examination of the combined genotypes of 5HT2C and COMT showed a 34.0% increased frequency of IGT or T2D (χ², p = 0.01).

Conclusions

The findings lend further support to the involvement of serotonin, noradrenaline and dopamine pathways on obesity and glucose homeostasis, in particular when combined genotype associations are explored.     

Obesity-Related Genomic Loci Are Associated with Type 2 Diabetes in a Han Chinese Population

Background and Aims

Obesity is a well-known risk factor for type 2 diabetes. Genome-wide association studies have identified a number of genetic loci associated with obesity. The aim of this study is to examine the contribution of obesity-related genomic loci to type 2 diabetes in a Chinese population.

Methods

We successfully genotyped 18 obesity-related single nucleotide polymorphisms among 5338 type 2 diabetic patients and 4663 controls. Both individual and joint effects of these single nucleotide polymorphisms on type 2 diabetes and quantitative glycemic traits (assessing β-cell function and insulin resistance) were analyzed using logistic and linear regression models, respectively.

Results

Two single nucleotide polymorphisms near MC4R and GNPDA2 genes were significantly associated with type 2 diabetes before adjusting for body mass index and waist circumference (OR (95% CI) = 1.14 (1.06, 1.22) for the A allele of rs12970134, P = 4.75×10⁻⁴; OR (95% CI) = 1.10 (1.03, 1.17) for the G allele of rs10938397, P = 4.54×10⁻³). When body mass index and waist circumference were further adjusted, the association of MC4R with type 2 diabetes remained significant (P = 1.81×10⁻²) and that of GNPDA2 was attenuated (P = 1.26×10⁻¹), suggesting the effect of the locus including GNPDA2 on type 2 diabetes may be mediated through obesity. Single nucleotide polymorphism rs2260000 within BAT2 was significantly associated with type 2 diabetes after adjusting for body mass index and waist circumference (P = 1.04×10⁻²). In addition, four single nucleotide polymorphisms (near or within SEC16B, BDNF, MAF and PRL genes) showed significant associations with quantitative glycemic traits in controls even after adjusting for body mass index and waist circumference (all P values<0.05).

Conclusions

This study indicates that obesity-related genomic loci were associated with type 2 diabetes and glycemic traits in the Han Chinese population.     

Insight in Genome-Wide Association of Metabolite Quantitative Traits by Exome Sequence Analyses

Metabolite quantitative traits carry great promise for epidemiological studies, and their genetic background has been addressed using Genome-Wide Association Studies (GWAS). Thus far, the role of less common variants has not been exhaustively studied. Here, we set out a GWAS for metabolite quantitative traits in serum, followed by exome sequence analysis to zoom in on putative causal variants in the associated genes. ¹H Nuclear Magnetic Resonance (¹H-NMR) spectroscopy experiments yielded successful quantification of 42 unique metabolites in 2,482 individuals from The Erasmus Rucphen Family (ERF) study. Heritability of metabolites were estimated by SOLAR. GWAS was performed by linear mixed models, using HapMap imputations. Based on physical vicinity and pathway analyses, candidate genes were screened for coding region variation using exome sequence data. Heritability estimates for metabolites ranged between 10% and 52%. GWAS replicated three known loci in the metabolome wide significance: CPS1 with glycine (P-value  = 1.27×10⁻³²), PRODH with proline (P-value  = 1.11×10⁻¹⁹), SLC16A9 with carnitine level (P-value  = 4.81×10⁻¹⁴) and uncovered a novel association between DMGDH and dimethyl-glycine (P-value  = 1.65×10⁻¹⁹) level. In addition, we found three novel, suggestively significant loci: TNP1 with pyruvate (P-value  = 1.26×10⁻⁸), KCNJ16 with 3-hydroxybutyrate (P-value  = 1.65×10⁻⁸) and 2p12 locus with valine (P-value  = 3.49×10⁻⁸). Exome sequence analysis identified potentially causal coding and regulatory variants located in the genes CPS1, KCNJ2 and PRODH, and revealed allelic heterogeneity for CPS1 and PRODH. Combined GWAS and exome analyses of metabolites detected by high-resolution ¹H-NMR is a robust approach to uncover metabolite quantitative trait loci (mQTL), and the likely causative variants in these loci. It is anticipated that insight in the genetics of intermediate phenotypes will provide additional insight into the genetics of complex traits.     

Functional Variants in NFKBIE and RTKN2 Involved in Activation of the NF-κB Pathway Are Associated with Rheumatoid Arthritis in Japanese

Rheumatoid arthritis is an autoimmune disease with a complex etiology, leading to inflammation of synovial tissue and joint destruction. Through a genome-wide association study (GWAS) and two replication studies in the Japanese population (7,907 cases and 35,362 controls), we identified two gene loci associated with rheumatoid arthritis susceptibility (NFKBIE at 6p21.1, rs2233434, odds ratio (OR) = 1.20, P = 1.3×10⁻¹⁵; RTKN2 at 10q21.2, rs3125734, OR = 1.20, P = 4.6×10⁻⁹). In addition to two functional non-synonymous SNPs in NFKBIE, we identified candidate causal SNPs with regulatory potential in NFKBIE and RTKN2 gene regions by integrating in silico analysis using public genome databases and subsequent in vitro analysis. Both of these genes are known to regulate the NF-κB pathway, and the risk alleles of the genes were implicated in the enhancement of NF-κB activity in our analyses. These results suggest that the NF-κB pathway plays a role in pathogenesis and would be a rational target for treatment of rheumatoid arthritis.