Genome-wide characterization of circulating metabolic biomarkers.

Genome-wide association analyses using high-throughput metabolomics platforms have led to novel insights into the biology of human metabolism1-7. This detailed knowledge of the genetic determinants of systemic metabolism has been pivotal for uncovering how genetic pathways influence biological mechanisms and complex diseases8-11. Here we present a genome-wide association study for 233 circulating metabolic traits quantified by nuclear magnetic resonance spectroscopy in up to 136,016 participants from 33 cohorts. We identify more than 400 independent loci and assign probable causal genes at two-thirds of these using manual curation of plausible biological candidates. We highlight the importance of sample and participant characteristics that can have significant effects on genetic associations. We use detailed metabolic profiling of lipoprotein- and lipid-associated variants to better characterize how known lipid loci and novel loci affect lipoprotein metabolism at a granular level. We demonstrate the translational utility of comprehensively phenotyped molecular data, characterizing the metabolic associations of intrahepatic cholestasis of pregnancy. Finally, we observe substantial genetic pleiotropy for multiple metabolic pathways and illustrate the importance of careful instrument selection in Mendelian randomization analysis, revealing a putative causal relationship between acetone and hypertension. Our publicly available results provide a foundational resource for the community to examine the role of metabolism across diverse diseases.

Genome-wide association meta-analysis identifies 48 risk variants and highlights the role of the stria vascularis in hearing loss.

Hearing loss is one of the top contributors to years lived with disability and is a risk factor for dementia. Molecular evidence on the cellular origins of hearing loss in humans is growing. Here, we performed a genome-wide association meta-analysis of clinically diagnosed and self-reported hearing impairment on 723,266 individuals and identified 48 significant loci, 10 of which are novel. A large proportion of associations comprised missense variants, half of which lie within known familial hearing loss loci. We used single-cell RNA-sequencing data from mouse cochlea and brain and mapped common-variant genomic results to spindle, root, and basal cells from the stria vascularis, a structure in the cochlea necessary for normal hearing. Our findings indicate the importance of the stria vascularis in the mechanism of hearing impairment, providing future paths for developing targets for therapeutic intervention in hearing loss.

Associations of carotid intima media thickness with gene expression in whole blood and genetically predicted gene expression across 48 tissues.

Carotid intima media thickness (cIMT) is a biomarker of subclinical atherosclerosis and a predictor of future cardiovascular events. Identifying associations between gene expression levels and cIMT may provide insight to atherosclerosis etiology. Here, we use two approaches to identify associations between mRNA levels and cIMT: differential gene expression analysis in whole blood and S-PrediXcan. We used microarrays to measure genome-wide whole blood mRNA levels of 5647 European individuals from four studies. We examined the association of mRNA levels with cIMT adjusted for various potential confounders. Significant associations were tested for replication in three studies totaling 3943 participants. Next, we applied S-PrediXcan to summary statistics from a cIMT genome-wide association study (GWAS) of 71 128 individuals to estimate the association between genetically determined mRNA levels and cIMT and replicated these analyses using S-PrediXcan on an independent GWAS on cIMT that included 22 179 individuals from the UK Biobank. mRNA levels of TNFAIP3, CEBPD and METRNL were inversely associated with cIMT, but these associations were not significant in the replication analysis. S-PrediXcan identified associations between cIMT and genetically determined mRNA levels for 36 genes, of which six were significant in the replication analysis, including TLN2, which had not been previously reported for cIMT. There was weak correlation between our results using differential gene expression analysis and S-PrediXcan. Differential expression analysis and S-PrediXcan represent complementary approaches for the discovery of associations between phenotypes and gene expression. Using these approaches, we prioritize TNFAIP3, CEBPD, METRNL and TLN2 as new candidate genes whose differential expression might modulate cIMT.

Plasma proteomic signature of fatty liver disease: The Rotterdam Study.

