Genetic Variants Associated with Serum Thyroid Stimulating Hormone (TSH) Levels in European Americans and African Americans from the eMERGE Network

Thyroid stimulating hormone (TSH) hormone levels are normally tightly regulated within an individual; thus, relatively small variations may indicate thyroid disease. Genome-wide association studies (GWAS) have identified variants in PDE8B and FOXE1 that are associated with TSH levels. However, prior studies lacked racial/ethnic diversity, limiting the generalization of these findings to individuals of non-European ethnicities. The Electronic Medical Records and Genomics (eMERGE) Network is a collaboration across institutions with biobanks linked to electronic medical records (EMRs). The eMERGE Network uses EMR-derived phenotypes to perform GWAS in diverse populations for a variety of phenotypes. In this report, we identified serum TSH levels from 4,501 European American and 351 African American euthyroid individuals in the eMERGE Network with existing GWAS data. Tests of association were performed using linear regression and adjusted for age, sex, body mass index (BMI), and principal components, assuming an additive genetic model. Our results replicate the known association of PDE8B with serum TSH levels in European Americans (rs2046045 p = 1.85×10⁻¹⁷, β = 0.09). FOXE1 variants, associated with hypothyroidism, were not genome-wide significant (rs10759944: p = 1.08×10⁻⁶, β = −0.05). No SNPs reached genome-wide significance in African Americans. However, multiple known associations with TSH levels in European ancestry were nominally significant in African Americans, including PDE8B (rs2046045 p = 0.03, β = −0.09), VEGFA (rs11755845 p = 0.01, β = −0.13), and NFIA (rs334699 p = 1.50×10⁻³, β = −0.17). We found little evidence that SNPs previously associated with other thyroid-related disorders were associated with serum TSH levels in this study. These results support the previously reported association between PDE8B and serum TSH levels in European Americans and emphasize the need for additional genetic studies in more diverse populations.     

Genome-Wide Association of Body Fat Distribution in African Ancestry Populations Suggests New Loci

Central obesity, measured by waist circumference (WC) or waist-hip ratio (WHR), is a marker of body fat distribution. Although obesity disproportionately affects minority populations, few studies have conducted genome-wide association study (GWAS) of fat distribution among those of predominantly African ancestry (AA). We performed GWAS of WC and WHR, adjusted and unadjusted for BMI, in up to 33,591 and 27,350 AA individuals, respectively. We identified loci associated with fat distribution in AA individuals using meta-analyses of GWA results for WC and WHR (stage 1). Overall, 25 SNPs with single genomic control (GC)-corrected p-values<5.0×10⁻⁶ were followed-up (stage 2) in AA with WC and with WHR. Additionally, we interrogated genomic regions of previously identified European ancestry (EA) WHR loci among AA. In joint analysis of association results including both Stage 1 and 2 cohorts, 2 SNPs demonstrated association, rs2075064 at LHX2, p = 2.24×10⁻⁸ for WC-adjusted-for-BMI, and rs6931262 at RREB1, p = 2.48×10⁻⁸ for WHR-adjusted-for-BMI. However, neither signal was genome-wide significant after double GC-correction (LHX2: p = 6.5×10⁻⁸; RREB1: p = 5.7×10⁻⁸). Six of fourteen previously reported loci for waist in EA populations were significant (p<0.05 divided by the number of independent SNPs within the region) in AA studied here (TBX15-WARS2, GRB14, ADAMTS9, LY86, RSPO3, ITPR2-SSPN). Further, we observed associations with metabolic traits: rs13389219 at GRB14 associated with HDL-cholesterol, triglycerides, and fasting insulin, and rs13060013 at ADAMTS9 with HDL-cholesterol and fasting insulin. Finally, we observed nominal evidence for sexual dimorphism, with stronger results in AA women at the GRB14 locus (p for interaction = 0.02). In conclusion, we identified two suggestive loci associated with fat distribution in AA populations in addition to confirming 6 loci previously identified in populations of EA. These findings reinforce the concept that there are fat distribution loci that are independent of generalized adiposity.     

Association of rs780094 in GCKR with Metabolic Traits and Incident Diabetes and Cardiovascular Disease: The ARIC Study

Objective

