Polygenic risk scores associate with blood pressure traits across the lifespan

AIMS

Hypertension is a major modifiable cause of morbidity and mortality that affects over 1 billion people worldwide. Blood pressure (BP) traits have a strong genetic component that can be quantified with polygenic risk scores (PRSs). To date, the performance of BP PRSs has mainly been assessed in adults, and less is known about polygenic hypertension risk in childhood.

METHODS AND RESULTS

Multiple PRSs for systolic BP (SBP), diastolic BP (DBP), and pulse pressure were developed using either genome-wide significant weights, pruning and thresholding, or Bayesian regression. Among 87 total PRSs, the top performer for each trait was applied in independent cohorts of children and adult to assess genotype-phenotype associations and disease risk across the lifespan. Differences between those with low (1st decile), average (2nd-9th decile), and high (10th decile) PRS emerge in the first years of life and are maintained throughout adulthood. These diverging BP trajectories also seem to affect cardiovascular and renal disease risk, with increased risk observed among those in the top decile and reduced risk among those in the bottom decile of the polygenic risk distribution compared with the rest of the population.

CONCLUSION

Genetic risk factors are associated with BP traits across the lifespan, beginning in the first years of life. Given the importance of exposure time in disease pathogenesis and the early rise in BP levels among those genetically susceptible, PRSs may help identify high-risk individuals prior to hypertension onset, facilitate primordial prevention, and reduce the burden of this public health challenge.

A genome-wide association study provides insights into the genetic etiology of 57 essential and non-essential trace elements in humans

Trace elements are important for human health but may exert toxic or adverse effects. Mechanisms of uptake, distribution, metabolism, and excretion are partly under genetic control but have not yet been extensively mapped. Here we report a comprehensive multi-element genome-wide association study of 57 essential and non-essential trace elements. We perform genome-wide association meta-analyses of 14 trace elements in up to 6564 Scandinavian whole blood samples, and genome-wide association studies of 43 trace elements in up to 2819 samples measured only in the Trøndelag Health Study (HUNT). We identify 11 novel genetic loci associated with blood concentrations of arsenic, cadmium, manganese, selenium, and zinc in genome-wide association meta-analyses. In HUNT, several genome-wide significant loci are also indicated for other trace elements. Using two-sample Mendelian randomization, we find several indications of weak to moderate effects on health outcomes, the most precise being a weak harmful effect of increased zinc on prostate cancer. However, independent validation is needed. Our current understanding of trace element-associated genetic variants may help establish consequences of trace elements on human health.

X-chromosome and kidney function: evidence from a multi-trait genetic analysis of 908,697 individuals reveals sex-specific and sex-differential findings in genes regulated by androgen response elements

X-chromosomal genetic variants are understudied but can yield valuable insights into sexually dimorphic human traits and diseases. We performed a sex-stratified cross-ancestry X-chromosome-wide association meta-analysis of seven kidney-related traits (n = 908,697), identifying 23 loci genome-wide significantly associated with two of the traits: 7 for uric acid and 16 for estimated glomerular filtration rate (eGFR), including four novel eGFR loci containing the functionally plausible prioritized genes ACSL4, CLDN2, TSPAN6 and the female-specific DRP2. Further, we identified five novel sex-interactions, comprising male-specific effects at FAM9B and AR/EDA2R, and three sex-differential findings with larger genetic effect sizes in males at DCAF12L1 and MST4 and larger effect sizes in females at HPRT1. All prioritized genes in loci showing significant sex-interactions were located next to androgen response elements (ARE). Five ARE genes showed sex-differential expressions. This study contributes new insights into sex-dimorphisms of kidney traits along with new prioritized gene targets for further molecular research.

GWAS meta-analysis of psoriasis identifies new susceptibility alleles impacting disease mechanisms and therapeutic targets

Psoriasis is a common, debilitating immune-mediated skin disease. Genetic studies have identified biological mechanisms of psoriasis risk, including those targeted by effective therapies. However, the genetic liability to psoriasis is not fully explained by variation at robustly identified risk loci. To move towards a saturation map of psoriasis susceptibility we meta-analysed 18 GWAS comprising 36,466 cases and 458,078 controls and identified 109 distinct psoriasis susceptibility loci, including 45 that have not been previously reported. These include susceptibility variants at loci in which the therapeutic targets IL17RA and AHR are encoded, and deleterious coding variants supporting potential new drug targets (including in STAP2, CPVL and POU2F3). We conducted a transcriptome-wide association study to identify regulatory effects of psoriasis susceptibility variants and cross-referenced these against single cell expression profiles in psoriasis-affected skin, highlighting roles for the transcriptional regulation of haematopoietic cell development and epigenetic modulation of interferon signalling in psoriasis pathobiology.

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.

Cluster Headache Genomewide Association Study and Meta-Analysis Identifies Eight Loci and Implicates Smoking as Causal Risk Factor