BACKGROUND AND AIMS: Fatty liver disease (FLD) is caused by excess fat in the liver, and its global prevalence exceeds 33%. The role of protein expression on the pathogenesis of FLD and accompanied fibrosis and its potential as a disease biomarker is currently not clear. Hence, we aimed to identify plasma proteomics associated with FLD and fibrosis using population-based data. APPROACH AND RESULTS: Blood samples were collected from 2578 participants from the population-based Rotterdam Study cohort. The proximity extension assay reliably measured plasma levels of 171 cardiometabolic and inflammatory-related proteins (Olink Proteomics). FLD was assessed by ultrasound, and fibrosis by transient elastography. Logistic regression models quantified the association of plasma proteomics with FLD and fibrosis. In addition, we aimed to validate our results in liver organoids. The cross-sectional analysis identified 27 proteins significantly associated with FLD surpassing the Bonferroni-corrected p <2.92×10 -4 . The strongest association was observed for FGF-21 (ß=0.45, p =1.07×10 -18 ) and carboxylesterase 1 (CES1) protein (ß=0.66, p =4.91×10 -40 ). Importantly, 15 of the 27 proteins significantly associated with FLD were also associated with liver fibrosis. Finally, consistent with plasma proteomic profiling, we found the expression levels of IL-18 receptor 1 (IL-18R1) and CES1 to be upregulated in an FLD model of 3-dimensional culture human liver organoids. CONCLUSIONS: Among the general population, several inflammatory and cardiometabolic plasma proteins were associated with FLD and fibrosis. Particularly, plasma levels of FGF-21, IL-18R1, and CES1 were largely dependent on the presence of FLD and fibrosis and may therefore be important in their pathogenesis.

The Rotterdam Study. Design update and major findings between 2020 and 2024.

The Rotterdam Study is a population-based cohort study, started in 1990 in the district of Ommoord in the city of Rotterdam, the Netherlands, with the aim to describe the prevalence and incidence, unravel the etiology, and identify targets for prediction, prevention or intervention of multifactorial diseases in mid-life and elderly. The study currently includes 17,931 participants (overall response rate 65%), aged 40 years and over, who are examined in-person every 3 to 5 years in a dedicated research facility, and who are followed-up continuously through automated linkage with health care providers, both regionally and nationally. Research within the Rotterdam Study is carried out along two axes. First, research lines are oriented around diseases and clinical conditions, which are reflective of medical specializations. Second, cross-cutting research lines transverse these clinical demarcations allowing for inter- and multidisciplinary research. These research lines generally reflect subdomains within epidemiology. This paper describes recent methodological updates and main findings from each of these research lines. Also, future perspective for coming years highlighted.

A functional genomics approach reveals suggestive quantitative trait loci associated with combined TLR4 and BCP crystal-induced inflammation and osteoarthritis.

OBJECTIVE: Basic calcium phosphate (BCP) crystals can activate the NLRP3 inflammasome and are potentially involved in the pathogenesis of osteoarthritis (OA). In order to elucidate relevant inflammatory mechanisms in OA, we used a functional genomics approach to assess genetic variation influencing BCP crystal-induced cytokine production. METHOD: Peripheral blood mononuclear cells (PBMCs) were isolated from healthy volunteers who were previously genotyped and stimulated with BCP crystals and/or lipopolysaccharide (LPS) after which cytokines release was assessed. Cytokine quantitative trait locus (cQTL) mapping was performed. For in vitro validation of the cQTL located in anoctamin 3 (ANO3), PBMCs were incubated with Tamoxifen and Benzbromarone prior to stimulation. Additionally, we performed co-localisation analysis of our top cQTLs with the most recent OA meta-analysis of genome-wide association studies (GWAS). RESULTS: We observed that BCP crystals and LPS synergistically induce IL-1ß in human PBMCs. cQTL analysis revealed several suggestive loci influencing cytokine release upon stimulation, among which are quantitative trait locus annotated to ANO3 and GLIS3. As functional validation, anoctamin inhibitors reduced IL-1ß release in PBMCs after stimulation. Co-localisation analysis showed that the GLIS3 locus was shared between LPS/BCP crystal-induced IL-1ß and genetic association with Knee OA. CONCLUSIONS: We identified and functionally validated a new locus, ANO3, associated with LPS/BCP crystal-induced inflammation in PBMCs. Moreover, the cQTL in the GLIS3 locus co-localises with the previously found locus associated with Knee OA, suggesting that this Knee OA locus might be explained through an inflammatory mechanism. These results form a basis for further exploration of inflammatory mechanisms in OA.

A systematic review of the sex differences in risk factors for knee osteoarthritis.

