Defining disease severity in atopic dermatitis and psoriasis for the application to biomarker research- an inter-disciplinary perspective

More severe atopic dermatitis (AD) and psoriasis are associated with a higher cumulative impact on quality of life, multimorbidity and healthcare costs. Proactive, early intervention in those most at risk of severe disease may reduce this cumulative burden and modify the disease trajectory to limit progression. The lack of reliable biomarkers for this at-risk group represents a barrier to such a paradigm shift in practice. To expedite discovery and validation, the BIOMAP consortium (Biomarkers in AD and Psoriasis, a large-scale European, inter-disciplinary research initiative) has curated clinical and molecular data across diverse study designs and sources including cross-sectional and cohort studies (small scale through to large multi-centre registries), clinical trials, electronic health records and large-scale population-based biobanks. We map all dataset disease severity instruments and measures to three key domains (symptoms, inflammatory activity and disease course), and describe important co-dependencies and relationships across variables and domains. We prioritise definitions for more severe disease with reference to international consensus, reference standards and/or expert opinion. Key factors to consider when analysing datasets across these diverse study types include explicit early consideration of biomarker purpose and clinical context, candidate biomarkers associated with disease severity at a point in time and over time and how they are related, taking the stage of biomarker development into account when selecting disease severity measures for analyses and, validating biomarker associations with disease severity outcomes using both physician- and patient-reported measures and across domains. The outputs from this exercise will ensure coherence and focus across the BIOMAP consortium so that mechanistic insights and biomarkers are clinically relevant, patient-centric and more generalisable to current and future research efforts.

Protein identification for stroke progression via Mendelian Randomization in Million Veteran Program and UK Biobank

Background

Individuals who have experienced a stroke, or transient ischemic attack, face a heightened risk of future cardiovascular events. Identification of genetic and molecular risk factors for subsequent cardiovascular outcomes may identify effective therapeutic targets to improve prognosis after an incident stroke.

Methods

We performed genome-wide association studies (GWAS) for subsequent major adverse cardiovascular events (MACE) (Ncases=51,929, Ncntrl=39,980) and subsequent arterial ischemic stroke (AIS) Ncases=45,120, Ncntrl=46,789) after first incident stroke within the Million Veteran Program and UK Biobank. We then used genetic variants associated with proteins (pQTLs) to determine the effect of 1,463 plasma protein abundances on subsequent MACE using Mendelian randomization (MR).

Results

Two variants were significantly associated with subsequent cardiovascular events: rs76472767 (OR=0.75, 95% CI = 0.64-0.85, p= 3.69×10-08) with subsequent AIS and rs13294166 (OR=1.52, 95% CI = 1.37-1.67, p=3.77×10-08) with subsequent MACE. Using MR, we identified 2 proteins with an effect on subsequent MACE after a stroke: CCL27 (effect OR= 0.77, 95% CI = 0.66-0.88, adj. p=0.05), and TNFRSF14 (effect OR=1.42, 95% CI = 1.24-1.60, adj. p=0.006). These proteins are not associated with incident AIS and are implicated to have a role in inflammation.

Conclusions

We found evidence that two proteins with little effect on incident stroke appear to influence subsequent MACE after incident AIS. These associations suggest that inflammation is a contributing factor to subsequent MACE outcomes after incident AIS and highlights potential novel targets.

Investigation of genetic determinants of cognitive change in later life

Cognitive decline is a major health concern and identification of genes that may serve as drug targets to slow decline is important to adequately support an aging population. Whilst genetic studies of cross-sectional cognition have been carried out, cognitive change is less well-understood. Here, using data from the TOMMORROW trial, we investigate genetic associations with cognitive change in a cognitively normal older cohort. We conducted a genome-wide association study of trajectories of repeated cognitive measures (using generalised estimating equation (GEE) modelling) and tested associations with polygenic risk scores (PRS) of potential risk factors. We identified two genetic variants associated with change in attention domain scores, rs534221751 (p = 1 × 10-8 with slope 1) and rs34743896 (p = 5 × 10-10 with slope 2), implicating NCAM2 and CRIPT/ATP6V1E2 genes, respectively. We also found evidence for the association between an education PRS and baseline cognition (at >65 years of age), particularly in the language domain. We demonstrate the feasibility of conducting GWAS of cognitive change using GEE modelling and our results suggest that there may be novel genetic associations for cognitive change that have not previously been associated with cross-sectional cognition. We also show the importance of the education PRS on cognition much later in life. These findings warrant further investigation and demonstrate the potential value of using trial data and trajectory modelling to identify genetic variants associated with cognitive change.

European and multi-ancestry genome-wide association meta-analysis of atopic dermatitis highlights importance of systemic immune regulation

Atopic dermatitis (AD) is a common inflammatory skin condition and prior genome-wide association studies (GWAS) have identified 71 associated loci. In the current study we conducted the largest AD GWAS to date (discovery N = 1,086,394, replication N = 3,604,027), combining previously reported cohorts with additional available data. We identified 81 loci (29 novel) in the European-only analysis (which all replicated in a separate European analysis) and 10 additional loci in the multi-ancestry analysis (3 novel). Eight variants from the multi-ancestry analysis replicated in at least one of the populations tested (European, Latino or African), while two may be specific to individuals of Japanese ancestry. AD loci showed enrichment for DNAse I hypersensitivity and eQTL associations in blood. At each locus we prioritised candidate genes by integrating multi-omic data. The implicated genes are predominantly in immune pathways of relevance to atopic inflammation and some offer drug repurposing opportunities.

Strategies to investigate and mitigate collider bias in genetic and Mendelian randomisation studies of disease progression

Genetic studies of disease progression can be used to identify factors that may influence survival or prognosis, which may differ from factors that influence on disease susceptibility. Studies of disease progression feed directly into therapeutics for disease, whereas studies of incidence inform prevention strategies. However, studies of disease progression are known to be affected by collider (also known as "index event") bias since the disease progression phenotype can only be observed for individuals who have the disease. This applies equally to observational and genetic studies, including genome-wide association studies and Mendelian randomisation (MR) analyses. In this paper, our aim is to review several statistical methods that can be used to detect and adjust for index event bias in studies of disease progression, and how they apply to genetic and MR studies using both individual- and summary-level data. Methods to detect the presence of index event bias include the use of negative controls, a comparison of associations between risk factors for incidence in individuals with and without the disease, and an inspection of Miami plots. Methods to adjust for the bias include inverse probability weighting (with individual-level data), or Slope-Hunter and Dudbridge et al.'s index event bias adjustment (when only summary-level data are available). We also outline two approaches for sensitivity analysis. We then illustrate how three methods to minimise bias can be used in practice with two applied examples. Our first example investigates the effects of blood lipid traits on mortality from coronary heart disease, while our second example investigates genetic associations with breast cancer mortality.

Epigenetic Regulation of F2RL3 Associates With Myocardial Infarction and Platelet Function

BACKGROUND

DNA hypomethylation at the F2RL3 (F2R like thrombin or trypsin receptor 3) locus has been associated with both smoking and atherosclerotic cardiovascular disease; whether these smoking-related associations form a pathway to disease is unknown. F2RL3 encodes protease-activated receptor 4, a potent thrombin receptor expressed on platelets. Given the role of thrombin in platelet activation and the role of thrombus formation in myocardial infarction, alterations to this biological pathway could be important for ischemic cardiovascular disease.

METHODS

We conducted multiple independent experiments to assess whether DNA hypomethylation at F2RL3 in response to smoking is associated with risk of myocardial infarction via changes to platelet reactivity. Using cohort data (N=3205), we explored the relationship between smoking, DNA hypomethylation at F2RL3, and myocardial infarction. We compared platelet reactivity in individuals with low versus high DNA methylation at F2RL3 (N=41). We used an in vitro model to explore the biological response of F2RL3 to cigarette smoke extract. Finally, a series of reporter constructs were used to investigate how differential methylation could impact F2RL3 gene expression.

RESULTS

