Recent publications

Lashin HMS, Nadkarni S, Oggero S, Jones HR, Knight JC, Hinds CJ, Perretti M. 2017. Microvesicle Subsets in Sepsis Due to Community Acquired Pneumonia Compared to Faecal Peritonitis. Shock, | Show Abstract | Read more

RATIONALE: Microvesicles (MV) act as a non-soluble means of inter-cellular communication, with effector roles in disease pathogenesis and potentially as biomarkers. Previously, we reported that neutrophil MV expressing alpha-2-macroglobulin (A2MG) are protective in experimental sepsis and associate with survival in a small cohort of patients with sepsis due to community acquired pneumonia (CAP). OBJECTIVES: To characterise MV profiles in sepsis due to CAP or faecal peritonitis (FP) and determine their relation to outcome. To investigate the effects of novel sepsis treatments (granulocyte-macrophage colony stimulating factor (GM-CSF) and interferon-υ (IFN-γ)) on MV production and functions in-vitro. METHODS: Flow cytometry analysis of MV identified the cell of origin and the proportion of A2MG expression in the plasma of patients with sepsis secondary to CAP (n = 60) or FP (n = 40) and compared to healthy volunteers (HV, n = 10). The association between MV subsets and outcome was examined. The ability of GM-CSF and IFN-γ on A2MG MV production from whole blood was examined together with the assessment of their effect on neutrophil and endothelial functions. RESULTS: Circulating cell-derived and A2MG MV were higher in CAP compared to FP and HV. A2MG MV were higher in survivors of CAP, but not in FP. GM-CSF and IFN-γ enhanced A2MG MV production, with these MV eliciting pathogen clearance in vitro. CONCLUSIONS: Plasma MV profiles vary according to the source of infection. A2MG MV are associated with survival in CAP but not FP. We propose specific MV subsets as novel biomarkers in sepsis and potential effector for some of the actions of experimental therapeutic interventions.

Burren OS, Rubio García A, Javierre BM, Rainbow DB, Cairns J, Cooper NJ, Lambourne JJ, Schofield E, Castro Dopico X, Ferreira RC et al. 2017. Chromosome contacts in activated T cells identify autoimmune disease candidate genes. Genome Biol, 18 (1), pp. 165. | Show Abstract | Read more

BACKGROUND: Autoimmune disease-associated variants are preferentially found in regulatory regions in immune cells, particularly CD4(+) T cells. Linking such regulatory regions to gene promoters in disease-relevant cell contexts facilitates identification of candidate disease genes. RESULTS: Within 4 h, activation of CD4(+) T cells invokes changes in histone modifications and enhancer RNA transcription that correspond to altered expression of the interacting genes identified by promoter capture Hi-C. By integrating promoter capture Hi-C data with genetic associations for five autoimmune diseases, we prioritised 245 candidate genes with a median distance from peak signal to prioritised gene of 153 kb. Just under half (108/245) prioritised genes related to activation-sensitive interactions. This included IL2RA, where allele-specific expression analyses were consistent with its interaction-mediated regulation, illustrating the utility of the approach. CONCLUSIONS: Our systematic experimental framework offers an alternative approach to candidate causal gene identification for variants with cell state-specific functional effects, with achievable sample sizes.

Grotz AK, Gloyn AL, Thomsen SK. 2017. Prioritising Causal Genes at Type 2 Diabetes Risk Loci. Curr Diab Rep, 17 (9), pp. 76. | Show Abstract | Read more

PURPOSE OF REVIEW: Genome-wide association studies (GWAS) for type 2 diabetes (T2D) risk have identified a large number of genetic loci associated with disease susceptibility. However, progress moving from association signals through causal genes to functional understanding has so far been slow, hindering clinical translation. This review discusses the benefits and limitations of emerging, unbiased approaches for prioritising causal genes at T2D risk loci. RECENT FINDINGS: Candidate causal genes can be identified by a number of different strategies that rely on genetic data, genomic annotations, and functional screening of selected genes. To overcome the limitations of each particular method, integration of multiple data sets is proving essential for establishing confidence in the prioritised genes. Previous studies have also highlighted the need to support these efforts through identification of causal variants and disease-relevant tissues. Prioritisation of causal genes at T2D risk loci by integrating complementary lines of evidence promises to accelerate our understanding of disease pathology and promote translation into new therapeutics.

Scicluna BP, van Vught LA, Zwinderman AH, Wiewel MA, Davenport EE, Burnham KL, Nürnberg P, Schultz MJ, Horn J, Cremer OL et al. 2017. Classification of patients with sepsis according to blood genomic endotype: a prospective cohort study. Lancet Respir Med, | Show Abstract | Read more

BACKGROUND: Host responses during sepsis are highly heterogeneous, which hampers the identification of patients at high risk of mortality and their selection for targeted therapies. In this study, we aimed to identify biologically relevant molecular endotypes in patients with sepsis. METHODS: This was a prospective observational cohort study that included consecutive patients admitted for sepsis to two intensive care units (ICUs) in the Netherlands between Jan 1, 2011, and July 20, 2012 (discovery and first validation cohorts) and patients admitted with sepsis due to community-acquired pneumonia to 29 ICUs in the UK (second validation cohort). We generated genome-wide blood gene expression profiles from admission samples and analysed them by unsupervised consensus clustering and machine learning. The primary objective of this study was to establish endotypes for patients with sepsis, and assess the association of these endotypes with clinical traits and survival outcomes. We also established candidate biomarkers for the endotypes to allow identification of patient endotypes in clinical practice. FINDINGS: The discovery cohort had 306 patients, the first validation cohort had 216, and the second validation cohort had 265 patients. Four molecular endotypes for sepsis, designated Mars1-4, were identified in the discovery cohort, and were associated with 28-day mortality (log-rank p=0·022). In the discovery cohort, the worst outcome was found for patients classified as having a Mars1 endotype, and at 28 days, 35 (39%) of 90 people with a Mars1 endotype had died (hazard ratio [HR] vs all other endotypes 1·86 [95% CI 1·21-2·86]; p=0·0045), compared with 23 (22%) of 105 people with a Mars2 endotype (HR 0·64 [0·40-1·04]; p=0·061), 16 (23%) of 71 people with a Mars3 endotype (HR 0·71 [0·41-1·22]; p=0·19), and 13 (33%) of 40 patients with a Mars4 endotype (HR 1·13 [0·63-2·04]; p=0·69). Analysis of the net reclassification improvement using a combined clinical and endotype model significantly improved risk prediction to 0·33 (0·09-0·58; p=0·008). A 140-gene expression signature reliably stratified patients with sepsis to the four endotypes in both the first and second validation cohorts. Only Mars1 was consistently significantly associated with 28-day mortality across the cohorts. To facilitate possible clinical use, a biomarker was derived for each endotype; BPGM and TAP2 reliably identified patients with a Mars1 endotype. INTERPRETATION: This study provides a method for the molecular classification of patients with sepsis to four different endotypes upon ICU admission. Detection of sepsis endotypes might assist in providing personalised patient management and in selection for trials. FUNDING: Center for Translational Molecular Medicine, Netherlands.

Pekalski ML, García AR, Ferreira RC, Rainbow DB, Smyth DJ, Mashar M, Brady J, Savinykh N, Dopico XC, Mahmood S et al. 2017. Neonatal and adult recent thymic emigrants produce IL-8 and express complement receptors CR1 and CR2. JCI Insight, 2 (16), | Show Abstract | Read more

The maintenance of peripheral naive T lymphocytes in humans is dependent on their homeostatic division, not continuing emigration from the thymus, which undergoes involution with age. However, postthymic maintenance of naive T cells is still poorly understood. Previously we reported that recent thymic emigrants (RTEs) are contained in CD31+CD25- naive T cells as defined by their levels of signal joint T cell receptor rearrangement excision circles (sjTRECs). Here, by differential gene expression analysis followed by protein expression and functional studies, we define that the naive T cells having divided the least since thymic emigration express complement receptors (CR1 and CR2) known to bind complement C3b- and C3d-decorated microbial products and, following activation, produce IL-8 (CXCL8), a major chemoattractant for neutrophils in bacterial defense. We also observed an IL-8-producing memory T cell subpopulation coexpressing CR1 and CR2 and with a gene expression signature resembling that of RTEs. The functions of CR1 and CR2 on T cells remain to be determined, but we note that CR2 is the receptor for Epstein-Barr virus, which is a cause of T cell lymphomas and a candidate environmental factor in autoimmune disease.

Milicic A, S Rollier C, Tang CK, Longley R, Hill AVS, Reyes-Sandoval A. 2017. Adjuvanting a viral vectored vaccine against pre-erythrocytic malaria. Sci Rep, 7 (1), pp. 7284. | Show Abstract | Read more

The majority of routinely given vaccines require two or three immunisations for full protective efficacy. Single dose vaccination has long been considered a key solution to improving the global immunisation coverage. Recent infectious disease outbreaks have further highlighted the need for vaccines that can achieve full efficacy after a single administration. Viral vectors are a potent immunisation platform, benefiting from intrinsic immuno-stimulatory features while retaining excellent safety profile through the use of non-replicating viruses. We investigated the scope for enhancing the protective efficacy of a single dose adenovirus-vectored malaria vaccine in a mouse model of malaria by co-administering it with vaccine adjuvants. Out of 11 adjuvants, only two, Abisco(®)-100 and CoVaccineHT(TM), enhanced vaccine efficacy and sterile protection following malaria challenge. The CoVaccineHT(TM) adjuvanted vaccine induced significantly higher proportion of antigen specific central memory CD8(+) cells, and both adjuvants resulted in increased proportion of CD8(+) T cells expressing the CD107a degranulation marker in the absence of IFNγ, TNFα and IL2 production. Our results show that the efficacy of vaccines designed to induce protective T cell responses can be positively modulated with chemical adjuvants and open the possibility of achieving full protection with a single dose immunisation.

Venkatraman N, Silman D, Folegatti PM, Hill AVS. 2017. Vaccines against Ebola virus. Vaccine, | Show Abstract | Read more

We have just witnessed the largest and most devastating outbreak of Ebola virus disease, which highlighted the urgent need for development of an efficacious vaccine that could be used to curtail future outbreaks. Prior to 2014, there had been limited impetus worldwide to develop a vaccine since the virus was first discovered in 1976. Though too many lives were lost during this outbreak, it resulted in the significantly accelerated clinical development of a number of candidate vaccines through an extraordinary collaborative global effort coordinated by the World Health Organisation (WHO) and involving a number of companies, trial centres, funders, global stakeholders and agencies. We have acquired substantial safety and immunogenicity data on a number of vaccines in Caucasian and African populations. The rapid pace of events led to the initiation of the landmark efficacy trial testing the rVSV-vectored vaccine, which showed high level efficacy in an outbreak setting when deployed using an innovative ring vaccination strategy. Though the Public Health Emergency of International Concern (PHEIC) declared by the WHO has now been lifted, the global scientific community faces numerous challenges ahead to ensure that there is a licensed, deployable vaccine available for use in future outbreaks for at least the Zaire and Sudan strains of Ebola virus. There remain several unanswered questions on the durability of protection, mechanistic immunological correlates and preferred deployment strategies. This review outlines a brief history of the development of Ebola vaccines, the significant progress made since the scale of the outbreak became apparent, some lessons learnt and how they could shape future development of vaccines and the management of similar outbreaks.

Censin JC, Nowak C, Cooper N, Bergsten P, Todd JA, Fall T. 2017. Childhood adiposity and risk of type 1 diabetes: A Mendelian randomization study. PLoS Med, 14 (8), pp. e1002362. | Show Abstract | Read more

BACKGROUND: The incidence of type 1 diabetes (T1D) is increasing globally. One hypothesis is that increasing childhood obesity rates may explain part of this increase, but, as T1D is rare, intervention studies are challenging to perform. The aim of this study was to assess this hypothesis with a Mendelian randomization approach that uses genetic variants as instrumental variables to test for causal associations. METHODS AND FINDINGS: We created a genetic instrument of 23 single nucleotide polymorphisms (SNPs) associated with childhood adiposity in children aged 2-10 years. Summary-level association results for these 23 SNPs with childhood-onset (<17 years) T1D were extracted from a meta-analysis of genome-wide association study with 5,913 T1D cases and 8,828 reference samples. Using inverse-variance weighted Mendelian randomization analysis, we found support for an effect of childhood adiposity on T1D risk (odds ratio 1.32, 95% CI 1.06-1.64 per standard deviation score in body mass index [SDS-BMI]). A sensitivity analysis provided evidence of horizontal pleiotropy bias (p = 0.04) diluting the estimates towards the null. We therefore applied Egger regression and multivariable Mendelian randomization methods to control for this type of bias and found evidence in support of a role of childhood adiposity in T1D (odds ratio in Egger regression, 2.76, 95% CI 1.40-5.44). Limitations of our study include that underlying genes and their mechanisms for most of the genetic variants included in the score are not known. Mendelian randomization requires large sample sizes, and power was limited to provide precise estimates. This research has been conducted using data from the Early Growth Genetics (EGG) Consortium, the Genetic Investigation of Anthropometric Traits (GIANT) Consortium, the Tobacco and Genetics (TAG) Consortium, and the Social Science Genetic Association Consortium (SSGAC), as well as meta-analysis results from a T1D genome-wide association study. CONCLUSIONS: This study provides genetic support for a link between childhood adiposity and T1D risk. Together with evidence from observational studies, our findings further emphasize the importance of measures to reduce the global epidemic of childhood obesity and encourage mechanistic studies.

Yates LR, Knappskog S, Wedge D, Farmery JHR, Gonzalez S, Martincorena I, Alexandrov LB, Van Loo P, Haugland HK, Lilleng PK et al. 2017. Genomic Evolution of Breast Cancer Metastasis and Relapse. Cancer Cell, 32 (2), pp. 169-184.e7. | Citations: 1 (Web of Science Lite) | Show Abstract | Read more

Patterns of genomic evolution between primary and metastatic breast cancer have not been studied in large numbers, despite patients with metastatic breast cancer having dismal survival. We sequenced whole genomes or a panel of 365 genes on 299 samples from 170 patients with locally relapsed or metastatic breast cancer. Several lines of analysis indicate that clones seeding metastasis or relapse disseminate late from primary tumors, but continue to acquire mutations, mostly accessing the same mutational processes active in the primary tumor. Most distant metastases acquired driver mutations not seen in the primary tumor, drawing from a wider repertoire of cancer genes than early drivers. These include a number of clinically actionable alterations and mutations inactivating SWI-SNF and JAK2-STAT3 pathways.

Buchanan J, Wordsworth S, Clifford R, Robbe P, Taylor JC, Schuh A, Knight SJL. 2017. Using Genomic Information to Guide Ibrutinib Treatment Decisions in Chronic Lymphocytic Leukaemia: A Cost-Effectiveness Analysis. Pharmacoeconomics, 35 (8), pp. 845-858. | Show Abstract | Read more

BACKGROUND: Genomic tests may improve the stratification of patients to receive new therapies in several disease areas. However, the use of expensive targeted therapies can impact on the cost effectiveness of these tests. This study presents an economic evaluation of genomic testing in chronic lymphocytic leukaemia in the context of the UK National Health Service. METHODS: Cost-effectiveness and cost-utility analyses (using life-years and quality-adjusted life-years) were undertaken from a National Health Service and societal perspective. Five strategies were evaluated across several age groups using Markov modelling: three strategies that reflected varying current genetic testing practice and two configurations of genomic testing (including ibrutinib treatment). RESULTS: Genomic testing strategies yielded the most life-years/quality-adjusted life-years per patient, but were not cost effective compared with a threshold of £30,000 per life-year/quality-adjusted life-year gained. Cost-effectiveness acceptability curves indicated that there was some uncertainty surrounding this result. A genomic testing strategy becomes the most cost-effective option if a higher end-of-life cost-effectiveness threshold of £50,000 is applied, if a societal costing perspective is considered in 25-year-old patients or if the cost of ibrutinib treatment falls. CONCLUSION: Stratifying patients with chronic lymphocytic leukaemia to targeted treatment using genomic testing improves health outcomes, but will likely only represent a cost-effective use of limited National Health Service resources if a higher cost-effectiveness threshold or societal costing perspective is applied, or if the price of ibrutinib treatment is reduced. This result may be broadly indicative of the likely cost effectiveness of other genomic tests that inform the stratification of patients to high cost-targeted therapies.

Cortes A, Dendrou CA, Motyer A, Jostins L, Vukcevic D, Dilthey A, Donnelly P, Leslie S, Fugger L, McVean G. 2017. Bayesian analysis of genetic association across tree-structured routine healthcare data in the UK Biobank. Nat Genet, 49 (9), pp. 1311-1318. | Citations: 1 (Web of Science Lite) | Show Abstract | Read more

Genetic discovery from the multitude of phenotypes extractable from routine healthcare data can transform understanding of the human phenome and accelerate progress toward precision medicine. However, a critical question when analyzing high-dimensional and heterogeneous data is how best to interrogate increasingly specific subphenotypes while retaining statistical power to detect genetic associations. Here we develop and employ a new Bayesian analysis framework that exploits the hierarchical structure of diagnosis classifications to analyze genetic variants against UK Biobank disease phenotypes derived from self-reporting and hospital episode statistics. Our method displays a more than 20% increase in power to detect genetic effects over other approaches and identifies new associations between classical human leukocyte antigen (HLA) alleles and common immune-mediated diseases (IMDs). By applying the approach to genetic risk scores (GRSs), we show the extent of genetic sharing among IMDs and expose differences in disease perception or diagnosis with potential clinical implications.

Vieira JM, Howard S, Villa Del Campo C, Bollini S, Dubé KN, Masters M, Barnette DN, Rohling M, Sun X, Hankins LE et al. 2017. BRG1-SWI/SNF-dependent regulation of the Wt1 transcriptional landscape mediates epicardial activity during heart development and disease. Nat Commun, 8 pp. 16034. | Show Abstract | Read more

Epicardium-derived cells (EPDCs) contribute cardiovascular cell types during development and in adulthood respond to Thymosin β4 (Tβ4) and myocardial infarction (MI) by reactivating a fetal gene programme to promote neovascularization and cardiomyogenesis. The mechanism for epicardial gene (re-)activation remains elusive. Here we reveal that BRG1, the essential ATPase subunit of the SWI/SNF chromatin-remodelling complex, is required for expression of Wilms' tumour 1 (Wt1), fetal EPDC activation and subsequent differentiation into coronary smooth muscle, and restores Wt1 activity upon MI. BRG1 physically interacts with Tβ4 and is recruited by CCAAT/enhancer-binding protein β (C/EBPβ) to discrete regulatory elements in the Wt1 locus. BRG1-Tβ4 co-operative binding promotes optimal transcription of Wt1 as the master regulator of embryonic EPDCs. Moreover, chromatin immunoprecipitation-sequencing reveals BRG1 binding at further key loci suggesting SWI/SNF activity across the fetal epicardial gene programme. These findings reveal essential functions for chromatin-remodelling in the activation of EPDCs during cardiovascular development and repair.