The minor T-allele of rs780094 in the glucokinase regulator gene (GCKR) associates with a number of metabolic traits including higher triglyceride levels and improved glycemic regulation in study populations of mostly European ancestry. Using data from the Atherosclerosis Risk in Communities (ARIC) Study, we sought to replicate these findings, examine them in a large population-based sample of African American study participants, and to investigate independent associations with other metabolic traits in order to determine if variation in GKCR contributes to their observed clustering. In addition, we examined the association of rs780094 with incident diabetes, coronary heart disease (CHD), and stroke over up mean follow-up times of 8, 15, and 15 years, respectively.

Research Design and Methods

Race-stratified analyses were conducted among 10,929 white and 3,960 black participants aged 45–64 at baseline assuming an additive genetic model and using linear and logistic regression and Cox proportional hazards models.

Results

Previous findings replicated among white participants in multivariable adjusted models: the T-allele of rs780094 was associated with lower fasting glucose (p = 10⁻⁷) and insulin levels (p = 10⁻⁶), lower insulin resistance (HOMA-IR, p = 10⁻⁹), less prevalent diabetes (p = 10⁻⁶), and higher CRP (p = 10⁻⁸), 2-h postprandial glucose (OGTT, p = 10⁻⁶), and triglyceride levels (p = 10⁻³¹). Moreover, the T-allele was independently associated with higher HDL cholesterol levels (p = 0.022), metabolic syndrome prevalence (p = 0.043), and lower beta-cell function measured as HOMA-B (p = 0.011). Among black participants, the T-allele was associated only with higher triglyceride levels (p = 0.004) and lower insulin levels (p = 0.002) and HOMA-IR (p = 0.013). Prospectively, the T-allele was associated with reduced incidence of diabetes (p = 10⁻⁴) among white participants, but not with incidence of CHD or stroke.

Conclusions

Our findings indicate rs780094 has independent associations with multiple metabolic traits as well as incident diabetes, but not incident CHD or stroke. The magnitude of association between the SNP and most traits was of lower magnitude among African American compared to white participants.     

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.     

Epigenome-Wide DNA Methylation in Hearing Ability: New Mechanisms for an Old Problem

Epigenetic regulation of gene expression has been shown to change over time and may be associated with environmental exposures in common complex traits. Age-related hearing impairment is a complex disorder, known to be heritable, with heritability estimates of 57–70%. Epigenetic regulation might explain the observed difference in age of onset and magnitude of hearing impairment with age. Epigenetic epidemiology studies using unrelated samples can be limited in their ability to detect small effects, and recent epigenetic findings in twins underscore the power of this well matched study design. We investigated the association between venous blood DNA methylation epigenome-wide and hearing ability. Pure-tone audiometry (PTA) and Illumina HumanMethylation array data were obtained from female twin volunteers enrolled in the TwinsUK register. Two study groups were explored: first, an epigenome-wide association scan (EWAS) was performed in a discovery sample (n = 115 subjects, age range: 47–83 years, Illumina 27 k array), then replication of the top ten associated probes from the discovery EWAS was attempted in a second unrelated sample (n = 203, age range: 41–86 years, Illumina 450 k array). Finally, a set of monozygotic (MZ) twin pairs (n = 21 pairs) within the discovery sample (Illumina 27 k array) was investigated in more detail in an MZ discordance analysis. Hearing ability was strongly associated with DNA methylation levels in the promoter regions of several genes, including TCF25 (cg01161216, p = 6.6×10⁻⁶), FGFR1 (cg15791248, p = 5.7×10⁻⁵) and POLE (cg18877514, p = 6.3×10⁻⁵). Replication of these results in a second sample confirmed the presence of differential methylation at TCF25 (p(replication) = 6×10⁻⁵) and POLE (p(replication) = 0.016). In the MZ discordance analysis, twins'' intrapair difference in hearing ability correlated with DNA methylation differences at ACP6 (cg01377755, r = −0.75, p = 1.2×10⁻⁴) and MEF2D (cg08156349, r = −0.75, p = 1.4×10⁻⁴). Examination of gene expression in skin, suggests an influence of differential methylation on expression, which may account for the variation in hearing ability with age.     

Genome-Wide Association Study of Lp-PLA2 Activity and Mass in the Framingham Heart Study

