Genome-wide meta-analysis identifies six novel loci associated with habitual coffee consumption

Coffee, a major dietary source of caffeine, is among the most widely consumed beverages in the world and has received considerable attention regarding health risks and benefits. We conducted a genome-wide (GW) meta-analysis of predominately regular-type coffee consumption (cups per day) among up to 91 462 coffee consumers of European ancestry with top single-nucleotide polymorphisms (SNPs) followed-up in ~30 062 and 7964 coffee consumers of European and African-American ancestry, respectively. Studies from both stages were combined in a trans-ethnic meta-analysis. Confirmed loci were examined for putative functional and biological relevance. Eight loci, including six novel loci, met GW significance (log10Bayes factor (BF)>5.64) with per-allele effect sizes of 0.03-0.14 cups per day. Six are located in or near genes potentially involved in pharmacokinetics (ABCG2, AHR, POR and CYP1A2) and pharmacodynamics (BDNF and SLC6A4) of caffeine. Two map to GCKR and MLXIPL genes related to metabolic traits but lacking known roles in coffee consumption. Enhancer and promoter histone marks populate the regions of many confirmed loci and several potential regulatory SNPs are highly correlated with the lead SNP of each. SNP alleles near GCKR, MLXIPL, BDNF and CYP1A2 that were associated with higher coffee consumption have previously been associated with smoking initiation, higher adiposity and fasting insulin and glucose but lower blood pressure and favorable lipid, inflammatory and liver enzyme profiles (P<5 × 10-8).Our genetic findings among European and African-American adults reinforce the role of caffeine in mediating habitual coffee consumption and may point to molecular mechanisms underlying inter-individual variability in pharmacological and health effects of coffee.

Mapping Bias Overestimates Reference Allele Frequencies at the HLA Genes in the 1000 Genomes Project Phase I Data.

Next-generation sequencing (NGS) technologies have become the standard for data generation in studies of population genomics, as the 1000 Genomes Project (1000G). However, these techniques are known to be problematic when applied to highly polymorphic genomic regions, such as the human leukocyte antigen (HLA) genes. Because accurate genotype calls and allele frequency estimations are crucial to population genomics analyses, it is important to assess the reliability of NGS data. Here, we evaluate the reliability of genotype calls and allele frequency estimates of the single-nucleotide polymorphisms (SNPs) reported by 1000G (phase I) at five HLA genes (HLA-A, -B, -C, -DRB1, and -DQB1). We take advantage of the availability of HLA Sanger sequencing of 930 of the 1092 1000G samples and use this as a gold standard to benchmark the 1000G data. We document that 18.6% of SNP genotype calls in HLA genes are incorrect and that allele frequencies are estimated with an error greater than ±0.1 at approximately 25% of the SNPs in HLA genes. We found a bias toward overestimation of reference allele frequency for the 1000G data, indicating mapping bias is an important cause of error in frequency estimation in this dataset. We provide a list of sites that have poor allele frequency estimates and discuss the outcomes of including those sites in different kinds of analyses. Because the HLA region is the most polymorphic in the human genome, our results provide insights into the challenges of using of NGS data at other genomic regions of high diversity.

Molecular Analysis of NOTCH2 in Patients with Primary Open-Angle Glaucoma.

BACKGROUND AND PURPOSE: Transgenic mice overexpressing Notch2 in the uvea exhibit a hyperplastic ciliary body leading to increased IOP and glaucoma. The aim of this study was to investigate the possible presence of NOTCH2 variants in patients with primary open-angle glaucoma (POAG). METHODS: We screened DNA samples from 130 patients with POAG for NOTCH2 variants by denaturing high-performance liquid chromatography after PCR amplification and validated our data by direct Sanger sequencing. RESULTS: No mutations were observed in the coding regions of NOTCH2 or in the splice sites. 19 known SNPs (single nucleotide polymorphisms) were detected. An SNP located in intron 24, c.[4005+45A>G], was seen in 28.5% of the patients (37/130 patients). As this SNP is reported to have a minor allele frequency of 7% in the 1000 genomes database, it could be associated with POAG. However, we evaluated its frequency in an ethnic-matched control group of 96 subjects unaffected by POAG and observed a frequency of 29%, indicating that it was not related to POAG. CONCLUSION: NOTCH2 seemed to be a good candidate for POAG as it is expressed in the anterior segment in the human eye. However, mutational analysis did not show any causative mutation. This study also shows that proper ethnic-matched control groups are essential in association studies and that values given in databases are sometimes misleading.