OBJECTIVE

The objective of this study was to aggregate data for the first genomewide association study meta-analysis of cluster headache, to identify genetic risk variants, and gain biological insights.

METHODS

A total of 4,777 cases (3,348 men and 1,429 women) with clinically diagnosed cluster headache were recruited from 10 European and 1 East Asian cohorts. We first performed an inverse-variance genomewide association meta-analysis of 4,043 cases and 21,729 controls of European ancestry. In a secondary trans-ancestry meta-analysis, we included 734 cases and 9,846 controls of East Asian ancestry. Candidate causal genes were prioritized by 5 complementary methods: expression quantitative trait loci, transcriptome-wide association, fine-mapping of causal gene sets, genetically driven DNA methylation, and effects on protein structure. Gene set and tissue enrichment analyses, genetic correlation, genetic risk score analysis, and Mendelian randomization were part of the downstream analyses.

RESULTS

The estimated single nucleotide polymorphism (SNP)-based heritability of cluster headache was 14.5%. We identified 9 independent signals in 7 genomewide significant loci in the primary meta-analysis, and one additional locus in the trans-ethnic meta-analysis. Five of the loci were previously known. The 20 genes prioritized as potentially causal for cluster headache showed enrichment to artery and brain tissue. Cluster headache was genetically correlated with cigarette smoking, risk-taking behavior, attention deficit hyperactivity disorder (ADHD), depression, and musculoskeletal pain. Mendelian randomization analysis indicated a causal effect of cigarette smoking intensity on cluster headache. Three of the identified loci were shared with migraine.

INTERPRETATION

This first genomewide association study meta-analysis gives clues to the biological basis of cluster headache and indicates that smoking is a causal risk factor. ANN NEUROL 2023;94:713-726.

Body Composition, Cardiometabolic Risk Factors and Comorbidities in Psoriasis and the Effect of HLA-C*06:02 Status: The HUNT Study, Norway

Abstract is missing (Short communication)

Genetic loci and prioritization of genes for kidney function decline derived from a meta-analysis of 62 longitudinal genome-wide association studies

Estimated glomerular filtration rate (eGFR) reflects kidney function. Progressive eGFR-decline can lead to kidney failure, necessitating dialysis or transplantation. Hundreds of loci from genome-wide association studies (GWAS) for eGFR help explain population cross section variability. Since the contribution of these or other loci to eGFR-decline remains largely unknown, we derived GWAS for annual eGFR-decline and meta-analyzed 62 longitudinal studies with eGFR assessed twice over time in all 343,339 individuals and in high-risk groups. We also explored different covariate adjustment. Twelve genome-wide significant independent variants for eGFR-decline unadjusted or adjusted for eGFR-baseline (11 novel, one known for this phenotype), including nine variants robustly associated across models were identified. All loci for eGFR-decline were known for cross-sectional eGFR and thus distinguished a subgroup of eGFR loci. Seven of the nine variants showed variant-by-age interaction on eGFR cross section (further about 350,000 individuals), which linked genetic associations for eGFR-decline with age-dependency of genetic cross-section associations. Clinically important were two to four-fold greater genetic effects on eGFR-decline in high-risk subgroups. Five variants associated also with chronic kidney disease progression mapped to genes with functional in-silico evidence (UMOD, SPATA7, GALNTL5, TPPP). An unfavorable versus favorable nine-variant genetic profile showed increased risk odds ratios of 1.35 for kidney failure (95% confidence intervals 1.03-1.77) and 1.27 for acute kidney injury (95% confidence intervals 1.08-1.50) in over 2000 cases each, with matched controls). Thus, we provide a large data resource, genetic loci, and prioritized genes for kidney function decline, which help inform drug development pipelines revealing important insights into the age-dependency of kidney function genetics.

Differential and shared genetic effects on kidney function between diabetic and non-diabetic individuals

Reduced glomerular filtration rate (GFR) can progress to kidney failure. Risk factors include genetics and diabetes mellitus (DM), but little is known about their interaction. We conducted genome-wide association meta-analyses for estimated GFR based on serum creatinine (eGFR), separately for individuals with or without DM (nDM = 178,691, nnoDM = 1,296,113). Our genome-wide searches identified (i) seven eGFR loci with significant DM/noDM-difference, (ii) four additional novel loci with suggestive difference and (iii) 28 further novel loci (including CUBN) by allowing for potential difference. GWAS on eGFR among DM individuals identified 2 known and 27 potentially responsible loci for diabetic kidney disease. Gene prioritization highlighted 18 genes that may inform reno-protective drug development. We highlight the existence of DM-only and noDM-only effects, which can inform about the target group, if respective genes are advanced as drug targets. Largely shared effects suggest that most drug interventions to alter eGFR should be effective in DM and noDM.

Risk of lower respiratory tract infections: a genome-wide association study with Mendelian randomization analysis in three independent European populations