Lipoprotein-associated phospholipase A₂ (Lp-PLA₂) is an emerging risk factor and therapeutic target for cardiovascular disease. The activity and mass of this enzyme are heritable traits, but major genetic determinants have not been explored in a systematic, genome-wide fashion. We carried out a genome-wide association study of Lp-PLA₂ activity and mass in 6,668 Caucasian subjects from the population-based Framingham Heart Study. Clinical data and genotypes from the Affymetrix 550K SNP array were obtained from the open-access Framingham SHARe project. Each polymorphism that passed quality control was tested for associations with Lp-PLA₂ activity and mass using linear mixed models implemented in the R statistical package, accounting for familial correlations, and controlling for age, sex, smoking, lipid-lowering-medication use, and cohort. For Lp-PLA₂ activity, polymorphisms at four independent loci reached genome-wide significance, including the APOE/APOC1 region on chromosome 19 (p = 6×10⁻²⁴); CELSR2/PSRC1 on chromosome 1 (p = 3×10⁻¹⁵); SCARB1 on chromosome 12 (p = 1×10⁻⁸) and ZNF259/BUD13 in the APOA5/APOA1 gene region on chromosome 11 (p = 4×10⁻⁸). All of these remained significant after accounting for associations with LDL cholesterol, HDL cholesterol, or triglycerides. For Lp-PLA₂ mass, 12 SNPs achieved genome-wide significance, all clustering in a region on chromosome 6p12.3 near the PLA2G7 gene. Our analyses demonstrate that genetic polymorphisms may contribute to inter-individual variation in Lp-PLA₂ activity and mass.     

Whole-Body Insulin Sensitivity Rather than Body-Mass-Index Determines Fasting and Post-Glucose-Load Growth Hormone Concentrations

Background

Obese, non-acromegalic persons show lower growth hormone (GH) concentrations at fasting and reduced GH nadir during an oral glucose tolerance test (OGTT). However, this finding has never been studied with regard to whole-body insulin-sensitivity as a possible regulator.

Methods

In this retrospective analysis, non-acromegalic (NonACRO, n = 161) and acromegalic (ACRO, n = 35), non-diabetic subjects were subdivided into insulin-sensitive (IS) and –resistant (IR) groups according to the Clamp-like Index (CLIX)-threshold of 5 mg·kg⁻¹·min⁻¹ from the OGTT.

Results

Non-acromegalic IS (CLIX: 8.8±0.4 mg·kg⁻¹·min⁻¹) persons with similar age and sex distribution, but lower (p<0.001) body-mass-index (BMI = 25±0 kg/m², 84% females, 56±1 years) had 59% and 70%, respectively, higher (p<0.03) fasting GH and OGTT GH area under the curve concentrations than IR (CLIX: 3.5±0.1 mg·kg⁻¹·min⁻¹, p<0.001) subjects (BMI = 29±1 kg/m², 73% females, 58±1 years). When comparing on average overweight non-acromegalic IS and IR with similar anthropometry (IS: BMI: 27±0 kg/m², 82% females, 58±2 years; IR: BMI: 27±0 kg/m², 71% females, 60±1 years), but different CLIX (IS: 8.7±0.9 vs. IR: 3.8±0.1 mg·kg⁻¹·min⁻¹, p<0.001), the results remained almost the same. In addition, when adjusted for OGTT-mediated glucose rise, GH fall was less pronounced in IR. In contrast, in acromegalic subjects, no difference was found between IS and IR patients with regard to fasting and post-glucose-load GH concentrations.

Conclusions

Circulating GH concentrations at fasting and during the OGTT are lower in non-acromegalic insulin-resistant subjects. This study seems the first to demonstrate that insulin sensitivity rather than body-mass modulates fasting and post-glucose-load GH concentrations in non-diabetic non–acromegalic subjects.     

SORL1 Is Genetically Associated with Late-Onset Alzheimer’s Disease in Japanese,Koreans and Caucasians

