Exploring polygenic-environment and residual-environment interactions for depressive symptoms within the UK Biobank.

Substantial advances have been made in identifying genetic contributions to depression, but little is known about how the effect of genes can be modulated by the environment, creating a gene-environment interaction. Using multivariate reaction norm models (MRNMs) within the UK Biobank (N = 61294-91644), we investigate whether the polygenic and residual variance components of depressive symptoms are modulated by 17 a priori selected covariate traits-12 environmental variables and 5 biomarkers. MRNMs, a mixed-effects modelling approach, provide unbiased polygenic-covariate interaction estimates for a quantitative trait by controlling for outcome-covariate correlations and residual-covariate interactions. A continuous depressive symptom variable was the outcome in 17 MRNMs-one for each covariate trait. Each MRNM had a fixed-effects model (fixed effects included the covariate trait, demographic variables, and principal components) and a random effects model (where polygenic-covariate and residual-covariate interactions are modelled). Of the 17 selected covariates, 11 significantly modulate deviations in depressive symptoms through the modelled interactions, but no single interaction explains a large proportion of phenotypic variation. Results are dominated by residual-covariate interactions, suggesting that covariate traits (including neuroticism, childhood trauma, and BMI) typically interact with unmodelled variables, rather than a genome-wide polygenic component, to influence depressive symptoms. Only average sleep duration has a polygenic-covariate interaction explaining a demonstrably nonzero proportion of the variability in depressive symptoms. This effect is small, accounting for only 1.22% (95% confidence interval: [0.54, 1.89]) of variation. The presence of an interaction highlights a specific focus for intervention, but the negative results here indicate a limited contribution from polygenic-environment interactions.

Genetic and clinical characteristics of treatment-resistant depression using primary care records in two UK cohorts.

Treatment-resistant depression (TRD) is a major contributor to the disability caused by major depressive disorder (MDD). Primary care electronic health records provide an easily accessible approach to investigate TRD clinical and genetic characteristics. MDD defined from primary care records in UK Biobank (UKB) and EXCEED studies was compared with other measures of depression and tested for association with MDD polygenic risk score (PRS). Using prescribing records, TRD was defined from at least two switches between antidepressant drugs, each prescribed for at least 6 weeks. Clinical-demographic characteristics, SNP-based heritability (h2SNP) and genetic overlap with psychiatric and non-psychiatric traits were compared in TRD and non-TRD MDD cases. In 230,096 and 8926 UKB and EXCEED participants with primary care data, respectively, the prevalence of MDD was 8.7% and 14.2%, of which 13.2% and 13.5% was TRD, respectively. In both cohorts, MDD defined from primary care records was strongly associated with MDD PRS, and in UKB it showed overlap of 71-88% with other MDD definitions. In UKB, TRD vs healthy controls and non-TRD vs healthy controls h2SNP was comparable (0.25 [SE = 0.04] and 0.19 [SE = 0.02], respectively). TRD vs non-TRD was positively associated with the PRS of attention deficit hyperactivity disorder, with lower socio-economic status, obesity, higher neuroticism and other unfavourable clinical characteristics. This study demonstrated that MDD and TRD can be reliably defined using primary care records and provides the first large scale population assessment of the genetic, clinical and demographic characteristics of TRD.

Shared genetic risk between eating disorder- and substance-use-related phenotypes: Evidence from genome-wide association studies.