Influence of MCHR2 and MCHR2-AS1 Genetic Polymorphisms on Body Mass Index in Psychiatric Patients and In Population-Based Subjects with Present or Past Atypical Depression.

Obesity development during psychotropic treatments represents a major health issue in psychiatry. Melanin-concentrating hormone receptor 2 (MCHR2) is a central receptor involved in energy homeostasis. MCHR2 shares its promoter region with MCHR2-AS1, a long antisense non-coding RNA. The aim of this study was to determine whether tagging single nucleotide polymorphisms (tSNPs) of MCHR2 and MCHR2-AS1 are associated with the body mass index (BMI) in the psychiatric and in the general population. The influence of MCHR2 and MCHR2-AS1 tSNPs on BMI was firstly investigated in a discovery psychiatric sample (n1 = 474). Positive results were tested for replication in two other psychiatric samples (n2 = 164, n3 = 178) and in two population-based samples (CoLaus, n4 = 5409; GIANT, n5 = 113809). In the discovery sample, TT carriers of rs7754794C>T had 1.08 kg/m2 (p = 0.04) lower BMI as compared to C-allele carriers. This observation was replicated in an independent psychiatric sample (-2.18 kg/m2; p = 0.009). The association of rs7754794C>T and BMI seemed stronger in subjects younger than 45 years (median of age). In the population-based sample, a moderate association was observed (-0.17 kg/m2; p = 0.02) among younger individuals (<45y). Interestingly, this association was totally driven by patients meeting lifetime criteria for atypical depression, i.e. major depressive episodes characterized by symptoms such as an increased appetite. Indeed, patients with atypical depression carrying rs7754794-TT had 1.17 kg/m2 (p = 0.04) lower BMI values as compared to C-allele carriers, the effect being stronger in younger individuals (-2.50 kg/m2; p = 0.03; interaction between rs7754794 and age: p-value = 0.08). This study provides new insights on the possible influence of MCHR2 and/or MCHR2-AS1 on obesity in psychiatric patients and on the pathophysiology of atypical depression.

Impact of common risk factors of fibrosis progression in chronic hepatitis C.

OBJECTIVE: The natural course of chronic hepatitis C varies widely. To improve the profiling of patients at risk of developing advanced liver disease, we assessed the relative contribution of factors for liver fibrosis progression in hepatitis C. DESIGN: We analysed 1461 patients with chronic hepatitis C with an estimated date of infection and at least one liver biopsy. Risk factors for accelerated fibrosis progression rate (FPR), defined as ≥0.13 Metavir fibrosis units per year, were identified by logistic regression. Examined factors included age at infection, sex, route of infection, HCV genotype, body mass index (BMI), significant alcohol drinking (≥20 g/day for ≥5 years), HIV coinfection and diabetes. In a subgroup of 575 patients, we assessed the impact of single nucleotide polymorphisms previously associated with fibrosis progression in genome-wide association studies. Results were expressed as attributable fraction (AF) of risk for accelerated FPR. RESULTS: Age at infection (AF 28.7%), sex (AF 8.2%), route of infection (AF 16.5%) and HCV genotype (AF 7.9%) contributed to accelerated FPR in the Swiss Hepatitis C Cohort Study, whereas significant alcohol drinking, anti-HIV, diabetes and BMI did not. In genotyped patients, variants at rs9380516 (TULP1), rs738409 (PNPLA3), rs4374383 (MERTK) (AF 19.2%) and rs910049 (major histocompatibility complex region) significantly added to the risk of accelerated FPR. Results were replicated in three additional independent cohorts, and a meta-analysis confirmed the role of age at infection, sex, route of infection, HCV genotype, rs738409, rs4374383 and rs910049 in accelerating FPR. CONCLUSIONS: Most factors accelerating liver fibrosis progression in chronic hepatitis C are unmodifiable.