To discover susceptibility genes of late-onset Alzheimer’s disease (LOAD), we conducted a 3-stage genome-wide association study (GWAS) using three populations: Japanese from the Japanese Genetic Consortium for Alzheimer Disease (JGSCAD), Koreans, and Caucasians from the Alzheimer Disease Genetic Consortium (ADGC). In Stage 1, we evaluated data for 5,877,918 genotyped and imputed SNPs in Japanese cases (n = 1,008) and controls (n = 1,016). Genome-wide significance was observed with 12 SNPs in the APOE region. Seven SNPs from other distinct regions with p-values <2×10⁻⁵ were genotyped in a second Japanese sample (885 cases, 985 controls), and evidence of association was confirmed for one SORL1 SNP (rs3781834, P = 7.33×10⁻⁷ in the combined sample). Subsequent analysis combining results for several SORL1 SNPs in the Japanese, Korean (339 cases, 1,129 controls) and Caucasians (11,840 AD cases, 10,931 controls) revealed genome wide significance with rs11218343 (P = 1.77×10⁻⁹) and rs3781834 (P = 1.04×10⁻⁸). SNPs in previously established AD loci in Caucasians showed strong evidence of association in Japanese including rs3851179 near PICALM (P = 1.71×10⁻⁵) and rs744373 near BIN1 (P = 1.39×10⁻⁴). The associated allele for each of these SNPs was the same as in Caucasians. These data demonstrate for the first time genome-wide significance of LOAD with SORL1 and confirm the role of other known loci for LOAD in Japanese. Our study highlights the importance of examining associations in multiple ethnic populations.     

FCGR2B and FCRLB Gene Polymorphisms Associated with IgA Nephropathy

Background

IgA nephropathy (IgAN) is a complex syndrome characterized by deposition of IgA and IgA containing immune complexes (ICs) composed of IgG and complement C3 proteins in the mesangial area of glomeruli. The low-affinity receptors for the Fc region of IgG (FcγRs) are involved in autoantibody/immune complex-induced organ injury as well as ICs clearance. The aim of the study was to associate multiple polymorphisms within FCGR gene locus with IgAN in a large Chinese cohort.

Patients and Methods

60 single nucleotide polymorphisms (SNPs) spanning a 400 kb range within FCGR gene locus were analyzed in 2100 DNA samples from patients with biopsy proven IgAN and healthy age- and sex-matched controls from the same population in Chinese.

Results

Among the 60 SNPs investigated, 15 gene polymorphisms within FCGR gene locus (25%) were associated with susceptibility to IgAN. The most significantly associated SNPs within individual genes were FCGR2B rs12118043 (p = 8.74*10⁻³, OR 0.76, 95% CI 0.62–0.93), and FCRLB rs4657093 (p = 2.28*10⁻³, OR 0.77, 95% CI 0.65–0.91). Both conditional analysis and linkage disequilibrium analysis suggested they were independent signals associated with IgAN. Associations between FCGR2B rs12118043 and proteinuria (p = 3.65×10⁻²) as well as gross hematuria (p = 4.53×10⁻²), between FCRLB rs4657093 and levels of serum creatinine (p = 2.67×10⁻²) as well as eGFR (p = 5.41*10⁻³) were also observed. Electronic cis-expression quantative trait loci analysis supported their possible functional significance, with protective genotypes correlating lower gene expressions.

Conclusion

Our data from genetic associations and expression associations revealed potentially pathogenic roles of Fc receptor gene polymorphisms in IgAN.     

Mechanistic Phenotypes: An Aggregative Phenotyping Strategy to Identify Disease Mechanisms Using GWAS Data

A single mutation can alter cellular and global homeostatic mechanisms and give rise to multiple clinical diseases. We hypothesized that these disease mechanisms could be identified using low minor allele frequency (MAF<0.1) non-synonymous SNPs (nsSNPs) associated with “mechanistic phenotypes”, comprised of collections of related diagnoses. We studied two mechanistic phenotypes: (1) thrombosis, evaluated in a population of 1,655 African Americans; and (2) four groupings of cancer diagnoses, evaluated in 3,009 white European Americans. We tested associations between nsSNPs represented on GWAS platforms and mechanistic phenotypes ascertained from electronic medical records (EMRs), and sought enrichment in functional ontologies across the top-ranked associations. We used a two-step analytic approach whereby nsSNPs were first sorted by the strength of their association with a phenotype. We tested associations using two reverse genetic models and standard additive and recessive models. In the second step, we employed a hypothesis-free ontological enrichment analysis using the sorted nsSNPs to identify functional mechanisms underlying the diagnoses comprising the mechanistic phenotypes. The thrombosis phenotype was solely associated with ontologies related to blood coagulation (Fisher''s p = 0.0001, FDR p = 0.03), driven by the F5, P2RY12 and F2RL2 genes. For the cancer phenotypes, the reverse genetics models were enriched in DNA repair functions (p = 2×10−5, FDR p = 0.03) (POLG/FANCI, SLX4/FANCP, XRCC1, BRCA1, FANCA, CHD1L) while the additive model showed enrichment related to chromatid segregation (p = 4×10−6, FDR p = 0.005) (KIF25, PINX1). We were able to replicate nsSNP associations for POLG/FANCI, BRCA1, FANCA and CHD1L in independent data sets. Mechanism-oriented phenotyping using collections of EMR-derived diagnoses can elucidate fundamental disease mechanisms.     