OBJECTIVES: Previous systematic reviews focused on the evidence of common risk factors for knee OA (KOA); however, the effect and strength of association between risk factors and KOA might be different between the two sexes. The aim of the present systematic review was to determine the current evidence on sex differences in the association between risk factors and KOA and their prevalence. METHODS: We searched the following electronic bibliographic databases: MEDLINE (PubMed), EMBASE and Web of Science. A methodological quality assessment was conducted independently by two researchers according to an adapted version of the standardized set of criteria known as the Newcastle-Ottawa Quality Assessment Scale (NOS). The NOS, a star system, was converted to three categories of quality. RESULTS: In total, 27 studies reported sex-specific risk estimates on several risk factors for KOA. Out of the 22 longitudinal cohort studies (except one nested case-control), 12 were of good quality and 10 were of fair quality. The five cross-sectional studies consisted of one of good, three of fair and one of poor quality. There was an indication of sex differences in risk factors leading to higher risk of KOA: high BMI, alcohol consumption, atherosclerosis and high vitamin E levels in women, and high physical activity, soft drink consumption and abdominal obesity in men. Knee injury, high blood pressure and low step rate seem to affect both women and men. CONCLUSION: More good quality studies are needed to assess sex differences in risk factors for KOA, especially for symptomatic/clinical OA.

Plasma proteomics identifies CRTAC1 as a biomarker for osteoarthritis severity and progression.

OBJECTIVES: The aim of this study was to identify biomarkers for radiographic OA severity and progression acting within the inflammation and metabolic pathways. METHODS: For 3517 Rotterdam Study participants, 184 plasma protein levels were measured using Olink inflammation and cardiometabolic panels. We studied associations with severity and progression of knee, hip and hand OA and a composite overall OA burden score by multivariable regression models, adjusting for age, sex, cell counts and BMI. RESULTS: We found 18 significantly associated proteins for overall OA burden, of which 5 stayed significant after multiple testing correction: circulating cartilage acidic protein 1 (CRTAC1), cartilage oligomeric matrix protein (COMP), thrombospondin 4, IL-18 receptor 1 (IL-18R1) and TNF ligand superfamily member 14. These proteins were also associated with progression of knee OA, with the exception of IL-18R1. The strongest association was found for the level of CRTAC1, with 1 s.d. increase in protein level resulting in an increase of 0.09 (95% CI 0.06, 0.12) in the overall OA Kellgren-Lawrence sum score (P = 2.9 × 10-8) in the model adjusted for age, sex, BMI and cell counts. This association was also present with the severity of OA in all three joints and progression of knee OA and was independent of BMI. We observed a stronger association for CRTAC1 with OA than for the well-known OA biomarker COMP. CONCLUSION: We identified several compelling biomarkers reflecting the overall OA burden and the increased risk for OA progression. CRTAC1 was the most compelling and robust biomarker for OA severity and progression. Such a biomarker may be used for disease monitoring.

Younger facial looks are associate with a lower likelihood of several age-related morbidities in the middle-aged to elderly.

BACKGROUND: Looking older for one's chronological age is associated with a higher mortality rate. Yet it remains unclear how perceived facial age relates to morbidity and the degree to which facial ageing reflects systemic ageing of the human body. OBJECTIVES: To investigate the association between ΔPA and age-related morbidities of different organ systems, where ΔPA represents the difference between perceived age (PA) and chronological age. METHODS: We performed a cross-sectional analysis on data from the Rotterdam Study, a population-based cohort study in the Netherlands. High-resolution facial photographs of 2679 men and women aged 51.5-87.8 years of European descent were used to assess PA. PA was estimated and scored in 5-year categories using these photographs by a panel of men and women who were blinded for chronological age and medical history. A linear mixed model was used to generate the mean PAs. The difference between the mean PA and chronological age was calculated (ΔPA), where a higher (positive) ΔPA means that the person looks younger for their age and a lower (negative) ΔPA that the person looks older. ΔPA was tested as a continuous variable for association with ageing-related morbidities including cardiovascular, pulmonary, ophthalmological, neurocognitive, renal, skeletal and auditory morbidities in separate regression analyses, adjusted for age and sex (model 1) and additionally for body mass index, smoking and sun exposure (model 2). RESULTS: We observed 5-year higher ΔPA (i.e. looking younger by 5 years for one's age) to be associated with less osteoporosis [odds ratio (OR) 0.76, 95% confidence interval (CI) 0.62-0.93], less chronic obstructive pulmonary disease (OR 0.85, 95% CI 0.77-0.95), less age-related hearing loss (model 2; B = -0.76, 95% CI -1.35 to -0.17) and fewer cataracts (OR 0.84, 95% CI 0.73-0.97), but with better global cognitive functioning (g-factor; model 2; B = 0.07, 95% CI 0.04-0.10). CONCLUSIONS: PA is associated with multiple morbidities and better cognitive function, suggesting that systemic ageing and cognitive ageing are, to an extent, externally visible in the human face.