Observationally, DNA methylation at F2RL3 mediated an estimated 34% of the smoking effect on increased risk of myocardial infarction. An association between methylation group (low/high) and platelet reactivity was observed in response to PAR4 (protease-activated receptor 4) stimulation. In cells, cigarette smoke extract exposure was associated with a 4.9% to 9.3% reduction in DNA methylation at F2RL3 and a corresponding 1.7-(95% CI, 1.2-2.4, P=0.04) fold increase in F2RL3 mRNA. Results from reporter assays suggest the exon 2 region of F2RL3 may help control gene expression.

CONCLUSIONS

Smoking-induced epigenetic DNA hypomethylation at F2RL3 appears to increase PAR4 expression with potential downstream consequences for platelet reactivity. Combined evidence here not only identifies F2RL3 DNA methylation as a possible contributory pathway from smoking to cardiovascular disease risk but from any feature potentially influencing F2RL3 regulation in a similar manner.

Imputation provides an opportunity to study filaggrin (FLG) null mutations in large population cohorts that lack bespoke genotyping

Background: Low frequency mutations within the filaggrin (FLG) gene are established genetic risk factors for atopic dermatitis. Studies of FLG have typically used sequencing or bespoke genotyping. Large-scale population cohorts with genome-wide imputed data offer powerful genetic analysis opportunities, but bespoke FLG genotyping is often not feasible in such studies. Therefore, we aimed to determine the quality of selected FLG null genotype data extracted from genome-wide imputed sources, focussing on UK population data. Methods: We compared the allele frequencies of three FLG null mutations (R501X, R2447X and S3247X) in directly genotyped and genome-wide imputed data in the ALSPAC cohort. Logistic regression analysis was used to test the association of atopic dermatitis with imputed and genotyped FLG null mutations in ALSPAC and UK Biobank to investigate the usefulness of imputed FLG data. Results: The three FLG null mutations appear to be well imputed in datasets that use the Haplotype Reference Consortium (HRC) for imputation (0.3% discordance compared with directly genotyped data). However, a greater proportion of null alleles failed imputation compared to wild-type alleles. Despite the calling of FLG mutations in imputed data being imperfect, they are still strongly associated with atopic dermatitis (p-values between 7x10-10 and 5x10-75 in UK Biobank). Conclusions: HRC imputed data appears to be adequate for UK population-based genetic analysis of selected FLG null mutations.

Using multivariable Mendelian randomization to estimate the causal effect of bone mineral density on osteoarthritis risk, independently of body mass index

OBJECTIVES

Observational analyses suggest that high bone mineral density (BMD) is a risk factor for osteoarthritis (OA); it is unclear whether this represents a causal effect or shared aetiology and whether these relationships are body mass index (BMI)-independent. We performed bidirectional Mendelian randomization (MR) to uncover the causal pathways between BMD, BMI and OA.

METHODS

One-sample (1S)MR estimates were generated by two-stage least-squares regression. Unweighted allele scores instrumented each exposure. Two-sample (2S)MR estimates were generated using inverse-variance weighted random-effects meta-analysis. Multivariable MR (MVMR), including BMD and BMI instruments in the same model, determined the BMI-independent causal pathway from BMD to OA. Latent causal variable (LCV) analysis, using weight-adjusted femoral neck (FN)-BMD and hip/knee OA summary statistics, determined whether genetic correlation explained the causal effect of BMD on OA.

RESULTS

1SMR provided strong evidence for a causal effect of BMD estimated from heel ultrasound (eBMD) on hip and knee OA {odds ratio [OR]hip = 1.28 [95% confidence interval (CI) = 1.05, 1.57], p = 0.02, ORknee = 1.40 [95% CI = 1.20, 1.63], p = 3 × 10-5, OR per standard deviation [SD] increase}. 2SMR effect sizes were consistent in direction. Results suggested that the causal pathways between eBMD and OA were bidirectional (βhip = 1.10 [95% CI = 0.36, 1.84], p = 0.003, βknee = 4.16 [95% CI = 2.74, 5.57], p = 8 × 10-9, β = SD increase per doubling in risk). MVMR identified a BMI-independent causal pathway between eBMD and hip/knee OA. LCV suggested that genetic correlation (i.e. shared genetic aetiology) did not fully explain the causal effects of BMD on hip/knee OA.

CONCLUSIONS

These results provide evidence for a BMI-independent causal effect of eBMD on OA. Despite evidence of bidirectional effects, the effect of BMD on OA did not appear to be fully explained by shared genetic aetiology, suggesting a direct action of bone on joint deterioration.

Rare variant analysis in eczema identifies exonic variants in DUSP1, NOTCH4 and SLC9A4

Previous genome-wide association studies revealed multiple common variants involved in eczema but the role of rare variants remains to be elucidated. Here, we investigate the role of rare variants in eczema susceptibility. We meta-analyze 21 study populations including 20,016 eczema cases and 380,433 controls. Rare variants are imputed with high accuracy using large population-based reference panels. We identify rare exonic variants in DUSP1, NOTCH4, and SLC9A4 to be associated with eczema. In DUSP1 and NOTCH4 missense variants are predicted to impact conserved functional domains. In addition, five novel common variants at SATB1-AS1/KCNH8, TRIB1/LINC00861, ZBTB1, TBX21/OSBPL7, and CSF2RB are discovered. While genes prioritized based on rare variants are significantly up-regulated in the skin, common variants point to immune cell function. Over 20% of the single nucleotide variant-based heritability is attributable to rare and low-frequency variants. The identified rare/low-frequency variants located in functional protein domains point to promising targets for novel therapeutic approaches to eczema.

Genetic studies of eczema to date have mostly explored common genetic variation. Here, the authors perform a large meta-analysis for common and rare variants and discover 8 loci associated with eczema. Over 20% of the heritability of the condition is attributable to rare variants.

Investigating the causal relationship between allergic disease and mental health

BACKGROUND

Observational studies have reported an association between allergic disease and mental health, but a causal relationship has not been established. Here, we use Mendelian randomization (MR) to investigate a possible causal relationship between atopic disease and mental health phenotypes.

METHODS

The observational relationship between allergic disease and mental health was investigated in UK Biobank. The direction of causality was investigated with bidirectional two-sample MR using summary-level data from published genome-wide association studies. A genetic instrument was derived from associated variants for a broad allergic disease phenotype to test for causal relationships with various mental health outcomes. We also investigated whether these relationships were specific to atopic dermatitis (AD), asthma or hayfever. Given the multiple testing burden, we applied a Bonferroni correction to use an individual test p-value threshold of .0016 (32 tests).

RESULTS

We found strong evidence of an observational association between the broad allergic disease phenotype and depression (OR_self-report =1.45, 95% CI: 1.41-1.50, p = 3.6 × 10^-130 ), anxiety (OR=1.25, 95% CI: 1.18-1.33, p = 6.5 × 10^-13 ), bipolar disorder (OR_self-report =1.29, 95% CI: 1.12-1.47, p = 2.8 × 10^-4 ) and neuroticism (β = 0.38, 95% CI: 0.36-0.41, p = 6.8 × 10^-166 ). Similar associations were found between asthma, AD, hayfever individually with the mental health phenotypes, although the associations between AD and hayfever with bipolar disorder were weaker. There was little evidence of causality in either direction (all p-values>.02).

CONCLUSION

Using MR, we were unable to replicate most of the phenotypic associations between allergic disease and mental health. Any causal effects we detected were considerably attenuated compared with the phenotypic association. This suggests that most comorbidity observed clinically is unlikely to be causal.

Author Correction: Investigating causality in the association between vitamin D status and self-reported tiredness

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

MendelVar: gene prioritization at GWAS loci using phenotypic enrichment of Mendelian disease genes

Motivation

Gene prioritization at human GWAS loci is challenging due to linkage-disequilibrium and long-range gene regulatory mechanisms. However, identifying the causal gene is crucial to enable identification of potential drug targets and better understanding of molecular mechanisms. Mapping GWAS traits to known phenotypically relevant Mendelian disease genes near a locus is a promising approach to gene prioritization.

Results