Hanssen LLP, Kassouf MT, Oudelaar AM, Biggs D, Preece C, Downes DJ, Gosden M, Sharpe JA, Sloane-Stanley JA, Hughes JR et al. 2017. Tissue-specific CTCF-cohesin-mediated chromatin architecture delimits enhancer interactions and function in vivo. Nat Cell Biol, 19 (8), pp. 952-961. | Citations: 1 (Web of Science Lite) | Show Abstract | Read more

The genome is organized via CTCF-cohesin-binding sites, which partition chromosomes into 1-5 megabase (Mb) topologically associated domains (TADs), and further into smaller sub-domains (sub-TADs). Here we examined in vivo an ∼80 kb sub-TAD, containing the mouse α-globin gene cluster, lying within a ∼1 Mb TAD. We find that the sub-TAD is flanked by predominantly convergent CTCF-cohesin sites that are ubiquitously bound by CTCF but only interact during erythropoiesis, defining a self-interacting erythroid compartment. Whereas the α-globin regulatory elements normally act solely on promoters downstream of the enhancers, removal of a conserved upstream CTCF-cohesin boundary extends the sub-TAD to adjacent upstream CTCF-cohesin-binding sites. The α-globin enhancers now interact with the flanking chromatin, upregulating expression of genes within this extended sub-TAD. Rather than acting solely as a barrier to chromatin modification, CTCF-cohesin boundaries in this sub-TAD delimit the region of chromatin to which enhancers have access and within which they interact with receptive promoters.

Wain LV, Vaez A, Jansen R, Joehanes R, van der Most PJ, Erzurumluoglu AM, O'Reilly PF, Cabrera CP, Warren HR, Rose LM et al. 2017. Novel Blood Pressure Locus and Gene Discovery Using Genome-Wide Association Study and Expression Data Sets From Blood and the Kidney. Hypertension, 70 (3), pp. E4-+. | Show Abstract | Read more

Elevated blood pressure is a major risk factor for cardiovascular disease and has a substantial genetic contribution. Genetic variation influencing blood pressure has the potential to identify new pharmacological targets for the treatment of hypertension. To discover additional novel blood pressure loci, we used 1000 Genomes Project-based imputation in 150 134 European ancestry individuals and sought significant evidence for independent replication in a further 228 245 individuals. We report 6 new signals of association in or near HSPB7, TNXB, LRP12, LOC283335, SEPT9, and AKT2, and provide new replication evidence for a further 2 signals in EBF2 and NFKBIA Combining large whole-blood gene expression resources totaling 12 607 individuals, we investigated all novel and previously reported signals and identified 48 genes with evidence for involvement in blood pressure regulation that are significant in multiple resources. Three novel kidney-specific signals were also detected. These robustly implicated genes may provide new leads for therapeutic innovation.


© Institute of Mathematical Statistics, 2017. The combination of genetic information with electronic patient records promises to provide a powerful new resource for understanding human disease and its treatment. Here we develop and apply a novel stochastic compartmental model to a large dataset on Clostridium difficile infection (CDI) in three Oxfordshire hospitals over a 2.5 year period which combines genetic information on 858 confirmed cases of CDI with a database of 750,000 patient records. C. difficile is a major cause of healthcare-associated diarrhoea and is responsible for substantial mortality and morbidity, with relatively little known about its biology or its transmission epidemiology. Bayesian analysis of our model, via Markov chain Monte Carlo, provides new information about the biology of CDI, including genetic heterogeneity in infectiousness across different sequence types, and evidence for ward contamination as a significant mode of transmission, and allows inferences about the contribution of particular individuals, wards or hospitals to transmission of the bacterium, and assessment of changes in these over time following changes in hospital practice. Our work demonstrates the value of using statistical modelling and computational inference on large-scale hospital patient databases and genetic data.

Zillikens MC, Demissie S, Hsu YH, Yerges-Armstrong LM, Chou WC, Stolk L, Livshits G, Broer L, Johnson T, Koller DL et al. 2017. Large meta-analysis of genome-wide association studies identifies five loci for lean body mass. Nat Commun, 8 (1), pp. 80. | Show Abstract | Read more

Lean body mass, consisting mostly of skeletal muscle, is important for healthy aging. We performed a genome-wide association study for whole body (20 cohorts of European ancestry with n = 38,292) and appendicular (arms and legs) lean body mass (n = 28,330) measured using dual energy X-ray absorptiometry or bioelectrical impedance analysis, adjusted for sex, age, height, and fat mass. Twenty-one single-nucleotide polymorphisms were significantly associated with lean body mass either genome wide (p < 5 × 10(-8)) or suggestively genome wide (p < 2.3 × 10(-6)). Replication in 63,475 (47,227 of European ancestry) individuals from 33 cohorts for whole body lean body mass and in 45,090 (42,360 of European ancestry) subjects from 25 cohorts for appendicular lean body mass was successful for five single-nucleotide polymorphisms in/near HSD17B11, VCAN, ADAMTSL3, IRS1, and FTO for total lean body mass and for three single-nucleotide polymorphisms in/near VCAN, ADAMTSL3, and IRS1 for appendicular lean body mass. Our findings provide new insight into the genetics of lean body mass.Lean body mass is a highly heritable trait and is associated with various health conditions. Here, Kiel and colleagues perform a meta-analysis of genome-wide association studies for whole body lean body mass and find five novel genetic loci to be significantly associated.

Kaur G, Gras S, Mobbs JI, Vivian JP, Cortes A, Barber T, Kuttikkatte SB, Jensen LT, Attfield KE, Dendrou CA et al. 2017. Structural and regulatory diversity shape HLA-C protein expression levels. Nat Commun, 8 pp. 15924. | Show Abstract | Read more

Expression of HLA-C varies widely across individuals in an allele-specific manner. This variation in expression can influence efficacy of the immune response, as shown for infectious and autoimmune diseases. MicroRNA binding partially influences differential HLA-C expression, but the additional contributing factors have remained undetermined. Here we use functional and structural analyses to demonstrate that HLA-C expression is modulated not just at the RNA level, but also at the protein level. Specifically, we show that variation in exons 2 and 3, which encode the α1/α2 domains, drives differential expression of HLA-C allomorphs at the cell surface by influencing the structure of the peptide-binding cleft and the diversity of peptides bound by the HLA-C molecules. Together with a phylogenetic analysis, these results highlight the diversity and long-term balancing selection of regulatory factors that modulate HLA-C expression.

Behjati S, Tarpey PS, Haase K, Ye H, Young MD, Alexandrov LB, Farndon SJ, Collord G, Wedge DC, Martincorena I et al. 2017. Recurrent mutation of IGF signalling genes and distinct patterns of genomic rearrangement in osteosarcoma. Nat Commun, 8 pp. 15936. | Citations: 1 (Web of Science Lite) | Show Abstract | Read more

Osteosarcoma is a primary malignancy of bone that affects children and adults. Here, we present the largest sequencing study of osteosarcoma to date, comprising 112 childhood and adult tumours encompassing all major histological subtypes. A key finding of our study is the identification of mutations in insulin-like growth factor (IGF) signalling genes in 8/112 (7%) of cases. We validate this observation using fluorescence in situ hybridization (FISH) in an additional 87 osteosarcomas, with IGF1 receptor (IGF1R) amplification observed in 14% of tumours. These findings may inform patient selection in future trials of IGF1R inhibitors in osteosarcoma. Analysing patterns of mutation, we identify distinct rearrangement profiles including a process characterized by chromothripsis and amplification. This process operates recurrently at discrete genomic regions and generates driver mutations. It may represent an age-independent mutational mechanism that contributes to the development of osteosarcoma in children and adults alike.

Saare M, Rekker K, Laisk-Podar T, Rahmioglu N, Zondervan K, Salumets A, Götte M, Peters M. 2017. Challenges in endometriosis miRNA studies - From tissue heterogeneity to disease specific miRNAs. Biochim Biophys Acta, 1863 (9), pp. 2282-2292. | Show Abstract | Read more

In order to uncover miRNA changes in endometriosis pathogenesis, both endometriotic lesions and endometrial biopsies, as well as stromal and epithelial cells isolated from these tissues have been investigated and a large number of dysregulated miRNAs have been reported. However, the concordance between the result of different studies has remained small. One potential explanation for limited overlap between the proposed disease-related miRNAs could be the heterogeneity in tissue composition, as some studies have compared highly heterogeneous whole-lesion biopsies with endometrial tissue, some have compared the endometrium from patients and controls, and some have used pure cell fractions isolated from lesions and endometrium. This review focuses on the results of published miRNA studies in endometriosis to reveal the potential impact of tissue heterogeneity on the discovery of disease-specific miRNA alterations in endometriosis. Additionally, functional studies that explore the roles of endometriosis-involved miRNAs are discussed.

Cebrian-Serrano A, Davies B. 2017. CRISPR-Cas orthologues and variants: optimizing the repertoire, specificity and delivery of genome engineering tools. Mamm Genome, 28 (7-8), pp. 247-261. | Show Abstract | Read more

Robust and cost-effective genome editing in a diverse array of cells and model organisms is now possible thanks to the discovery of the RNA-guided endonucleases of the CRISPR-Cas system. The commonly used Cas9 of Streptococcus pyogenes shows high levels of activity but, depending on the application, has been associated with some shortcomings. Firstly, the enzyme has been shown to cause mutagenesis at genomic sequences resembling the target sequence. Secondly, the stringent requirement for a specific motif adjacent to the selected target site can limit the target range of this enzyme. Lastly, the physical size of Cas9 challenges the efficient delivery of genomic engineering tools based on this enzyme as viral particles for potential therapeutic applications. Related and parallel strategies have been employed to address these issues. Taking advantage of the wealth of structural information that is becoming available for CRISPR-Cas effector proteins, Cas9 has been redesigned by mutagenizing key residues contributing to activity and target recognition. The protein has also been shortened and redesigned into component subunits in an attempt to facilitate its efficient delivery. Furthermore, the CRISPR-Cas toolbox has been expanded by exploring the properties of Cas9 orthologues and other related effector proteins from diverse bacterial species, some of which exhibit different target site specificities and reduced molecular size. It is hoped that the improvements in accuracy, target range and efficiency of delivery will facilitate the therapeutic application of these site-specific nucleases.

Kukushkina V, Modhukur V, Suhorutšenko M, Peters M, Mägi R, Rahmioglu N, Velthut-Meikas A, Altmäe S, Esteban FJ, Vilo J et al. 2017. DNA methylation changes in endometrium and correlation with gene expression during the transition from pre-receptive to receptive phase. Sci Rep, 7 (1), pp. 3916. | Show Abstract | Read more

The inner uterine lining (endometrium) is a unique tissue going through remarkable changes each menstrual cycle. Endometrium has its characteristic DNA methylation profile, although not much is known about the endometrial methylome changes throughout the menstrual cycle. The impact of methylome changes on gene expression and thereby on the function of the tissue, including establishing receptivity to implanting embryo, is also unclear. Therefore, this study used genome-wide technologies to characterize the methylome and the correlation between DNA methylation and gene expression in endometrial biopsies collected from 17 healthy fertile-aged women from pre-receptive and receptive phase within one menstrual cycle. Our study showed that the overall methylome remains relatively stable during this stage of the menstrual cycle, with small-scale changes affecting 5% of the studied CpG sites (22,272 out of studied 437,022 CpGs, FDR < 0.05). Of differentially methylated CpG sites with the largest absolute changes in methylation level, approximately 30% correlated with gene expression measured by RNA sequencing, with negative correlations being more common in 5' UTR and positive correlations in the gene 'Body' region. According to our results, extracellular matrix organization and immune response are the pathways most affected by methylation changes during the transition from pre-receptive to receptive phase.

Leffler EM, Band G, Busby GBJ, Kivinen K, Le QS, Clarke GM, Bojang KA, Conway DJ, Jallow M, Sisay-Joof F et al. 2017. Resistance to malaria through structural variation of red blood cell invasion receptors SCIENCE, 356 (6343), | Citations: 2 (Web of Science Lite) | Read more

Payne RO, Silk SE, Elias SC, Milne KH, Rawlinson TA, Llewellyn D, Shakri AR, Jin J, Labbé GM, Edwards NJ et al. 2017. Human vaccination against Plasmodium vivax Duffy-binding protein induces strain-transcending antibodies. JCI Insight, 2 (12), | Show Abstract | Read more

BACKGROUND: Plasmodium vivax is the most widespread human malaria geographically; however, no effective vaccine exists. Red blood cell invasion by the P. vivax merozoite depends on an interaction between the Duffy antigen receptor for chemokines (DARC) and region II of the parasite's Duffy-binding protein (PvDBP_RII). Naturally acquired binding-inhibitory antibodies against this interaction associate with clinical immunity, but it is unknown whether these responses can be induced by human vaccination. METHODS: Safety and immunogenicity of replication-deficient chimpanzee adenovirus serotype 63 (ChAd63) and modified vaccinia virus Ankara (MVA) viral vectored vaccines targeting PvDBP_RII (Salvador I strain) were assessed in an open-label dose-escalation phase Ia study in 24 healthy UK adults. Vaccines were delivered by the intramuscular route in a ChAd63-MVA heterologous prime-boost regimen using an 8-week interval. RESULTS: Both vaccines were well tolerated and demonstrated a favorable safety profile in malaria-naive adults. PvDBP_RII-specific ex-vivo IFN-γ T cell, antibody-secreting cell, memory B cell, and serum IgG responses were observed after the MVA boost immunization. Vaccine-induced antibodies inhibited the binding of vaccine homologous and heterologous variants of recombinant PvDBP_RII to the DARC receptor, with median 50% binding-inhibition titers greater than 1:100. CONCLUSION: We have demonstrated for the first time to our knowledge that strain-transcending antibodies can be induced against the PvDBP_RII antigen by vaccination in humans. These vaccine candidates warrant further clinical evaluation of efficacy against the blood-stage P. vivax parasite. TRIAL REGISTRATION: NCT01816113. FUNDING: Support was provided by the UK Medical Research Council, UK National Institute of Health Research Oxford Biomedical Research Centre, and the Wellcome Trust.

Nolte IM, Munoz ML, Tragante V, Amare AT, Jansen R, Vaez A, von der Heyde B, Avery CL, Bis JC, Dierckx B et al. 2017. Genetic loci associated with heart rate variability and their effects on cardiac disease risk. Nat Commun, 8 pp. 15805. | Citations: 1 (Web of Science Lite) | Show Abstract | Read more

Reduced cardiac vagal control reflected in low heart rate variability (HRV) is associated with greater risks for cardiac morbidity and mortality. In two-stage meta-analyses of genome-wide association studies for three HRV traits in up to 53,174 individuals of European ancestry, we detect 17 genome-wide significant SNPs in eight loci. HRV SNPs tag non-synonymous SNPs (in NDUFA11 and KIAA1755), expression quantitative trait loci (eQTLs) (influencing GNG11, RGS6 and NEO1), or are located in genes preferentially expressed in the sinoatrial node (GNG11, RGS6 and HCN4). Genetic risk scores account for 0.9 to 2.6% of the HRV variance. Significant genetic correlation is found for HRV with heart rate (-0.74<rg<-0.55) and blood pressure (-0.35<rg<-0.20). These findings provide clinically relevant biological insight into heritable variation in vagal heart rhythm regulation, with a key role for genetic variants (GNG11, RGS6) that influence G-protein heterotrimer action in GIRK-channel induced pacemaker membrane hyperpolarization.

Shungin D, Deng WQ, Varga TV, Luan J, Mihailov E, Metspalu A, GIANT Consortium, Morris AP, Forouhi NG, Lindgren C et al. 2017. Ranking and characterization of established BMI and lipid associated loci as candidates for gene-environment interactions. PLoS Genet, 13 (6), pp. e1006812. | Show Abstract | Read more

Phenotypic variance heterogeneity across genotypes at a single nucleotide polymorphism (SNP) may reflect underlying gene-environment (G×E) or gene-gene interactions. We modeled variance heterogeneity for blood lipids and BMI in up to 44,211 participants and investigated relationships between variance effects (Pv), G×E interaction effects (with smoking and physical activity), and marginal genetic effects (Pm). Correlations between Pv and Pm were stronger for SNPs with established marginal effects (Spearman's ρ = 0.401 for triglycerides, and ρ = 0.236 for BMI) compared to all SNPs. When Pv and Pm were compared for all pruned SNPs, only BMI was statistically significant (Spearman's ρ = 0.010). Overall, SNPs with established marginal effects were overrepresented in the nominally significant part of the Pv distribution (Pbinomial <0.05). SNPs from the top 1% of the Pm distribution for BMI had more significant Pv values (PMann-Whitney = 1.46×10-5), and the odds ratio of SNPs with nominally significant (<0.05) Pm and Pv was 1.33 (95% CI: 1.12, 1.57) for BMI. Moreover, BMI SNPs with nominally significant G×E interaction P-values (Pint<0.05) were enriched with nominally significant Pv values (Pbinomial = 8.63×10-9 and 8.52×10-7 for SNP × smoking and SNP × physical activity, respectively). We conclude that some loci with strong marginal effects may be good candidates for G×E, and variance-based prioritization can be used to identify them.

Muller J, Parizotto E, Antrobus R, Francis J, Bunce C, Stranks A, Nichols M, McClain M, Hill AVS, Ramasamy A, Gilbert SC. 2017. Development of an objective gene expression panel as an alternative to self-reported symptom scores in human influenza challenge trials. J Transl Med, 15 (1), pp. 134. | Show Abstract | Read more

BACKGROUND: Influenza challenge trials are important for vaccine efficacy testing. Currently, disease severity is determined by self-reported scores to a list of symptoms which can be highly subjective. A more objective measure would allow for improved data analysis. METHODS: Twenty-one volunteers participated in an influenza challenge trial. We calculated the daily sum of scores (DSS) for a list of 16 influenza symptoms. Whole blood collected at baseline and 24, 48, 72 and 96 h post challenge was profiled on Illumina HT12v4 microarrays. Changes in gene expression most strongly correlated with DSS were selected to train a Random Forest model and tested on two independent test sets consisting of 41 individuals profiled on a different microarray platform and 33 volunteers assayed by qRT-PCR. RESULTS: 1456 probes are significantly associated with DSS at 1% false discovery rate. We selected 19 genes with the largest fold change to train a random forest model. We observed good concordance between predicted and actual scores in the first test set (r = 0.57; RMSE = -16.1%) with the greatest agreement achieved on samples collected approximately 72 h post challenge. Therefore, we assayed samples collected at baseline and 72 h post challenge in the second test set by qRT-PCR and observed good concordance (r = 0.81; RMSE = -36.1%). CONCLUSIONS: We developed a 19-gene qRT-PCR panel to predict DSS, validated on two independent datasets. A transcriptomics based panel could provide a more objective measure of symptom scoring in future influenza challenge studies. Trial registration Samples were obtained from a clinical trial with the Identifier: NCT02014870, first registered on December 5, 2013.

