The role of plasma concentrations and drug characteristics of beta-blockers in fall risk of older persons.

Beta-blocker usage is inconsistently associated with increased fall risk in the literature. However, due to age-related changes and interindividual heterogeneity in pharmacokinetics and dynamics, it is difficult to predict which older adults are more at risk for falls. Therefore, we wanted to explore whether elevated plasma concentrations of selective and nonselective beta-blockers are associated with an increased risk of falls in older beta-blocker users. To answer our research question, we analyzed samples of selective (metoprolol, n = 316) and nonselective beta-blockers (sotalol, timolol, propranolol, and carvedilol, n = 179) users from the B-PROOF cohort. The associations between the beta-blocker concentration and time to first fall were assessed using Cox proportional hazard models. Change of concentration over time in relation to fall risk was assessed with logistic regression models. Models were adjusted for potential confounders. Our results showed that above the median concentration of metoprolol was associated with an increased fall risk (HR 1.55 [1.11-2.16], p = .01). No association was found for nonselective beta-blocker concentrations. Also, changes in concentration over time were not associated with increased fall risk. To conclude, metoprolol plasma concentrations were associated with an increased risk of falls in metoprolol users while no associations were found for nonselective beta-blockers users. This might be caused by a decreased ß1-selectivity in high plasma concentrations. In the future, beta-blocker concentrations could potentially help clinicians estimate fall risk in older beta-blockers users and personalize treatment.

Novel Genetic Determinants of Dental Maturation in Children.

Dental occlusion requires harmonious development of teeth, jaws, and other elements of the craniofacial complex, which are regulated by environmental and genetic factors. We performed the first genome-wide association study (GWAS) on dental development (DD) using the Demirjian radiographic method. Radiographic assessments from participants of the Generation R Study (primary study population, N1 = 2,793; mean age of 9.8 y) were correlated with ~30 million genetic variants while adjusting for age, sex, and genomic principal components (proxy for population stratification). Variants associated with DD at genome-wide significant level (P < 5 × 10-8) mapped to 16q12.2 (IRX5) (lead variant rs3922616, B = 0.16; P = 2.2 × 10-8). We used Fisher's combined probability tests weighted by sample size to perform a meta-analysis (N = 14,805) combining radiographic DD at a mean age of 9.8 y from Generation R with data from a previous GWAS (N2 = 12,012) on number of teeth (NT) in infants used as proxy of DD at a mean age of 9.8 y (including the ALSPAC and NFBC1966). This GWAS meta-analysis revealed 3 novel loci mapping to 7p15.3 (IGF2BP3: P = 3.2 × 10-8), 14q13.3 (PAX9: P = 1.9 × 10-8), and 16q12.2 (IRX5: P = 1.2 × 10-9) and validated 8 previously reported NT loci. A polygenic allele score constructed from these 11 loci was associated with radiographic DD in an independent Generation R set of children (N = 703; B = 0.05, P = 0.004). Furthermore, profiling of the identified genes across an atlas of murine and human stem cells observed expression in the cells involved in the formation of bone and/or dental tissues (>0.3 frequency per kilobase of transcript per million mapped reads), likely reflecting functional specialization. Our findings provide biological insight into the polygenic architecture of the pediatric dental maturation process.

Candidate genetic variants and antidepressant-related fall risk in middle-aged and older adults.