We present MendelVar, a comprehensive tool that integrates knowledge from four databases on Mendelian disease genes with enrichment testing for a range of associated functional annotations such as Human Phenotype Ontology, Disease Ontology and variants from ClinVar. This open web-based platform enables users to strengthen the case for causal importance of phenotypically matched candidate genes at GWAS loci. We demonstrate the use of MendelVar in post-GWAS gene annotation for type 1 diabetes, type 2 diabetes, blood lipids and atopic dermatitis.

Availability and implementation

MendelVar is freely available at https://mendelvar.mrcieu.ac.uk

Supplementary information

Supplementary data are available at Bioinformatics online.

Triangulating molecular evidence to prioritize candidate causal genes at established atopic dermatitis loci

Genome-wide association studies for atopic dermatitis (AD) have identified 25 reproducible loci. We attempt to prioritize candidate causal genes at these loci using extensive molecular resources compiled into a bioinformatics pipeline. We identified a list of 103 molecular resources for AD aetiology, including expression, protein and DNA methylation QTL datasets in skin or immune-relevant tissues which were tested for overlap with GWAS signals. This was combined with functional annotation using regulatory variant prediction, and features such as promoter-enhancer interactions, expression studies and variant fine-mapping. For each gene at each locus, we condensed the evidence into a prioritization score. Across the investigated loci, we detected significant enrichment of genes with adaptive immune regulatory function and epidermal barrier formation among the top prioritized genes. At 8 loci, we were able to prioritize a single candidate gene (IL6R, ADO, PRR5L, IL7R, ETS1, INPP5D, MDM1, TRAF3). In addition, at 6 of the 25 loci, our analysis prioritizes less familiar candidates (SLC22A5, IL2RA, MDM1, DEXI, ADO, STMN3). Our analysis provides support for previously implicated genes at several AD GWAS loci, as well as evidence for plausible additional candidates at others, which may represent potential targets for drug discovery.

Mendelian randomization provides evidence for a causal effect of higher serum IGF-1 concentration on risk of hip and knee osteoarthritis

Objectives

How insulin-like growth factor-1 (IGF-1) is related to OA is not well understood. We determined relationships between IGF-1 and hospital-diagnosed hand, hip and knee OA in UK Biobank, using Mendelian randomization (MR) to determine causality.

Methods

Serum IGF-1 was assessed by chemiluminescent immunoassay. OA was determined using Hospital Episode Statistics. One-sample MR (1SMR) was performed using two-stage least-squares regression, with an unweighted IGF-1 genetic risk score as an instrument. Multivariable MR included BMI as an additional exposure (instrumented by BMI genetic risk score). MR analyses were adjusted for sex, genotyping chip and principal components. We then performed two-sample MR (2SMR) using summary statistics from Cohorts for Heart and Aging Research in Genetic Epidemiology (CHARGE) (IGF-1, N = 30 884) and the recent genome-wide association study meta-analysis (N = 455 221) of UK Biobank and Arthritis Research UK OA Genetics (arcOGEN).

Results

A total of 332 092 adults in UK Biobank had complete data. Their mean (s.d.) age was 56.5 (8.0) years and 54% were female. IGF-1 was observationally related to a reduced odds of hand OA [odds ratio per doubling = 0.87 (95% CI 0.82, 0.93)], and an increased odds of hip OA [1.04 (1.01, 1.07)], but was unrelated to knee OA [0.99 (0.96, 1.01)]. Using 1SMR, we found strong evidence for an increased risk of hip [odds ratio per s.d. increase = 1.57 (1.21, 2.01)] and knee [1.30 (1.07, 1.58)] OA with increasing IGF-1 concentration. By contrast, we found no evidence for a causal effect of IGF-1 concentration on hand OA [0.98 (0.57, 1.70)]. Results were consistent when estimated using 2SMR and in multivariable MR analyses accounting for BMI.

Conclusion

We have found evidence that increased serum IGF-1 is causally related to higher risk of hip and knee OA.

Osteoarthritis: Insights Offered by the Study of Bone Mass Genetics

PURPOSE OF REVIEW

This paper reviews how bone genetics has contributed to our understanding of the pathogenesis of osteoarthritis. As well as identifying specific genetic mechanisms involved in osteoporosis which also contribute to osteoarthritis, we review whether bone mineral density (BMD) plays a causal role in OA development.

RECENT FINDINGS

We examined whether those genetically predisposed to elevated BMD are at increased risk of developing OA, using our high bone mass (HBM) cohort. HBM individuals were found to have a greater prevalence of OA compared with family controls and greater development of radiographic features of OA over 8 years, with predominantly osteophytic OA. Initial Mendelian randomisation analysis provided additional support for a causal effect of increased BMD on increased OA risk. In contrast, more recent investigation estimates this relationship to be bi-directional. However, both these findings could be explained instead by shared biological pathways. Pathways which contribute to BMD appear to play an important role in OA development, likely reflecting shared common mechanisms as opposed to a causal effect of raised BMD on OA. Studies in HBM individuals suggest this reflects an important role of mechanisms involved in bone formation in OA development; however further work is required to establish whether the same applies to more common forms of OA within the general population.

The influence of adult hip shape genetic variants on adolescent hip shape: Findings from a population-based DXA study

OBJECTIVE

Hip shape is a well-recognized risk factor for hip osteoarthritis (OA) and hip fracture. We aimed to investigate whether the genetic variants known to be associated with adult hip shape were also associated with adolescent hip shape.

METHODS

Hip DXA scans, obtained in offspring from the Avon Longitudinal Study of Parents and Children (ALSPAC) at two time points (mean ages 13.8 and 17.8 years), were used to quantify hip morphology using a 53-point Statistical Shape Model (SSM). Principal component analysis was used to generate hip shape modes (HSMs). Genetic variants which had previously shown genome-wide significant association with specific HSMs in adults were tested for association with the same HSMs in adolescents (at each timepoint separately) using SNPTEST v2.

RESULTS

Complete genotypic and phenotypic data were available for 3550 and 3175 individuals at 14 and 18 years, respectively. The strongest evidence for association with adolescent hip shape was for a variant located near SOX9 (rs2158915) with consistent effects across both time points for HSM1 (age 14: beta -0.05, p = 9.9 × 10-8; age 18: -0.05, p = 3.3 × 10-6) and HSM5 (age 14: beta -0.07, p = 1.6 × 10-4; age 18: -0.1, p = 2.7 × 10-6). There was also strong evidence of association between rs10743612 (near PTHLH) and HSM1 (age 14: 0.05, p = 1.1 × 10-5; age 18: 0.04, p = 0.003) and between rs6537291 (near HHIP) and HSM2 (age 14: -0.06, p = 0.001; age 18: -0.07, p = 0.001) across both time points. The genes with the strongest associations with hip shape in adolescents, (SOX9, PTHLH and HHIP) are known to be involved in endochondral bone formation. HSM1 indicates narrower aspect ratio of the upper femur, whereas both HSM2 and HSM5 reflect variation in the femoral head size and femoral neck width, features previously found to be related to the risk of OA in later life. The SOX9 locus has previously been found to associate with increased risk of hip fracture.

CONCLUSION

In conclusion, variants implicated in endochondral bone formation appear to consistently influence hip shape between adolescence and adulthood, including those aspects related to risk of hip OA and/or fracture in later life.

Increased development of radiographic hip osteoarthritis in individuals with high bone mass: a prospective cohort study

BACKGROUND

Individuals with high bone mass (HBM) have a greater odds of prevalent radiographic hip osteoarthritis (OA), reflecting an association with bone-forming OA sub-phenotypes (e.g. osteophytosis, subchondral sclerosis). As the role of bone mineral density (BMD) in hip OA progression is unclear, we aimed to determine if individuals with HBM have increased incidence and/or progression of bone-forming OA sub-phenotypes.

METHODS

We analysed an adult cohort with and without HBM (L1 and/or total hip BMD Z-score > + 3.2) with pelvic radiographs collected at baseline and 8-year follow-up. Sub-phenotypes were graded using the OARSI atlas. Superior/inferior acetabular/femoral osteophyte and medial/superior joint space narrowing (JSN) grades were summed and Δosteophyte and ΔJSN derived. Pain and functional limitations were quantified using the WOMAC questionnaire. Associations between HBM status and change in OA sub-phenotypes were determined using multivariable linear/logistic regression, adjusting for age, sex, height, total body fat mass, follow-up time and baseline sub-phenotype grade. Generalised estimating equations accounted for individual-level clustering.