Towards sex-specific osteoarthritis risk models: evaluation of risk factors for knee osteoarthritis in males and females.

OBJECTIVES: The aim of this study was to identify sex-specific prevalence and strength of risk factors for the incidence of radiographic knee OA (incRKOA). METHODS: Our study population consisted of 10 958 Rotterdam Study participants free of knee OA in one or both knees at baseline. One thousand and sixty-four participants developed RKOA after a median follow-up time of 9.6 years. We estimated the association between each available risk factor and incRKOA using sex stratified multivariate regression models with generalized estimating equations. Subsequently, we statistically tested sex differences between risk estimates and calculated the population attributable fractions (PAFs) for modifiable risk factors. RESULTS: The prevalence of the investigated risk factors was, in general, higher in women compared with men, except that alcohol intake and smoking were higher in men and high BMI showed equal prevalence. We found significantly different risk estimates between men and women: high level of physical activity [relative risk (RR) 1.76 (95% CI: 1.29-2.40)] or a Kellgren and Lawrence score 1 at baseline [RR 5.48 (95% CI: 4.51-6.65)] was higher in men. Among borderline significantly different risk estimates was BMI ≥27, associated with higher risk for incRKOA in women [RR 2.00 (95% CI: 1.74-2.31)]. The PAF for higher BMI was 25.6% in women and 19.3% in men. CONCLUSION: We found sex-specific differences in both presence and relative risk of several risk factors for incRKOA. Especially BMI, a modifiable risk factor, impacts women more strongly than men. These risk factors can be used in the development of personalized prevention strategies and in building sex-specific prediction tools to identify high risk profile patients.

Shorter periconception maternal telomere length and the risk of congenital cardiac outflow defects in the offspring.

BACKGROUND: Congenital cardiac outflow defects (COD) are the largest group of congenital heart defects, with ventricular septal defect (VSD) as the most prevalent phenotype. Increased maternal age, excessive oxidative stress and inflammation are involved in the pathophysiology of COD and enhance telomere length (TL) shortening. We investigated the association between periconception maternal TL and the risk of having a child with COD. METHODS: From a multicentre case-control trial, 306 case mothers of a child with COD and 424 control mothers of a child without a congenital malformation were selected. Relative TL was measured by qPCR. Multivariable logistic regression was used to compute crude and adjusted odds ratios, per standard deviation decrease, between maternal T/S ratio and COD and VSD risk. Adjustments were made for maternal age. Additional adjustments were made in a second model. RESULTS: Shorter maternal relative TL was significantly associated with an OR of 1.29 (95% CI 1.04-1.61), p = .02, for the risk of VSD in offspring, which remained significant after an adjustment for maternal age (adjOR 1.25(95% CI 1.01-1.55), p = .04). No association between maternal TL and the risk of overall COD in offspring was observed. CONCLUSION: Shorter maternal relative TL is associated with an approximately 1.3-OR for the risk, per SD in relative TL shortening, of VSD in the offspring. These findings need further confirmation in other studies on the predictive value of maternal TL.

The interrelatedness of chronic cough and chronic pain.