OBJECTIVE

Lower respiratory tract infections (LRTIs) are a leading cause of morbidity and mortality worldwide. Few studies have previously investigated genetic susceptibility and potential risk factors for LRTI.

METHODS

We used data from the UK Biobank, Trøndelag Health Study (HUNT), and FinnGen to conduct a genome-wide association study (GWAS). Cases were subjects hospitalized with LRTI, and controls were subjects with no such hospitalization. We conducted stratification and interaction analyses to evaluate whether the genetic effect of LRTI differed by sex or smoking. Mendelian randomization (MR) analyses were conducted to identify the unconfounded relationship between cardiometabolic risk factors and LRTI.

RESULTS

A total of 25 320 cases and 575 294 controls were included. The 15q25.1 locus reached genome-wide significance in the meta-analysis (rs10519203: OR 0.94, p 3.87e-11). The protective effect of effect allele of rs10519203 was present among smokers (OR 0.90, 95%CI 0.87-0.92, p 1.38e-15) but not among never-smokers (OR 1.01, 95%CI 0.97-1.06, p 5.20e-01). In MR analyses, we found that increasing body mass index (OR 1.31, 95%CI 1.24-1.40, p 3.78e-18), lifetime smoking (OR 2.83, 95%CI 2.34-3.42, p 6.56e-27), and systolic blood pressure robustly increased the risk of LRTIs (OR 1.11, 95%CI 1.02-1.22, p 1.48e-02).

CONCLUSION

A region in 15q25.1 was strongly associated with LRTI susceptibility. Reduction in the prevalence of smoking, overweight, obesity, and hypertension may reduce the disease burden of LRTIs.

Genome-wide risk prediction of common diseases across ancestries in one million people

Polygenic risk scores (PRS) measure genetic disease susceptibility by combining risk effects across the genome. For coronary artery disease (CAD), type 2 diabetes (T2D), and breast and prostate cancer, we performed cross-ancestry evaluation of genome-wide PRSs in six biobanks in Europe, the United States, and Asia. We studied transferability of these highly polygenic, genome-wide PRSs across global ancestries, within European populations with different health-care systems, and local population substructures in a population isolate. All four PRSs had similar accuracy across European and Asian populations, with poorer transferability in the smaller group of individuals of African ancestry. The PRSs had highly similar effect sizes in different populations of European ancestry, and in early- and late-settlement regions with different recent population bottlenecks in Finland. Comparing genome-wide PRSs to PRSs containing a smaller number of variants, the highly polygenic, genome-wide PRSs generally displayed higher effect sizes and better transferability across global ancestries. Our findings indicate that in the populations investigated, the current genome-wide polygenic scores for common diseases have potential for clinical utility within different health-care settings for individuals of European ancestry, but that the utility in individuals of African ancestry is currently much lower.

Genome-wide analysis of 102,084 migraine cases identifies 123 risk loci and subtype-specific risk alleles

Migraine affects over a billion individuals worldwide but its genetic underpinning remains largely unknown. Here, we performed a genome-wide association study of 102,084 migraine cases and 771,257 controls and identified 123 loci, of which 86 are previously unknown. These loci provide an opportunity to evaluate shared and distinct genetic components in the two main migraine subtypes: migraine with aura and migraine without aura. Stratification of the risk loci using 29,679 cases with subtype information indicated three risk variants that seem specific for migraine with aura (in HMOX2, CACNA1A and MPPED2), two that seem specific for migraine without aura (near SPINK2 and near FECH) and nine that increase susceptibility for migraine regardless of subtype. The new risk loci include genes encoding recent migraine-specific drug targets, namely calcitonin gene-related peptide (CALCA/CALCB) and serotonin 1F receptor (HTR1F). Overall, genomic annotations among migraine-associated variants were enriched in both vascular and central nervous system tissue/cell types, supporting unequivocally that neurovascular mechanisms underlie migraine pathophysiology.

Trans-ethnic Mendelian-randomization study reveals causal relationships between cardiometabolic factors and chronic kidney disease

BACKGROUND

This study was to systematically test whether previously reported risk factors for chronic kidney disease (CKD) are causally related to CKD in European and East Asian ancestries using Mendelian randomization.

METHODS

A total of 45 risk factors with genetic data in European ancestry and 17 risk factors in East Asian participants were identified as exposures from PubMed. We defined the CKD by clinical diagnosis or by estimated glomerular filtration rate of <60 ml/min/1.73 m2. Ultimately, 51 672 CKD cases and 958 102 controls of European ancestry from CKDGen, UK Biobank and HUNT, and 13 093 CKD cases and 238 118 controls of East Asian ancestry from Biobank Japan, China Kadoorie Biobank and Japan-Kidney-Biobank/ToMMo were included.

RESULTS