Insights into the Genetic Architecture of Early Stage Age-Related Macular Degeneration: A Genome-Wide Association Study Meta-Analysis

Genetic factors explain a majority of risk variance for age-related macular degeneration (AMD). While genome-wide association studies (GWAS) for late AMD implicate genes in complement, inflammatory and lipid pathways, the genetic architecture of early AMD has been relatively under studied. We conducted a GWAS meta-analysis of early AMD, including 4,089 individuals with prevalent signs of early AMD (soft drusen and/or retinal pigment epithelial changes) and 20,453 individuals without these signs. For various published late AMD risk loci, we also compared effect sizes between early and late AMD using an additional 484 individuals with prevalent late AMD. GWAS meta-analysis confirmed previously reported association of variants at the complement factor H (CFH) (peak P = 1.5×10⁻³¹) and age-related maculopathy susceptibility 2 (ARMS2) (P = 4.3×10⁻²⁴) loci, and suggested Apolipoprotein E (ApoE) polymorphisms (rs2075650; P = 1.1×10⁻⁶) associated with early AMD. Other possible loci that did not reach GWAS significance included variants in the zinc finger protein gene GLI3 (rs2049622; P = 8.9×10⁻⁶) and upstream of GLI2 (rs6721654; P = 6.5×10⁻⁶), encoding retinal Sonic hedgehog signalling regulators, and in the tyrosinase (TYR) gene (rs621313; P = 3.5×10⁻⁶), involved in melanin biosynthesis. For a range of published, late AMD risk loci, estimated effect sizes were significantly lower for early than late AMD. This study confirms the involvement of multiple established AMD risk variants in early AMD, but suggests weaker genetic effects on the risk of early AMD relative to late AMD. Several biological processes were suggested to be potentially specific for early AMD, including pathways regulating RPE cell melanin content and signalling pathways potentially involved in retinal regeneration, generating hypotheses for further investigation.     

A Genome-Wide Association Study of Circulating Galectin-3

Galectin-3 is a lectin involved in fibrosis, inflammation and proliferation. Increased circulating levels of galectin-3 have been associated with various diseases, including cancer, immunological disorders, and cardiovascular disease. To enhance our knowledge on galectin-3 biology we performed the first genome-wide association study (GWAS) using the Illumina HumanCytoSNP-12 array imputed with the HapMap 2 CEU panel on plasma galectin-3 levels in 3,776 subjects and follow-up genotyping in an additional 3,516 subjects. We identified 2 genome wide significant loci associated with plasma galectin-3 levels. One locus harbours the LGALS3 gene (rs2274273; P = 2.35×10⁻¹⁸⁸) and the other locus the ABO gene (rs644234; P = 3.65×10⁻⁴⁷). The variance explained by the LGALS3 locus was 25.6% and by the ABO locus 3.8% and jointly they explained 29.2%. Rs2274273 lies in high linkage disequilibrium with two non-synonymous SNPs (rs4644; r² = 1.0, and rs4652; r² = 0.91) and wet lab follow-up genotyping revealed that both are strongly associated with galectin-3 levels (rs4644; P = 4.97×10⁻⁴⁶⁵ and rs4652 P = 1.50×10⁻⁴²¹) and were also associated with LGALS3 gene-expression. The origins of our associations should be further validated by means of functional experiments.     

Oligoclonal Band Status in Scandinavian Multiple Sclerosis Patients Is Associated with Specific Genetic Risk Alleles