RESULTS

Of 136 individuals, 62% had HBM at baseline, 72% were female and mean (SD) age was 59 (10) years. HBM was positively associated with both Δosteophytes and ΔJSN (adjusted mean grade differences between individuals with and without HBM β_osteophyte = 0.30 [0.01, 0.58], p = 0.019 and β_JSN = 0.10 [0.01, 0.18], p = 0.019). Incident subchondral sclerosis was rare. HBM individuals had higher WOMAC hip functional limitation scores (β = 8.3 [0.7, 15.98], p = 0.032).

CONCLUSIONS

HBM is associated with the worsening of hip osteophytes and JSN over an average of 8 years, as well as increased hip pain and functional limitation.

Genome-wide association study identifies 48 common genetic variants associated with handedness

Handedness has been extensively studied because of its relationship with language and the over-representation of left-handers in some neurodevelopmental disorders. Using data from the UK Biobank, 23andMe and the International Handedness Consortium, we conducted a genome-wide association meta-analysis of handedness (N = 1,766,671). We found 41 loci associated (P < 5 × 10-8) with left-handedness and 7 associated with ambidexterity. Tissue-enrichment analysis implicated the CNS in the aetiology of handedness. Pathways including regulation of microtubules and brain morphology were also highlighted. We found suggestive positive genetic correlations between left-handedness and neuropsychiatric traits, including schizophrenia and bipolar disorder. Furthermore, the genetic correlation between left-handedness and ambidexterity is low (rG = 0.26), which implies that these traits are largely influenced by different genetic mechanisms. Our findings suggest that handedness is highly polygenic and that the genetic variants that predispose to left-handedness may underlie part of the association with some psychiatric disorders.

Genetic Analyses of Common Infections in the Avon Longitudinal Study of Parents and Children Cohort

The burden of infections on an individual and public health is profound. Many observational studies have shown a link between infections and the pathogenesis of disease; however a greater understanding of the role of host genetics is essential. Children from the longitudinal birth cohort, the Avon Longitudinal Study of Parents and Children, had 14 antibodies measured in plasma at age 7: Alpha-casein protein, beta-casein protein, cytomegalovirus, Epstein-Barr virus, feline herpes virus, Helicobacter pylori, herpes simplex virus 1, influenza virus subtype H1N1, influenza virus subtype H3N2, measles virus, Saccharomyces cerevisiae, Theiler's virus, Toxoplasma gondii, and SAG1 protein domain, a surface antigen of Toxoplasma gondii measured for greater precision. We performed genome-wide association analyses of antibody levels against these 14 infections (N = 357 - 5010) and identified three genome-wide signals (P < 5×10^-8), two associated with measles virus antibodies and one with Toxoplasma gondii antibodies. In an association analysis focused on the human leukocyte antigen (HLA) region of the genome, we further detected 15 HLA alleles at a two-digit resolution and 23 HLA alleles at a four-digit resolution associated with five antibodies, with eight HLA alleles associated with Epstein-Barr virus antibodies showing strong evidence of replication in UK Biobank. We discuss how our findings from antibody levels complement other studies using self-reported phenotypes in understanding the architecture of host genetics related to infections.

Individuals with high bone mass have increased progression of radiographic and clinical features of knee osteoarthritis

OBJECTIVE

High bone mass (HBM) is associated with an increased prevalence of radiographic knee OA (kOA), characterized by osteophytosis. We aimed to determine if progression of radiographic kOA, and its sub-phenotypes, is increased in HBM and whether observed changes are clinically relevant.

DESIGN

A cohort with and without HBM (L1 and/or total hip bone mineral density Z-score≥+3.2) had knee radiographs collected at baseline and 8-year follow-up. Sub-phenotypes were graded using the OARSI atlas. Medial/lateral tibial/femoral osteophyte and medial/lateral joint space narrowing (JSN) grades were summed and Δosteophytes, ΔJSN derived. Pain, function and stiffness were quantified using the WOMAC questionnaire. Associations between HBM status and sub-phenotype progression were determined using multivariable linear/poisson regression, adjusting for age, sex, height, baseline sub-phenotype grade, menopause, education and total body fat mass (TBFM). Generalized estimating equations accounted for individual-level clustering.

RESULTS

169 individuals had repeated radiographs, providing 330 knee images; 63% had HBM, 73% were female, mean (SD) age was 58 (12) years. Whilst HBM was not clearly associated with overall Kellgren-Lawrence measured progression (RR = 1.55 [0.56.4.32]), HBM was positively associated with both Δosteophytes and ΔJSN individually (adjusted mean differences between individuals with and without HBM 0.45 [0.01.0.89] and 0.15 [0.01.0.29], respectively). HBM individuals had higher WOMAC knee pain scores (β = 7.42 [1.17.13.66]), largely explained by adjustment for osteophyte score (58% attenuated) rather than JSN (30% attenuated) or TBFM (16% attenuated). The same pattern was observed for symptomatic stiffness and functional limitation.

CONCLUSIONS

HBM is associated with osteophyte progression, which appears to contribute to increased reported pain, stiffness and functional loss.

What is the evidence for interactions between filaggrin null mutations and environmental exposures in the aetiology of atopic dermatitis? A systematic review

BACKGROUND

Epidemiological studies indicate that gene-environment interactions play a role in atopic dermatitis (AD).

OBJECTIVES

To review the evidence for gene-environment interactions in AD aetiology, focusing on filaggrin (FLG) loss-of-function mutations.

METHODS

A systematic search from inception to September 2018 in Embase, MEDLINE and BIOSIS was performed. Search terms included all synonyms for AD and filaggrin/FLG; any genetic or epidemiological study design using any statistical methods were included. Quality assessment using criteria modified from guidance (ROBINS-I and Human Genome Epidemiology Network) for nonrandomized and genetic studies was completed, including consideration of power. Heterogeneity of study design and analyses precluded the use of meta-analysis.

RESULTS

Of 1817 papers identified, 12 studies fulfilled the inclusion criteria required and performed formal interaction testing. There was some evidence for FLG-environment interactions in six of the studies (P-value for interaction ≤ 0·05), including early-life cat ownership, older siblings, water hardness, phthalate exposure, higher urinary phthalate metabolite levels (which all increased AD risk additional to FLG null genotype) and prolonged breastfeeding (which decreased AD risk in the context of FLG null genotype). Major limitations of published studies were the low numbers of individuals (ranging from five to 94) with AD and FLG loss-of-function mutations and exposure to specific environmental factors, and variation in exposure definitions.

CONCLUSIONS

Evidence on FLG-environment interactions in AD aetiology is limited. However, many of the studies lacked large enough sample sizes to assess these interactions fully. Further research is needed with larger sample sizes and clearly defined exposure assessment. Linked Comment: Park and Seo. Br J Dermatol 2020; 183:411.

Mendelian randomization analysis does not support causal associations of birth weight with hypertension risk and blood pressure in adulthood

Epidemiology studies suggested that low birthweight was associated with a higher risk of hypertension in later life. However, little is known about the causality of such associations. In our study, we evaluated the causal association of low birthweight with adulthood hypertension following a standard analytic protocol using the study-level data of 183,433 participants from 60 studies (CHARGE-BIG consortium), as well as that with blood pressure using publicly available summary-level genome-wide association data from EGG consortium of 153,781 participants, ICBP consortium and UK Biobank cohort together of 757,601 participants. We used seven SNPs as the instrumental variable in the study-level analysis and 47 SNPs in the summary-level analysis. In the study-level analyses, decreased birthweight was associated with a higher risk of hypertension in adults (the odds ratio per 1 standard deviation (SD) lower birthweight, 1.22; 95% CI 1.16 to 1.28), while no association was found between genetically instrumented birthweight and hypertension risk (instrumental odds ratio for causal effect per 1 SD lower birthweight, 0.97; 95% CI 0.68 to 1.41). Such results were consistent with that from the summary-level analyses, where the genetically determined low birthweight was not associated with blood pressure measurements either. One SD lower genetically determined birthweight was not associated with systolic blood pressure (β = - 0.76, 95% CI - 2.45 to 1.08 mmHg), 0.06 mmHg lower diastolic blood pressure (β = - 0.06, 95% CI - 0.93 to 0.87 mmHg), or pulse pressure (β = - 0.65, 95% CI - 1.38 to 0.69 mmHg, all p > 0.05). Our findings suggest that the inverse association of birthweight with hypertension risk from observational studies was not supported by large Mendelian randomization analyses.

