Search

With aging skeletal muscle fibers undergo repeating cycles of denervation and reinnervation. In approximately the 8th decade of life reinnervation no longer keeps pace, resulting in the accumulation of persistently denervated muscle fibers that in turn cause an acceleration of muscle dysfunction. The significance of denervation in important clinical outcomes with aging is poorly studied. The Study of Muscle, Mobility, and Aging (SOMMA) is a large cohort study with the primary objective to assess how aging muscle biology impacts clinically important traits. Using transcriptomics data from vastus lateralis muscle biopsies in 575 participants we have selected 49 denervation-responsive genes to provide insights to the burden of denervation in SOMMA, to test the hypothesis that greater expression of denervation-responsive genes negatively associates with SOMMA participant traits that included time to walk 400 meters, fitness (VO2peak), maximal mitochondrial respiration, muscle mass and volume, and leg muscle strength and power. Consistent with our hypothesis, increased transcript levels of: a calciumdependent intercellular adhesion glycoprotein (CDH15), acetylcholine receptor subunits (CHRNA1, CHRND, CHRNE), a glycoprotein promoting reinnervation (NCAM1), a transcription factor regulating aspects of muscle organization (RUNX1), and a sodium channel (SCN5A) were each negatively associated with at least 3 of these traits. VO2peak and maximal respiration had the strongest negative associations with 15 and 19 denervation-responsive genes, respectively. In conclusion, the abundance of denervationresponsive gene transcripts is a significant determinant of muscle and mobility outcomes in aging humans, supporting the imperative to identify new treatment strategies to restore innervation in advanced age.

BackgroundHip minimum joint space width (mJSW) provides a proxy for cartilage thickness. This study aimed to conduct a genome-wide association study (GWAS) of mJSW to (i) identify new genetic determinants of mJSW and (ii) identify which mJSW loci convey hip osteoarthritis (HOA) risk and would therefore be of therapeutic interest.MethodsGWAS meta-analysis of hip mJSW derived from plain X-rays and DXA was performed, stratified by sex and adjusted for age and ancestry principal components. Mendelian randomisation (MR) and cluster analyses were used to examine causal effect of mJSW on HOA.Findings50,745 individuals were included in the meta-analysis. 42 SNPs, which mapped to 39 loci, were identified. Mendelian randomisation (MR) revealed little evidence of a causal effect of mJSW on HOA (ORIVW 0.98 [95% CI 0.82-1.18]). However, MR-Clust analysis suggested the null MR estimates reflected the net effect of two distinct causal mechanisms cancelling each other out, one of which was protective, whereas the other increased HOA susceptibility. For the latter mechanism, all loci were positively associated with height, suggesting mechanisms leading to greater height and mJSW increase the risk of HOA in later life.InterpretationsOne group of mJSW loci reduce HOA risk via increased mJSW, suggesting possible utility as targets for chondroprotective therapies. The second group of mJSW loci increased HOA risk, despite increasing mJSW, but were also positively related to height, suggesting they contribute to HOA risk via a growth-related mechanism.FundingPrimarily funded by the Medical Research Council and Wellcome Trust.

Study objectivesPeriodic limb movement in sleep is a common sleep phenotype characterized by repetitive leg movements that occur during or before sleep. We conducted a genome-wide association study (GWAS) of periodic limb movements in sleep (PLMS) using a joint analysis (i.e., discovery, replication, and joint meta-analysis) of four cohorts (MrOS, the Wisconsin Sleep Cohort Study, HypnoLaus, and MESA), comprised of 6843 total subjects.MethodsThe MrOS study and Wisconsin Sleep Cohort Study (N = 1745 cases) were used for discovery. Replication in the HypnoLaus and MESA cohorts (1002 cases) preceded joint meta-analysis. We also performed LD score regression, estimated heritability, and computed genetic correlations between potentially associated traits such as restless leg syndrome (RLS) and insomnia. The causality and direction of the relationships between PLMS and RLS was evaluated using Mendelian randomization.ResultsWe found 2 independent loci were significantly associated with PLMS: rs113851554 (p = 3.51 × 10-12, β = 0.486), an SNP located in a putative regulatory element of intron eight of MEIS1 (2p14); and rs9369062 (p = 3.06 × 10-22, β = 0.2093), a SNP located in the intron region of BTBD9 (6p12); both of which were also lead signals in RLS GWAS. PLMS is genetically correlated with insomnia, risk of stroke, and RLS, but not with iron deficiency. Pleiotropy adjusted Mendelian randomization analysis identified a causal effect of RLS on PLMS.ConclusionsBecause PLMS is more common than RLS, PLMS may have multiple causes and additional studies are needed to further validate these findings.

Genetic association studies of many heritable traits resulting from physiological testing often have modest sample sizes due to the cost and burden of the required phenotyping. This reduces statistical power and limits discovery of multiple genetic associations. We present a strategy to leverage pleiotropy between traits to both discover new loci and to provide mechanistic hypotheses of the underlying pathophysiology. Specifically, we combine a colocalization test with a locus-level test of pleiotropy. In simulations, we show that this approach is highly selective for identifying true pleiotropy driven by the same causative variant, thereby improves the chance to replicate the associations in underpowered validation cohorts and leads to higher interpretability. Here, as an exemplar, we use Obstructive Sleep Apnea (OSA), a common disorder diagnosed using overnight multi-channel physiological testing. We leverage pleiotropy with relevant cellular and cardio-metabolic phenotypes and gene expression traits to map new risk loci in an underpowered OSA GWAS. We identify several pleiotropic loci harboring suggestive associations to OSA and genome-wide significant associations to other traits, and show that their OSA association replicates in independent cohorts of diverse ancestries. By investigating pleiotropic loci, our strategy allows proposing new hypotheses about OSA pathobiology across many physiological layers. For example, we identify and replicate the pleiotropy across the plateletcrit, OSA and an eQTL of DNA primase subunit 1 (PRIM1) in immune cells. We find suggestive links between OSA, a measure of lung function (FEV1/FVC), and an eQTL of matrix metallopeptidase 15 (MMP15) in lung tissue. We also link a previously known genome-wide significant peak for OSA in the hexokinase 1 (HK1) locus to hematocrit and other red blood cell related traits. Thus, the analysis of pleiotropic associations has the potential to assemble diverse phenotypes into a chain of mechanistic hypotheses that provide insight into the pathogenesis of complex human diseases.