BACKGROUND: Antidepressant use has been associated with increased fall risk. Antidepressant-related adverse drug reactions (e.g. orthostatic hypotension) depend partly on genetic variation. We hypothesized that candidate genetic polymorphisms are associated with fall risk in older antidepressant users. METHODS: The association between antidepressant use and falls was cross-sectionally investigated in a cohort of Dutch older adults by logistic regression analyses. In case of significant interaction product term of antidepressant use and candidate polymorphism, the association between the variant genotype and fall risk was assessed within antidepressant users and the association between antidepressant use and fall risk was investigated stratified per genotype. Secondly, a look-up of the candidate genes was performed in an existing genome-wide association study on drug-related falls in antidepressant users within the UK Biobank. In antidepressant users, genetic associations for our candidate polymorphisms for fall history were investigated. RESULTS: In antidepressant users(n = 566), for rs28371725 (CYP2D6*41) fall risk was decreased in TC/variant allele carriers compared to CC/non-variant allele carriers (OR = 0.45, 95% CI 0.26-0.80). Concerning rs1057910 (CYP2C9*3), fall risk was increased in CA/variant allele carriers compared to AA/non-variant allele carriers (OR = 1.95, 95% CI 1.17-3.27). Regarding, rs1045642 (ABCB1), fall risk was increased in AG/variant allele carriers compared to GG/non-variant allele carriers (OR = 1.69, 95% CI 1.07-2.69). Concerning the ABCB1-haplotype (rs1045642/rs1128503), fall risk was increased in AA-AA/variant allele carriers compared to GG-GG/non-variant allele carriers (OR = 1.86, 95% CI 1.05-3.29). In the UK Biobank, in antidepressant users(n = 34,000) T/variant-allele of rs28371725 (CYP2D*41) was associated with increased fall risk (OR = 1.06, 95% CI 1.01-1.12). G/non-variant-allele of rs4244285 (CY2C19*2) was associated with decreased risk (OR = 0.96, 95% CI 0.92-1.00). CONCLUSION: This is the first study showing that certain genetic variants modify antidepressant-related fall risk. The results were not always consistent across the studies and should be validated in a study with a prospective design. However, pharmacogenetics might have value in antidepressant (de)prescribing in falls prevention.

Genome-wide analysis of mitochondrial DNA copy number reveals loci implicated in nucleotide metabolism, platelet activation, and megakaryocyte proliferation.

Mitochondrial DNA copy number (mtDNA-CN) measured from blood specimens is a minimally invasive marker of mitochondrial function that exhibits both inter-individual and intercellular variation. To identify genes involved in regulating mitochondrial function, we performed a genome-wide association study (GWAS) in 465,809 White individuals from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium and the UK Biobank (UKB). We identified 133 SNPs with statistically significant, independent effects associated with mtDNA-CN across 100 loci. A combination of fine-mapping, variant annotation, and co-localization analyses was used to prioritize genes within each of the 133 independent sites. Putative causal genes were enriched for known mitochondrial DNA depletion syndromes (p = 3.09 × 10-15) and the gene ontology (GO) terms for mtDNA metabolism (p = 1.43 × 10-8) and mtDNA replication (p = 1.2 × 10-7). A clustering approach leveraged pleiotropy between mtDNA-CN associated SNPs and 41 mtDNA-CN associated phenotypes to identify functional domains, revealing three distinct groups, including platelet activation, megakaryocyte proliferation, and mtDNA metabolism. Finally, using mitochondrial SNPs, we establish causal relationships between mitochondrial function and a variety of blood cell-related traits, kidney function, liver function and overall (p = 0.044) and non-cancer mortality (p = 6.56 × 10-4).

TCERG1L allelic variation is associated with cisplatin-induced hearing loss in childhood cancer, a PanCareLIFE study.

In children with cancer, the heterogeneity in ototoxicity occurrence after similar treatment suggests a role for genetic susceptibility. Using a genome-wide association study (GWAS) approach, we identified a genetic variant in TCERG1L (rs893507) to be associated with hearing loss in 390 non-cranial irradiated, cisplatin-treated children with cancer. These results were replicated in two independent, similarly treated cohorts (n = 192 and 188, respectively) (combined cohort: P = 5.3 × 10-10, OR 3.11, 95% CI 2.2-4.5). Modulating TCERG1L expression in cultured human cells revealed significantly altered cellular responses to cisplatin-induced cytokine secretion and toxicity. These results contribute to insights into the genetic and pathophysiological basis of cisplatin-induced ototoxicity.