The effect of pubertal timing, as reflected by height tempo, on proximal femur shape: Findings from a population-based study in adolescents

OBJECTIVE

To examine the relationship between pubertal timing (using measures of height tempo) and proximal femur shape in a large adolescent cohort.

METHODS

Hip DXA scans were obtained in offspring from the Avon Longitudinal Study of Parents and Children. To quantify hip morphology, the images were analyzed using Shape software based on a 53-point statistical shape model and independent modes of variation (hip shape mode (HSM) scores) for each image were generated. Height tempo (which corresponds to age at peak height velocity (aPHV)) was estimated from serial height measurements collected between age 5-20 years. Multivariable linear regression was used to examine cross-sectional associations between height tempo and the top ten HSMs at age 14 and 18, adjusting for sex and fat mass index (FMI).

RESULTS

Complete outcome and covariate data were available from 3827 and 3507 participants at age 14 and 18 years, respectively. Mean aPHV was 13.5 and 11.8 years for males and females, respectively. At age 14, height tempo was associated with a majority of modes, except for HSM4 and there was strong evidence of interaction by sex. In males, all modes showed evidence of an association with tempo, independent of FMI, with the strongest observed for HSM8 (adjusted β 0.38 (0.33, 0.43) p = 4.1 × 10^-50). Compared with males, the associations were generally weaker in females, with the strongest effect observed for HSM8 (adjusted β 0.10 (0.05, 0.14) p = 1.6 × 10^-5). The overall effect of later pubertal timing on proximal femur shape in males was a narrower femoral neck and larger superolateral head, whereas in females these changes were hard to discern. When assessed at age 18, there was little relationship between tempo and proximal femur shape in either sex.

CONCLUSION

Our results indicate that significant changes in hip shape occur during puberty, including aspects of shape which may be related to future risk of hip OA and/or fracture. However, puberty timing per se does not appear to exert long lasting effects on proximal femur shape.

A Rare Mutation in SMAD9 Associated With High Bone Mass Identifies the SMAD-Dependent BMP Signaling Pathway as a Potential Anabolic Target for Osteoporosis

Novel anabolic drug targets are needed to treat osteoporosis. Having established a large national cohort with unexplained high bone mass (HBM), we aimed to identify a novel monogenic cause of HBM and provide insight into a regulatory pathway potentially amenable to therapeutic intervention. We investigated a pedigree with unexplained HBM in whom previous sequencing had excluded known causes of monogenic HBM. Whole exome sequencing identified a rare (minor allele frequency 0.0023), highly evolutionarily conserved missense mutation in SMAD9 (c.65T>C, p.Leu22Pro) segregating with HBM in this autosomal dominant family. The same mutation was identified in another two unrelated individuals both with HBM. In silico protein modeling predicts the mutation severely disrupts the MH1 DNA-binding domain of SMAD9. Affected individuals have bone mineral density (BMD) Z-scores +3 to +5, mandible enlargement, a broad frame, torus palatinus/mandibularis, pes planus, increased shoe size, and a tendency to sink when swimming. Peripheral quantitative computed tomography (pQCT) measurement demonstrates increased trabecular volumetric BMD and increased cortical thickness conferring greater predicted bone strength; bone turnover markers are low/normal. Notably, fractures and nerve compression are not found. Both genome-wide and gene-based association testing involving estimated BMD measured at the heel in 362,924 white British subjects from the UK Biobank Study showed strong associations with SMAD9 (P_GWAS = 6 × 10^-16 ; P_GENE = 8 × 10^-17 ). Furthermore, we found Smad9 to be highly expressed in both murine cortical bone-derived osteocytes and skeletal elements of zebrafish larvae. Our findings support SMAD9 as a novel HBM gene and a potential novel osteoanabolic target for osteoporosis therapeutics. SMAD9 is thought to inhibit bone morphogenetic protein (BMP)-dependent target gene transcription to reduce osteoblast activity. Thus, we hypothesize SMAD9 c.65T>C is a loss-of-function mutation reducing BMP inhibition. Lowering SMAD9 as a potential novel anabolic mechanism for osteoporosis therapeutics warrants further investigation. © 2019 The Authors. Journal of Bone and Mineral Research published by American Society for Bone and Mineral Research.

Metabolomics analysis in adults with high bone mass identifies a relationship between bone resorption and circulating citrate which replicates in the general population

OBJECTIVE

Bone turnover, which regulates bone mass, may exert metabolic consequences, particularly on markers of glucose metabolism and adiposity. To better understand these relationships, we examined cross-sectional associations between bone turnover markers (BTMs) and metabolic traits in a population with high bone mass (HBM, BMD Z-score ≥+3.2).

DESIGN

β-C-terminal telopeptide of type-I collagen (β-CTX), procollagen type-1 amino-terminal propeptide (P1NP) and osteocalcin were assessed by electrochemiluminescence immunoassays. Metabolic traits, including lipids and glycolysis-related metabolites, were measured using nuclear magnetic resonance spectroscopy. Associations of BTMs with metabolic traits were assessed using generalized estimating equation linear regression, accounting for within-family correlation, adjusting for potential confounders (age, sex, height, weight, menopause, bisphosphonate and oral glucocorticoid use).

RESULTS

A total of 198 adults with HBM had complete data, mean [SD] age 61.6 [13.7] years; 77% were female. Of 23 summary metabolic traits, citrate was positively related to all BTMs: adjusted ββ-CTX  = 0.050 (95% CI 0.024, 0.076), P = 1.71 × 10-4 , βosteocalcin  = 6.54 × 10-4 (1.87 × 10-4 , 0.001), P = .006 and βP1NP  = 2.40 × 10-4 (6.49 × 10-5 , 4.14 × 10-4 ), P = .007 (β = increase in citrate (mmol/L) per 1 µg/L BTM increase). Inverse relationships of β-CTX (β = -0.276 [-0.434, -0.118], P = 6.03 × 10-4 ) and osteocalcin (-0.004 [-0.007, -0.001], P = .020) with triglycerides were also identified. We explored the generalizability of these associations in 3664 perimenopausal women (age 47.9 [4.4] years) from a UK family cohort. We confirmed a positive, albeit lower magnitude, association between β-CTX and citrate (adjusted βwomen  = 0.020 [0.013, 0.026], P = 1.95 × 10-9 ) and an inverse association of similar magnitude between β-CTX and triglycerides (β = -0.354 [-0.471, -0.237], P = 3.03 × 10-9 ).

CONCLUSIONS

Bone resorption is positively related to circulating citrate and inversely related to triglycerides. Further studies are justified to determine whether plasma citrate or triglyceride concentrations are altered by factors known to modulate bone resorption, such as bisphosphonates.

Novel genetic loci affecting facial shape variation in humans

The human face represents a combined set of highly heritable phenotypes, but knowledge on its genetic architecture remains limited, despite the relevance for various fields. A series of genome-wide association studies on 78 facial shape phenotypes quantified from 3-dimensional facial images of 10,115 Europeans identified 24 genetic loci reaching study-wide suggestive association (p < 5 × 10-8), among which 17 were previously unreported. A follow-up multi-ethnic study in additional 7917 individuals confirmed 10 loci including six unreported ones (padjusted < 2.1 × 10-3). A global map of derived polygenic face scores assembled facial features in major continental groups consistent with anthropological knowledge. Analyses of epigenomic datasets from cranial neural crest cells revealed abundant cis-regulatory activities at the face-associated genetic loci. Luciferase reporter assays in neural crest progenitor cells highlighted enhancer activities of several face-associated DNA variants. These results substantially advance our understanding of the genetic basis underlying human facial variation and provide candidates for future in-vivo functional studies.