Burns A, Alsolami R, Becq J, Timbs A, Bruce D, Robbe P, Vavoulis D, Cabes M, Dreau H, Taylor J et al. 2017. Whole-genome sequencing of chronic lymphocytic leukaemia reveals distinct differences in the mutational landscape between IgHV(mut) and IgHV(unmut) subgroups. Leukemia, | Show Abstract | Read more

Chronic lymphocytic leukaemia (CLL) consists of two biologically and clinically distinct subtypes defined by the abundance of somatic hypermutation (SHM) affecting the Ig variable heavy-chain locus (IgHV). The molecular mechanisms underlying these subtypes are incompletely understood. Here, we present a comprehensive whole-genome sequencing analysis of somatically acquired genetic events from 46 CLL patients, including a systematic comparison of coding and non-coding single-nucleotide variants, copy number variants and structural variants, regions of kataegis and mutation signatures between IgHV(mut) and IgHV(unmut) subtypes. We demonstrate that one-quarter of non-coding mutations in regions of kataegis outside the Ig loci are located in genes relevant to CLL. We show that non-coding mutations in ATM may negatively impact on ATM expression and find non-coding and regulatory region mutations in TCL1A, and in IgHV(unmut) CLL in IKZF3, SAMHD1,PAX5 and BIRC3. Finally, we show that IgHV(unmut) CLL is dominated by coding mutations in driver genes and an aging signature, whereas IgHV(mut) CLL has a high incidence of promoter and enhancer mutations caused by aberrant activation-induced cytidine deaminase activity. Taken together, our data support the hypothesis that differences in clinical outcome and biological characteristics between the two subgroups might reflect differences in mutation distribution, incidence and distinct underlying mutagenic mechanisms.Leukemia advance online publication, 27 June 2017; doi:10.1038/leu.2017.177.

Rolo J, Worning P, Nielsen JB, Bowden R, Bouchami O, Damborg P, Guardabassi L, Perreten V, Tomasz A, Westh H et al. 2017. Evolutionary Origin of the Staphylococcal Cassette Chromosome mec (SCCmec). Antimicrob Agents Chemother, 61 (6), pp. e02302-16-e02302-16. | Show Abstract | Read more

Several lines of evidence indicate that the most primitive staphylococcal species, those of the Staphylococcus sciuri group, were involved in the first stages of evolution of the staphylococcal cassette chromosome mec (SCCmec), the genetic element carrying the β-lactam resistance gene mecA However, many steps are still missing from this evolutionary history. In particular, it is not known how mecA was incorporated into the mobile element SCC prior to dissemination among Staphylococcus aureus and other pathogenic staphylococcal species. To gain insights into the possible contribution of several species of the Staphylococcus sciuri group to the assembly of SCCmec, we sequenced the genomes of 106 isolates, comprising S. sciuri (n = 76), Staphylococcus vitulinus (n = 18), and Staphylococcus fleurettii (n = 12) from animal and human sources, and characterized the native location of mecA and the SCC insertion site by using a variety of comparative genomic approaches. Moreover, we performed a single nucleotide polymorphism (SNP) analysis of the genomes in order to understand SCCmec evolution in relation to phylogeny. We found that each of three species of the S. sciuri group contributed to the evolution of SCCmec: S. vitulinus and S. fleurettii contributed to the assembly of the mec complex, and S. sciuri most likely provided the mobile element in which mecA was later incorporated. We hypothesize that an ancestral SCCmec III cassette (an element carried by one of the most epidemic methicillin-resistant S. aureus clones) originated in S. sciuri possibly by a recombination event in a human host or a human-created environment and later was transferred to S. aureus.

Christophersen IE, Rienstra M, Roselli C, Yin X, Geelhoed B, Barnard J, Lin H, Arking DE, Smith AV, Albert CM et al. 2017. Large-scale analyses of common and rare variants identify 12 new loci associated with atrial fibrillation. Nat Genet, 49 (6), pp. 946-952. | Citations: 4 (Web of Science Lite) | Show Abstract | Read more

Atrial fibrillation affects more than 33 million people worldwide and increases the risk of stroke, heart failure, and death. Fourteen genetic loci have been associated with atrial fibrillation in European and Asian ancestry groups. To further define the genetic basis of atrial fibrillation, we performed large-scale, trans-ancestry meta-analyses of common and rare variant association studies. The genome-wide association studies (GWAS) included 17,931 individuals with atrial fibrillation and 115,142 referents; the exome-wide association studies (ExWAS) and rare variant association studies (RVAS) involved 22,346 cases and 132,086 referents. We identified 12 new genetic loci that exceeded genome-wide significance, implicating genes involved in cardiac electrical and structural remodeling. Our results nearly double the number of known genetic loci for atrial fibrillation, provide insights into the molecular basis of atrial fibrillation, and may facilitate the identification of new potential targets for drug discovery.

Thomsen SK, Gloyn AL. 2017. Human genetics as a model for target validation: finding new therapies for diabetes. Diabetologia, 60 (6), pp. 960-970. | Citations: 1 (Web of Science Lite) | Show Abstract | Read more

Type 2 diabetes is a global epidemic with major effects on healthcare expenditure and quality of life. Currently available treatments are inadequate for the prevention of comorbidities, yet progress towards new therapies remains slow. A major barrier is the insufficiency of traditional preclinical models for predicting drug efficacy and safety. Human genetics offers a complementary model to assess causal mechanisms for target validation. Genetic perturbations are 'experiments of nature' that provide a uniquely relevant window into the long-term effects of modulating specific targets. Here, we show that genetic discoveries over the past decades have accurately predicted (now known) therapeutic mechanisms for type 2 diabetes. These findings highlight the potential for use of human genetic variation for prospective target validation, and establish a framework for future applications. Studies into rare, monogenic forms of diabetes have also provided proof-of-principle for precision medicine, and the applicability of this paradigm to complex disease is discussed. Finally, we highlight some of the limitations that are relevant to the use of genome-wide association studies (GWAS) in the search for new therapies for diabetes. A key outstanding challenge is the translation of GWAS signals into disease biology and we outline possible solutions for tackling this experimental bottleneck.

Ferreira RC, Rainbow DB, Rubio García A, Pekalski ML, Porter L, Oliveira JJ, Waldron-Lynch F, Wicker LS, Todd JA. 2017. In-depth immunophenotyping data of IL-6R on the human peripheral regulatory T cell (Treg) compartment. Data Brief, 12 pp. 676-691. | Show Abstract | Read more

We provide in this paper a detailed characterization of the human peripheral CD4(+) CD127(low)CD25(+) regulatory T cell (Treg) compartment, with a particular emphasis in defining the population expressing higher levels of the IL-6 receptor (IL-6R). We provide a description of the phenotype of this population by assessing both the surface expression by flow cytometry as well as their transcriptional profile and functional features. In addition, we also present functional data describing the responsiveness of these subsets to IL-6 signalling in vitro and to IL-2 in vivo. The data presented in this paper support the research article "Human IL-6R(hi)TIGIT(-) CD4(+)CD127(low)CD25(+) T cells display potent in vitro suppressive capacity and a distinct Th17 profile" (Ferreira RC et al., 2017; doi: 10.1016/j.clim.2017.03.002) [1].

Scott RA, Scott LJ, Mägi R, Marullo L, Gaulton KJ, Kaakinen M, Pervjakova N, Pers TH, Johnson AD, Eicher JD et al. 2017. An Expanded Genome-Wide Association Study of Type 2 Diabetes in Europeans. Diabetes, pp. db161253-db161253. | Show Abstract | Read more

To characterise type 2 diabetes (T2D) associated variation across the allele frequency spectrum, we conducted a meta-analysis of genome-wide association data from 26,676 T2D cases and 132,532 controls of European ancestry after imputation using the 1000 Genomes multi-ethnic reference panel. Promising association signals were followed-up in additional data sets (of 14,545 or 7,397 T2D cases and 38,994 or 71,604 controls). We identified 13 novel T2D-associated loci (p<5×10(-8)), including variants near the GLP2R, GIP, and HLA-DQA1 genes. Our analysis brought the total number of independent T2D associations to 128 distinct signals at 113 loci. Despite substantially increased sample size and more complete coverage of low-frequency variation, all novel associations were driven by common SNVs. Credible sets of potentially causal variants were generally larger than those based on imputation with earlier reference panels, consistent with resolution of causal signals to common risk haplotypes. Stratification of T2D-associated loci based on T2D-related quantitative trait associations revealed tissue-specific enrichment of regulatory annotations in pancreatic islet enhancers for loci influencing insulin secretion, and in adipocytes, monocytes and hepatocytes for insulin action-associated loci. These findings highlight the predominant role played by common variants of modest effect and the diversity of biological mechanisms influencing T2D pathophysiology.

Vavoulis DV, Taylor JC, Schuh A. 2017. Hierarchical probabilistic models for multiple gene/variant associations based on next-generation sequencing data. Bioinformatics, | Show Abstract | Read more

Motivation: The identification of genetic variants influencing gene expression (known as expression quantitative trait loci or eQTLs ) is important in unravelling the genetic basis of complex traits. Detecting multiple eQTLs simultaneously in a population based on paired DNA-seq and RNA-seq assays employs two competing types of models: models which rely on appropriate transformations of RNA-seq data (and are powered by a mature mathematical theory), or count-based models , which represent digital gene expression explicitly, thus rendering such transformations unnecessary. The latter constitutes an immensely popular methodology, which is however plagued by mathematical intractability. Results: We develop tractable count-based models, which are amenable to efficient estimation through the introduction of latent variables and the appropriate application of recent statistical theory in a sparse Bayesian modelling framework. Furthermore, we examine several transformation methods for RNA-seq read counts and we introduce arcsin , logit and Laplace smoothing as preprocessing steps for transformation-based models. Using natural and carefully simulated data from the 1000 Genomes and gEUVADIS projects, we benchmark both approaches under a variety of scenarios, including the presence of noise and violation of basic model assumptions. We demonstrate that an arcsin transformation of Laplace-smoothed data is at least as good as state-of-the-art models, particularly at small samples. Furthermore, we show that an over-dispersed Poisson model is comparable to the celebrated Negative Binomial, but much easier to estimate. These results provide strong support for transformation-based vs. count-based (particularly Negative-Binomial-based) models for eQTL mapping. Availability: All methods are implemented in the free software eQTLseq : Contact: Supplementary information: Supplementary data are available at Bioinformatics online.

Sapkota Y, Steinthorsdottir V, Morris AP, Fassbender A, Rahmioglu N, De Vivo I, Buring JE, Zhang F, Edwards TL, Jones S et al. 2017. Meta-analysis identifies five novel loci associated with endometriosis highlighting key genes involved in hormone metabolism. Nat Commun, 8 pp. 15539. | Citations: 1 (Scopus) | Show Abstract | Read more

Endometriosis is a heritable hormone-dependent gynecological disorder, associated with severe pelvic pain and reduced fertility; however, its molecular mechanisms remain largely unknown. Here we perform a meta-analysis of 11 genome-wide association case-control data sets, totalling 17,045 endometriosis cases and 191,596 controls. In addition to replicating previously reported loci, we identify five novel loci significantly associated with endometriosis risk (P<5 × 10(-8)), implicating genes involved in sex steroid hormone pathways (FN1, CCDC170, ESR1, SYNE1 and FSHB). Conditional analysis identified five secondary association signals, including two at the ESR1 locus, resulting in 19 independent single nucleotide polymorphisms (SNPs) robustly associated with endometriosis, which together explain up to 5.19% of variance in endometriosis. These results highlight novel variants in or near specific genes with important roles in sex steroid hormone signalling and function, and offer unique opportunities for more targeted functional research efforts.

Gemmell LC, Webster KE, Kirtley S, Vincent K, Zondervan KT, Becker CM. 2017. The management of menopause in women with a history of endometriosis: a systematic review. Hum Reprod Update, 23 (4), pp. 481-500. | Show Abstract | Read more

BACKGROUND: Endometriosis is typically regarded as a premenopausal disease, resolving after natural or iatrogenic menopause due to declining oestrogen levels. Nonetheless, case reports over the years have highlighted the incidence of recurrent postmenopausal endometriosis. It is now clear that both recurrence and malignant transformation of endometriotic foci can occur in the postmenopausal period. Postmenopausal women are commonly treated with hormone replacement therapy (HRT) to treat climacteric symptoms and prevent bone loss; however, HRT may reactivate endometriosis and stimulate malignant transformation in women with a history of endometriosis. Given the uncertain risks of initiating HRT, it is difficult to determine the best menopausal management for this group of women. OBJECTIVE AND RATIONAL: The aim of this study was to systematically review the existing literature on management of menopausal symptoms in women with a history of endometriosis. We also aimed to evaluate the published literature on the risks associated with HRT in these women, and details regarding optimal formulations and timing (i.e. initiation and duration) of HRT. SEARCH METHODS: Four electronic databases (MEDLINE via OVID, Embase via OVID, PsycINFO via OVID and CINAHL via EbscoHost) were searched from database inception until June 2016, using a combination of relevant controlled vocabulary terms and free-text terms related to 'menopause' and 'endometriosis'. Inclusion criteria were: menopausal women with a history of endometriosis and menopausal treatment including HRT or other preparations. Case reports/series, observational studies and clinical trials were included. Narrative review articles, organizational guidelines and conference abstracts were excluded, as were studies that did not report on any form of menopausal management. Articles were assessed for risk of bias and quality using GRADE criteria. OUTCOMES: We present a synthesis of the existing case reports of endometriosis recurrence or malignant transformation in women undergoing treatment for menopausal symptoms. We highlight common presenting symptoms, potential risk factors and outcomes amongst the studies. Sparse high-quality evidence was identified, with few observational studies and only two randomized controlled trials. Given this paucity of data, no definitive conclusions can be drawn concerning risk. WIDER IMPLICATIONS: Due to the lack of high-quality studies, it remains unclear how to advise women with a history of endometriosis regarding the management of menopausal symptoms. The absolute risk of disease recurrence and malignant transformation cannot be quantified, and the impact of HRT use on these outcomes is not known. Multicentre randomized trials or large observational studies are urgently needed to inform clinicians and patients alike.

Parks T, Mirabel MM, Kado J, Auckland K, Nowak J, Rautanen A, Mentzer AJ, Marijon E, Jouven X, Perman ML et al. 2017. Association between a common immunoglobulin heavy chain allele and rheumatic heart disease risk in Oceania. Nat Commun, 8 pp. 14946. | Show Abstract | Read more

The indigenous populations of the South Pacific experience a high burden of rheumatic heart disease (RHD). Here we report a genome-wide association study (GWAS) of RHD susceptibility in 2,852 individuals recruited in eight Oceanian countries. Stratifying by ancestry, we analysed genotyped and imputed variants in Melanesians (607 cases and 1,229 controls) before follow-up of suggestive loci in three further ancestral groups: Polynesians, South Asians and Mixed or other populations (totalling 399 cases and 617 controls). We identify a novel susceptibility signal in the immunoglobulin heavy chain (IGH) locus centring on a haplotype of nonsynonymous variants in the IGHV4-61 gene segment corresponding to the IGHV4-61*02 allele. We show each copy of IGHV4-61*02 is associated with a 1.4-fold increase in the risk of RHD (odds ratio 1.43, 95% confidence intervals 1.27-1.61, P=4.1 × 10(-9)). These findings provide new insight into the role of germline variation in the IGH locus in disease susceptibility.

Cappuccini F, Pollock E, Stribbling S, Hill AVS, Redchenko I. 2017. 5T4 oncofoetal glycoprotein: an old target for a novel prostate cancer immunotherapy. Oncotarget, 8 (29), pp. 47474-47489. | Show Abstract | Read more

The tumour-associated antigen 5T4 is an attractive target for cancer immunotherapy. However to date, reported 5T4-specific cellular immune responses induced by various immunisation platforms have been largely weak or non-existent. In the present study, we have evaluated a heterologous prime boost regime based on the simian adenovirus ChAdOx1 and modified vaccinia virus Ankara (MVA) expressing 5T4 for immunogenicity and tumour protective efficacy in a mouse cancer model. Vaccination-induced immune responses were strong, durable and attributable primarily to CD8+ T cells. By comparison, homologous MVA vaccination regimen did not induce detectable 5T4-specific T cell responses. ChAdOx1-MVA vaccinated mice were completely protected against subsequent B16 melanoma challenge, but in therapeutic settings this regime was only modestly effective in delaying tumour outgrowth. Concomitant delivery of the vaccine with monoclonal antibodies (mAbs) targeting immune checkpoint regulators LAG-3, PD-1 or PD-L1 demonstrated that the combination of vaccine with anti PD-1 mAb could significantly delay tumour growth and increase overall survival of tumour-bearing mice. Our findings support a translation of the combinatorial approach based on the heterologous ChAdOx1-MVA vaccination platform with immune checkpoint blockade into the clinic for the treatment of 5T4-positive tumours such as prostate, renal, colorectal, gastric, ovarian, lung cancer and mesothelioma.

Wright DJ, Day FR, Kerrison ND, Zink F, Cardona A, Sulem P, Thompson DJ, Sigurjonsdottir S, Gudbjartsson DF, Helgason A et al. 2017. Genetic variants associated with mosaic Y chromosome loss highlight cell cycle genes and overlap with cancer susceptibility. Nat Genet, 49 (5), pp. 674-679. | Citations: 2 (Web of Science Lite) | Show Abstract | Read more

The Y chromosome is frequently lost in hematopoietic cells, which represents the most common somatic alteration in men. However, the mechanisms that regulate mosaic loss of chromosome Y (mLOY), and its clinical relevance, are unknown. We used genotype-array-intensity data and sequence reads from 85,542 men to identify 19 genomic regions (P < 5 × 10(-8)) that are associated with mLOY. Cumulatively, these loci also predicted X chromosome loss in women (n = 96,123; P = 4 × 10(-6)). Additional epigenome-wide methylation analyses using whole blood highlighted 36 differentially methylated sites associated with mLOY. The genes identified converge on aspects of cell proliferation and cell cycle regulation, including DNA synthesis (NPAT), DNA damage response (ATM), mitosis (PMF1, CENPN and MAD1L1) and apoptosis (TP53). We highlight the shared genetic architecture between mLOY and cancer susceptibility, in addition to inferring a causal effect of smoking on mLOY. Collectively, our results demonstrate that genotype-array-intensity data enables a measure of cell cycle efficiency at population scale and identifies genes implicated in aneuploidy, genome instability and cancer susceptibility.