Eight risk factors showed reliable evidence of causal effects on CKD in Europeans, including genetically predicted body mass index (BMI), hypertension, systolic blood pressure, high-density lipoprotein cholesterol, apolipoprotein A-I, lipoprotein(a), type 2 diabetes (T2D) and nephrolithiasis. In East Asians, BMI, T2D and nephrolithiasis showed evidence of causality on CKD. In two independent replication analyses, we observed that increased hypertension risk showed reliable evidence of a causal effect on increasing CKD risk in Europeans but in contrast showed a null effect in East Asians. Although liability to T2D showed consistent effects on CKD, the effects of glycaemic phenotypes on CKD were weak. Non-linear Mendelian randomization indicated a threshold relationship between genetically predicted BMI and CKD, with increased risk at BMI of >25 kg/m2.

CONCLUSIONS

Eight cardiometabolic risk factors showed causal effects on CKD in Europeans and three of them showed causality in East Asians, providing insights into the design of future interventions to reduce the burden of CKD.

Genome-Wide Association Study of 2,093 Cases With Idiopathic Polyneuropathy and 445,256 Controls Identifies First Susceptibility Loci

Background: About one third of patients with chronic polyneuropathy have no obvious underlying etiology and are classified as having idiopathic polyneuropathy. The lack of knowledge about pathomechanisms and predisposing factors limits the development of effective prevention and treatment for these patients. We report the first genome-wide association study (GWAS) of idiopathic polyneuropathy.

Methods: Cases with idiopathic polyneuropathy and healthy controls were identified by linkage to hospital records. We performed genome-wide association studies using genetic data from two large population-based health studies, the Trøndelag Health Study (HUNT, 1,147 cases and 62,204 controls) and UK Biobank (UKB, 946 cases and 383,052 controls). In a two-stage analysis design, we first treated HUNT as a discovery cohort and UK Biobank as a replication cohort. Secondly, we combined the two studies in a meta-analysis. Downstream analyses included genetic correlation to other traits and diseases. We specifically examined previously reported risk loci, and genes known to cause hereditary polyneuropathy.

Results: No replicable risk loci were identified in the discovery-replication stage, in line with the limited predicted power of this approach. When combined in a meta-analysis, two independent loci reached statistical significance (rs7294354 in B4GALNT3, P-value 4.51 × 10⁻⁸) and (rs147738081 near NR5A2, P-value 4.75 × 10⁻⁸). Idiopathic polyneuropathy genetically correlated with several anthropometric measures, most pronounced for height, and with several pain-related traits.

Conclusions: In this first GWAS of idiopathic polyneuropathy we identify two risk-loci that indicate possible pathogenetic mechanisms. Future collaborative efforts are needed to replicate and expand on these findings.

Genome-wide association study of cardiac troponin I in the general population

Circulating cardiac troponin proteins are associated with structural heart disease and predict incident cardiovascular disease in the general population. However, the genetic contribution to cardiac troponin I (cTnI) concentrations and its causal effect on cardiovascular phenotypes are unclear. We combine data from two large population-based studies, the Trøndelag Health Study and the Generation Scotland Scottish Family Health Study, and perform a genome-wide association study of high-sensitivity cTnI concentrations with 48 115 individuals. We further use two-sample Mendelian randomization to investigate the causal effects of circulating cTnI on acute myocardial infarction (AMI) and heart failure (HF). We identified 12 genetic loci (8 novel) associated with cTnI concentrations. Associated protein-altering variants highlighted putative functional genes: CAND2, HABP2, ANO5, APOH, FHOD3, TNFAIP2, KLKB1 and LMAN1. Phenome-wide association tests in 1688 phecodes and 83 continuous traits in UK Biobank showed associations between a genetic risk score for cTnI and cardiac arrhythmias, metabolic and anthropometric measures. Using two-sample Mendelian randomization, we confirmed the non-causal role of cTnI in AMI (5948 cases, 355 246 controls). We found indications for a causal role of cTnI in HF (47 309 cases and 930 014 controls), but this was not supported by secondary analyses using left ventricular mass as outcome (18 257 individuals). Our findings clarify the biology underlying the heritable contribution to circulating cTnI and support cTnI as a non-causal biomarker for AMI in the general population. Using genetically informed methods for causal inference helps inform the role and value of measuring cTnI in the general population.

Deciphering osteoarthritis genetics across 826,690 individuals from 9 populations

Osteoarthritis affects over 300 million people worldwide. Here, we conduct a genome-wide association study meta-analysis across 826,690 individuals (177,517 with osteoarthritis) and identify 100 independently associated risk variants across 11 osteoarthritis phenotypes, 52 of which have not been associated with the disease before. We report thumb and spine osteoarthritis risk variants and identify differences in genetic effects between weight-bearing and non-weight-bearing joints. We identify sex-specific and early age-at-onset osteoarthritis risk loci. We integrate functional genomics data from primary patient tissues (including articular cartilage, subchondral bone, and osteophytic cartilage) and identify high-confidence effector genes. We provide evidence for genetic correlation with phenotypes related to pain, the main disease symptom, and identify likely causal genes linked to neuronal processes. Our results provide insights into key molecular players in disease processes and highlight attractive drug targets to accelerate translation.