Prevalence of and risk factors for extended-spectrum beta-lactamase genes carriership in a population-based cohort of middle-aged and elderly.

INTRODUCTION: Increasing resistance to beta-lactam antibiotics is an alarming development worldwide. Fecal carriership of TEM, SHV, CTX-M and CMY was studied in a community-dwelling population of middle-aged and elderly individuals. PATIENTS AND METHODS: Feces was obtained from individuals of the Rotterdam Study. Carriership of the TEM, SHV, CTX-M and CMY genes was determined using real-time polymerase chain reaction (qPCR). Possible associations were investigated between carriership of these genes and several risk factors, such as the use of antimicrobial drugs, diabetes mellitus, protein pump inhibitor (PPI) use, travelling, the composition of the gut microbiota, and intake of certain foods. RESULTS: The most prevalent gene was TEM (53.0%), followed by SHV (18.4%), CTX-M (5.4%) and CMY (3.6%). Use of penicillins with extended spectrum was associated with TEM carriership, whereas use of macrolides and lincosamides was associated with TEM and SHV carriership. Interestingly, use of PPIs was associated with a higher prevalence of carriership of TEM, SHV and CMY (TEM: odds ratio [OR] 1.34; 95% confidence interval [CI] 1.05-1.77; SHV: OR 2.17; 95%CI 1.55-2.87; CMY: OR 2.26; 95%CI 1.23-4.11). Furthermore, associations were found between the richness and composition of the gut microbiota and TEM and SHV carriership. CONCLUSIONS: The prevalence of carriership of TEM was substantial, but the prevalence of carriership of the extended-spectrum ß-lactamase gene, CTX-M and the AmpC ß-lactamase gene, CMY was relatively low in this community-dwelling, population-based cohort. The composition of the microbiota might play a role in the retention of resistance genes, but future studies are necessary to further elucidate this relationship.

Age-dependent sex differences in calcium and phosphate homeostasis.

BACKGROUND: Sex differences in calcium and phosphate have been observed. We aimed to assess a relation with age. METHODS: We used the laboratory values of serum calcium, phosphate and albumin from three different samples ( 2005, 2010 and 2014 years) using the hospital information system of Erasmus MC, Rotterdam. The samples were divided into three age groups: 1-17, 18-44 and ≥45 years. Sex differences in calcium and phosphate were analyzed using ANCOVA, adjusting for age and serum albumin. Furthermore, sex by age interactions were determined and we analyzed differences between age groups stratified by sex. RESULTS: In all three samples there was a significant sex × age interaction for serum calcium and phosphate, whose levels were significantly higher in women compared to men above 45 years. No sex differences in the younger age groups were found. In men, serum calcium and phosphate levels were highest in the youngest age group compared to age groups of 18-44 and ≥45 years. In women, serum calcium levels were significantly higher in the age group 1-17 and the age group ≥45 years compared to the 18-44 years age group. In women, serum phosphate was different between the three different age groups with highest level in the group 1-17 years and lowest in the group 18-44 years. CONCLUSION: There are age- dependent sex differences in serum calcium and phosphate. Furthermore, we found differences in serum calcium and phosphate between different age groups. Underlying mechanisms for these age- and sex- differences are not yet fully elucidated.

Serum fatty acid chain length associates with prevalent symptomatic end-stage osteoarthritis, independent of BMI.