Purshouse K, Schuh A, Fairfax BP, Knight S, Antoniou P, Dreau H, Popitsch N, Gatter K, Roberts I, Browning L et al. 2017. Whole-genome sequencing identifies homozygous BRCA2 deletion guiding treatment in dedifferentiated prostate cancer. Cold Spring Harb Mol Case Stud, 3 (3), pp. a001362. | Show Abstract | Read more

Whole-genome sequencing (WGS) has transformed the understanding of the genetic drivers of cancer and is increasingly being used in cancer medicine to identify personalized therapies. Here we describe a case in which the application of WGS identified a tumoral BRCA2 deletion in a patient with aggressive dedifferentiated prostate cancer that was repeat-biopsied after disease progression. This would not have been detected by standard BRCA testing, and it led to additional treatment with a maintenance poly ADP ribose polymerase (PARP) inhibitor following platinum-based chemotherapy. This case demonstrates that repeat biopsy upon disease progression and application of WGS to tumor samples has meaningful clinical utility and the potential to transform outcomes in patients with cancer.

Windholz J, Kovacs P, Schlicke M, Franke C, Mahajan A, Morris AP, Lemke JR, Klammt J, Kiess W, Schöneberg T et al. 2017. Copy number variations in "classical" obesity candidate genes are not frequently associated with severe early-onset obesity in children. J Pediatr Endocrinol Metab, 30 (5), pp. 507-515. | Show Abstract | Read more

BACKGROUND: Obesity is genetically heterogeneous and highly heritable, although polymorphisms explain the phenotype in only a small proportion of obese children. We investigated the presence of copy number variations (CNVs) in "classical" genes known to be associated with (monogenic) early-onset obesity in children. METHODS: In 194 obese Caucasian children selected for early-onset and severe obesity from our obesity cohort we screened for deletions and/or duplications by multiplex ligation-dependent probe amplification reaction (MLPA). As we found one MLPA probe to interfere with a polymorphism in SIM1 we investigated its association with obesity and other phenotypic traits in our extended cohort of 2305 children. RESULTS: In the selected subset of most severely obese children, we did not find CNV with MLPA in POMC, LEP, LEPR, MC4R, MC3R or MC2R genes. However, one SIM1 probe located at exon 9 gave signals suggestive for SIM1 insufficiency in 52 patients. Polymerase chain reaction (PCR) analysis identified this as a false positive result due to interference with single nucleotide polymorphism (SNP) rs3734354/rs3734355. We, therefore, investigated for associations of this polymorphism with obesity and metabolic traits in our extended cohort. We found rs3734354/rs3734355 to be associated with body mass index-standard deviation score (BMI-SDS) (p = 0.003), but not with parameters of insulin metabolism, blood pressure or food intake. CONCLUSIONS: In our modest sample of severely obese children, we were unable to find CNVs in well-established monogenic obesity genes. Nevertheless, we found an association of rs3734354 in SIM1 with obesity of early-onset type in children, although not with obesity-related traits.

Wild PS, Felix JF, Schillert A, Teumer A, Chen MH, Leening MJG, Völker U, Großmann V, Brody JA, Irvin MR et al. 2017. Large-scale genome-wide analysis identifies genetic variants associated with cardiac structure and function. J Clin Invest, 127 (5), pp. 1798-1812. | Show Abstract | Read more

BACKGROUND: Understanding the genetic architecture of cardiac structure and function may help to prevent and treat heart disease. This investigation sought to identify common genetic variations associated with inter-individual variability in cardiac structure and function. METHODS: A GWAS meta-analysis of echocardiographic traits was performed, including 46,533 individuals from 30 studies (EchoGen consortium). The analysis included 16 traits of left ventricular (LV) structure, and systolic and diastolic function. RESULTS: The discovery analysis included 21 cohorts for structural and systolic function traits (n = 32,212) and 17 cohorts for diastolic function traits (n = 21,852). Replication was performed in 5 cohorts (n = 14,321) and 6 cohorts (n = 16,308), respectively. Besides 5 previously reported loci, the combined meta-analysis identified 10 additional genome-wide significant SNPs: rs12541595 near MTSS1 and rs10774625 in ATXN2 for LV end-diastolic internal dimension; rs806322 near KCNRG, rs4765663 in CACNA1C, rs6702619 near PALMD, rs7127129 in TMEM16A, rs11207426 near FGGY, rs17608766 in GOSR2, and rs17696696 in CFDP1 for aortic root diameter; and rs12440869 in IQCH for Doppler transmitral A-wave peak velocity. Findings were in part validated in other cohorts and in GWAS of related disease traits. The genetic loci showed associations with putative signaling pathways, and with gene expression in whole blood, monocytes, and myocardial tissue. CONCLUSION: The additional genetic loci identified in this large meta-analysis of cardiac structure and function provide insights into the underlying genetic architecture of cardiac structure and warrant follow-up in future functional studies. FUNDING: For detailed information per study, see Acknowledgments.

Justice AE, Winkler TW, Feitosa MF, Graff M, Fisher VA, Young K, Barata L, Deng X, Czajkowski J, Hadley D et al. 2017. Genome-wide meta-analysis of 241,258 adults accounting for smoking behaviour identifies novel loci for obesity traits. Nat Commun, 8 pp. 14977. | Citations: 2 (Scopus) | Show Abstract | Read more

Few genome-wide association studies (GWAS) account for environmental exposures, like smoking, potentially impacting the overall trait variance when investigating the genetic contribution to obesity-related traits. Here, we use GWAS data from 51,080 current smokers and 190,178 nonsmokers (87% European descent) to identify loci influencing BMI and central adiposity, measured as waist circumference and waist-to-hip ratio both adjusted for BMI. We identify 23 novel genetic loci, and 9 loci with convincing evidence of gene-smoking interaction (GxSMK) on obesity-related traits. We show consistent direction of effect for all identified loci and significance for 18 novel and for 5 interaction loci in an independent study sample. These loci highlight novel biological functions, including response to oxidative stress, addictive behaviour, and regulatory functions emphasizing the importance of accounting for environment in genetic analyses. Our results suggest that tobacco smoking may alter the genetic susceptibility to overall adiposity and body fat distribution.

Collins KA, Snaith R, Cottingham MG, Gilbert SC, Hill AVS. 2017. Enhancing protective immunity to malaria with a highly immunogenic virus-like particle vaccine. Sci Rep, 7 pp. 46621. | Show Abstract | Read more

The leading malaria vaccine in development is the circumsporozoite protein (CSP)-based particle vaccine, RTS,S, which targets the pre-erythrocytic stage of Plasmodium falciparum infection. It induces modest levels of protective efficacy, thought to be mediated primarily by CSP-specific antibodies. We aimed to enhance vaccine efficacy by generating a more immunogenic CSP-based particle vaccine and therefore developed a next-generation RTS,S-like vaccine, called R21. The major improvement is that in contrast to RTS,S, R21 particles are formed from a single CSP-hepatitis B surface antigen (HBsAg) fusion protein, and this leads to a vaccine composed of a much higher proportion of CSP than in RTS,S. We demonstrate that in BALB/c mice R21 is immunogenic at very low doses and when administered with the adjuvants Abisco-100 and Matrix-M it elicits sterile protection against transgenic sporozoite challenge. Concurrent induction of potent cellular and humoral immune responses was also achieved by combining R21 with TRAP-based viral vectors and protective efficacy was significantly enhanced. In addition, in contrast to RTS,S, only a minimal antibody response to the HBsAg carrier was induced. These studies identify an anti-sporozoite vaccine component that may improve upon the current leading malaria vaccine RTS,S. R21 is now under evaluation in Phase 1/2a clinical trials.

Salman AM, Montoya-Díaz E, West H, Lall A, Atcheson E, Lopez-Camacho C, Ramesar J, Bauza K, Collins KA, Brod F et al. 2017. Rational development of a protective P. vivax vaccine evaluated with transgenic rodent parasite challenge models. Sci Rep, 7 pp. 46482. | Citations: 1 (Scopus) | Show Abstract | Read more

Development of a protective and broadly-acting vaccine against the most widely distributed human malaria parasite, Plasmodium vivax, will be a major step towards malaria elimination. However, a P. vivax vaccine has remained elusive by the scarcity of pre-clinical models to test protective efficacy and support further clinical trials. In this study, we report the development of a highly protective CSP-based P. vivax vaccine, a virus-like particle (VLP) known as Rv21, able to provide 100% sterile protection against a stringent sporozoite challenge in rodent models to malaria, where IgG2a antibodies were associated with protection in absence of detectable PvCSP-specific T cell responses. Additionally, we generated two novel transgenic rodent P. berghei parasite lines, where the P. berghei csp gene coding sequence has been replaced with either full-length P. vivax VK210 or the allelic VK247 csp that additionally express GFP-Luciferase. Efficacy of Rv21 surpassed viral-vectored vaccination using ChAd63 and MVA. We show for the first time that a chimeric VK210/247 antigen can elicit high level cross-protection against parasites expressing either CSP allele, which provide accessible and affordable models suitable to support the development of P. vivax vaccines candidates. Rv21 is progressing to GMP production and has entered a path towards clinical evaluation.

Vicente JL, Clarkson CS, Caputo B, Gomes B, Pombi M, Sousa CA, Antao T, Dinis J, Bottà G, Mancini E et al. 2017. Massive introgression drives species radiation at the range limit of Anopheles gambiae. Sci Rep, 7 pp. 46451. | Show Abstract | Read more

Impacts of introgressive hybridisation may range from genomic erosion and species collapse to rapid adaptation and speciation but opportunities to study these dynamics are rare. We investigated the extent, causes and consequences of a hybrid zone between Anopheles coluzzii and Anopheles gambiae in Guinea-Bissau, where high hybridisation rates appear to be stable at least since the 1990s. Anopheles gambiae was genetically partitioned into inland and coastal subpopulations, separated by a central region dominated by A. coluzzii. Surprisingly, whole genome sequencing revealed that the coastal region harbours a hybrid form characterised by an A. gambiae-like sex chromosome and massive introgression of A. coluzzii autosomal alleles. Local selection on chromosomal inversions may play a role in this process, suggesting potential for spatiotemporal stability of the coastal hybrid form and providing resilience against introgression of medically-important loci and traits, found to be more prevalent in inland A. gambiae.

Ansari MA, Pedergnana V, L C Ip C, Magri A, Von Delft A, Bonsall D, Chaturvedi N, Bartha I, Smith D, Nicholson G et al. 2017. Genome-to-genome analysis highlights the effect of the human innate and adaptive immune systems on the hepatitis C virus. Nat Genet, 49 (5), pp. 666-673. | Citations: 2 (Web of Science Lite) | Show Abstract | Read more

Outcomes of hepatitis C virus (HCV) infection and treatment depend on viral and host genetic factors. Here we use human genome-wide genotyping arrays and new whole-genome HCV viral sequencing technologies to perform a systematic genome-to-genome study of 542 individuals who were chronically infected with HCV, predominantly genotype 3. We show that both alleles of genes encoding human leukocyte antigen molecules and genes encoding components of the interferon lambda innate immune system drive viral polymorphism. Additionally, we show that IFNL4 genotypes determine HCV viral load through a mechanism dependent on a specific amino acid residue in the HCV NS5A protein. These findings highlight the interplay between the innate immune system and the viral genome in HCV control.

Roberts AR, Vecellio M, Cortes A, Knight JC, Cohen CJ, Wordsworth BP. 2017. Investigation of a possible extended risk haplotype in the IL23R region associated with ankylosing spondylitis. Genes Immun, 18 (2), pp. 105-108. | Show Abstract | Read more

The IL23R region on chromosome 1 exhibits complex associations with ankylosing spondylitis (AS). We used publicly available epigenomic information and historical genetic association data to identify a putative regulatory element (PRE) in the intergenic region between IL23R and IL12RB2, which includes two single-nucleotide polymorphisms (SNPs) independently associated with AS-rs924080 (P=2 × 10(-3)) and rs11578380 (P=2 × 10(-4)). In luciferase reporter assays, this PRE showed silencer activity (P<0.001). Haplotype and conditional analysis of 4230 historical AS cases and 9700 controls revealed a possible AS-associated extended haplotype, including the PRE and risk variants at three SNPs (rs11209026, rs11209032 and rs924080), but excluding the rs11578380 risk variant. However, the rs924080 association was absent after conditioning on the primary association with rs11209032, which, in contrast, was robust to conditioning on all other AS-associated SNPs in this region (P<2 × 10(-8)). The role of this putative silencer on some IL23R extended haplotypes therefore remains unclear.

Ford T, Wenden C, Mbekeani A, Dally L, Cox JH, Morin M, Winstone N, Hill AV, Gilmour J, Ewer KJ. 2017. Cryopreservation-related loss of antigen-specific IFNγ producing CD4(+) T-cells can skew immunogenicity data in vaccine trials: Lessons from a malaria vaccine trial substudy. Vaccine, 35 (15), pp. 1898-1906. | Show Abstract | Read more

Ex vivo functional immunoassays such as ELISpot and intracellular cytokine staining (ICS) by flow cytometry are crucial tools in vaccine development both in the identification of novel immunogenic targets and in the immunological assessment of samples from clinical trials. Cryopreservation and subsequent thawing of PBMCs via validated processes has become a mainstay of clinical trials due to processing restrictions inherent in the disparate location and capacity of trial centres, and also in the need to standardize biological assays at central testing facilities. Logistical and financial requirement to batch process samples from multiple study timepoints are also key. We used ELISpot and ICS assays to assess antigen-specific immunogenicity in blood samples taken from subjects enrolled in a phase II malaria heterologous prime-boost vaccine trial and showed that the freeze thaw process can result in a 3-5-fold reduction of malaria antigen-specific IFNγ-producing CD3(+)CD4(+) effector populations from PBMC samples taken post vaccination. We have also demonstrated that peptide responsive CD8(+) T cells are relatively unaffected, as well as CD4(+) T cell populations that do not produce IFNγ. These findings contribute to a growing body of data that could be consolidated and synthesised as guidelines for clinical trials with the aim of increasing the efficiency of vaccine development pipelines.

Folkersen L, Fauman E, Sabater-Lleal M, Strawbridge RJ, Frånberg M, Sennblad B, Baldassarre D, Veglia F, Humphries SE, Rauramaa R et al. 2017. Mapping of 79 loci for 83 plasma protein biomarkers in cardiovascular disease. PLoS Genet, 13 (4), pp. e1006706. | Show Abstract | Read more

Recent advances in highly multiplexed immunoassays have allowed systematic large-scale measurement of hundreds of plasma proteins in large cohort studies. In combination with genotyping, such studies offer the prospect to 1) identify mechanisms involved with regulation of protein expression in plasma, and 2) determine whether the plasma proteins are likely to be causally implicated in disease. We report here the results of genome-wide association (GWA) studies of 83 proteins considered relevant to cardiovascular disease (CVD), measured in 3,394 individuals with multiple CVD risk factors. We identified 79 genome-wide significant (p<5e-8) association signals, 55 of which replicated at P<0.0007 in separate validation studies (n = 2,639 individuals). Using automated text mining, manual curation, and network-based methods incorporating information on expression quantitative trait loci (eQTL), we propose plausible causal mechanisms for 25 trans-acting loci, including a potential post-translational regulation of stem cell factor by matrix metalloproteinase 9 and receptor-ligand pairs such as RANK-RANK ligand. Using public GWA study data, we further evaluate all 79 loci for their causal effect on coronary artery disease, and highlight several potentially causal associations. Overall, a majority of the plasma proteins studied showed evidence of regulation at the genetic level. Our results enable future studies of the causal architecture of human disease, which in turn should aid discovery of new drug targets.

Rolo J, Worning P, Boye Nielsen J, Sobral R, Bowden R, Bouchami O, Damborg P, Guardabassi L, Perreten V, Westh H et al. 2017. Evidence for the evolutionary steps leading to mecA-mediated β-lactam resistance in staphylococci. PLoS Genet, 13 (4), pp. e1006674. | Citations: 2 (Web of Science Lite) | Show Abstract | Read more

The epidemiologically most important mechanism of antibiotic resistance in Staphylococcus aureus is associated with mecA-an acquired gene encoding an extra penicillin-binding protein (PBP2a) with low affinity to virtually all β-lactams. The introduction of mecA into the S. aureus chromosome has led to the emergence of methicillin-resistant S. aureus (MRSA) pandemics, responsible for high rates of mortality worldwide. Nonetheless, little is known regarding the origin and evolution of mecA. Different mecA homologues have been identified in species belonging to the Staphylococcus sciuri group representing the most primitive staphylococci. In this study we aimed to identify evolutionary steps linking these mecA precursors to the β-lactam resistance gene mecA and the resistance phenotype. We sequenced genomes of 106 S. sciuri, S. vitulinus and S. fleurettii strains and determined their oxacillin susceptibility profiles. Single-nucleotide polymorphism (SNP) analysis of the core genome was performed to assess the genetic relatedness of the isolates. Phylogenetic analysis of the mecA gene homologues and promoters was achieved through nucleotide/amino acid sequence alignments and mutation rates were estimated using a Bayesian analysis. Furthermore, the predicted structure of mecA homologue-encoded PBPs of oxacillin-susceptible and -resistant strains were compared. We showed for the first time that oxacillin resistance in the S. sciuri group has emerged multiple times and by a variety of different mechanisms. Development of resistance occurred through several steps including structural diversification of the non-binding domain of native PBPs; changes in the promoters of mecA homologues; acquisition of SCCmec and adaptation of the bacterial genetic background. Moreover, our results suggest that it was exposure to β-lactams in human-created environments that has driven evolution of native PBPs towards a resistance determinant. The evolution of β-lactam resistance in staphylococci highlights the numerous resources available to bacteria to adapt to the selective pressure of antibiotics.

Ju YS, Martincorena I, Gerstung M, Petljak M, Alexandrov LB, Rahbari R, Wedge DC, Davies HR, Ramakrishna M, Fullam A et al. 2017. Somatic mutations reveal asymmetric cellular dynamics in the early human embryo. Nature, 543 (7647), pp. 714-718. | Citations: 4 (Scopus) | Show Abstract | Read more