Protein-coding variants contribute to the risk of atopic dermatitis and skin-specific gene expression

BACKGROUND

Fifteen percent of atopic dermatitis (AD) liability-scale heritability could be attributed to 31 susceptibility loci identified by using genome-wide association studies, with only 3 of them (IL13, IL-6 receptor [IL6R], and filaggrin [FLG]) resolved to protein-coding variants.

OBJECTIVE

We examined whether a significant portion of unexplained AD heritability is further explained by low-frequency and rare variants in the gene-coding sequence.

METHODS

We evaluated common, low-frequency, and rare protein-coding variants using exome chip and replication genotype data of 15,574 patients and 377,839 control subjects combined with whole-transcriptome data on lesional, nonlesional, and healthy skin samples of 27 patients and 38 control subjects.

RESULTS

An additional 12.56% (SE, 0.74%) of AD heritability is explained by rare protein-coding variation. We identified docking protein 2 (DOK2) and CD200 receptor 1 (CD200R1) as novel genome-wide significant susceptibility genes. Rare coding variants associated with AD are further enriched in 5 genes (IL-4 receptor [IL4R], IL13, Janus kinase 1 [JAK1], JAK2, and tyrosine kinase 2 [TYK2]) of the IL13 pathway, all of which are targets for novel systemic AD therapeutics. Multiomics-based network and RNA sequencing analysis revealed DOK2 as a central hub interacting with, among others, CD200R1, IL6R, and signal transducer and activator of transcription 3 (STAT3). Multitissue gene expression profile analysis for 53 tissue types from the Genotype-Tissue Expression project showed that disease-associated protein-coding variants exert their greatest effect in skin tissues.

CONCLUSION

Our discoveries highlight a major role of rare coding variants in AD acting independently of common variants. Further extensive functional studies are required to detect all potential causal variants and to specify the contribution of the novel susceptibility genes DOK2 and CD200R1 to overall disease susceptibility.

Differential associations of allergic disease genetic variants with developmental profiles of eczema, wheeze and rhinitis

BACKGROUND

Allergic diseases (eczema, wheeze and rhinitis) in children often present as heterogeneous phenotypes. Understanding genetic associations of specific patterns of symptoms might facilitate understanding of the underlying biological mechanisms.

OBJECTIVE

To examine associations between allergic disease-related variants identified in a recent genome-wide association study and latent classes of allergic diseases (LCADs) in two population-based birth cohorts.

METHODS

Eight previously defined LCADs between birth and 11 years: "No disease," "Atopic march," "Persistent eczema and wheeze," "Persistent eczema with later-onset rhinitis," "Persistent wheeze with later-onset rhinitis," "Transient wheeze," "Eczema only" and "Rhinitis only" were used as the study outcome. Weighted multinomial logistic regression was used to estimate associations between 135 SNPs (and a polygenic risk score, PRS) and LCADs among 6345 individuals from The Avon Longitudinal Study of Parents and Children (ALSPAC). Heterogeneity across LCADs was assessed before and after Bonferroni correction. Results were replicated in Manchester Asthma and Allergy Study (MAAS) (n = 896) and pooled in a meta-analysis.

RESULTS

We found strong evidence for differential genetic associations across the LCADs; pooled PRS heterogeneity P-value = 3.3 × 10-14 , excluding "no disease" class. The associations between the PRS and LCADs in MAAS were remarkably similar to ALSPAC. Two SNPs (a protein-truncating variant in FLG and a SNP within an intron of GSDMB) had evidence for differential association (pooled P-values ≤ 0.006). The FLG locus was differentially associated across LCADs that included eczema, with stronger associations for LCADs with comorbid wheeze and rhinitis. The GSDMB locus in contrast was equally associated across LCADs that included wheeze.

CONCLUSIONS AND CLINICAL RELEVANCE

We have shown complex, but distinct patterns of genetic associations with LCADs, suggesting that heterogeneous mechanisms underlie individual disease trajectories. Establishing the combination of allergic diseases with which each genetic variant is associated may inform therapeutic development and/or predictive modelling.

Describing the application of statistical shape modelling to DXA images to quantify the shape of the proximal femur at ages 14 and 18 years in the Avon Longitudinal Study of Parents and Children

Bones are complex objects with considerable variation in the shape and structure often attributed to anatomical, environmental or genetic differences. In addition, bone shape has been of interest in relation to its associations with disease processes. Hip shape is an important determinant of hip osteoarthritis and osteoporotic hip fracture; however, its quantification is difficult. While previous studies largely focused on individual geometrical indices of hip geometry such as neck-shaft angle or femoral neck width, statistical shape modelling offers the means to quantify the entire contour of the proximal femur, including lesser trochanter and acetabular eyebrow. We describe the derivation of independent modes of variation (hip shape mode scores) to characterise variation in hip shape from dual-energy X-ray absorptiometry (DXA) images in the Avon Longitudinal Study of Parents and Children (ALSPAC) offspring, using statistical shape modelling. ALSPAC is a rich source of phenotypic and genotypic data which provides a unique opportunity to investigate the environmental and genetic influences on hip shape in adolescence, as well as comparison with adult hip shape.

Clinical onset of atopic eczema: Results from 2 nationally representative British birth cohorts followed through midlife

Background

Atopic eczema onset is described primarily in early childhood, and the frequency and characteristics of adult-onset disease remain controversial.

Objective

We sought to determine the proportion of subjects who report atopic eczema symptoms between birth and midadulthood and to examine demographic, immunologic, and genetic factors associated with period of symptom onset.

Methods

We conducted a longitudinal study using data from 2 nationally representative community-based birth cohorts from the United Kingdom: the British Cohort Studies 1958 and 1970. Subjects were followed from birth through age 42 to 50 years. The primary outcome was the age period of self-reported atopic eczema symptom onset based on repeated measures of self-reported atopic eczema at each survey wave.

Results

The annual period prevalence of atopic eczema ranged from 5% to 15% in 2 cohorts of more than 17,000 participants each followed from birth through middle age. There was no clear trend in prevalence by age, and among adults reporting active atopic eczema during a given year, only 38% had symptom onset reported in childhood. When compared with subjects whose eczema started in childhood, those with adult-onset disease were more likely to be women, from Scotland or Northern England, of lower childhood socioeconomic group, smokers in adulthood, and less likely to have a history of asthma. In a subanalysis using data from the 1958 cohort only, genetic mutations previously associated with atopic eczema, including filaggrin-null mutations, and allergen-specific IgE were more common among those with childhood-onset disease.

Conclusion

Rates of self-reported atopic eczema remain high after childhood, and adult-onset atopic eczema has different risk factor associations than childhood-onset eczema.

EMSY expression affects multiple components of the skin barrier with relevance to atopic dermatitis

Background

Atopic dermatitis (AD) is a common, complex, and highly heritable inflammatory skin disease. Genome-wide association studies offer opportunities to identify molecular targets for drug development. A risk locus on chromosome 11q13.5 lies between 2 candidate genes, EMSY and LRRC32 (leucine-rich repeat-containing 32) but the functional mechanisms affecting risk of AD remain unclear.

Objectives

We sought to apply a combination of genomic and molecular analytic techniques to investigate which genes are responsible for genetic risk at this locus and to define mechanisms contributing to atopic skin disease.

Methods

We used interrogation of available genomic and chromosome conformation data in keratinocytes, small interfering RNA (siRNA)–mediated knockdown in skin organotypic culture and functional assessment of barrier parameters, mass spectrometric global proteomic analysis and quantitative lipid analysis, electron microscopy of organotypic skin, and immunohistochemistry of human skin samples.

Results

Genomic data indicate active promoters in the genome-wide association study locus and upstream of EMSY; EMSY, LRRC32, and intergenic variants all appear to be within a single topologically associating domain. siRNA-knockdown of EMSY in organotypic culture leads to enhanced development of barrier function, reflecting increased expression of structural and functional proteins, including filaggrin and filaggrin-2, as well as long-chain ceramides. Conversely, overexpression of EMSY in keratinocytes leads to a reduction in markers of barrier formation. Skin biopsy samples from patients with AD show greater EMSY staining in the nucleus, which is consistent with an increased functional effect of this transcriptional control protein.

