Genome-wide association study in Han Chinese identifies three novel loci for human height

Human height is a complex genetic trait with high heritability but discovery efforts in Asian populations are limited. We carried out a meta-analysis of genome-wide association studies (GWAS) for height in 6,534 subjects with in silico replication of 1,881 subjects in Han Chinese. We identified three novel loci reaching the genome-wide significance threshold (P < 5 × 10^?8), which mapped in or near ZNF638 (rs12612930, P = 2.02 × 10^?10), MAML2 (rs11021504, P = 7.81 × 10^?9), and C18orf12 (rs11082671, P = 1.87 × 10^?8). We also confirmed two loci previously reported in European populations including CS (rs3816804, P = 2.63 × 10^?9) and CYP19A1 (rs3751599, P = 4.80 × 10^?10). In addition, we provided evidence supporting 35 SNPs identified by previous GWAS (P < 0.05). Our study provides new insights into the genetic determination of biological regulation of human height.     

Genome-wide scan revealed that polymorphisms in the PNPLA3, SAMM50, and PARVB genes are associated with development and progression of nonalcoholic fatty liver disease in Japan

We examined the genetic background of nonalcoholic fatty liver disease (NAFLD) in the Japanese population, by performing a genome-wide association study (GWAS). For GWAS, 392 Japanese NAFLD subjects and 934 control individuals were analyzed. For replication studies, 172 NAFLD and 1,012 control subjects were monitored. After quality control, 261,540 single-nucleotide polymorphisms (SNPs) in autosomal chromosomes were analyzed using a trend test. Association analysis was also performed using multiple logistic regression analysis using genotypes, age, gender and body mass index (BMI) as independent variables. Multiple linear regression analyses were performed to evaluate allelic effect of significant SNPs on biochemical traits and histological parameters adjusted by age, gender, and BMI. Rs738409 in the PNPLA3 gene was most strongly associated with NAFLD after adjustment (P = 6.8 × 10^?14, OR = 2.05). Rs2896019, and rs381062 in the PNPLA3 gene, rs738491, rs3761472, and rs2143571 in the SAMM50 gene, rs6006473, rs5764455, and rs6006611 in the PARVB gene had also significant P values (<2.0 × 10^?10) and high odds ratios (1.84–2.02). These SNPs were found to be in the same linkage disequilibrium block and were associated with decreased serum triglycerides and increased aspartate aminotransferase (AST) and alanine aminotransferase (ALT) in NAFLD patients. These SNPs were associated with steatosis grade and NAFLD activity score (NAS). Rs738409, rs2896019, rs738491, rs6006473, rs5764455, and rs6006611 were associated with fibrosis. Polymorphisms in the SAMM50 and PARVB genes in addition to those in the PNPLA3 gene were observed to be associated with the development and progression of NAFLD.     

Genetic susceptibility to obesity and diet intakes: association and interaction analyses in the Malmö Diet and Cancer Study

Gene–environment interactions need to be studied to better understand the obesity. We aimed at determining whether genetic susceptibility to obesity associates with diet intake levels and whether diet intakes modify the genetic susceptibility. In 29,480 subjects of the population-based Malmö Diet and Cancer Study (MDCS), we first assessed association between 16 genome-wide association studies identified obesity-related single-nucleotide polymorphisms (SNPs) with body mass index (BMI) and associated traits. We then conducted association analyses between a genetic risk score (GRS) comprising of 13 replicated SNPs and the individual SNPs, and relative dietary intakes of fat, carbohydrates, protein, fiber and total energy intake, as well as interaction analyses on BMI and associated traits among 26,107 nondiabetic MDCS participants. GRS associated strongly with increased BMI (P = 3.6 × 10^?34), fat mass (P = 6.3 × 10^?28) and fat-free mass (P = 1.3 × 10^?24). Higher GRS associated with lower total energy intake (P = 0.001) and higher intake of fiber (P = 2.3 × 10^?4). No significant interactions were observed between GRS and the studied dietary intakes on BMI or related traits. Of the individual SNPs, after correcting for multiple comparisons, NEGR1 rs2815752 associated with diet intakes and BDNF rs4923461 showed interaction with protein intake on BMI. In conclusion, our study does not provide evidence for a major role for macronutrient-, fiber- or total energy intake levels in modifying genetic susceptibility to obesity measured as GRS. However, our data suggest that the number of risk alleles as well as some of the individual obesity loci may have a role in regulation of food and energy intake and that some individual loci may interact with diet.     