Somatic cells acquire mutations throughout the course of an individual's life. Mutations occurring early in embryogenesis are often present in a substantial proportion of, but not all, cells in postnatal humans and thus have particular characteristics and effects. Depending on their location in the genome and the proportion of cells they are present in, these mosaic mutations can cause a wide range of genetic disease syndromes and predispose carriers to cancer. They have a high chance of being transmitted to offspring as de novo germline mutations and, in principle, can provide insights into early human embryonic cell lineages and their contributions to adult tissues. Although it is known that gross chromosomal abnormalities are remarkably common in early human embryos, our understanding of early embryonic somatic mutations is very limited. Here we use whole-genome sequences of normal blood from 241 adults to identify 163 early embryonic mutations. We estimate that approximately three base substitution mutations occur per cell per cell-doubling event in early human embryogenesis and these are mainly attributable to two known mutational signatures. We used the mutations to reconstruct developmental lineages of adult cells and demonstrate that the two daughter cells of many early embryonic cell-doubling events contribute asymmetrically to adult blood at an approximately 2:1 ratio. This study therefore provides insights into the mutation rates, mutational processes and developmental outcomes of cell dynamics that operate during early human embryogenesis.

Reilly SN, Liu X, Carnicer R, Recalde A, Muszkiewicz A, Jayaram R, Carena MC, Wijesurendra R, Stefanini M, Surdo NC et al. 2016. Up-regulation of miR-31 in human atrial fibrillation begets the arrhythmia by depleting dystrophin and neuronal nitric oxide synthase. Sci Transl Med, 8 (340), pp. 340ra74. | Citations: 4 (Web of Science Lite) | Show Abstract | Read more

Atrial fibrillation (AF) is a growing public health burden, and its treatment remains a challenge. AF leads to electrical remodeling of the atria, which in turn promotes AF maintenance and resistance to treatment. Although remodeling has long been a therapeutic target in AF, its causes remain poorly understood. We show that atrial-specific up-regulation of microRNA-31 (miR-31) in goat and human AF depletes neuronal nitric oxide synthase (nNOS) by accelerating mRNA decay and alters nNOS subcellular localization by repressing dystrophin translation. By shortening action potential duration and abolishing rate-dependent adaptation of the action potential duration, miR-31 overexpression and/or disruption of nNOS signaling recapitulates features of AF-induced remodeling and significantly increases AF inducibility in mice in vivo. By contrast, silencing miR-31 in atrial myocytes from patients with AF restores dystrophin and nNOS and normalizes action potential duration and its rate dependency. These findings identify atrial-specific up-regulation of miR-31 in human AF as a key mechanism causing atrial dystrophin and nNOS depletion, which in turn contributes to the atrial phenotype begetting this arrhythmia. miR-31 may therefore represent a potential therapeutic target in AF.

Manning A, Highland HM, Gasser J, Sim X, Tukiainen T, Fontanillas P, Grarup N, Rivas MA, Mahajan A, Locke AE et al. 2017. A Low-Frequency Inactivating AKT2 Variant Enriched in the Finnish Population Is Associated With Fasting Insulin Levels and Type 2 Diabetes Risk. Diabetes, 66 (7), pp. 2019-2032. | Show Abstract | Read more

To identify novel coding association signals and facilitate characterization of mechanisms influencing glycemic traits and type 2 diabetes risk, we analyzed 109,215 variants derived from exome array genotyping together with an additional 390,225 variants from exome sequence in up to 39,339 normoglycemic individuals from five ancestry groups. We identified a novel association between the coding variant (p.Pro50Thr) in AKT2 and fasting plasma insulin (FI), a gene in which rare fully penetrant mutations are causal for monogenic glycemic disorders. The low-frequency allele is associated with a 12% increase in FI levels. This variant is present at 1.1% frequency in Finns but virtually absent in individuals from other ancestries. Carriers of the FI-increasing allele had increased 2-h insulin values, decreased insulin sensitivity, and increased risk of type 2 diabetes (odds ratio 1.05). In cellular studies, the AKT2-Thr50 protein exhibited a partial loss of function. We extend the allelic spectrum for coding variants in AKT2 associated with disorders of glucose homeostasis and demonstrate bidirectional effects of variants within the pleckstrin homology domain of AKT2.

Pagnamenta AT, Murakami Y, Taylor JM, Anzilotti C, Howard MF, Miller V, Johnson DS, Tadros S, Mansour S, Temple IK et al. 2017. Analysis of exome data for 4293 trios suggests GPI-anchor biogenesis defects are a rare cause of developmental disorders. Eur J Hum Genet, 25 (6), pp. 669-679. | Show Abstract | Read more

Over 150 different proteins attach to the plasma membrane using glycosylphosphatidylinositol (GPI) anchors. Mutations in 18 genes that encode components of GPI-anchor biogenesis result in a phenotypic spectrum that includes learning disability, epilepsy, microcephaly, congenital malformations and mild dysmorphic features. To determine the incidence of GPI-anchor defects, we analysed the exome data from 4293 parent-child trios recruited to the Deciphering Developmental Disorders (DDD) study. All probands recruited had a neurodevelopmental disorder. We searched for variants in 31 genes linked to GPI-anchor biogenesis and detected rare biallelic variants in PGAP3, PIGN, PIGT (n=2), PIGO and PIGL, providing a likely diagnosis for six families. In five families, the variants were in a compound heterozygous configuration while in a consanguineous Afghani kindred, a homozygous c.709G>C; p.(E237Q) variant in PIGT was identified within 10-12 Mb of autozygosity. Validation and segregation analysis was performed using Sanger sequencing. Across the six families, five siblings were available for testing and in all cases variants co-segregated consistent with them being causative. In four families, abnormal alkaline phosphatase results were observed in the direction expected. FACS analysis of knockout HEK293 cells that had been transfected with wild-type or mutant cDNA constructs demonstrated that the variants in PIGN, PIGT and PIGO all led to reduced activity. Splicing assays, performed using leucocyte RNA, showed that a c.336-2A>G variant in PIGL resulted in exon skipping and p.D113fs*2. Our results strengthen recently reported disease associations, suggest that defective GPI-anchor biogenesis may explain ~0.15% of individuals with developmental disorders and highlight the benefits of data sharing.

Ormondroyd E, Mackley MP, Blair E, Craft J, Knight JC, Taylor J, Taylor JC, Wilkie AO, Watkins H. 2017. Insights from early experience of a Rare Disease Genomic Medicine Multidisciplinary Team: a qualitative study. Eur J Hum Genet, 25 (6), pp. 680-686. | Show Abstract | Read more

Whole-exome/whole-genome sequencing (WES/WGS) has the potential to enhance genetic diagnosis of rare disease, and is increasingly becoming part of routine clinical care in mainstream medicine. Effective translation will require ongoing efforts in a number of areas including: selection of appropriate patients, provision of effective consent, pre- and post-test genetic counselling, improving variant interpretation algorithms and practices, and management of secondary findings including those found incidentally and those actively sought. Allied to this is the need for an effective education programme for all members of clinical teams involved in care of patients with rare disease, as well as to maintain public confidence in the use of these technologies. We established a Genomic Medicine Multidisciplinary Team (GM-MDT) in 2014 to build on the experiences of earlier successful research-based WES/WGS studies, to address these needs and to review results including pertinent and secondary findings. Here we report on a qualitative study of decision-making in the GM-MDT combined with analysis of semi-structured interviews with GM-MDT members. Study findings show that members appreciate the clinical and scientific diversity of the GM-MDT and value it for education and oversight. To date, discussions have focussed on case selection including the extent and interpretation of clinical and family history information required to establish likely monogenic aetiology and inheritance model. Achieving a balance between effective use of WES/WGS - prioritising cases in a diverse and highly complex patient population where WES/WGS will be tractable - and meeting the recruitment targets of a large project is considered challenging.

Monahan KJ, Alsina D, Bach S, Buchanan J, Burn J, Clark S, Dawson P, De Souza B, Din FV, Dolwani S et al. 2017. Urgent improvements needed to diagnose and manage Lynch syndrome. BMJ, 356 pp. j1388. | Citations: 1 (Web of Science Lite) | Read more

Hoehn KB, Lunter G, Pybus OG. 2017. A Phylogenetic Codon Substitution Model for Antibody Lineages. Genetics, 206 (1), pp. 417-427. | Show Abstract | Read more

Phylogenetic methods have shown promise in understanding the development of broadly neutralizing antibody lineages (bNAbs). However, the mutational process that generates these lineages, somatic hypermutation, is biased by hotspot motifs which violates important assumptions in most phylogenetic substitution models. Here, we develop a modified GY94-type substitution model that partially accounts for this context dependency while preserving independence of sites during calculation. This model shows a substantially better fit to three well-characterized bNAb lineages than the standard GY94 model. We also demonstrate how our model can be used to test hypotheses concerning the roles of different hotspot and coldspot motifs in the evolution of B-cell lineages. Further, we explore the consequences of the idea that the number of hotspot motifs, and perhaps the mutation rate in general, is expected to decay over time in individual bNAb lineages.

Irshad S, Bansal M, Guarnieri P, Davis H, Al Haj Zen A, Baran B, Pinna CMA, Rahman H, Biswas S, Bardella C et al. 2017. Bone morphogenetic protein and Notch signalling crosstalk in poor-prognosis, mesenchymal-subtype colorectal cancer. J Pathol, 242 (2), pp. 178-192. | Show Abstract | Read more

The functional role of bone morphogenetic protein (BMP) signalling in colorectal cancer (CRC) is poorly defined, with contradictory results in cancer cell line models reflecting the inherent difficulties of assessing a signalling pathway that is context-dependent and subject to genetic constraints. By assessing the transcriptional response of a diploid human colonic epithelial cell line to BMP ligand stimulation, we generated a prognostic BMP signalling signature, which was applied to multiple CRC datasets to investigate BMP heterogeneity across CRC molecular subtypes. We linked BMP and Notch signalling pathway activity and function in human colonic epithelial cells, and normal and neoplastic tissue. BMP induced Notch through a γ-secretase-independent interaction, regulated by the SMAD proteins. In homeostasis, BMP/Notch co-localization was restricted to cells at the top of the intestinal crypt, with more widespread interaction in some human CRC samples. BMP signalling was downregulated in the majority of CRCs, but was conserved specifically in mesenchymal-subtype tumours, where it interacts with Notch to induce an epithelial-mesenchymal transition (EMT) phenotype. In intestinal homeostasis, BMP-Notch pathway crosstalk is restricted to differentiating cells through stringent pathway segregation. Conserved BMP activity and loss of signalling stringency in mesenchymal-subtype tumours promotes a synergistic BMP-Notch interaction, and this correlates with poor patient prognosis. BMP signalling heterogeneity across CRC subtypes and cell lines can account for previous experimental contradictions. Crosstalk between the BMP and Notch pathways will render mesenchymal-subtype CRC insensitive to γ-secretase inhibition unless BMP activation is concomitantly addressed. © 2017 The Authors. Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

Hu Y, Tanaka T, Zhu J, Guan W, Wu JHY, Psaty BM, McKnight B, King IB, Sun Q, Richard M et al. 2017. Discovery and fine-mapping of loci associated with MUFAs through trans-ethnic meta-analysis in Chinese and European populations. J Lipid Res, 58 (5), pp. 974-981. | Show Abstract | Read more

MUFAs are unsaturated FAs with one double bond and are derived from endogenous synthesis and dietary intake. Accumulating evidence has suggested that plasma and erythrocyte MUFA levels are associated with cardiometabolic disorders, including CVD, T2D, and metabolic syndrome (MS). Previous genome-wide association studies (GWASs) have identified seven loci for plasma and erythrocyte palmitoleic and oleic acid levels in populations of European origin. To identify additional MUFA-associated loci and the potential functional variant at each locus, we performed ethnic-specific GWAS meta-analyses and trans-ethnic meta-analyses in more than 15,000 participants of Chinese and European ancestry. We identified novel genome-wide significant associations for vaccenic acid at FADS1/2 and PKD2L1 [log10(Bayes factor) ≥ 8.07] and for gondoic acid at FADS1/2 and GCKR [log10(Bayes factor) ≥ 6.22], and also observed improved fine-mapping resolutions at FADS1/2 and GCKR loci. The greatest improvement was observed at GCKR, where the number of variants in the 99% credible set was reduced from 16 (covering 94.8 kb) to 5 (covering 19.6 kb, including a missense variant rs1260326) after trans-ethnic meta-analysis. We also confirmed the previously reported associations of PKD2L1, FADS1/2, GCKR, and HIF1AN with palmitoleic acid and of FADS1/2 and LPCAT3 with oleic acid in the Chinese-specific GWAS and the trans-ethnic meta-analyses. Pathway-based analyses suggested that the identified loci were in unsaturated FA metabolism and signaling pathways. Our findings provide novel insight into the genetic basis relevant to MUFA metabolism and biology.

Campbell KR, Yau C. 2017. Probabilistic modeling of bifurcations in single-cell gene expression data using a Bayesian mixture of factor analyzers. Wellcome Open Res, 2 pp. 19. | Show Abstract | Read more

Modeling bifurcations in single-cell transcriptomics data has become an increasingly popular field of research. Several methods have been proposed to infer bifurcation structure from such data, but all rely on heuristic non-probabilistic inference. Here we propose the first generative, fully probabilistic model for such inference based on a Bayesian hierarchical mixture of factor analyzers. Our model exhibits competitive performance on large datasets despite implementing full Markov-Chain Monte Carlo sampling, and its unique hierarchical prior structure enables automatic determination of genes driving the bifurcation process. We additionally propose an Empirical-Bayes like extension that deals with the high levels of zero-inflation in single-cell RNA-seq data and quantify when such models are useful. We apply or model to both real and simulated single-cell gene expression data and compare the results to existing pseudotime methods. Finally, we discuss both the merits and weaknesses of such a unified, probabilistic approach in the context practical bioinformatics analyses.

Shimelis H, Mesman RLS, Von Nicolai C, Ehlen A, Guidugli L, Martin C, Calléja FMGR, Meeks H, Hallberg E, Hinton J et al. 2017. BRCA2 Hypomorphic Missense Variants Confer Moderate Risks of Breast Cancer. Cancer Res, 77 (11), pp. 2789-2799. | Citations: 1 (Web of Science Lite) | Show Abstract | Read more

Breast cancer risks conferred by many germline missense variants in the BRCA1 and BRCA2 genes, often referred to as variants of uncertain significance (VUS), have not been established. In this study, associations between 19 BRCA1 and 33 BRCA2 missense substitution variants and breast cancer risk were investigated through a breast cancer case-control study using genotyping data from 38 studies of predominantly European ancestry (41,890 cases and 41,607 controls) and nine studies of Asian ancestry (6,269 cases and 6,624 controls). The BRCA2 c.9104A>C, p.Tyr3035Ser (OR = 2.52; P = 0.04), and BRCA1 c.5096G>A, p.Arg1699Gln (OR = 4.29; P = 0.009) variant were associated with moderately increased risks of breast cancer among Europeans, whereas BRCA2 c.7522G>A, p.Gly2508Ser (OR = 2.68; P = 0.004), and c.8187G>T, p.Lys2729Asn (OR = 1.4; P = 0.004) were associated with moderate and low risks of breast cancer among Asians. Functional characterization of the BRCA2 variants using four quantitative assays showed reduced BRCA2 activity for p.Tyr3035Ser compared with wild-type. Overall, our results show how BRCA2 missense variants that influence protein function can confer clinically relevant, moderately increased risks of breast cancer, with potential implications for risk management guidelines in women with these specific variants. Cancer Res; 77(11); 2789-99. ©2017 AACR.

Ferreira RC, Rainbow DB, Rubio García A, Pekalski ML, Porter L, Oliveira JJ, Waldron-Lynch F, Wicker LS, Todd JA. 2017. Human IL-6R(hi)TIGIT(-) CD4(+)CD127(low)CD25(+) T cells display potent in vitro suppressive capacity and a distinct Th17 profile. Clin Immunol, 179 pp. 25-39. | Citations: 1 (Scopus) | Show Abstract | Read more

To date many clinical studies aim to increase the number and/or fitness of CD4(+)CD127(low)CD25(+) regulatory T cells (Tregs) in vivo to harness their regulatory potential in the context of treating autoimmune disease. Here, we sought to define the phenotype and function of Tregs expressing the highest levels of IL-6 receptor (IL-6R). We have identified a population of CD4(+)CD127(low)CD25(+) TIGIT(-) T cells distinguished by their elevated IL-6R expression that lacked expression of HELIOS, showed higher CTLA-4 expression, and displayed increased suppressive capacity compared to IL-6R(hi)TIGIT(+) Tregs. IL-6R(hi)TIGIT(-) CD127(low)CD25(+) T cells contained a majority of cells demethylated at FOXP3 and displayed a Th17 transcriptional signature, including RORC (RORγt) and the capacity of producing both pro- and anti-inflammatory cytokines, such as IL-17, IL-22 and IL-10. We propose that in vivo, in the presence of IL-6-associated inflammation, the suppressive function of CD4(+)CD127(low)CD25(+) FOXP3(+)IL-6R(hi)TIGIT(-) T cells is temporarily disarmed allowing further activation of the effector functions and potential pathogenic tissue damage.

Votintseva AA, Bradley P, Pankhurst L, Del Ojo Elias C, Loose M, Nilgiriwala K, Chatterjee A, Smith EG, Sanderson N, Walker TM et al. 2017. Same-Day Diagnostic and Surveillance Data for Tuberculosis via Whole-Genome Sequencing of Direct Respiratory Samples. J Clin Microbiol, 55 (5), pp. 1285-1298. | Citations: 7 (Scopus) | Show Abstract | Read more

Routine full characterization of Mycobacterium tuberculosis is culture based, taking many weeks. Whole-genome sequencing (WGS) can generate antibiotic susceptibility profiles to inform treatment, augmented with strain information for global surveillance; such data could be transformative if provided at or near the point of care. We demonstrate a low-cost method of DNA extraction directly from patient samples for M. tuberculosis WGS. We initially evaluated the method by using the Illumina MiSeq sequencer (40 smear-positive respiratory samples obtained after routine clinical testing and 27 matched liquid cultures). M. tuberculosis was identified in all 39 samples from which DNA was successfully extracted. Sufficient data for antibiotic susceptibility prediction were obtained from 24 (62%) samples; all results were concordant with reference laboratory phenotypes. Phylogenetic placement was concordant between direct and cultured samples. With Illumina MiSeq/MiniSeq, the workflow from patient sample to results can be completed in 44/16 h at a reagent cost of £96/£198 per sample. We then employed a nonspecific PCR-based library preparation method for sequencing on an Oxford Nanopore Technologies MinION sequencer. We applied this to cultured Mycobacterium bovis strain BCG DNA and to combined culture-negative sputum DNA and BCG DNA. For flow cell version R9.4, the estimated turnaround time from patient to identification of BCG, detection of pyrazinamide resistance, and phylogenetic placement was 7.5 h, with full susceptibility results 5 h later. Antibiotic susceptibility predictions were fully concordant. A critical advantage of MinION is the ability to continue sequencing until sufficient coverage is obtained, providing a potential solution to the problem of variable amounts of M. tuberculosis DNA in direct samples.