Rationale: Obstructive sleep apnea (OSA) is a common disorder associated with increased risk for cardiovascular disease, diabetes, and premature mortality. There is strong clinical and epidemiologic evidence supporting the importance of genetic factors influencing OSA but limited data implicating specific genes. Objectives: To search for rare variants contributing to OSA severity. Methods: Leveraging high-depth genomic sequencing data from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program and imputed genotype data from multiple population-based studies, we performed linkage analysis in the CFS (Cleveland Family Study), followed by multistage gene-based association analyses in independent cohorts for apnea-hypopnea index (AHI) in a total of 7,708 individuals of European ancestry. Measurements and Main Results: Linkage analysis in the CFS identified a suggestive linkage peak on chromosome 7q31 (LOD = 2.31). Gene-based analysis identified 21 noncoding rare variants in CAV1 (Caveolin-1) associated with lower AHI after accounting for multiple comparisons (P = 7.4 × 10-8). These noncoding variants together significantly contributed to the linkage evidence (P

Chronic increases in pro-inflammatory cytokines in older adults, known as inflammaging, are an important risk factor for morbidity and mortality in the aging population. It has been suggested that circadian disruption may play a role in chronic inflammation, but there has been limited study that investigated the overall profile of 24-hour rest-activity rhythms in relation to inflammation using longitudinal data. In the Outcomes of Sleep Disorders in Older Men Study, we applied the extended cosine model to derive multiple rest-activity rhythm characteristics using multiday actigraphy, and examined their associations with 6 inflammatory markers (ie, C-reactive protein [CRP], interleukin 6 [IL-6], tumor necrosis factor alpha [TNF-α], tumor necrosis factor alpha soluble receptor II [TNF-α-sRII], interleukin-1β [IL-1β], interferon gamma [IFN-γ]) measured from fasting blood. We assessed both the cross-sectional association between rest-activity rhythms and inflammatory markers measured at baseline, and the prospective association between baseline rest-activity rhythms and changes in inflammatory markers over 3.5 years of follow-up. We found that multiple rest-activity characteristics, including lower amplitude and relative amplitude, and decreased overall rhythmicity, were associated with higher levels of CRP, IL-6, TNF-α, and TNF-α-sRII, but not IL-1β and IFN-γ at baseline. Moreover, the lowest quartile of these 3 rest-activity characteristics was associated with an approximately 2-fold increase in the odds of having elevated inflammation (ie, having 3 or more markers in the highest quartile) at baseline. However, we found little evidence supporting a relationship between rest-activity rhythm characteristics and changes in inflammatory markers. Future studies should clarify the dynamic relationship between rest-activity rhythms and inflammation in different populations, and evaluate the effects of improving rest-activity profiles on inflammation and related disease outcomes.

BackgroundSleep-disordered breathing is a common disorder associated with significant morbidity. The genetic architecture of sleep-disordered breathing remains poorly understood. Through the NHLBI Trans-Omics for Precision Medicine (TOPMed) program, we performed the first whole-genome sequence analysis of sleep-disordered breathing.MethodsThe study sample was comprised of 7988 individuals of diverse ancestry. Common-variant and pathway analyses included an additional 13,257 individuals. We examined five complementary traits describing different aspects of sleep-disordered breathing: the apnea-hypopnea index, average oxyhemoglobin desaturation per event, average and minimum oxyhemoglobin saturation across the sleep episode, and the percentage of sleep with oxyhemoglobin saturation

Chronic increases in pro-inflammatory cytokines in older adults, known as inflammaging, is an important risk factor for morbidity and mortality in the aging population. It has been suggested that circadian disruption may play a role in chronic inflammation, but there has been limited study that investigated the overall profile of 24-hour rest-activity rhythms in relation to inflammation using longitudinal data. In the Outcomes of Sleep Disorders in Older Men Study, we applied the extended cosine model to derive multiple rest-activity rhythm characteristics using multi-day actigraphy, and examined their associations with six inflammatory markers (i.e., CRP, IL-6, TNF-α, TNF-α-sRII, IL-1 β, IFN-γ) measured from fasting blood. We assessed both the cross-sectional association between rest-activity rhythms and inflammatory markers measured at baseline, and the prospective association between baseline rest-activity rhythms and changes in in inflammatory markers over 3.5 years of follow up. We found that multiple rest-activity characteristics, including lower amplitude and relative amplitude, and decreased overall rhythmicity, were associated with higher levels of CRP, IL-6, TNF-α, and TNF-α-sRII, but not IL-1β and IFN-γ at baseline. Moreover, the lowest quartile of these three rest-activity characteristics was associated with an approximately two-fold increase in the odds of having elevated inflammation (i.e. having three or more markers in the highest quartile) at baseline. However, we found little evidence supporting a relationship between rest-activity rhythm characteristics and changes in inflammatory markers. Future studies should clarify the dynamic relationship between rest-activity rhythms and inflammation in different populations, and evaluate the effects of improving rest-activity profiles on inflammation and related disease outcomes.

signatureSearch is an R/Bioconductor package that integrates a suite of existing and novel algorithms into an analysis environment for gene expression signature (GES) searching combined with functional enrichment analysis (FEA) and visualization methods to facilitate the interpretation of the search results. In a typical GES search (GESS), a query GES is searched against a database of GESs obtained from large numbers of measurements, such as different genetic backgrounds, disease states and drug perturbations. Database matches sharing correlated signatures with the query indicate related cellular responses frequently governed by connected mechanisms, such as drugs mimicking the expression responses of a disease. To identify which processes are predominantly modulated in the GESS results, we developed specialized FEA methods combined with drug-target network visualization tools. The provided analysis tools are useful for studying the effects of genetic, chemical and environmental perturbations on biological systems, as well as searching single cell GES databases to identify novel network connections or cell types. The signatureSearch software is unique in that it provides access to an integrated environment for GESS/FEA routines that includes several novel search and enrichment methods, efficient data structures, and access to pre-built GES databases, and allowing users to work with custom databases.

ObjectiveDisruption of rest-activity rhythms is cross-sectionally associated with metabolic disorders, including type 2 diabetes, yet it remains unclear whether it predicts impaired glucose metabolism and homeostasis. The aim of this study is to examine the cross-sectional and prospective associations between rest-activity rhythm characteristics and glycemic measures in a cohort of older men.Research design and methodsBaseline rest-activity rhythms were derived from actigraphy with use of extended cosine model analysis. With subjects fasting, glucose, insulin, and HOMA of insulin resistance (HOMA-IR) were measured from blood at baseline and after ∼3.5 years. Type 2 diabetes was defined based on self-report, medication use, and fasting glucose.ResultsIn the cross-sectional analysis (n = 2,450), lower 24-h amplitude-to-mesor ratio (i.e., mean activity-adjusted rhythm amplitude) and reduced overall rhythmicity were associated with higher fasting insulin and HOMA-IR (all P trend < 0.0001), indicating increased insulin resistance. The odds of baseline type 2 diabetes were significantly higher among those in the lowest quartile of amplitude (Q1) (odds ratio [OR]Q1 vs. Q4 1.63 [95% CI 1.14, 2.30]) and late acrophase group (ORlate vs. normal 1.46 [95% CI 1.04, 2.04]). In the prospective analysis (n = 861), multiple rest-activity characteristics predicted a two- to threefold increase in type 2 diabetes risk, including a lower amplitude (ORQ1 vs. Q4 3.81 [95% CI 1.45, 10.00]) and amplitude-to-mesor ratio (OR 2.79 [95% CI 1.10, 7.07]), reduced overall rhythmicity (OR 3.49 [95% CI 1.34, 9.10]), and a late acrophase (OR 2.44 [1.09, 5.47]).ConclusionsRest-activity rhythm characteristics are associated with impaired glycemic metabolism and homeostasis and higher risk of incident type 2 diabetes.

Abstract Sleep disturbances and physical inactivity have been associated with chronic inflammation, an important risk factor for cognitive decline in the aging population. However most previous studies focused on the cross-sectional relationships between sleep and physical activity and inflammation. In the Outcomes of Sleep Disorders in Older Men (MrOS Sleep) study, we studied both the cross-sectional and prospective associations between characteristics of 24-hour rest-activity rhythms measured by actigraphy and inflammation index measured by multiple circulating markers. In cross-sectional analysis, a lower amplitude is associated with elevated inflammation (Odds ratio Q4 vs Q1 (95% Confidence interval): 1.65 (1.22, 2.24)). In prospective analysis, an earlier acrophase (

Human longevity is heritable, but genome-wide association (GWA) studies have had limited success. Here, we perform two meta-analyses of GWA studies of a rigorous longevity phenotype definition including 11,262/3484 cases surviving at or beyond the age corresponding to the 90th/99th survival percentile, respectively, and 25,483 controls whose age at death or at last contact was at or below the age corresponding to the 60th survival percentile. Consistent with previous reports, rs429358 (apolipoprotein E (ApoE) ε4) is associated with lower odds of surviving to the 90th and 99th percentile age, while rs7412 (ApoE ε2) shows the opposite. Moreover, rs7676745, located near GPR78, associates with lower odds of surviving to the 90th percentile age. Gene-level association analysis reveals a role for tissue-specific expression of multiple genes in longevity. Finally, genetic correlation of the longevity GWA results with that of several disease-related phenotypes points to a shared genetic architecture between health and longevity.

Obstructive sleep apnea (OSA) is a common disorder associated with increased risk of cardiovascular disease and mortality. Its prevalence and severity vary across ancestral background. Although OSA traits are heritable, few genetic associations have been identified. To identify genetic regions associated with OSA and improve statistical power, we applied admixture mapping on three primary OSA traits [the apnea hypopnea index (AHI), overnight average oxyhemoglobin saturation (SaO2) and percentage time SaO2

We aimed to report the first genomewide association study (GWAS) meta-analysis of dual-energy X-ray absorptiometry (DXA)-derived hip shape, which is thought to be related to the risk of both hip osteoarthritis and hip fracture. Ten hip shape modes (HSMs) were derived by statistical shape modeling using SHAPE software, from hip DXA scans in the Avon Longitudinal Study of Parents and Children (ALSPAC; adult females), TwinsUK (mixed sex), Framingham Osteoporosis Study (FOS; mixed), Osteoporotic Fractures in Men study (MrOS), and Study of Osteoporotic Fractures (SOF; females) (total N = 15,934). Associations were adjusted for age, sex, and ancestry. Five genomewide significant (p  0.5) were identified, which intersected with open chromatin sites as detected by ATAC-seq performed on embryonic mouse proximal femora. In conclusion, we identified eight SNPs independently associated with hip shape, most of which were associated with height and/or mapped close to endochondral bone formation genes, consistent with a contribution of processes involved in limb growth to hip shape and pathological sequelae. These findings raise the possibility that genetic studies of hip shape might help in understanding potential pathways involved in hip osteoarthritis and hip fracture. © 2018 The Authors. Journal of Bone and Mineral Research Published by Wiley Periodicals, Inc.

BackgroundLean body mass (LM) plays an important role in mobility and metabolic function. We previously identified five loci associated with LM adjusted for fat mass in kilograms. Such an adjustment may reduce the power to identify genetic signals having an association with both lean mass and fat mass.ObjectivesTo determine the impact of different fat mass adjustments on genetic architecture of LM and identify additional LM loci.MethodsWe performed genome-wide association analyses for whole-body LM (20 cohorts of European ancestry with n = 38,292) measured using dual-energy X-ray absorptiometry) or bioelectrical impedance analysis, adjusted for sex, age, age2, and height with or without fat mass adjustments (Model 1 no fat adjustment; Model 2 adjustment for fat mass as a percentage of body mass; Model 3 adjustment for fat mass in kilograms).ResultsSeven single-nucleotide polymorphisms (SNPs) in separate loci, including one novel LM locus (TNRC6B), were successfully replicated in an additional 47,227 individuals from 29 cohorts. Based on the strengths of the associations in Model 1 vs Model 3, we divided the LM loci into those with an effect on both lean mass and fat mass in the same direction and refer to those as "sumo wrestler" loci (FTO and MC4R). In contrast, loci with an impact specifically on LM were termed "body builder" loci (VCAN and ADAMTSL3). Using existing available genome-wide association study databases, LM increasing alleles of SNPs in sumo wrestler loci were associated with an adverse metabolic profile, whereas LM increasing alleles of SNPs in "body builder" loci were associated with metabolic protection.ConclusionsIn conclusion, we identified one novel LM locus (TNRC6B). Our results suggest that a genetically determined increase in lean mass might exert either harmful or protective effects on metabolic traits, depending on its relation to fat mass.

BackgroundThe electrocardiographically quantified QRS duration measures ventricular depolarization and conduction. QRS prolongation has been associated with poor heart failure prognosis and cardiovascular mortality, including sudden death. While previous genome-wide association studies (GWAS) have identified 32 QRS SNPs across 26 loci among European, African, and Asian-descent populations, the genetics of QRS among Hispanics/Latinos has not been previously explored.MethodsWe performed a GWAS of QRS duration among Hispanic/Latino ancestry populations (n = 15,124) from four studies using 1000 Genomes imputed genotype data (adjusted for age, sex, global ancestry, clinical and study-specific covariates). Study-specific results were combined using fixed-effects, inverse variance-weighted meta-analysis.ResultsWe identified six loci associated with QRS (P

Sleep disordered breathing (SDB)-related overnight hypoxemia is associated with cardiometabolic disease and other comorbidities. Understanding the genetic bases for variations in nocturnal hypoxemia may help understand mechanisms influencing oxygenation and SDB-related mortality. We conducted genome-wide association tests across 10 cohorts and 4 populations to identify genetic variants associated with three correlated measures of overnight oxyhemoglobin saturation: average and minimum oxyhemoglobin saturation during sleep and the percent of sleep with oxyhemoglobin saturation under 90%. The discovery sample consisted of 8,326 individuals. Variants with p < 1 × 10(-6) were analyzed in a replication group of 14,410 individuals. We identified 3 significantly associated regions, including 2 regions in multi-ethnic analyses (2q12, 10q22). SNPs in the 2q12 region associated with minimum SpO2 (rs78136548 p = 2.70 × 10(-10)). SNPs at 10q22 were associated with all three traits including average SpO2 (rs72805692 p = 4.58 × 10(-8)). SNPs in both regions were associated in over 20,000 individuals and are supported by prior associations or functional evidence. Four additional significant regions were detected in secondary sex-stratified and combined discovery and replication analyses, including a region overlapping Reelin, a known marker of respiratory complex neurons.These are the first genome-wide significant findings reported for oxyhemoglobin saturation during sleep, a phenotype of high clinical interest. Our replicated associations with HK1 and IL18R1 suggest that variants in inflammatory pathways, such as the biologically-plausible NLRP3 inflammasome, may contribute to nocturnal hypoxemia.

Back pain is the #1 cause of years lived with disability worldwide, yet surprisingly little is known regarding the biology underlying this symptom. We conducted a genome-wide association study (GWAS) meta-analysis of chronic back pain (CBP). Adults of European ancestry were included from 15 cohorts in the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium, and from the UK Biobank interim data release. CBP cases were defined as those reporting back pain present for ≥3-6 months; non-cases were included as comparisons ("controls"). Each cohort conducted genotyping using commercially available arrays followed by imputation. GWAS used logistic regression models with additive genetic effects, adjusting for age, sex, study-specific covariates, and population substructure. The threshold for genome-wide significance in the fixed-effect inverse-variance weighted meta-analysis was p

BACKGROUND/RATIONALE:Experimental studies support the role of neurotransmitter genes in dementia risk, but human studies utilizing single variants in candidate genes have had limited success. METHODS:We used the gene-based testing program Versatile Gene-based Association Study to assess whether aggregate variation across 6 neurotransmitter pathways influences risk of cognitive decline in 8159 cognitively normal elderly (≥65 years old) adults from 3 community-based cohorts. RESULTS:Common genetic variation in GNG4 and KCNQ2 was associated with cognitive decline. In human brain tissue data sets, both GNG4 and KCNQ2 show higher expression in hippocampus relative to other brain regions; GNG4 expression decreases with advancing age. Both GNG4 and KCNQ2 show highest expression in fetal astrocytes. CONCLUSION:Genetic variation analyses and gene expression data suggest that GNG4 and KCNQ2 may be associated with cognitive decline in normal aging. Gene-based testing of neurotransmitter pathways may confirm and reveal novel risk genes in future studies of healthy cognitive aging.

Obstructive sleep apnea (OSA) is a common heritable disorder displaying marked sexual dimorphism in disease prevalence and progression. Previous genetic association studies have identified a few genetic loci associated with OSA and related quantitative traits, but they have only focused on single ethnic groups, and a large proportion of the heritability remains unexplained. The apnea-hypopnea index (AHI) is a commonly used quantitative measure characterizing OSA severity. Because OSA differs by sex, and the pathophysiology of obstructive events differ in rapid eye movement (REM) and non-REM (NREM) sleep, we hypothesized that additional genetic association signals would be identified by analyzing the NREM/REM-specific AHI and by conducting sex-specific analyses in multiethnic samples. We performed genome-wide association tests for up to 19,733 participants of African, Asian, European, and Hispanic/Latino American ancestry in 7 studies. We identified rs12936587 on chromosome 17 as a possible quantitative trait locus for NREM AHI in men (N = 6,737; P = 1.7 × 10-8) but not in women (P = 0.77). The association with NREM AHI was replicated in a physiological research study (N = 67; P = 0.047). This locus overlapping the RAI1 gene and encompassing genes PEMT1, SREBF1, and RASD1 was previously reported to be associated with coronary artery disease, lipid metabolism, and implicated in Potocki-Lupski syndrome and Smith-Magenis syndrome, which are characterized by abnormal sleep phenotypes. We also identified gene-by-sex interactions in suggestive association regions, suggesting that genetic variants for AHI appear to vary by sex, consistent with the clinical observations of strong sexual dimorphism.

Bone mineral density (BMD) assessed by DXA is used to evaluate bone health. In children, total body (TB) measurements are commonly used; in older individuals, BMD at the lumbar spine (LS) and femoral neck (FN) is used to diagnose osteoporosis. To date, genetic variants in more than 60 loci have been identified as associated with BMD. To investigate the genetic determinants of TB-BMD variation along the life course and test for age-specific effects, we performed a meta-analysis of 30 genome-wide association studies (GWASs) of TB-BMD including 66,628 individuals overall and divided across five age strata, each spanning 15 years. We identified variants associated with TB-BMD at 80 loci, of which 36 have not been previously identified; overall, they explain approximately 10% of the TB-BMD variance when combining all age groups and influence the risk of fracture. Pathway and enrichment analysis of the association signals showed clustering within gene sets implicated in the regulation of cell growth and SMAD proteins, overexpressed in the musculoskeletal system, and enriched in enhancer and promoter regions. These findings reveal TB-BMD as a relevant trait for genetic studies of osteoporosis, enabling the identification of variants and pathways influencing different bone compartments. Only variants in ESR1 and close proximity to RANKL showed a clear effect dependency on age. This most likely indicates that the majority of genetic variants identified influence BMD early in life and that their effect can be captured throughout the life course.

Electrocardiographic PR interval measures atrio-ventricular depolarization and conduction, and abnormal PR interval is a risk factor for atrial fibrillation and heart block. Our genome-wide association study of over 92,000 European-descent individuals identifies 44 PR interval loci (34 novel). Examination of these loci reveals known and previously not-yet-reported biological processes involved in cardiac atrial electrical activity. Genes in these loci are over-represented in cardiac disease processes including heart block and atrial fibrillation. Variants in over half of the 44 loci were associated with atrial or blood transcript expression levels, or were in high linkage disequilibrium with missense variants. Six additional loci were identified either by meta-analysis of ~105,000 African and European-descent individuals and/or by pleiotropic analyses combining PR interval with heart rate, QRS interval, and atrial fibrillation. These findings implicate developmental pathways, and identify transcription factors, ion-channel genes, and cell-junction/cell-signaling proteins in atrio-ventricular conduction, identifying potential targets for drug development.

A correction to this article has been published and is linked from the HTML version of this article.

Poor bone quality contributes to bone fragility in diabetes, aging, and osteogenesis imperfecta. However, the mechanisms controlling bone quality are not well understood, contributing to the current lack of strategies to diagnose or treat bone quality deficits. Transforming growth factor beta (TGF-β) signaling is a crucial mechanism known to regulate the material quality of bone, but its cellular target in this regulation is unknown. Studies showing that osteocytes directly remodel their perilacunar/canalicular matrix led us to hypothesize that TGF-β controls bone quality through perilacunar/canalicular remodeling (PLR). Using inhibitors and mice with an osteocyte-intrinsic defect in TGF-β signaling (TβRIIocy-/-), we show that TGF-β regulates PLR in a cell-intrinsic manner to control bone quality. Altogether, this study emphasizes that osteocytes are key in executing the biological control of bone quality through PLR, thereby highlighting the fundamental role of osteocyte-mediated PLR in bone homeostasis and fragility.

Lean body mass, consisting mostly of skeletal muscle, is important for healthy aging. We performed a genome-wide association study for whole body (20 cohorts of European ancestry with n = 38,292) and appendicular (arms and legs) lean body mass (n = 28,330) measured using dual energy X-ray absorptiometry or bioelectrical impedance analysis, adjusted for sex, age, height, and fat mass. Twenty-one single-nucleotide polymorphisms were significantly associated with lean body mass either genome wide (p

Type 2 diabetes (T2D) affects Latinos at twice the rate seen in populations of European descent. We recently identified a risk haplotype spanning SLC16A11 that explains ∼20% of the increased T2D prevalence in Mexico. Here, through genetic fine-mapping, we define a set of tightly linked variants likely to contain the causal allele(s). We show that variants on the T2D-associated haplotype have two distinct effects: (1) decreasing SLC16A11 expression in liver and (2) disrupting a key interaction with basigin, thereby reducing cell-surface localization. Both independent mechanisms reduce SLC16A11 function and suggest SLC16A11 is the causal gene at this locus. To gain insight into how SLC16A11 disruption impacts T2D risk, we demonstrate that SLC16A11 is a proton-coupled monocarboxylate transporter and that genetic perturbation of SLC16A11 induces changes in fatty acid and lipid metabolism that are associated with increased T2D risk. Our findings suggest that increasing SLC16A11 function could be therapeutically beneficial for T2D. VIDEO ABSTRACT.

BackgroundIncreased heart rate and a prolonged QT interval are important risk factors for cardiovascular morbidity and mortality, and can be influenced by the use of various medications, including tricyclic/tetracyclic antidepressants (TCAs). We aim to identify genetic loci that modify the association between TCA use and RR and QT intervals.Methods and resultsWe conducted race/ethnic-specific genome-wide interaction analyses (with HapMap phase II imputed reference panel imputation) of TCAs and resting RR and QT intervals in cohorts of European (n=45 706; n=1417 TCA users), African (n=10 235; n=296 TCA users) and Hispanic/Latino (n=13 808; n=147 TCA users) ancestry, adjusted for clinical covariates. Among the populations of European ancestry, two genome-wide significant loci were identified for RR interval: rs6737205 in BRE (β=56.3, pinteraction=3.9e-9) and rs9830388 in UBE2E2 (β=25.2, pinteraction=1.7e-8). In Hispanic/Latino cohorts, rs2291477 in TGFBR3 significantly modified the association between TCAs and QT intervals (β=9.3, pinteraction=2.55e-8). In the meta-analyses of the other ethnicities, these loci either were excluded from the meta-analyses (as part of quality control), or their effects did not reach the level of nominal statistical significance (pinteraction>0.05). No new variants were identified in these ethnicities. No additional loci were identified after inverse-variance-weighted meta-analysis of the three ancestries.ConclusionsAmong Europeans, TCA interactions with variants in BRE and UBE2E2 were identified in relation to RR intervals. Among Hispanic/Latinos, variants in TGFBR3 modified the relation between TCAs and QT intervals. Future studies are required to confirm our results.

This study presents an analysis of the small molecule bioactivity profiles across large quantities of diverse protein families represented in PubChem BioAssay. We compared the bioactivity profiles of FDA approved drugs to non-FDA approved compounds, and report several distinct patterns characteristic of the approved drugs. We found that a large fraction of the previously reported higher target promiscuity among FDA approved compounds, compared to non-FDA approved bioactives, was frequently due to cross-reactivity within rather than across protein families. We identified 804 potentially novel protein target candidates for FDA approved drugs, as well as 901 potentially novel target candidates with active non-FDA approved compounds, but no FDA approved drugs with activity against these targets. We also identified 486348 potentially novel compounds active against the same targets as FDA approved drugs, as well as 153402 potentially novel compounds active against targets without active FDA approved drugs. By quantifying the agreement among replicated screens, we estimated that more than half of these novel outcomes are reproducible. Using biclustering, we identified many dense clusters of FDA approved drugs with enriched activity against a common set of protein targets. We also report the distribution of compound promiscuity using a Bayesian statistical model, and report the sensitivity and specificity of two common methods for identifying promiscuous compounds. Aggregator assays exhibited greater accuracy in identifying highly promiscuous compounds, while PAINS substructures were able to identify a much larger set of "middle range" promiscuous compounds. Additionally, we report a large number of promiscuous compounds not identified as aggregators or PAINS. In summary, the results of this study represent a rich reference for selecting novel drug and target protein candidates, as well as for eliminating candidate compounds with unselective activities.

BackgroundIn addition to lowering low density lipoprotein cholesterol (LDL-C), statin therapy also raises high density lipoprotein cholesterol (HDL-C) levels. Inter-individual variation in HDL-C response to statins may be partially explained by genetic variation.Methods and resultsWe performed a meta-analysis of genome-wide association studies (GWAS) to identify variants with an effect on statin-induced high density lipoprotein cholesterol (HDL-C) changes. The 123 most promising signals with p

BackgroundOsteoporosis is a major public health problem associated with excess disability and mortality. It is estimated that 50-70% of the variation in osteoporotic fracture risk is attributable to genetic factors. The purpose of this hypothesis-generating study was to identify possible genetic determinants of fracture among African American (AA) women in a GWAS meta-analysis.MethodsData on clinical fractures (all fractures except fingers, toes, face, skull or sternum) were analyzed among AA female participants in the Women's Health Initiative (WHI) (N = 8155), Cardiovascular Health Study (CHS) (N = 504), BioVU (N = 704), Health ABC (N = 651), and the Johnston County Osteoarthritis Project (JoCoOA) (N = 291). Affymetrix (WHI) and Illumina (Health ABC, JoCoOA, BioVU, CHS) GWAS panels were used for genotyping, and a 1:1 ratio of YRI:CEU HapMap haplotypes was used as an imputation reference panel. We used Cox proportional hazard models or logistic regression to evaluate the association of ~ 2.5 million SNPs with fracture risk, adjusting for ancestry, age, and geographic region where applicable. We conducted a fixed-effects, inverse variance-weighted meta-analysis. Genome-wide significance was set at P

The electrocardiographic QRS duration, a measure of ventricular depolarization and conduction, is associated with cardiovascular mortality. While single nucleotide polymorphisms (SNPs) associated with QRS duration have been identified at 22 loci in populations of European descent, the genetic architecture of QRS duration in non-European populations is largely unknown. We therefore performed a genome-wide association study (GWAS) meta-analysis of QRS duration in 13,031 African Americans from ten cohorts and a transethnic GWAS meta-analysis with additional results from populations of European descent. In the African American GWAS, a single genome-wide significant SNP association was identified (rs3922844, P = 4 × 10-14) in intron 16 of SCN5A, a voltage-gated cardiac sodium channel gene. The QRS-prolonging rs3922844 C allele was also associated with decreased SCN5A RNA expression in human atrial tissue (P = 1.1 × 10-4). High density genotyping revealed that the SCN5A association region in African Americans was confined to intron 16. Transethnic GWAS meta-analysis identified novel SNP associations on chromosome 18 in MYL12A (rs1662342, P = 4.9 × 10-8) and chromosome 1 near CD1E and SPTA1 (rs7547997, P = 7.9 × 10-9). The 22 QRS loci previously identified in populations of European descent were enriched for significant SNP associations with QRS duration in African Americans (P = 9.9 × 10-7), and index SNP associations in or near SCN5A, SCN10A, CDKN1A, NFIA, HAND1, TBX5 and SETBP1 replicated in African Americans. In summary, rs3922844 was associated with QRS duration and SCN5A expression, two novel QRS loci were identified using transethnic meta-analysis, and a significant proportion of QRS-SNP associations discovered in populations of European descent were transferable to African Americans when adequate power was achieved.

Decline in muscle strength with aging is an important predictor of health trajectory in the elderly. Several factors, including genetics, are proposed contributors to variability in muscle strength. To identify genetic contributors to muscle strength, a meta-analysis of genomewide association studies of handgrip was conducted. Grip strength was measured using a handheld dynamometer in 27 581 individuals of European descent over 65 years of age from 14 cohort studies. Genomewide association analysis was conducted on ~2.7 million imputed and genotyped variants (SNPs). Replication of the most significant findings was conducted using data from 6393 individuals from three cohorts. GWAS of lower body strength was also characterized in a subset of cohorts. Two genomewide significant (P-value< 5 × 10(-8) ) and 39 suggestive (P-value< 5 × 10(-5) ) associations were observed from meta-analysis of the discovery cohorts. After meta-analysis with replication cohorts, genomewide significant association was observed for rs752045 on chromosome 8 (β = 0.47, SE = 0.08, P-value = 5.20 × 10(-10) ). This SNP is mapped to an intergenic region and is located within an accessible chromatin region (DNase hypersensitivity site) in skeletal muscle myotubes differentiated from the human skeletal muscle myoblasts cell line. This locus alters a binding motif of the CCAAT/enhancer-binding protein-β (CEBPB) that is implicated in muscle repair mechanisms. GWAS of lower body strength did not yield significant results. A common genetic variant in a chromosomal region that regulates myotube differentiation and muscle repair may contribute to variability in grip strength in the elderly. Further studies are needed to uncover the mechanisms that link this genetic variant with muscle strength.

RationaleObstructive sleep apnea is a common disorder associated with increased risk for cardiovascular disease, diabetes, and premature mortality. Although there is strong clinical and epidemiologic evidence supporting the importance of genetic factors in influencing obstructive sleep apnea, its genetic basis is still largely unknown. Prior genetic studies focused on traits defined using the apnea-hypopnea index, which contains limited information on potentially important genetically determined physiologic factors, such as propensity for hypoxemia and respiratory arousability.ObjectivesTo define novel obstructive sleep apnea genetic risk loci for obstructive sleep apnea, we conducted genome-wide association studies of quantitative traits in Hispanic/Latino Americans from three cohorts.MethodsGenome-wide data from as many as 12,558 participants in the Hispanic Community Health Study/Study of Latinos, Multi-Ethnic Study of Atherosclerosis, and Starr County Health Studies population-based cohorts were metaanalyzed for association with the apnea-hypopnea index, average oxygen saturation during sleep, and average respiratory event duration.Measurements and main resultsTwo novel loci were identified at genome-level significance (rs11691765, GPR83, P = 1.90 × 10-8 for the apnea-hypopnea index, and rs35424364; C6ORF183/CCDC162P, P = 4.88 × 10-8 for respiratory event duration) and seven additional loci were identified with suggestive significance (P

Knowledge of the genetic basis of the type 2 diabetes (T2D)-related quantitative traits fasting glucose (FG) and insulin (FI) in African ancestry (AA) individuals has been limited. In non-diabetic subjects of AA (n = 20,209) and European ancestry (EA; n = 57,292), we performed trans-ethnic (AA+EA) fine-mapping of 54 established EA FG or FI loci with detailed functional annotation, assessed their relevance in AA individuals, and sought previously undescribed loci through trans-ethnic (AA+EA) meta-analysis. We narrowed credible sets of variants driving association signals for 22/54 EA-associated loci; 18/22 credible sets overlapped with active islet-specific enhancers or transcription factor (TF) binding sites, and 21/22 contained at least one TF motif. Of the 54 EA-associated loci, 23 were shared between EA and AA. Replication with an additional 10,096 AA individuals identified two previously undescribed FI loci, chrX FAM133A (rs213676) and chr5 PELO (rs6450057). Trans-ethnic analyses with regulatory annotation illuminate the genetic architecture of glycemic traits and suggest gene regulation as a target to advance precision medicine for T2D. Our approach to utilize state-of-the-art functional annotation and implement trans-ethnic association analysis for discovery and fine-mapping offers a framework for further follow-up and characterization of GWAS signals of complex trait loci.

Single nucleotide polymorphisms (SNPs) in and near ABCA7, BIN1, CASS4, CD2AP, CD33, CELF1, CLU, complement receptor 1 (CR1), EPHA1, EXOC3L2, FERMT2, HLA cluster (DRB5-DQA), INPP5D, MEF2C, MS4A cluster (MS4A3-MS4A6E), NME8, PICALM, PTK2B, SLC24A4, SORL1, and ZCWPW1 have been associated with Alzheimer's disease (AD) in large meta-analyses. We aimed to determine whether established AD-associated genes are associated with longitudinal cognitive decline by examining aggregate variation across these gene regions. In two single-sex cohorts of older, community-dwelling adults, we examined the association between SNPs in previously implicated gene regions and cognitive decline (age-adjusted person-specific cognitive slopes) using a Sequence Kernel Association Test (SKAT). In regions which showed aggregate significance, we examined the univariate association between individual SNPs in the region and cognitive decline. Only two of the original AD-associated SNPs were significantly associated with cognitive decline in our cohorts. We identified significant aggregate-level associations between cognitive decline and the gene regions BIN1, CD33, CELF1, CR1, HLA cluster, and MEF2C in the all-female cohort and significant associations with ABCA7, HLA cluster, MS4A6E, PICALM, PTK2B, SLC24A4, and SORL1 in the all-male cohort. We also identified a block of eight correlated SNPs in CD33 and several blocks of correlated SNPs in CELF1 that were significantly associated with cognitive decline in univariate analysis in the all-female cohort.

Sleep duration is implicated in the etiologies of chronic diseases and premature mortality. However, the genetic basis for sleep duration is poorly defined. We sought to identify novel genetic components influencing sleep duration in a multi-ethnic sample. Meta-analyses were conducted of genetic associations with self-reported, habitual sleep duration from seven Candidate Gene Association Resource (CARe) cohorts of over 25 000 individuals of African, Asian, European and Hispanic American ancestry. All individuals were genotyped for ∼50 000 SNPs from 2000 candidate heart, lung, blood and sleep genes. African-Americans had additional genome-wide genotypes. Four cohorts provided replication. A SNP (rs17601612) in the dopamine D2 receptor gene (DRD2) was significantly associated with sleep duration (P = 9.8 × 10(-7)). Conditional analysis identified a second DRD2 signal with opposite effects on sleep duration. In exploratory analysis, suggestive association was observed for rs17601612 with polysomnographically determined sleep latency (P = 0.002). The lead DRD2 signal was recently identified in a schizophrenia GWAS, and a genetic risk score of 11 additional schizophrenia GWAS loci genotyped on the IBC array was also associated with longer sleep duration (P = 0.03). These findings support a role for DRD2 in influencing sleep duration. Our work motivates future pharmocogenetics research on alerting agents such as caffeine and modafinil that interact with the dopaminergic pathway and further investigation of genetic overlap between sleep and neuro-psychiatric traits.

The extent to which low-frequency (minor allele frequency (MAF) between 1-5%) and rare (MAF ≤ 1%) variants contribute to complex traits and disease in the general population is mainly unknown. Bone mineral density (BMD) is highly heritable, a major predictor of osteoporotic fractures, and has been previously associated with common genetic variants, as well as rare, population-specific, coding variants. Here we identify novel non-coding genetic variants with large effects on BMD (ntotal = 53,236) and fracture (ntotal = 508,253) in individuals of European ancestry from the general population. Associations for BMD were derived from whole-genome sequencing (n = 2,882 from UK10K (ref. 10); a population-based genome sequencing consortium), whole-exome sequencing (n = 3,549), deep imputation of genotyped samples using a combined UK10K/1000 Genomes reference panel (n = 26,534), and de novo replication genotyping (n = 20,271). We identified a low-frequency non-coding variant near a novel locus, EN1, with an effect size fourfold larger than the mean of previously reported common variants for lumbar spine BMD (rs11692564(T), MAF = 1.6%, replication effect size = +0.20 s.d., Pmeta = 2 × 10(-14)), which was also associated with a decreased risk of fracture (odds ratio = 0.85; P = 2 × 10(-11); ncases = 98,742 and ncontrols = 409,511). Using an En1(cre/flox) mouse model, we observed that conditional loss of En1 results in low bone mass, probably as a consequence of high bone turnover. We also identified a novel low-frequency non-coding variant with large effects on BMD near WNT16 (rs148771817(T), MAF = 1.2%, replication effect size = +0.41 s.d., Pmeta = 1 × 10(-11)). In general, there was an excess of association signals arising from deleterious coding and conserved non-coding variants. These findings provide evidence that low-frequency non-coding variants have large effects on BMD and fracture, thereby providing rationale for whole-genome sequencing and improved imputation reference panels to study the genetic architecture of complex traits and disease in the general population.

BackgroundDepressive symptoms are common in older adults and associated with poor outcomes. Although circadian genes have been implicated in depression, the relationship between circadian genes and depressive symptoms in older adults is unclear.MethodsA cross-sectional genetic association study of 529 single nucleotide polymorphisms (SNPs) representing 30 candidate circadian genes was performed in two population-based cohorts: the Osteoporotic Fractures in Men Study (MrOS; N=270, age: 76.58±5.61 years) and the Study of Osteoporotic Fractures (SOF) in women (N=1740, 84.05±3.53 years) and a meta-analysis was performed. Depressive symptoms were assessed with the Geriatric Depression Scale categorizing participants as having none-few symptoms (0-2), some depressive symptoms (>2 to

ObjectivesTo identify genetic associations with hip osteoarthritis (HOA), we performed a meta-analysis of genome-wide association studies (GWAS) of HOA.MethodsThe GWAS meta-analysis included approximately 2.5 million imputed HapMap single nucleotide polymorphisms (SNPs). HOA cases and controls defined radiographically and by total hip replacement were selected from the Osteoporotic Fractures in Men (MrOS) Study and the Study of Osteoporotic Fractures (SOF) (654 cases and 4697 controls, combined). Replication of genome-wide significant SNP associations (p ≤5×10(-8)) was examined in five studies (3243 cases and 6891 controls, combined). Functional studies were performed using in vitro models of chondrogenesis and osteogenesis.ResultsThe A allele of rs788748, located 65 kb upstream of the IGFBP3 gene, was associated with lower HOA odds at the genome-wide significance level in the discovery stage (OR 0.71, p=2×10(-8)). The association replicated in five studies (OR 0.92, p=0.020), but the joint analysis of discovery and replication results was not genome-wide significant (p=1×10(-6)). In separate study populations, the rs788748 A allele was also associated with lower circulating IGFBP3 protein levels (p=4×10(-13)), suggesting that this SNP or a variant in linkage disequilibrium could be an IGFBP3 regulatory variant. Results from functional studies were consistent with association results. Chondrocyte hypertrophy, a deleterious event in OA pathogenesis, was largely prevented upon IGFBP3 knockdown in chondrocytes. Furthermore, IGFBP3 overexpression induced cartilage catabolism and osteogenic differentiation.ConclusionsResults from GWAS and functional studies provided suggestive links between IGFBP3 and HOA.