Assessment of association of rs2200733 on chromosome 4q25 with atrial fibrillation and ischemic stroke in a Chinese Han population

Atrial fibrillation (AF) is the most common arrhythmia in the clinical setting and an independent risk factor for stroke. Approximately 10 million Chinese people are affected by AF, but the genetic basis is largely unknown. A recent genome-wide association study in Iceland identified association between SNP rs2200733 on 4q25 and AF; however, many independent replication studies are essential to unequivocally validate this association. To assess the association between rs2200733 and AF as well as that between rs2200733 and ischemic stroke in a mainland Chinese Han population, we carried out case–control association studies with 383 AF patients versus 851 non-AF controls and 811 ischemic stroke patients versus 688 non-stroke controls. Highly significant association was detected between rs2200733 and AF in a Chinese Han population (allelic P = 3.7 × 10^?11 with OR = 1.81; genotypic P = 4.1 × 10^?12 with a dominant model). When the AF cases were divided into lone AF (32.6%) and other types of AF (67.4%), significantly stronger association was found with lone AF (OR = 2.40, P = 1.3 × 10^?9 compared to OR = 1.59, P = 6.2 × 10^?7 for other types of AF; P = 0.02 for two ORs). No significant association was found between rs2200733 and ischemic stroke. Our results suggest that SNP rs2200733 confers a highly significant risk of AF, but not ischemic stroke, in a more representative Chinese Han population in the mainland China.     

Genetic variants in the JAK1 gene confer higher risk of Behcet’s disease with ocular involvement in Han Chinese

Recent surveys have identified SLC22A4, SLC22A5, RUNX1, JAK1 as susceptibility genes for various immune-related diseases. An association study was performed in 738 Behcet’s patients with ocular involvement and 1,873 controls using the iPLEX system method. The first-stage study for 30 SNPs showed that SNPs rs2780815, rs310241, rs3790532 in JAK1 were associated with Behcet’s disease in Han Chinese (Pc_{(Bonferroni correction)} = 0.022–7.7 × 10^?3). The G allele and AA genotype of SNP rs2834643 in RUNX1 (Pc = 0.041–1.75 × 10^?3), but none of the other SNPs, were associated with Behcet’s disease. Haplotype analysis for the SLC22A4, SLC22A5 genes showed an increased tendency for AGTCTGCCGC frequency in patients compared with controls; however, the significance was lost after Bonferroni correction (P = 0.004, Pc > 0.05). Subsequently, we further replicated the significantly associated SNPs using another independent cohort. Replication and combining studies showed that three SNPs rs2780815, rs310241, rs3790532 in JAK1, but not SNP rs2834643 in RUNX1, were consistently associated with Behcet’s disease (replication: Pc = 0.012–9.60 × 10^?4; combining: Pc = 0.030–1.90 × 10^?4). SNPs rs2780815, rs310241, rs3790532 were estimated to confer a population attributable risk of 35.0, 28.0, 27.0 %, respectively. We found a strong association between HLA-B51 with Behcet’s disease in Chinese Han population (P = 1.35 × 10^?73; OR = 5.15; 95 % CI 4.28–6.19). GMDR analysis showed that no gene–gene interaction was detectable between JAK1 and HLA-B51. Logistic analysis indicated that the JAK1 gene was an independent risk factor for Behcet’s disease (P > 0.05). Real-time PCR analysis showed that no difference on the expression of JAK1 in PBMCs or LPS-stimulated PBMCs between individuals with the different rs1762780815 genotypes studied (P > 0.05). In conclusion, this study suggests that JAK1, but not SLC22A4, SLC22A5 and RUNX1, contributes to the genetic susceptibility to Behcet’s disease with ocular involvement.     

Association, interaction, and replication analysis of genes encoding serotonin transporter and 5-HT3 receptor subunits A and B in alcohol dependence