Eating disorders and substance use disorders frequently co-occur. Twin studies reveal shared genetic variance between liabilities to eating disorders and substance use, with the strongest associations between symptoms of bulimia nervosa and problem alcohol use (genetic correlation [rg ], twin-based = 0.23-0.53). We estimated the genetic correlation between eating disorder and substance use and disorder phenotypes using data from genome-wide association studies (GWAS). Four eating disorder phenotypes (anorexia nervosa [AN], AN with binge eating, AN without binge eating, and a bulimia nervosa factor score), and eight substance-use-related phenotypes (drinks per week, alcohol use disorder [AUD], smoking initiation, current smoking, cigarettes per day, nicotine dependence, cannabis initiation, and cannabis use disorder) from eight studies were included. Significant genetic correlations were adjusted for variants associated with major depressive disorder and schizophrenia. Total study sample sizes per phenotype ranged from ~2400 to ~537 000 individuals. We used linkage disequilibrium score regression to calculate single nucleotide polymorphism-based genetic correlations between eating disorder- and substance-use-related phenotypes. Significant positive genetic associations emerged between AUD and AN (rg = 0.18; false discovery rate q = 0.0006), cannabis initiation and AN (rg = 0.23; q < 0.0001), and cannabis initiation and AN with binge eating (rg = 0.27; q = 0.0016). Conversely, significant negative genetic correlations were observed between three nondiagnostic smoking phenotypes (smoking initiation, current smoking, and cigarettes per day) and AN without binge eating (rgs = -0.19 to -0.23; qs < 0.04). The genetic correlation between AUD and AN was no longer significant after co-varying for major depressive disorder loci. The patterns of association between eating disorder- and substance-use-related phenotypes highlights the potentially complex and substance-specific relationships among these behaviors.

Genetic fine mapping of systemic lupus erythematosus MHC associations in Europeans and African Americans.

Genetic variation within the major histocompatibility complex (MHC) contributes substantial risk for systemic lupus erythematosus, but high gene density, extreme polymorphism and extensive linkage disequilibrium (LD) have made fine mapping challenging. To address the problem, we compared two association techniques in two ancestrally diverse populations, African Americans (AAs) and Europeans (EURs). We observed a greater number of Human Leucocyte Antigen (HLA) alleles in AA consistent with the elevated level of recombination in this population. In EUR we observed 50 different A-C-B-DRB1-DQA-DQB multilocus haplotype sequences per hundred individuals; in the AA sample, these multilocus haplotypes were twice as common compared to Europeans. We also observed a strong narrow class II signal in AA as opposed to the long-range LD observed in EUR that includes class I alleles. We performed a Bayesian model choice of the classical HLA alleles and a frequentist analysis that combined both single nucleotide polymorphisms (SNPs) and classical HLA alleles. Both analyses converged on a similar subset of risk HLA alleles: in EUR HLA- B*08:01 + B*18:01 + (DRB1*15:01 frequentist only) + DQA*01:02 + DQB*02:01 + DRB3*02 and in AA HLA-C*17:01 + B*08:01 + DRB1*15:03 + (DQA*01:02 frequentist only) + DQA*02:01 + DQA*05:01+ DQA*05:05 + DQB*03:19 + DQB*02:02. We observed two additional independent SNP associations in both populations: EUR rs146903072 and rs501480; AA rs389883 and rs114118665. The DR2 serotype was best explained by DRB1*15:03 + DQA*01:02 in AA and by DRB1*15:01 + DQA*01:02 in EUR. The DR3 serotype was best explained by DQA*05:01 in AA and by DQB*02:01 in EUR. Despite some differences in underlying HLA allele risk models in EUR and AA, SNP signals across the extended MHC showed remarkable similarity and significant concordance in direction of effect for risk-associated variants.

Genetic and Inflammatory Biomarkers Classify Small Intestine Inflammation in Asymptomatic First-degree Relatives of Patients With Crohn's Disease.

BACKGROUND & AIMS: Relatives of individuals with Crohn's disease (CD) carry CD-associated genetic variants and are often exposed to environmental factors that increase their risk for this disease. We aimed to estimate the utility of genotype, smoking status, family history, and biomarkers can calculate risk in asymptomatic first-degree relatives of patients with CD. METHODS: We recruited 480 healthy first-degree relatives (full siblings, offspring or parents) of patients with CD through the Guy's and St Thomas' NHS Foundation Trust and from members of Crohn's and Colitis, United Kingdom. DNA samples were genotyped using the Immunochip. We calculated a risk score for 454 participants, based on 72 genetic variants associated with CD, family history, and smoking history. Participants were assigned to highest and lowest risk score quartiles. We assessed pre-symptomatic inflammation by capsule endoscopy and measured 22 markers of inflammation in stool and serum samples (reference standard). Two machine-learning classifiers (elastic net and random forest) were used to assess the ability of the risk factors and biomarkers to identify participants with small intestinal inflammation in the same dataset. RESULTS: The machine-learning classifiers identified participants with pre-symptomatic intestinal inflammation: elastic net (area under the curve, 0.80; 95% CI, 0.62-0.98) and random forest (area under the curve, 0.87; 95% CI, 0.75-1.00). The elastic net method identified 3 variables that can be used to calculate odds for intestinal inflammation: combined family history of CD (odds ratio, 1.31), genetic risk score (odds ratio, 1.14), and fecal calprotectin (odds ratio, 1.04). These same 3 variables were among the 5 factors associated with intestinal inflammation in the random forest model. CONCLUSION: Using machine learning classifiers, we found that genetic variants associated with CD, family history, and fecal calprotectin together identify individuals with pre-symptomatic intestinal inflammation who are therefore at risk for CD. A tool for detecting people at risk for CD before they develop symptoms would help identify the individuals most likely to benefit from early intervention.