Since chronic cough has common neurobiological mechanisms and pathophysiology with chronic pain, both clinical disorders might be interrelated. Hence, we examined the association between chronic cough and chronic pain in adult subjects in the Rotterdam Study, a large prospective population-based cohort study.Using a standardised questionnaire, chronic pain was defined as pain lasting up to 6 months and grouped into a frequency of weekly/monthly or daily pain. Chronic cough was described as daily coughing for at least 3 months duration. The longitudinal and cross-sectional associations were investigated bi-directionally.Of 7141 subjects in the study, 54% (n=3888) reported chronic pain at baseline. The co-prevalence of daily chronic pain and chronic cough was 4.4%. Chronic cough was more prevalent in subjects with daily and weekly/monthly chronic pain compared with those without chronic pain (13.8% and 10.3% versus 8.2%; p<0.001). After adjustment for potential confounders, prevalent chronic pain was significantly associated with incident chronic cough (OR 1.47, 95% CI 1.08-1.99). The association remained significant in subjects with daily chronic pain (OR 1.49, 95% CI 1.06-2.11) with a similar effect estimate, albeit non-significant in those with weekly/monthly chronic pain (OR 1.43, 95% CI 0.98-2.10). After adjustment for covariables, subjects with chronic cough had a significant risk of developing chronic pain (OR 1.63, 95% CI 1.02-2.62) compared with those without chronic cough.Chronic cough and chronic pain confer risk on each other among adult subjects, indicating that both conditions might share common risk factors and/or pathophysiologic mechanisms.

Vitamin K antagonist anticoagulant usage is associated with increased incidence and progression of osteoarthritis.

OBJECTIVES: Vitamin K is hypothesised to play a role in osteoarthritis (OA) pathogenesis through effects on vitamin K-dependent bone and cartilage proteins, and therefore may represent a modifiable risk factor. A genetic variant in a vitamin K-dependent protein that is an essential inhibitor for cartilage calcification, matrix Gla protein (MGP), was associated with an increased risk for OA. Vitamin K antagonist anticoagulants (VKAs), such as warfarin and acenocoumarol, act as anticoagulants through inhibition of vitamin K-dependent blood coagulation proteins. VKAs likely also affect the functioning of other vitamin K-dependent proteins such as MGP. METHODS: We investigated the effect of acenocoumarol usage on progression and incidence of radiographic OA in 3494 participants of the Rotterdam Study cohort. We also examined the effect of MGP and VKORC1 single nucleotide variants on this association. RESULTS: Acenocoumarol usage was associated with an increased risk of OA incidence and progression (OR=2.50, 95% CI=1.94-3.20), both for knee (OR=2.34, 95% CI=1.67-3.22) and hip OA (OR=2.74, 95% CI=1.82-4.11). Among acenocoumarol users, carriers of the high VKORC1(BB) expression haplotype together with the MGP OA risk allele (rs1800801-T) had an increased risk of OA incidence and progression (OR=4.18, 95% CI=2.69-6.50), while this relationship was not present in non-users of that group (OR=1.01, 95% CI=0.78-1.33). CONCLUSIONS: These findings support the importance of vitamin K and vitamin K-dependent proteins, as MGP, in the pathogenesis of OA. Additionally, these results may have direct implications for the clinical prevention of OA, supporting the consideration of direct oral anticoagulants in favour of VKAs.

Genome-wide association study identifies RNF123 locus as associated with chronic widespread musculoskeletal pain.

BACKGROUND AND OBJECTIVES: Chronic widespread musculoskeletal pain (CWP) is a symptom of fibromyalgia and a complex trait with poorly understood pathogenesis. CWP is heritable (48%-54%), but its genetic architecture is unknown and candidate gene studies have produced inconsistent results. We conducted a genome-wide association study to get insight into the genetic background of CWP. METHODS: Northern Europeans from UK Biobank comprising 6914 cases reporting pain all over the body lasting >3 months and 242 929 controls were studied. Replication of three independent genome-wide significant single nucleotide polymorphisms was attempted in six independent European cohorts (n=43 080; cases=14 177). Genetic correlations with risk factors, tissue specificity and colocalisation were examined. RESULTS: Three genome-wide significant loci were identified (rs1491985, rs10490825, rs165599) residing within the genes Ring Finger Protein 123 (RNF123), ATPase secretory pathway Ca2+transporting 1 (ATP2C1) and catechol-O-methyltransferase (COMT). The RNF123 locus was replicated (meta-analysis p=0.0002), the ATP2C1 locus showed suggestive association (p=0.0227) and the COMT locus was not replicated. Partial genetic correlation between CWP and depressive symptoms, body mass index, age of first birth and years of schooling were identified. Tissue specificity and colocalisation analysis highlight the relevance of skeletal muscle in CWP. CONCLUSIONS: We report a novel association of RNF123 locus and a suggestive association of ATP2C1 locus with CWP. Both loci are consistent with a role of calcium regulation in CWP. The association with COMT, one of the most studied genes in chronic pain field, was not confirmed in the replication analysis.