Dentro SC, Wedge DC, Van Loo P. 2017. Principles of Reconstructing the Subclonal Architecture of Cancers. Cold Spring Harb Perspect Med, 7 (8), pp. a026625-a026625. | Show Abstract | Read more

Most cancers evolve from a single founder cell through a series of clonal expansions that are driven by somatic mutations. These clonal expansions can lead to several coexisting subclones sharing subsets of mutations. Analysis of massively parallel sequencing data can infer a tumor's subclonal composition through the identification of populations of cells with shared mutations. We describe the principles that underlie subclonal reconstruction through single nucleotide variants (SNVs) or copy number alterations (CNAs) from bulk or single-cell sequencing. These principles include estimating the fraction of tumor cells for SNVs and CNAs, performing clustering of SNVs from single- and multisample cases, and single-cell sequencing. The application of subclonal reconstruction methods is providing key insights into tumor evolution, identifying subclonal driver mutations, patterns of parallel evolution and differences in mutational signatures between cellular populations, and characterizing the mechanisms of therapy resistance, spread, and metastasis.

Wain LV, Shrine N, Artigas MS, Erzurumluoglu AM, Noyvert B, Bossini-Castillo L, Obeidat M, Henry AP, Portelli MA, Hall RJ et al. 2017. Genome-wide association analyses for lung function and chronic obstructive pulmonary disease identify new loci and potential druggable targets. Nat Genet, 49 (3), pp. 416-425. | Citations: 4 (Scopus) | Show Abstract | Read more

Chronic obstructive pulmonary disease (COPD) is characterized by reduced lung function and is the third leading cause of death globally. Through genome-wide association discovery in 48,943 individuals, selected from extremes of the lung function distribution in UK Biobank, and follow-up in 95,375 individuals, we increased the yield of independent signals for lung function from 54 to 97. A genetic risk score was associated with COPD susceptibility (odds ratio per 1 s.d. of the risk score (∼6 alleles) (95% confidence interval) = 1.24 (1.20-1.27), P = 5.05 × 10(-49)), and we observed a 3.7-fold difference in COPD risk between individuals in the highest and lowest genetic risk score deciles in UK Biobank. The 97 signals show enrichment in genes for development, elastic fibers and epigenetic regulation pathways. We highlight targets for drugs and compounds in development for COPD and asthma (genes in the inositol phosphate metabolism pathway and CHRM3) and describe targets for potential drug repositioning from other clinical indications.

Duffy CW, Ba H, Assefa S, Ahouidi AD, Deh YB, Tandia A, Kirsebom FCM, Kwiatkowski DP, Conway DJ. 2017. Population genetic structure and adaptation of malaria parasites on the edge of endemic distribution Molecular Ecology, 26 (11), pp. 2880-2894. | Read more

Kasela S, Kisand K, Tserel L, Kaleviste E, Remm A, Fischer K, Esko T, Westra HJ, Fairfax BP, Makino S et al. 2017. Pathogenic implications for autoimmune mechanisms derived by comparative eQTL analysis of CD4+ versus CD8+ T cells. PLoS Genet, 13 (3), pp. e1006643. | Show Abstract | Read more

Inappropriate activation or inadequate regulation of CD4+ and CD8+ T cells may contribute to the initiation and progression of multiple autoimmune and inflammatory diseases. Studies on disease-associated genetic polymorphisms have highlighted the importance of biological context for many regulatory variants, which is particularly relevant in understanding the genetic regulation of the immune system and its cellular phenotypes. Here we show cell type-specific regulation of transcript levels of genes associated with several autoimmune diseases in CD4+ and CD8+ T cells including a trans-acting regulatory locus at chr12q13.2 containing the rs1131017 SNP in the RPS26 gene. Most remarkably, we identify a common missense variant in IL27, associated with type 1 diabetes that results in decreased functional activity of the protein and reduced expression levels of downstream IRF1 and STAT1 in CD4+ T cells only. Altogether, our results indicate that eQTL mapping in purified T cells provides novel functional insights into polymorphisms and pathways associated with autoimmune diseases.

Goh C, Knight JC. 2017. Enhanced understanding of the host-pathogen interaction in sepsis: new opportunities for omic approaches. Lancet Respir Med, 5 (3), pp. 212-223. | Show Abstract | Read more

Progress in sepsis research has been severely hampered by a heterogeneous disease phenotype, limiting the interpretation of clinical trials and the development of effective therapeutic interventions. Application of omics-based methodologies is advancing understanding of the dysregulated host immune response to infection in sepsis. However, the frequently elusive nature of the infecting organism in sepsis has limited efforts to understand the effect of disease heterogeneity involving the pathogen. Recent advances in nucleic acid sequencing-based pathogen analysis provide the opportunity for more accurate and comprehensive microbiological diagnosis. In this Review, we explore how better understanding of the host-pathogen interaction can substantially enhance, and in turn benefit from, current and future application of omics-based approaches to understand the host response in sepsis. We illustrate this using recent work accounting for heterogeneity involving the pathogen. We propose that there is a timely opportunity to further resolve sepsis heterogeneity by considering host-pathogen interactions, enabling progress towards a precision medicine approach.

Wang L, Ko ER, Gilchrist JJ, Pittman KJ, Rautanen A, Pirinen M, Thompson JW, Dubois LG, Langley RJ, Jaslow SL et al. 2017. Human genetic and metabolite variation reveals that methylthioadenosine is a prognostic biomarker and an inflammatory regulator in sepsis. Sci Adv, 3 (3), pp. e1602096. | Citations: 1 (Web of Science Lite) | Show Abstract | Read more

Sepsis is a deleterious inflammatory response to infection with high mortality. Reliable sepsis biomarkers could improve diagnosis, prognosis, and treatment. Integration of human genetics, patient metabolite and cytokine measurements, and testing in a mouse model demonstrate that the methionine salvage pathway is a regulator of sepsis that can accurately predict prognosis in patients. Pathway-based genome-wide association analysis of nontyphoidal Salmonella bacteremia showed a strong enrichment for single-nucleotide polymorphisms near the components of the methionine salvage pathway. Measurement of the pathway's substrate, methylthioadenosine (MTA), in two cohorts of sepsis patients demonstrated increased plasma MTA in nonsurvivors. Plasma MTA was correlated with levels of inflammatory cytokines, indicating that elevated MTA marks a subset of patients with excessive inflammation. A machine-learning model combining MTA and other variables yielded approximately 80% accuracy (area under the curve) in predicting death. Furthermore, mice infected with Salmonella had prolonged survival when MTA was administered before infection, suggesting that manipulating MTA levels could regulate the severity of the inflammatory response. Our results demonstrate how combining genetic data, biomolecule measurements, and animal models can shape our understanding of disease and lead to new biomarkers for patient stratification and potential therapeutic targeting.

Franquinho F, Nogueira-Rodrigues J, Duarte JM, Esteves SS, Carter-Su C, Monaco AP, Molnár Z, Velayos-Baeza A, Brites P, Sousa MM. 2017. The Dyslexia-susceptibility Protein KIAA0319 Inhibits Axon Growth Through Smad2 Signaling. Cereb Cortex, 27 (3), pp. 1732-1747. | Citations: 1 (Scopus) | Show Abstract | Read more

KIAA0319 is a transmembrane protein associated with dyslexia with a presumed role in neuronal migration. Here we show that KIAA0319 expression is not restricted to the brain but also occurs in sensory and spinal cord neurons, increasing from early postnatal stages to adulthood and being downregulated by injury. This suggested that KIAA0319 participates in functions unrelated to neuronal migration. Supporting this hypothesis, overexpression of KIAA0319 repressed axon growth in hippocampal and dorsal root ganglia neurons; the intracellular domain of KIAA0319 was sufficient to elicit this effect. A similar inhibitory effect was observed in vivo as axon regeneration was impaired after transduction of sensory neurons with KIAA0319. Conversely, the deletion of Kiaa0319 in neurons increased neurite outgrowth in vitro and improved axon regeneration in vivo. At the mechanistic level, KIAA0319 engaged the JAK2-SH2B1 pathway to activate Smad2, which played a central role in KIAA0319-mediated repression of axon growth. In summary, we establish KIAA0319 as a novel player in axon growth and regeneration with the ability to repress the intrinsic growth potential of axons. This study describes a novel regulatory mechanism operating during peripheral nervous system and central nervous system axon growth, and offers novel targets for the development of effective therapies to promote axon regeneration.

Murray S, Pascoe B, Méric G, Mageiros L, Yahara K, Hitchings MD, Friedmann Y, Wilkinson TS, Gormley FJ, Mack D et al. 2017. Recombination-Mediated Host Adaptation by Avian Staphylococcus aureus. Genome Biol Evol, 9 (4), pp. 830-842. | Show Abstract | Read more

Staphylococcus aureus are globally disseminated among farmed chickens causing skeletal muscle infections, dermatitis, and septicaemia. The emergence of poultry-associated lineages has involved zoonotic transmission from humans to chickens but questions remain about the specific adaptations that promote proliferation of chicken pathogens. We characterized genetic variation in a population of genome-sequenced S. aureus isolates of poultry and human origin. Genealogical analysis identified a dominant poultry-associated sequence cluster within the CC5 clonal complex. Poultry and human CC5 isolates were significantly distinct from each other and more recombination events were detected in the poultry isolates. We identified 44 recombination events in 33 genes along the branch extending to the poultry-specific CC5 cluster, and 47 genes were found more often in CC5 poultry isolates compared with those from humans. Many of these gene sequences were common in chicken isolates from other clonal complexes suggesting horizontal gene transfer among poultry associated lineages. Consistent with functional predictions for putative poultry-associated genes, poultry isolates showed enhanced growth at 42 °C and greater erythrocyte lysis on chicken blood agar in comparison with human isolates. By combining phenotype information with evolutionary analyses of staphylococcal genomes, we provide evidence of adaptation, following a human-to-poultry host transition. This has important implications for the emergence and dissemination of new pathogenic clones associated with modern agriculture.

Lieberman S, Walsh T, Schechter M, Adar T, Goldin E, Beeri R, Sharon N, Baris H, Ben Avi L, Half E et al. 2017. Features of Patients With Hereditary Mixed Polyposis Syndrome Caused by Duplication of GREM1 and Implications for Screening and Surveillance. Gastroenterology, 152 (8), pp. 1876-1880.e1. | Citations: 1 (Web of Science Lite) | Show Abstract | Read more

Hereditary mixed polyposis syndrome is a rare colon cancer predisposition syndrome caused by a duplication of a noncoding sequence near the gremlin 1, DAN family BMP antagonist gene (GREM1) originally described in Ashkenazi Jews. Few families with GREM1 duplications have been described, so there are many questions about detection and management. We report 4 extended families with the duplication near GREM1 previously found in Ashkenazi Jews; 3 families were identified at cancer genetic clinics in Israel and 1 family was identified in a cohort of patients with familial colorectal cancer. Their clinical features include extracolonic tumors, onset of polyps in adolescence, and rapid progression of some polyps to advanced adenomas. One family met diagnostic criteria for Lynch syndrome. Expansion of the hereditary mixed polyposis syndrome phenotype can inform surveillance strategies for carriers of GREM1 duplications.

Telomeres Mendelian Randomization Collaboration, Haycock PC, Burgess S, Nounu A, Zheng J, Okoli GN, Bowden J, Wade KH, Timpson NJ, Evans DM et al. 2017. Association Between Telomere Length and Risk of Cancer and Non-Neoplastic Diseases: A Mendelian Randomization Study. JAMA Oncol, 3 (5), pp. 636-651. | Citations: 10 (Scopus) | Show Abstract | Read more

Importance: The causal direction and magnitude of the association between telomere length and incidence of cancer and non-neoplastic diseases is uncertain owing to the susceptibility of observational studies to confounding and reverse causation. Objective: To conduct a Mendelian randomization study, using germline genetic variants as instrumental variables, to appraise the causal relevance of telomere length for risk of cancer and non-neoplastic diseases. Data Sources: Genomewide association studies (GWAS) published up to January 15, 2015. Study Selection: GWAS of noncommunicable diseases that assayed germline genetic variation and did not select cohort or control participants on the basis of preexisting diseases. Of 163 GWAS of noncommunicable diseases identified, summary data from 103 were available. Data Extraction and Synthesis: Summary association statistics for single nucleotide polymorphisms (SNPs) that are strongly associated with telomere length in the general population. Main Outcomes and Measures: Odds ratios (ORs) and 95% confidence intervals (CIs) for disease per standard deviation (SD) higher telomere length due to germline genetic variation. Results: Summary data were available for 35 cancers and 48 non-neoplastic diseases, corresponding to 420 081 cases (median cases, 2526 per disease) and 1 093 105 controls (median, 6789 per disease). Increased telomere length due to germline genetic variation was generally associated with increased risk for site-specific cancers. The strongest associations (ORs [95% CIs] per 1-SD change in genetically increased telomere length) were observed for glioma, 5.27 (3.15-8.81); serous low-malignant-potential ovarian cancer, 4.35 (2.39-7.94); lung adenocarcinoma, 3.19 (2.40-4.22); neuroblastoma, 2.98 (1.92-4.62); bladder cancer, 2.19 (1.32-3.66); melanoma, 1.87 (1.55-2.26); testicular cancer, 1.76 (1.02-3.04); kidney cancer, 1.55 (1.08-2.23); and endometrial cancer, 1.31 (1.07-1.61). Associations were stronger for rarer cancers and at tissue sites with lower rates of stem cell division. There was generally little evidence of association between genetically increased telomere length and risk of psychiatric, autoimmune, inflammatory, diabetic, and other non-neoplastic diseases, except for coronary heart disease (OR, 0.78 [95% CI, 0.67-0.90]), abdominal aortic aneurysm (OR, 0.63 [95% CI, 0.49-0.81]), celiac disease (OR, 0.42 [95% CI, 0.28-0.61]) and interstitial lung disease (OR, 0.09 [95% CI, 0.05-0.15]). Conclusions and Relevance: It is likely that longer telomeres increase risk for several cancers but reduce risk for some non-neoplastic diseases, including cardiovascular diseases.

Perilla JR, Zhao G, Lu M, Ning J, Hou G, Byeon IL, Gronenborn AM, Polenova T, Zhang P. 2017. CryoEM Structure Refinement by Integrating NMR Chemical Shifts with Molecular Dynamics Simulations. J Phys Chem B, 121 (15), pp. 3853-3863. | Citations: 2 (Scopus) | Show Abstract | Read more

Single particle cryoEM has emerged as a powerful method for structure determination of proteins and complexes, complementing X-ray crystallography and NMR spectroscopy. Yet, for many systems, the resolution of cryoEM density map has been limited to 4-6 Å, which only allows for resolving bulky amino acids side chains, thus hindering accurate model building from the density map. On the other hand, experimental chemical shifts (CS) from solution and solid state MAS NMR spectra provide atomic level data for each amino acid within a molecule or a complex; however, structure determination of large complexes and assemblies based on NMR data alone remains challenging. Here, we present a novel integrated strategy to combine the highly complementary experimental data from cryoEM and NMR computationally by molecular dynamics simulations to derive an atomistic model, which is not attainable by either approach alone. We use the HIV-1 capsid protein (CA) C-terminal domain as well as the large capsid assembly to demonstrate the feasibility of this approach, termed NMR CS-biased cryoEM structure refinement.

Webb TR, Erdmann J, Stirrups KE, Stitziel NO, Masca NG, Jansen H, Kanoni S, Nelson CP, Ferrario PG, König IR et al. 2017. Systematic Evaluation of Pleiotropy Identifies 6 Further Loci Associated With Coronary Artery Disease. J Am Coll Cardiol, 69 (7), pp. 823-836. | Citations: 1 (Scopus) | Show Abstract | Read more

BACKGROUND: Genome-wide association studies have so far identified 56 loci associated with risk of coronary artery disease (CAD). Many CAD loci show pleiotropy; that is, they are also associated with other diseases or traits. OBJECTIVES: This study sought to systematically test if genetic variants identified for non-CAD diseases/traits also associate with CAD and to undertake a comprehensive analysis of the extent of pleiotropy of all CAD loci. METHODS: In discovery analyses involving 42,335 CAD cases and 78,240 control subjects we tested the association of 29,383 common (minor allele frequency >5%) single nucleotide polymorphisms available on the exome array, which included a substantial proportion of known or suspected single nucleotide polymorphisms associated with common diseases or traits as of 2011. Suggestive association signals were replicated in an additional 30,533 cases and 42,530 control subjects. To evaluate pleiotropy, we tested CAD loci for association with cardiovascular risk factors (lipid traits, blood pressure phenotypes, body mass index, diabetes, and smoking behavior), as well as with other diseases/traits through interrogation of currently available genome-wide association study catalogs. RESULTS: We identified 6 new loci associated with CAD at genome-wide significance: on 2q37 (KCNJ13-GIGYF2), 6p21 (C2), 11p15 (MRVI1-CTR9), 12q13 (LRP1), 12q24 (SCARB1), and 16q13 (CETP). Risk allele frequencies ranged from 0.15 to 0.86, and odds ratio per copy of the risk allele ranged from 1.04 to 1.09. Of 62 new and known CAD loci, 24 (38.7%) showed statistical association with a traditional cardiovascular risk factor, with some showing multiple associations, and 29 (47%) showed associations at p < 1 × 10(-4) with a range of other diseases/traits. CONCLUSIONS: We identified 6 loci associated with CAD at genome-wide significance. Several CAD loci show substantial pleiotropy, which may help us understand the mechanisms by which these loci affect CAD risk.

Kumar D, Puan KJ, Andiappan AK, Lee B, Westerlaken GH, Haase D, Melchiotti R, Li Z, Yusof N, Lum J et al. 2017. A functional SNP associated with atopic dermatitis controls cell type-specific methylation of the VSTM1 gene locus. Genome Med, 9 (1), pp. 18. | Show Abstract | Read more