Genomics of body fat percentage may contribute to sex bias in anorexia nervosa.

Anorexia nervosa (AN) occurs nine times more often in females than in males. Although environmental factors likely play a role, the reasons for this imbalanced sex ratio remain unresolved. AN displays high genetic correlations with anthropometric and metabolic traits. Given sex differences in body composition, we investigated the possible metabolic underpinnings of female propensity for AN. We conducted sex-specific GWAS in a healthy and medication-free subsample of the UK Biobank (n = 155,961), identifying 77 genome-wide significant loci associated with body fat percentage (BF%) and 174 with fat-free mass (FFM). Partitioned heritability analysis showed an enrichment for central nervous tissue-associated genes for BF%, which was more prominent in females than males. Genetic correlations of BF% and FFM with the largest GWAS of AN by the Psychiatric Genomics Consortium were estimated to explore shared genomics. The genetic correlations of BF%male and BF%female with AN differed significantly from each other (p < .0001, δ = -0.17), suggesting that the female preponderance in AN may, in part, be explained by sex-specific anthropometric and metabolic genetic factors increasing liability to AN.

Genetic variation at 16q24.2 is associated with small vessel stroke.

OBJECTIVE: Genome-wide association studies (GWAS) have been successful at identifying associations with stroke and stroke subtypes, but have not yet identified any associations solely with small vessel stroke (SVS). SVS comprises one quarter of all ischemic stroke and is a major manifestation of cerebral small vessel disease, the primary cause of vascular cognitive impairment. Studies across neurological traits have shown that younger-onset cases have an increased genetic burden. We leveraged this increased genetic burden by performing an age-at-onset informed GWAS meta-analysis, including a large younger-onset SVS population, to identify novel associations with stroke. METHODS: We used a three-stage age-at-onset informed GWAS to identify novel genetic variants associated with stroke. On identifying a novel locus associated with SVS, we assessed its influence on other small vessel disease phenotypes, as well as on messenger RNA (mRNA) expression of nearby genes, and on DNA methylation of nearby CpG sites in whole blood and in the fetal brain. RESULTS: We identified an association with SVS in 4,203 cases and 50,728 controls on chromosome 16q24.2 (odds ratio [OR; 95% confidence interval {CI}] = 1.16 [1.10-1.22]; p = 3.2 × 10-9 ). The lead single-nucleotide polymorphism (rs12445022) was also associated with cerebral white matter hyperintensities (OR [95% CI] = 1.10 [1.05-1.16]; p = 5.3 × 10-5 ; N = 3,670), but not intracerebral hemorrhage (OR [95% CI] = 0.97 [0.84-1.12]; p = 0.71; 1,545 cases, 1,481 controls). rs12445022 is associated with mRNA expression of ZCCHC14 in arterial tissues (p = 9.4 × 10-7 ) and DNA methylation at probe cg16596957 in whole blood (p = 5.3 × 10-6 ). INTERPRETATION: 16q24.2 is associated with SVS. Associations of the locus with expression of ZCCHC14 and DNA methylation suggest the locus acts through changes to regulatory elements. Ann Neurol 2017;81:383-394.

Trajectories of Callous-Unemotional Traits in Childhood Predict Different Forms of Peer Victimization in Adolescence.

