Recent publications

Sweeney TE, Perumal TM, Henao R, Nichols M, Howrylak JA, Choi AM, Bermejo-Martin JF, Almansa R, Tamayo E, Davenport EE et al. 2018. A community approach to mortality prediction in sepsis via gene expression analysis. Nat Commun, 9 (1), pp. 694. | Show Abstract | Read more

Improved risk stratification and prognosis prediction in sepsis is a critical unmet need. Clinical severity scores and available assays such as blood lactate reflect global illness severity with suboptimal performance, and do not specifically reveal the underlying dysregulation of sepsis. Here, we present prognostic models for 30-day mortality generated independently by three scientific groups by using 12 discovery cohorts containing transcriptomic data collected from primarily community-onset sepsis patients. Predictive performance is validated in five cohorts of community-onset sepsis patients in which the models show summary AUROCs ranging from 0.765-0.89. Similar performance is observed in four cohorts of hospital-acquired sepsis. Combining the new gene-expression-based prognostic models with prior clinical severity scores leads to significant improvement in prediction of 30-day mortality as measured via AUROC and net reclassification improvement index These models provide an opportunity to develop molecular bedside tests that may improve risk stratification and mortality prediction in patients with sepsis.

Smith D, Magri A, Bonsall D, Ip CL, Trebes A, Brown A, Piazza P, Bowden R, Nguyen D, STOP-HCV Consortium et al. 2018. Resistance analysis of genotype 3 HCV indicates subtypes inherently resistant to NS5A inhibitors. Hepatology, | Show Abstract | Read more

BACKGROUND AND AIMS: HCV genotype-3 (gt3) is highly prevalent globally, with non-gt3a subtypes common in Southeast Asia. Resistance-associated substitutions (RASs) have been shown to play a role in treatment failure. However, the role of RASs in gt3 is not well understood. We report the prevalence of RASs in a cohort of directly acting antivirals (DAA) treatment-naïve gt3 infected patients, including those with rarer subtypes and evaluate the effect of these RAS on DAAs in-vitro. METHODS: Baseline samples from 496 gt3 patients enrolled in the BOSON clinical trial were analysed by next-generation sequencing after probe-based enrichment for HCV. Whole viral genomes were analysed for the presence of RASs to approved DAAs. The resistance phenotype of RASs in combination to daclatasvir, velpatasvir, pibrentasvir, elbasvir and sofosbuvir was measured using the S52 ΔN gt3a replicon model. RESULTS: The NS5A A30K and Y93H substitutions were the most common at 8.8% (n=44) and 12.2% (n=61) respectively and showed a 10-fold and 11-fold increase in EC50 for daclatasvir compared to the unmodified replicon. Paired RASs (A30K + L31M and A30K + Y93H) were identified in 18 patients (9 of each pair); these combinations were shown to be highly resistant to daclatasvir, velpatasvir, elbasvir and pibrentasvir. The A30K + L31M combination was found in all gt3b and gt3g samples. CONCLUSION: Our study reveals high frequencies of RASs to NS5A Inhibitors in gt3 HCV. The paired A30K + L31M substitutions occur in all patients with gt3b and gt3g virus. In vitro analysis suggests that these subtypes may be inherently resistant to all approved NS5A Inhibitors for gt3 HCV. This article is protected by copyright. All rights reserved.

Zheng Y, Sethi R, Mangala LS, Taylor C, Goldsmith J, Wang M, Masuda K, Karaminejadranjbar M, Mannion D, Miranda F et al. 2018. Tuning microtubule dynamics to enhance cancer therapy by modulating FER-mediated CRMP2 phosphorylation. Nat Commun, 9 (1), pp. 476. | Show Abstract | Read more

Though used widely in cancer therapy, paclitaxel only elicits a response in a fraction of patients. A strong determinant of paclitaxel tumor response is the state of microtubule dynamic instability. However, whether the manipulation of this physiological process can be controlled to enhance paclitaxel response has not been tested. Here, we show a previously unrecognized role of the microtubule-associated protein CRMP2 in inducing microtubule bundling through its carboxy terminus. This activity is significantly decreased when the FER tyrosine kinase phosphorylates CRMP2 at Y479 and Y499. The crystal structures of wild-type CRMP2 and CRMP2-Y479E reveal how mimicking phosphorylation prevents tetramerization of CRMP2. Depletion of FER or reducing its catalytic activity using sub-therapeutic doses of inhibitors increases paclitaxel-induced microtubule stability and cytotoxicity in ovarian cancer cells and in vivo. This work provides a rationale for inhibiting FER-mediated CRMP2 phosphorylation to enhance paclitaxel on-target activity for cancer therapy.

Attar M, Sharma E, Li S, Bryer C, Cubitt L, Broxholme J, Lockstone H, Kinchen J, Simmons A, Piazza P et al. 2018. A practical solution for preserving single cells for RNA sequencing. Sci Rep, 8 (1), pp. 2151. | Show Abstract | Read more

The design and implementation of single-cell experiments is often limited by their requirement for fresh starting material. We have adapted a method for histological tissue fixation using dithio-bis(succinimidyl propionate) (DSP), or Lomant's Reagent, to stabilise cell samples for single-cell transcriptomic applications. DSP is a reversible cross-linker of free amine groups that has previously been shown to preserve tissue integrity for histology while maintaining RNA integrity and yield in bulk RNA extractions. Although RNA-seq data from DSP-fixed single cells appears to be prone to characteristic artefacts, such as slightly reduced yield of cDNA and a detectable 3' bias in comparison with fresh cells, cell preservation using DSP does not appear to substantially reduce RNA complexity at the gene level. In addition, there is evidence that instantaneous fixation of cells can reduce inter-cell technical variability. The ability of DSP-fixed cells to retain commonly used dyes, such as propidium iodide, enables the tracking of experimental sub-populations and the recording of cell viability at the point of fixation. Preserving cells using DSP will remove several barriers in the staging of single-cell experiments, including the transport of samples and the scheduling of shared equipment for downstream single-cell isolation and processing.

Robbe P, Popitsch N, Knight SJL, Antoniou P, Becq J, He M, Kanapin A, Samsonova A, Vavoulis DV, Ross MT et al. 2018. Clinical whole-genome sequencing from routine formalin-fixed, paraffin-embedded specimens: pilot study for the 100,000 Genomes Project. Genet Med, | Show Abstract | Read more

PurposeFresh-frozen (FF) tissue is the optimal source of DNA for whole-genome sequencing (WGS) of cancer patients. However, it is not always available, limiting the widespread application of WGS in clinical practice. We explored the viability of using formalin-fixed, paraffin-embedded (FFPE) tissues, available routinely for cancer patients, as a source of DNA for clinical WGS.MethodsWe conducted a prospective study using DNAs from matched FF, FFPE, and peripheral blood germ-line specimens collected from 52 cancer patients (156 samples) following routine diagnostic protocols. We compared somatic variants detected in FFPE and matching FF samples.ResultsWe found the single-nucleotide variant agreement reached 71% across the genome and somatic copy-number alterations (CNAs) detection from FFPE samples was suboptimal (0.44 median correlation with FF) due to nonuniform coverage. CNA detection was improved significantly with lower reverse crosslinking temperature in FFPE DNA extraction (80 °C or 65 °C depending on the methods). Our final data showed somatic variant detection from FFPE for clinical decision making is possible. We detected 98% of clinically actionable variants (including 30/31 CNAs).ConclusionWe present the first prospective WGS study of cancer patients using FFPE specimens collected in a routine clinical environment proving WGS can be applied in the clinic.GENETICS in MEDICINE advance online publication, 1 February 2018; doi:10.1038/gim.2017.241.

Kvikstad EM, Piazza P, Taylor JC, Lunter G. 2018. A high throughput screen for active human transposable elements. BMC Genomics, 19 (1), pp. 115. | Show Abstract | Read more

BACKGROUND: Transposable elements (TEs) are mobile genetic sequences that randomly propagate within their host's genome. This mobility has the potential to affect gene transcription and cause disease. However, TEs are technically challenging to identify, which complicates efforts to assess the impact of TE insertions on disease. Here we present a targeted sequencing protocol and computational pipeline to identify polymorphic and novel TE insertions using next-generation sequencing: TE-NGS. The method simultaneously targets the three subfamilies that are responsible for the majority of recent TE activity (L1HS, AluYa5/8, and AluYb8/9) thereby obviating the need for multiple experiments and reducing the amount of input material required. RESULTS: Here we describe the laboratory protocol and detection algorithm, and a benchmark experiment for the reference genome NA12878. We demonstrate a substantial enrichment for on-target fragments, and high sensitivity and precision to both reference and NA12878-specific insertions. We report 17 previously unreported loci for this individual which are supported by orthogonal long-read evidence, and we identify 1470 polymorphic and novel TEs in 12 additional samples that were previously undocumented in databases of insertion polymorphisms. CONCLUSIONS: We anticipate that future applications of TE-NGS alongside exome sequencing of patients with sporadic disease will reduce the number of unresolved cases, and improve estimates of the contribution of TEs to human genetic disease.

Amato R, Pearson RD, Almagro-Garcia J, Amaratunga C, Lim P, Suon S, Sreng S, Drury E, Stalker J, Miotto O et al. 2018. Origins of the current outbreak of multidrug-resistant malaria in southeast Asia: a retrospective genetic study. Lancet Infect Dis, | Show Abstract | Read more

BACKGROUND: Antimalarial resistance is rapidly spreading across parts of southeast Asia where dihydroartemisinin-piperaquine is used as first-line treatment for Plasmodium falciparum malaria. The first published reports about resistance to antimalarial drugs came from western Cambodia in 2013. Here, we analyse genetic changes in the P falciparum population of western Cambodia in the 6 years before those reports. METHODS: We analysed genome sequence data on 1492 P falciparum samples from 11 locations across southeast Asia, including 464 samples collected in western Cambodia between 2007 and 2013. Different epidemiological origins of resistance were identified by haplotypic analysis of the kelch13 artemisinin resistance locus and the plasmepsin 2-3 piperaquine resistance locus. FINDINGS: We identified more than 30 independent origins of artemisinin resistance, of which the KEL1 lineage accounted for 140 (91%) of 154 parasites resistant to dihydroartemisinin-piperaquine. In 2008, KEL1 combined with PLA1, the major lineage associated with piperaquine resistance. By 2013, the KEL1/PLA1 co-lineage had reached a frequency of 63% (24/38) in western Cambodia and had spread to northern Cambodia. INTERPRETATION: The KEL1/PLA1 co-lineage emerged in the same year that dihydroartemisinin-piperaquine became the first-line antimalarial drug in western Cambodia and spread rapidly thereafter, displacing other artemisinin-resistant parasite lineages. These findings have important implications for management of the global health risk associated with the current outbreak of multidrug-resistant malaria in southeast Asia. FUNDING: Wellcome Trust, Bill & Melinda Gates Foundation, Medical Research Council, UK Department for International Development, and the Intramural Research Program of the National Institute of Allergy and Infectious Diseases.