On the basis of the converging evidence showing regulation of drinking behavior by 5-HT_3AB receptors and the serotonin transporter, we hypothesized that the interactive effects of genetic variations in the genes HTR3A, HTR3B, and SLC6A4 confer greater susceptibility to alcohol dependence (AD) than do their effects individually. We examined the associations of AD with 22 SNPs across HTR3A, HTR3B, and two functional variants in SLC6A4 in 500 AD and 280 healthy control individuals of European descent. We found that the alleles of the low-frequency SNPs rs33940208:T in HTR3A and rs2276305:A in HTR3B were inversely and nominally significantly associated with AD with odds ratio (OR) and 95 % confidence interval of 0.212 and 0.073, 0.616 (P = 0.004) and 0.261 and 0.088, 0.777 (P = 0.016), respectively. Further, our gene-by-gene interaction analysis revealed that two four-variant models that differed by only one SNP carried a risk for AD (empirical P < 1 × 10^?6 for prediction accuracy of the two models based on 10⁶ permutations). Subsequent analysis of these two interaction models revealed an OR of 2.71 and 2.80, respectively, for AD (P < 0.001) in carriers of genotype combinations 5′-HTTLPR:LL/LS(SLC6A4)–rs1042173:TT/TG(SLC6A4)–rs1176744:AC(HTR3B)–rs3782025:AG(HTR3B) and 5′-HTTLPR:LL/LS(SLC6A4)–rs10160548:GT/TT(HTR3A)–rs1176744:AC(HTR3B)–rs3782025:AG(HTR3B). Combining all five genotypes resulted in an OR of 3.095 (P = 2.0 × 10^?4) for AD. Inspired by these findings, we conducted the analysis in an independent sample, OZ-ALC-GWAS (N = 6699), obtained from the NIH dbGAP database, which confirmed the findings, not only for all three risk genotype combinations (Z = 4.384, P = 1.0 × 10^?5; Z = 3.155, P = 1.6 × 10^?3; and Z = 3.389, P = 7.0 × 10^?4, respectively), but also protective effects for rs33940208:T (χ ² = 3.316, P = 0.0686) and rs2276305:A (χ ² = 7.224, P = 0.007). These findings reveal significant interactive effects among variants in SLC6A4–HTR3A–HTR3B affecting AD. Further studies are needed to confirm these findings and characterize the molecular mechanisms underlying these effects.     

Associations between interferon regulatory factor 5 polymorphisms and rheumatoid arthritis: a meta-analysis

The aim of this study was to determine whether interferon regulatory factor 5 (IRF5) polymorphisms confers susceptibility to rheumatoid arthritis (RA) in populations with different ethnicities. We searched the literature using the Pubmed and Embase databases and conducted meta-analyses on associations between the four IRF5 polymorphisms (rs2004640, rs729302, rs752637, and rs2280714) and RA susceptibility, using fixed and random effects models. A total of 12 comparison studies were considered in this meta-analysis, which in total involved 7,916 RA patients and 6,452 controls, and eight European, three Asian, and one Argentinean population. Meta-analysis showed an association between the minor allele of rs2004640 and RA in all subjects (odds ratio [OR] = 0.928, 95 % confidence interval [CI] = 0.865–0.996, P = 0.037). After stratification by ethnicity, analysis indicated that the minor allele was significantly associated with RA in Europeans (OR = 0.889, 95 % CI = 0.839–0.941, P = 5.03 × 10^?6), but not in Asians (OR = 1.057, 95 % CI = 0.978–1.144, P = 0.164). A direct comparison between anti-citrullinated peptide antibody-positive and -negative patients revealed no difference of the frequency of the rs2004640 minor allele (OR = 1.047, 95 % CI = 0.813–1.348, P = 0.724). Meta-analysis identified a significant association between RA and the minor allele of the rs729302 polymorphism in the overall population (OR = 0.896, 95 % CI = 0.826–0.972, P = 0.009) and in Asians (OR = 0.862, 95 % CI = 0.795–0.935, P = 3.50 × 10^?5), but not in Europeans (OR = 0.951, 95 % CI = 0.877–1.031, P = 0.225). Meta-analysis showed an association between the minor allele of rs752637 and RA in Europeans (OR = 0.858, 95 % CI = 0.789–0.932, P = 3.03 × 10^?5), but not in Asians (OR = 1.035, 95 % CI = 0.918–1.168, P = 0.572). No association was found between the rs2280714 polymorphism and RA susceptibility. This meta-analysis confirms that the IRF5 rs2004640, rs729302 and rs752637 polymorphisms are associated with RA susceptibility in different ethnic groups, especially in Europeans and Asians, but further study of this association is required in other ethnic groups.     

Genetic variants in microRNA and microRNA biogenesis pathway genes and breast cancer risk among women of African ancestry