Callous-unemotional (CU) traits (e.g., lack of empathy and guilt) differentiate a group of children at particularly high risk for engaging in aggressive behavior, notably bullying. However, little is known about whether youths with CU traits are at risk for being victimized by their peers. We examined the associations between trajectories of CU traits in childhood (between 7 and 12 years old) and peer victimization in adolescence (14 years old). The participants were drawn from the Twins Early Development Study, a longitudinal population-based study of twins born in England and in Wales. The trajectories of CU traits (i.e., stable high, increasing, decreasing and stable low) were identified through general growth mixture modeling. Four forms of peer victimization were considered: physical victimization, verbal victimization, social manipulation, and attacks on property. We found that youths with stable high levels, increasing levels, and decreasing levels of CU traits in childhood had higher levels of physical victimization in adolescence, not explained by other predictors at age 7 (e.g., conduct problems). Youths with increasing levels of CU traits, compared with the ones with stable low levels, also had higher levels of verbal victimization, social manipulation, and attacks on property. Our findings highlight the importance of distinct trajectories of CU traits in accounting for the experience of different forms of peer victimization. Youths with CU traits may benefit from bullying prevention programs, as they are likely to be the targets of peer victimization.

Genetic Factors Influencing Coagulation Factor XIII B-Subunit Contribute to Risk of Ischemic Stroke.

BACKGROUND AND PURPOSE: Abnormal coagulation has been implicated in the pathogenesis of ischemic stroke, but how this association is mediated and whether it differs between ischemic stroke subtypes is unknown. We determined the shared genetic risk between 14 coagulation factors and ischemic stroke and its subtypes. METHODS: Using genome-wide association study results for 14 coagulation factors from the population-based TwinsUK sample (N≈2000 for each factor), meta-analysis results from the METASTROKE consortium ischemic stroke genome-wide association study (12 389 cases, 62 004 controls), and genotype data for 9520 individuals from the WTCCC2 ischemic stroke study (3548 cases, 5972 controls-the largest METASTROKE subsample), we explored shared genetic risk for coagulation and stroke. We performed three analyses: (1) a test for excess concordance (or discordance) in single nucleotide polymorphism effect direction across coagulation and stroke, (2) an estimation of the joint effect of multiple coagulation-associated single nucleotide polymorphisms in stroke, and (3) an evaluation of common genetic risk between coagulation and stroke. RESULTS: One coagulation factor, factor XIII subunit B (FXIIIB), showed consistent effects in the concordance analysis, the estimation of polygenic risk, and the validation with genotype data, with associations specific to the cardioembolic stroke subtype. Effect directions for FXIIIB-associated single nucleotide polymorphisms were significantly discordant with cardioembolic disease (smallest P=5.7×10(-04)); the joint effect of FXIIIB-associated single nucleotide polymorphisms was significantly predictive of ischemic stroke (smallest P=1.8×10(-04)) and the cardioembolic subtype (smallest P=1.7×10(-04)). We found substantial negative genetic covariation between FXIIIB and ischemic stroke (rG=-0.71, P=0.01) and the cardioembolic subtype (rG=-0.80, P=0.03). CONCLUSIONS: Genetic markers associated with low FXIIIB levels increase risk of ischemic stroke cardioembolic subtype.

Why do spatial abilities predict mathematical performance?

Spatial ability predicts performance in mathematics and eventual expertise in science, technology and engineering. Spatial skills have also been shown to rely on neuronal networks partially shared with mathematics. Understanding the nature of this association can inform educational practices and intervention for mathematical underperformance. Using data on two aspects of spatial ability and three domains of mathematical ability from 4174 pairs of 12-year-old twins, we examined the relative genetic and environmental contributions to variation in spatial ability and to its relationship with different aspects of mathematics. Environmental effects explained most of the variation in spatial ability (~70%) and in mathematical ability (~60%) at this age, and the effects were the same for boys and girls. Genetic factors explained about 60% of the observed relationship between spatial ability and mathematics, with a substantial portion of the relationship explained by common environmental influences (26% and 14% by shared and non-shared environments respectively). These findings call for further research aimed at identifying specific environmental mediators of the spatial-mathematics relationship.