Flannick J, Fuchsberger C, Mahajan A, Teslovich TM, Agarwala V, Gaulton KJ, Caulkins L, Koesterer R, Ma C, Moutsianas L et al. 2018. Erratum: Sequence data and association statistics from 12,940 type 2 diabetes cases and controls. Sci Data, 5 pp. 180002. | Show Abstract | Read more

This corrects the article DOI: 10.1038/sdata.2017.179.

Kotecha D, Breithardt G, Camm AJ, Lip GYH, Schotten U, Ahlsson A, Arnar D, Atar D, Auricchio A, Bax J et al. 2018. Integrating new approaches to atrial fibrillation management: the 6th AFNET/EHRA Consensus Conference. Europace, | Show Abstract | Read more

There are major challenges ahead for clinicians treating patients with atrial fibrillation (AF). The population with AF is expected to expand considerably and yet, apart from anticoagulation, therapies used in AF have not been shown to consistently impact on mortality or reduce adverse cardiovascular events. New approaches to AF management, including the use of novel technologies and structured, integrated care, have the potential to enhance clinical phenotyping or result in better treatment selection and stratified therapy. Here, we report the outcomes of the 6th Consensus Conference of the Atrial Fibrillation Network (AFNET) and the European Heart Rhythm Association (EHRA), held at the European Society of Cardiology Heart House in Sophia Antipolis, France, 17-19 January 2017. Sixty-two global specialists in AF and 13 industry partners met to develop innovative solutions based on new approaches to screening and diagnosis, enhancing integration of AF care, developing clinical pathways for treating complex patients, improving stroke prevention strategies, and better patient selection for heart rate and rhythm control. Ultimately, these approaches can lead to better outcomes for patients with AF.

Jun G, Manning A, Almeida M, Zawistowski M, Wood AR, Teslovich TM, Fuchsberger C, Feng S, Cingolani P, Gaulton KJ et al. 2018. Evaluating the contribution of rare variants to type 2 diabetes and related traits using pedigrees. Proc Natl Acad Sci U S A, 115 (2), pp. 379-384. | Show Abstract | Read more

A major challenge in evaluating the contribution of rare variants to complex disease is identifying enough copies of the rare alleles to permit informative statistical analysis. To investigate the contribution of rare variants to the risk of type 2 diabetes (T2D) and related traits, we performed deep whole-genome analysis of 1,034 members of 20 large Mexican-American families with high prevalence of T2D. If rare variants of large effect accounted for much of the diabetes risk in these families, our experiment was powered to detect association. Using gene expression data on 21,677 transcripts for 643 pedigree members, we identified evidence for large-effect rare-variant cis-expression quantitative trait loci that could not be detected in population studies, validating our approach. However, we did not identify any rare variants of large effect associated with T2D, or the related traits of fasting glucose and insulin, suggesting that large-effect rare variants account for only a modest fraction of the genetic risk of these traits in this sample of families. Reliable identification of large-effect rare variants will require larger samples of extended pedigrees or different study designs that further enrich for such variants.

Turcot V, Lu Y, Highland HM, Schurmann C, Justice AE, Fine RS, Bradfield JP, Esko T, Giri A, Graff M et al. 2018. Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity. Nat Genet, 50 (1), pp. 26-41. | Show Abstract | Read more

Genome-wide association studies (GWAS) have identified >250 loci for body mass index (BMI), implicating pathways related to neuronal biology. Most GWAS loci represent clusters of common, noncoding variants from which pinpointing causal genes remains challenging. Here we combined data from 718,734 individuals to discover rare and low-frequency (minor allele frequency (MAF) < 5%) coding variants associated with BMI. We identified 14 coding variants in 13 genes, of which 8 variants were in genes (ZBTB7B, ACHE, RAPGEF3, RAB21, ZFHX3, ENTPD6, ZFR2 and ZNF169) newly implicated in human obesity, 2 variants were in genes (MC4R and KSR2) previously observed to be mutated in extreme obesity and 2 variants were in GIPR. The effect sizes of rare variants are ~10 times larger than those of common variants, with the largest effect observed in carriers of an MC4R mutation introducing a stop codon (p.Tyr35Ter, MAF = 0.01%), who weighed ~7 kg more than non-carriers. Pathway analyses based on the variants associated with BMI confirm enrichment of neuronal genes and provide new evidence for adipocyte and energy expenditure biology, widening the potential of genetically supported therapeutic targets in obesity.

Flannick J, Fuchsberger C, Mahajan A, Teslovich TM, Agarwala V, Gaulton KJ, Caulkins L, Koesterer R, Ma C, Moutsianas L et al. 2017. Sequence data and association statistics from 12,940 type 2 diabetes cases and controls. Sci Data, 4 pp. 170179. | Show Abstract | Read more

To investigate the genetic basis of type 2 diabetes (T2D) to high resolution, the GoT2D and T2D-GENES consortia catalogued variation from whole-genome sequencing of 2,657 European individuals and exome sequencing of 12,940 individuals of multiple ancestries. Over 27M SNPs, indels, and structural variants were identified, including 99% of low-frequency (minor allele frequency [MAF] 0.1-5%) non-coding variants in the whole-genome sequenced individuals and 99.7% of low-frequency coding variants in the whole-exome sequenced individuals. Each variant was tested for association with T2D in the sequenced individuals, and, to increase power, most were tested in larger numbers of individuals (>80% of low-frequency coding variants in ~82 K Europeans via the exome chip, and ~90% of low-frequency non-coding variants in ~44 K Europeans via genotype imputation). The variants, genotypes, and association statistics from these analyses provide the largest reference to date of human genetic information relevant to T2D, for use in activities such as T2D-focused genotype imputation, functional characterization of variants or genes, and other novel analyses to detect associations between sequence variation and T2D.

Baker A-M, Huang W, Wang X-MM, Jansen M, Ma X-J, Kim J, Anderson CM, Wu X, Pan L, Su N et al. 2017. Robust RNA-based in situ mutation detection delineates colorectal cancer subclonal evolution. Nat Commun, 8 (1), pp. 1998. | Show Abstract | Read more

Intra-tumor heterogeneity (ITH) is a major underlying cause of therapy resistance and disease recurrence, and is a read-out of tumor growth. Current genetic ITH analysis methods do not preserve spatial context and may not detect rare subclones. Here, we address these shortfalls by developing and validating BaseScope-a novel mutation-specific RNA in situ hybridization assay. We target common point mutations in the BRAF, KRAS and PIK3CA oncogenes in archival colorectal cancer samples to precisely map the spatial and morphological context of mutant subclones. Computational modeling suggests that subclones must arise sufficiently early, or carry a considerable fitness advantage, to form large or spatially disparate subclones. Examples of putative treatment-resistant cells isolated in small topographical areas are observed. The BaseScope assay represents a significant technical advance for in situ mutation detection that provides new insight into tumor evolution, and could have ramifications for selecting patients for treatment.

Rodríguez-Galán A, Salman AM, Bowyer G, Collins KA, Longley RJ, Brod F, Ulaszewska M, Ewer KJ, Janse CJ, Khan SM et al. 2017. An in vitro assay to measure antibody-mediated inhibition of P. berghei sporozoite invasion against P. falciparum antigens. Sci Rep, 7 (1), pp. 17011. | Show Abstract | Read more

A large research effort is currently underway to find an effective and affordable malaria vaccine. Tools that enable the rapid evaluation of protective immune responses are essential to vaccine development as they can provide selection criteria to rank order vaccine candidates. In this study we have revisited the Inhibition of Sporozoite Invasion (ISI) assay to assess the ability of antibodies to inhibit sporozoite infection of hepatocytes. By using GFP expressing sporozoites of the rodent parasite P. berghei we are able to robustly quantify parasite infection of hepatocyte cell lines by flow cytometry. In conjunction with recently produced transgenic P. berghei parasites that express P. falciparum sporozoite antigens, we have been able to use this assay to measure antibody mediated inhibition of sporozoite invasion against one of the lead malaria antigens P. falciparum CSP. By combining chimeric rodent parasites expressing P. falciparum antigens and a flow cytometric readout of infection, we are able to robustly assess vaccine-induced antibodies, from mice, rhesus macaques and human clinical trials, for their functional ability to block sporozoite invasion of hepatocytes.

Ewer K, Sebastian S, Spencer AJ, Gilbert S, Hill AVS, Lambe T. 2017. Chimpanzee adenoviral vectors as vaccines for outbreak pathogens. Hum Vaccin Immunother, 13 (12), pp. 3020-3032. | Show Abstract | Read more

The 2014-15 Ebola outbreak in West Africa highlighted the potential for large disease outbreaks caused by emerging pathogens and has generated considerable focus on preparedness for future epidemics. Here we discuss drivers, strategies and practical considerations for developing vaccines against outbreak pathogens. Chimpanzee adenoviral (ChAd) vectors have been developed as vaccine candidates for multiple infectious diseases and prostate cancer. ChAd vectors are safe and induce antigen-specific cellular and humoral immunity in all age groups, as well as circumventing the problem of pre-existing immunity encountered with human Ad vectors. For these reasons, such viral vectors provide an attractive platform for stockpiling vaccines for emergency deployment in response to a threatened outbreak of an emerging pathogen. Work is already underway to develop vaccines against a number of other outbreak pathogens and we will also review progress on these approaches here, particularly for Lassa fever, Nipah and MERS.