Conclusion

Our findings demonstrate an important role for EMSY in transcriptional regulation and skin barrier formation, supporting EMSY inhibition as a therapeutic approach.

Bi-allelic Loss-of-Function CACNA1B Mutations in Progressive Epilepsy-Dyskinesia

The occurrence of non-epileptic hyperkinetic movements in the context of developmental epileptic encephalopathies is an increasingly recognized phenomenon. Identification of causative mutations provides an important insight into common pathogenic mechanisms that cause both seizures and abnormal motor control. We report bi-allelic loss-of-function CACNA1B variants in six children from three unrelated families whose affected members present with a complex and progressive neurological syndrome. All affected individuals presented with epileptic encephalopathy, severe neurodevelopmental delay (often with regression), and a hyperkinetic movement disorder. Additional neurological features included postnatal microcephaly and hypotonia. Five children died in childhood or adolescence (mean age of death: 9 years), mainly as a result of secondary respiratory complications. CACNA1B encodes the pore-forming subunit of the pre-synaptic neuronal voltage-gated calcium channel Ca_v2.2/N-type, crucial for SNARE-mediated neurotransmission, particularly in the early postnatal period. Bi-allelic loss-of-function variants in CACNA1B are predicted to cause disruption of Ca²⁺ influx, leading to impaired synaptic neurotransmission. The resultant effect on neuronal function is likely to be important in the development of involuntary movements and epilepsy. Overall, our findings provide further evidence for the key role of Ca_v2.2 in normal human neurodevelopment.

Maternal and fetal genetic effects on birth weight and their relevance to cardio-metabolic risk factors

Birth weight variation is influenced by fetal and maternal genetic and non-genetic factors, and has been reproducibly associated with future cardio-metabolic health outcomes. In expanded genome-wide association analyses of own birth weight (n = 321,223) and offspring birth weight (n = 230,069 mothers), we identified 190 independent association signals (129 of which are novel). We used structural equation modeling to decompose the contributions of direct fetal and indirect maternal genetic effects, then applied Mendelian randomization to illuminate causal pathways. For example, both indirect maternal and direct fetal genetic effects drive the observational relationship between lower birth weight and higher later blood pressure: maternal blood pressure-raising alleles reduce offspring birth weight, but only direct fetal effects of these alleles, once inherited, increase later offspring blood pressure. Using maternal birth weight-lowering genotypes to proxy for an adverse intrauterine environment provided no evidence that it causally raises offspring blood pressure, indicating that the inverse birth weight-blood pressure association is attributable to genetic effects, and not to intrauterine programming.

Genome wide analysis for mouth ulcers identifies associations at immune regulatory loci

Mouth ulcers are the most common ulcerative condition and encompass several clinical diagnoses, including recurrent aphthous stomatitis (RAS). Despite previous evidence for heritability, it is not clear which specific genetic loci are implicated in RAS. In this genome-wide association study (n = 461,106) heritability is estimated at 8.2% (95% CI: 6.4%, 9.9%). This study finds 97 variants which alter the odds of developing non-specific mouth ulcers and replicate these in an independent cohort (n = 355,744) (lead variant after meta-analysis: rs76830965, near IL12A, OR 0.72 (95% CI: 0.71, 0.73); P = 4.4e−483). Additional effect estimates from three independent cohorts with more specific phenotyping and specific study characteristics support many of these findings. In silico functional analyses provide evidence for a role of T cell regulation in the aetiology of mouth ulcers. These results provide novel insight into the pathogenesis of a common, important condition.

Oral ulcerations are sores of the mucous membrane of the mouth and highly prevalent in the population. Here, in a genome-wide association study, the authors identify 97 loci associated with mouth ulcers highlighting genes involved in T cell-mediated immunity and T_H1 responses.

Investigating causality in the association between vitamin D status and self-reported tiredness

Self-reported tiredness or low energy, often referred to as fatigue, has been linked to low levels of circulating 25-hydroxyvitamin D (25OHD), a biomarker of vitamin D status. Although it is uncertain if the association is causal, fatigue is a common indication for testing, and correcting, low 25OHD-levels. We used two-sample Mendelian randomization to test for genetic evidence of a causal association between low 25OHD-levels and fatigue. Genetic-25OHD associations were estimated from the largest genome-wide association study of vitamin D to date, and genetic-fatigue associations were estimated in 327,478 individuals of European descent in UK Biobank, of whom 19,526 (5.96%) reported fatigue (tiredness or low energy nearly every day over the past two weeks). Using seven genome-wide significant 25OHD-reducing genetic variants, there was little evidence for a causal effect of 25OHD on fatigue (odds ratio for fatigue was 1.05 with 95% confidence interval of 0.87–1.27 per 1-SD decrease in log-transformed 25OHD). There was also little evidence of association between any individual 25OHD-reducing variant and fatigue. Our results suggest that a clinically relevant protective effect of 25OHD-levels on fatigue is unlikely. Therefore, vitamin D supplementation of the general population to raise 25OHD-levels is not likely to be useful in preventing fatigue.

Using statistical shape modelling of DXA images to quantify the shape of the proximal femur at ages 14 and 18 years in the Avon Longitudinal Study of Parents and Children

Hip shape is an important determinant of hip osteoarthritis and osteoporotic hip fracture; however, little is known about its development in childhood and adolescence. While previous studies largely focused on individual geometrical indices of hip geometry such as neck-shaft angle or femoral neck width, statistical shape modelling offers the means to quantify the entire contour of the proximal femur, including lesser trochanter and acetabular eyebrow. We describe the derivation of independent modes of variation (hip shape mode scores) to characterise variation in hip shape from dual-energy X-ray absorptiometry (DXA) images in the Avon Longitudinal Study of Parents and Children (ALSPAC) offspring, using statistical shape modelling. ALSPAC is a rich source of phenotypic and genotypic data which provides a unique opportunity to investigate the environmental and genetic influences on hip shape in adolescence, as well as comparison with adult hip shape.