Higher body mass index (BMI) is associated with osteoarthritis (OA) in both weight-bearing and non-weight-bearing joints, suggesting a link between OA and poor metabolic health beyond mechanical loading. This risk may be influenced by systemic factors accompanying BMI. Fluctuations in concentrations of metabolites may mark or even contribute to development of OA. This study explores the association of metabolites with radiographic knee/hip OA prevalence and progression. A 1H-NMR-metabolomics assay was performed on plasma samples of 1564 cases for prevalent OA and 2,125 controls collected from the Rotterdam Study, CHECK, GARP/NORREF and LUMC-arthroplasty cohorts. OA prevalence and 5 to 10 year progression was assessed by means of Kellgren-Lawrence (KL) score and the OARSI-atlas. End-stage knee/hip OA (TJA) was defined as indication for arthroplasty surgery. Controls did not have OA at baseline or follow-up. Principal component analysis of 227 metabolites demonstrated 23 factors, of which 19 remained interpretable after quality-control. Associations of factor scores with OA definitions were investigated with logistic regression. Fatty acids chain length (FALen), which was included in two factors which associated with TJA, was individually associated with both overall OA as well as TJA. Increased Fatty Acid chain Length is associated with OA.

Long-term effects of antimicrobial drugs on the composition of the human gut microbiota.

INTRODUCTION: Antimicrobial drugs are known to have effects on the human gut microbiota. We studied the long-term temporal relationship between several antimicrobial drug groups and the composition of the human gut microbiota determined in feces samples. METHODS: Feces samples were obtained from a community-dwelling cohort of middle-aged and elderly individuals (Rotterdam Study). Bacterial DNA was isolated and sequenced using V3/V4 16 S ribosomal RNA sequencing (Illumina MiSeq). The time between the last prescription of several antimicrobial drug groups and the day of sampling was categorized into 0-12, 12-24, 24-48 and >48 months. The effects of the antimicrobial drug groups on the Shannon alpha-diversity (diversity), the Bray-Curtis beta-diversity (community structure), the Firmicutes/Bacteroidetes (F/B) ratio and individual genera were determined. RESULTS: We studied the gut microbiota of 1413 individuals (57.5% female, median age 62.6 years). The alpha-diversity was significantly lower up to 4 years after prescriptions of macrolides and lincosamides. It was also lower in the first year after the use of beta-lactams. The community structure (beta-diversity) of the microbiota was significantly different up to 4 years for macrolides and lincosamides, the first year for beta-lactams and at least the first year for quinolones. For the F/B ratio, drugs with a high anaerobic activity shifted the ratio toward Firmicutes in the first year whereas other antimicrobial drugs shifted the ratio toward Bacteroidetes. CONCLUSION: Use of antimicrobial drugs is associated with a shift in the composition of the gut microbiota.These effects differ in strength and duration, depending on the antimicrobial drug group used. These findings should be considered when prescribing antimicrobial drugs.

Genome-wide meta-analysis of variant-by-diuretic interactions as modulators of lipid traits in persons of European and African ancestry.

Hypertension (HTN) is a significant risk factor for cardiovascular morbidity and mortality. Metabolic abnormalities, including adverse cholesterol and triglycerides (TG) profiles, are frequent comorbid findings with HTN and contribute to cardiovascular disease. Diuretics, which are used to treat HTN and heart failure, have been associated with worsening of fasting lipid concentrations. Genome-wide meta-analyses with 39,710 European-ancestry (EA) individuals and 9925 African-ancestry (AA) individuals were performed to identify genetic variants that modify the effect of loop or thiazide diuretic use on blood lipid concentrations. Both longitudinal and cross sectional data were used to compute cohort-specific interaction results, which were then combined through meta-analysis in each ancestry. These ancestry-specific results were further combined through trans-ancestry meta-analysis. Analysis of EA data identified two genome-wide significant (p < 5 × 10-8) loci with single nucleotide variant (SNV)-loop diuretic interaction on TG concentrations (including COL11A1). Analysis of AA data identified one genome-wide significant locus adjacent to BMP2 with SNV-loop diuretic interaction on TG concentrations. Trans-ancestry analysis strengthened evidence of association for SNV-loop diuretic interaction at two loci (KIAA1217 and BAALC). There were few significant SNV-thiazide diuretic interaction associations on TG concentrations and for either diuretic on cholesterol concentrations. Several promising loci were identified that may implicate biologic pathways that contribute to adverse metabolic side effects from diuretic therapy.