MicroRNAs (miRNA) regulate breast biology by binding to specific RNA sequences, leading to RNA degradation and inhibition of translation of their target genes. While germline genetic variations may disrupt some of these interactions between miRNAs and their targets, studies assessing the relationship between genetic variations in the miRNA network and breast cancer risk are still limited, particularly among women of African ancestry. We systematically put together a list of 822 and 10,468 genetic variants among primary miRNA sequences and 38 genes in the miRNA biogenesis pathway, respectively; and examined their association with breast cancer risk in the ROOT consortium which includes women of African ancestry. Findings were replicated in an independent consortium. Logistic regression was used to estimate the odds ratio (OR) and 95 % confidence intervals (CI). For overall breast cancer risk, three single-nucleotide polymorphisms (SNPs) in miRNA biogenesis genes DROSHA rs78393591 (OR = 0.69, 95 % CI: 0.55–0.88, P = 0.003), ESR1 rs523736 (OR = 0.88, 95 % CI: 0.82–0.95, P = 3.99 × 10^?4), and ZCCHC11 rs114101502 (OR = 1.33, 95 % CI: 1.11–1.59, P = 0.002), and one SNP in primary miRNA sequence (rs116159732 in miR-6826, OR = 0.74, 95 % CI: 0.63–0.89, P = 0.001) were found to have significant associations in both discovery and validation phases. In a subgroup analysis, two SNPs were associated with risk of estrogen receptor (ER)-negative breast cancer, and three SNPs were associated with risk of ER-positive breast cancer. Several variants in miRNA and miRNA biogenesis pathway genes were associated with breast cancer risk. Risk associations varied by ER status, suggesting potential new mechanisms in etiology.     

Common variants in the ATP2B1 gene are associated with hypertension and arterial stiffness in Chinese population

Genome-wide association studies have identified the ATP2B1 gene associated with blood pressure (BP), the evidence from large scale Chinese population was still rare. We performed the current replication study to test the association of the ATP2B1 gene and hypertension and BP in two unrelated Chinese cohorts including 2,831 unrelated individuals with hypertension and 1,987 controls in total. We also examined the influences of the ATP2B1 gene on the arterial stiffness through evaluation of carotid-femoral pulse wave velocities (cf-PWV) in 164 untreated hypertensives. The major findings of this study were that four loci––rs10858911, rs2681472, rs17249754 and rs1401982––associated with any or all of four traits: hypertension (P = 0.001–4.6E–05; odds ratio, 0.83–0.87), systolic BP (P = 0.003–0.004), diastolic BP (P = 0.002–0.003) and cf-PWV (P = 0.002–0.004). All the comparisons were adjusted for sex, age, age² and body mass index. We validated the association of the ATP2B1 gene and susceptibility to hypertension, BP traits and cf-PWV in Chinese population. In addition, further genetic and functional research was warranted to elucidate the concrete locus in the ATP2B1 gene that influenced the manifestation of BP and vascular function.     

Association of four insulin resistance genes with type 2 diabetes mellitus and hypertension in the Chinese Han population

Insulin resistance plays an important role in the development of type 2 diabetes mellitus (T2DM) and hypertension. The purpose of the present study was to evaluate the association between four insulin resistance genes (ADIPOQ, LEPR, RETN, and TRIB3) and both T2DM and hypertension. A total of 768 Han Chinese subjects were recruited into this study, including 188 cases who had T2DM alone, 223 cases who had hypertension alone, 181 cases with both T2DM and hypertension, and 176 control subjects with neither T2DM nor hypertension. Twenty-three tag SNPs in four insulin resistance genes were genotyped and analyzed for association with T2DM and hypertension. One intron SNP (rs13306519) in LEPR and one 3′UTR SNP (rs1063537) in ADIPOQ demonstrated a significant association with T2DM (P = 0.024 and 0.014 respectively). Another intron SNP (rs12037879) in LEPR and a promoter region SNP (rs266729) in ADIPOQ were significantly associated with hypertension (P = 0.041 and 0.042, respectively). These associations survived the permutation test (P = 0.023, 0.018, 0.026, and 0.035, respectively). These associations were still found to be significant in the additive model after adjusting for potential confounding factors including age, sex, BMI, HDL, LDL, total cholesterol, and triglyceride levels (P = 0.024, 0.016, 0.04, and 0.043, respectively). No other gene variants were found to be significantly associated with T2DM or hypertension (P > 0.05). None of the studied gene variants were found to be significantly associated with T2DM+ hypertension (P > 0.05). A significant interaction was observed between two SNPs rs13306519 in LEPR and rs266729 in ADIPOQ for T2DM (P_int = 0.012, OR_int = 2.67) and hypertension (P_int = 0.0041, OR_int = 2.23). These findings suggest that variants in ADIPOQ and LEPR are risk factors for T2DM and hypertension in the Chinese population and that variants in RETN and TRIB3 are not major risk factors for these diseases.     