BACKGROUND: Expression quantitative trait loci (eQTL) databases represent a valuable resource to link disease-associated SNPs to specific candidate genes whose gene expression is significantly modulated by the SNP under investigation. We previously identified signal inhibitory receptor on leukocytes-1 (SIRL-1) as a powerful regulator of human innate immune cell function. While it is constitutively high expressed on neutrophils, on monocytes the SIRL-1 surface expression varies strongly between individuals. The underlying mechanism of regulation, its genetic control as well as potential clinical implications had not been explored yet. METHODS: Whole blood eQTL data of a Chinese cohort was used to identify SNPs regulating the expression of VSTM1, the gene encoding SIRL-1. The genotype effect was validated by flow cytometry (cell surface expression), correlated with electrophoretic mobility shift assay (EMSA), chromatin immunoprecipitation (ChIP) and bisulfite sequencing (C-methylation) and its functional impact studied the inhibition of reactive oxygen species (ROS). RESULTS: We found a significant association of a single CpG-SNP, rs612529T/C, located in the promoter of VSTM1. Through flow cytometry analysis we confirmed that primarily in the monocytes the protein level of SIRL-1 is strongly associated with genotype of this SNP. In monocytes, the T allele of this SNP facilitates binding of the transcription factors YY1 and PU.1, of which the latter has been recently shown to act as docking site for modifiers of DNA methylation. In line with this notion rs612529T associates with a complete demethylation of the VSTM1 promoter correlating with the allele-specific upregulation of SIRL-1 expression. In monocytes, this upregulation strongly impacts the IgA-induced production of ROS by these cells. Through targeted association analysis we found a significant Meta P value of 1.14 × 10(-6) for rs612529 for association to atopic dermatitis (AD). CONCLUSION: Low expression of SIRL-1 on monocytes is associated with an increased risk for the manifestation of an inflammatory skin disease. It thus underlines the role of both the cell subset and this inhibitory immune receptor in maintaining immune homeostasis in the skin. Notably, the genetic regulation is achieved by a single CpG-SNP, which controls the overall methylation state of the promoter gene segment.

Cai N, Bigdeli TB, Kretzschmar WW, Li Y, Liang J, Hu J, Peterson RE, Bacanu S, Webb BT, Riley B et al. 2017. 11,670 whole-genome sequences representative of the Han Chinese population from the CONVERGE project. Sci Data, 4 pp. 170011. | Citations: 2 (Web of Science Lite) | Show Abstract | Read more

The China, Oxford and Virginia Commonwealth University Experimental Research on Genetic Epidemiology (CONVERGE) project on Major Depressive Disorder (MDD) sequenced 11,670 female Han Chinese at low-coverage (1.7X), providing the first large-scale whole genome sequencing resource representative of the largest ethnic group in the world. Samples are collected from 58 hospitals from 23 provinces around China. We are able to call 22 million high quality single nucleotide polymorphisms (SNP) from the nuclear genome, representing the largest SNP call set from an East Asian population to date. We use these variants for imputation of genotypes across all samples, and this has allowed us to perform a successful genome wide association study (GWAS) on MDD. The utility of these data can be extended to studies of genetic ancestry in the Han Chinese and evolutionary genetics when integrated with data from other populations. Molecular phenotypes, such as copy number variations and structural variations can be detected, quantified and analysed in similar ways.

Uimari O, Rahmioglu N, Nyholt DR, Vincent K, Missmer SA, Becker C, Morris AP, Montgomery GW, Zondervan KT. 2017. Genome-wide genetic analyses highlight mitogen-activated protein kinase (MAPK) signaling in the pathogenesis of endometriosis. Hum Reprod, 32 (4), pp. 780-793. | Citations: 3 (Web of Science Lite) | Show Abstract | Read more

Study question: Do genome-wide association study (GWAS) data for endometriosis provide insight into novel biological pathways associated with its pathogenesis? Summary answer: GWAS analysis uncovered multiple pathways that are statistically enriched for genetic association signals, analysis of Stage A disease highlighted a novel variant in MAP3K4, while top pathways significantly associated with all endometriosis and Stage A disease included several mitogen-activated protein kinase (MAPK)-related pathways. What is known already: Endometriosis is a complex disease with an estimated heritability of 50%. To date, GWAS revealed 10 genomic regions associated with endometriosis, explaining <4% of heritability, while half of the heritability is estimated to be due to common risk variants. Pathway analyses combine the evidence of single variants into gene-based measures, leveraging the aggregate effect of variants in genes and uncovering biological pathways involved in disease pathogenesis. Study design size, duration: Pathway analysis was conducted utilizing the International Endogene Consortium GWAS data, comprising 3194 surgically confirmed endometriosis cases and 7060 controls of European ancestry with genotype data imputed up to 1000 Genomes Phase three reference panel. GWAS was performed for all endometriosis cases and for Stage A (revised American Fertility Society (rAFS) I/II, n = 1686) and B (rAFS III/IV, n = 1364) cases separately. The identified significant pathways were compared with pathways previously investigated in the literature through candidate association studies. Participants/materials, setting, methods: The most comprehensive biological pathway databases, MSigDB (including BioCarta, KEGG, PID, SA, SIG, ST and GO) and PANTHER were utilized to test for enrichment of genetic variants associated with endometriosis. Statistical enrichment analysis was performed using the MAGENTA (Meta-Analysis Gene-set Enrichment of variaNT Associations) software. Main results and the role of chance: The first genome-wide association analysis for Stage A endometriosis revealed a novel locus, rs144240142 (P = 6.45 × 10-8, OR = 1.71, 95% CI = 1.23-2.37), an intronic single-nucleotide polymorphism (SNP) within MAP3K4. This SNP was not associated with Stage B disease (P = 0.086). MAP3K4 was also shown to be differentially expressed in eutopic endometrium between Stage A endometriosis cases and controls (P = 3.8 × 10-4), but not with Stage B disease (P = 0.26). A total of 14 pathways enriched with genetic endometriosis associations were identified (false discovery rate (FDR)-P < 0.05). The pathways associated with any endometriosis were Grb2-Sos provides linkage to MAPK signaling for integrins pathway (P = 2.8 × 10-5, FDR-P = 3.0 × 10-3), Wnt signaling (P = 0.026, FDR-P = 0.026) and p130Cas linkage to MAPK signaling for integrins pathway (P = 6.0 × 10-4, FDR-P = 0.029); with Stage A endometriosis: extracellular signal-regulated kinase (ERK)1 ERK2 MAPK (P = 5.0 × 10-4, FDR-P = 5.0 × 10-4) and with Stage B endometriosis: two overlapping pathways that related to extracellular matrix biology-Core matrisome (P = 1.4 × 10-3, FDR-P = 0.013) and ECM glycoproteins (P = 1.8 × 10-3, FDR-P = 7.1 × 10-3). Genes arising from endometriosis candidate gene studies performed to date were enriched for Interleukin signaling pathway (P = 2.3 × 10-12), Apoptosis signaling pathway (P = 9.7 × 10-9) and Gonadotropin releasing hormone receptor pathway (P = 1.2 × 10-6); however, these pathways did not feature in the results based on GWAS data. Large scale data: Not applicable. Limitations, reasons for caution: The analysis is restricted to (i) variants in/near genes that can be assigned to pathways, excluding intergenic variants; (ii) the gene-based pathway definition as registered in the databases; (iii) women of European ancestry. Wider implications of the findings: The top ranked pathways associated with overall and Stage A endometriosis in particular involve integrin-mediated MAPK activation and intracellular ERK/MAPK acting downstream in the MAPK cascade, both acting in the control of cell division, gene expression, cell movement and survival. Other top enriched pathways in Stage B disease include ECM glycoprotein pathways important for extracellular structure and biochemical support. The results highlight the need for increased efforts to understand the functional role of these pathways in endometriosis pathogenesis, including the investigation of the biological effects of the genetic variants on downstream molecular processes in tissue relevant to endometriosis. Additionally, our results offer further support for the hypothesis of at least partially distinct causal pathophysiology for minimal/mild (rAFS I/II) vs. moderate/severe (rAFS III/IV) endometriosis. Study funding/competing interest(s): The genome-wide association data and Wellcome Trust Case Control Consortium (WTCCC) were generated through funding from the Wellcome Trust (WT084766/Z/08/Z, 076113 and 085475) and the National Health and Medical Research Council (NHMRC) of Australia (241944, 339462, 389927, 389875, 389891, 389892, 389938, 443036, 442915, 442981, 496610, 496739, 552485 and 552498). N.R. was funded by a grant from the Medical Research Council UK (MR/K011480/1). A.P.M. is a Wellcome Trust Senior Fellow in Basic Biomedical Science (grant WT098017). All authors declare there are no conflicts of interest.

Alves E, Salman AM, Leoratti F, Lopez-Camacho C, Viveros-Sandoval ME, Lall A, El-Turabi A, Bachmann MF, Hill AV, Janse CJ et al. 2017. Evaluation of Plasmodium vivax Cell-Traversal Protein for Ookinetes and Sporozoites as a Preerythrocytic P. vivax Vaccine. Clin Vaccine Immunol, 24 (4), pp. e00501-16-e00501-16. | Citations: 1 (Web of Science Lite) | Show Abstract | Read more

Four different vaccine platforms, each targeting the human malaria parasite Plasmodium vivax cell-traversal protein for ookinetes and sporozoites (PvCelTOS), were generated and assessed for protective efficacy. These platforms consisted of a recombinant chimpanzee adenoviral vector 63 (ChAd63) expressing PvCelTOS (Ad), a recombinant modified vaccinia virus Ankara expressing PvCelTOS (MVA), PvCelTOS conjugated to bacteriophage Qβ virus-like particles (VLPs), and a recombinant PvCelTOS protein expressed in eukaryotic HEK293T cells (protein). Inbred BALB/c mice and outbred CD-1 mice were immunized using the following prime-boost regimens: Ad-MVA, Ad-VLPs, and Ad-protein. Protective efficacy against sporozoite challenge was assessed after immunization using a novel chimeric rodent Plasmodium berghei parasite (Pb-PvCelTOS). This chimeric parasite expresses P. vivax CelTOS in place of the endogenous P. berghei CelTOS and produces fully infectious sporozoites. A single Ad immunization in BALB/c and CD-1 mice induced anti-PvCelTOS antibodies which were boosted efficiently using MVA, VLP, or protein immunization. PvCelTOS-specific gamma interferon- and tumor necrosis factor alpha-producing CD8(+) T cells were induced at high frequencies by all prime-boost regimens in BALB/c mice but not in CD-1 mice; in CD-1 mice, they were only marginally increased after boosting with MVA. Despite the induction of anti-PvCelTOS antibodies and PvCelTOS-specific CD8(+) T-cell responses, only low levels of protective efficacy against challenge with Pb-PvCelTOS sporozoites were obtained using any immunization strategy. In BALB/c mice, no immunization regimens provided significant protection against a Pb-PvCelTOS chimeric sporozoite challenge. In CD-1 mice, modest protective efficacy against challenge with chimeric P. berghei sporozoites expressing either PvCelTOS or P. falciparum CelTOS was observed using the Ad-protein vaccination regimen.

McCarthy MI. 2017. Painting a new picture of personalised medicine for diabetes. Diabetologia, pp. 1-7. | Show Abstract | Read more

The current focus on delivery of personalised (or precision) medicine reflects the expectation that developments in genomics, imaging and other domains will extend our diagnostic and prognostic capabilities, and enable more effective targeting of current and future preventative and therapeutic options. The clinical benefits of this approach are already being realised in rare diseases and cancer but the impact on management of complex diseases, such as type 2 diabetes, remains limited. This may reflect reliance on inappropriate models of disease architecture, based around rare, high-impact genetic and environmental exposures that are poorly suited to our emerging understanding of type 2 diabetes. This review proposes an alternative 'palette' model, centred on a molecular taxonomy that focuses on positioning an individual with respect to the major pathophysiological processes that contribute to diabetes risk and progression. This model anticipates that many individuals with diabetes will have multiple parallel defects that affect several of these processes. One corollary of this model is that research efforts should, at least initially, be targeted towards identifying and characterising individuals whose adverse metabolic trajectory is dominated by perturbation in a restricted set of processes.

Rutledge GG, Böhme U, Sanders M, Reid AJ, Cotton JA, Maiga-Ascofare O, Djimdé AA, Apinjoh TO, Amenga-Etego L, Manske M et al. 2017. Plasmodium malariae and P. ovale genomes provide insights into malaria parasite evolution. Nature, 542 (7639), pp. 101-104. | Citations: 7 (Web of Science Lite) | Show Abstract | Read more

Elucidation of the evolutionary history and interrelatedness of Plasmodium species that infect humans has been hampered by a lack of genetic information for three human-infective species: P. malariae and two P. ovale species (P. o. curtisi and P. o. wallikeri). These species are prevalent across most regions in which malaria is endemic and are often undetectable by light microscopy, rendering their study in human populations difficult. The exact evolutionary relationship of these species to the other human-infective species has been contested. Using a new reference genome for P. malariae and a manually curated draft P. o. curtisi genome, we are now able to accurately place these species within the Plasmodium phylogeny. Sequencing of a P. malariae relative that infects chimpanzees reveals similar signatures of selection in the P. malariae lineage to another Plasmodium lineage shown to be capable of colonization of both human and chimpanzee hosts. Molecular dating suggests that these host adaptations occurred over similar evolutionary timescales. In addition to the core genome that is conserved between species, differences in gene content can be linked to their specific biology. The genome suggests that P. malariae expresses a family of heterodimeric proteins on its surface that have structural similarities to a protein crucial for invasion of red blood cells. The data presented here provide insight into the evolution of the Plasmodium genus as a whole.

Leedham SJ. 2017. MAP(K)ing the Path to Stem Cell Quiescence and the Elusive Enteroendocrine Cell Cell Stem Cell, 20 (2), pp. 153-154. | Show Abstract | Read more

© 2017 The existence and interaction of proliferating and quiescent intestinal stem cells have been debated since their discovery in the 1970s. In this issue of Cell Stem Cell, using murine intestinal organoids, Basak et al. (2017) induce stem cell quiescence by selective inhibition of EGF/MAPK signaling and define culture conditions that direct differentiation to the enteroendocrine lineage.

Bliss CM, Drammeh A, Bowyer G, Sanou GS, Jagne YJ, Ouedraogo O, Edwards NJ, Tarama C, Ouedraogo N, Ouedraogo M et al. 2017. Viral Vector Malaria Vaccines Induce High-Level T Cell and Antibody Responses in West African Children and Infants. Mol Ther, 25 (2), pp. 547-559. | Show Abstract | Read more

Heterologous prime-boosting with viral vectors encoding the pre-erythrocytic antigen thrombospondin-related adhesion protein fused to a multiple epitope string (ME-TRAP) induces CD8(+) T cell-mediated immunity to malaria sporozoite challenge in European malaria-naive and Kenyan semi-immune adults. This approach has yet to be evaluated in children and infants. We assessed this vaccine strategy among 138 Gambian and Burkinabe children in four cohorts: 2- to 6-year olds in The Gambia, 5- to 17-month-olds in Burkina Faso, and 5- to 12-month-olds and 10-week-olds in The Gambia. We assessed induction of cellular immunity, taking into account the distinctive hematological status of young infants, and characterized the antibody response to vaccination. T cell responses peaked 7 days after boosting with modified vaccinia virus Ankara (MVA), with highest responses in infants aged 10 weeks at priming. Incorporating lymphocyte count into the calculation of T cell responses facilitated a more physiologically relevant comparison of cellular immunity across different age groups. Both CD8(+) and CD4(+) T cells secreted cytokines. Induced antibodies were up to 20-fold higher in all groups compared with Gambian and United Kingdom (UK) adults, with comparable or higher avidity. This immunization regimen elicited strong immune responses, particularly in young infants, supporting future evaluation of efficacy in this key target age group for a malaria vaccine.

Cook JP, Mahajan A, Morris AP. 2017. Guidance for the utility of linear models in meta-analysis of genetic association studies of binary phenotypes. Eur J Hum Genet, 25 (2), pp. 240-245. | Show Abstract | Read more

Linear mixed models are increasingly used for the analysis of genome-wide association studies (GWAS) of binary phenotypes because they can efficiently and robustly account for population stratification and relatedness through inclusion of random effects for a genetic relationship matrix. However, the utility of linear (mixed) models in the context of meta-analysis of GWAS of binary phenotypes has not been previously explored. In this investigation, we present simulations to compare the performance of linear and logistic regression models under alternative weighting schemes in a fixed-effects meta-analysis framework, considering designs that incorporate variable case-control imbalance, confounding factors and population stratification. Our results demonstrate that linear models can be used for meta-analysis of GWAS of binary phenotypes, without loss of power, even in the presence of extreme case-control imbalance, provided that one of the following schemes is used: (i) effective sample size weighting of Z-scores or (ii) inverse-variance weighting of allelic effect sizes after conversion onto the log-odds scale. Our conclusions thus provide essential recommendations for the development of robust protocols for meta-analysis of binary phenotypes with linear models.

Liley J, Todd JA, Wallace C. 2017. A method for identifying genetic heterogeneity within phenotypically defined disease subgroups. Nat Genet, 49 (2), pp. 310-316. | Show Abstract | Read more

Many common diseases show wide phenotypic variation. We present a statistical method for determining whether phenotypically defined subgroups of disease cases represent different genetic architectures, in which disease-associated variants have different effect sizes in two subgroups. Our method models the genome-wide distributions of genetic association statistics with mixture Gaussians. We apply a global test without requiring explicit identification of disease-associated variants, thus maximizing power in comparison to standard variant-by-variant subgroup analysis. Where evidence for genetic subgrouping is found, we present methods for post hoc identification of the contributing genetic variants. We demonstrate the method on a range of simulated and test data sets, for which expected results are already known. We investigate subgroups of individuals with type 1 diabetes (T1D) defined by autoantibody positivity, establishing evidence for differential genetic architecture with positivity for thyroid-peroxidase-specific antibody, driven generally by variants in known T1D-associated genomic regions.

McCarthy MI. 2017. Genetics of T2DM in 2016: Biological and translational insights from T2DM genetics. Nat Rev Endocrinol, 13 (2), pp. 71-72. | Read more

Haghikia A, Dendrou CA, Schneider R, Grüter T, Postert T, Matzke M, Stephanik H, Fugger L, Gold R. 2017. Severe B-cell-mediated CNS disease secondary to alemtuzumab therapy. Lancet Neurol, 16 (2), pp. 104-106. | Citations: 2 (Web of Science Lite) | Read more

Hamblin A, Wordsworth S, Fermont JM, Page S, Kaur K, Camps C, Kaisaki P, Gupta A, Talbot D, Middleton M et al. 2017. Clinical applicability and cost of a 46-gene panel for genomic analysis of solid tumours: Retrospective validation and prospective audit in the UK National Health Service. PLoS Med, 14 (2), pp. e1002230. | Show Abstract | Read more