The effect of antimicrobial drug use on the composition of the genitourinary microbiota in an elderly population.

BACKGROUND: The urinary tract is inhabited by a diversity of microorganisms, known as the genitourinary microbiota. Here, we investigated the association between the use of antimicrobial drugs and the composition of the genitourinary microbiota. RESULTS: Clean-catch urinary samples were collected from 27 participants of the Rotterdam Study. Bacterial DNA was extracted and the 16S ribosomal RNA gene variable regions V3 and V4 were analyzed using Illumina sequencing. 23 of the 27 participants were included in the analysis. The population consisted of 10 men and 13 women with a mean age of 75 ± 3 years. The time between the last prescription of an antimicrobial drug and sampling was determined and categorized. The use of antimicrobial drugs prior to urine sampling was associated with statistically significant differences in the beta-diversity of the genitourinary microbiota. No association was found between antimicrobial drug use and the alpha-diversity of the genitourinary microbiota. Operational Taxonomic Units (OTUs) that were lowest in participants who used antimicrobial drug belonged to Lactobacillus and Finegoldia. In contrast, an OTU belonging to the genus Parabacteroides had higher abundances. Also, an OTU belonging to the species E.coli was higher in the participants who used antimicrobial drugs. CONCLUSION: Prior use of antimicrobial drugs is associated with a different composition of the genitourinary microbiota. Our results might indicate a persisting effect of antimicrobial drugs on the composition of the microbiota, but reverse causality cannot be ruled out. Future studies are needed to differentiate between two possibilities. Genitourinary dysbiosis could be the result of antimicrobial drug use or genitourinary dysbiosis could be a risk factor for urinary tract infections resulting in increased use of antimicrobial drugs. This may have important implications for treatment and prevention of (recurrent) UTIs.

Non-genetic and genetic predictors of a superficial first basal cell carcinoma.

BACKGROUND: Several observational studies have suggested differences in the risk factor profile between patients with superficial basal cell carcinomas (BCCs) and non-superficial BCCs. OBJECTIVE: To test the reproducibility of previous study findings and to find new genetic and non-genetic predictors for patients with a superficial first BCC. METHODS: A total of 14.628 participants of northwestern European descent aged 45 years or older from a prospective population-based cohort study (Rotterdam Study) were linked with the Dutch Pathology Registry (PALGA) of whom 1528 were identified as BCC patients. After exclusion, 948 eligible BCC patients remained for further non-genetic analyses and 1014 for genetic analyses. We included 11 phenotypic, environmental and tumour-specific characteristics, and 20 candidate single nucleotide polymorphisms (SNP) as potential predictors for patients with a superficial first BCC. We performed binary logistic multivariable regression analyses. RESULTS: We found that patients with a superficial first BCC were significantly younger, almost two times more often female and 12-18 times more likely to have their BCC on the trunk or extremities than patients with a non-superficial first BCC. One SNP (rs12203592), mapped to IRF4, looked promising (OR 1.83, 95% CI 1.13-2.97, P-value <0.05), but after adjustment for multiple testing, no significant differences in genetic make-up between superficial BCC and non-superficial BCC patients were found. CONCLUSION: We conclude that patients with a superficial BCC differ from non-superficial BCC patients with respect to environmental factors (tumour localization as a proxy for UVR exposure) and phenotypic characteristics (age and sex), but we found no difference in genotype. As superficial BCC patients develop their first BCCs at a younger age, they could be at higher lifetime risk for subsequent skin cancers and therefore be an important group for secondary prevention.

The influence of genetic variation on late toxicities in childhood cancer survivors: A review.