Liu DJ, Peloso GM, Yu H, Butterworth AS, Wang X, Mahajan A, Saleheen D, Emdin C, Alam D, Alves AC et al. 2017. Exome-wide association study of plasma lipids in >300,000 individuals. Nat Genet, 49 (12), pp. 1758-1766. | Citations: 2 (Scopus) | Show Abstract | Read more

We screened variants on an exome-focused genotyping array in >300,000 participants (replication in >280,000 participants) and identified 444 independent variants in 250 loci significantly associated with total cholesterol (TC), high-density-lipoprotein cholesterol (HDL-C), low-density-lipoprotein cholesterol (LDL-C), and/or triglycerides (TG). At two loci (JAK2 and A1CF), experimental analysis in mice showed lipid changes consistent with the human data. We also found that: (i) beta-thalassemia trait carriers displayed lower TC and were protected from coronary artery disease (CAD); (ii) excluding the CETP locus, there was not a predictable relationship between plasma HDL-C and risk for age-related macular degeneration; (iii) only some mechanisms of lowering LDL-C appeared to increase risk for type 2 diabetes (T2D); and (iv) TG-lowering alleles involved in hepatic production of TG-rich lipoproteins (TM6SF2 and PNPLA3) tracked with higher liver fat, higher risk for T2D, and lower risk for CAD, whereas TG-lowering alleles involved in peripheral lipolysis (LPL and ANGPTL4) had no effect on liver fat but decreased risks for both T2D and CAD.

Conti V, Carabalona A, Pallesi-Pocachard E, Leventer RJ, Schaller F, Parrini E, Deparis AA, Watrin F, Buhler E, Novara F et al. 2017. A Novel Strategy Combining Array-CGH, Whole-exome Sequencing and In Utero Electroporation in Rodents to Identify Causative Genes for Brain Malformations. J Vis Exp, 2017 (130), | Show Abstract | Read more

Birth defects that involve the cerebral cortex - also known as malformations of cortical development (MCD) - are important causes of intellectual disability and account for 20-40% of drug-resistant epilepsy in childhood. High-resolution brain imaging has facilitated in vivo identification of a large group of MCD phenotypes. Despite the advances in brain imaging, genomic analysis and generation of animal models, a straightforward workflow to systematically prioritize candidate genes and to test functional effects of putative mutations is missing. To overcome this problem, an experimental strategy enabling the identification of novel causative genes for MCD was developed and validated. This strategy is based on identifying candidate genomic regions or genes via array-CGH or whole-exome sequencing and characterizing the effects of their inactivation or of overexpression of specific mutations in developing rodent brains via in utero electroporation. This approach led to the identification of the C6orf70 gene, encoding for a putative vesicular protein, to the pathogenesis of periventricular nodular heterotopia, a MCD caused by defective neuronal migration.

Anopheles gambiae 1000 Genomes Consortium, Data analysis group, Partner working group, Sample collections—Angola:, Burkina Faso:, Cameroon:, Gabon:, Guinea:, Guinea-Bissau:, Kenya: et al. 2017. Genetic diversity of the African malaria vector Anopheles gambiae. Nature, 552 (7683), pp. 96-100. | Citations: 1 (Web of Science Lite) | Show Abstract | Read more

The sustainability of malaria control in Africa is threatened by the rise of insecticide resistance in Anopheles mosquitoes, which transmit the disease. To gain a deeper understanding of how mosquito populations are evolving, here we sequenced the genomes of 765 specimens of Anopheles gambiae and Anopheles coluzzii sampled from 15 locations across Africa, and identified over 50 million single nucleotide polymorphisms within the accessible genome. These data revealed complex population structure and patterns of gene flow, with evidence of ancient expansions, recent bottlenecks, and local variation in effective population size. Strong signals of recent selection were observed in insecticide-resistance genes, with several sweeps spreading over large geographical distances and between species. The design of new tools for mosquito control using gene-drive systems will need to take account of high levels of genetic diversity in natural mosquito populations.

Wang M, Quinn CM, Perilla JR, Zhang H, Shirra R, Hou G, Byeon I-J, Suiter CL, Ablan S, Urano E et al. 2017. Quenching protein dynamics interferes with HIV capsid maturation. Nat Commun, 8 (1), pp. 1779. | Show Abstract | Read more

Maturation of HIV-1 particles encompasses a complex morphological transformation of Gag via an orchestrated series of proteolytic cleavage events. A longstanding question concerns the structure of the C-terminal region of CA and the peptide SP1 (CA-SP1), which represents an intermediate during maturation of the HIV-1 virus. By integrating NMR, cryo-EM, and molecular dynamics simulations, we show that in CA-SP1 tubes assembled in vitro, which represent the features of an intermediate assembly state during maturation, the SP1 peptide exists in a dynamic helix-coil equilibrium, and that the addition of the maturation inhibitors Bevirimat and DFH-055 causes stabilization of a helical form of SP1. Moreover, the maturation-arresting SP1 mutation T8I also induces helical structure in SP1 and further global dynamical and conformational changes in CA. Overall, our results show that dynamics of CA and SP1 are critical for orderly HIV-1 maturation and that small molecules can inhibit maturation by perturbing molecular motions.

Burns A, Alsolami R, Becq J, Stamatopoulos B, Timbs A, Bruce D, Robbe P, Vavoulis D, Clifford R, Cabes M et al. 2018. Whole-genome sequencing of chronic lymphocytic leukaemia reveals distinct differences in the mutational landscape between IgHVmut and IgHVunmut subgroups. Leukemia, 32 (2), pp. 573. | Show Abstract | Read more

This corrects the article DOI: 10.1038/leu.2017.177.

Al-Mossawi MH, Chen L, Fang H, Ridley A, de Wit J, Yager N, Hammitzsch A, Pulyakhina I, Fairfax BP, Simone D et al. 2017. Unique transcriptome signatures and GM-CSF expression in lymphocytes from patients with spondyloarthritis. Nat Commun, 8 (1), pp. 1510. | Show Abstract | Read more

Spondyloarthritis encompasses a group of common inflammatory diseases thought to be driven by IL-17A-secreting type-17 lymphocytes. Here we show increased numbers of GM-CSF-producing CD4 and CD8 lymphocytes in the blood and joints of patients with spondyloarthritis, and increased numbers of IL-17A+GM-CSF+ double-producing CD4, CD8, γδ and NK cells. GM-CSF production in CD4 T cells occurs both independently and in combination with classical Th1 and Th17 cytokines. Type 3 innate lymphoid cells producing predominantly GM-CSF are expanded in synovial tissues from patients with spondyloarthritis. GM-CSF+CD4+ cells, isolated using a triple cytokine capture approach, have a specific transcriptional signature. Both GM-CSF+ and IL-17A+GM-CSF+ double-producing CD4 T cells express increased levels of GPR65, a proton-sensing receptor associated with spondyloarthritis in genome-wide association studies and pathogenicity in murine inflammatory disease models. Silencing GPR65 in primary CD4 T cells reduces GM-CSF production. GM-CSF and GPR65 may thus serve as targets for therapeutic intervention of spondyloarthritis.

Latva-Rasku A, Honka M-J, Stančáková A, Koistinen HA, Kuusisto J, Guan L, Manning AK, Stringham H, Gloyn AL, Lindgren CM et al. 2018. A Partial Loss-of-Function Variant in AKT2 Is Associated With Reduced Insulin-Mediated Glucose Uptake in Multiple Insulin-Sensitive Tissues: A Genotype-Based Callback Positron Emission Tomography Study. Diabetes, 67 (2), pp. 334-342. | Show Abstract | Read more

Rare fully penetrant mutations in AKT2 are an established cause of monogenic disorders of glucose metabolism. Recently, a novel partial loss-of-function AKT2 coding variant (p.Pro50Thr) was identified that is nearly specific to Finns (frequency 1.1%), with the low-frequency allele associated with an increase in fasting plasma insulin level and risk of type 2 diabetes. The effects of the p.Pro50Thr AKT2 variant (p.P50T/AKT2) on insulin-stimulated glucose uptake (GU) in the whole body and in different tissues have not previously been investigated. We identified carriers (N = 20) and matched noncarriers (N = 25) for this allele in the population-based Metabolic Syndrome in Men (METSIM)study and invited these individuals back for positron emission tomography study with [18F]-fluorodeoxyglucose during euglycemic hyperinsulinemia. When we compared p.P50T/AKT2 carriers to noncarriers, we found a 39.4% reduction in whole-body GU (P = 0.006) and a 55.6% increase in the rate of endogenous glucose production (P = 0.038). We found significant reductions in GU in multiple tissues-skeletal muscle (36.4%), liver (16.1%), brown adipose (29.7%), and bone marrow (32.9%)-and increases of 16.8-19.1% in seven tested brain regions. These data demonstrate that the p.P50T substitution of AKT2 influences insulin-mediated GU in multiple insulin-sensitive tissues and may explain, at least in part, the increased risk of type 2 diabetes in p.P50T/AKT2 carriers.

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

A correction to this article has been published and is linked from the HTML version of this article.

Cassidy CK, Himes BA, Luthey-Schulten Z, Zhang P. 2017. CryoEM-based hybrid modeling approaches for structure determination. Curr Opin Microbiol, 43 pp. 14-23. | Show Abstract | Read more

Recent advances in cryo-electron microscopy (cryoEM) have dramatically improved the resolutions at which vitrified biological specimens can be studied, revealing new structural and mechanistic insights over a broad range of spatial scales. Bolstered by these advances, much effort has been directed toward the development of hybrid modeling methodologies for the construction and refinement of high-fidelity atomistic models from cryoEM data. In this brief review, we will survey the key elements of cryoEM-based hybrid modeling, providing an overview of available computational tools and strategies as well as several recent applications.

Payne RO, Silk SE, Elias SC, Miura K, Diouf A, Galaway F, de Graaf H, Brendish NJ, Poulton ID, Griffiths OJ et al. 2017. Human vaccination against RH5 induces neutralizing antimalarial antibodies that inhibit RH5 invasion complex interactions. JCI Insight, 2 (21), | Show Abstract | Read more