The combined impact of metabolic gene polymorphisms on elite endurance athlete status and related phenotypes

Endurance performance is a complex phenotype subject to the influence of both environmental and genetic factors. Although the last decade has seen a variety of specific genetic factors proposed, many in metabolic pathways, each is likely to make a limited contribution to an ‘elite’ phenotype: it seems more likely that such status depends on the simultaneous presence of multiple such variants. The aim of the study was to investigate individually and in combination the association of common metabolic gene polymorphisms with endurance athlete status, the proportion of slow-twitch muscle fibers and maximal oxygen consumption. A total of 1,423 Russian athletes and 1,132 controls were genotyped for 15 gene polymorphisms, of which most were previously reported to be associated with athlete status or related intermediate phenotypes. Muscle fiber composition of m. vastus lateralis in 45 healthy men was determined by immunohistochemistry. Maximal oxygen consumption of 50 male rowers of national competitive standard was determined during an incremental test to exhaustion on a rowing ergometer. Ten ‘endurance alleles’ (NFATC4 Gly160, PPARA rs4253778 G, PPARD rs2016520 C, PPARGC1A Gly482, PPARGC1B 203Pro, PPP3R1 promoter 5I, TFAM 12Thr, UCP2 55Val, UCP3 rs1800849 T and VEGFA rs2010963 C) were first identified showing discrete associations with elite endurance athlete status. Next, to assess the combined impact of all 10 gene polymorphisms, all athletes were classified according to the number of ‘endurance’ alleles they possessed. The proportion of subjects with a high (≥9) number of ‘endurance’ alleles was greater in the best endurance athletes compared with controls (85.7 vs. 37.8%, P = 7.6 × 10^?6). The number of ‘endurance’ alleles was shown to be positively correlated (r = 0.50; P = 4.0 × 10^?4) with the proportion of fatigue-resistant slow-twitch fibers, and with maximal oxygen consumption (r = 0.46; P = 7.0 × 10^?4). These data suggest that the likelihood of becoming an elite endurance athlete depends on the carriage of a high number of endurance-related alleles.     

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

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

Analysis of coding variants identified from exome sequencing resources for association with diabetic and non-diabetic nephropathy in African Americans

Prior studies have identified common genetic variants influencing diabetic and non-diabetic nephropathy, diseases which disproportionately affect African Americans. Recently, exome sequencing techniques have facilitated identification of coding variants on a genome-wide basis in large samples. Exonic variants in known or suspected end-stage kidney disease (ESKD) or nephropathy genes can be tested for their ability to identify association either singly or in combination with known associated common variants. Coding variants in genes with prior evidence for association with ESKD or nephropathy were identified in the NHLBI-ESP GO database and genotyped in 5,045 African Americans (3,324 cases with type 2 diabetes associated nephropathy [T2D-ESKD] or non-T2D ESKD, and 1,721 controls) and 1,465 European Americans (568 T2D-ESKD cases and 897 controls). Logistic regression analyses were performed to assess association, with admixture and APOL1 risk status incorporated as covariates. Ten of 31 SNPs were associated in African Americans; four replicated in European Americans. In African Americans, SNPs in OR2L8, OR2AK2, C6orf167 (MMS22L), LIMK2, APOL3, APOL2, and APOL1 were nominally associated (P = 1.8 × 10^?4–0.044). Haplotype analysis of common and coding variants increased evidence of association at the OR2L13 and APOL1 loci (P = 6.2 × 10^?5 and 4.6 × 10^?5, respectively). SNPs replicating in European Americans were in OR2AK2, LIMK2, and APOL2 (P = 0.0010-0.037). Meta-analyses highlighted four SNPs associated in T2D-ESKD and all-cause ESKD. Results from this study suggest a role for coding variants in the development of diabetic, non-diabetic, and/or all-cause ESKD in African Americans and/or European Americans.     

A functional variant of pre-miRNA-196a2 confers risk for Behcet’s disease but not for Vogt–Koyanagi–Harada syndrome or AAU in ankylosing spondylitis