INTRODUCTION: The variability in late toxicities among childhood cancer survivors (CCS) is only partially explained by treatment and baseline patient characteristics. Inter-individual variability in the association between treatment exposure and risk of late toxicity suggests that genetic variation possibly modifies this association. We reviewed the available literature on genetic susceptibility of late toxicity after childhood cancer treatment related to components of metabolic syndrome, bone mineral density, gonadal impairment and hearing impairment. METHODS: A systematic literature search was performed, using Embase, Cochrane Library, Google Scholar, MEDLINE, and Web of Science databases. Eligible publications included all English language reports of candidate gene studies and genome wide association studies (GWAS) that aimed to identify genetic risk factors associated with the four late toxicities, defined as toxicity present after end of treatment. RESULTS: Twenty-seven articles were identified, including 26 candidate gene studies: metabolic syndrome (n = 6); BMD (n = 6); gonadal impairment (n = 2); hearing impairment (n = 12) and one GWAS (metabolic syndrome). Eighty percent of the genetic studies on late toxicity after childhood cancer had relatively small sample sizes (n < 200), leading to insufficient power, and lacked adjustment for multiple comparisons. Only four (4/26 = 15%) candidate gene studies had their findings validated in independent replication cohorts as part of their own report. CONCLUSION: Genetic susceptibility associations are not consistent or not replicated and therefore, currently no evidence-based recommendations can be made for hearing impairment, gonadal impairment, bone mineral density impairment and metabolic syndrome in CCS. To advance knowledge related to genetic variation influencing late toxicities among CCS, future studies need adequate power, independent cohorts for replication, harmonization of disease outcomes and sample collections, and (international) collaboration.

Development of a prediction model for future risk of radiographic hip osteoarthritis.

OBJECTIVE: To develop and validate a prognostic model for incident radiologic hip osteoarthritis (HOA) and determine the value of previously identified predictive factors. DESIGN: We first validated previously reported predictive factors for HOA by performing univariate and multivariate analyses for all predictors in three large prospective cohorts (total sample size of 4548 with 653 incident cases). The prognostic model was developed in 2327 individuals followed for 10 years from the Rotterdam Study-I (RS-I) cohort. External validation of the model was tested on discrimination in two other cohorts: RS-II (n = 1435) and the Cohort Hip and Cohort Knee (CHECK) study (n = 786). RESULTS: From the total number of 28 previously reported predictive factors, we were able to replicate 13 factors, while 15 factors were not significantly predictive in a meta-analysis of the three cohorts. The basic model including the demographic, questionnaire, and clinical examination variables (area under the receiver-operating characteristic curve (AUC) = 0.67) or genetic markers (AUC = 0.55) or urinary C-terminal cross-linked telopeptide of type II collagen (uCTX-II) levels (AUC = 0.67) alone were poor predictors of HOA in all cohorts. Imaging factors showed the highest predictive value for the development of HOA (AUC = 0.74). Addition of imaging variables to the basic model led to substantial improvement in the discriminative ability of the model (AUC = 0.78) compared with uCTX-II (AUC = 0.74) or genetic markers (AUC = 0.68). Applying external validation, similar results were observed in the RS-II and the CHECK cohort. CONCLUSIONS: The developed prediction model included demographic, a limited number of questionnaire, and imaging risk factors seems promising for prediction of HOA.

Integration of Murine and Human Studies for Mapping Periodontitis Susceptibility.