The development of a highly effective vaccine remains a key strategic goal to aid the control and eventual eradication of Plasmodium falciparum malaria. In recent years, the reticulocyte-binding protein homolog 5 (RH5) has emerged as the most promising blood-stage P. falciparum candidate antigen to date, capable of conferring protection against stringent challenge in Aotus monkeys. We report on the first clinical trial to our knowledge to assess the RH5 antigen - a dose-escalation phase Ia study in 24 healthy, malaria-naive adult volunteers. We utilized established viral vectors, the replication-deficient chimpanzee adenovirus serotype 63 (ChAd63), and the attenuated orthopoxvirus modified vaccinia virus Ankara (MVA), encoding RH5 from the 3D7 clone of P. falciparum. Vaccines were administered i.m. in a heterologous prime-boost regimen using an 8-week interval and were well tolerated. Vaccine-induced anti-RH5 serum antibodies exhibited cross-strain functional growth inhibition activity (GIA) in vitro, targeted linear and conformational epitopes within RH5, and inhibited key interactions within the RH5 invasion complex. This is the first time to our knowledge that substantial RH5-specific responses have been induced by immunization in humans, with levels greatly exceeding the serum antibody responses observed in African adults following years of natural malaria exposure. These data support the progression of RH5-based vaccines to human efficacy testing.

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, 66 (11), pp. 2888-2902. | Citations: 6 (Scopus) | Show Abstract | Read more

To characterize 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 case and 132,532 control subjects of European ancestry after imputation using the 1000 Genomes multiethnic reference panel. Promising association signals were followed up in additional data sets (of 14,545 or 7,397 T2D case and 38,994 or 71,604 control subjects). We identified 13 novel T2D-associated loci (P< 5 × 10-8), including variants near theGLP2R,GIP, andHLA-DQA1genes. 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 single nucleotide variants. 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.

Apinjoh TO, Mugri RN, Miotto O, Chi HF, Tata RB, Anchang-Kimbi JK, Fon EM, Tangoh DA, Nyingchu RV, Jacob C et al. 2017. Molecular markers for artemisinin and partner drug resistance in natural Plasmodium falciparum populations following increased insecticide treated net coverage along the slope of mount Cameroon: cross-sectional study. Infect Dis Poverty, 6 (1), pp. 136. | Show Abstract | Read more

BACKGROUND: Drug resistance is one of the greatest challenges of malaria control programmes, with the monitoring of parasite resistance to artemisinins or to Artemisinin Combination Therapy (ACT) partner drugs critical to elimination efforts. Markers of resistance to a wide panel of antimalarials were assessed in natural parasite populations from southwestern Cameroon. METHODS: Individuals with asymptomatic parasitaemia or uncomplicated malaria were enrolled through cross-sectional surveys from May 2013 to March 2014 along the slope of mount Cameroon. Plasmodium falciparum malaria parasitaemic blood, screened by light microscopy, was depleted of leucocytes using CF11 cellulose columns and the parasite genotype ascertained by sequencing on the Illumina HiSeq platform. RESULTS: A total of 259 participants were enrolled in this study from three different altitudes. While some alleles associated with drug resistance in pfdhfr, pfmdr1 and pfcrt were highly prevalent, less than 3% of all samples carried mutations in the pfkelch13 gene, none of which were amongst those associated with slow artemisinin parasite clearance rates in Southeast Asia. The most prevalent haplotypes were triple mutants Pfdhfr I51R59N108I164(99%), pfcrt- C72V73I74E75T76(47.3%), and single mutants PfdhpsS436G437K540A581A613(69%) and Pfmdr1 N86F184D1246(53.2%). CONCLUSIONS: The predominance of the Pf pfcrt CVIET and Pf dhfr IRN triple mutant parasites and absence of pfkelch13 resistance alleles suggest that the amodiaquine and pyrimethamine components of AS-AQ and SP may no longer be effective in their role while chloroquine resistance still persists in southwestern Cameroon.

Matalliotakis M, Zervou MI, Matalliotaki C, Rahmioglu N, Koumantakis G, Kalogiannidis I, Prapas I, Zondervan K, Spandidos DA, Matalliotakis I, Goulielmos GN. 2017. The role of gene polymorphisms in endometriosis. Mol Med Rep, 16 (5), pp. 5881-5886. | Citations: 1 (Scopus) | Show Abstract | Read more

Endometriosis is a benign gynecologic disorder, affecting up to 10% of women, characterized by the presence of functional endometrial tissue at ectopic positions generally within the peritoneum. It is a heritable condition influenced by multiple genetic and environmental factors, with an overall heritability estimated at approximately 50%. In this study, we investigated whether single nucleotide polymorphisms (SNPs) rs7521902, rs10859871 and rs11031006, mapping to WNT4, VEZT and FSHB genetic loci, respectively, are associated with risk for endometriosis in a Greek population. This study included 166 women with histologically confirmed endometriosis diagnosed through surgery and 150 normal controls. Genotyping of the rs7521902, rs10859871 and rs11031006 SNPs was performed with Taqman primer/probe sets. A significant association was detected with the AC genotype of rs7521902 (WNT4) in patients with stage III and IV disease only. Evidence for association with endometriosis was also found for the AC genotype of the rs10859871 of VEZT. Notably, a significant difference in the distribution of the AG genotype and the minor allele A of FSHB rs11031006 SNP was found between the endometriosis patients and controls. We found a genetic association between rs11031006 (FSHB) SNP and endometriosis. WNT4 and VEZT genes constitute the most consistently associated genes with endometriosis. In the present study, an association of rs7521902 (WNT4) and rs10859871 (VEZT) was confirmed in women with endometriosis at the genotypic but not the allelic level.

Strawbridge RJ, Silveira A, Hoed MD, Gustafsson S, Luan J, Rybin D, Dupuis J, Li-Gao R, Kavousi M, Dehghan A et al. 2017. Identification of a novel proinsulin-associated SNP and demonstration that proinsulin is unlikely to be a causal factor in subclinical vascular remodelling using Mendelian randomisation. Atherosclerosis, 266 pp. 196-204. | Show Abstract | Read more

BACKGROUND AND AIMS: Increased proinsulin relative to insulin levels have been associated with subclinical atherosclerosis (measured by carotid intima-media thickness (cIMT)) and are predictive of future cardiovascular disease (CVD), independently of established risk factors. The mechanisms linking proinsulin to atherosclerosis and CVD are unclear. A genome-wide meta-analysis has identified nine loci associated with circulating proinsulin levels. Using proinsulin-associated SNPs, we set out to use a Mendelian randomisation approach to test the hypothesis that proinsulin plays a causal role in subclinical vascular remodelling. METHODS: We studied the high CVD-risk IMPROVE cohort (n = 3345), which has detailed biochemical phenotyping and repeated, state-of-the-art, high-resolution carotid ultrasound examinations. Genotyping was performed using Illumina Cardio-Metabo and Immuno arrays, which include reported proinsulin-associated loci. Participants with type 2 diabetes (n = 904) were omitted from the analysis. Linear regression was used to identify proinsulin-associated genetic variants. RESULTS: We identified a proinsulin locus on chromosome 15 (rs8029765) and replicated it in data from 20,003 additional individuals. An 11-SNP score, including the previously identified and the chromosome 15 proinsulin-associated loci, was significantly and negatively associated with baseline IMTmean and IMTmax (the primary cIMT phenotypes) but not with progression measures. However, MR-Eggers refuted any significant effect of the proinsulin-associated 11-SNP score, and a non-pleiotropic SNP score of three variants (including rs8029765) demonstrated no effect on baseline or progression cIMT measures. CONCLUSIONS: We identified a novel proinsulin-associated locus and demonstrated that whilst proinsulin levels are associated with cIMT measures, proinsulin per se is unlikely to have a causative effect on cIMT.

Thurner M, Shenhav L, Wesolowska-Andersen A, Bennett AJ, Barrett A, Gloyn AL, McCarthy MI, Beer NL, Efrat S. 2017. Genes Associated with Pancreas Development and Function Maintain Open Chromatin in iPSCs Generated from Human Pancreatic Beta Cells. Stem Cell Reports, 9 (5), pp. 1395-1405. | Show Abstract | Read more

Current in vitro islet differentiation protocols suffer from heterogeneity and low efficiency. Induced pluripotent stem cells (iPSCs) derived from pancreatic beta cells (BiPSCs) preferentially differentiate toward endocrine pancreas-like cells versus those from fibroblasts (FiPSCs). We interrogated genome-wide open chromatin in BiPSCs and FiPSCs via ATAC-seq and identified ∼8.3k significant, differential open chromatin sites (DOCS) between the two iPSC subtypes (false discovery rate [FDR] < 0.05). DOCS where chromatin was more accessible in BiPSCs (Bi-DOCS) were significantly enriched for known regulators of endodermal development, including bivalent and weak enhancers, and FOXA2 binding sites (FDR < 0.05). Bi-DOCS were associated with genes related to pancreas development and beta-cell function, including transcription factors mutated in monogenic diabetes (PDX1, NKX2-2, HNF1A; FDR < 0.05). Moreover, Bi-DOCS correlated with enhanced gene expression in BiPSC-derived definitive endoderm and pancreatic progenitor cells. Bi-DOCS therefore highlight genes and pathways governing islet-lineage commitment, which can be exploited for differentiation protocol optimization, diabetes disease modeling, and therapeutic purposes.

Cheng J, Demeulemeester J, Wedge DC, Vollan HKM, Pitt JJ, Russnes HG, Pandey BP, Nilsen G, Nord S, Bignell GR et al. 2017. Pan-cancer analysis of homozygous deletions in primary tumours uncovers rare tumour suppressors. Nat Commun, 8 (1), pp. 1221. | Show Abstract | Read more

Homozygous deletions are rare in cancers and often target tumour suppressor genes. Here, we build a compendium of 2218 primary tumours across 12 human cancer types and systematically screen for homozygous deletions, aiming to identify rare tumour suppressors. Our analysis defines 96 genomic regions recurrently targeted by homozygous deletions. These recurrent homozygous deletions occur either over tumour suppressors or over fragile sites, regions of increased genomic instability. We construct a statistical model that separates fragile sites from regions showing signatures of positive selection for homozygous deletions and identify candidate tumour suppressors within those regions. We find 16 established tumour suppressors and propose 27 candidate tumour suppressors. Several of these genes (including MGMT, RAD17, and USP44) show prior evidence of a tumour suppressive function. Other candidate tumour suppressors, such as MAFTRR, KIAA1551, and IGF2BP2, are novel. Our study demonstrates how rare tumour suppressors can be identified through copy number meta-analysis.