Does phylogeographic structure relate to climatic niche divergence? A test using maritime pine (Pinus pinaster Ait.)

Aim To disentangle the effects of environmental and geographical processes driving phylogenetic distances among clades of maritime pine (Pinus pinaster). To assess the implications for conservation management of combining molecular information with species distribution models (SDMs; which predict species distribution based on known occurrence records and on environmental variables). Location Western Mediterranean Basin and European Atlantic coast. Methods We undertook two cluster analyses for eight genetically defined pine clades based on climatic niche and genetic similarities. We assessed niche similarity by means of a principal component analysis and Schoener's D metric. To calculate genetic similarity, we used the unweighted pair group method with arithmetic mean based on Nei's distance using 266 single nucleotide polymorphisms. We then assessed the contribution of environmental and geographical distances to phylogenetic distance by means of Mantel regression with variance partitioning. Finally, we compared the projection obtained from SDMs fitted from the species level (SDMsp) and composed from the eight clade-level models (SDMcm). Results Genetically and environmentally defined clusters were identical. Environmental and geographical distances explained 12.6% of the phylogenetic distance variation and, overall, geographical and environmental overlap among clades was low. Large differences were detected between SDMsp and SDMcm (57.75% of disagreement in the areas predicted as suitable). Main conclusions The genetic structure within the maritime pine subspecies complex is primarily a consequence of its demographic history, as seen by the high proportion of unexplained variation in phylogenetic distances. Nevertheless, our results highlight the contribution of local environmental adaptation in shaping the lower-order, phylogeographical distribution patterns and spatial genetic structure of maritime pine: (1) genetically and environmentally defined clusters are consistent, and (2) environment, rather than geography, explained a higher proportion of variation in phylogenetic distance. SDMs, key tools in conservation management, better characterize the fundamental niche of the species when they include molecular information.

Buried in the Middle but Guilty: Intronic Mutations in the TCIRG1 Gene Cause Human Autosomal Recessive Osteopetrosis.

Autosomal recessive osteopetrosis (ARO) is a rare genetic bone disease with genotypic and phenotypic heterogeneity, sometimes translating into delayed diagnosis and treatment. In particular, cases of intermediate severity often constitute a diagnostic challenge and represent good candidates for exome sequencing. Here, we describe the tortuous path to identification of the molecular defect in two siblings, in which osteopetrosis diagnosed in early childhood followed a milder course, allowing them to reach the adult age in relatively good conditions with no specific therapy. No clearly pathogenic mutation was identified either with standard amplification and resequencing protocols or with exome sequencing analysis. While evaluating the possible impact of a 3'UTR variant on the TCIRG1 expression, we found a novel single nucleotide change buried in the middle of intron 15 of the TCIRG1 gene, about 150 nucleotides away from the closest canonical splice site. By sequencing a number of independent cDNA clones covering exons 14 to 17, we demonstrated that this mutation reduced splicing efficiency but did not completely abrogate the production of the normal transcript. Prompted by this finding, we sequenced the same genomic region in 33 patients from our unresolved ARO cohort and found three additional novel single nucleotide changes in a similar location and with a predicted disruptive effect on splicing, further confirmed in one of them at the transcript level. Overall, we identified an intronic region in TCIRG1 that seems to be particularly prone to splicing mutations, allowing the production of a small amount of protein sufficient to reduce the severity of the phenotype usually associated with TCIRG1 defects. On this basis, we would recommend including TCIRG1 not only in the molecular work-up of severe infantile osteopetrosis but also in intermediate cases and carefully evaluating the possible effects of intronic changes. © 2015 American Society for Bone and Mineral Research.

Association between parental history and genetic risk scores for coronary heart disease prediction: The population-based CoLaus study.