Periodontitis is one of the most common inflammatory human diseases with a strong genetic component. Due to the limited sample size of available periodontitis cohorts and the underlying trait heterogeneity, genome-wide association studies (GWASs) of chronic periodontitis (CP) have largely been unsuccessful in identifying common susceptibility factors. A combination of quantitative trait loci (QTL) mapping in mice with association studies in humans has the potential to discover novel risk loci. To this end, we assessed alveolar bone loss in response to experimental periodontal infection in 25 lines (286 mice) from the Collaborative Cross (CC) mouse population using micro-computed tomography (µCT) analysis. The orthologous human chromosomal regions of the significant QTL were analyzed for association using imputed genotype data (OmniExpress BeadChip arrays) derived from case-control samples of aggressive periodontitis (AgP; 896 cases, 7,104 controls) and chronic periodontitis (CP; 2,746 cases, 1,864 controls) of northwest European and European American descent, respectively. In the mouse genome, QTL mapping revealed 2 significant loci (-log P = 5.3; false discovery rate = 0.06) on chromosomes 1 ( Perio3) and 14 ( Perio4). The mapping resolution ranged from ~1.5 to 3 Mb. Perio3 overlaps with a previously reported QTL associated with residual bone volume in F2 cross and includes the murine gene Ccdc121. Its human orthologue showed previously a nominal significant association with CP in humans. Use of variation data from the genomes of the CC founder strains further refined the QTL and suggested 7 candidate genes ( CAPN8, DUSP23, PCDH17, SNORA17, PCDH9, LECT1, and LECT2). We found no evidence of association of these candidates with the human orthologues. In conclusion, the CC populations enabled mapping of confined QTL that confer susceptibility to alveolar bone loss in mice and larger human phenotype-genotype samples and additional expression data from gingival tissues are likely required to identify true positive signals.

A DNA methylation biomarker of alcohol consumption.

The lack of reliable measures of alcohol intake is a major obstacle to the diagnosis and treatment of alcohol-related diseases. Epigenetic modifications such as DNA methylation may provide novel biomarkers of alcohol use. To examine this possibility, we performed an epigenome-wide association study of methylation of cytosine-phosphate-guanine dinucleotide (CpG) sites in relation to alcohol intake in 13 population-based cohorts (ntotal=13 317; 54% women; mean age across cohorts 42-76 years) using whole blood (9643 European and 2423 African ancestries) or monocyte-derived DNA (588 European, 263 African and 400 Hispanic ancestry) samples. We performed meta-analysis and variable selection in whole-blood samples of people of European ancestry (n=6926) and identified 144 CpGs that provided substantial discrimination (area under the curve=0.90-0.99) for current heavy alcohol intake (⩾42 g per day in men and ⩾28 g per day in women) in four replication cohorts. The ancestry-stratified meta-analysis in whole blood identified 328 (9643 European ancestry samples) and 165 (2423 African ancestry samples) alcohol-related CpGs at Bonferroni-adjusted P<1 × 10-7. Analysis of the monocyte-derived DNA (n=1251) identified 62 alcohol-related CpGs at P<1 × 10-7. In whole-blood samples of people of European ancestry, we detected differential methylation in two neurotransmitter receptor genes, the γ-Aminobutyric acid-A receptor delta and γ-aminobutyric acid B receptor subunit 1; their differential methylation was associated with expression levels of a number of genes involved in immune function. In conclusion, we have identified a robust alcohol-related DNA methylation signature and shown the potential utility of DNA methylation as a clinically useful diagnostic test to detect current heavy alcohol consumption.

A rare missense variant in RCL1 segregates with depression in extended families.

Depression is the most prevalent psychiatric disorder with a complex and elusive etiology that is moderately heritable. Identification of genes would greatly facilitate the elucidation of the biological mechanisms underlying depression, however, its complex etiology has proved to be a major bottleneck in the identification of its genetic risk factors, especially in genome-wide association-like studies. In this study, we exploit the properties of a genetic isolate and its family-based structure to explore whether relatively rare exonic variants influence the burden of depressive symptoms in families. Using a multistep approach involving linkage and haplotype analyses followed by exome sequencing in the Erasmus Rucphen Family (ERF) study, we identified a rare (minor allele frequency (MAF)=1%) missense c.1114C>T mutation (rs115482041) in the RCL1 gene segregating with depression across multiple generations. Rs115482041 showed significant association with depressive symptoms (N=2393, ßT-allele=2.33, P-value=1 × 10-4) and explained 2.9% of the estimated genetic variance of depressive symptoms (22%) in ERF. Despite being twice as rare (MAF<0.5%), c.1114C>T showed similar effect and significant association with depressive symptoms in samples from the independent population-based Rotterdam study (N=1604, ßT-allele=3.60, P-value=3 × 10-2). A comparison of RCL1 expression in human and mouse brain revealed a striking co-localization of RCL1 with the layer 1 interlaminar subclass of astrocytes found exclusively in higher-order primates. Our findings identify RCL1 as a novel candidate gene for depression and offer insights into mechanisms through which RCL1 may be relevant for depression.