A genome-wide association study with 1,126,563 individuals identifies new risk loci for Alzheimer's disease

Late-onset Alzheimer's disease is a prevalent age-related polygenic disease that accounts for 50-70% of dementia cases. Currently, only a fraction of the genetic variants underlying Alzheimer's disease have been identified. Here we show that increased sample sizes allowed identification of seven previously unidentified genetic loci contributing to Alzheimer's disease. This study highlights microglia, immune cells and protein catabolism as relevant to late-onset Alzheimer's disease, while identifying and prioritizing previously unidentified genes of potential interest. We anticipate that these results can be included in larger meta-analyses of Alzheimer's disease to identify further genetic variants that contribute to Alzheimer's pathology.

Discovery and prioritization of variants and genes for kidney function in >1.2 million individuals

Genes underneath signals from genome-wide association studies (GWAS) for kidney function are promising targets for functional studies, but prioritizing variants and genes is challenging. By GWAS meta-analysis for creatinine-based estimated glomerular filtration rate (eGFR) from the Chronic Kidney Disease Genetics Consortium and UK Biobank (n = 1,201,909), we expand the number of eGFRcrea loci (424 loci, 201 novel; 9.8% eGFRcrea variance explained by 634 independent signal variants). Our increased sample size in fine-mapping (n = 1,004,040, European) more than doubles the number of signals with resolved fine-mapping (99% credible sets down to 1 variant for 44 signals, ≤5 variants for 138 signals). Cystatin-based eGFR and/or blood urea nitrogen association support 348 loci (n = 460,826 and 852,678, respectively). Our customizable tool for Gene PrioritiSation reveals 23 compelling genes including mechanistic insights and enables navigation through genes and variants likely relevant for kidney function in human to help select targets for experimental follow-up.

Genetic Susceptibility Loci in Genomewide Association Study of Cluster Headache

Objective

Identifying common genetic variants that confer genetic risk for cluster headache.

Methods

We conducted a case–control study in the Dutch Leiden University Cluster headache neuro‐Analysis program (LUCA) study population (n = 840) and unselected controls from the Netherlands Epidemiology of Obesity Study (NEO; n = 1,457). Replication was performed in a Norwegian sample of 144 cases from the Trondheim Cluster headache sample and 1,800 controls from the Nord‐Trøndelag Health Survey (HUNT). Gene set and tissue enrichment analyses, blood cell‐derived RNA‐sequencing of genes around the risk loci and linkage disequilibrium score regression were part of the downstream analyses.

Results

An association was found with cluster headache for 4 independent loci (r² < 0.1) with genomewide significance (p < 5 × 10⁻⁸), rs11579212 (odds ratio [OR] = 1.51, 95% confidence interval [CI] = 1.33–1.72 near RP11‐815 M8.1), rs6541998 (OR = 1.53, 95% CI = 1.37–1.74 near MERTK), rs10184573 (OR = 1.43, 95% CI = 1.26–1.61 near AC093590.1), and rs2499799 (OR = 0.62, 95% CI = 0.54–0.73 near UFL1/FHL5), collectively explaining 7.2% of the variance of cluster headache. SNPs rs11579212, rs10184573, and rs976357, as proxy SNP for rs2499799 (r² = 1.0), replicated in the Norwegian sample (p < 0.05). Gene‐based mapping yielded ASZ1 as possible fifth locus. RNA‐sequencing indicated differential expression of POLR1B and TMEM87B in cluster headache patients.

Interpretation

This genomewide association study (GWAS) identified and replicated genetic risk loci for cluster headache with effect sizes larger than those typically seen in complex genetic disorders. ANN NEUROL 2021;90:203–216

Review: Mitigating the risks posed by intensification in livestock production: the examples of antimicrobial resistance and zoonoses

Major shifts in how animals are bred, raised and slaughtered are involved in the intensification of livestock systems. Globally, these changes have produced major increases in access to protein-rich foods with high levels of micronutrients. Yet the intensification of livestock systems generates numerous externalities including environmental degradation, zoonotic disease transmission and the emergence of antimicrobial resistance (AMR) genes. Where the process of intensification is most advanced, the expertise, institutions and regulations required to manage these externalities have developed over time, often in response to hard lessons, crises and challenges to public health. By exploring the drivers of intensification, the foci of future intensification can be identified. Low- and middle-income (LMICs) countries are likely to experience significant intensification in livestock production in the near future; however, the lessons learned elsewhere are not being transferred rapidly enough to develop risk mitigation capacity in these settings. At present, fragmentary approaches to address these problems present an incomplete picture of livestock populations, antimicrobial use, and disease risks in LMIC settings. A worldwide improvement in evidence-based zoonotic disease and AMR management within intensifying livestock production systems demands better information on the burden of livestock-associated disease, antimicrobial use and resistance and resources allocated to mitigation.