The presence of oligoclonal bands (OCB) in cerebrospinal fluid (CSF) is a typical finding in multiple sclerosis (MS). We applied data from Norwegian, Swedish and Danish (i.e. Scandinavian) MS patients from a genome-wide association study (GWAS) to search for genetic differences in MS relating to OCB status. GWAS data was compared in 1367 OCB positive and 161 OCB negative Scandinavian MS patients, and nine of the most associated SNPs were genotyped for replication in 3403 Scandinavian MS patients. HLA-DRB1 genotypes were analyzed in a subset of the OCB positive (n = 2781) and OCB negative (n = 292) MS patients and compared to 890 healthy controls. Results from the genome-wide analyses showed that single nucleotide polymorphisms (SNPs) from the HLA complex and six other loci were associated to OCB status. In SNPs selected for replication, combined analyses showed genome-wide significant association for two SNPs in the HLA complex; rs3129871 (p = 5.7×10⁻¹⁵) and rs3817963 (p = 5.7×10⁻¹⁰) correlating with the HLA-DRB1*15 and the HLA-DRB1*04 alleles, respectively. We also found suggestive association to one SNP in the Calsyntenin-2 gene (p = 8.83×10⁻⁷). In HLA-DRB1 analyses HLA-DRB1*15∶01 was a stronger risk factor for OCB positive than OCB negative MS, whereas HLA-DRB1*04∶04 was associated with increased risk of OCB negative MS and reduced risk of OCB positive MS. Protective effects of HLA-DRB1*01∶01 and HLA-DRB1*07∶01 were detected in both groups. The groups were different with regard to age at onset (AAO), MS outcome measures and gender. This study confirms both shared and distinct genetic risk for MS subtypes in the Scandinavian population defined by OCB status and indicates different clinical characteristics between the groups. This suggests differences in disease mechanisms between OCB negative and OCB positive MS with implications for patient management, which need to be further studied.     

Association of Variants at UMOD with Chronic Kidney Disease and Kidney Stones—Role of Age and Comorbid Diseases

Chronic kidney disease (CKD) is a worldwide public health problem that is associated with substantial morbidity and mortality. To search for sequence variants that associate with CKD, we conducted a genome-wide association study (GWAS) that included a total of 3,203 Icelandic cases and 38,782 controls. We observed an association between CKD and a variant with 80% population frequency, rs4293393-T, positioned next to the UMOD gene (GeneID: 7369) on chromosome 16p12 (OR = 1.25, P = 4.1×10⁻¹⁰). This gene encodes uromodulin (Tamm-Horsfall protein), the most abundant protein in mammalian urine. The variant also associates significantly with serum creatinine concentration (SCr) in Icelandic subjects (N = 24,635, P = 1.3×10⁻²³) but not in a smaller set of healthy Dutch controls (N = 1,819, P = 0.39). Our findings validate the association between the UMOD variant and both CKD and SCr recently discovered in a large GWAS. In the Icelandic dataset, we demonstrate that the effect on SCr increases substantially with both age (P = 3.0×10⁻¹⁷) and number of comorbid diseases (P = 0.008). The association with CKD is also stronger in the older age groups. These results suggest that the UMOD variant may influence the adaptation of the kidney to age-related risk factors of kidney disease such as hypertension and diabetes. The variant also associates with serum urea (P = 1.0×10⁻⁶), uric acid (P = 0.0064), and suggestively with gout. In contrast to CKD, the UMOD variant confers protection against kidney stones when studied in 3,617 Icelandic and Dutch kidney stone cases and 43,201 controls (OR = 0.88, P = 5.7×10⁻⁵).     

A Genome-Wide Association Study of Total Serum and Mite-Specific IgEs in Asthma Patients

Immunoglobulin E (IgE) is one of the central players in asthma and allergic diseases. Although the serum IgE level, a useful endophenotype, is generally increased in patients with asthma, genetic factors influencing IgE regulation in asthma are still not fully understood. To identify the genetic variations associated with total serum and mite-specific IgEs in asthmatics, a genome-wide association study (GWAS) of 657,366 single nucleotide polymorphisms (SNPs) was performed in 877 Korean asthmatics. This study found that several new genes might be associated with total IgE in asthmatics, such as CRIM1 (rs848512, P = 1.18×10⁻⁶; rs711254, P = 6.73×10⁻⁶), ZNF71 (rs10404342, P = 7.60×10⁻⁶), TLN1 (rs4879926, P = 7.74×10⁻⁶), and SYNPO2 (rs1472066, P = 8.36×10⁻⁶; rs1038770, P = 8.66×10⁻⁶). Regarding the association of specific IgE to house dust mites, it was observed that intergenic SNPs nearby to OPRK1 and LOC730217 might be associated with Dermatophagoides pteronyssinus (D.p.) and Dermatophagoides farinae (D.f.) in asthmatics, respectively. In further pathway analysis, the phosphatidylinositol signaling system and adherens junction pathways were estimated to play a role in the regulation of total IgE levels in asthma. Although functional evaluations and replications of these results in other populations are needed, this GWAS of serum IgE in asthmatics could facilitate improved understanding of the role of the newly identified genetic variants in asthma and its related phenotypes.     