BACKGROUND: Single gene tests to predict whether cancers respond to specific targeted therapies are performed increasingly often. Advances in sequencing technology, collectively referred to as next generation sequencing (NGS), mean the entire cancer genome or parts of it can now be sequenced at speed with increased depth and sensitivity. However, translation of NGS into routine cancer care has been slow. Healthcare stakeholders are unclear about the clinical utility of NGS and are concerned it could be an expensive addition to cancer diagnostics, rather than an affordable alternative to single gene testing. METHODS AND FINDINGS: We validated a 46-gene hotspot cancer panel assay allowing multiple gene testing from small diagnostic biopsies. From 1 January 2013 to 31 December 2013, solid tumour samples (including non-small-cell lung carcinoma [NSCLC], colorectal carcinoma, and melanoma) were sequenced in the context of the UK National Health Service from 351 consecutively submitted prospective cases for which treating clinicians thought the patient had potential to benefit from more extensive genetic analysis. Following histological assessment, tumour-rich regions of formalin-fixed paraffin-embedded (FFPE) sections underwent macrodissection, DNA extraction, NGS, and analysis using a pipeline centred on Torrent Suite software. With a median turnaround time of seven working days, an integrated clinical report was produced indicating the variants detected, including those with potential diagnostic, prognostic, therapeutic, or clinical trial entry implications. Accompanying phenotypic data were collected, and a detailed cost analysis of the panel compared with single gene testing was undertaken to assess affordability for routine patient care. Panel sequencing was successful for 97% (342/351) of tumour samples in the prospective cohort and showed 100% concordance with known mutations (detected using cobas assays). At least one mutation was identified in 87% (296/342) of tumours. A locally actionable mutation (i.e., available targeted treatment or clinical trial) was identified in 122/351 patients (35%). Forty patients received targeted treatment, in 22/40 (55%) cases solely due to use of the panel. Examination of published data on the potential efficacy of targeted therapies showed theoretically actionable mutations (i.e., mutations for which targeted treatment was potentially appropriate) in 66% (71/107) and 39% (41/105) of melanoma and NSCLC patients, respectively. At a cost of £339 (US$449) per patient, the panel was less expensive locally than performing more than two or three single gene tests. Study limitations include the use of FFPE samples, which do not always provide high-quality DNA, and the use of "real world" data: submission of cases for sequencing did not always follow clinical guidelines, meaning that when mutations were detected, patients were not always eligible for targeted treatments on clinical grounds. CONCLUSIONS: This study demonstrates that more extensive tumour sequencing can identify mutations that could improve clinical decision-making in routine cancer care, potentially improving patient outcomes, at an affordable level for healthcare providers.

Macintyre G, Van Loo P, Corcoran NM, Wedge DC, Markowetz F, Hovens CM. 2017. How Subclonal Modeling Is Changing the Metastatic Paradigm. Clin Cancer Res, 23 (3), pp. 630-635. | Citations: 1 (Web of Science Lite) | Show Abstract | Read more

A concerted effort to sequence matched primary and metastatic tumors is vastly improving our ability to understand metastasis in humans. Compelling evidence has emerged that supports the existence of diverse and surprising metastatic patterns. Enhancing these efforts is a new class of algorithms that facilitate high-resolution subclonal modeling of metastatic spread. Here we summarize how subclonal models of metastasis are influencing the metastatic paradigm. Clin Cancer Res; 23(3); 630-5. ©2016 AACR.

Warren HR, Evangelou E, Cabrera CP, Gao H, Ren M, Mifsud B, Ntalla I, Surendran P, Liu C, Cook JP et al. 2017. Genome-wide association analysis identifies novel blood pressure loci and offers biological insights into cardiovascular risk. Nat Genet, 49 (3), pp. 403-415. | Citations: 10 (Scopus) | Show Abstract | Read more

Elevated blood pressure is the leading heritable risk factor for cardiovascular disease worldwide. We report genetic association of blood pressure (systolic, diastolic, pulse pressure) among UK Biobank participants of European ancestry with independent replication in other cohorts, and robust validation of 107 independent loci. We also identify new independent variants at 11 previously reported blood pressure loci. In combination with results from a range of in silico functional analyses and wet bench experiments, our findings highlight new biological pathways for blood pressure regulation enriched for genes expressed in vascular tissues and identify potential therapeutic targets for hypertension. Results from genetic risk score models raise the possibility of a precision medicine approach through early lifestyle intervention to offset the impact of blood pressure-raising genetic variants on future cardiovascular disease risk.

McCarthy MI, MacArthur DG. 2017. Human disease genomics: from variants to biology. Genome Biol, 18 (1), pp. 20. | Show Abstract | Read more

We summarize the remarkable progress that has been made in the identification and functional characterization of DNA sequence variants associated with disease.

Chang HH, Worby CJ, Yeka A, Nankabirwa J, Kamya MR, Staedke SG, Dorsey G, Murphy M, Neafsey DE, Jeffreys AE et al. 2017. THE REAL McCOIL: A method for the concurrent estimation of the complexity of infection and SNP allele frequency for malaria parasites. PLoS Comput Biol, 13 (1), pp. e1005348. | Citations: 1 (Web of Science Lite) | Show Abstract | Read more

As many malaria-endemic countries move towards elimination of Plasmodium falciparum, the most virulent human malaria parasite, effective tools for monitoring malaria epidemiology are urgent priorities. P. falciparum population genetic approaches offer promising tools for understanding transmission and spread of the disease, but a high prevalence of multi-clone or polygenomic infections can render estimation of even the most basic parameters, such as allele frequencies, challenging. A previous method, COIL, was developed to estimate complexity of infection (COI) from single nucleotide polymorphism (SNP) data, but relies on monogenomic infections to estimate allele frequencies or requires external allele frequency data which may not available. Estimates limited to monogenomic infections may not be representative, however, and when the average COI is high, they can be difficult or impossible to obtain. Therefore, we developed THE REAL McCOIL, Turning HEterozygous SNP data into Robust Estimates of ALelle frequency, via Markov chain Monte Carlo, and Complexity Of Infection using Likelihood, to incorporate polygenomic samples and simultaneously estimate allele frequency and COI. This approach was tested via simulations then applied to SNP data from cross-sectional surveys performed in three Ugandan sites with varying malaria transmission. We show that THE REAL McCOIL consistently outperforms COIL on simulated data, particularly when most infections are polygenomic. Using field data we show that, unlike with COIL, we can distinguish epidemiologically relevant differences in COI between and within these sites. Surprisingly, for example, we estimated high average COI in a peri-urban subregion with lower transmission intensity, suggesting that many of these cases were imported from surrounding regions with higher transmission intensity. THE REAL McCOIL therefore provides a robust tool for understanding the molecular epidemiology of malaria across transmission settings.

Liston A, Todd JA, Lagou V. 2017. Beta-Cell Fragility As a Common Underlying Risk Factor in Type 1 and Type 2 Diabetes. Trends Mol Med, 23 (2), pp. 181-194. | Citations: 2 (Scopus) | Show Abstract | Read more

Type 1 and type 2 diabetes are distinct clinical entities primarily driven by autoimmunity and metabolic dysfunction, respectively. However, there is a growing appreciation that they may share an etiopathological factor, namely the role of variation in beta-cell sensitivity to stress factors. Increased sensitivity increases the risk of beta-cell death or insulin secretion dysfunction. The beta-cell fragility model proposes that this variation contributes to the risk of developing either type 1 or type 2 diabetes, in the presence of immunological and/or metabolic stress factors. Therapeutics that increase the resistance of beta cells to these factors and decreasing fragility may constitute a new class of anti-diabetogenics, with potential use across both diseases.

Vonk JJ, Yeshaw WM, Pinto F, Faber AI, Lahaye LL, Kanon B, van der Zwaag M, Velayos-Baeza A, Freire R, van IJzendoorn SC et al. 2017. Drosophila Vps13 Is Required for Protein Homeostasis in the Brain. PLoS One, 12 (1), pp. e0170106. | Show Abstract | Read more

Chorea-Acanthocytosis is a rare, neurodegenerative disorder characterized by progressive loss of locomotor and cognitive function. It is caused by loss of function mutations in the Vacuolar Protein Sorting 13A (VPS13A) gene, which is conserved from yeast to human. The consequences of VPS13A dysfunction in the nervous system are still largely unspecified. In order to study the consequences of VPS13A protein dysfunction in the ageing central nervous system we characterized a Drosophila melanogaster Vps13 mutant line. The Drosophila Vps13 gene encoded a protein of similar size as human VPS13A. Our data suggest that Vps13 is a peripheral membrane protein located to endosomal membranes and enriched in the fly head. Vps13 mutant flies showed a shortened life span and age associated neurodegeneration. Vps13 mutant flies were sensitive to proteotoxic stress and accumulated ubiquitylated proteins. Levels of Ref(2)P, the Drosophila orthologue of p62, were increased and protein aggregates accumulated in the central nervous system. Overexpression of the human Vps13A protein in the mutant flies partly rescued apparent phenotypes. This suggests a functional conservation of human VPS13A and Drosophila Vps13. Our results demonstrate that Vps13 is essential to maintain protein homeostasis in the larval and adult Drosophila brain. Drosophila Vps13 mutants are suitable to investigate the function of Vps13 in the brain, to identify genetic enhancers and suppressors and to screen for potential therapeutic targets for Chorea-Acanthocytosis.

Neville MJ, Lee W, Humburg P, Wong D, Barnardo M, Karpe F, Knight JC. 2017. High resolution HLA haplotyping by imputation for a British population bioresource. Hum Immunol, 78 (3), pp. 242-251. | Citations: 1 (Web of Science Lite) | Show Abstract | Read more

This study aimed to establish the occurrence and frequency of HLA alleles and haplotypes for a healthy British Caucasian population bioresource from Oxfordshire. We present the results of imputation from HLA SNP genotyping data using SNP2HLA for 5553 individuals from Oxford Biobank, defining one- and two-field alleles together with amino acid polymorphisms. We show that this achieves a high level of accuracy with validation using sequence-specific primer amplification PCR. We define six- and eight-locus HLA haplotypes for this population by Bayesian methods implemented using PHASE. We determine patterns of linkage disequilibrium and recombination for these individuals involving classical HLA loci and show how analysis within a haplotype block structure may be more tractable for imputed data. Our findings contribute to knowledge of HLA diversity in healthy populations and further validate future large-scale use of HLA imputation as an informative approach in population bioresources.

Carrat GR, Hu M, Nguyen-Tu MS, Chabosseau P, Gaulton KJ, van de Bunt M, Siddiq A, Falchi M, Thurner M, Canouil M et al. 2017. Decreased STARD10 Expression Is Associated with Defective Insulin Secretion in Humans and Mice. Am J Hum Genet, 100 (2), pp. 238-256. | Show Abstract | Read more

Genetic variants near ARAP1 (CENTD2) and STARD10 influence type 2 diabetes (T2D) risk. The risk alleles impair glucose-induced insulin secretion and, paradoxically but characteristically, are associated with decreased proinsulin:insulin ratios, indicating improved proinsulin conversion. Neither the identity of the causal variants nor the gene(s) through which risk is conferred have been firmly established. Whereas ARAP1 encodes a GTPase activating protein, STARD10 is a member of the steroidogenic acute regulatory protein (StAR)-related lipid transfer protein family. By integrating genetic fine-mapping and epigenomic annotation data and performing promoter-reporter and chromatin conformational capture (3C) studies in β cell lines, we localize the causal variant(s) at this locus to a 5 kb region that overlaps a stretch-enhancer active in islets. This region contains several highly correlated T2D-risk variants, including the rs140130268 indel. Expression QTL analysis of islet transcriptomes from three independent subject groups demonstrated that T2D-risk allele carriers displayed reduced levels of STARD10 mRNA, with no concomitant change in ARAP1 mRNA levels. Correspondingly, β-cell-selective deletion of StarD10 in mice led to impaired glucose-stimulated Ca(2+) dynamics and insulin secretion and recapitulated the pattern of improved proinsulin processing observed at the human GWAS signal. Conversely, overexpression of StarD10 in the adult β cell improved glucose tolerance in high fat-fed animals. In contrast, manipulation of Arap1 in β cells had no impact on insulin secretion or proinsulin conversion in mice. This convergence of human and murine data provides compelling evidence that the T2D risk associated with variation at this locus is mediated through reduction in STARD10 expression in the β cell.

Sandor C, Robertson P, Lang C, Heger A, Booth H, Vowles J, Witty L, Bowden R, Hu M, Cowley SA et al. 2017. Transcriptomic profiling of purified patient-derived dopamine neurons identifies convergent perturbations and therapeutics for Parkinson's disease. Hum Mol Genet, 26 (3), pp. 552-566. | Citations: 2 (Web of Science Lite) | Show Abstract | Read more

While induced pluripotent stem cell (iPSC) technologies enable the study of inaccessible patient cell types, cellular heterogeneity can confound the comparison of gene expression profiles between iPSC-derived cell lines. Here, we purified iPSC-derived human dopaminergic neurons (DaNs) using the intracellular marker, tyrosine hydroxylase. Once purified, the transcriptomic profiles of iPSC-derived DaNs appear remarkably similar to profiles obtained from mature post-mortem DaNs. Comparison of the profiles of purified iPSC-derived DaNs derived from Parkinson's disease (PD) patients carrying LRRK2 G2019S variants to controls identified significant functional convergence amongst differentially-expressed (DE) genes. The PD LRRK2-G2019S associated profile was positively matched with expression changes induced by the Parkinsonian neurotoxin rotenone and opposed by those induced by clioquinol, a compound with demonstrated therapeutic efficacy in multiple PD models. No functional convergence amongst DE genes was observed following a similar comparison using non-purified iPSC-derived DaN-containing populations, with cellular heterogeneity appearing a greater confound than genotypic background.

de Castro IJ, Budzak J, Di Giacinto ML, Ligammari L, Gokhan E, Spanos C, Moralli D, Richardson C, de Las Heras JI, Salatino S et al. 2017. Repo-Man/PP1 regulates heterochromatin formation in interphase. Nat Commun, 8 pp. 14048. | Citations: 1 (Web of Science Lite) | Show Abstract | Read more

Repo-Man is a protein phosphatase 1 (PP1) targeting subunit that regulates mitotic progression and chromatin remodelling. After mitosis, Repo-Man/PP1 remains associated with chromatin but its function in interphase is not known. Here we show that Repo-Man, via Nup153, is enriched on condensed chromatin at the nuclear periphery and at the edge of the nucleopore basket. Repo-Man/PP1 regulates the formation of heterochromatin, dephosphorylates H3S28 and it is necessary and sufficient for heterochromatin protein 1 binding and H3K27me3 recruitment. Using a novel proteogenomic approach, we show that Repo-Man is enriched at subtelomeric regions together with H2AZ and H3.3 and that depletion of Repo-Man alters the peripheral localization of a subset of these regions and alleviates repression of some polycomb telomeric genes. This study shows a role for a mitotic phosphatase in the regulation of the epigenetic landscape and gene expression in interphase.

Spencer AJ, Longley RJ, Gola A, Ulaszewska M, Lambe T, Hill AV. 2017. The Threshold of Protection from Liver-Stage Malaria Relies on a Fine Balance between the Number of Infected Hepatocytes and Effector CD8(+) T Cells Present in the Liver. J Immunol, 198 (5), pp. 2006-2016. | Citations: 1 (Scopus) | Show Abstract | Read more

Since the demonstration of sterile protection afforded by injection of irradiated sporozoites, CD8(+) T cells have been shown to play a significant role in protection from liver-stage malaria. This is, however, dependent on the presence of an extremely high number of circulating effector cells, thought to be necessary to scan, locate, and kill infected hepatocytes in the short time that parasites are present in the liver. We used an adoptive transfer model to elucidate the kinetics of the effector CD8(+) T cell response in the liver following Plasmodium berghei sporozoite challenge. Although effector CD8(+) T cells require <24 h to find, locate, and kill infected hepatocytes, active migration of Ag-specific CD8(+) T cells into the liver was not observed during the 2-d liver stage of infection, as divided cells were only detected from day 3 postchallenge. However, the percentage of donor cells recruited into division was shown to indicate the level of Ag presentation from infected hepatocytes. By titrating the number of transferred Ag-specific effector CD8(+) T cells and sporozoites, we demonstrate that achieving protection toward liver-stage malaria is reliant on CD8(+) T cells being able to locate infected hepatocytes, resulting in a protection threshold dependent on a fine balance between the number of infected hepatocytes and CD8(+) T cells present in the liver. With such a fine balance determining protection, achieving a high number of CD8(+) T cells will be critical to the success of a cell-mediated vaccine against liver-stage malaria.

Mägi R, Suleimanov YV, Clarke GM, Kaakinen M, Fischer K, Prokopenko I, Morris AP. 2017. SCOPA and META-SCOPA: software for the analysis and aggregation of genome-wide association studies of multiple correlated phenotypes. BMC Bioinformatics, 18 (1), pp. 25. | Citations: 1 (Scopus) | Show Abstract | Read more

BACKGROUND: Genome-wide association studies (GWAS) of single nucleotide polymorphisms (SNPs) have been successful in identifying loci contributing genetic effects to a wide range of complex human diseases and quantitative traits. The traditional approach to GWAS analysis is to consider each phenotype separately, despite the fact that many diseases and quantitative traits are correlated with each other, and often measured in the same sample of individuals. Multivariate analyses of correlated phenotypes have been demonstrated, by simulation, to increase power to detect association with SNPs, and thus may enable improved detection of novel loci contributing to diseases and quantitative traits. RESULTS: We have developed the SCOPA software to enable GWAS analysis of multiple correlated phenotypes. The software implements "reverse regression" methodology, which treats the genotype of an individual at a SNP as the outcome and the phenotypes as predictors in a general linear model. SCOPA can be applied to quantitative traits and categorical phenotypes, and can accommodate imputed genotypes under a dosage model. The accompanying META-SCOPA software enables meta-analysis of association summary statistics from SCOPA across GWAS. Application of SCOPA to two GWAS of high-and low-density lipoprotein cholesterol, triglycerides and body mass index, and subsequent meta-analysis with META-SCOPA, highlighted stronger association signals than univariate phenotype analysis at established lipid and obesity loci. The META-SCOPA meta-analysis also revealed a novel signal of association at genome-wide significance for triglycerides mapping to GPC5 (lead SNP rs71427535, p = 1.1x10(-8)), which has not been reported in previous large-scale GWAS of lipid traits. CONCLUSIONS: The SCOPA and META-SCOPA software enable discovery and dissection of multiple phenotype association signals through implementation of a powerful reverse regression approach.

Total publications on this page: 100

Total citations for publications on this page: 95