Mitochondrial genome-wide association study of migraine - the HUNT Study

BACKGROUND

Variation in mitochondrial DNA (mtDNA) has been indicated in migraine pathogenesis, but genetic studies to date have focused on candidate variants, with sparse findings. We aimed to perform the first mitochondrial genome-wide association study of migraine, examining both single variants and mitochondrial haplogroups.

METHODS

In total, 71,860 participants from the population-based Nord-Trøndelag Health Study were genotyped. We excluded samples not passing quality control for nuclear genotypes, in addition to samples with low call rate and closely maternally related. We analysed 775 mitochondrial DNA variants in 4021 migraine cases and 14,288 headache-free controls, using logistic regression. In addition, we analysed 3831 cases and 13,584 controls who could be reliably assigned to a mitochondrial haplogroup. Lastly, we attempted to replicate previously reported mitochondrial DNA candidate variants.

RESULTS

Neither of the mitochondrial variants or haplogroups were associated with migraine. In addition, none of the previously reported mtDNA candidate variants replicated in our data.

CONCLUSIONS

Our findings do not support a major role of mitochondrial genetic variation in migraine pathophysiology, but a larger sample is needed to detect rare variants and future studies should also examine heteroplasmic variation, epigenetic changes and copy-number variation.

Household socio-economic position and individual infectious disease risk in rural Kenya

The importance of household socio-economic position (SEP) in shaping individual infectious disease risk is increasingly recognised, particularly in low income settings. However, few studies have measured the extent to which this association is consistent for the range of pathogens that are typically endemic among the rural poor in the tropics. This cross-sectional study assessed the relationship between SEP and human infection within a single community in western Kenya using a set of pathogens with diverse transmission routes. The relationships between household SEP and individual infection with Plasmodium falciparum, hookworm (Ancylostoma duodenale and/or Necator americanus), Entamoeba histolytica/dispar, Mycobacterium tuberculosis, and HIV, and co-infections between hookworm, P. falciparum and E. histolytica/dispar, were assessed using multivariable logistic and multinomial regression. Individuals in households with the lowest SEP were at greatest risk of infection with P. falciparum, hookworm and E. histolytica/dispar, as well as co-infection with each pathogen. Infection with M. tuberculosis, by contrast, was most likely in individuals living in households with the highest SEP. There was no evidence of a relationship between individual HIV infection and household SEP. We demonstrate the existence of a household socio-economic gradient within a rural farming community in Kenya which impacts upon individual infectious disease risk. Structural adjustments that seek to reduce poverty, and therefore the socio-economic inequalities that exist in this community, would be expected to substantially reduce overall infectious disease burden. However, policy makers and researchers should be aware that heterogeneous relationships can exist between household SEP and infection risk for different pathogens in low income settings.

Environmental predictors of bovine Eimeria infection in western Kenya

Eimeriosis is caused by a protozoan infection affecting most domestic animal species. Outbreaks in cattle are associated with various environmental factors in temperate climates but limited work has been done in tropical settings. The objective of this work was to determine the prevalence and environmental factors associated with bovine Eimeria spp. infection in a mixed farming area of western Kenya. A total of 983 cattle were sampled from 226 cattle-keeping households. Faecal samples were collected directly from the rectum via digital extraction and analysed for the presence of Eimeria spp. infection using the MacMaster technique. Individual and household level predictors of infection were explored using mixed effects logistic regression. The prevalence of individual animal Eimeria infection was 32.8% (95% CI 29.9-35.9). A positive linear relationship was found between risk of Eimeria infection and increasing temperature (OR = 1.4, 95% CI 1.06-1.86) and distance to areas at risk of flooding (OR = 1.49, 95% CI 1.17-1.91). There was weak evidence of non-linear relationship between Eimeria infection and the proportion of the area around a household that was classified as swamp (OR = 1.12, 95% CI 0.87-1.44; OR (quadratic term) = 0.85, 95% CI 0.73-1.00), and the sand content of the soil (OR = 1.18, 95% CI 0.91-1.53; OR (quadratic term) = 1.1, 95% CI 0.99-1.23). The risk of animal Eimeria spp. infection is influenced by a number of climatic and soil-associated conditions.

The importance of p53 pathway genetics in inherited and somatic cancer genomes