Altemose N, Noor N, Bitoun E, Tumian A, Imbeault M, Chapman JR, Aricescu AR, Myers SR. 2017. A map of human PRDM9 binding provides evidence for novel behaviors of PRDM9 and other zinc-finger proteins in meiosis. Elife, 6 | Show Abstract | Read more

PRDM9 binding localizes almost all meiotic recombination sites in humans and mice. However, most PRDM9-bound loci do not become recombination hotspots. To explore factors that affect binding and subsequent recombination outcomes, we mapped human PRDM9 binding sites in a transfected human cell line and measured PRDM9-induced histone modifications. These data reveal varied DNA-binding modalities of PRDM9. We also find that human PRDM9 frequently binds promoters, despite their low recombination rates, and it can activate expression of a small number of genes including CTCFL and VCX. Furthermore, we identify specific sequence motifs that predict consistent, localized meiotic recombination suppression around a subset of PRDM9 binding sites. These motifs strongly associate with KRAB-ZNF protein binding, TRIM28 recruitment, and specific histone modifications. Finally, we demonstrate that, in addition to binding DNA, PRDM9's zinc fingers also mediate its multimerization, and we show that a pair of highly diverged alleles preferentially form homo-multimers.

Ormondroyd E, Mackley MP, Blair E, Craft J, Knight JC, Taylor JC, Taylor J, Watkins H. 2017. "Not pathogenic until proven otherwise": perspectives of UK clinical genomics professionals toward secondary findings in context of a Genomic Medicine Multidisciplinary Team and the 100,000 Genomes Project. Genet Med, | Show Abstract | Read more

PurposeApproaches to secondary findings in genome sequencing (GS) are unresolved. In the United Kingdom, GS is now routinely available through the 100,000 Genomes Project, which offers participants feedback of limited secondary findings.MethodsIn Oxford, a Genomic Medicine Multidisciplinary Team (GM-MDT) governs local access to GS, and reviews findings. Semistructured interviews were conducted with 19 GM-MDT members to explore perspectives on secondary findings.ResultsWhile enthusiastic about GS for diagnosing rare disease, members question the rationale for genome screening largely because of lack of evidence for clinical utility and limited justification for use of resources. Members' views are drawn from diverse experiences; they feel a strong sense of responsibility to act in participants' best interests. The capacity to return limited secondary findings should be enabled, but members favor a cautious approach that is responsive to accumulating evidence. Informed participant choice is considered critical, yet challenging. Discrimination of variants is considered essential, and requiring of specialist input and consensus. Multiple areas requiring enhanced engagement and education are identified, i.e., for patients, the public, and health-care professionals; at present, mainstreaming of genomics may be premature.ConclusionUK experts believe that evidence to inform policy toward secondary findings is lacking, arguing for caution.Genetics in medicine advance online publication, 26 October 2017; doi:10.1038/gim.2017.157.

Davies B, Brown LA, Cais O, Watson J, Clayton AJ, Chang VT, Biggs D, Preece C, Hernandez-Pliego P, Krohn J et al. 2017. A point mutation in the ion conduction pore of AMPA receptor GRIA3 causes dramatically perturbed sleep patterns as well as intellectual disability. Hum Mol Genet, 26 (20), pp. 3869-3882. | Show Abstract | Read more

The discovery of genetic variants influencing sleep patterns can shed light on the physiological processes underlying sleep. As part of a large clinical sequencing project, WGS500, we sequenced a family in which the two male children had severe developmental delay and a dramatically disturbed sleep-wake cycle, with very long wake and sleep durations, reaching up to 106-h awake and 48-h asleep. The most likely causal variant identified was a novel missense variant in the X-linked GRIA3 gene, which has been implicated in intellectual disability. GRIA3 encodes GluA3, a subunit of AMPA-type ionotropic glutamate receptors (AMPARs). The mutation (A653T) falls within the highly conserved transmembrane domain of the ion channel gate, immediately adjacent to the analogous residue in the Grid2 (glutamate receptor) gene, which is mutated in the mouse neurobehavioral mutant, Lurcher. In vitro, the GRIA3(A653T) mutation stabilizes the channel in a closed conformation, in contrast to Lurcher. We introduced the orthologous mutation into a mouse strain by CRISPR-Cas9 mutagenesis and found that hemizygous mutants displayed significant differences in the structure of their activity and sleep compared to wild-type littermates. Typically, mice are polyphasic, exhibiting multiple sleep bouts of sleep several minutes long within a 24-h period. The Gria3A653T mouse showed significantly fewer brief bouts of activity and sleep than the wild-types. Furthermore, Gria3A653T mice showed enhanced period lengthening under constant light compared to wild-type mice, suggesting an increased sensitivity to light. Our results suggest a role for GluA3 channel activity in the regulation of sleep behavior in both mice and humans.

Devoy A, Kalmar B, Stewart M, Park H, Burke B, Noy SJ, Redhead Y, Humphrey J, Lo K, Jaeger J et al. 2017. Humanized mutant FUS drives progressive motor neuron degeneration without aggregation in 'FUSDelta14' knockin mice. Brain, 140 (11), pp. 2797-2805. | Show Abstract | Read more

Mutations in FUS are causative for amyotrophic lateral sclerosis with a dominant mode of inheritance. In trying to model FUS-amyotrophic lateral sclerosis (ALS) in mouse it is clear that FUS is dosage-sensitive and effects arise from overexpression per se in transgenic strains. Novel models are required that maintain physiological levels of FUS expression and that recapitulate the human disease-with progressive loss of motor neurons in heterozygous animals. Here, we describe a new humanized FUS-ALS mouse with a frameshift mutation, which fulfils both criteria: the FUS Delta14 mouse. Heterozygous animals express mutant humanized FUS protein at physiological levels and have adult onset progressive motor neuron loss and denervation of neuromuscular junctions. Additionally, we generated a novel antibody to the unique human frameshift peptide epitope, allowing specific identification of mutant FUS only. Using our new FUSDelta14 ALS mouse-antibody system we show that neurodegeneration occurs in the absence of FUS protein aggregation. FUS mislocalization increases as disease progresses, and mutant FUS accumulates at the rough endoplasmic reticulum. Further, transcriptomic analyses show progressive changes in ribosomal protein levels and mitochondrial function as early disease stages are initiated. Thus, our new physiological mouse model has provided novel insight into the early pathogenesis of FUS-ALS.

Inshaw JRJ, Walker NM, Wallace C, Bottolo L, Todd JA. 2018. The chromosome 6q22.33 region is associated with age at diagnosis of type 1 diabetes and disease risk in those diagnosed under 5 years of age. Diabetologia, 61 (1), pp. 147-157. | Show Abstract | Read more

AIMS/HYPOTHESIS: The genetic risk of type 1 diabetes has been extensively studied. However, the genetic determinants of age at diagnosis (AAD) of type 1 diabetes remain relatively unexplained. Identification of AAD genes and pathways could provide insight into the earliest events in the disease process. METHODS: Using ImmunoChip data from 15,696 cases, we aimed to identify regions in the genome associated with AAD. RESULTS: Two regions were convincingly associated with AAD (p < 5 × 10-8): the MHC on 6p21, and 6q22.33. Fine-mapping of 6q22.33 identified two AAD-associated haplotypes in the region nearest to the genes encoding protein tyrosine phosphatase receptor kappa (PTPRK) and thymocyte-expressed molecule involved in selection (THEMIS). We examined the susceptibility to type 1 diabetes at these SNPs by performing a meta-analysis including 19,510 control participants. Although these SNPs were not associated with type 1 diabetes overall (p > 0.001), the SNP most associated with AAD, rs72975913, was associated with susceptibility to type 1 diabetes in those individuals diagnosed at less than 5 years old (p = 2.3 × 10-9). CONCLUSION/INTERPRETATION: PTPRK and its neighbour THEMIS are required for early development of the thymus, which we can assume influences the initiation of autoimmunity. Non-HLA genes may only be detectable as risk factors for the disease in individuals diagnosed under the age 5 years because, after that period of immune development, their role in disease susceptibility has become redundant.

Vavoulis DV, Taylor JC, Schuh A. 2017. Hierarchical probabilistic models for multiple gene/variant associations based on next-generation sequencing data. Bioinformatics, 33 (19), pp. 3058-3064. | 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 versus count-based (particularly Negative-Binomial-based) models for eQTL mapping. Availability and implementation: All methods are implemented in the free software eQTLseq: https://github.com/dvav/eQTLseq. Contact: dimitris.vavoulis@well.ox.ac.uk. Supplementary information: Supplementary data are available at Bioinformatics online.

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, 5 (10), pp. 816-826. | Citations: 4 (Scopus) | 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.

Kraja AT, Cook JP, Warren HR, Surendran P, Liu C, Evangelou E, Manning AK, Grarup N, Drenos F, Sim X et al. 2017. New Blood Pressure-Associated Loci Identified in Meta-Analyses of 475 000 Individuals. Circ Cardiovasc Genet, 10 (5), pp. e001778-e001778. | Citations: 1 (Web of Science Lite) | Show Abstract | Read more

BACKGROUND: Genome-wide association studies have recently identified >400 loci that harbor DNA sequence variants that influence blood pressure (BP). Our earlier studies identified and validated 56 single nucleotide variants (SNVs) associated with BP from meta-analyses of exome chip genotype data. An additional 100 variants yielded suggestive evidence of association. METHODS AND RESULTS: Here, we augment the sample with 140 886 European individuals from the UK Biobank, in whom 77 of the 100 suggestive SNVs were available for association analysis with systolic BP or diastolic BP or pulse pressure. We performed 2 meta-analyses, one in individuals of European, South Asian, African, and Hispanic descent (pan-ancestry, ≈475 000), and the other in the subset of individuals of European descent (≈423 000). Twenty-one SNVs were genome-wide significant (P<5×10-8) for BP, of which 4 are new BP loci: rs9678851 (missense, SLC4A1AP), rs7437940 (AFAP1), rs13303 (missense, STAB1), and rs1055144 (7p15.2). In addition, we identified a potentially independent novel BP-associated SNV, rs3416322 (missense, SYNPO2L) at a known locus, uncorrelated with the previously reported SNVs. Two SNVs are associated with expression levels of nearby genes, and SNVs at 3 loci are associated with other traits. One SNV with a minor allele frequency <0.01, (rs3025380 at DBH) was genome-wide significant. CONCLUSIONS: We report 4 novel loci associated with BP regulation, and 1 independent variant at an established BP locus. This analysis highlights several candidate genes with variation that alter protein function or gene expression for potential follow-up.