Understanding the science-learning environment: A genetically sensitive approach.

Previous studies have shown that environmental influences on school science performance increase in importance from primary to secondary school. Here we assess for the first time the relationship between the science-learning environment and science performance using a genetically sensitive approach to investigate the aetiology of this link. 3000 pairs of 14-year-old twins from the UK Twins Early Development Study reported on their experiences of the science-learning environment and were assessed for their performance in science using a web-based test of scientific enquiry. Multivariate twin analyses were used to investigate the genetic and environmental links between environment and outcome. The most surprising result was that the science-learning environment was almost as heritable (43%) as performance on the science test (50%), and showed negligible shared environmental influence (3%). Genetic links explained most (56%) of the association between learning environment and science outcome, indicating gene-environment correlation.

Chaotic homes and children's disruptive behavior: a longitudinal cross-lagged twin study.

Chaotic home lives are correlated with behavior problems in children. In the study reported here, we tested whether there was a cross-lagged relation between children's experience of chaos and their disruptive behaviors (conduct problems and hyperactivity-inattention). Using genetically informative models, we then tested for the first time whether the influence of household chaos on disruptive behavior was environmentally mediated and whether genetic influences on children's disruptive behaviors accounted for the heritability of household chaos. We measured children's perceptions of household chaos and parents' ratings of children's disruptive behavior at ages 9 and 12 in a sample of 6,286 twin pairs from the Twins Early Development Study (TEDS). There was a phenotypic cross-lagged relation between children's experiences of household chaos and their disruptive behavior. In genetically informative models, we found that the effect of household chaos on subsequent disruptive behavior was environmentally mediated. However, genetic influences on disruptive behavior did not explain why household chaos was heritable.

Chaotic homes and school achievement: a twin study.

BACKGROUND: Chaotic homes predict poor school performance. Given that it is known that genes affect both children's experience of household chaos and their school achievement, to what extent is the relationship between high levels of noise and environmental confusion in the home, and children's school performance, mediated by heritable child effects? This is the first study to explore the genetic and environmental pathways between household chaos and academic performance. METHOD: Children's perceptions of family chaos at ages 9 and 12 and their school performance at age 12 were assessed in more than 2,300 twin pairs. The use of child-specific measures in a multivariate genetic analysis made it possible to investigate the genetic and environmental origins of the covariation between children's experience of chaos in the home and their school achievement. RESULTS: Children's experience of family chaos and their school achievement were significantly correlated in the expected negative direction (r = -.26). As expected, shared environmental factors explained a large proportion (63%) of the association. However, genetic factors accounted for a significant proportion (37%) of the association between children's experience of household chaos and their school performance. CONCLUSIONS: The association between chaotic homes and poor performance in school, previously assumed to be entirely environmental in origin, is in fact partly genetic. How children's home environment affects their academic achievement is not simply in the direction environment → child → outcome. Instead, genetic factors that influence children's experience of the disordered home environment also affect how well they do at school. The relationship between the child, their environment and their performance at school is complex: both genetic and environmental factors play a role.

The genetic case for cardiorespiratory fitness as a clinical vital sign and the routine prescription of physical activity in healthcare.