Identification of Novel Loci Associated With Hip Shape: A Meta-Analysis of Genomewide Association Studies

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 < 5 × 10-9 , adjusted for 10 independent outcomes) single-nucleotide polymorphisms (SNPs) were associated with HSM1, and three SNPs with HSM2. One SNP, in high linkage disequilibrium with rs2158915 associated with HSM1, was associated with HSM5 at genomewide significance. In a look-up of previous GWASs, three of the identified SNPs were associated with hip osteoarthritis, one with hip fracture, and five with height. Seven SNPs were within 200 kb of genes involved in endochondral bone formation, namely SOX9, PTHrP, RUNX1, NKX3-2, FGFR4, DICER1, and HHIP. The SNP adjacent to DICER1 also showed osteoblast cis-regulatory activity of GSC, in which mutations have previously been reported to cause hip dysplasia. For three of the lead SNPs, SNPs in high LD (r2  > 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.

Eleven loci with new reproducible genetic associations with allergic disease risk

BACKGROUND

A recent genome-wide association study (GWAS) identified 99 loci that contain genetic risk variants shared between asthma, hay fever, and eczema. Many more risk loci shared between these common allergic diseases remain to be discovered, which could point to new therapeutic opportunities.

OBJECTIVE

We sought to identify novel risk loci shared between asthma, hay fever, and eczema by applying a gene-based test of association to results from a published GWAS that included data from 360,838 subjects.

METHODS

We used approximate conditional analysis to adjust the results from the published GWAS for the effects of the top risk variants identified in that study. We then analyzed the adjusted GWAS results with the EUGENE gene-based approach, which combines evidence for association with disease risk across regulatory variants identified in different tissues. Novel gene-based associations were followed up in an independent sample of 233,898 subjects from the UK Biobank study.

RESULTS

Of the 19,432 genes tested, 30 had a significant gene-based association at a Bonferroni-corrected P value of 2.5 × 10-6. Of these, 20 were also significantly associated (P < .05/30 = .0016) with disease risk in the replication sample, including 19 that were located in 11 loci not reported to contain allergy risk variants in previous GWASs. Among these were 9 genes with a known function that is directly relevant to allergic disease: FOSL2, VPRBP, IPCEF1, PRR5L, NCF4, APOBR, IL27, ATXN2L, and LAT. For 4 genes (eg, ATXN2L), a genetically determined decrease in gene expression was associated with decreased allergy risk, and therefore drugs that inhibit gene expression or function are predicted to ameliorate disease symptoms. The opposite directional effect was observed for 14 genes, including IL27, a cytokine known to suppress TH2 responses.

CONCLUSION

Using a gene-based approach, we identified 11 risk loci for allergic disease that were not reported in previous GWASs. Functional studies that investigate the contribution of the 19 associated genes to the pathophysiology of allergic disease and assess their therapeutic potential are warranted.

Association of alcohol consumption with allergic disease and asthma: a multi-centre Mendelian randomization analysis

AIMS

To use the rs1229984 variant associated with alcohol consumption as an instrument for alcohol consumption to test the causality of the association of alcohol consumption with hay fever, asthma, allergic sensitization and serum total immunoglobulin (Ig)E.

DESIGN

Observational and Mendelian randomization analyses using genetic variants as unbiased markers of exposure to estimate causal effects, subject to certain assumptions.

SETTING

Europe.

PARTICIPANTS

We included a total of 466 434 people aged 15-82 years from 17 population-based studies conducted from 1997 to 2015.

MEASUREMENTS

The rs1229984 (ADH1B) was genotyped; alcohol consumption, hay fever and asthma were self-reported. Specific and total IgE were measured from serum samples.

FINDINGS

Observational analyses showed that ever-drinking versus non-drinking, but not amount of alcohol intake, was positively associated with hay fever and inversely associated with asthma but not with allergic sensitization or serum total immunoglobulin (Ig)E. However, Mendelian randomization analyses did not suggest that the observational associations are causal. The causal odds ratio (OR) per genetically assessed unit of alcohol/week was an OR = 0.907 [95% confidence interval (CI) = 0.806, 1.019; P = 0.101] for hay fever, an OR = 0.897 (95% CI = 0.790, 1.019; P = 0.095) for asthma, an OR = 0.971 (95% CI =  0.804, 1.174; P = 0.763) for allergic sensitization and a 4.7% change (95% CI = -5.5%, 14.9%; P = 0.366) for total IgE.

CONCLUSIONS

In observational analyses, ever-drinking versus not drinking was positively associated with hay fever and negatively associated with asthma. However, the Mendelian randomization results were not consistent with these associations being causal.

Apparent latent structure within the UK Biobank sample has implications for epidemiological analysis

Large studies use genotype data to discover genetic contributions to complex traits and infer relationships between those traits. Co-incident geographical variation in genotypes and health traits can bias these analyses. Here we show that single genetic variants and genetic scores composed of multiple variants are associated with birth location within UK Biobank and that geographic structure in genotype data cannot be accounted for using routine adjustment for study centre and principal components derived from genotype data. We find that major health outcomes appear geographically structured and that coincident structure in health outcomes and genotype data can yield biased associations. Understanding and accounting for this phenomenon will be important when making inference from genotype data in large studies.

Genome-wide association study of offspring birth weight in 86 577 women identifies five novel loci and highlights maternal genetic effects that are independent of fetal genetics

Genome-wide association studies of birth weight have focused on fetal genetics, whereas relatively little is known about the role of maternal genetic variation. We aimed to identify maternal genetic variants associated with birth weight that could highlight potentially relevant maternal determinants of fetal growth. We meta-analysed data on up to 8.7 million SNPs in up to 86 577 women of European descent from the Early Growth Genetics (EGG) Consortium and the UK Biobank. We used structural equation modelling (SEM) and analyses of mother–child pairs to quantify the separate maternal and fetal genetic effects. Maternal SNPs at 10 loci (MTNR1B, HMGA2, SH2B3, KCNAB1, L3MBTL3, GCK, EBF1, TCF7L2, ACTL9, CYP3A7) were associated with offspring birth weight at P < 5 × 10⁻⁸. In SEM analyses, at least 7 of the 10 associations were consistent with effects of the maternal genotype acting via the intrauterine environment, rather than via effects of shared alleles with the fetus. Variants, or correlated proxies, at many of the loci had been previously associated with adult traits, including fasting glucose (MTNR1B, GCK and TCF7L2) and sex hormone levels (CYP3A7), and one (EBF1) with gestational duration. The identified associations indicate that genetic effects on maternal glucose, cytochrome P450 activity and gestational duration, and potentially on maternal blood pressure and immune function, are relevant for fetal growth. Further characterization of these associations in mechanistic and causal analyses will enhance understanding of the potentially modifiable maternal determinants of fetal growth, with the goal of reducing the morbidity and mortality associated with low and high birth weights.

Evidence of a causal relationship between body mass index and psoriasis: A mendelian randomization study

BACKGROUND

Psoriasis is a common inflammatory skin disease that has been reported to be associated with obesity. We aimed to investigate a possible causal relationship between body mass index (BMI) and psoriasis.

METHODS AND FINDINGS

Following a review of published epidemiological evidence of the association between obesity and psoriasis, mendelian randomization (MR) was used to test for a causal relationship with BMI. We used a genetic instrument comprising 97 single-nucleotide polymorphisms (SNPs) associated with BMI as a proxy for BMI (expected to be much less confounded than measured BMI). One-sample MR was conducted using individual-level data (396,495 individuals) from the UK Biobank and the Nord-Trøndelag Health Study (HUNT), Norway. Two-sample MR was performed with summary-level data (356,926 individuals) from published BMI and psoriasis genome-wide association studies (GWASs). The one-sample and two-sample MR estimates were meta-analysed using a fixed-effect model. To test for a potential reverse causal effect, MR analysis with genetic instruments comprising variants from recent genome-wide analyses for psoriasis were used to test whether genetic risk for this skin disease has a causal effect on BMI. Published observational data showed an association of higher BMI with psoriasis. A mean difference in BMI of 1.26 kg/m2 (95% CI 1.02-1.51) between psoriasis cases and controls was observed in adults, while a 1.55 kg/m2 mean difference (95% CI 1.13-1.98) was observed in children. The observational association was confirmed in UK Biobank and HUNT data sets. Overall, a 1 kg/m2 increase in BMI was associated with 4% higher odds of psoriasis (meta-analysis odds ratio [OR] = 1.04; 95% CI 1.03-1.04; P = 1.73 × 10(-60)). MR analyses provided evidence that higher BMI causally increases the odds of psoriasis (by 9% per 1 unit increase in BMI; OR = 1.09 (1.06-1.12) per 1 kg/m2; P = 4.67 × 10(-9)). In contrast, MR estimates gave little support to a possible causal effect of psoriasis genetic risk on BMI (0.004 kg/m2 change in BMI per doubling odds of psoriasis (-0.003 to 0.011). Limitations of our study include possible misreporting of psoriasis by patients, as well as potential misdiagnosis by clinicians. In addition, there is also limited ethnic variation in the cohorts studied.

CONCLUSIONS

Our study, using genetic variants as instrumental variables for BMI, provides evidence that higher BMI leads to a higher risk of psoriasis. This supports the prioritization of therapies and lifestyle interventions aimed at controlling weight for the prevention or treatment of this common skin disease. Mechanistic studies are required to improve understanding of this relationship.

Using SITAR (SuperImposition by Translation and Rotation) to estimate age at peak height velocity in Avon Longitudinal Study of Parents and Children

Puberty is a time of substantial biological and psychological changes. One of the hallmarks of puberty is a rapid growth spurt, however its timing varies between individuals. The impact of pubertal timing on later health outcomes has been of interest in life course epidemiology, however its measurement can be challenging. Age at peak height velocity (aPHV) offers an objective measure of pubertal timing without having to rely on physical examination or self-report. We describe the derivation of aPHV estimates in Avon Longitudinal Study of Parents and Children (ALSPAC) offspring, using Superimposition by Translation And Rotation (SITAR) mixed effects growth curve analysis. ALSPAC is a rich source of phenotypic and genotypic data and given the importance of pubertal timing for later health outcomes, these data offer an opportunity to explore the determinants and consequences of aPHV.