Wijesurendra RS, Liu A, Ferreira VM, Neubauer S, Casadei B. 2017. Histopathological and Immunological Characteristics of Tachycardia-Induced Cardiomyopathy. J Am Coll Cardiol, 70 (13), pp. 1686. | Read more

Camacho N, Van Loo P, Edwards S, Kay JD, Matthews L, Haase K, Clark J, Dennis N, Thomas S, Kremeyer B et al. 2017. Appraising the relevance of DNA copy number loss and gain in prostate cancer using whole genome DNA sequence data. PLoS Genet, 13 (9), pp. e1007001. | Show Abstract | Read more

A variety of models have been proposed to explain regions of recurrent somatic copy number alteration (SCNA) in human cancer. Our study employs Whole Genome DNA Sequence (WGS) data from tumor samples (n = 103) to comprehensively assess the role of the Knudson two hit genetic model in SCNA generation in prostate cancer. 64 recurrent regions of loss and gain were detected, of which 28 were novel, including regions of loss with more than 15% frequency at Chr4p15.2-p15.1 (15.53%), Chr6q27 (16.50%) and Chr18q12.3 (17.48%). Comprehensive mutation screens of genes, lincRNA encoding sequences, control regions and conserved domains within SCNAs demonstrated that a two-hit genetic model was supported in only a minor proportion of recurrent SCNA losses examined (15/40). We found that recurrent breakpoints and regions of inversion often occur within Knudson model SCNAs, leading to the identification of ZNF292 as a target gene for the deletion at 6q14.3-q15 and NKX3.1 as a two-hit target at 8p21.3-p21.2. The importance of alterations of lincRNA sequences was illustrated by the identification of a novel mutational hotspot at the KCCAT42, FENDRR, CAT1886 and STCAT2 loci at the 16q23.1-q24.3 loss. Our data confirm that the burden of SCNAs is predictive of biochemical recurrence, define nine individual regions that are associated with relapse, and highlight the possible importance of ion channel and G-protein coupled-receptor (GPCR) pathways in cancer development. We concluded that a two-hit genetic model accounts for about one third of SCNA indicating that mechanisms, such haploinsufficiency and epigenetic inactivation, account for the remaining SCNA losses.

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, pp. 1-1. | 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.

Dendrou CA, Fugger L. 2017. Immunomodulation in multiple sclerosis: promises and pitfalls. Curr Opin Immunol, 49 pp. 37-43. | Show Abstract | Read more

Multiple sclerosis (MS) afflicts about 2.5 million people globally and poses a major personal and socioeconomic burden. The recognition of MS as an inflammatory disease, characterized by infiltration of immune cells into the central nervous system, has spurred research into the autoimmune etiology of the condition and has provided the rationale for its treatment through immunomodulation. Experience with immunotherapies in MS to date has suggested a disparity between the observed immune cell infiltration and the progressive loss of neurons. However, recent clinical efforts are providing new insights into progressive MS that once again place the immune system at center stage. This article reviews the main mechanisms of MS immunopathogenesis, and the benefits, risks and challenges of immunomodulatory treatments for the disease.

Alvarez FJD, He S, Perilla JR, Jang S, Schulten K, Engelman AN, Scheres SHW, Zhang P. 2017. CryoEM structure of MxB reveals a novel oligomerization interface critical for HIV restriction. Sci Adv, 3 (9), pp. e1701264. | Show Abstract | Read more

Human dynamin-like, interferon-induced myxovirus resistance 2 (Mx2 or MxB) is a potent HIV-1 inhibitor. Antiviral activity requires both the amino-terminal region of MxB and protein oligomerization, each of which has eluded structural determination due to difficulties in protein preparation. We report that maltose binding protein-fused, full-length wild-type MxB purifies as oligomers and further self-assembles into helical arrays in physiological salt. Guanosine triphosphate (GTP), but not guanosine diphosphate, binding results in array disassembly, whereas subsequent GTP hydrolysis allows its reformation. Using cryo-electron microscopy (cryoEM), we determined the MxB assembly structure at 4.6 Å resolution, representing the first near-atomic resolution structure in the mammalian dynamin superfamily. The structure revealed previously described and novel MxB assembly interfaces. Mutational analyses demonstrated a critical role for one of the novel interfaces in HIV-1 restriction.

Wheeler E, Leong A, Liu C-T, Hivert M-F, Strawbridge RJ, Podmore C, Li M, Yao J, Sim X, Hong J et al. 2017. Impact of common genetic determinants of Hemoglobin A1c on type 2 diabetes risk and diagnosis in ancestrally diverse populations: A transethnic genome-wide meta-analysis. PLoS Med, 14 (9), pp. e1002383. | Citations: 5 (Scopus) | Show Abstract | Read more

BACKGROUND: Glycated hemoglobin (HbA1c) is used to diagnose type 2 diabetes (T2D) and assess glycemic control in patients with diabetes. Previous genome-wide association studies (GWAS) have identified 18 HbA1c-associated genetic variants. These variants proved to be classifiable by their likely biological action as erythrocytic (also associated with erythrocyte traits) or glycemic (associated with other glucose-related traits). In this study, we tested the hypotheses that, in a very large scale GWAS, we would identify more genetic variants associated with HbA1c and that HbA1c variants implicated in erythrocytic biology would affect the diagnostic accuracy of HbA1c. We therefore expanded the number of HbA1c-associated loci and tested the effect of genetic risk-scores comprised of erythrocytic or glycemic variants on incident diabetes prediction and on prevalent diabetes screening performance. Throughout this multiancestry study, we kept a focus on interancestry differences in HbA1c genetics performance that might influence race-ancestry differences in health outcomes. METHODS & FINDINGS: Using genome-wide association meta-analyses in up to 159,940 individuals from 82 cohorts of European, African, East Asian, and South Asian ancestry, we identified 60 common genetic variants associated with HbA1c. We classified variants as implicated in glycemic, erythrocytic, or unclassified biology and tested whether additive genetic scores of erythrocytic variants (GS-E) or glycemic variants (GS-G) were associated with higher T2D incidence in multiethnic longitudinal cohorts (N = 33,241). Nineteen glycemic and 22 erythrocytic variants were associated with HbA1c at genome-wide significance. GS-G was associated with higher T2D risk (incidence OR = 1.05, 95% CI 1.04-1.06, per HbA1c-raising allele, p = 3 × 10-29); whereas GS-E was not (OR = 1.00, 95% CI 0.99-1.01, p = 0.60). In Europeans and Asians, erythrocytic variants in aggregate had only modest effects on the diagnostic accuracy of HbA1c. Yet, in African Americans, the X-linked G6PD G202A variant (T-allele frequency 11%) was associated with an absolute decrease in HbA1c of 0.81%-units (95% CI 0.66-0.96) per allele in hemizygous men, and 0.68%-units (95% CI 0.38-0.97) in homozygous women. The G6PD variant may cause approximately 2% (N = 0.65 million, 95% CI 0.55-0.74) of African American adults with T2D to remain undiagnosed when screened with HbA1c. Limitations include the smaller sample sizes for non-European ancestries and the inability to classify approximately one-third of the variants. Further studies in large multiethnic cohorts with HbA1c, glycemic, and erythrocytic traits are required to better determine the biological action of the unclassified variants. CONCLUSIONS: As G6PD deficiency can be clinically silent until illness strikes, we recommend investigation of the possible benefits of screening for the G6PD genotype along with using HbA1c to diagnose T2D in populations of African ancestry or groups where G6PD deficiency is common. Screening with direct glucose measurements, or genetically-informed HbA1c diagnostic thresholds in people with G6PD deficiency, may be required to avoid missed or delayed diagnoses.

Burren OS, Rubio García A, Javierre B-M, 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.

Cutts A, Venn O, Dilthey A, Gupta A, Vavoulis D, Dreau H, Middleton M, McVean G, Taylor JC, Schuh A. 2017. Characterisation of the changing genomic landscape of metastatic melanoma using cell free DNA. NPJ Genom Med, 2 (1), pp. 25. | Show Abstract | Read more

Cancer is characterised by complex somatically acquired genetic aberrations that manifest as intra-tumour and inter-tumour genetic heterogeneity and can lead to treatment resistance. In this case study, we characterise the genome-wide somatic mutation dynamics in a metastatic melanoma patient during therapy using low-input (50 ng) PCR-free whole genome sequencing of cell-free DNA from pre-treatment and post-relapse blood samples. We identify de novo tumour-specific somatic mutations from cell-free DNA, while the sequence context of single nucleotide variants showed the characteristic UV-damage mutation signature of melanoma. To investigate the behaviour of individual somatic mutations during proto-oncogene B-Raf -targeted and immune checkpoint inhibition, amplicon-based deep sequencing was used to verify and track frequencies of 212 single nucleotide variants at 10 distinct time points over 13 months of treatment. Under checkpoint inhibition therapy, we observed an increase in mutant allele frequencies indicating progression on therapy 88 days before clinical determination of non-response positron emission tomogrophy-computed tomography. We also revealed mutations from whole genome sequencing of cell-free DNA that were not present in the tissue biopsy, but that later contributed to relapse. Our findings have potential clinical applications where high quality tumour-tissue derived DNA is not available.

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. | Citations: 1 (Scopus) | 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.

Grotz AK, Gloyn AL, Thomsen SK. 2017. Prioritising Causal Genes at Type 2 Diabetes Risk Loci. Curr Diab Rep, 17 (9), pp. 76. | Citations: 1 (Web of Science Lite) | 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.

Mägi R, Horikoshi M, Sofer T, Mahajan A, Kitajima H, Franceschini N, McCarthy MI, COGENT-Kidney Consortium, T2D-GENES Consortium, Morris AP. 2017. Trans-ethnic meta-regression of genome-wide association studies accounting for ancestry increases power for discovery and improves fine-mapping resolution. Hum Mol Genet, 26 (18), pp. 3639-3650. | Show Abstract | Read more