BACKGROUND: Cardiorespiratory fitness (CRF) and physical activity (PA) are well-established predictors of morbidity and all-cause mortality. However, CRF is not routinely measured and PA not routinely prescribed as part of standard healthcare. The American Heart Association (AHA) recently presented a scientific case for the inclusion of CRF as a clinical vital sign based on epidemiological and clinical observation. Here, we leverage genetic data in the UK Biobank (UKB) to strengthen the case for CRF as a vital sign and make a case for the prescription of PA. METHODS: We derived two CRF measures from the heart rate data collected during a submaximal cycle ramp test: CRF-vo2max, an estimate of the participants' maximum volume of oxygen uptake, per kilogram of body weight, per minute; and CRF-slope, an estimate of the rate of increase of heart rate during exercise. Average PA over a 7-day period was derived from a wrist-worn activity tracker. After quality control, 70,783 participants had data on the two derived CRF measures, and 89,683 had PA data. We performed genome-wide association study (GWAS) analyses by sex, and post-GWAS techniques to understand genetic architecture of the traits and prioritise functional genes for follow-up. RESULTS: We found strong evidence that genetic variants associated with CRF and PA influenced genetic expression in a relatively small set of genes in the heart, artery, lung, skeletal muscle and adipose tissue. These functionally relevant genes were enriched among genes known to be associated with coronary artery disease (CAD), type 2 diabetes (T2D) and Alzheimer's disease (three of the top 10 causes of death in high-income countries) as well as Parkinson's disease, pulmonary fibrosis, and blood pressure, heart rate, and respiratory phenotypes. Genetic variation associated with lower CRF and PA was also correlated with several disease risk factors (including greater body mass index, body fat and multiple obesity phenotypes); a typical T2D profile (including higher insulin resistance, higher fasting glucose, impaired beta-cell function, hyperglycaemia, hypertriglyceridemia); increased risk for CAD and T2D; and a shorter lifespan. CONCLUSIONS: Genetics supports three decades of evidence for the inclusion of CRF as a clinical vital sign. Given the genetic, clinical and epidemiological evidence linking CRF and PA to increased morbidity and mortality, regular measurement of CRF as a marker of health and routine prescription of PA could be a prudent strategy to support public health.

Multiple measures of depression to enhance validity of major depressive disorder in the UK Biobank.

BACKGROUND: The UK Biobank contains data with varying degrees of reliability and completeness for assessing depression. A third of participants completed a Mental Health Questionnaire (MHQ) containing the gold-standard Composite International Diagnostic Interview (CIDI) criteria for assessing mental health disorders. AIMS: To investigate whether multiple observations of depression from sources other than the MHQ can enhance the validity of major depressive disorder (MDD). METHOD: In participants who did not complete the MHQ, we calculated the number of other depression measures endorsed, for example from hospital episode statistics and interview data. We compared cases defined this way with CIDI-defined cases for several estimates: the variance explained by polygenic risk scores (PRS), area under the curve attributable to PRS, single nucleotide polymorphisms (SNPs)-based heritability and genetic correlations with summary statistics from the Psychiatric Genomics Consortium MDD genome-wide association study. RESULTS: The strength of the genetic contribution increased with the number of measures endorsed. For example, SNP-based heritability increased from 7% in participants who endorsed only one measure of depression, to 21% in those who endorsed four or five measures of depression. The strength of the genetic contribution to cases defined by at least two measures approximated that for CIDI-defined cases. Most genetic correlations between UK Biobank and the Psychiatric Genomics Consortium MDD study exceeded 0.7, but there was variability between pairwise comparisons. CONCLUSIONS: Multiple measures of depression can serve as a reliable approximation for case status where the CIDI measure is not available, indicating sample size can be optimised using the entire suite of UK Biobank data.

In search of genes associated with risk for psychopathic tendencies in children: a two-stage genome-wide association study of pooled DNA.

BACKGROUND: Quantitative genetic data from our group indicates that antisocial behaviour (AB) is strongly heritable when coupled with psychopathic, callous-unemotional (CU) personality traits. We have also demonstrated that the genetic influences for AB and CU overlap considerably. We conducted a genome-wide association scan that capitalises on these findings in an attempt to identify quantitative trait loci (QTLs) that may increase risk for psychopathic tendencies (AB+/CU+). METHODS: Teacher ratings at age 7 were used to screen 8374 twins with available DNA samples for individuals that were high vs. low on both AB and CU. In Stage 1, we screened for allele frequency differences in 642,432 autosomal single-nucleotide polymorphisms (SNPs) using the Affymetrix 6.0 GeneChip with pooled DNA for high-scoring (AB+/CU+) versus low-scoring children (N = approximately 300/group). In Stage 2, we tested the 3000 most strongly associated SNPs from Stage 1 for association in the same direction in a second sample of high- versus low-scoring children from the same twin study (18% co-twins). RESULTS: Using allele frequencies estimated from pooled DNA, we found suggestive evidence for enrichment of association in the second stage of our two-stage genome-wide association design and focus on reporting the 30 top-ranking SNPs nominally associated with psychopathic tendencies. These SNPs include neurodevelopmental genes such as ROBO2. CONCLUSIONS: Although none of the SNPs reached genome-wide statistical significance we have generated a list of SNPs that are potentially associated with psychopathic tendencies, which we believe warrant verification and replication in large independent and clinical samples.