Genome-Wide Association Study Implicates Chromosome 9q21.31 as a Susceptibility Locus for Asthma in Mexican Children

Many candidate genes have been studied for asthma, but replication has varied. Novel candidate genes have been identified for various complex diseases using genome-wide association studies (GWASs). We conducted a GWAS in 492 Mexican children with asthma, predominantly atopic by skin prick test, and their parents using the Illumina HumanHap 550 K BeadChip to identify novel genetic variation for childhood asthma. The 520,767 autosomal single nucleotide polymorphisms (SNPs) passing quality control were tested for association with childhood asthma using log-linear regression with a log-additive risk model. Eleven of the most significantly associated GWAS SNPs were tested for replication in an independent study of 177 Mexican case–parent trios with childhood-onset asthma and atopy using log-linear analysis. The chromosome 9q21.31 SNP rs2378383 (p = 7.10×10⁻⁶ in the GWAS), located upstream of transducin-like enhancer of split 4 (TLE4), gave a p-value of 0.03 and the same direction and magnitude of association in the replication study (combined p = 6.79×10⁻⁷). Ancestry analysis on chromosome 9q supported an inverse association between the rs2378383 minor allele (G) and childhood asthma. This work identifies chromosome 9q21.31 as a novel susceptibility locus for childhood asthma in Mexicans. Further, analysis of genome-wide expression data in 51 human tissues from the Novartis Research Foundation showed that median GWAS significance levels for SNPs in genes expressed in the lung differed most significantly from genes not expressed in the lung when compared to 50 other tissues, supporting the biological plausibility of our overall GWAS findings and the multigenic etiology of childhood asthma.     

Genome-Wide Association Study Identifies Major Loci for Carcass Weight on BTA14 in Hanwoo (Korean Cattle)

This genome-wide association study (GWAS) was conducted to identify major loci that are significantly associated with carcass weight, and their effects, in order to provide increased understanding of the genetic architecture of carcass weight in Hanwoo. This genome-wide association study identified one major chromosome region ranging from 23 Mb to 25 Mb on chromosome 14 as being associated with carcass weight in Hanwoo. Significant Bonferroni-corrected genome-wide associations (P<1.52×10⁻⁶) were detected for 6 Single Nucleotide Polymorphic (SNP) loci for carcass weight on chromosome 14. The most significant SNP was BTB-01280026 (P = 4.02×10⁻¹¹), located in the 25 Mb region on Bos taurus autosome 14 (BTA14). The other 5 significant SNPs were Hapmap27934-BTC-065223 (P = 4.04×10⁻¹¹) in 25.2 Mb, BTB-01143580 (P = 6.35×10⁻¹¹) in 24.3 Mb, Hapmap30932-BTC-011225 (P = 5.92×10⁻¹⁰) in 24.8 Mb, Hapmap27112-BTC-063342 (P = 5.18×10⁻⁹) in 25.4 Mb, and Hapmap24414-BTC-073009 (P = 7.38×10⁻⁸) in 25.4 Mb, all on BTA 14. One SNP (BTB-01143580; P = 6.35×10⁻¹¹) lies independently from the other 5 SNPs. The 5 SNPs that lie together showed a large Linkage disequilibrium (LD) block (block size of 553 kb) with LD coefficients ranging from 0.53 to 0.89 within the block. The most significant SNPs accounted for 6.73% to 10.55% of additive genetic variance, which is quite a large proportion of the total additive genetic variance. The most significant SNP (BTB-01280026; P = 4.02×10⁻¹¹) had 16.96 kg of allele substitution effect, and the second most significant SNP (Hapmap27934-BTC-065223; P = 4.04×10⁻¹¹) had 18.06 kg of effect on carcass weight, which correspond to 44% and 47%, respectively, of the phenotypic standard deviation for carcass weight in Hanwoo cattle. Our results demonstrated that carcass weight was affected by a major Quantitative Trait Locus (QTL) with a large effect and by many SNPs with small effects that are normally distributed.