BACKGROUND AND AIMS: Parental history (PH) and genetic risk scores (GRSs) are separately associated with coronary heart disease (CHD), but evidence regarding their combined effects is lacking. We aimed to evaluate the joint associations and predictive ability of PH and GRSs for incident CHD. METHODS: Data for 4283 Caucasians were obtained from the population-based CoLaus Study, over median follow-up time of 5.6 years. CHD was defined as incident myocardial infarction, angina, percutaneous coronary revascularization or bypass grafting. Single nucleotide polymorphisms for CHD identified by genome-wide association studies were used to construct unweighted and weighted versions of three GRSs, comprising of 38, 53 and 153 SNPs respectively. RESULTS: PH was associated with higher values of all weighted GRSs. After adjustment for age, sex, smoking, diabetes, systolic blood pressure, low and high density lipoprotein cholesterol, PH was significantly associated with CHD [HR 2.61, 95% CI (1.47-4.66)] and further adjustment for GRSs did not change this estimate. Similarly, one standard deviation change of the weighted 153-SNPs GRS was significantly associated with CHD [HR 1.50, 95% CI (1.26-1.80)] and remained so, after further adjustment for PH. The weighted, 153-SNPs GRS, but not PH, modestly improved discrimination [(C-index improvement, 0.016), p = 0.048] and reclassification [(NRI improvement, 8.6%), p = 0.027] beyond cardiovascular risk factors. After including both the GRS and PH, model performance improved further [(C-index improvement, 0.022), p = 0.006]. CONCLUSION: After adjustment for cardiovascular risk factors, PH and a weighted, polygenic GRS were jointly associated with CHD and provided additive information for coronary events prediction.

Impact of a cis-associated gene expression SNP on chromosome 20q11.22 on bipolar disorder susceptibility, hippocampal structure and cognitive performance.

BackgroundBipolar disorder is a highly heritable polygenic disorder. Recent enrichment analyses suggest that there may be true risk variants for bipolar disorder in the expression quantitative trait loci (eQTL) in the brain.AimsWe sought to assess the impact of eQTL variants on bipolar disorder risk by combining data from both bipolar disorder genome-wide association studies (GWAS) and brain eQTL.MethodTo detect single nucleotide polymorphisms (SNPs) that influence expression levels of genes associated with bipolar disorder, we jointly analysed data from a bipolar disorder GWAS (7481 cases and 9250 controls) and a genome-wide brain (cortical) eQTL (193 healthy controls) using a Bayesian statistical method, with independent follow-up replications. The identified risk SNP was then further tested for association with hippocampal volume (n = 5775) and cognitive performance (n = 342) among healthy individuals.ResultsIntegrative analysis revealed a significant association between a brain eQTL rs6088662 on chromosome 20q11.22 and bipolar disorder (log Bayes factor = 5.48; bipolar disorder P = 5.85×10(-5)). Follow-up studies across multiple independent samples confirmed the association of the risk SNP (rs6088662) with gene expression and bipolar disorder susceptibility (P = 3.54×10(-8)). Further exploratory analysis revealed that rs6088662 is also associated with hippocampal volume and cognitive performance in healthy individuals.ConclusionsOur findings suggest that 20q11.22 is likely a risk region for bipolar disorder; they also highlight the informative value of integrating functional annotation of genetic variants for gene expression in advancing our understanding of the biological basis underlying complex disorders, such as bipolar disorder.

Genome-wide Association Studies Identify Genetic Loci Associated With Albuminuria in Diabetes.

Elevated concentrations of albumin in the urine, albuminuria, are a hallmark of diabetic kidney disease and are associated with an increased risk for end-stage renal disease and cardiovascular events. To gain insight into the pathophysiological mechanisms underlying albuminuria, we conducted meta-analyses of genome-wide association studies and independent replication in up to 5,825 individuals of European ancestry with diabetes and up to 46,061 without diabetes, followed by functional studies. Known associations of variants in CUBN, encoding cubilin, with the urinary albumin-to-creatinine ratio (UACR) were confirmed in the overall sample (P = 2.4 × 10(-10)). Gene-by-diabetes interactions were detected and confirmed for variants in HS6ST1 and near RAB38/CTSC. Single nucleotide polymorphisms at these loci demonstrated a genetic effect on UACR in individuals with but not without diabetes. The change in the average UACR per minor allele was 21% for HS6ST1 (P = 6.3 × 10(-7)) and 13% for RAB38/CTSC (P = 5.8 × 10(-7)). Experiments using streptozotocin-induced diabetic Rab38 knockout and control rats showed higher urinary albumin concentrations and reduced amounts of megalin and cubilin at the proximal tubule cell surface in Rab38 knockout versus control rats. Relative expression of RAB38 was higher in tubuli of patients with diabetic kidney disease compared with control subjects. The loci identified here confirm known pathways and highlight novel pathways influencing albuminuria.