Decades of research have shown that mutations in the p53 stress response pathway affect the incidence of diverse cancers more than mutations in other pathways. However, most evidence is limited to somatic mutations and rare inherited mutations. Using newly abundant genomic data, we demonstrate that commonly inherited genetic variants in the p53 pathway also affect the incidence of a broad range of cancers more than variants in other pathways. The cancer-associated single nucleotide polymorphisms (SNPs) of the p53 pathway have strikingly similar genetic characteristics to well-studied p53 pathway cancer-causing somatic mutations. Our results enable insights into p53-mediated tumour suppression in humans and into p53 pathway-based cancer surveillance and treatment strategies.

Rapid Resolution of Donor-Derived Glomerular Fibrin Thrombi After Deceased Donor Kidney Transplantation

The aim of this study was to determine the clinical and histologic outcomes related to transplanting kidneys from deceased donors with glomerular fibrin thrombi (GFT). We included all cases transplanted between October 2003 and October 2014 that had either a preimplantation biopsy or an immediate postreperfusion biopsy showing GFT. The study cohort included 61 recipients (9.9%) with GFT and 557 in the control group without GFT. Delayed graft function occurred in 49% of the GFT group and 39% in the control group (p = 0.14). Serum creatinine at 1, 4, and 12 months and estimated GFR at 12 months were similar in the two groups. Estimated 1-year graft survival was 93.2% in the GFT group and 95.1% in the control group (p = 0.22 by log-rank). Fifty-two of the 61 patients in the GFT group (85%) had a 1-month protocol biopsy, and only two biopsies (4%) showed residual focal glomerular thrombi. At the 1-year protocol biopsy, the prevalence of moderate to severe interstitial fibrosis and tubular atrophy was 24% in the GFT group and 30% in the control group (p = 0.42). We concluded that GFT resolves rapidly after transplantation and that transplanting selected kidneys from deceased donors with GFT is a safe practice.

Circular RNAs are depleted of polymorphisms at microRNA binding sites

Motivation: Circular RNAs (circRNAs) are an abundant class of highly stable RNAs that can affect gene regulation by binding and preventing microRNAs (miRNAs) from regulating their messenger RNA (mRNA) targets. Mammals have thousands of circRNAs with predicted miRNA binding sites, but only two circRNAs have been verified as being actual miRNA sponges. As it is unclear whether these thousands of predicted miRNA binding sites are functional, we investigated whether miRNA seed sites within human circRNAs are under selective pressure.

Results: Using SNP data from the 1000 Genomes Project, we found a significant decrease in SNP density at miRNA seed sites compared with flanking sequences and random sites. This decrease was similar to that of miRNA seed sites in 3' untranslated regions, suggesting that many of the predicted miRNA binding sites in circRNAs are functional and under similar selective pressure as miRNA binding sites in mRNAs.

Contact:pal.satrom@ntnu.no

Supplementary information:Supplementary data are available at Bioinformatics online.

The eGenVar data management system--cataloguing and sharing sensitive data and metadata for the life sciences

Systematic data management and controlled data sharing aim at increasing reproducibility, reducing redundancy in work, and providing a way to efficiently locate complementing or contradicting information. One method of achieving this is collecting data in a central repository or in a location that is part of a federated system and providing interfaces to the data. However, certain data, such as data from biobanks or clinical studies, may, for legal and privacy reasons, often not be stored in public repositories. Instead, we describe a metadata cataloguing system and a software suite for reporting the presence of data from the life sciences domain. The system stores three types of metadata: file information, file provenance and data lineage, and content descriptions. Our software suite includes both graphical and command line interfaces that allow users to report and tag files with these different metadata types. Importantly, the files remain in their original locations with their existing access-control mechanisms in place, while our system provides descriptions of their contents and relationships. Our system and software suite thereby provide a common framework for cataloguing and sharing both public and private data.

Database URL:http://bigr.medisin.ntnu.no/data/eGenVar/

Single nucleotide polymorphisms can create alternative polyadenylation signals and affect gene expression through loss of microRNA-regulation

Alternative polyadenylation (APA) can for example occur when a protein-coding gene has several polyadenylation (polyA) signals in its last exon, resulting in messenger RNAs (mRNAs) with different 3′ untranslated region (UTR) lengths. Different 3′UTR lengths can give different microRNA (miRNA) regulation such that shortened transcripts have increased expression. The APA process is part of human cells' natural regulatory processes, but APA also seems to play an important role in many human diseases. Although altered APA in disease can have many causes, we reasoned that mutations in DNA elements that are important for the polyA process, such as the polyA signal and the downstream GU-rich region, can be one important mechanism. To test this hypothesis, we identified single nucleotide polymorphisms (SNPs) that can create or disrupt APA signals (APA-SNPs). By using a data-integrative approach, we show that APA-SNPs can affect 3′UTR length, miRNA regulation, and mRNA expression—both between homozygote individuals and within heterozygote individuals. Furthermore, we show that a significant fraction of the alleles that cause APA are strongly and positively linked with alleles found by genome-wide studies to be associated with disease. Our results confirm that APA-SNPs can give altered gene regulation and that APA alleles that give shortened transcripts and increased gene expression can be important hereditary causes for disease.