This study aimed to investigate the predisposition of common pre-miRNA SNPs with Behcet’s disease (BD), Vogt–Koyanagi–Harada (VKH) syndrome and acute anterior uveitis (AAU) associated with ankylosing spondylitis (AS). A two-stage association study was carried out in 859 BD, 400 VKH syndrome, 209 AAU⁺AS⁺ patients and 1,685 controls all belonging to a Chinese Han population. Genotyping, the expression of miR-196a and Bach1 (the target gene of miR-196a), cell proliferation, cytokine production were examined by PCR–RFLP, real-time PCR, CCK8 and ELISA. In the first stage study, the results showed significantly increased frequencies of the miR-196a2/rs11614913 TT genotype and T allele in BD patients (adjusted P _c = 0.024, OR = 1.63; adjusted P _c = 5.4 × 10^?3, OR = 1.45, respectively). However, no significant association of the tested SNPs with VKH and AAU⁺AS⁺ patients was observed. The second stage and combined studies confirmed the association of rs11614913 with BD (TT genotype: adjusted P _c = 6×10^?5, OR = 1.53; T allele: adjusted P _c = 8×10^?6, OR = 1.35; CC genotype: adjusted P _c = 0.024, OR = 0.68). A stratified analysis showed an association of the rs11614913 TT genotype and T allele with the arthritis subgroup of BD (P _c = 5.3 × 10^?3, OR = 1.89; P _c = 0.015, OR = 1.56, respectively). Functional experiments showed a decreased miR-196a expression, an increased Bach1 expression and an increased production of IL-1β and MCP-1 in TT cases compared to CC cases (P = 0.023, P = 0.0073, P = 0.012, P = 0.002, respectively). This study shows that a functional variant of miR-196a2 confers risk for BD but not for VKH syndrome or AAU⁺AS⁺ by modulating the miR-196a gene expression and by regulating pro-inflammatory IL-1β and MCP-1 production.     

A Genome-wide Expression Association Analysis Identifies Genes and Pathways Associated with Amyotrophic Lateral Sclerosis

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease with strong genetic components. To identity novel risk variants for ALS, utilizing the latest genome-wide association studies (GWAS) and eQTL study data, we conducted a genome-wide expression association analysis by summary data-based Mendelian randomization (SMR) method. Summary data were derived from a large-scale GWAS of ALS, involving 12577 cases and 23475 controls. The eQTL annotation dataset included 923,021 cis-eQTL for 14,329 genes and 4732 trans-eQTL for 2612 genes. Genome-wide single gene expression association analysis was conducted by SMR software. To identify ALS-associated biological pathways, the SMR analysis results were further subjected to gene set enrichment analysis (GSEA). SMR single gene analysis identified one significant and four suggestive genes associated with ALS, including C9ORF72 (P value = 7.08 × 10^?6), NT5C3L (P value = 1.33 × 10^?5), GGNBP2 (P value = 1.81 × 10^?5), ZNHIT3(P value = 2.94 × 10^?5), and KIAA1600(P value = 9.97 × 10^?5). GSEA identified 7 significant biological pathways, such as PEROXISOME (empirical P value = 0.006), GLYCOLYSIS_GLUCONEOGENESIS (empirical P value = 0.043), and ARACHIDONIC_ACID_ METABOLISM (empirical P value = 0.040). Our study provides novel clues for the genetic mechanism studies of ALS.     

Identification of methylation quantitative trait loci (mQTLs) influencing promoter DNA methylation of alcohol dependence risk genes