The nature (and nurture) of children's perceptions of family chaos.

Chaos in the home is a key environment in cognitive and behavioral development. However, we show that children's experience of home chaos is partly genetically mediated. We assessed children's perceptions of household chaos at ages 9 and 12 in 2337 pairs of twins. Using child-specific reports allowed us to use structural equation modeling to explore the genetic and environmental etiology of children's perceptions of chaos. We found that these perceptions are significantly heritable (22%), with the remainder explained by environmental influences. Finding that genes influence children's experience of chaotic environments has far-reaching implications for how we conceptualize the family home and its impact on cognitive and behavioral development.

Genome-wide association study of MRI markers of cerebral small vessel disease in 42,310 participants.

Cerebral small vessel disease is a major cause of stroke and dementia, but its genetic basis is incompletely understood. We perform a genetic study of three MRI markers of the disease in UK Biobank imaging data and other sources: white matter hyperintensities (N = 42,310), fractional anisotropy (N = 17,663) and mean diffusivity (N = 17,467). Our aim is to better understand the disease pathophysiology. Across the three traits, we identify 31 loci, of which 21 were previously unreported. We perform a transcriptome-wide association study to identify associations with gene expression in relevant tissues, identifying 66 associated genes across the three traits. This genetic study provides insights into the understanding of the biological mechanisms underlying small vessel disease.

Classical Human Leukocyte Antigen Alleles and C4 Haplotypes Are Not Significantly Associated With Depression.

BACKGROUND: The prevalence of depression is higher in individuals with autoimmune diseases, but the mechanisms underlying the observed comorbidities are unknown. Shared genetic etiology is a plausible explanation for the overlap, and in this study we tested whether genetic variation in the major histocompatibility complex (MHC), which is associated with risk for autoimmune diseases, is also associated with risk for depression. METHODS: We fine-mapped the classical MHC (chr6: 29.6-33.1 Mb), imputing 216 human leukocyte antigen (HLA) alleles and 4 complement component 4 (C4) haplotypes in studies from the Psychiatric Genomics Consortium Major Depressive Disorder Working Group and the UK Biobank. The total sample size was 45,149 depression cases and 86,698 controls. We tested for association between depression status and imputed MHC variants, applying both a region-wide significance threshold (3.9 × 10-6) and a candidate threshold (1.6 × 10-4). RESULTS: No HLA alleles or C4 haplotypes were associated with depression at the region-wide threshold. HLA-B*08:01 was associated with modest protection for depression at the candidate threshold for testing in HLA genes in the meta-analysis (odds ratio = 0.98, 95% confidence interval = 0.97-0.99). CONCLUSIONS: We found no evidence that an increased risk for depression was conferred by HLA alleles, which play a major role in the genetic susceptibility to autoimmune diseases, or C4 haplotypes, which are strongly associated with schizophrenia. These results suggest that any HLA or C4 variants associated with depression either are rare or have very modest effect sizes.

ukbtools: An R package to manage and query UK Biobank data.

INTRODUCTION: The UK Biobank (UKB) is a resource that includes detailed health-related data on about 500,000 individuals and is available to the research community. However, several obstacles limit immediate analysis of the data: data files vary in format, may be very large, and have numerical codes for column names. RESULTS: ukbtools removes all the upfront data wrangling required to get a single dataset for statistical analysis. All associated data files are merged into a single dataset with descriptive column names. The package also provides tools to assist in quality control by exploring the primary demographics of subsets of participants; query of disease diagnoses for one or more individuals, and estimating disease frequency relative to a reference variable; and to retrieve genetic metadata. CONCLUSION: Having a dataset with meaningful variable names, a set of UKB-specific exploratory data analysis tools, disease query functions, and a set of helper functions to explore and write genetic metadata to file, will rapidly enable UKB users to undertake their research.