Variants in DNA that affect gene expression—so-called regulatory variants—are thought to play important roles in common complex diseases, such as cancer. In contrast to variants in protein-coding regions, regulatory variants do not affect protein sequence and function. Instead, regulatory variants affect the amount of protein produced. The 3′ untranslated region (UTR) is one gene region that is critically important for gene regulation; cancers for example, often express genes with shortened 3′UTRs that, compared with full-length 3′UTRs, have higher and more stable expression levels. We have investigated one kind of regulatory variant that can affect the 3′UTR length and thereby cause disease. We identified several such variants in different genes and found that these variants affected the genes' expression. Some of these variants were also strongly linked with known markers for disease, suggesting that these regulatory variants are important hereditary causes for disease.

Inferring causative variants in microRNA target sites

MicroRNAs (miRNAs) regulate genes post transcription by pairing with messenger RNA (mRNA). Variants such as single nucleotide polymorphisms (SNPs) in miRNA regulatory regions might result in altered protein levels and disease. Genome-wide association studies (GWAS) aim at identifying genomic regions that contain variants associated with disease, but lack tools for finding causative variants. We present a computational tool that can help identifying SNPs associated with diseases, by focusing on SNPs affecting miRNA-regulation of genes. The tool predicts the effects of SNPs in miRNA target sites and uses linkage disequilibrium to map these miRNA-related variants to SNPs of interest in GWAS. We compared our predicted SNP effects in miRNA target sites with measured SNP effects from allelic imbalance sequencing. Our predictions fit measured effects better than effects based on differences in free energy or differences of TargetScan context scores. We also used our tool to analyse data from published breast cancer and Parkinson's disease GWAS and significant trait-associated SNPs from the NHGRI GWAS Catalog. A database of predicted SNP effects is available at http://www.bigr.medisin.ntnu.no/mirsnpscore/. The database is based on haplotype data from the CEU HapMap population and miRNAs from miRBase 16.0.

A risk variant in an miR-125b binding site in BMPR1B is associated with breast cancer pathogenesis

MicroRNAs regulate diverse cellular processes and play an integral role in cancer pathogenesis. Genomic variation within miRNA target sites may therefore be important sources for genetic differences in cancer risk. To investigate this possibility, we mapped HapMap single nucleotide polymorphisms (SNP) to putative miRNA recognition sites within genes dysregulated in estrogen receptor-stratified breast tumors and used local linkage disequilibrium patterns to identify high-ranking SNPs in the Cancer Genetic Markers of Susceptibility (CGEMS) breast cancer genome-wide association study for further testing. Two SNPs, rs1970801 and rs11097457, scoring in the top 100 from the CGEMS study, were in strong linkage disequilibrium with rs1434536, an SNP that resides within a miR-125b target site in the 3' untranslated region of the bone morphogenic receptor type 1B (BMPR1B) gene encoding a transmembrane serine/threonine kinase. We validated the CGEMS association findings for rs1970801 in an independent cohort of admixture-corrected cases identified from families with multiple case histories. Subsequent association testing of rs1434536 for these cases and CGEMS controls with imputed genotypes supported the association. Furthermore, luciferase reporter assays and overexpression of miR-125b-mimics combined with quantitative reverse transcription-PCR showed that BMPR1B transcript is a direct target of miR-125b and that miR-125b differentially regulates the C and T alleles of rs1434536. These results suggest that allele-specific regulation of BMPR1B by miR-125b explains the observed disease risk. Our approach is general and can help identify and explain the mechanisms behind disease association for alleles that affect miRNA regulation.

Effects of ultraviolet irradiation on the rate and sequence of DNA replication in synchronized Chinese hamster cells

The effects of ultraviolet light (UV) irradiation on the rate of DNA replication in synchronized Chinese hamster ovary (CHO) cells were investigated. A technique for measuring semiconservative DNA replication was employed that involved growing the cells in medium containing 5-bromodeoxyuridine and subsequently determining the amount of DNA that acquired hybrid buoyant density in CsCl density gradients. One of the advantages of this technique was that it allowed a characterization of the extent of DNA replication as well as rate after irradiation. It was found that while there was a dose-dependent reduction in the rate of DNA replication following UV-irradiation, doses of up to 10 J/m2 (which produce many dimers per replication) did not prevent the ultimate replication of the entire genome. Hence, we conclude that dimers cannot be absolute blocks to DNA replication. In order to account for the total genome replication observed, a mechanism must exist that allows genome replication between dimers. The degree of reduction in the rate of replication by UV was the same whether the cells were irradiated at the G1-S boundary or 1 h into S-phase. Previous work had shown that cells in early S-phase are considerably more sensitive to UV than cells at the G1-S boundary. Experiments specifically designed to test for reiterative replication showed that UV does not induce a second round of DNA replication within the same S-phase.