Large-scale pharmacogenomic study of sulfonylureas and the QT, JT and QRS intervals: CHARGE Pharmacogenomics Working Group.

Sulfonylureas, a commonly used class of medication used to treat type 2 diabetes, have been associated with an increased risk of cardiovascular disease. Their effects on QT interval duration and related electrocardiographic phenotypes are potential mechanisms for this adverse effect. In 11 ethnically diverse cohorts that included 71 857 European, African-American and Hispanic/Latino ancestry individuals with repeated measures of medication use and electrocardiogram (ECG) measurements, we conducted a pharmacogenomic genome-wide association study of sulfonylurea use and three ECG phenotypes: QT, JT and QRS intervals. In ancestry-specific meta-analyses, eight novel pharmacogenomic loci met the threshold for genome-wide significance (P<5 × 10-8), and a pharmacokinetic variant in CYP2C9 (rs1057910) that has been associated with sulfonylurea-related treatment effects and other adverse drug reactions in previous studies was replicated. Additional research is needed to replicate the novel findings and to understand their biological basis.

Association between subclinical thyroid dysfunction and change in bone mineral density in prospective cohorts.

BACKGROUND: Subclinical hyperthyroidism (SHyper) has been associated with increased risk of hip and other fractures, but the linking mechanisms remain unclear. OBJECTIVE: To investigate the association between subclinical thyroid dysfunction and bone loss. METHODS: Individual participant data analysis was performed after a systematic literature search in MEDLINE/EMBASE (1946-2016). Two reviewers independently screened and selected prospective cohorts providing baseline thyroid status and serial bone mineral density (BMD) measurements. We classified thyroid status as euthyroidism (thyroid-stimulating hormone [TSH] 0.45-4.49 mIU/L), SHyper (TSH < 0.45 mIU/L) and subclinical hypothyroidism (SHypo, TSH ≥ 4.50-19.99 mIU/L) both with normal free thyroxine levels. Our primary outcome was annualized percentage BMD change (%ΔBMD) from serial dual X-ray absorptiometry scans of the femoral neck, total hip and lumbar spine, obtained from multivariable regression in a random-effects two-step approach. RESULTS: Amongst 5458 individuals (median age 72 years, 49.1% women) from six prospective cohorts, 451 (8.3%) had SHypo and 284 (5.2%) had SHyper. During 36 569 person-years of follow-up, those with SHyper had a greater annual bone loss at the femoral neck versus euthyroidism: %ΔBMD = -0.18 (95% CI: -0.34, -0.02; I2 = 0%), with a nonstatistically significant pattern at the total hip: %ΔBMD = -0.14 (95% CI: -0.38, 0.10; I2 = 53%), but not at the lumbar spine: %ΔBMD = 0.03 (95% CI: -0.30, 0.36; I2 = 25%); especially participants with TSH < 0.10 mIU/L showed an increased bone loss in the femoral neck (%Δ BMD = -0.59; [95% CI: -0.99, -0.19]) and total hip region (%ΔBMD = -0.46 [95% CI: -1.05, -0.13]). In contrast, SHypo was not associated with bone loss at any site. CONCLUSION: Amongst adults, SHyper was associated with increased femoral neck bone loss, potentially contributing to the increased fracture risk.