Trans-ethnic meta-analysis of genome-wide association studies (GWAS) across diverse populations can increase power to detect complex trait loci when the underlying causal variants are shared between ancestry groups. However, heterogeneity in allelic effects between GWAS at these loci can occur that is correlated with ancestry. Here, a novel approach is presented to detect SNP association and quantify the extent of heterogeneity in allelic effects that is correlated with ancestry. We employ trans-ethnic meta-regression to model allelic effects as a function of axes of genetic variation, derived from a matrix of mean pairwise allele frequency differences between GWAS, and implemented in the MR-MEGA software. Through detailed simulations, we demonstrate increased power to detect association for MR-MEGA over fixed- and random-effects meta-analysis across a range of scenarios of heterogeneity in allelic effects between ethnic groups. We also demonstrate improved fine-mapping resolution, in loci containing a single causal variant, compared to these meta-analysis approaches and PAINTOR, and equivalent performance to MANTRA at reduced computational cost. Application of MR-MEGA to trans-ethnic GWAS of kidney function in 71,461 individuals indicates stronger signals of association than fixed-effects meta-analysis when heterogeneity in allelic effects is correlated with ancestry. Application of MR-MEGA to fine-mapping four type 2 diabetes susceptibility loci in 22,086 cases and 42,539 controls highlights: (i) strong evidence for heterogeneity in allelic effects that is correlated with ancestry only at the index SNP for the association signal at the CDKAL1 locus; and (ii) 99% credible sets with six or fewer variants for five distinct association signals.

Mitchell C, Ploem C, Chico V, Ormondroyd E, Hall A, Wallace S, Fay M, Goodwin D, Bell J, Phillips S et al. 2017. Exploring the potential duty of care in clinical genomics under UK law. Med Law Int, 17 (3), pp. 158-182. | Show Abstract | Read more

Genome-wide sequencing technologies are beginning to be used in projects that have both clinical diagnostic and research components. The clinical application of this technology, which generates a huge amount of information of varying diagnostic certainty, involves addressing a number of challenges to establish appropriate standards. In this article, we explore the way that UK law may respond to three of these key challenges and could establish new legal duties in relation to feedback of findings that are unrelated to the presenting condition (secondary, additional or incidental findings); duties towards genetic relatives as well as the patient and duties on the part of researchers and professionals who do not have direct contact with patients. When considering these issues, the courts will take account of European and international comparisons, developing guidance and relevant ethical, social and policy factors. The UK courts will also be strongly influenced by precedent set in case law.

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. | Citations: 1 (Web of Science Lite) | 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 CoVaccineHTTM, enhanced vaccine efficacy and sterile protection following malaria challenge. The CoVaccineHTTM 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, | Citations: 2 (Scopus) | 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.

Burnham KL, Davenport EE, Radhakrishnan J, Humburg P, Gordon AC, Hutton P, Svoren-Jabalera E, Garrard C, Hill AVS, Hinds CJ, Knight JC. 2017. Shared and Distinct Aspects of the Sepsis Transcriptomic Response to Fecal Peritonitis and Pneumonia. Am J Respir Crit Care Med, 196 (3), pp. 328-339. | Citations: 5 (Web of Science Lite) | Show Abstract | Read more