Genome-wide association of phenotypes based on clustering patterns of hand osteoarthritis identify WNT9A as novel osteoarthritis gene.

BACKGROUND: Despite recent advances in the understanding of the genetic architecture of osteoarthritis (OA), only two genetic loci have been identified for OA of the hand, in part explained by the complexity of the different hand joints and heterogeneity of OA pathology. METHODS: We used data from the Rotterdam Study (RSI, RSII and RSIII) to create three hand OA phenotypes based on clustering patterns of radiographic OA severity to increase power in our modest discovery genome-wide association studies in the RS (n=8700), and sought replication in an independent cohort, the Framingham Heart Study (n=1203). We used multiple approaches that leverage different levels of information and functional data to further investigate the underlying biological mechanisms and candidate genes for replicated loci. We also attempted to replicate known OA loci at other joint sites, including the hips and knees. RESULTS: We found two novel genome-wide significant loci for OA in the thumb joints. We identified WNT9A as a possible novel causal gene involved in OA pathogenesis. Furthermore, several previously identified genetic loci for OA seem to confer risk for OA across multiple joints: TGFa, RUNX2, COL27A1, ASTN2, IL11 and GDF5 loci. CONCLUSIONS: We identified a robust novel genetic locus for hand OA on chromosome 1, of which WNT9A is the most likely causal gene. In addition, multiple genetic loci were identified to be associated with OA across multiple joints. Our study confirms the potential for novel insight into the genetic architecture of OA by using biologically meaningful stratified phenotypes.

Mendelian randomization study on vitamin D levels and osteoarthritis risk: a concise report.

OBJECTIVE: The role of vitamin D in OA is unclear and previous epidemiological studies have provided inconsistent results. We conducted a two-sample Mendelian randomization (MR) study to investigate the causal relationship between genetically determined serum vitamin D levels and hip/knee OA. METHODS: Six single-nucleotide polymorphisms (SNPs) associated with vitamin D levels in the Study of Underlying Genetic Determinants of Vitamin D and Highly Related Traits Consortium were selected as instrumental variables. Summary statistics of the SNPs effects on OA were derived from the Iceland and UK Biobank, comprising 23 877 knee OA cases, 17 151 hip OA cases and >562 000 controls. The control samples match the OA cases in age, sex and county of origin. RESULTS: The MR analyses showed no causal association between genetically determined vitamin D levels and knee OA [odds ratio (OR) 1.03 (95% CI 0.84, 1.26)] or hip OA [OR 1.06 (95% CI 0.83, 1.35)]. CONCLUSION: Genetic variations associated with low vitamin D serum levels are not associated with increased risk of hip or knee OA in community-dwelling older adults, suggesting that vitamin D levels are not causally linked to OA. It is therefore unlikely that vitamin D supplementation protects against hip or knee OA.

Somatic TARDBP variants as a cause of semantic dementia.

The aetiology of late-onset neurodegenerative diseases is largely unknown. Here we investigated whether de novo somatic variants for semantic dementia can be detected, thereby arguing for a more general role of somatic variants in neurodegenerative disease. Semantic dementia is characterized by a non-familial occurrence, early onset (<65 years), focal temporal atrophy and TDP-43 pathology. To test whether somatic variants in neural progenitor cells during brain development might lead to semantic dementia, we compared deep exome sequencing data of DNA derived from brain and blood of 16 semantic dementia cases. Somatic variants observed in brain tissue and absent in blood were validated using amplicon sequencing and digital PCR. We identified two variants in exon one of the TARDBP gene (L41F and R42H) at low level (1-3%) in cortical regions and in dentate gyrus in two semantic dementia brains, respectively. The pathogenicity of both variants is supported by demonstrating impaired splicing regulation of TDP-43 and by altered subcellular localization of the mutant TDP-43 protein. These findings indicate that somatic variants may cause semantic dementia as a non-hereditary neurodegenerative disease, which might be exemplary for other late-onset neurodegenerative disorders.