Interaction of DNA methylation and sequence variants that are methylation quantitative trait loci (mQTLs) may influence susceptibility to diseases such as alcohol dependence (AD). We used genome-wide genotype data from 268 African Americans (AAs: 129 AD cases and 139 controls) and 143 European Americans (EAs: 129 AD cases and 14 controls) to identify mQTLs that were associated with promoter CpGs in 82 AD risk genes. 282 significant mQTL–CpG pairs (9.9 × 10^?100 ≤ P _nominal ≤ 7.7 × 10^?8) in AAs and 313 significant mQTL–CpG pairs (2.7 × 10^?53 ≤ P _nominal ≤ 9.9 × 10^?8) in EAs were identified [i.e., mQTL–CpG associations survived multiple-testing correction, q values (false discovery rate) ≤ 0.05]. The most significant mQTL was rs1800759, which was strongly associated with CpG cg12011299 in both AAs (P _nominal = 9.9 × 10^?100; q = 6.7 × 10^?91) and EAs (P _nominal = 2.7 × 10^?53; q = 1.4 × 10^?44). Rs1800759 (previously known to be associated to AD) and CpG cg12011299 (distance: 37 bp) are both located in alcohol dehydrogenase (ADH) 4 gene (ADH4) promoter region. In general, the strength of association between mQTLs and CpGs was inversely correlated with the distance between them. Association was also influenced by race and AD. Additionally, 48.3 % of the mQTLs identified in AAs and 65.6 % of the mQTLs identified in EAs were predicted to be expression QTLs. Three mQTLs (rs2173201, rs4147542, and rs4147541 in ADH1B-AHD1C gene cluster region) found in AAs were previously identified by our genome-wide association studies as being significantly associated with AD in AAs. Thus, DNA methylation, which can be influenced by sequence variants and is implicated in gene expression regulation, appears to at least partially underlie the association of genetic variation with AD.     

The majority of the genetic risk for Paget’s disease of bone is explained by genetic variants close to the CSF1, OPTN, TM7SF4, and TNFRSF11A genes

Paget’s disease of bone (PDB) is one of the most frequent metabolic bone disorders (1–5%), next to osteoporosis, affecting individuals above age 55. Sequestosome1 mutations explain a part of the PDB patients, but still the disease pathogenesis in the remaining PDB patients is largely unknown. Therefore, association studies investigating the relationship between genetic polymorphisms and sporadic PDB have been performed to find the genetic risk variants. Previously such studies indicated a role of the OPG and RANK gene. The latter was recently confirmed in a genome-wide association study (GWAS) which also indicated the involvement of chromosomal regions harbouring the CSF1 and OPTN gene. In this study, we sought to replicate these findings in a Belgian and a Dutch population. Similar significant results were obtained for the single nucleotide polymorphisms and the haplotypes. The most significant results are found in the CSF1 gene region, followed by the OPTN and TNFRSF11A gene region (p values ranging from 1.3 × 10^?4 to 3.8 × 10^?8, OR = 1.523–1.858). We next obtained significant association with a polymorphism from the chromosomal region around the TM7SF4 gene (p = 2.7 × 10^?3, OR = 1.427), encoding DC-STAMP which did not reach genome-wide significance in the GWAS, but based on its function in osteoclasts it can be considered a strong candidate gene. After meta-analysis with the GWAS data, p values ranged between 2.6 × 10^?4 and 8.8 × 10^?32. The calculated cumulative population attributable risk of these four loci turned out to be about 67% in our two populations, indicating that most of the genetic risk for PDB is coming from genetic variants close to these four genes.     

A genome-wide association identified the common genetic variants influence disease severity in β0-thalassemia/hemoglobin E

β-Thalassemia/HbE disease is clinically variable. In searching for genetic factors modifying the disease severity, patients were selected based on their disease severities, and a genome-wide association study (GWAS) was performed. Genotyping was conducted with the Illumina Human 610-Quad BeadChips array using DNAs from 618 Thai β⁰-thalassemia/HbE patients who were classified as 383 severe and 235 mild phenotypes by a validated scoring system. Twenty-three SNPs in three independent genes/regions were identified as being significantly associated with the disease severity. The highest association was observed with SNPs in the β-globin gene cluster (chr.11p15), and rs2071348 of the HBBP1 gene revealed the most significant association [P = 2.96 × 10^?13, odds ratio (OR) = 4.33 (95% confidence interval (CI), 2.74–6.84)]. The second was identified in the intergenic region between the HBS1L and MYB genes (chr.6q23), among which rs9376092 was the most significant [P = 2.36 × 10^?10, OR = 3.07 (95% CI, 2.16–4.38)]. The third region was located in the BCL11A gene (chr.2p16.1), and rs766432 showed the most significant association [P = 5.87 × 10^?10, OR = 3.06 (95% CI, 2.15–4.37)]. All three loci were replicated in an independent cohort of 174 Indonesian patients. The associations to fetal hemoglobin levels were also observed with SNPs on these three regions. Our data indicate that several genetic loci act in concert to influence HbF levels of β⁰-thalassemia/HbE patients. This study revealed that all the three reported loci and the α-globin gene locus are the best and common predictors of the disease severity in β-thalassemia.