RATIONALE: Heterogeneity in the septic response has hindered efforts to understand pathophysiology and develop targeted therapies. Source of infection, with different causative organisms and temporal changes, might influence this heterogeneity. OBJECTIVES: To investigate individual and temporal variations in the transcriptomic response to sepsis due to fecal peritonitis, and to compare these with the same parameters in community-acquired pneumonia. METHODS: We performed genome-wide gene expression profiling in peripheral blood leukocytes of adult patients admitted to intensive care with sepsis due to fecal peritonitis (n = 117) or community-acquired pneumonia (n = 126), and of control subjects without sepsis (n = 10). MEASUREMENTS AND MAIN RESULTS: A substantial portion of the transcribed genome (18%) was differentially expressed compared with that of control subjects, independent of source of infection, with eukaryotic initiation factor 2 signaling being the most enriched canonical pathway. We identified two sepsis response signature (SRS) subgroups in fecal peritonitis associated with early mortality (P = 0.01; hazard ratio, 4.78). We defined gene sets predictive of SRS group, and serial sampling demonstrated that subgroup membership is dynamic during intensive care unit admission. We found that SRS is the major predictor of transcriptomic variation; a small number of genes (n = 263) were differentially regulated according to the source of infection, enriched for IFN signaling and antigen presentation. We define temporal changes in gene expression from disease onset involving phagosome formation as well as natural killer cell and IL-3 signaling. CONCLUSIONS: The majority of the sepsis transcriptomic response is independent of the source of infection and includes signatures reflecting immune response state and prognosis. A modest number of genes show evidence of specificity. Our findings highlight opportunities for patient stratification and precision medicine in sepsis.

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. | Citations: 2 (Web of Science Lite) | 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.

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: 4 (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.

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. | Citations: 1 (Scopus) | 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: 2 (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.

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-+. | Citations: 4 (Scopus) | 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.

Komdeur FL, Prins TM, van de Wall S, Plat A, Wisman GBA, Hollema H, Daemen T, Church DN, de Bruyn M, Nijman HW. 2017. CD103+ tumor-infiltrating lymphocytes are tumor-reactive intraepithelial CD8+ T cells associated with prognostic benefit and therapy response in cervical cancer. Oncoimmunology, 6 (9), pp. e1338230. | Show Abstract | Read more

Human papilloma virus (HPV)-induced cervical cancer constitutively expresses viral E6/E7 oncoproteins and is an excellent target for T cell-based immunotherapy. However, not all tumor-infiltrating T cells confer equal benefit to patients, with epithelial T cells being superior to stromal T cells. To assess whether the epithelial T cell biomarker CD103 could specifically discriminate the beneficial antitumor T cells, association of CD103 with clinicopathological variables and outcome was analyzed in the TCGA cervical cancer data set (n = 304) and by immunohistochemistry (IHC) in an independent cohort (n = 460). Localization of CD103+ cells in the tumor was assessed by immunofluorescence. Furthermore, use of CD103 as a response biomarker was assessed in an in vivo E6/E7+ tumor model. Our results show that CD103 gene expression was strongly correlated with cytotoxic T cell markers (e.g. CD8/GZMB/PD1) in the TCGA series. In line with this, CD103+ cells in the IHC series co-expressed CD8 and were preferentially located in cervical tumor epithelium. High CD103+ cell infiltration was strongly associated with an improved prognosis in both series, and appeared to be a better predictor of outcome than CD8. Interestingly, the prognostic benefit of CD103 in both series seemed limited to patients receiving radiotherapy. In a preclinical mouse model, HPV E6/E7-targeted therapeutic vaccination in combination with radiotherapy increased the intratumoral number of CD103+ CD8+ T cells, providing a potential mechanistic basis for our results. In conclusion, CD103 is a promising marker for rapid assessment of tumor-reactive T cell infiltration of cervical cancers and a promising response biomarker for E6/E7-targeted immunotherapy.

Cule M, Donnelly P. 2017. STOCHASTIC MODELLING AND INFERENCE IN ELECTRONIC HOSPITAL DATABASES FOR THE SPREAD OF INFECTIONS: CLOSTRIDIUM DIFFICILE TRANSMISSION IN OXFORDSHIRE HOSPITALS 2007-2010 ANNALS OF APPLIED STATISTICS, 11 (2), pp. 655-679. | Show Abstract | Read more

© 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 evi dence 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 Y-H, Yerges-Armstrong LM, Chou W-C, 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. | Citations: 3 (Web of Science Lite) | 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.

Titsias MK, Yau C. 2017. The Hamming Ball Sampler JOURNAL OF THE AMERICAN STATISTICAL ASSOCIATION, 112 (520), pp. 1598-1611. | Show Abstract | Read more

© 2017 The Author(s). Published with license by Taylor & Francis © 2017, © Michalis K. Titsias and Christopher Yau. We introduce the Hamming ball sampler, a novel Markov chain Monte Carlo algorithm, for efficient inference in statistical models involving high-dimensional discrete state spaces. The sampling scheme uses an auxiliary variable construction that adaptively truncates the model space allowing iterative exploration of the full model space. The approach generalizes conventional Gibbs sampling schemes for discrete spaces and provides an intuitive means for user-controlled balance between statistical efficiency and computational tractability. We illustrate the generic utility of our sampling algorithm through application to a range of statistical models. Supplementary materials for this article are available online.

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. | Citations: 1 (Web of Science Lite) | 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.

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. | Citations: 2 (Web of Science Lite) | 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.

Hu Z, Yau C, Ahmed AA. 2017. A pan-cancer genome-wide analysis reveals tumour dependencies by induction of nonsense-mediated decay. Nat Commun, 8 pp. 15943. | Citations: 3 (Scopus) | Show Abstract | Read more

Nonsense-mediated decay (NMD) eliminates transcripts with premature termination codons. Although NMD-induced loss-of-function has been shown to contribute to the genesis of particular cancers, its global functional consequence in tumours has not been characterized. Here we develop an algorithm to predict NMD and apply it on somatic mutations reported in The Cancer Genome Atlas. We identify more than 73 K mutations that are predicted to elicit NMD (NMD-elicit). NMD-elicit mutations in tumour suppressor genes (TSGs) are associated with significant reduction in gene expression. We discover cancer-specific NMD-elicit signatures in TSGs and cancer-associated genes. Our analysis reveals a previously unrecognized dependence of hypermutated tumours on hypofunction of genes that are involved in chromatin remodelling and translation. Half of hypermutated stomach adenocarcinomas are associated with NMD-elicit mutations of the translation initiators LARP4B and EIF5B. Our results unravel strong therapeutic opportunities by targeting tumour dependencies on NMD-elicit mutations.

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: 3 (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.

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.

Wang J, Mouradov D, Wang X, Jorissen RN, Chambers MC, Zimmerman LJ, Vasaikar S, Love CG, Li S, Lowes K et al. 2017. Colorectal Cancer Cell Line Proteomes Are Representative of Primary Tumors and Predict Drug Sensitivity. Gastroenterology, 153 (4), pp. 1082-1095. | Citations: 1 (Web of Science Lite) | Show Abstract | Read more

BACKGROUND AND AIMS: Proteomics holds promise for individualizing cancer treatment. We analyzed to what extent the proteomic landscape of human colorectal cancer (CRC) is maintained in established CRC cell lines and the utility of proteomics for predicting therapeutic responses. METHODS: Proteomic and transcriptomic analyses were performed on 44 CRC cell lines, compared against primary CRCs (n=95) and normal tissues (n=60), and integrated with genomic and drug sensitivity data. RESULTS: Cell lines mirrored the proteomic aberrations of primary tumors, in particular for intrinsic programs. Tumor relationships of protein expression with DNA copy number aberrations and signatures of post-transcriptional regulation were recapitulated in cell lines. The 5 proteomic subtypes previously identified in tumors were represented among cell lines. Nonetheless, systematic differences between cell line and tumor proteomes were apparent, attributable to stroma, extrinsic signaling, and growth conditions. Contribution of tumor stroma obscured signatures of DNA mismatch repair identified in cell lines with a hypermutation phenotype. Global proteomic data showed improved utility for predicting both known drug-target relationships and overall drug sensitivity as compared with genomic or transcriptomic measurements. Inhibition of targetable proteins associated with drug responses further identified corresponding synergistic or antagonistic drug combinations. Our data provide evidence for CRC proteomic subtype-specific drug responses. CONCLUSIONS: Proteomes of established CRC cell line are representative of primary tumors. Proteomic data tend to exhibit improved prediction of drug sensitivity as compared with genomic and transcriptomic profiles. Our integrative proteogenomic analysis highlights the potential of proteome profiling to inform personalized cancer medicine.

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: Clinicaltrials.gov 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: 4 (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. | Citations: 1 (Scopus) | 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. | Citations: 1 (Web of Science Lite) | 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 ClinicalTrials.gov Identifier: NCT02014870, first registered on December 5, 2013.

Manivanh L, Pierret A, Rattanavong S, Kounnavongsa O, Buisson Y, Elliott I, Maeght J-L, Xayyathip K, Silisouk J, Vongsouvath M et al. 2017. Burkholderia pseudomallei in a lowland rice paddy: seasonal changes and influence of soil depth and physico-chemical properties. Sci Rep, 7 (1), pp. 3031. | Show Abstract | Read more

Melioidosis, a severe infection with the environmental bacterium Burkholderia pseudomallei, is being recognised increasingly frequently. What determines its uneven distribution within endemic areas is poorly understood. We cultured soil from a rice field in Laos for B. pseudomallei at different depths on 4 occasions over a 13-month period. We also measured physical and chemical parameters in order to identify associated characteristics. Overall, 195 of 653 samples (29.7%) yielded B. pseudomallei. A higher prevalence of B. pseudomallei was found at soil depths greater than the 30 cm currently recommended for B. pseudomallei environmental sampling. B. pseudomallei was associated with a high soil water content and low total nitrogen, carbon and organic matter content. Our results suggested that a sampling grid of 25 five metre square quadrats (i.e. 25 × 25 m) should be sufficient to detect B. pseudomallei at a given location if samples are taken at a soil depth of at least 60 cm. However, culture of B. pseudomallei in environmental samples is difficult and liable to variation. Future studies should both rely on molecular approaches and address the micro-heterogeneity of soil when investigating physico-chemical associations with the presence of B. pseudomallei.

Burns A, Alsolami R, Becq J, Stamatopoulos B, Timbs A, Bruce D, Robbe P, Vavoulis D, Clifford R, Cabes M et al. 2018. Whole-genome sequencing of chronic lymphocytic leukaemia reveals distinct differences in the mutational landscape between IgHVmutand IgHVunmutsubgroups. Leukemia, 32 (2), pp. 332-342. | 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 IgHVmutand IgHVunmutsubtypes. 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 IgHVunmutCLL in IKZF3, SAMHD1,PAX5 and BIRC3. Finally, we show that IgHVunmutCLL is dominated by coding mutations in driver genes and an aging signature, whereas IgHVmutCLL 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.

Irshad S, Bansal M, Guarnieri P, Davis H, Zen AAH, 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 JOURNAL OF PATHOLOGY, 242 (2), pp. 178-192. | Citations: 2 (Scopus) | Show Abstract | Read more

© 2017 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland. 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.

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: 3 (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.

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.

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), pp. eaam6393-eaam6393. | Citations: 7 (Scopus) | Show Abstract | Read more

The malaria parasitePlasmodium falciparuminvades human red blood cells by a series of interactions between host and parasite surface proteins. By analyzing genome sequence data from human populations, including 1269 individuals from sub-Saharan Africa, we identify a diverse array of large copy-number variants affecting the host invasion receptor genesGYPAandGYPBWe find that a nearby association with severe malaria is explained by a complex structural rearrangement involving the loss ofGYPBand gain of twoGYPB-Ahybrid genes, which encode a serologically distinct blood group antigen known as Dantu. This variant reduces the risk of severe malaria by 40% and has recently increased in frequency in parts of Kenya, yet it appears to be absent from west Africa. These findings link structural variation of red blood cell invasion receptors with natural resistance to severe malaria.

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: 15 (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.

Alharbi NK, Padron-Regalado E, Thompson CP, Kupke A, Wells D, Sloan MA, Grehan K, Temperton N, Lambe T, Warimwe G et al. 2017. ChAdOx1 and MVA based vaccine candidates against MERS-CoV elicit neutralising antibodies and cellular immune responses in mice. Vaccine, 35 (30), pp. 3780-3788. | Citations: 3 (Web of Science Lite) | Show Abstract | Read more

The Middle East respiratory syndrome coronavirus (MERS-CoV) has infected more than 1900 humans, since 2012. The syndrome ranges from asymptomatic and mild cases to severe pneumonia and death. The virus is believed to be circulating in dromedary camels without notable symptoms since the 1980s. Therefore, dromedary camels are considered the only animal source of infection. Neither antiviral drugs nor vaccines are approved for veterinary or medical use despite active research on this area. Here, we developed four vaccine candidates against MERS-CoV based on ChAdOx1 and MVA viral vectors, two candidates per vector. All vaccines contained the full-length spike gene of MERS-CoV; ChAdOx1 MERS vaccines were produced with or without the leader sequence of the human tissue plasminogen activator gene (tPA) where MVA MERS vaccines were produced with tPA, but either the mH5 or F11 promoter driving expression of the spike gene. All vaccine candidates were evaluated in a mouse model in prime only or prime-boost regimens. ChAdOx1 MERS with tPA induced higher neutralising antibodies than ChAdOx1 MERS without tPA. A single dose of ChAdOx1 MERS with tPA elicited cellular immune responses as well as neutralising antibodies that were boosted to a significantly higher level by MVA MERS. The humoral immunogenicity of a single dose of ChAdOx1 MERS with tPA was equivalent to two doses of MVA MERS (also with tPA). MVA MERS with mH5 or F11 promoter induced similar antibody levels; however, F11 promoter enhanced the cellular immunogenicity of MVA MERS to significantly higher magnitudes. In conclusion, our study showed that MERS-CoV vaccine candidates could be optimized by utilising different viral vectors, various genetic designs of the vectors, or different regimens to increase immunogenicity. ChAdOx1 and MVA vectored vaccines have been safely evaluated in camels and humans and these MERS vaccine candidates should now be tested in camels and in clinical trials.

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+CD127lowCD25+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 signallingin vitroand to IL-2in vivo. The data presented in this paper support the research article "Human IL-6RhiTIGIT-CD4+CD127lowCD25+T cells display potentin vitrosuppressive capacity and a distinct Th17 profile" (Ferreira RC et al., 2017; doi: 10.1016/j.clim.2017.03.002) [1].

Salaga M, Leedham S. 2017. Young GI angle: Effective data presentation. United European Gastroenterol J, 5 (4), pp. 601-602. | Read more

Holmes CC, Walker SG. 2017. Assigning a value to a power likelihood in a general Bayesian model BIOMETRIKA, 104 (2), pp. 497-503. | Citations: 1 (Web of Science Lite) | Show Abstract | Read more

© 2017 Biometrika Trust. Bayesian robustness under model misspecification is a current area of active research. Among recent ideas is that of raising the likelihood function to a power. In this paper we discuss the choice of appropriate power and provide examples.

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: 8 (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. | Citations: 1 (Web of Science Lite) | 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. | Citations: 2 (Web of Science Lite) | 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.

Muranen TA, Greco D, Blomqvist C, Aittomäki K, Khan S, Hogervorst F, Verhoef S, Pharoah PDP, Dunning AM, Shah M et al. 2017. Genetic modifiers of CHEK2*1100delC-associated breast cancer risk. Genet Med, 19 (5), pp. 599-603. | Show Abstract | Read more

PURPOSE: CHEK2*1100delC is a founder variant in European populations that confers a two- to threefold increased risk of breast cancer (BC). Epidemiologic and family studies have suggested that the risk associated with CHEK2*1100delC is modified by other genetic factors in a multiplicative fashion. We have investigated this empirically using data from the Breast Cancer Association Consortium (BCAC). METHODS: Using genotype data from 39,139 (624 1100delC carriers) BC patients and 40,063 (224) healthy controls from 32 BCAC studies, we analyzed the combined risk effects of CHEK2*1100delC and 77 common variants in terms of a polygenic risk score (PRS) and pairwise interaction. RESULTS: The PRS conferred odds ratios (OR) of 1.59 (95% CI: 1.21-2.09) per standard deviation for BC for CHEK2*1100delC carriers and 1.58 (1.55-1.62) for noncarriers. No evidence of deviation from the multiplicative model was found. The OR for the highest quintile of the PRS was 2.03 (0.86-4.78) for CHEK2*1100delC carriers, placing them in the high risk category according to UK NICE guidelines. The OR for the lowest quintile was 0.52 (0.16-1.74), indicating a lifetime risk close to the population average. CONCLUSION: Our results confirm the multiplicative nature of risk effects conferred by CHEK2*1100delC and the common susceptibility variants. Furthermore, the PRS could identify carriers at a high lifetime risk for clinical actions.Genet Med advance online publication 06 October 2016.

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: 25 (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.

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: 4 (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.

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: 4 (Scopus) | 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.

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 homozygousBRCA2deletion 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 tumoralBRCA2deletion 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 M-H, 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. | Citations: 2 (Scopus) | 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.

Total publications on this page: 100

Total citations for publications on this page: 143