Prevalence and determinants of periodic limb movements in the general population.

OBJECTIVE: Periodic limb movements during sleep (PLMS) are sleep phenomena characterized by periodic episodes of repetitive stereotyped limb movements. The aim of this study was to describe the prevalence and determinants of PLMS in a middle to older aged general population. METHODS: Data from 2,162 subjects (51.2% women, mean age = 58.4 ± 11.1 years) participating in a population-based study (HypnoLaus, Lausanne, Switzerland) were collected. Assessments included laboratory tests, sociodemographic data, personal and treatment history, and full polysomnography at home. PLMS index (PLMSI) was determined, and PLMSI > 15/h was considered as significant. RESULTS: Prevalence of PLMSI > 15/h was 28.6% (31.3% in men, 26% in women). Compared to subjects with PLMSI ≤ 15/h, subjects with PLMSI > 15/h were older (p < 0.001), were predominantly males (p = 0.007), had a higher proportion of restless legs syndrome (RLS; p < 0.001), had a higher body mass index (p = 0.001), and had a lower mean glomerular filtration rate (p < 0.001). Subjects with PLMSI > 15/h also had a higher prevalence of diabetes, hypertension, and beta-blocker or hypnotic treatments. The prevalence of antidepressant use was higher, but not statistically significant (p = 0.07). Single nucleotide polymorphisms (SNPs) within BTBD9 (rs3923809), TOX3 (rs3104788), and MEIS1 (rs2300478) genes were significantly associated with PLSMI > 15/h. Conversely, mean hemoglobin and ferritin levels were similar in both groups. In the multivariate analysis, age, male gender, antidepressant intake, RLS, and rs3923809, rs3104788, and rs2300478 SNPs were independently associated with PLMSI > 15/h. INTERPRETATION: PLMS are highly prevalent in our middle-aged European population. Age, male gender, RLS, antidepressant treatment, and specific BTBD9, TOX3, and MEIS1 SNP distribution are independent predictors of PLMSI > 15/h. ANN NEUROL 2016;79:464-474.

Hybridization Capture Using RAD Probes (hyRAD), a New Tool for Performing Genomic Analyses on Collection Specimens.

In the recent years, many protocols aimed at reproducibly sequencing reduced-genome subsets in non-model organisms have been published. Among them, RAD-sequencing is one of the most widely used. It relies on digesting DNA with specific restriction enzymes and performing size selection on the resulting fragments. Despite its acknowledged utility, this method is of limited use with degraded DNA samples, such as those isolated from museum specimens, as these samples are less likely to harbor fragments long enough to comprise two restriction sites making possible ligation of the adapter sequences (in the case of double-digest RAD) or performing size selection of the resulting fragments (in the case of single-digest RAD). Here, we address these limitations by presenting a novel method called hybridization RAD (hyRAD). In this approach, biotinylated RAD fragments, covering a random fraction of the genome, are used as baits for capturing homologous fragments from genomic shotgun sequencing libraries. This simple and cost-effective approach allows sequencing of orthologous loci even from highly degraded DNA samples, opening new avenues of research in the field of museum genomics. Not relying on the restriction site presence, it improves among-sample loci coverage. In a trial study, hyRAD allowed us to obtain a large set of orthologous loci from fresh and museum samples from a non-model butterfly species, with a high proportion of single nucleotide polymorphisms present in all eight analyzed specimens, including 58-year-old museum samples. The utility of the method was further validated using 49 museum and fresh samples of a Palearctic grasshopper species for which the spatial genetic structure was previously assessed using mtDNA amplicons. The application of the method is eventually discussed in a wider context. As it does not rely on the restriction site presence, it is therefore not sensitive to among-sample loci polymorphisms in the restriction sites that usually causes loci dropout. This should enable the application of hyRAD to analyses at broader evolutionary scales.