Genome-wide meta-analysis of 158,000 individuals of European ancestry identifies three loci associated with chronic back pain.

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<5×10(-8). Suggestive (p<5×10(-7)) and genome-wide significant (p<5×10(-8)) variants were carried forward for replication or further investigation in the remaining UK Biobank participants not included in the discovery sample. The discovery sample comprised 158,025 individuals, including 29,531 CBP cases. A genome-wide significant association was found for the intronic variant rs12310519 in SOX5 (OR 1.08, p = 7.2×10(-10)). This was subsequently replicated in 283,752 UK Biobank participants not included in the discovery sample, including 50,915 cases (OR 1.06, p = 5.3×10(-11)), and exceeded genome-wide significance in joint meta-analysis (OR 1.07, p = 4.5×10(-19)). We found suggestive associations at three other loci in the discovery sample, two of which exceeded genome-wide significance in joint meta-analysis: an intergenic variant, rs7833174, located between CCDC26 and GSDMC (OR 1.05, p = 4.4×10(-13)), and an intronic variant, rs4384683, in DCC (OR 0.97, p = 2.4×10(-10)). In this first reported meta-analysis of GWAS for CBP, we identified and replicated a genetic locus associated with CBP (SOX5). We also identified 2 other loci that reached genome-wide significance in a 2-stage joint meta-analysis (CCDC26/GSDMC and DCC).

COPD GWAS variant at 19q13.2 in relation with DNA methylation and gene expression.

Chronic obstructive pulmonary disease (COPD) is among the major health burdens in adults. While cigarette smoking is the leading risk factor, a growing number of genetic variations have been discovered to influence disease susceptibility. Epigenetic modifications may mediate the response of the genome to smoking and regulate gene expression. Chromosome 19q13.2 region is associated with both smoking and COPD, yet its functional role is unclear. Our study aimed to determine whether rs7937 (RAB4B, EGLN2), a top genetic variant in 19q13.2 region identified in genome-wide association studies of COPD, is associated with differential DNA methylation in blood (N = 1490) and gene expression in blood (N = 721) and lungs (N = 1087). We combined genetic and epigenetic data from the Rotterdam Study (RS) to perform the epigenome-wide association analysis of rs7937. Further, we used genetic and transcriptomic data from blood (RS) and from lung tissue (Lung expression quantitative trait loci mapping study), to perform the transcriptome-wide association study of rs7937. Rs7937 was significantly (FDR < 0.05) and consistently associated with differential DNA methylation in blood at 4 CpG sites in cis, independent of smoking. One methylation site (cg11298343-EGLN2) was also associated with COPD (P = 0.001). Additionally, rs7937 was associated with gene expression levels in blood in cis (EGLN2), 42% mediated through cg11298343, and in lung tissue, in cis and trans (NUMBL, EGLN2, DNMT3A, LOC101929709 and PAK2). Our results suggest that changes of DNA methylation and gene expression may be intermediate steps between genetic variants and COPD, but further causal studies in lung tissue should confirm this hypothesis.

Genome-wide association study identifies nine novel loci for 2D:4D finger ratio, a putative retrospective biomarker of testosterone exposure in utero.

The ratio of the length of the index finger to that of the ring finger (2D:4D) is sexually dimorphic and is commonly used as a non-invasive biomarker of prenatal androgen exposure. Most association studies of 2D:4D ratio with a diverse range of sex-specific traits have typically involved small sample sizes and have been difficult to replicate, raising questions around the utility and precise meaning of the measure. In the largest genome-wide association meta-analysis of 2D:4D ratio to date (N = 15 661, with replication N = 75 821), we identified 11 loci (9 novel) explaining 3.8% of the variance in mean 2D:4D ratio. We also found weak evidence for association (ß = 0.06; P = 0.02) between 2D:4D ratio and sensitivity to testosterone [length of the CAG microsatellite repeat in the androgen receptor (AR) gene] in females only. Furthermore, genetic variants associated with (adult) testosterone levels and/or sex hormone-binding globulin were not associated with 2D:4D ratio in our sample. Although we were unable to find strong evidence from our genetic study to support the hypothesis that 2D:4D ratio is a direct biomarker of prenatal exposure to androgens in healthy individuals, our findings do not explicitly exclude this possibility, and pathways involving testosterone may become apparent as the size of the discovery sample increases further. Our findings provide new insight into the underlying biology shaping 2D:4D variation in the general population.