Recent publications

Diogo D, Tian C, Franklin CS, Alanne-Kinnunen M, March M, Spencer CCA, Vangjeli C, Weale ME, Mattsson H, Kilpeläinen E et al. 2018. Phenome-wide association studies across large population cohorts support drug target validation. Nat Commun, 9 (1), pp. 4285. | Show Abstract | Read more

Phenome-wide association studies (PheWAS) have been proposed as a possible aid in drug development through elucidating mechanisms of action, identifying alternative indications, or predicting adverse drug events (ADEs). Here, we select 25 single nucleotide polymorphisms (SNPs) linked through genome-wide association studies (GWAS) to 19 candidate drug targets for common disease indications. We interrogate these SNPs by PheWAS in four large cohorts with extensive health information (23andMe, UK Biobank, FINRISK, CHOP) for association with 1683 binary endpoints in up to 697,815 individuals and conduct meta-analyses for 145 mapped disease endpoints. Our analyses replicate 75% of known GWAS associations (P < 0.05) and identify nine study-wide significant novel associations (of 71 with FDR < 0.1). We describe associations that may predict ADEs, e.g., acne, high cholesterol, gout, and gallstones with rs738409 (p.I148M) in PNPLA3 and asthma with rs1990760 (p.T946A) in IFIH1. Our results demonstrate PheWAS as a powerful addition to the toolkit for drug discovery.

Young AI, Wauthier FL, Donnelly P. 2018. Identifying loci affecting trait variability and detecting interactions in genome-wide association studies. Nat Genet, | Show Abstract | Read more

Identification of genetic variants with effects on trait variability can provide insights into the biological mechanisms that control variation and can identify potential interactions. We propose a two-degree-of-freedom test for jointly testing mean and variance effects to identify such variants. We implement the test in a linear mixed model, for which we provide an efficient algorithm and software. To focus on biologically interesting settings, we develop a test for dispersion effects, that is, variance effects not driven solely by mean effects when the trait distribution is non-normal. We apply our approach to body mass index in the subsample of the UK Biobank population with British ancestry (n ~408,000) and show that our approach can increase the power to detect associated loci. We identify and replicate novel associations with significant variance effects that cannot be explained by the non-normality of body mass index, and we provide suggestive evidence for a connection between leptin levels and body mass index variability.

Bycroft C, Freeman C, Petkova D, Band G, Elliott LT, Sharp K, Motyer A, Vukcevic D, Delaneau O, O'Connell J et al. 2018. The UK Biobank resource with deep phenotyping and genomic data. Nature, 562 (7726), pp. 203-209. | Citations: 1 (Scopus) | Show Abstract | Read more

The UK Biobank project is a prospective cohort study with deep genetic and phenotypic data collected on approximately 500,000 individuals from across the United Kingdom, aged between 40 and 69 at recruitment. The open resource is unique in its size and scope. A rich variety of phenotypic and health-related information is available on each participant, including biological measurements, lifestyle indicators, biomarkers in blood and urine, and imaging of the body and brain. Follow-up information is provided by linking health and medical records. Genome-wide genotype data have been collected on all participants, providing many opportunities for the discovery of new genetic associations and the genetic bases of complex traits. Here we describe the centralized analysis of the genetic data, including genotype quality, properties of population structure and relatedness of the genetic data, and efficient phasing and genotype imputation that increases the number of testable variants to around 96 million. Classical allelic variation at 11 human leukocyte antigen genes was imputed, resulting in the recovery of signals with known associations between human leukocyte antigen alleles and many diseases.

Elliott LT, Sharp K, Alfaro-Almagro F, Shi S, Miller KL, Douaud G, Marchini J, Smith SM. 2018. Genome-wide association studies of brain imaging phenotypes in UK Biobank. Nature, 562 (7726), pp. 210-216. | Show Abstract | Read more

The genetic architecture of brain structure and function is largely unknown. To investigate this, we carried out genome-wide association studies of 3,144 functional and structural brain imaging phenotypes from UK Biobank (discovery dataset 8,428 subjects). Here we show that many of these phenotypes are heritable. We identify 148 clusters of associations between single nucleotide polymorphisms and imaging phenotypes that replicate at P < 0.05, when we would expect 21 to replicate by chance. Notable significant, interpretable associations include: iron transport and storage genes, related to magnetic susceptibility of subcortical brain tissue; extracellular matrix and epidermal growth factor genes, associated with white matter micro-structure and lesions; genes that regulate mid-line axon development, associated with organization of the pontine crossing tract; and overall 17 genes involved in development, pathway signalling and plasticity. Our results provide insights into the genetic architecture of the brain that are relevant to neurological and psychiatric disorders, brain development and ageing.

Halbroth BR, Sebastian S, Poyntz HC, Bregu M, Cottingham MG, Hill AVS, Spencer AJ. 2018. Development of a Molecular Adjuvant to Enhance Antigen-Specific CD8+ T Cell Responses. Sci Rep, 8 (1), pp. 15020. | Show Abstract | Read more

Despite promising progress in malaria vaccine development, an efficacious subunit vaccine against P. falciparum remains to be licensed and deployed. This study aimed to improve on the immunogenicity of the leading liver-stage vaccine candidate (ChAd63-MVA ME-TRAP), known to confer protection by eliciting high levels of antigen-specific CD8+ T cells. We previously showed fusion of ME-TRAP to the human MHC class II invariant chain (Ii) could enhance CD8+ T cell responses in non-human primates, but did not progress to clinical testing due to potential risk of auto-immunity by vaccination of humans with a self-antigen. Initial immunogenicity analyses of ME-TRAP fused to subdomains of the Ii showed that the Ii transmembrane domain alone can enhance CD8+ T cell responses. Subsequently, truncated Ii sequences with low homology to human Ii were developed and shown to enhance CD8+ T cell responses. By systematically mutating the TM domain sequence, multimerization of the Ii chain was shown to be important for immune enhancement. We subsequently identified several proteins from a variety of microbial pathogens with similar characteristics, that also enhance the CD8+ T cell response and could therefore be used in viral vector vaccines when potent cell mediated immunity is required.

Seelig E, Howlett J, Porter L, Truman L, Heywood J, Kennet J, Arbon EL, Anselmiova K, Walker NM, Atkar R et al. 2018. The DILfrequency study is an adaptive trial to identify optimal IL-2 dosing in patients with type 1 diabetes. JCI Insight, 3 (19), | Show Abstract | Read more

BACKGROUND: Type 1 diabetes (T1D) results from loss of immune regulation, leading to the development of autoimmunity to pancreatic β cells, involving autoreactive T effector cells (Teffs). Tregs, which prevent autoimmunity, require IL-2 for maintenance of immunosuppressive functions. Using a response-adaptive design, we aimed to determine the optimal regimen of aldesleukin (recombinant human IL-2) to physiologically enhance Tregs while limiting expansion of Teffs. METHODS: DILfrequency is a nonrandomized, open-label, response-adaptive study of participants, aged 18-70 years, with T1D. The initial learning phase allocated 12 participants to 6 different predefined regimens. Then, 3 cohorts of 8 participants were sequentially allocated dose frequencies, based on repeated interim analyses of all accumulated trial data. The coprimary endpoints were percentage change in Tregs and Teffs and CD25 (α subunit of the IL-2 receptor) expression by Tregs, from baseline to steady state. RESULTS: Thirty-eight participants were enrolled, with thirty-six completing treatment. The optimal regimen to maintain a steady-state increase in Tregs of 30% and CD25 expression of 25% without Teff expansion is 0.26 × 106 IU/m2 (95% CI -0.007 to 0.485) every 3 days. Tregs and CD25 were dose-frequency responsive, Teffs were not. The commonest adverse event was injection site reaction (464 of 694 events). CONCLUSIONS: Using a response-adaptive design, aldesleukin treatment can be optimized. Our methodology can generally be employed to immediately access proof of mechanism, thereby leading to more efficient and safe drug development. TRIAL REGISTRATION: International Standard Randomised Controlled Trial Number Register, ISRCTN40319192; ClinicalTrials.gov, NCT02265809. FUNDING: Sir Jules Thorn Trust, the Swiss National Science Foundation, Wellcome, JDRF, and NIHR Cambridge Biomedical Research Centre.

Gordon D, Dafinca R, Scaber J, Alegre-Abarrategui J, Farrimond L, Scott C, Biggs D, Kent L, Oliver PL, Davies B et al. 2018. Single-copy expression of an amyotrophic lateral sclerosis-linked TDP-43 mutation (M337V) in BAC transgenic mice leads to altered stress granule dynamics and progressive motor dysfunction. Neurobiol Dis, 121 pp. 148-162. | Show Abstract | Read more

Mutations in the gene encoding the RNA-binding protein TDP-43 cause amyotrophic lateral sclerosis (ALS), clinically and pathologically indistinguishable from the majority of 'sporadic' cases of ALS, establishing altered TDP-43 function and distribution as a primary mechanism of neurodegeneration. Transgenic mouse models in which TDP-43 is overexpressed only partially recapitulate the key cellular pathology of human ALS, but may also lead to non-specific toxicity. To avoid the potentially confounding effects of overexpression, and to maintain regulated spatio-temporal and cell-specific expression, we generated mice in which an 80 kb genomic fragment containing the intact human TDP-43 locus (either TDP-43WT or TDP-43M337V) and its regulatory regions was integrated into the Rosa26 (Gt(ROSA26)Sor) locus in a single copy. At 3 months of age, TDP-43M337V mice are phenotypically normal but by around 6 months develop progressive motor function deficits associated with loss of neuromuscular junction integrity, leading to a reduced lifespan. RNA sequencing shows that widespread mis-splicing is absent prior to the development of a motor phenotype, though differential expression analysis reveals a distinct transcriptional profile in pre-symptomatic TDP-43M337V spinal cords. Despite the presence of clear motor abnormalities, there was no evidence of TDP-43 cytoplasmic aggregation in vivo at any timepoint. In primary embryonic spinal motor neurons and in embryonic stem cell (ESC)-derived motor neurons, mutant TDP-43 undergoes cytoplasmic mislocalisation, and is associated with altered stress granule assembly and dynamics. Overall, this mouse model provides evidence that ALS may arise through acquired TDP-43 toxicity associated with defective stress granule function. The normal phenotype until 6 months of age can facilitate the study of early pathways underlying ALS.

Gola A, Silman D, Walters AA, Sridhar S, Uderhardt S, Salman AM, Halbroth BR, Bellamy D, Bowyer G, Powlson J et al. 2018. Prime and target immunization protects against liver-stage malaria in mice. Sci Transl Med, 10 (460), | Show Abstract | Read more

Despite recent advances in treatment and vector control, malaria is still a leading cause of death, emphasizing the need for an effective vaccine. The malaria life cycle can be subdivided into three stages: the invasion and growth within liver hepatocytes (pre-erythrocytic stage), the blood stage (erythrocytic stage), and, finally, the sexual stage (occurring within the mosquito vector). Antigen (Ag)-specific CD8+ T cells are effectively induced by heterologous prime-boost viral vector immunization and known to correlate with liver-stage protection. However, liver-stage malaria vaccines have struggled to generate and maintain the high numbers of Plasmodium-specific circulating T cells necessary to confer sterile protection. We describe an alternative "prime and target" vaccination strategy aimed specifically at inducing high numbers of tissue-resident memory T cells present in the liver at the time of hepatic infection. This approach bypasses the need for very high numbers of circulating T cells and markedly increases the efficacy of subunit immunization against liver-stage malaria with clinically relevant Ags and clinically tested viral vectors in murine challenge models. Translation to clinical use has begun, with encouraging results from a pilot safety and feasibility trial of intravenous chimpanzee adenovirus vaccination in humans. This work highlights the value of a prime-target approach for immunization against malaria and suggests that this strategy may represent a more general approach for prophylaxis or immunotherapy of other liver infections and diseases.

Kinchen J, Chen HH, Parikh K, Antanaviciute A, Jagielowicz M, Fawkner-Corbett D, Ashley N, Cubitt L, Mellado-Gomez E, Attar M et al. 2018. Structural Remodeling of the Human Colonic Mesenchyme in Inflammatory Bowel Disease. Cell, 175 (2), pp. 372-386.e17. | Show Abstract | Read more

Intestinal mesenchymal cells play essential roles in epithelial homeostasis, matrix remodeling, immunity, and inflammation. But the extent of heterogeneity within the colonic mesenchyme in these processes remains unknown. Using unbiased single-cell profiling of over 16,500 colonic mesenchymal cells, we reveal four subsets of fibroblasts expressing divergent transcriptional regulators and functional pathways, in addition to pericytes and myofibroblasts. We identified a niche population located in proximity to epithelial crypts expressing SOX6, F3 (CD142), and WNT genes essential for colonic epithelial stem cell function. In colitis, we observed dysregulation of this niche and emergence of an activated mesenchymal population. This subset expressed TNF superfamily member 14 (TNFSF14), fibroblastic reticular cell-associated genes, IL-33, and Lysyl oxidases. Further, it induced factors that impaired epithelial proliferation and maturation and contributed to oxidative stress and disease severity in vivo. Our work defines how the colonic mesenchyme remodels to fuel inflammation and barrier dysfunction in IBD.

Leonavicius K, Royer C, Preece C, Davies B, Biggins JS, Srinivas S. 2018. Mechanics of mouse blastocyst hatching revealed by a hydrogel-based microdeformation assay. Proc Natl Acad Sci U S A, 115 (41), pp. 10375-10380. | Show Abstract | Read more

Mammalian embryos are surrounded by an acellular shell, the zona pellucida. Hatching out of the zona is crucial for implantation and continued development of the embryo. Clinically, problems in hatching can contribute to failure in assisted reproductive intervention. Although hatching is fundamentally a mechanical process, due to limitations in methodology most studies focus on its biochemical properties. To understand the role of mechanical forces in hatching, we developed a hydrogel deformation-based method and analytical approach for measuring pressure in cyst-like tissues. Using this approach, we found that, in cultured blastocysts, pressure increased linearly, with intermittent falls. Inhibition of Na/K-ATPase led to a dosage-dependent reduction in blastocyst cavity pressure, consistent with its requirement for cavity formation. Reducing blastocyst pressure reduced the probability of hatching, highlighting the importance of mechanical forces in hatching. These measurements allowed us to infer details of microphysiology such as osmolarity, ion and water transport kinetics across the trophectoderm, and zona stiffness, allowing us to model the embryo as a thin-shell pressure vessel. We applied this technique to test whether cryopreservation, a process commonly used in assisted reproductive technology (ART), leads to alteration of the embryo and found that thawed embryos generated significantly lower pressure than fresh embryos, a previously unknown effect of cryopreservation. We show that reduced pressure is linked to delayed hatching. Our approach can be used to optimize in vitro fertilization (IVF) using precise measurement of embryo microphysiology. It is also applicable to other biological systems involving cavity formation, providing an approach for measuring forces in diverse contexts.

Westra H-J, Martínez-Bonet M, Onengut-Gumuscu S, Lee A, Luo Y, Teslovich N, Worthington J, Martin J, Huizinga T, Klareskog L et al. 2018. Fine-mapping and functional studies highlight potential causal variants for rheumatoid arthritis and type 1 diabetes. Nat Genet, 50 (10), pp. 1366-1374. | Show Abstract | Read more

To define potentially causal variants for autoimmune disease, we fine-mapped1,2 76 rheumatoid arthritis (11,475 cases, 15,870 controls)3 and type 1 diabetes loci (9,334 cases, 11,111 controls)4. After sequencing 799 1-kilobase regulatory (H3K4me3) regions within these loci in 568 individuals, we observed accurate imputation for 89% of common variants. We defined credible sets of ≤5 causal variants at 5 rheumatoid arthritis and 10 type 1 diabetes loci. We identified potentially causal missense variants at DNASE1L3, PTPN22, SH2B3, and TYK2, and noncoding variants at MEG3, CD28-CTLA4, and IL2RA. We also identified potential candidate causal variants at SIRPG and TNFAIP3. Using functional assays, we confirmed allele-specific protein binding and differential enhancer activity for three variants: the CD28-CTLA4 rs117701653 SNP, MEG3 rs34552516 indel, and TNFAIP3 rs35926684 indel.

Muriuki JM, Mentzer AJ, Kimita W, Ndungu FM, Macharia AW, Webb EL, Lule SA, Morovat A, Hill AVS, Bejon P et al. 2018. Iron Status and Associated Malaria Risk Among African Children. Clin Infect Dis, | Show Abstract | Read more

Background: It remains unclear whether improving iron status increases malaria risk, and few studies have looked at the effect of host iron status on subsequent malaria infection. We therefore aimed to determine whether a child's iron status influences their subsequent risk of malaria infection in sub-Saharan Africa. Methods: We assayed iron and inflammatory biomarkers from community-based cohorts of 1309 Kenyan and 1374 Ugandan children aged 0 - 7 years and conducted prospective surveillance for episodes of malaria. Poisson regression models were fitted to determine the effect of iron status on the incidence rate ratio of malaria using longitudinal data covering a period of 6 months. Models were adjusted for age, sex, parasitemia, inflammation and study site. Results: At baseline, the prevalence of iron deficiency (ID) was 36.9% and 34.6% in Kenyan and Ugandan children, respectively. Iron deficiency anemia (IDA) affected 23.6% of Kenyan and 17.6% of Ugandan children. Malaria risk was lower in children with ID (IRR = 0.7; 95% CI: 0.6, 0.8; p<0.001) and IDA (IRR = 0.7; 95% CI: 0.6, 0.9; p=0.006). Low transferrin saturation (<10%) was similarly associated with lower malaria risk (IRR = 0.8; 95% CI: 0.6, 0.9; p=0.016). However, variation in hepcidin, soluble transferrin receptors (sTfR) and hemoglobin / anemia was not associated with altered malaria risk. Conclusions: ID appears to protect against malaria infection in African children when defined using ferritin and transferrin saturation, but not when defined by hepcidin, sTfR or hemoglobin. Further research is required to determine causality.

Volpato V, Smith J, Sandor C, Ried JS, Baud A, Handel A, Newey SE, Wessely F, Attar M, Whiteley E et al. 2018. Reproducibility of Molecular Phenotypes after Long-Term Differentiation to Human iPSC-Derived Neurons: A Multi-Site Omics Study. Stem Cell Reports, 11 (4), pp. 897-911. | Show Abstract | Read more

Reproducibility in molecular and cellular studies is fundamental to scientific discovery. To establish the reproducibility of a well-defined long-term neuronal differentiation protocol, we repeated the cellular and molecular comparison of the same two iPSC lines across five distinct laboratories. Despite uncovering acceptable variability within individual laboratories, we detect poor cross-site reproducibility of the differential gene expression signature between these two lines. Factor analysis identifies the laboratory as the largest source of variation along with several variation-inflating confounders such as passaging effects and progenitor storage. Single-cell transcriptomics shows substantial cellular heterogeneity underlying inter-laboratory variability and being responsible for biases in differential gene expression inference. Factor analysis-based normalization of the combined dataset can remove the nuisance technical effects, enabling the execution of robust hypothesis-generating studies. Our study shows that multi-center collaborations can expose systematic biases and identify critical factors to be standardized when publishing novel protocols, contributing to increased cross-site reproducibility.

Hill A, Beitelshees M, Pfeifer B, Jones CH. 2018. Standardizing pneumococcal biofilm release to PncO expression, a predictive measurement of virulence Infection and Immunity, 86 (9), pp. IAI.00494-18-IAI.00494-18. | Show Abstract | Read more

© 2018 American Society for Microbiology. A critical component in clinical trials for vaccines against pneumococcal disease is the establishment of robust preclinical models and clinical correlates of protection, which, in the case of the causative bacterial agent Streptococcus pneumoniae, include standard sepsis/pneumonia mouse models and opsonophagocytic activity (OPA), respectively. Despite broad usage, these gold-standard measures are ill equipped to evaluate nontraditional antigens that target virulence factors beyond capsular polysaccharides and/or proteins not associated with colonization or routine growth. These assays are further complicated by observed inconsistencies in the expression of target protein antigens and in the quantity of usable bacteria provided from respective growth processes. In an effort to overcome these issues, we performed an extensive optimization study of the critical steps in a bacterial biofilm dispersion model (termed "the biofilm model") to identify conditions that yield the greatest quantity of released pneumococci displaying a consistent virulence phenotype. Using this knowledge, we developed a secondary immune absorbance assay to provide immediate insight into the phenotypic state of bacteria conditioned using the biofilm model. Specifically, positive correlations between the expression of PncO (a key virulence-associated protein antigen) and immune absorbance (R2= 0.96), capsule shedding, and OPA assay titers were translated into a predictive readout of virulence in sepsis and pneumonia challenge models. These results present a methodology for generating consistent lots of virulent bacteria to standardize inputs in preclinical and clinical models for testing vaccines against biofilm-associated bacteria.

Duyvesteyn HME, Kotecha A, Ginn HM, Hecksel CW, Beale EV, de Haas F, Evans G, Zhang P, Chiu W, Stuart DI. 2018. Machining protein microcrystals for structure determination by electron diffraction. Proc Natl Acad Sci U S A, 115 (38), pp. 9569-9573. | Show Abstract | Read more

We demonstrate that ion-beam milling of frozen, hydrated protein crystals to thin lamella preserves the crystal lattice to near-atomic resolution. This provides a vehicle for protein structure determination, bridging the crystal size gap between the nanometer scale of conventional electron diffraction and micron scale of synchrotron microfocus beamlines. The demonstration that atomic information can be retained suggests that milling could provide such detail on sections cut from vitrified cells.

Cross W, Kovac M, Mustonen V, Temko D, Davis H, Baker A-M, Biswas S, Arnold R, Chegwidden L, Gatenbee C et al. 2018. The evolutionary landscape of colorectal tumorigenesis. Nat Ecol Evol, 2 (10), pp. 1661-1672. | Show Abstract | Read more

The evolutionary events that cause colorectal adenomas (benign) to progress to carcinomas (malignant) remain largely undetermined. Using multi-region genome and exome sequencing of 24 benign and malignant colorectal tumours, we investigate the evolutionary fitness landscape occupied by these neoplasms. Unlike carcinomas, advanced adenomas frequently harbour sub-clonal driver mutations-considered to be functionally important in the carcinogenic process-that have not swept to fixation, and have relatively high genetic heterogeneity. Carcinomas are distinguished from adenomas by widespread aneusomies that are usually clonal and often accrue in a 'punctuated' fashion. We conclude that adenomas evolve across an undulating fitness landscape, whereas carcinomas occupy a sharper fitness peak, probably owing to stabilizing selection.

Wang DR, Agosto-Pérez FJ, Chebotarov D, Shi Y, Marchini J, Fitzgerald M, McNally KL, Alexandrov N, McCouch SR. 2018. An imputation platform to enhance integration of rice genetic resources. Nat Commun, 9 (1), pp. 3519. | Show Abstract | Read more

As sequencing and genotyping technologies evolve, crop genetics researchers accumulate increasing numbers of genomic data sets from various genotyping platforms on different germplasm panels. Imputation is an effective approach to increase marker density of existing data sets toward the goal of integrating resources for downstream applications. While a number of imputation software packages are available, the limitations to utilization for the rice community include high computational demand and lack of a reference panel. To address these challenges, we develop the Rice Imputation Server, a publicly available web application leveraging genetic information from a globally diverse rice reference panel assembled here. This resource allows researchers to benefit from increased marker density without needing to perform imputation on their own machines. We demonstrate improvements that imputed data provide to rice genome-wide association (GWA) results of grain amylose content and show that the major functional nucleotide polymorphism is tagged only in the imputed data set.

Lunter G. 2018. Haplotype matching in large cohorts using the Li and Stephens model. Bioinformatics, | Show Abstract | Read more

Motivation: The Li and Stephens model, which approximates the coalescent describing the pattern of variation in a population, underpins a range of key tools and results in genetics. Although highly efficient compared to the coalescent, standard implementations of this model still cannot deal with the very large reference cohorts that are starting to become available, and practical implementations use heuristics to achieve reasonable runtimes. Results: Here I describe a new, exact algorithm ("fastLS") that implements the Li and Stephens model and achieves runtimes independent of the size of the reference cohort. Key to achieving this runtime is the use of the Burrows-Wheeler transform, allowing the algorithm to efficiently identify partial haplotype matches across a cohort. I show that the proposed data structure is very similar to, and generalizes, Durbin's positional Burrows-Wheeler transform.

Bolli N, Maura F, Minvielle S, Gloznik D, Szalat R, Fullam A, Martincorena I, Dawson KJ, Samur MK, Zamora J et al. 2018. Genomic patterns of progression in smoldering multiple myeloma NATURE COMMUNICATIONS, 9 (1), | Show Abstract | Read more

© 2018, The Author(s). We analyzed whole genomes of unique paired samples from smoldering multiple myeloma (SMM) patients progressing to multiple myeloma (MM). We report that the genomic landscape, including mutational profile and structural rearrangements at the smoldering stage is very similar to MM. Paired sample analysis shows two different patterns of progression: a “static progression model”, where the subclonal architecture is retained as the disease progressed to MM suggesting that progression solely reflects the time needed to accumulate a sufficient disease burden; and a “spontaneous evolution model”, where a change in the subclonal composition is observed. We also observe that activation-induced cytidine deaminase plays a major role in shaping the mutational landscape of early subclinical phases, while progression is driven by APOBEC cytidine deaminases. These results provide a unique insight into myelomagenesis with potential implications for the definition of smoldering disease and timing of treatment initiation.

Bonsall D, Black S, Howe AY, Chase R, Ingravallo P, Pak I, Brown A, Smith DA, Bowden R, Barnes E. 2018. Characterization of hepatitis C virus resistance to grazoprevir reveals complex patterns of mutations following on-treatment breakthrough that are not observed at relapse. Infect Drug Resist, 11 pp. 1119-1135. | Show Abstract | Read more

Purpose: A detailed analysis of hepatitis C virus (HCV) resistance-associated substitutions (RASs) is required to understand why people fail direct-acting antiviral therapies. This study was conducted to assess RASs in an analysis of 2 trials evaluating the second-generation NS3/4A protease inhibitor grazoprevir (GZR) in combination with peginterferon/ribavirin. Patients and methods: From a total of 113 participants with HCV genotype 1 infection, RASs were evaluated in 25 patients who relapsed and 6 patients with on-treatment virologic breakthrough using consensus Sanger and clonal sequence analysis of NS3/NS4a genes, with in vitro testing of replicon mutants. Next-generation sequencing (NGS) was used in a subset of participants to assess minority variants and the evolution of the whole viral genome. Results: Baseline RASs did not predict treatment failure. Relapse was most commonly associated with RASs at D168, although additional RASs (Y56, R155 and A156) were also detected, particularly in participants with on-treatment breakthrough. Treatment-emergent RASs usually reverted to wild-type (WT), suggesting these mutations were associated with a negative fitness cost (confirmed using in vitro assays). NGS was the most sensitive assay for the detection of minor variants. Significant viral sequence divergence (up to 5.9% codons) was observed across whole genomes in association with the acquisition and reversion of RASs. Conclusion: Relapse with GZR and peginterferon/ribavirin is commonly associated with single RASs in NS3 that generally revert to WT, while breakthrough follows more complex patterns of viral resistance. NGS suggests that large diverse pools of viral quasispecies that emerge with RASs facilitate rapid viral evolution.

Small KS, Todorčević M, Civelek M, El-Sayed Moustafa JS, Wang X, Simon MM, Fernandez-Tajes J, Mahajan A, Horikoshi M, Hugill A et al. 2018. Author Correction: Regulatory variants at KLF14 influence type 2 diabetes risk via a female-specific effect on adipocyte size and body composition. Nat Genet, 50 (9), pp. 1342. | Show Abstract | Read more

In the version of this article originally published, minus signs were missing from the three β-values for BMI given in Table 1. The errors have been corrected in the HTML and PDF versions of the article.

Varlakhanova NV, Alvarez FJD, Brady TM, Tornabene BA, Hosford CJ, Chappie JS, Zhang P, Ford MGJ. 2018. Structures of the fungal dynamin-related protein Vps1 reveal a unique, open helical architecture. J Cell Biol, 217 (10), pp. 3608-3624. | Citations: 1 (European Pubmed Central) | Show Abstract | Read more

Dynamin-related proteins (DRPs) are large multidomain GTPases required for diverse membrane-remodeling events. DRPs self-assemble into helical structures, but how these structures are tailored to their cellular targets remains unclear. We demonstrate that the fungal DRP Vps1 primarily localizes to and functions at the endosomal compartment. We present crystal structures of a Vps1 GTPase-bundle signaling element (BSE) fusion in different nucleotide states to capture GTP hydrolysis intermediates and concomitant conformational changes. Using cryoEM, we determined the structure of full-length GMPPCP-bound Vps1. The Vps1 helix is more open and flexible than that of dynamin. This is due to further opening of the BSEs away from the GTPase domains. A novel interface between adjacent GTPase domains forms in Vps1 instead of the contacts between the BSE and adjacent stalks and GTPase domains as seen in dynamin. Disruption of this interface abolishes Vps1 function in vivo. Hence, Vps1 exhibits a unique helical architecture, highlighting structural flexibilities of DRP self-assembly.

Thomsen SK, Raimondo A, Hastoy B, Sengupta S, Dai X-Q, Bautista A, Censin J, Payne AJ, Umapathysivam MM, Spigelman AF et al. 2018. Type 2 diabetes risk alleles in PAM impact insulin release from human pancreatic β-cells. Nat Genet, 50 (8), pp. 1122-1131. | Show Abstract | Read more

The molecular mechanisms underpinning susceptibility loci for type 2 diabetes (T2D) remain poorly understood. Coding variants in peptidylglycine α-amidating monooxygenase (PAM) are associated with both T2D risk and insulinogenic index. Here, we demonstrate that the T2D risk alleles impact negatively on overall PAM activity via defects in expression and catalytic function. PAM deficiency results in reduced insulin content and altered dynamics of insulin secretion in a human β-cell model and primary islets from cadaveric donors. Thus, our results demonstrate a role for PAM in β-cell function, and establish molecular mechanisms for T2D risk alleles at this locus.

Ramond C, Beydag-Tasöz BS, Azad A, van de Bunt M, Petersen MBK, Beer NL, Glaser N, Berthault C, Gloyn AL, Hansson M et al. 2018. Understanding human fetal pancreas development using subpopulation sorting, RNA sequencing and single-cell profiling. Development, 145 (16), | Citations: 1 (European Pubmed Central) | Show Abstract | Read more

To decipher the populations of cells present in the human fetal pancreas and their lineage relationships, we developed strategies to isolate pancreatic progenitors, endocrine progenitors and endocrine cells. Transcriptome analysis of the individual populations revealed a large degree of conservation among vertebrates in the drivers of gene expression changes that occur at different steps of differentiation, although notably, sometimes, different members of the same gene family are expressed. The transcriptome analysis establishes a resource to identify novel genes and pathways involved in human pancreas development. Single-cell profiling further captured intermediate stages of differentiation and enabled us to decipher the sequence of transcriptional events occurring during human endocrine differentiation. Furthermore, we evaluate how well individual pancreatic cells derived in vitro from human pluripotent stem cells mirror the natural process occurring in human fetuses. This comparison uncovers a few differences at the progenitor steps, a convergence at the steps of endocrine induction, and the current inability to fully resolve endocrine cell subtypes in vitro.

Bowyer G, Rampling T, Powlson J, Morter R, Wright D, Hill AVS, Ewer KJ. 2018. Activation-induced Markers Detect Vaccine-Specific CD4⁺ T Cell Responses Not Measured by Assays Conventionally Used in Clinical Trials. Vaccines (Basel), 6 (3), | Show Abstract | Read more

Immunogenicity of T cell-inducing vaccines, such as viral vectors or DNA vaccines and Bacillus Calmette-Guérin (BCG), are frequently assessed by cytokine-based approaches. While these are sensitive methods that have shown correlates of protection in various vaccine studies, they only identify a small proportion of the vaccine-specific T cell response. Responses to vaccination are likely to be heterogeneous, particularly when comparing prime and boost or assessing vaccine performance across diverse populations. Activation-induced markers (AIM) can provide a broader view of the total antigen-specific T cell response to enable a more comprehensive evaluation of vaccine immunogenicity. We tested an AIM assay for the detection of vaccine-specific CD4⁺ and CD8⁺ T cell responses in healthy UK adults vaccinated with viral vectored Ebola vaccine candidates, ChAd3-EBO-Z and MVA-EBO-Z. We used the markers, CD25, CD134 (OX40), CD274 (PDL1), and CD107a, to sensitively identify vaccine-responsive T cells. We compared the use of OX40⁺CD25⁺ and OX40⁺PDL1⁺ in CD4⁺ T cells and OX40⁺CD25⁺ and CD25⁺CD107a⁺ in CD8⁺ T cells for their sensitivity, specificity, and associations with other measures of vaccine immunogenicity. We show that activation-induced markers can be used as an additional method of demonstrating vaccine immunogenicity, providing a broader picture of the global T cell response to vaccination.

Li X, Francies HE, Secrier M, Perner J, Miremadi A, Galeano-Dalmau N, Barendt WJ, Letchford L, Leyden GM, Goffin EK et al. 2018. Organoid cultures recapitulate esophageal adenocarcinoma heterogeneity providing a model for clonality studies and precision therapeutics. Nat Commun, 9 (1), pp. 2983. | Show Abstract | Read more

Esophageal adenocarcinoma (EAC) incidence is increasing while 5-year survival rates remain less than 15%. A lack of experimental models has hampered progress. We have generated clinically annotated EAC organoid cultures that recapitulate the morphology, genomic, and transcriptomic landscape of the primary tumor including point mutations, copy number alterations, and mutational signatures. Karyotyping of organoid cultures has confirmed polyclonality reflecting the clonal architecture of the primary tumor. Furthermore, subclones underwent clonal selection associated with driver gene status. Medium throughput drug sensitivity testing demonstrates the potential of targeting receptor tyrosine kinases and downstream mediators. EAC organoid cultures provide a pre-clinical tool for studies of clonal evolution and precision therapeutics.

Oliveira JJ, Karrar S, Rainbow DB, Pinder CL, Clarke P, Rubio García A, Al-Assar O, Burling K, Morris S, Stratton R et al. 2018. The plasma biomarker soluble SIGLEC-1 is associated with the type I interferon transcriptional signature, ethnic background and renal disease in systemic lupus erythematosus. Arthritis Res Ther, 20 (1), pp. 152. | Show Abstract | Read more

BACKGROUND: The molecular heterogeneity of autoimmune and inflammatory diseases has been one of the main obstacles to the development of safe and specific therapeutic options. Here, we evaluated the diagnostic and clinical value of a robust, inexpensive, immunoassay detecting the circulating soluble form of the monocyte-specific surface receptor sialic acid binding Ig-like lectin 1 (sSIGLEC-1). METHODS: We developed an immunoassay to measure sSIGLEC-1 in small volumes of plasma/serum from systemic lupus erythematosus (SLE) patients (n = 75) and healthy donors (n = 504). Samples from systemic sclerosis patients (n = 99) were studied as an autoimmune control. We investigated the correlation between sSIGLEC-1 and both monocyte surface SIGLEC-1 and type I interferon-regulated gene (IRG) expression. Associations of sSIGLEC-1 with clinical features were evaluated in an independent cohort of SLE patients (n = 656). RESULTS: Plasma concentrations of sSIGLEC-1 strongly correlated with expression of SIGLEC-1 on the surface of blood monocytes and with IRG expression in SLE patients. We found ancestry-related differences in sSIGLEC-1 concentrations in SLE patients, with patients of non-European ancestry showing higher levels compared to patients of European ancestry. Higher sSIGLEC-1 concentrations were associated with lower serum complement component 3 and increased frequency of renal complications in European patients, but not with the SLE Disease Activity Index clinical score. CONCLUSIONS: Our sSIGLEC-1 immunoassay provides a specific and easily assayed marker for monocyte-macrophage activation, and interferonopathy in SLE and other diseases. Further studies can extend its clinical associations and its potential use to stratify patients and as a secondary endpoint in clinical trials.

Bowyer G, Grobbelaar A, Rampling T, Venkatraman N, Morelle D, Ballou RW, Hill AVS, Ewer KJ. 2018. CXCR3+ T Follicular Helper Cells Induced by Co-Administration of RTS,S/AS01B and Viral-Vectored Vaccines Are Associated With Reduced Immunogenicity and Efficacy Against Malaria. Front Immunol, 9 (JUL), pp. 1660. | Show Abstract | Read more

A malaria vaccine strategy targeting multiple lifecycle stages may be required to achieve a high level of efficacy. In two Phase IIa clinical trials, we tested immunogenicity and efficacy of RTS,S/AS01B administered alone, in a staggered regimen with viral-vectored vaccines or co-administered with viral-vectored vaccines. RTS,S/AS01B induces high titers of antibody against sporozoites and viral-vectored vaccines ChAd63 ME-TRAP and MVA ME-TRAP induce potent T cell responses against infected hepatocytes. By combining these two strategies, we aimed to improve efficacy by inducing immune responses targeting multiple parasite antigens. Vaccination with RTS,S/AS01B alone or in a staggered regimen with viral vectors produced strong immune responses and demonstrated high levels of protection against controlled human malaria infection. However, concomitant administration of these vaccines significantly reduced humoral immunogenicity and protective efficacy. Strong Th1-biased cytokine responses induced by MVA ME-TRAP were associated with a skew in circulating T follicular helper cells toward a CXCR3+ phenotype and a reduction in antibody quantity and quality. This study illustrates that while a multistage-targeting vaccine strategy could provide high-level efficacy, the regimen design will require careful optimization.

Ghezraoui H, Oliveira C, Becker JR, Bilham K, Moralli D, Anzilotti C, Fischer R, Deobagkar-Lele M, Sanchiz-Calvo M, Fueyo-Marcos E et al. 2018. 53BP1 cooperation with the REV7-shieldin complex underpins DNA structure-specific NHEJ. Nature, 560 (7716), pp. 122-127. | Citations: 4 (European Pubmed Central) | Show Abstract | Read more

53BP1 governs a specialized, context-specific branch of the classical non-homologous end joining DNA double-strand break repair pathway. Mice lacking 53bp1 (also known as Trp53bp1) are immunodeficient owing to a complete loss of immunoglobulin class-switch recombination1,2, and reduced fidelity of long-range V(D)J recombination3. The 53BP1-dependent pathway is also responsible for pathological joining events at dysfunctional telomeres4, and its unrestricted activity in Brca1-deficient cellular and tumour models causes genomic instability and oncogenesis5-7. Cells that lack core non-homologous end joining proteins are profoundly radiosensitive8, unlike 53BP1-deficient cells9,10, which suggests that 53BP1 and its co-factors act on specific DNA substrates. Here we show that 53BP1 cooperates with its downstream effector protein REV7 to promote non-homologous end joining during class-switch recombination, but REV7 is not required for 53BP1-dependent V(D)J recombination. We identify shieldin-a four-subunit putative single-stranded DNA-binding complex comprising REV7, c20orf196 (SHLD1), FAM35A (SHLD2) and FLJ26957 (SHLD3)-as the factor that explains this specificity. Shieldin is essential for REV7-dependent DNA end-protection and non-homologous end joining during class-switch recombination, and supports toxic non-homologous end joining in Brca1-deficient cells, yet is dispensable for REV7-dependent interstrand cross-link repair. The 53BP1 pathway therefore comprises distinct double-strand break repair activities within chromatin and single-stranded DNA compartments, which explains both the immunological differences between 53bp1- and Rev7- deficient mice and the context specificity of the pathway.

Ehsan M, Kelly M, Hooper C, Yavari A, Beglov J, Bellahcene M, Ghataorhe K, Poloni G, Goel A, Kyriakou T et al. 2018. Mutant Muscle LIM Protein C58G causes cardiomyopathy through protein depletion. J Mol Cell Cardiol, 121 pp. 287-296. | Show Abstract | Read more

Cysteine and glycine rich protein 3 (CSRP3) encodes Muscle LIM Protein (MLP), a well-established disease gene for Hypertrophic Cardiomyopathy (HCM). MLP, in contrast to the proteins encoded by the other recognised HCM disease genes, is non-sarcomeric, and has important signalling functions in cardiomyocytes. To gain insight into the disease mechanisms involved, we generated a knock-in mouse (KI) model, carrying the well documented HCM-causing CSRP3 mutation C58G. In vivo phenotyping of homozygous KI/KI mice revealed a robust cardiomyopathy phenotype with diastolic and systolic left ventricular dysfunction, which was supported by increased heart weight measurements. Transcriptome analysis by RNA-seq identified activation of pro-fibrotic signalling, induction of the fetal gene programme and activation of markers of hypertrophic signalling in these hearts. Further ex vivo analyses validated the activation of these pathways at transcript and protein level. Intriguingly, the abundance of MLP decreased in KI/KI mice by 80% and in KI/+ mice by 50%. Protein depletion was also observed in cellular studies for two further HCM-causing CSRP3 mutations (L44P and S54R/E55G). We show that MLP depletion is caused by proteasome action. Moreover, MLP C58G interacts with Bag3 and results in a proteotoxic response in the homozygous knock-in mice, as shown by induction of Bag3 and associated heat shock proteins. In conclusion, the newly generated mouse model provides insights into the underlying disease mechanisms of cardiomyopathy caused by mutations in the non-sarcomeric protein MLP. Furthermore, our cellular experiments suggest that protein depletion and proteasomal overload also play a role in other HCM-causing CSPR3 mutations that we investigated, indicating that reduced levels of functional MLP may be a common mechanism for HCM-causing CSPR3 mutations.

Shafrir AL, Farland LV, Shah DK, Harris HR, Kvaskoff M, Zondervan K, Missmer SA. 2018. Risk for and consequences of endometriosis: A critical epidemiologic review. Best Pract Res Clin Obstet Gynaecol, 51 pp. 1-15. | Show Abstract | Read more

Endometriosis affects approximately 10% of women of reproductive age. Characteristics robustly associated with a greater risk for endometriosis include early age at menarche, short menstrual cycle length, and lean body size, whereas greater parity has been associated with a lower risk. Relationships with other potential characteristics including physical activity, dietary factors, and lactation have been less consistent, partially because of the need for rigorous data collection and a longitudinal study design. Critical methodologic complexities include the need for a clear case definition; valid selection of comparison/control groups; and consideration of diagnostic bias and reverse causation when exploring demographic characteristics, medical history, and lifestyle factors. Reviewers and editors must demand a detailed description of rigorous methods to facilitate comparison and replication to advance our understanding of endometriosis.

Dennison SM, Reichartz M, Seaton KE, Dutta S, Wille-Reece U, Hill AVS, Ewer KJ, Rountree W, Sarzotti-Kelsoe M, Ozaki DA et al. 2018. Qualified Biolayer Interferometry Avidity Measurements Distinguish the Heterogeneity of Antibody Interactions with Plasmodium falciparum Circumsporozoite Protein Antigens. J Immunol, 201 (4), pp. 1315-1326. | Show Abstract | Read more

Ab avidity is a measure of the overall strength of Ab-Ag interactions and hence is important for understanding the functional efficiency of Abs. In vaccine evaluations, Ab avidity measurements can provide insights into immune correlates of protection and generate hypotheses regarding mechanisms of protection to improve vaccine design to achieve higher levels of efficacy. The commonly used Ab avidity assays require the use of chaotropic reagents to measure avidity index. In this study, using real-time detection of Ab-Ag binding by biolayer interferometry (BLI) technique, we have developed a qualified assay for measuring avidity of vaccine-induced Abs specific for Plasmodium falciparum circumsporozoite protein (CSP) Ags. Human mAb derived from plasmablasts of recipients of RTS,S/AS01 (RTS,S), the most advanced malaria vaccine candidate, were used in the assay development to measure Ag-specific binding responses and rate constants of association and dissociation. The optimized BLI binding assay demonstrated 1) good precision (percentage of coefficient of variation <20), 2) high specificity, 3) a lower limit of detection and quantitation in the 0.3-3.3 nM range, and 4) a range of linearity up to 50-100 nM for the CSP Ags tested. Analysis of polyclonal sera of malaria vaccinees demonstrated the suitability of this method to distinguish among vaccinees and rank Ab responses by avidity. These results demonstrate that precise, specific, and sensitive BLI measurements of Ab avidity in polyclonal sera from malaria vaccinees can map Ab response heterogeneity and potentially help to determine the role of Ab avidity as an immune correlate of protection for vaccines.

Baker A-M, Cross W, Curtius K, Al Bakir I, Choi C-HR, Davis HL, Temko D, Biswas S, Martinez P, Williams MJ et al. 2018. Evolutionary history of human colitis-associated colorectal cancer. Gut, pp. gutjnl-2018-316191-gutjnl-2018-316191. | Show Abstract | Read more

OBJECTIVE: IBD confers an increased lifetime risk of developing colorectal cancer (CRC), and colitis-associated CRC (CA-CRC) is molecularly distinct from sporadic CRC (S-CRC). Here we have dissected the evolutionary history of CA-CRC using multiregion sequencing. DESIGN: Exome sequencing was performed on fresh-frozen multiple regions of carcinoma, adjacent non-cancerous mucosa and blood from 12 patients with CA-CRC (n=55 exomes), and key variants were validated with orthogonal methods. Genome-wide copy number profiling was performed using single nucleotide polymorphism arrays and low-pass whole genome sequencing on archival non-dysplastic mucosa (n=9), low-grade dysplasia (LGD; n=30), high-grade dysplasia (HGD; n=13), mixed LGD/HGD (n=7) and CA-CRC (n=19). Phylogenetic trees were reconstructed, and evolutionary analysis used to reveal the temporal sequence of events leading to CA-CRC. RESULTS: 10/12 tumours were microsatellite stable with a median mutation burden of 3.0 single nucleotide alterations (SNA) per Mb, ~20% higher than S-CRC (2.5 SNAs/Mb), and consistent with elevated ageing-associated mutational processes. Non-dysplastic mucosa had considerable mutation burden (median 47 SNAs), including mutations shared with the neighbouring CA-CRC, indicating a precancer mutational field. CA-CRCs were often near triploid (40%) or near tetraploid (20%) and phylogenetic analysis revealed that copy number alterations (CNAs) began to accrue in non-dysplastic bowel, but the LGD/HGD transition often involved a punctuated 'catastrophic' CNA increase. CONCLUSIONS: Evolutionary genomic analysis revealed precancer clones bearing extensive SNAs and CNAs, with progression to cancer involving a dramatic accrual of CNAs at HGD. Detection of the cancerised field is an encouraging prospect for surveillance, but punctuated evolution may limit the window for early detection.

Guidi LG, Holloway ZG, Arnoult C, Ray PF, Monaco AP, Molnár Z, Velayos-Baeza A. 2018. AU040320 deficiency leads to disruption of acrosome biogenesis and infertility in homozygous mutant mice. Sci Rep, 8 (1), pp. 10379. | Show Abstract | Read more

Study of knockout (KO) mice has helped understand the link between many genes/proteins and human diseases. Identification of infertile KO mice provides valuable tools to characterize the molecular mechanisms underlying gamete formation. The KIAA0319L gene has been described to have a putative association with dyslexia; surprisingly, we observed that homozygous KO males for AU040320, KIAA0319L ortholog, are infertile and present a globozoospermia-like phenotype. Mutant spermatozoa are mostly immotile and display a malformed roundish head with no acrosome. In round spermatids, proacrosomal vesicles accumulate close to the acroplaxome but fail to coalesce into a single acrosomal vesicle. In wild-type mice AU040320 localises to the trans-Golgi-Network of germ cells but cannot be detected in mature acrosomes. Our results suggest AU040320 may be necessary for the normal formation of proacrosomal vesicles or the recruitment of cargo proteins required for downstream events leading to acrosomal fusion. Mutations in KIAA0319L could lead to human infertility; we screened for KIAA0319L mutations in a selected cohort of globozoospermia patients in which no genetic abnormalities have been previously identified, but detected no pathogenic changes in this particular cohort.

Atcheson E, Bauza K, Salman AM, Alves E, Blight J, Viveros-Sandoval ME, Janse CJ, Khan SM, Hill AVS, Reyes-Sandoval A. 2018. Tailoring a Plasmodium vivax Vaccine To Enhance Efficacy through a Combination of a CSP Virus-Like Particle and TRAP Viral Vectors. Infect Immun, 86 (9), pp. IAI.00114-18-IAI.00114-18. | Citations: 1 (Scopus) | Show Abstract | Read more

Vivax malaria remains one of the most serious and neglected tropical diseases, with 132 to 391 million clinical cases per year and 2.5 billion people at risk of infection. A vaccine against Plasmodium vivax could have more impact than any other intervention, and the use of a vaccine targeting multiple antigens may result in higher efficacy against sporozoite infection than targeting a single antigen. Here, two leading P. vivax preerythrocytic vaccine candidate antigens, the P. vivax circumsporozoite protein (PvCSP) and the thrombospondin-related adhesion protein (PvTRAP) were delivered as a combined vaccine. This strategy provided a dose-sparing effect, with 100% sterile protection in mice using doses that individually conferred low or no protection, as with the unadjuvanted antigens PvTRAP (0%) and PvCSP (50%), and reached protection similar to that of adjuvanted components. Efficacy against malaria infection was assessed using a new mouse challenge model consisting of a double-transgenic Plasmodium berghei parasite simultaneously expressing PvCSP and PvTRAP used in mice immunized with the virus-like particle (VLP) Rv21 previously reported to induce high efficacy in mice using Matrix-M adjuvant, while PvTRAP was concomitantly administered in chimpanzee adenovirus and modified vaccinia virus Ankara (MVA) vectors (viral-vectored TRAP, or vvTRAP) to support effective induction of T cells. We examined immunity elicited by these vaccines in the context of two adjuvants approved for human use (AddaVax and Matrix-M). Matrix-M supported the highest anti-PvCSP antibody titers when combined with Rv21, and, interestingly, mixing PvCSP Rv21 and PvTRAP viral vectors enhanced immunity to malaria over levels provided by single vaccines.

Wedden S, Miller K, Frayling IM, Thomas T, Chefani A, Miller K, Hamblin A, Taylor JC, D'Arrigo C. 2018. Colorectal Cancer Stratification in the Routine Clinical Pathway: A District General Hospital Experience. Appl Immunohistochem Mol Morphol, pp. 1-1. | Show Abstract | Read more

Colorectal cancer (CRC) has many subtypes with different prognoses and response to treatment. Patients must be characterized to access the most appropriate treatment and improve outcomes. An increasing number of biomarkers are required for characterization but are not in routine use. We investigated whether CRC can be stratified routinely within a small district general hospital to inform clinical decision making at local multidisciplinary team meeting/tumor board level. We evaluated mismatch repair (MMR) and EGFR signaling pathways using predominantly in-house immunohistochemical (IHC) tests (MSH2, MSH6, MLH1, PMS2, BRAF-V600E, Her2, PTEN, cMET) as well as send away PCR/NGS tests (NRAS, KRAS, and BRAF). We demonstrated that many of the tests required for personalized treatment of CRC can be done locally and timely. Send away tests need to be requested shortly after cut-up and this needs to be firmly established in the tissue pathways for the results to be considered at multidisciplinary team meeting/tumor board. We have shown that MMR IHC combined with BRAFV600E IHC is practical and easy to perform in a small district general hospital, has full concordance with DNA-based tests and satisfies the latest NICE requirements for the identification of potential Lynch syndrome patients. We provide a framework for the interpretation and presentation of test results. It is a practical classification that clinical pathologists can use to communicate effectively with the clinical team. It is broadly based on molecular subtyping, firmly focused on treatment decisions and dependent on the panel of molecular tests currently available.This is an open-access article distributed under the terms of the Creative Commons Attribution-Non Commercial-No Derivatives License 4.0 (CCBY-NC-ND), where it is permissible to download and share the work provided it is properly cited. The work cannot be changed in any way or used commercially without permission from the journal. http://creativecommons.org/licenses/by-nc-nd/4.0/.

Auburn S, Benavente ED, Miotto O, Pearson RD, Amato R, Grigg MJ, Barber BE, William T, Handayuni I, Marfurt J et al. 2018. Genomic analysis of a pre-elimination Malaysian Plasmodium vivax population reveals selective pressures and changing transmission dynamics. Nat Commun, 9 (1), pp. 2585. | Show Abstract | Read more

The incidence of Plasmodium vivax infection has declined markedly in Malaysia over the past decade despite evidence of high-grade chloroquine resistance. Here we investigate the genetic changes in a P. vivax population approaching elimination in 51 isolates from Sabah, Malaysia and compare these with data from 104 isolates from Thailand and 104 isolates from Indonesia. Sabah displays extensive population structure, mirroring that previously seen with the emergence of artemisinin-resistant P. falciparum founder populations in Cambodia. Fifty-four percent of the Sabah isolates have identical genomes, consistent with a rapid clonal expansion. Across Sabah, there is a high prevalence of loci known to be associated with antimalarial drug resistance. Measures of differentiation between the three countries reveal several gene regions under putative selection in Sabah. Our findings highlight important factors pertinent to parasite resurgence and molecular cues that can be used to monitor low-endemic populations at the end stages of P. vivax elimination.

Botuyan MV, Cui G, Drané P, Oliveira C, Detappe A, Brault ME, Parnandi N, Chaubey S, Thompson JR, Bragantini B et al. 2018. Mechanism of 53BP1 activity regulation by RNA-binding TIRR and a designer protein. Nat Struct Mol Biol, 25 (7), pp. 591-600. | Citations: 1 (European Pubmed Central) | Show Abstract | Read more

Dynamic protein interaction networks such as DNA double-strand break (DSB) signaling are modulated by post-translational modifications. The DNA repair factor 53BP1 is a rare example of a protein whose post-translational modification-binding function can be switched on and off. 53BP1 is recruited to DSBs by recognizing histone lysine methylation within chromatin, an activity directly inhibited by the 53BP1-binding protein TIRR. X-ray crystal structures of TIRR and a designer protein bound to 53BP1 now reveal a unique regulatory mechanism in which an intricate binding area centered on an essential TIRR arginine residue blocks the methylated-chromatin-binding surface of 53BP1. A 53BP1 separation-of-function mutation that abolishes TIRR-mediated regulation in cells renders 53BP1 hyperactive in response to DSBs, highlighting the key inhibitory function of TIRR. This 53BP1 inhibition is relieved by TIRR-interacting RNA molecules, providing proof-of-principle of RNA-triggered 53BP1 recruitment to DSBs.

Van Gool IC, Rayner E, Osse EM, Nout RA, Creutzberg CL, Tomlinson IPM, Church DN, Smit VTHBM, de Wind N, Bosse T, Drost M. 2018. Adjuvant Treatment for POLE Proofreading Domain-Mutant Cancers: Sensitivity to Radiotherapy, Chemotherapy, and Nucleoside Analogues. Clin Cancer Res, 24 (13), pp. 3197-3203. | Show Abstract | Read more

Purpose: Pathogenic POLE proofreading domain mutations are found in many malignancies where they are associated with ultramutation and favorable prognosis. The extent to which this prognosis depends on their sensitivity to adjuvant treatment is unknown, as is the optimal therapy for advanced-staged or recurrent POLE-mutant cancers.Experimental Design: We examined the recurrence-free survival of women with POLE-mutant and POLE-wild-type endometrial cancers (EC) in the observation arm of the randomized PORTEC-1 endometrial cancer trial (N = 245 patients with stage I endometrial cancer for analysis). Sensitivity to radiotherapy and selected chemotherapeutics was compared between Pole-mutant mouse-derived embryonic stem (mES) cells, generated using CRISPR-Cas9 (Pole mutations D275A/E275A, and cancer-associated P286R, S297F, V411L) and isogenic wild-type cell lines.Results: In the observation arm of the PORTEC-1 trial (N = 245), women with POLE-mutant endometrial cancers (N = 16) had an improved recurrence-free survival (10-year recurrence-free survival 100% vs. 80.1% for POLE-wild-type; HR, 0.143; 95% confidence interval, 0.001-0.996; P = 0.049). Pole mutations did not increase sensitivity to radiotherapy nor to chemotherapeutics in mES cells. In contrast, Pole-mutant cells displayed significantly increased sensitivity to cytarabine and fludarabine (IC50Pole P286R-mutant vs. wild-type: 0.05 vs. 0.17 μmol/L for cytarabine, 4.62 vs. 11.1 μmol/L for fludarabine; P < 0.001 for both comparisons).Conclusions: The favorable prognosis of POLE-mutant cancers cannot be explained by increased sensitivity to currently used adjuvant treatments. These results support studies exploring minimization of adjuvant therapy for early-stage POLE-mutant cancers, including endometrial and colorectal cancers. Conversely, POLE mutations result in hypersensitivity to nucleoside analogues, suggesting the use of these compounds as a potentially effective targeted treatment for advanced-stage POLE-mutant cancers. Clin Cancer Res; 24(13); 3197-203. ©2018 AACR.

Temko D, Van Gool IC, Rayner E, Glaire M, Makino S, Brown M, Chegwidden L, Palles C, Depreeuw J, Beggs A et al. 2018. Somatic POLE exonuclease domain mutations are early events in sporadic endometrial and colorectal carcinogenesis, determining driver mutational landscape, clonal neoantigen burden and immune response. J Pathol, 245 (3), pp. 283-296. | Show Abstract | Read more

Genomic instability, which is a hallmark of cancer, is generally thought to occur in the middle to late stages of tumourigenesis, following the acquisition of permissive molecular aberrations such as TP53 mutation or whole genome doubling. Tumours with somatic POLE exonuclease domain mutations are notable for their extreme genomic instability (their mutation burden is among the highest in human cancer), distinct mutational signature, lymphocytic infiltrate, and excellent prognosis. To what extent these characteristics are determined by the timing of POLE mutations in oncogenesis is unknown. Here, we have shown that pathogenic POLE mutations are detectable in non-malignant precursors of endometrial and colorectal cancer. Using genome and exome sequencing, we found that multiple driver mutations in POLE-mutant cancers show the characteristic POLE mutational signature, including those in genes conventionally regarded as initiators of tumourigenesis. In POLE-mutant cancers, the proportion of monoclonal predicted neoantigens was similar to that in other cancers, but the absolute number was much greater. We also found that the prominent CD8+ T-cell infiltrate present in POLE-mutant cancers was evident in their precursor lesions. Collectively, these data indicate that somatic POLE mutations are early, quite possibly initiating, events in the endometrial and colorectal cancers in which they occur. The resulting early onset of genomic instability may account for the striking immune response and excellent prognosis of these tumours, as well as their early presentation. © 2018 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

López-Camacho C, Abbink P, Larocca RA, Dejnirattisai W, Boyd M, Badamchi-Zadeh A, Wallace ZR, Doig J, Velazquez RS, Neto RDL et al. 2018. Rational Zika vaccine design via the modulation of antigen membrane anchors in chimpanzee adenoviral vectors. Nat Commun, 9 (1), pp. 2441. | Show Abstract | Read more

Zika virus (ZIKV) emerged on a global scale and no licensed vaccine ensures long-lasting anti-ZIKV immunity. Here we report the design and comparative evaluation of four replication-deficient chimpanzee adenoviral (ChAdOx1) ZIKV vaccine candidates comprising the addition or deletion of precursor membrane (prM) and envelope, with or without its transmembrane domain (TM). A single, non-adjuvanted vaccination of ChAdOx1 ZIKV vaccines elicits suitable levels of protective responses in mice challenged with ZIKV. ChAdOx1 prME ∆TM encoding prM and envelope without TM provides 100% protection, as well as long-lasting anti-envelope immune responses and no evidence of in vitro antibody-dependent enhancement to dengue virus. Deletion of prM and addition of TM reduces protective efficacy and yields lower anti-envelope responses. Our finding that immunity against ZIKV can be enhanced by modulating antigen membrane anchoring highlights important parameters in the design of viral vectored ZIKV vaccines to support further clinical assessments.

Campbell KR, Yau C. 2018. Uncovering pseudotemporal trajectories with covariates from single cell and bulk expression data. Nat Commun, 9 (1), pp. 2442. | Show Abstract | Read more

Pseudotime algorithms can be employed to extract latent temporal information from cross-sectional data sets allowing dynamic biological processes to be studied in situations where the collection of time series data is challenging or prohibitive. Computational techniques have arisen from single-cell 'omics and cancer modelling where pseudotime can be used to learn about cellular differentiation or tumour progression. However, methods to date typically implicitly assume homogeneous genetic, phenotypic or environmental backgrounds, which becomes limiting as data sets grow in size and complexity. We describe a novel statistical framework that learns how pseudotime trajectories can be modulated through covariates that encode such factors. We apply this model to both single-cell and bulk gene expression data sets and show that the approach can recover known and novel covariate-pseudotime interaction effects. This hybrid regression-latent variable model framework extends pseudotemporal modelling from its most prevalent area of single cell genomics to wider applications.

Inshaw JRJ, Cutler AJ, Burren OS, Stefana MI, Todd JA. 2018. Approaches and advances in the genetic causes of autoimmune disease and their implications. Nat Immunol, 19 (7), pp. 674-684. | Citations: 1 (European Pubmed Central) | Show Abstract | Read more

Genome-wide association studies are transformative in revealing the polygenetic basis of common diseases, with autoimmune diseases leading the charge. Although the field is just over 10 years old, advances in understanding the underlying mechanistic pathways of these conditions, which result from a dense multifactorial blend of genetic, developmental and environmental factors, have already been informative, including insights into therapeutic possibilities. Nevertheless, the challenge of identifying the actual causal genes and pathways and their biological effects on altering disease risk remains for many identified susceptibility regions. It is this fundamental knowledge that will underpin the revolution in patient stratification, the discovery of therapeutic targets and clinical trial design in the next 20 years. Here we outline recent advances in analytical and phenotyping approaches and the emergence of large cohorts with standardized gene-expression data and other phenotypic data that are fueling a bounty of discovery and improved understanding of human physiology.

Lenglet M, Robriquet F, Schwarz K, Camps C, Couturier A, Hoogewijs D, Buffet A, Knight SJL, Gad S, Couvé S et al. 2018. Identification of a new VHL exon and complex splicing alterations in familial erythrocytosis or von Hippel-Lindau disease. Blood, 132 (5), pp. 469-483. | Show Abstract | Read more

Chuvash polycythemia is an autosomal recessive form of erythrocytosis associated with a homozygous p.Arg200Trp mutation in the von Hippel-Lindau (VHL) gene. Since this discovery, additional VHL mutations have been identified in patients with congenital erythrocytosis, in a homozygous or compound-heterozygous state. VHL is a major tumor suppressor gene, mutations in which were first described in patients presenting with VHL disease, which is characterized by the development of highly vascularized tumors. Here, we identify a new VHL cryptic exon (termed E1') deep in intron 1 that is naturally expressed in many tissues. More importantly, we identify mutations in E1' in 7 families with erythrocytosis (1 homozygous case and 6 compound-heterozygous cases with a mutation in E1' in addition to a mutation in VHL coding sequences) and in 1 large family with typical VHL disease but without any alteration in the other VHL exons. In this study, we show that the mutations induced a dysregulation of VHL splicing with excessive retention of E1' and were associated with a downregulation of VHL protein expression. In addition, we demonstrate a pathogenic role for synonymous mutations in VHL exon 2 that altered splicing through E2-skipping in 5 families with erythrocytosis or VHL disease. In all the studied cases, the mutations differentially affected splicing, correlating with phenotype severity. This study demonstrates that cryptic exon retention and exon skipping are new VHL alterations and reveals a novel complex splicing regulation of the VHL gene. These findings open new avenues for diagnosis and research regarding the VHL-related hypoxia-signaling pathway.

Batty EM, Chaemchuen S, Blacksell S, Richards AL, Paris D, Bowden R, Chan C, Lachumanan R, Day N, Donnelly P et al. 2018. Long-read whole genome sequencing and comparative analysis of six strains of the human pathogen Orientia tsutsugamushi. PLoS Negl Trop Dis, 12 (6), pp. e0006566. | Show Abstract | Read more

BACKGROUND: Orientia tsutsugamushi is a clinically important but neglected obligate intracellular bacterial pathogen of the Rickettsiaceae family that causes the potentially life-threatening human disease scrub typhus. In contrast to the genome reduction seen in many obligate intracellular bacteria, early genetic studies of Orientia have revealed one of the most repetitive bacterial genomes sequenced to date. The dramatic expansion of mobile elements has hampered efforts to generate complete genome sequences using short read sequencing methodologies, and consequently there have been few studies of the comparative genomics of this neglected species. RESULTS: We report new high-quality genomes of O. tsutsugamushi, generated using PacBio single molecule long read sequencing, for six strains: Karp, Kato, Gilliam, TA686, UT76 and UT176. In comparative genomics analyses of these strains together with existing reference genomes from Ikeda and Boryong strains, we identify a relatively small core genome of 657 genes, grouped into core gene islands and separated by repeat regions, and use the core genes to infer the first whole-genome phylogeny of Orientia. CONCLUSIONS: Complete assemblies of multiple Orientia genomes verify initial suggestions that these are remarkable organisms. They have larger genomes compared with most other Rickettsiaceae, with widespread amplification of repeat elements and massive chromosomal rearrangements between strains. At the gene level, Orientia has a relatively small set of universally conserved genes, similar to other obligate intracellular bacteria, and the relative expansion in genome size can be accounted for by gene duplication and repeat amplification. Our study demonstrates the utility of long read sequencing to investigate complex bacterial genomes and characterise genomic variation.

Jewson J, Smith JQ, Holmes C. 2018. Principles of Bayesian Inference Using General Divergence Criteria ENTROPY, 20 (6), pp. 442-442. | Show Abstract | Read more

© 2018 by the authors. When it is acknowledged that all candidate parameterised statistical models are misspecified relative to the data generating process, the decision maker (DM) must currently concern themselves with inference for the parameter value minimising the Kullback-Leibler (KL)-divergence between the model and this process (Walker, 2013). However, it has long been known that minimising the KL-divergence places a large weight on correctly capturing the tails of the sample distribution. As a result, the DM is required to worry about the robustness of their model to tail misspecifications if they want to conduct principled inference. In this paper we alleviate these concerns for the DM. We advance recent methodological developments in general Bayesian updating (Bissiri, Holmes & Walker, 2016) to propose a statistically well principled Bayesian updating of beliefs targeting the minimisation of more general divergence criteria. We improve both the motivation and the statistical foundations of existing Bayesian minimum divergence estimation (Hooker & Vidyashankar, 2014; Ghosh & Basu, 2016), allowing the well principled Bayesian to target predictions from the model that are close to the genuine model in terms of some alternative divergence measure to the KL-divergence. Our principled formulation allows us to consider a broader range of divergences than have previously been considered. In fact, we argue defining the divergence measure forms an important, subjective part of any statistical analysis, and aim to provide some decision theoretic rational for this selection. We illustrate how targeting alternative divergence measures can impact the conclusions of simple inference tasks, and discuss then how our methods might apply to more complicated, high dimensional models.

Thézé J, Li T, du Plessis L, Bouquet J, Kraemer MUG, Somasekar S, Yu G, de Cesare M, Balmaseda A, Kuan G et al. 2018. Genomic Epidemiology Reconstructs the Introduction and Spread of Zika Virus in Central America and Mexico. Cell Host Microbe, 23 (6), pp. 855-864.e7. | Citations: 1 (European Pubmed Central) | Show Abstract | Read more

The Zika virus (ZIKV) epidemic in the Americas established ZIKV as a major public health threat and uncovered its association with severe diseases, including microcephaly. However, genetic epidemiology in some at-risk regions, particularly Central America and Mexico, remains limited. We report 61 ZIKV genomes from this region, generated using metagenomic sequencing with ZIKV-specific enrichment, and combine phylogenetic, epidemiological, and environmental data to reconstruct ZIKV transmission. These analyses revealed multiple independent ZIKV introductions to Central America and Mexico. One introduction, likely from Brazil via Honduras, led to most infections and the undetected spread of ZIKV through the region from late 2014. Multiple lines of evidence indicate biannual peaks of ZIKV transmission in the region, likely driven by varying local environmental conditions for mosquito vectors and herd immunity. The spatial and temporal heterogeneity of ZIKV transmission in Central America and Mexico challenges arbovirus surveillance and disease control measures.

Ternette N, Olde Nordkamp MJM, Müller J, Anderson AP, Nicastri A, Hill AVS, Kessler BM, Li D. 2018. Immunopeptidomic Profiling of HLA-A2-Positive Triple Negative Breast Cancer Identifies Potential Immunotherapy Target Antigens. Proteomics, 18 (12), pp. e1700465. | Citations: 1 (European Pubmed Central) | Show Abstract | Read more

The recent development in immune checkpoint inhibitors and chimeric antigen receptor (CAR) T-cells in the treatment of cancer has not only demonstrated the potency of utilizing T-cell reactivity for cancer therapy, but has also highlighted the need for developing new approaches to discover targets suitable for such novel therapeutics. Here we analyzed the immunopeptidomes of six HLA-A2-positive triple negative breast cancer (TNBC) samples by nano-ultra performance liquid chromatography tandem mass spectrometry (nUPLC-MS2 ). Immunopeptidomic profiling identified a total of 19 675 peptides from tumor and adjacent normal tissue and 130 of the peptides were found to have higher abundance in tumor than in normal tissues. To determine potential therapeutic target proteins, we calculated the average tumor-associated cohort coverage (aTaCC) that represents the percentage coverage of each protein in this cohort by peptides that had higher tumoral abundance. Cofilin-1 (CFL-1), interleukin-32 (IL-32), proliferating cell nuclear antigen (PCNA), syntenin-1 (SDCBP), and ribophorin-2 (RPN-2) were found to have the highest aTaCC scores. We propose that these antigens could be evaluated further for their potential as targets in breast cancer immunotherapy and the small cohort immunopeptidomics analysis technique could be used in a wide spectrum of target discovery. Data are available via ProteomeXchange with identifier PXD009738.

Sun BB, Maranville JC, Peters JE, Stacey D, Staley JR, Blackshaw J, Burgess S, Jiang T, Paige E, Surendran P et al. 2018. Genomic atlas of the human plasma proteome. Nature, 558 (7708), pp. 73-79. | Citations: 6 (Scopus) | Show Abstract | Read more

Although plasma proteins have important roles in biological processes and are the direct targets of many drugs, the genetic factors that control inter-individual variation in plasma protein levels are not well understood. Here we characterize the genetic architecture of the human plasma proteome in healthy blood donors from the INTERVAL study. We identify 1,927 genetic associations with 1,478 proteins, a fourfold increase on existing knowledge, including trans associations for 1,104 proteins. To understand the consequences of perturbations in plasma protein levels, we apply an integrated approach that links genetic variation with biological pathway, disease, and drug databases. We show that protein quantitative trait loci overlap with gene expression quantitative trait loci, as well as with disease-associated loci, and find evidence that protein biomarkers have causal roles in disease using Mendelian randomization analysis. By linking genetic factors to diseases via specific proteins, our analyses highlight potential therapeutic targets, opportunities for matching existing drugs with new disease indications, and potential safety concerns for drugs under development.

Laisk T, Kukushkina V, Palmer D, Laber S, Chen C-Y, Ferreira T, Rahmioglu N, Zondervan K, Becker C, Smoller J et al. 2018. GWAS meta-analysis highlights the hypothalamic-pituitary-gonadal axis (HPG axis) in the genetic regulation of menstrual cycle length | Show Abstract | Read more

The normal menstrual cycle requires a delicate interplay between the hypothalamus, pituitary, and ovary. Therefore, its length is an important indicator of female reproductive health. Menstrual cycle length has been shown to be partially controlled by genetic factors, especially in the follicle stimulating hormone beta-subunit (FSHB) locus. GWAS meta-analysis of menstrual cycle length in 44,871 women of European ancestry confirmed the previously observed association with the FSHB locus and identified four additional novel signals in, or near, the GNRH1, PGR, NR5A2 and INS-IGF2 genes. These findings confirm the role of the hypothalamic-pituitary-gonadal axis in the genetic regulation of menstrual cycle length, but also highlight potential novel local regulatory mechanisms, such as those mediated by IGF2.

Willetts M, Hollowell S, Aslett L, Holmes C, Doherty A. 2018. Statistical machine learning of sleep and physical activity phenotypes from sensor data in 96,220 UK Biobank participants. Sci Rep, 8 (1), pp. 7961. | Show Abstract | Read more

Current public health guidelines on physical activity and sleep duration are limited by a reliance on subjective self-reported evidence. Using data from simple wrist-worn activity monitors, we developed a tailored machine learning model, using balanced random forests with Hidden Markov Models, to reliably detect a number of activity modes. We show that physical activity and sleep behaviours can be classified with 87% accuracy in 159,504 minutes of recorded free-living behaviours from 132 adults. These trained models can be used to infer fine resolution activity patterns at the population scale in 96,220 participants. For example, we find that men spend more time in both low- and high- intensity behaviours, while women spend more time in mixed behaviours. Walking time is highest in spring and sleep time lowest during the summer. This work opens the possibility of future public health guidelines informed by the health consequences associated with specific, objectively measured, physical activity and sleep behaviours.

Fernandez-Rozadilla C, Kartsonaki C, Woolley C, McClellan M, Whittington D, Horgan G, Leedham S, Kriaucionis S, East JE, Tomlinson I. 2018. Author Correction: Telomere length and genetics are independent colorectal tumour risk factors in an evaluation of biomarkers in normal bowel. Br J Cancer, 118 (12), pp. 1683. | Show Abstract | Read more

Since the publication of this paper, the authors noticed that James E. East was assigned to the incorrect affiliation. The affiliation information is provided correctly, above.

Binks S, Varley J, Lee W, Makuch M, Elliott K, Gelfand JM, Jacob S, Leite MI, Maddison P, Chen M et al. 2018. Distinct HLA associations of LGI1 and CASPR2-antibody diseases. Brain, 141 pp. 2263-2271. | Show Abstract | Read more

The recent biochemical distinction between antibodies against leucine-rich, glioma-inactivated-1 (LGI1), contactin-associated protein-2 (CASPR2) and intracellular epitopes of voltage-gated potassium-channels (VGKCs) demands aetiological explanations. Given established associations between human leucocyte antigen (HLA) alleles and adverse drug reactions, and our clinical observation of frequent adverse drugs reactions in patients with LGI1 antibodies, we compared HLA alleles between healthy controls (n = 5553) and 111 Caucasian patients with VGKC-complex autoantibodies. In patients with LGI1 antibodies (n = 68), HLA-DRB1*07:01 was strongly represented [odds ratio = 27.6 (95% confidence interval 12.9-72.2), P = 4.1 × 10-26]. In contrast, patients with CASPR2 antibodies (n = 31) showed over-representation of HLA-DRB1*11:01 [odds ratio = 9.4 (95% confidence interval 4.6-19.3), P = 5.7 × 10-6]. Other allelic associations for patients with LGI1 antibodies reflected linkage, and significant haplotypic associations included HLA-DRB1*07:01-DQA1*02:01-DQB1*02:02, by comparison to DRB1*11:01-DQA1*05:01-DQB1*03:01 in CASPR2-antibody patients. Conditional analysis in LGI1-antibody patients resolved further independent class I and II associations. By comparison, patients with both LGI1 and CASPR2 antibodies (n = 3) carried yet another complement of HLA variants, and patients with intracellular VGKC antibodies (n = 9) lacked significant HLA associations. Within LGI1- or CASPR2-antibody patients, HLA associations did not correlate with clinical features. In silico predictions identified unique CASPR2- and LGI1-derived peptides potentially presented by the respective over-represented HLA molecules. These highly significant HLA associations dichotomize the underlying immunology in patients with LGI1 or CASPR2 antibodies, and inform T cell specificities and cellular interactions at disease initiation.

Gallagher C, Makinen N, Harris H, Uimari O, Cook J, Shigesi N, Rahmioglu N, Ferreira T, Velez-Edwards D, Edwards T et al. 2018. Genome-wide association analysis identifies 27 novel loci associated with uterine leiomyomata revealing common genetic origins with endometriosis | Show Abstract | Read more

Uterine leiomyomata (UL), also known as uterine fibroids, are the most common neoplasms of the reproductive tract and the primary cause for hysterectomy, leading to considerable impact on women's lives as well as high economic burden. Genetic epidemiologic studies indicate that heritable risk factors contribute to UL pathogenesis. Previous genome-wide association studies (GWAS) identified five loci associated with UL at genome-wide significance (P < 5 x 10 -8 ). We conducted GWAS meta-analysis in 20,406 cases and 223,918 female controls of white European ancestry, identifying 24 genome-wide significant independent loci; 17 replicated in an unrelated cohort of 15,068 additional cases and 43,587 female controls. Aggregation of discovery and replication studies (35,474 cases and 267,505 female controls) revealed six additional significant loci. Interestingly, four of the 17 loci identified and replicated in these analyses have also been associated with risk for endometriosis, another common gynecologic disorder. These findings increase our understanding of the biological mechanisms underlying UL development, and suggest overlapping genetic origins with endometriosis.

Barazas M, Annunziato S, Pettitt SJ, de Krijger I, Ghezraoui H, Roobol SJ, Lutz C, Frankum J, Song FF, Brough R et al. 2018. The CST Complex Mediates End Protection at Double-Strand Breaks and Promotes PARP Inhibitor Sensitivity in BRCA1-Deficient Cells. Cell Rep, 23 (7), pp. 2107-2118. | Citations: 6 (Web of Science Lite) | Show Abstract | Read more

Selective elimination of BRCA1-deficient cells by inhibitors of poly(ADP-ribose) polymerase (PARP) is a prime example of the concept of synthetic lethality in cancer therapy. This interaction is counteracted by the restoration of BRCA1-independent homologous recombination through loss of factors such as 53BP1, RIF1, and REV7/MAD2L2, which inhibit end resection of DNA double-strand breaks (DSBs). To identify additional factors involved in this process, we performed CRISPR/SpCas9-based loss-of-function screens and selected for factors that confer PARP inhibitor (PARPi) resistance in BRCA1-deficient cells. Loss of members of the CTC1-STN1-TEN1 (CST) complex were found to cause PARPi resistance in BRCA1-deficient cells in vitro and in vivo. We show that CTC1 depletion results in the restoration of end resection and that the CST complex may act downstream of 53BP1/RIF1. These data suggest that, in addition to its role in protecting telomeres, the CST complex also contributes to protecting DSBs from end resection.

Todd JA. 2018. Evidence that UBASH3 is a causal gene for type 1 diabetes. Eur J Hum Genet, 26 (7), pp. 925-927. | Citations: 1 (European Pubmed Central) | Read more

Posth C, Naegele K, Colleran H, Valentin F, Bedford S, Kami KW, Shing R, Buckley H, Kinaston R, Walworth M et al. 2018. Language continuity despite population replacement in Remote Oceania NATURE ECOLOGY & EVOLUTION, 2 (4), pp. 731-740. | Read more

Grace C, Clarke R, Goel A, Farrall M, Watkins H, Hopewell JC. 2018. Lack of genetic support for shared aetiology of Coronary Artery Disease and Late-onset Alzheimer's disease. Sci Rep, 8 (1), pp. 7102. | Show Abstract | Read more

Epidemiological studies suggest a positive association between coronary artery disease (CAD) and late-onset Alzheimer's disease (LOAD). This large-scale genetic study brings together 'big data' resources to examine the causal impact of genetic determinants of CAD on risk of LOAD. A two-sample Mendelian randomization approach was adopted to estimate the causal effect of CAD on risk of LOAD using summary data from 60,801 CAD cases from CARDIoGRAMplusC4D and 17,008 LOAD cases from the IGAP Consortium. Additional analyses assessed the independent relevance of genetic associations at the APOE locus for both CAD and LOAD. Higher genetically determined risk of CAD was associated with a slightly higher risk of LOAD (Odds Ratio (OR) per log-odds unit of CAD [95% CI]: 1.07 [1.01-1.15]; p = 0.027). However, after exclusion of the APOE locus, the estimate of the causal effect of CAD for LOAD was attenuated and no longer significant (OR 0.94 [0.88-1.01]; p = 0.072). This Mendelian randomization study indicates that the APOE locus is the chief determinant of shared genetic architecture between CAD and LOAD, and suggests a lack of causal relevance of CAD for risk of LOAD after exclusion of APOE.

Coughlan L, Sridhar S, Payne R, Edmans M, Milicic A, Venkatraman N, Lugonja B, Clifton L, Qi C, Folegatti PM et al. 2018. Corrigendum to "Heterologous Two-dose Vaccination with Simian Adenovirus and Poxvirus Vectors Elicits Long-lasting Cellular Immunity to Influenza Virus A in Healthy Adults" [EBioMedicine 29 (2018) 146-154]. EBioMedicine, 31 pp. 321. | Show Abstract | Read more

© 2018 The Authors The authors wish to point out that L. Coughlan and S. Sridhar were both at the Jenner Institute, University of Oxford, OX3 7DQ, UK, when the work for the paper was completed.

Osgood JA, Knight JC. 2018. Translating GWAS in rheumatic disease: approaches to establishing mechanism and function for genetic associations with ankylosing spondylitis. Brief Funct Genomics, 17 (5), pp. 308-318. | Show Abstract | Read more

Ankylosing spondylitis (AS) is a highly heritable chronic inflammatory arthritis characterized by osteoproliferation, fusion of affected joints and systemic manifestations. Many disease associations for AS have been reported through genome-wide association studies; however, identifying modulated genes and functional mechanism remains challenging. This review summarizes current genetic associations involving AS and describes strategic approaches for functional follow-up of disease-associated variants. Fine mapping using methods leveraging Bayesian approaches are outlined. Evidence highlighting the importance of context specificity for regulatory variants is reviewed, noting current evidence in AS for the relevant cell and tissue type to conduct such analyses. Technological advances for understanding the regulatory landscape within which functional variants may act are discussed using exemplars. Approaches include defining regulatory elements based on chromatin accessibility, effects of variants on genes at a distance through evidence of physical interactions (chromatin conformation capture), expression quantitative trait loci mapping and single-cell methodologies. Opportunities for mechanistic studies to investigate the function of specific variants, regulatory elements and genes enabled by genome editing using clustered regularly interspaced short palindromic repeats/Cas9 are also described. Further progress in our understanding of the genetics of AS through functional genomic and epigenomic approaches offers new opportunities to understand mechanism and develop innovative treatments.

Mackley MP, Blair E, Parker M, Taylor JC, Watkins H, Ormondroyd E. 2018. Views of rare disease participants in a UK whole-genome sequencing study towards secondary findings: a qualitative study. Eur J Hum Genet, 26 (5), pp. 652-659. | Show Abstract | Read more

With large-scale genome sequencing initiatives underway, vast amounts of genomic data are being generated. Results-including secondary findings (SF)-are being returned, although policies around generation and management remain inconsistent. In order to inform relevant policy, it is essential that the views of stakeholders be considered-including participants who have made decisions about SF since the wider debate began. We conducted semi-structured interviews with sixteen rare disease patients and parents enroled in genome sequencing to explore views towards SF. Informed by extensive contact with the healthcare system, interviewees demonstrated high levels of understanding of genetic testing and held pragmatic views: many are content not knowing SF. Interviewees expressed trust in the system and healthcare providers, as well as an appreciation of limited resources; acknowledging existing disease burden, many preferred to focus on their primary condition. Many demonstrated an expectation for recontact and assumed the possibility of later access to initially declined SF. In the absence of such an infrastructure, it is important that responsibilities for recontact are delineated, expectations are addressed, and the long-term impact of decisions is made clear during consent. In addition, some interviewees demonstrated fluid views towards SF, and suggestions were made that perceptions may be influenced by family history. Further research into the changing desirability of SF and behavioural impact of disclosure are needed, and the development and introduction of mechanisms to respond to changes in patient views should be considered.

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

In the published version of this paper, the name of author Emanuele Di Angelantonio was misspelled. This error has now been corrected in the HTML and PDF versions of the article.

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

In the version of this article originally published, one of the two authors with the name Wei Zhao was omitted from the author list and the affiliations for both authors were assigned to the single Wei Zhao in the author list. In addition, the ORCID for Wei Zhao (Department of Biostatistics and Epidemiology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA) was incorrectly assigned to author Wei Zhou. The errors have been corrected in the HTML and PDF versions of the article.

Opi DH, Swann O, Macharia A, Uyoga S, Band G, Ndila CM, Harrison EM, Thera MA, Kone AK, Diallo DA et al. 2018. Two complement receptor one alleles have opposing associations with cerebral malaria and interact with α+thalassaemia. Elife, 7 | Citations: 1 (European Pubmed Central) | Show Abstract | Read more

Malaria has been a major driving force in the evolution of the human genome. In sub-Saharan African populations, two neighbouring polymorphisms in the Complement Receptor One (CR1) gene, named Sl2 and McCb, occur at high frequencies, consistent with selection by malaria. Previous studies have been inconclusive. Using a large case-control study of severe malaria in Kenyan children and statistical models adjusted for confounders, we estimate the relationship between Sl2 and McCb and malaria phenotypes, and find they have opposing associations. The Sl2 polymorphism is associated with markedly reduced odds of cerebral malaria and death, while the McCb polymorphism is associated with increased odds of cerebral malaria. We also identify an apparent interaction between Sl2 and α+thalassaemia, with the protective association of Sl2 greatest in children with normal α-globin. The complex relationship between these three mutations may explain previous conflicting findings, highlighting the importance of considering genetic interactions in disease-association studies.

Clarke R, Valdes-Marquez E, Hill M, Gordon J, Farrall M, Hamsten A, Watkins H, Hopewell JC. 2018. Plasma cytokines and risk of coronary heart disease in the PROCARDIS study. Open Heart, 5 (1), pp. e000807. | Show Abstract | Read more

Objective: The aims of the study were to examine the associations of plasma levels of five cytokines (interleukin (IL)-6, IL-5, interferon-gamma (IFN-γ), tumour necrosis factor-alpha (TNF-α) and IL-6 receptor (IL-6R)) and C reactive protein (CRP) with risk of coronary heart disease (CHD). Methods: In a case-control study of 931 CHD cases and 974 controls, logistic regression was used to estimate the OR and 95% CI of CHD for extreme thirds of biomarkers after adjustment for established risk factors. Sensitivity analyses were conducted in non-statin and in non-aspirin users. Results: Plasma levels of CRP were moderately correlated with IL-6 (r=0.45) in controls, but more weakly correlated with other cytokines. Likewise, all other cytokines were only weakly correlated with each other. After adjustment for established risk factors, the ORs (95% CI) for CHD comparing extreme thirds of cytokine levels (defined in controls) were 2.53 (1.86 to 3.43) for IL-6, 1.46 (1.11 to 1.93) for IL-5 and 1.46 (1.09 to 1.95) for IFN-γ, respectively. However, neither TNF-α, IL-6R nor CRP was significantly associated with CHD. After further adjustment for the associated cytokines, only IL-5 (1.34; 1.00 to 1.80) and IL-6 (2.39; 1.73 to 3.30) remained significantly associated with CHD. The risk associations of cytokines in non-users of statins or aspirin were comparable with the overall population. Conclusions: This study confirmed the importance of IL-6 as the most strongly associated cytokine with CHD risk, but also demonstrated novel and independent associations of IL-5 with CHD that warrant further investigation using larger panels of cytokines.

Gupta SK, Wesolowska-Andersen A, Ringgaard AK, Jaiswal H, Song L, Hastoy B, Ingvorsen C, Taheri-Ghahfarokhi A, Magnusson B, Maresca M et al. 2018. NKX6.1 induced pluripotent stem cell reporter lines for isolation and analysis of functionally relevant neuronal and pancreas populations. Stem Cell Res, 29 pp. 220-231. | Show Abstract | Read more

Recent studies have reported significant advances in the differentiation of human pluripotent stem cells to clinically relevant cell types such as the insulin producing beta-like cells and motor neurons. However, many of the current differentiation protocols lead to heterogeneous cell cultures containing cell types other than the targeted cell fate. Genetically modified human pluripotent stem cells reporting the expression of specific genes are of great value for differentiation protocol optimization and for the purification of relevant cell populations from heterogeneous cell cultures. Here we present the generation of human induced pluripotent stem cell (iPSC) lines with a GFP reporter inserted in the endogenous NKX6.1 locus. Characterization of the reporter lines demonstrated faithful GFP labelling of NKX6.1 expression during pancreas and motor neuron differentiation. Cell sorting and gene expression profiling by RNA sequencing revealed that NKX6.1-positive cells from pancreatic differentiations closely resemble human beta cells. Furthermore, functional characterization of the isolated cells demonstrated that glucose-stimulated insulin secretion is mainly confined to the NKX6.1-positive cells. We expect that the NKX6.1-GFP iPSC lines and the results presented here will contribute to the further refinement of differentiation protocols and characterization of hPSC-derived beta cells and motor neurons for disease modelling and cell replacement therapies.

Perez-Alcantara M, Honoré C, Wesolowska-Andersen A, Gloyn AL, McCarthy MI, Hansson M, Beer NL, van de Bunt M. 2018. Patterns of differential gene expression in a cellular model of human islet development, and relationship to type 2 diabetes predisposition. Diabetologia, 61 (7), pp. 1614-1622. | Show Abstract | Read more

AIMS/HYPOTHESIS: Most type 2 diabetes-associated genetic variants identified via genome-wide association studies (GWASs) appear to act via the pancreatic islet. Observed defects in insulin secretion could result from an impact of these variants on islet development and/or the function of mature islets. Most functional studies have focused on the latter, given limitations regarding access to human fetal islet tissue. Capitalising upon advances in in vitro differentiation, we characterised the transcriptomes of human induced pluripotent stem cell (iPSC) lines differentiated along the pancreatic endocrine lineage, and explored the contribution of altered islet development to the pathogenesis of type 2 diabetes. METHODS: We performed whole-transcriptome RNA sequencing of human iPSC lines from three independent donors, at baseline and at seven subsequent stages during in vitro islet differentiation. Differentially expressed genes (q < 0.01, log2 fold change [FC] > 1) were assigned to the stages at which they were most markedly upregulated. We used these data to characterise upstream transcription factors directing different stages of development, and to explore the relationship between RNA expression profiles and genes mapping to type 2 diabetes GWAS signals. RESULTS: We identified 9409 differentially expressed genes across all stages, including many known markers of islet development. Integration of differential expression data with information on transcription factor motifs highlighted the potential contribution of REST to islet development. Over 70% of genes mapping within type 2 diabetes-associated credible intervals showed peak differential expression during islet development, and type 2 diabetes GWAS loci of largest effect (including TCF7L2; log2FC = 1.2; q = 8.5 × 10-10) were notably enriched in genes differentially expressed at the posterior foregut stage (q = 0.002), as calculated by gene set enrichment analyses. In a complementary analysis of enrichment, genes differentially expressed in the final, beta-like cell stage of in vitro differentiation were significantly enriched (hypergeometric test, permuted p value <0.05) for genes within the credible intervals of type 2 diabetes GWAS loci. CONCLUSIONS/INTERPRETATION: The present study characterises RNA expression profiles during human islet differentiation, identifies potential transcriptional regulators of the differentiation process, and suggests that the inherited predisposition to type 2 diabetes is partly mediated through modulation of islet development. DATA AVAILABILITY: Sequence data for this study has been deposited at the European Genome-phenome Archive (EGA), under accession number EGAS00001002721.

Thorne L, Nalwoga A, Mentzer AJ, de Rougemont A, Hosmillo M, Webb E, Nampiija M, Muhwezi A, Carstensen T, Gurdasani D et al. 2018. The First Norovirus Longitudinal Seroepidemiological Study From Sub-Saharan Africa Reveals High Seroprevalence of Diverse Genotypes Associated With Host Susceptibility Factors. J Infect Dis, 218 (5), pp. 716-725. | Citations: 1 (European Pubmed Central) | Show Abstract | Read more

Background: Human noroviruses (HuNoVs) are a prominent cause of gastroenteritis, yet fundamental questions remain regarding epidemiology, diversity, and immunity in sub-Saharan African children. We investigated HuNoV seroprevalence and genetic and sociodemographic risk factors in Ugandan children. Methods: We randomly screened 797 participants of a longitudinal birth cohort (Entebbe, EMaBS) and 378 from a cross-sectional survey (rural Lake Victoria, LaVIISWA), for antibodies against HuNoV genotypes by ELISA. We used linear regression modeling to test for associations between HuNoV antibody levels and sociodemographic factors, and with the human susceptibility rs601338 FUT2 secretor SNP and histo-blood group antigens (A/B/O). Results: Of EMaBS participants, 76.6% were seropositive by age 1, rising to 94.5% by age 2 years. Seroprevalence in 1 year olds of the rural LaVIISWA survey was even higher (95%). In the birth cohort, 99% of seropositive 2 year olds had responses to multiple HuNoV genotypes. We identified associations between secretor status and genogroup GII antibody levels (GII.4 P = 3.1 × 10-52), as well as ABO and GI (GI.2 P = 2.1 × 10-12). Conclusions: HuNoVs are highly prevalent in Ugandan children, indicating a substantial burden of diarrhea-associated morbidity with recurrent infections. Public health interventions, including vaccination, and increased surveillance are urgently needed.

Lannagan TRM, Lee YK, Wang T, Roper J, Bettington ML, Fennell L, Vrbanac L, Jonavicius L, Somashekar R, Gieniec K et al. 2018. Genetic editing of colonic organoids provides a molecularly distinct and orthotopic preclinical model of serrated carcinogenesis. Gut, pp. gutjnl-2017-315920-gutjnl-2017-315920. | Citations: 1 (European Pubmed Central) | Show Abstract | Read more

OBJECTIVE: Serrated colorectal cancer (CRC) accounts for approximately 25% of cases and includes tumours that are among the most treatment resistant and with worst outcomes. This CRC subtype is associated with activating mutations in the mitogen-activated kinase pathway gene, BRAF, and epigenetic modifications termed the CpG Island Methylator Phenotype, leading to epigenetic silencing of key tumour suppressor genes. It is still not clear which (epi-)genetic changes are most important in neoplastic progression and we begin to address this knowledge gap herein. DESIGN: We use organoid culture combined with CRISPR/Cas9 genome engineering to sequentially introduce genetic alterations associated with serrated CRC and which regulate the stem cell niche, senescence and DNA mismatch repair. RESULTS: Targeted biallelic gene alterations were verified by DNA sequencing. Organoid growth in the absence of niche factors was assessed, as well as analysis of downstream molecular pathway activity. Orthotopic engraftment of complex organoid lines, but not BrafV600E alone, quickly generated adenocarcinoma in vivo with serrated features consistent with human disease. Loss of the essential DNA mismatch repair enzyme, Mlh1, led to microsatellite instability. Sphingolipid metabolism genes are differentially regulated in both our mouse models of serrated CRC and human CRC, with key members of this pathway having prognostic significance in the human setting. CONCLUSION: We generate rapid, complex models of serrated CRC to determine the contribution of specific genetic alterations to carcinogenesis. Analysis of our models alongside patient data has led to the identification of a potential susceptibility for this tumour type.

Wedge DC, Gundem G, Mitchell T, Woodcock DJ, Martincorena I, Ghori M, Zamora J, Butler A, Whitaker H, Kote-Jarai Z et al. 2018. Sequencing of prostate cancers identifies new cancer genes, routes of progression and drug targets. Nat Genet, 50 (5), pp. 682-692. | Citations: 3 (Scopus) | Show Abstract | Read more

Prostate cancer represents a substantial clinical challenge because it is difficult to predict outcome and advanced disease is often fatal. We sequenced the whole genomes of 112 primary and metastatic prostate cancer samples. From joint analysis of these cancers with those from previous studies (930 cancers in total), we found evidence for 22 previously unidentified putative driver genes harboring coding mutations, as well as evidence for NEAT1 and FOXA1 acting as drivers through noncoding mutations. Through the temporal dissection of aberrations, we identified driver mutations specifically associated with steps in the progression of prostate cancer, establishing, for example, loss of CHD1 and BRCA2 as early events in cancer development of ETS fusion-negative cancers. Computational chemogenomic (canSAR) analysis of prostate cancer mutations identified 11 targets of approved drugs, 7 targets of investigational drugs, and 62 targets of compounds that may be active and should be considered candidates for future clinical trials.

Pavić T, Juszczak A, Pape Medvidović E, Burrows C, Šekerija M, Bennett AJ, Ćuća Knežević J, Gloyn AL, Lauc G, McCarthy MI et al. 2018. Maturity onset diabetes of the young due to HNF1A variants in Croatia. Biochem Med (Zagreb), 28 (2), pp. 020703. | Show Abstract | Read more

Introduction: Maturity onset diabetes of the young due to HNF1A mutations (HNF1A-MODY) is the most frequent form of monogenic diabetes in adults. It is often misdiagnosed as type 1 or type 2 diabetes, but establishing genetic diagnosis is important, as treatment differs from the common types of diabetes. HNF1A-MODY has not been investigated in Croatia before due to limited access to genetic testing. In this study we aimed to describe the characteristics of young adults diagnosed with diabetes before the age of 45 years, who have rare HNF1A allele variants, and estimate the prevalence of HNF1A-MODY in Croatia. Materials and methods: We recruited 477 C-peptide positive and beta cell antibody negative subjects through the Croatian Diabetes Registry. HNF1A was sequenced for all participants and systematic assessment of the variants found was performed. The prevalence of HNF1A-MODY was calculated in the study group and results extrapolated to estimate the proportion of diabetic individuals with HNF1A-MODY in Croatia and the population prevalence. Results: Our study identified 13 individuals harbouring rare HNF1A allelic variants. After systematic assessment, 8 were assigned a diagnosis of HNF1A-MODY. Two individuals were able to discontinue insulin treatment following the diagnosis. We estimated that HNF1A-MODY in Croatia has a prevalence of 66 (95% CI 61 - 72) cases per million. Conclusions: The estimated prevalence of HNF1A-MODY in Croatia is similar to that reported in other European countries. Finding cases lead to important treatment changes for patients. This strongly supports the introduction of diagnostic genetic testing for monogenic diabetes in Croatia.

Reichold M, Klootwijk ED, Reinders J, Otto EA, Milani M, Broeker C, Laing C, Wiesner J, Devi S, Zhou W et al. 2018. Glycine Amidinotransferase (GATM), Renal Fanconi Syndrome, and Kidney Failure. J Am Soc Nephrol, 29 (7), pp. 1849-1858. | Citations: 3 (Web of Science Lite) | Show Abstract | Read more

Background For many patients with kidney failure, the cause and underlying defect remain unknown. Here, we describe a novel mechanism of a genetic order characterized by renal Fanconi syndrome and kidney failure.Methods We clinically and genetically characterized members of five families with autosomal dominant renal Fanconi syndrome and kidney failure. We performed genome-wide linkage analysis, sequencing, and expression studies in kidney biopsy specimens and renal cells along with knockout mouse studies and evaluations of mitochondrial morphology and function. Structural studies examined the effects of recognized mutations.Results The renal disease in these patients resulted from monoallelic mutations in the gene encoding glycine amidinotransferase (GATM), a renal proximal tubular enzyme in the creatine biosynthetic pathway that is otherwise associated with a recessive disorder of creatine deficiency. In silico analysis showed that the particular GATM mutations, identified in 28 members of the five families, create an additional interaction interface within the GATM protein and likely cause the linear aggregation of GATM observed in patient biopsy specimens and cultured proximal tubule cells. GATM aggregates-containing mitochondria were elongated and associated with increased ROS production, activation of the NLRP3 inflammasome, enhanced expression of the profibrotic cytokine IL-18, and increased cell death.Conclusions In this novel genetic disorder, fully penetrant heterozygous missense mutations in GATM trigger intramitochondrial fibrillary deposition of GATM and lead to elongated and abnormal mitochondria. We speculate that this renal proximal tubular mitochondrial pathology initiates a response from the inflammasome, with subsequent development of kidney fibrosis.

Ning J, Zhong Z, Fischer DK, Harris G, Watkins SC, Ambrose Z, Zhang P. 2018. Truncated CPSF6 Forms Higher-Order Complexes That Bind and Disrupt HIV-1 Capsid. J Virol, 92 (13), | Show Abstract | Read more

Cleavage and polyadenylation specificity factor 6 (CPSF6) is a human protein that binds HIV-1 capsid and mediates nuclear transport and integration targeting of HIV-1 preintegration complexes. Truncation of the protein at its C-terminal nuclear-targeting arginine/serine-rich (RS) domain produces a protein, CPSF6-358, that potently inhibits HIV-1 infection by targeting the capsid and inhibiting nuclear entry. To understand the molecular mechanism behind this restriction, the interaction between CPSF6-358 and HIV-1 capsid was characterized using in vitro and in vivo assays. Purified CPSF6-358 protein formed oligomers and bound in vitro-assembled wild-type (WT) capsid protein (CA) tubes, but not CA tubes containing a mutation in the putative binding site of CPSF6. Intriguingly, binding of CPSF6-358 oligomers to WT CA tubes physically disrupted the tubular assemblies into small fragments. Furthermore, fixed- and live-cell imaging showed that stably expressed CPSF6-358 forms cytoplasmic puncta upon WT HIV-1 infection and leads to capsid permeabilization. These events did not occur when the HIV-1 capsid contained a mutation known to prevent CPSF6 binding, nor did they occur in the presence of a small-molecule inhibitor of capsid binding to CPSF6-358. Together, our in vitro biochemical and transmission electron microscopy data and in vivo intracellular imaging results provide the first direct evidence for an oligomeric nature of CPSF6-358 and suggest a plausible mechanism for restriction of HIV-1 infection by CPSF6-358.IMPORTANCE After entry into cells, the HIV-1 capsid, which contains the viral genome, interacts with numerous host cell factors to facilitate crucial events required for replication, including uncoating. One such host cell factor, called CPSF6, is predominantly located in the cell nucleus and interacts with HIV-1 capsid. The interaction between CA and CPSF6 is critical during HIV-1 replication in vivo Truncation of CPSF6 leads to its localization to the cell cytoplasm and inhibition of HIV-1 infection. Here, we determined that truncated CPSF6 protein forms large higher-order complexes that bind directly to HIV-1 capsid, leading to its disruption. Truncated CPSF6 expression in cells leads to premature capsid uncoating that is detrimental to HIV-1 infection. Our study provides the first direct evidence for an oligomeric nature of truncated CPSF6 and insights into the highly regulated process of HIV-1 capsid uncoating.

Berlanga-Taylor AJ, Plant K, Dahl A, Lau E, Hill M, Sims D, Heger A, Emberson J, Armitage J, Clarke R, Knight JC. 2018. Genomic Response to Vitamin D Supplementation in the Setting of a Randomized, Placebo-Controlled Trial. EBioMedicine, 31 pp. 133-142. | Citations: 1 (Web of Science Lite) | Show Abstract | Read more

BACKGROUND: Vitamin D deficiency has been associated with multiple diseases, but the causal relevance and underlying processes are not fully understood. Elucidating the mechanisms of action of drug treatments in humans is challenging, but application of functional genomic approaches in randomized trials may afford an opportunity to systematically assess molecular responses. METHODS: In the Biochemical Efficacy and Safety Trial of Vitamin D (BEST-D), a double-blind, placebo-controlled, dose-finding, randomized clinical trial, 305 community-dwelling individuals aged over 65 years were randomly allocated to treatment with vitamin D3 4000 IU, 2000 IU or placebo daily for 12 months. Genome-wide genotypes at baseline, and transcriptome and plasma levels of cytokines (IFN-γ, IL-10, IL-8, IL-6 and TNF-α) at baseline and after 12 months, were measured. The trial had >90% power to detect 1.2-fold changes in gene expression. FINDINGS: Allocation to vitamin D for 12-months was associated with 2-fold higher plasma levels of 25-hydroxy-vitamin D (25[OH]D, 4000 IU regimen), but had no significant effect on whole-blood gene expression (FDR < 5%) or on plasma levels of cytokines compared with placebo. In pre-specified analysis, rs7041 (intron variant, GC) had a significant effect on circulating levels of 25(OH)D in the low dose, but not in the placebo or high dose vitamin D regimen. A gene expression quantitative trait locus analysis (eQTL) demonstrated evidence of 31,568 cis-eQTLs (unique SNP-probe pairs) among individuals at baseline and 34,254 after supplementation for 12 months (any dose). No significant associations involving vitamin D supplementation response eQTLs were found. INTERPRETATION: We performed a comprehensive functional genomics and molecular analysis of vitamin D supplementation in a randomized, placebo-controlled trial. Although this study was limited to mostly Caucasian individuals aged over 65 years, the results differ from many previous studies and do not support a strong effect of vitamin D on long-term transcriptomic changes in blood or on plasma cytokine levels. The trial demonstrates the feasibility of applying functional genomic and genetic approaches in randomized trials to assess molecular and individual level responses. KEY RESULT: Supplementation with high-dose vitamin D in older people for 12 months in a randomized, placebo-controlled trial had no significant effect on gene expression or on plasma concentrations of selected cytokines. TRIAL REGISTRATION: SRCTN registry (Number 07034656) and the European Clinical Trials Database (EudraCT Number 2011-005763-24).

Brown SDM, Holmes CC, Mallon A-M, Meehan TF, Smedley D, Wells S. 2018. High-throughput mouse phenomics for characterizing mammalian gene function. Nat Rev Genet, 19 (6), pp. 357-370. | Citations: 2 (Web of Science Lite) | Show Abstract | Read more

We are entering a new era of mouse phenomics, driven by large-scale and economical generation of mouse mutants coupled with increasingly sophisticated and comprehensive phenotyping. These studies are generating large, multidimensional gene-phenotype data sets, which are shedding new light on the mammalian genome landscape and revealing many hitherto unknown features of mammalian gene function. Moreover, these phenome resources provide a wealth of disease models and can be integrated with human genomics data as a powerful approach for the interpretation of human genetic variation and its relationship to disease. In the future, the development of novel phenotyping platforms allied to improved computational approaches, including machine learning, for the analysis of phenotype data will continue to enhance our ability to develop a comprehensive and powerful model of mammalian gene-phenotype space.

Schuh A, Dreau H, Knight SJL, Ridout K, Mizani T, Vavoulis D, Colling R, Antoniou P, Kvikstad EM, Pentony MM et al. 2018. Clinically actionable mutation profiles in patients with cancer identified by whole-genome sequencing. Cold Spring Harb Mol Case Stud, 4 (2), pp. a002279-a002279. | Show Abstract | Read more

Next-generation sequencing (NGS) efforts have established catalogs of mutations relevant to cancer development. However, the clinical utility of this information remains largely unexplored. Here, we present the results of the first eight patients recruited into a clinical whole-genome sequencing (WGS) program in the United Kingdom. We performed PCR-free WGS of fresh frozen tumors and germline DNA at 75× and 30×, respectively, using the HiSeq2500 HTv4. Subtracted tumor VCFs and paired germlines were subjected to comprehensive analysis of coding and noncoding regions, integration of germline with somatically acquired variants, and global mutation signatures and pathway analyses. Results were classified into tiers and presented to a multidisciplinary tumor board. WGS results helped to clarify an uncertain histopathological diagnosis in one case, led to informed or supported prognosis in two cases, leading to de-escalation of therapy in one, and indicated potential treatments in all eight. Overall 26 different tier 1 potentially clinically actionable findings were identified using WGS compared with six SNVs/indels using routine targeted NGS. These initial results demonstrate the potential of WGS to inform future diagnosis, prognosis, and treatment choice in cancer and justify the systematic evaluation of the clinical utility of WGS in larger cohorts of patients with cancer.

Painter JN, O'Mara TA, Morris AP, Cheng THT, Gorman M, Martin L, Hodson S, Jones A, Martin NG, Gordon S et al. 2018. Genetic overlap between endometriosis and endometrial cancer: evidence from cross-disease genetic correlation and GWAS meta-analyses. Cancer Med, 7 (5), pp. 1978-1987. | Citations: 2 (European Pubmed Central) | Show Abstract | Read more

Epidemiological, biological, and molecular data suggest links between endometriosis and endometrial cancer, with recent epidemiological studies providing evidence for an association between a previous diagnosis of endometriosis and risk of endometrial cancer. We used genetic data as an alternative approach to investigate shared biological etiology of these two diseases. Genetic correlation analysis of summary level statistics from genomewide association studies (GWAS) using LD Score regression revealed moderate but significant genetic correlation (rg  = 0.23, P = 9.3 × 10-3 ), and SNP effect concordance analysis provided evidence for significant SNP pleiotropy (P = 6.0 × 10-3 ) and concordance in effect direction (P = 2.0 × 10-3 ) between the two diseases. Cross-disease GWAS meta-analysis highlighted 13 distinct loci associated at P ≤ 10-5 with both endometriosis and endometrial cancer, with one locus (SNP rs2475335) located within PTPRD associated at a genomewide significant level (P = 4.9 × 10-8 , OR = 1.11, 95% CI = 1.07-1.15). PTPRD acts in the STAT3 pathway, which has been implicated in both endometriosis and endometrial cancer. This study demonstrates the value of cross-disease genetic analysis to support epidemiological observations and to identify biological pathways of relevance to multiple diseases.

Mahajan A, Wessel J, Willems SM, Zhao W, Robertson NR, Chu AY, Gan W, Kitajima H, Taliun D, Rayner NW et al. 2018. Refining the accuracy of validated target identification through coding variant fine-mapping in type 2 diabetes. Nat Genet, 50 (4), pp. 559-571. | Citations: 11 (Scopus) | Show Abstract | Read more

We aggregated coding variant data for 81,412 type 2 diabetes cases and 370,832 controls of diverse ancestry, identifying 40 coding variant association signals (P < 2.2 × 10-7); of these, 16 map outside known risk-associated loci. We make two important observations. First, only five of these signals are driven by low-frequency variants: even for these, effect sizes are modest (odds ratio ≤1.29). Second, when we used large-scale genome-wide association data to fine-map the associated variants in their regional context, accounting for the global enrichment of complex trait associations in coding sequence, compelling evidence for coding variant causality was obtained for only 16 signals. At 13 others, the associated coding variants clearly represent 'false leads' with potential to generate erroneous mechanistic inference. Coding variant associations offer a direct route to biological insight for complex diseases and identification of validated therapeutic targets; however, appropriate mechanistic inference requires careful specification of their causal contribution to disease predisposition.

Small KS, Todorčević M, Civelek M, El-Sayed Moustafa JS, Wang X, Simon MM, Fernandez-Tajes J, Mahajan A, Horikoshi M, Hugill A et al. 2018. Regulatory variants at KLF14 influence type 2 diabetes risk via a female-specific effect on adipocyte size and body composition. Nat Genet, 50 (4), pp. 572-580. | Citations: 6 (Web of Science Lite) | Show Abstract | Read more

Individual risk of type 2 diabetes (T2D) is modified by perturbations to the mass, distribution and function of adipose tissue. To investigate the mechanisms underlying these associations, we explored the molecular, cellular and whole-body effects of T2D-associated alleles near KLF14. We show that KLF14 diabetes-risk alleles act in adipose tissue to reduce KLF14 expression and modulate, in trans, the expression of 385 genes. We demonstrate, in human cellular studies, that reduced KLF14 expression increases pre-adipocyte proliferation but disrupts lipogenesis, and in mice, that adipose tissue-specific deletion of Klf14 partially recapitulates the human phenotype of insulin resistance, dyslipidemia and T2D. We show that carriers of the KLF14 T2D risk allele shift body fat from gynoid stores to abdominal stores and display a marked increase in adipocyte cell size, and that these effects on fat distribution, and the T2D association, are female specific. The metabolic risk associated with variation at this imprinted locus depends on the sex both of the subject and of the parent from whom the risk allele derives.

Behr M, Holmes C, Munk A. 2018. MULTISCALE BLIND SOURCE SEPARATION ANNALS OF STATISTICS, 46 (2), pp. 711-744. | Citations: 1 (Scopus) | Show Abstract | Read more

© Institute of Mathematical Statistics, 2018. We provide a new methodology for statistical recovery of single linear mixtures of piecewise constant signals (sources) with unknown mixing weights and change points in a multiscale fashion. We show exact recovery within an ε-neighborhood of the mixture when the sources take only values in a known finite alphabet. Based on this we provide the SLAM (Separates Linear Alphabet Mixtures) estimators for the mixing weights and sources. For Gaussian error, we obtain uniform confidence sets and optimal rates (up to log-factors) for all quantities. SLAM is efficiently computed as a nonconvex optimization problem by a dynamic program tailored to the finite alphabet assumption. Its performance is investigated in a simulation study. Finally, it is applied to assign copy-number aberrations from genetic sequencing data to different clones and to estimate their proportions.

Elliott I, Batty E, Ming D, Robinson M, Nawtaisong P, de Cesare M, Newton P, Bowden R. 2018. Oxford Nanopore MinION sequencing enables rapid whole-genome assembly of Rickettsia typhi in a resource-limited setting | Show Abstract | Read more

The infrastructure challenges and costs of next-generation sequencing have been largely overcome, for many sequencing applications, by Oxford Nanopore Technologies portable MinION sequencer. However the question remains open whether MinION-based bacterial whole-genome sequencing (WGS) is by itself sufficient for the accurate assessment of phylogenetic and epidemiological relationships between isolates and whether such tasks can be undertaken in resource-limited settings. To investigate this question, we sequenced the genome of an isolate of Rickettsia typhi, an important and neglected cause of fever across much of the tropics and subtropics, for which only three genomic sequences previously existed. We prepared and sequenced libraries on a MinION in Vientiane, Lao PDR using v9.5 chemistry and in parallel we sequenced the same isolate on the Illumina platform in a genomics laboratory in the UK. The MinION sequence reads yielded a single contiguous assembly, in which the addition of Illumina data revealed 226 base-substitution and 5,856 in/del errors. The combined assembly represents the first complete genome sequence of a human R. typhi isolate collected in the last 50 years and differed from the genomes of existing strains collected over a 90-year time period at very few sites, and with no re-arrangements. Filtering based on the known error profile of MinION data improved the accuracy of the Nanopore-only assembly. However, the frequency of false-positive errors remained greater than true sequence divergence from recorded sequences. While Nanopore-only sequencing cannot yet recover phylogenetic signal in R. typhi, such an approach may be applicable for more diverse organisms.

Gilchrist JJ, Mentzer AJ, Rautanen A, Pirinen M, Mwarumba S, Njuguna P, Mturi N, Wellcome Trust Case-Control Consortium 2, Kenyan Bacteraemia Study Group, Williams TN et al. 2018. Genetic variation in VAC14 is associated with bacteremia secondary to diverse pathogens in African children. Proc Natl Acad Sci U S A, 115 (16), pp. E3601-E3603. | Citations: 1 (European Pubmed Central) | Read more

Pagnamenta AT, Murakami Y, Anzilotti C, Titheradge H, Oates AJ, Morton J, DDD Study, Kinoshita T, Kini U, Taylor JC. 2018. A homozygous variant disrupting the PIGH start-codon is associated with developmental delay, epilepsy, and microcephaly. Hum Mutat, 39 (6), pp. 822-826. | Citations: 2 (European Pubmed Central) | Show Abstract | Read more

Defective glycosylphosphatidylinositol (GPI)-anchor biogenesis can cause a spectrum of predominantly neurological problems. For eight genes critical to this biological process, disease associations are not yet reported. Scanning exomes from 7,833 parent-child trios and 1,792 singletons from the DDD study for biallelic variants in this gene-set uncovered a rare PIGH variant in a boy with epilepsy, microcephaly, and behavioral difficulties. Although only 2/2 reads harbored this c.1A > T transversion, the presence of ∼25 Mb autozygosity at this locus implied homozygosity, which was confirmed using Sanger sequencing. A similarly-affected sister was also homozygous. FACS analysis of PIGH-deficient CHO cells indicated that cDNAs with c.1A > T could not efficiently restore expression of GPI-APs. Truncation of PIGH protein was consistent with the utilization of an in-frame start-site at codon 63. In summary, we describe siblings harboring a homozygous c.1A > T variant resulting in defective GPI-anchor biogenesis and highlight the importance of exploring low-coverage variants within autozygous regions.

Jiang J, Thalamuthu A, Ho JE, Mahajan A, Ek WE, Brown DA, Breit SN, Wang TJ, Gyllensten U, Chen M-H et al. 2018. A Meta-Analysis of Genome-Wide Association Studies of Growth Differentiation Factor-15 Concentration in Blood. Front Genet, 9 (MAR), pp. 97. | Show Abstract | Read more

Blood levels of growth differentiation factor-15 (GDF-15), also known as macrophage inhibitory cytokine-1 (MIC-1), have been associated with various pathological processes and diseases, including cardiovascular disease and cancer. Prior studies suggest genetic factors play a role in regulating blood MIC-1/GDF-15 concentration. In the current study, we conducted the largest genome-wide association study (GWAS) to date using a sample of ∼5,400 community-based Caucasian participants, to determine the genetic variants associated with MIC-1/GDF-15 blood concentration. Conditional and joint (COJO), gene-based association, and gene-set enrichment analyses were also carried out to identify novel loci, genes, and pathways. Consistent with prior results, a locus on chromosome 19, which includes nine single nucleotide polymorphisms (SNPs) (top SNP, rs888663, p = 1.690 × 10-35), was significantly associated with blood MIC-1/GDF-15 concentration, and explained 21.47% of its variance. COJO analysis showed evidence for two independent signals within this locus. Gene-based analysis confirmed the chromosome 19 locus association and in addition, a putative locus on chromosome 1. Gene-set enrichment analyses showed that the"COPI-mediated anterograde transport" gene-set was associated with MIC-1/GDF15 blood concentration with marginal significance after FDR correction (p = 0.067). In conclusion, a locus on chromosome 19 was associated with MIC-1/GDF-15 blood concentration with genome-wide significance, with evidence for a new locus (chromosome 1). Future studies using independent cohorts are needed to confirm the observed associations especially for the chromosomes 1 locus, and to further investigate and identify the causal SNPs that contribute to MIC-1/GDF-15 levels.

Saare M, Krigul KL, Laisk-Podar T, Ponandai-Srinivasan S, Rahmioglu N, Lalit Kumar PG, Zondervan K, Salumets A, Peters M. 2018. DNA methylation alterations-potential cause of endometriosis pathogenesis or a reflection of tissue heterogeneity? Biol Reprod, 99 (2), pp. 273-282. | Citations: 1 (European Pubmed Central) | Show Abstract | Read more

Alterations in the DNA methylation pattern of endometriotic lesions and endometrium of endometriosis patients have been proposed as one potential factor accompanying the endometriosis development. Although many differentially methylated genes have been associated with the pathogenesis of this disease, the overlap between the results of different studies has remained small. Among other potential confounders, the impact of tissue heterogeneity on the outcome of DNA methylation studies should be considered, as tissues are mixtures of different cell types with their own specific DNA methylation signatures. This review focuses on the results of DNA methylation studies in endometriosis from the cellular heterogeneity perspective. We consider both the studies using highly heterogeneous whole-lesion biopsies and endometrial tissue, as well as pure cell fractions isolated from lesions and endometrium to understand the potential impact of the cellular composition to the results of endometriosis DNA methylation studies. Also, future perspectives on how to diminish the impact of tissue heterogeneity in similar studies are provided.

Martin HC, Batty EM, Hussin J, Westall P, Daish T, Kolomyjec S, Piazza P, Bowden R, Hawkins M, Grant T et al. 2018. Insights into Platypus Population Structure and History from Whole-Genome Sequencing. Mol Biol Evol, 35 (5), pp. 1238-1252. | Citations: 1 (European Pubmed Central) | Show Abstract | Read more

The platypus is an egg-laying mammal which, alongside the echidna, occupies a unique place in the mammalian phylogenetic tree. Despite widespread interest in its unusual biology, little is known about its population structure or recent evolutionary history. To provide new insights into the dispersal and demographic history of this iconic species, we sequenced the genomes of 57 platypuses from across the whole species range in eastern mainland Australia and Tasmania. Using a highly improved reference genome, we called over 6.7 M SNPs, providing an informative genetic data set for population analyses. Our results show very strong population structure in the platypus, with our sampling locations corresponding to discrete groupings between which there is no evidence for recent gene flow. Genome-wide data allowed us to establish that 28 of the 57 sampled individuals had at least a third-degree relative among other samples from the same river, often taken at different times. Taking advantage of a sampled family quartet, we estimated the de novo mutation rate in the platypus at 7.0 × 10-9/bp/generation (95% CI 4.1 × 10-9-1.2 × 10-8/bp/generation). We estimated effective population sizes of ancestral populations and haplotype sharing between current groupings, and found evidence for bottlenecks and long-term population decline in multiple regions, and early divergence between populations in different regions. This study demonstrates the power of whole-genome sequencing for studying natural populations of an evolutionarily important species.

Chen X, Gustafsson S, Whitington T, Borné Y, Lorentzen E, Sun J, Almgren P, Su J, Karlsson R, Song J et al. 2018. A genome-wide association study of IgM antibody against phosphorylcholine: shared genetics and phenotypic relationship to chronic lymphocytic leukemia. Hum Mol Genet, 27 (10), pp. 1809-1818. | Show Abstract | Read more

Phosphorylcholine (PC) is an epitope on oxidized low-density lipoprotein (oxLDL), apoptotic cells and several pathogens like Streptococcus pneumoniae. Immunoglobulin M against PC (IgM anti-PC) has the ability to inhibit uptake of oxLDL by macrophages and increase clearance of apoptotic cells. From our genome-wide association studies (GWASs) in four European-ancestry cohorts, six single nucleotide polymorphisms (SNPs) in 11q24.1 were discovered (in 3002 individuals) and replicated (in 646 individuals) to be associated with serum level of IgM anti-PC (the leading SNP rs35923643-G, combined β = 0.19, 95% confidence interval 0.13-0.24, P = 4.3 × 10-11). The haplotype tagged by rs35923643-G (or its proxy SNP rs735665-A) is also known as the top risk allele for chronic lymphocytic leukemia (CLL), and a main increasing allele for general IgM. By using summary GWAS results of IgM anti-PC and CLL in the polygenic risk score (PRS) analysis, PRS on the basis of IgM anti-PC risk alleles positively associated with CLL risk (explained 0.6% of CLL variance, P = 1.2 × 10-15). Functional prediction suggested that rs35923643-G might impede the binding of Runt-related transcription factor 3, a tumor suppressor playing a central role in the immune regulation of cancers. Contrary to the expectations from the shared genetics between IgM anti-PC and CLL, an inverse relationship at the phenotypic level was found in a nested case-control study (30 CLL cases with 90 age- and sex-matched controls), potentially reflecting reverse causation. The suggested function of the top variant as well as the phenotypic association between IgM anti-PC and CLL risk needs replication and motivates further studies.

Ramamurthy N, Marchi E, Ansari MA, Pedergnana V, Mclean A, Hudson E, STOP HCV consortium, Bowden R, Spencer CCA, Barnes E, Klenerman P. 2018. Impact of IFNL4 genotype on Interferon-stimulated Gene Expression during DAA therapy for Hepatitis C. Hepatology, 68 (3), pp. 859-871. | Show Abstract | Read more

New directly acting antivirals (DAAs) provide very high cure rates in most patients infected by hepatitis C virus (HCV). However, some patient groups have been relatively harder to treat including those with cirrhosis or infected with HCV genotype 3. In the recent BOSON trial, genotype 3, cirrhotic patients receiving a 16 week course of sofosbuvir and ribavirin had a sustained virologic response rate (SVR) of around 50%. In cirrhotic patients, IFNL4 CC genotype was significantly associated with SVR. This genotype was also associated with a lower interferon-stimulated gene (ISG) signature in peripheral blood and in liver at baseline. Unexpectedly, patients with the CC genotype showed a dynamic increase in ISG expression between weeks 4 and 16 of DAA therapy, while the reverse was true for non-CC patients. These data provide an important dynamic link between host genotype and phenotype in HCV therapy also potentially relevant to naturally acquired infection. This article is protected by copyright. All rights reserved.

Malik R, Chauhan G, Traylor M, Sargurupremraj M, Okada Y, Mishra A, Rutten-Jacobs L, Giese A-K, van der Laan SW, Gretarsdottir S et al. 2018. Multiancestry genome-wide association study of 520,000 subjects identifies 32 loci associated with stroke and stroke subtypes. Nat Genet, 50 (4), pp. 524-537. | Citations: 12 (Web of Science Lite) | Show Abstract | Read more

Stroke has multiple etiologies, but the underlying genes and pathways are largely unknown. We conducted a multiancestry genome-wide-association meta-analysis in 521,612 individuals (67,162 cases and 454,450 controls) and discovered 22 new stroke risk loci, bringing the total to 32. We further found shared genetic variation with related vascular traits, including blood pressure, cardiac traits, and venous thromboembolism, at individual loci (n = 18), and using genetic risk scores and linkage-disequilibrium-score regression. Several loci exhibited distinct association and pleiotropy patterns for etiological stroke subtypes. Eleven new susceptibility loci indicate mechanisms not previously implicated in stroke pathophysiology, with prioritization of risk variants and genes accomplished through bioinformatics analyses using extensive functional datasets. Stroke risk loci were significantly enriched in drug targets for antithrombotic therapy.

Gilchrist JJ, Rautanen A, Fairfax BP, Mills TC, Naranbhai V, Trochet H, Pirinen M, Muthumbi E, Mwarumba S, Njuguna P et al. 2018. Risk of nontyphoidal Salmonella bacteraemia in African children is modified by STAT4. Nat Commun, 9 (1), pp. 1014. | Citations: 1 (Scopus) | Show Abstract | Read more

Nontyphoidal Salmonella (NTS) is a major cause of bacteraemia in Africa. The disease typically affects HIV-infected individuals and young children, causing substantial morbidity and mortality. Here we present a genome-wide association study (180 cases, 2677 controls) and replication analysis of NTS bacteraemia in Kenyan and Malawian children. We identify a locus in STAT4, rs13390936, associated with NTS bacteraemia. rs13390936 is a context-specific expression quantitative trait locus for STAT4 RNA expression, and individuals carrying the NTS-risk genotype demonstrate decreased interferon-γ (IFNγ) production in stimulated natural killer cells, and decreased circulating IFNγ concentrations during acute NTS bacteraemia. The NTS-risk allele at rs13390936 is associated with protection against a range of autoimmune diseases. These data implicate interleukin-12-dependent IFNγ-mediated immunity as a determinant of invasive NTS disease in African children, and highlight the shared genetic architecture of infectious and autoimmune disease.

Méric G, Mageiros L, Pascoe B, Woodcock DJ, Mourkas E, Lamble S, Bowden R, Jolley KA, Raymond B, Sheppard SK. 2018. Lineage-specific plasmid acquisition and the evolution of specialized pathogens in Bacillus thuringiensis and the Bacillus cereus group. Mol Ecol, 27 (7), pp. 1524-1540. | Citations: 1 (European Pubmed Central) | Show Abstract | Read more

Bacterial plasmids can vary from small selfish genetic elements to large autonomous replicons that constitute a significant proportion of total cellular DNA. By conferring novel function to the cell, plasmids may facilitate evolution but their mobility may be opposed by co-evolutionary relationships with chromosomes or encouraged via the infectious sharing of genes encoding public goods. Here, we explore these hypotheses through large-scale examination of the association between plasmids and chromosomal DNA in the phenotypically diverse Bacillus cereus group. This complex group is rich in plasmids, many of which encode essential virulence factors (Cry toxins) that are known public goods. We characterized population genomic structure, gene content and plasmid distribution to investigate the role of mobile elements in diversification. We analysed coding sequence within the core and accessory genome of 190 B. cereus group isolates, including 23 novel sequences and genes from 410 reference plasmid genomes. While cry genes were widely distributed, those with invertebrate toxicity were predominantly associated with one sequence cluster (clade 2) and phenotypically defined Bacillus thuringiensis. Cry toxin plasmids in clade 2 showed evidence of recent horizontal transfer and variable gene content, a pattern of plasmid segregation consistent with transfer during infectious cooperation. Nevertheless, comparison between clades suggests that co-evolutionary interactions may drive association between plasmids and chromosomes and limit wider transfer of key virulence traits. Proliferation of successful plasmid and chromosome combinations is a feature of specialized pathogens with characteristic niches (Bacillus anthracis, B. thuringiensis) and has occurred multiple times in the B. cereus group.

Lee W, Plant K, Humburg P, Knight JC. 2018. AltHapAlignR: improved accuracy of RNA-seq analyses through the use of alternative haplotypes. Bioinformatics, 34 (14), pp. 2401-2408. | Show Abstract | Read more

Motivation: Reliance on mapping to a single reference haplotype currently limits accurate estimation of allele or haplotype-specific expression using RNA-sequencing, notably in highly polymorphic regions such as the major histocompatibility complex. Results: We present AltHapAlignR, a method incorporating alternate reference haplotypes to generate gene- and haplotype-level estimates of transcript abundance for any genomic region where such information is available. We validate using simulated and experimental data to quantify input allelic ratios for major histocompatibility complex haplotypes, demonstrating significantly improved correlation with ground truth estimates of gene counts compared to standard single reference mapping. We apply AltHapAlignR to RNA-seq data from 462 individuals, showing how significant underestimation of expression of the majority of classical human leukocyte antigen genes using conventional mapping can be corrected using AltHapAlignR to allow more accurate quantification of gene expression for individual alleles and haplotypes. Availability and implementation: Source code freely available at https://github.com/jknightlab/AltHapAlignR. Supplementary information: Supplementary data are available at Bioinformatics online.

Matalliotakis M, Zervou MI, Eliopoulos E, Matalliotaki C, Rahmioglu N, Kalogiannidis I, Zondervan K, Spandidos DA, Matalliotakis I, Goulielmos GN. 2018. The role of IL‑16 gene polymorphisms in endometriosis. Int J Mol Med, 41 (3), pp. 1469-1476. | Show Abstract | Read more

Endometriosis is one of the most common gynecological diseases affecting up to 10% of the female population of childbearing age and a major cause of pain and infertility. It is influenced by multiple genetic, epigenetic and environmental factors. Interleukin‑16 (IL‑16) is a proinflammatory cytokine playing a pivotal role in many inflammatory and autoimmune diseases as well as in the pathogenesis of endometriosis. The aim of the present study was to evaluate the association of two IL‑16 gene single nucleotide polymorphisms (SNPs), rs4072111 and rs11556218, with the risk of endometriosis in women from Greece as well as to gain insight about the structural consequences of these two exonic SNPs regarding development of the disease. A total of 159 women with endometriosis (stages I‑IV) hospitalized for endometriosis, diagnosed by laparoscopic intervention and histologically confirmed, and 146 normal controls were recruited and genotyped. Subjects were genotyped using a polymerase chain reaction restriction fragment length polymorphism (PCR‑RFLP) strategy. A significant association was detected regarding the GG and GT genotype as well as 'G' allele of rs11556218 in patients with endometriosis. The rs4072111 SNP of the IL‑16 gene was not found to be associated with an increased susceptibility to endometriosis either for all patients (stages I‑IV) or for stage III and IV of the disease only. Our results demonstrated that rs11556218 is associated with endometriosis in Greek women, probably by resulting in the aberrant expression of IL‑16, as suggested by the bioinformatics analysis conducted on the SNP‑derived protein sequences, which indicated a possible association between mutation and functional modification of Pro‑IL‑16.

Coughlan L, Sridhar S, Payne R, Edmans M, Milicic A, Venkatraman N, Lugonja B, Clifton L, Qi C, Folegatti PM et al. 2018. Heterologous Two-Dose Vaccination with Simian Adenovirus and Poxvirus Vectors Elicits Long-Lasting Cellular Immunity to Influenza Virus A in Healthy Adults. EBioMedicine, 29 pp. 146-154. | Citations: 6 (Scopus) | Show Abstract | Read more

BACKGROUND: T-cell responses against highly conserved influenza antigens have been previously associated with protection. However, these immune responses are poorly maintained following recovery from influenza infection and are not boosted by inactivated influenza vaccines. We have previously demonstrated the safety and immunogenicity of two viral vectored vaccines, modified vaccinia virus Ankara (MVA) and the chimpanzee adenovirus ChAdOx1 expressing conserved influenza virus antigens, nucleoprotein (NP) and matrix protein-1 (M1). We now report on the safety and long-term immunogenicity of multiple combination regimes of these vaccines in young and older adults. METHODS: We conducted a Phase I open-label, randomized, multi-center study in 49 subjects aged 18-46years and 24 subjects aged 50years or over. Following vaccination, adverse events were recorded and the kinetics of the T cell response determined at multiple time points for up to 18months. FINDINGS: Both vaccines were well tolerated. A two dose heterologous vaccination regimen significantly increased the magnitude of pre-existing T-cell responses to NP and M1 after both doses in young and older adults. The fold-increase and peak immune responses after a single MVA-NP+M1 vaccination was significantly higher compared to ChAdOx1 NP+M1. In a mixed regression model, T-cell responses over 18months were significantly higher following the two dose vaccination regimen of MVA/ChAdOx1 NP+M1. INTERPRETATION: A two dose heterologous vaccination regimen of MVA/ChAdOx1 NP+M1 was safe and immunogenic in young and older adults, offering a promising vaccination strategy for inducing long-term broadly cross-reactive protection against influenza A. FUNDING SOURCE: Medical Research Council UK, NIHR BMRC Oxford.

Yoshida K, Iyori M, Blagborough AM, Salman AM, Dulal P, Sala KA, Yamamoto DS, Khan SM, Janse CJ, Biswas S et al. 2018. Adenovirus-prime and baculovirus-boost heterologous immunization achieves sterile protection against malaria sporozoite challenge in a murine model. Sci Rep, 8 (1), pp. 3896. | Citations: 1 (Scopus) | Show Abstract | Read more

With the increasing prevalence of artemisinin-resistant malaria parasites, a highly efficacious and durable vaccine for malaria is urgently required. We have developed an experimental virus-vectored vaccine platform based on an envelope-modified baculovirus dual-expression system (emBDES). Here, we show a conceptually new vaccine platform based on an adenovirus-prime/emBDES-boost heterologous immunization regimen expressing the Plasmodium falciparum circumsporozoite protein (PfCSP). A human adenovirus 5-prime/emBDES-boost heterologous immunization regimen consistently achieved higher sterile protection against transgenic P. berghei sporozoites expressing PfCSP after a mosquito-bite challenge than reverse-ordered or homologous immunization. This high protective efficacy was also achieved with a chimpanzee adenovirus 63-prime/emBDES-boost heterologous immunization regimen against an intravenous sporozoite challenge. Thus, we show that the adenovirus-prime/emBDES-boost heterologous immunization regimen confers sterile protection against sporozoite challenge by two individual routes, providing a promising new malaria vaccine platform for future clinical use.

Lipson M, Skoglund P, Spriggs M, Valentin F, Bedford S, Shing R, Buckley H, Phillip I, Ward GK, Mallick S et al. 2018. Population Turnover in Remote Oceania Shortly after Initial Settlement. Curr Biol, 28 (7), pp. 1157-1165.e7. | Citations: 5 (Scopus) | Show Abstract | Read more

Ancient DNA from Vanuatu and Tonga dating to about 2,900-2,600 years ago (before present, BP) has revealed that the "First Remote Oceanians" associated with the Lapita archaeological culture were directly descended from the population that, beginning around 5000 BP, spread Austronesian languages from Taiwan to the Philippines, western Melanesia, and eventually Remote Oceania. Thus, ancestors of the First Remote Oceanians must have passed by the Papuan-ancestry populations they encountered in New Guinea, the Bismarck Archipelago, and the Solomon Islands with minimal admixture [1]. However, all present-day populations in Near and Remote Oceania harbor >25% Papuan ancestry, implying that additional eastward migration must have occurred. We generated genome-wide data for 14 ancient individuals from Efate and Epi Islands in Vanuatu from 2900-150 BP, as well as 185 present-day individuals from 18 islands. We find that people of almost entirely Papuan ancestry arrived in Vanuatu by around 2300 BP, most likely reflecting migrations a few hundred years earlier at the end of the Lapita period, when there is also evidence of changes in skeletal morphology and cessation of long-distance trade between Near and Remote Oceania [2, 3]. Papuan ancestry was subsequently diluted through admixture but remains at least 80%-90% in most islands. Through a fine-grained analysis of ancestry profiles, we show that the Papuan ancestry in Vanuatu derives from the Bismarck Archipelago rather than the geographically closer Solomon Islands. However, the Papuan ancestry in Polynesia-the most remote Pacific islands-derives from different sources, documenting a third stream of migration from Near to Remote Oceania.

Posth C, Nägele K, Colleran H, Valentin F, Bedford S, Kami KW, Shing R, Buckley H, Kinaston R, Walworth M et al. 2018. Language continuity despite population replacement in Remote Oceania. Nat Ecol Evol, 2 (4), pp. 731-740. | Citations: 2 (Scopus) | Show Abstract | Read more

Recent genomic analyses show that the earliest peoples reaching Remote Oceania-associated with Austronesian-speaking Lapita culture-were almost completely East Asian, without detectable Papuan ancestry. However, Papuan-related genetic ancestry is found across present-day Pacific populations, indicating that peoples from Near Oceania have played a significant, but largely unknown, ancestral role. Here, new genome-wide data from 19 ancient South Pacific individuals provide direct evidence of a so-far undescribed Papuan expansion into Remote Oceania starting ~2,500 yr BP, far earlier than previously estimated and supporting a model from historical linguistics. New genome-wide data from 27 contemporary ni-Vanuatu demonstrate a subsequent and almost complete replacement of Lapita-Austronesian by Near Oceanian ancestry. Despite this massive demographic change, incoming Papuan languages did not replace Austronesian languages. Population replacement with language continuity is extremely rare-if not unprecedented-in human history. Our analyses show that rather than one large-scale event, the process was incremental and complex, with repeated migrations and sex-biased admixture with peoples from the Bismarck Archipelago.

Johnston JJ, van der Smagt JJ, Rosenfeld JA, Pagnamenta AT, Alswaid A, Baker EH, Blair E, Borck G, Brinkmann J, Craigen W et al. 2018. Autosomal recessive Noonan syndrome associated with biallelic LZTR1 variants. Genet Med, | Citations: 3 (European Pubmed Central) | Show Abstract | Read more

PurposeTo characterize the molecular genetics of autosomal recessive Noonan syndrome.MethodsFamilies underwent phenotyping for features of Noonan syndrome in children and their parents. Two multiplex families underwent linkage analysis. Exome, genome, or multigene panel sequencing was used to identify variants. The molecular consequences of observed splice variants were evaluated by reverse-transcription polymerase chain reaction.ResultsTwelve families with a total of 23 affected children with features of Noonan syndrome were evaluated. The phenotypic range included mildly affected patients, but it was lethal in some, with cardiac disease and leukemia. All of the parents were unaffected. Linkage analysis using a recessive model supported a candidate region in chromosome 22q11, which includes LZTR1, previously shown to harbor mutations in patients with Noonan syndrome inherited in a dominant pattern. Sequencing analyses of 21 live-born patients and a stillbirth identified biallelic pathogenic variants in LZTR1, including putative loss-of-function, missense, and canonical and noncanonical splicing variants in the affected children, with heterozygous, clinically unaffected parents and heterozygous or normal genotypes in unaffected siblings.ConclusionThese clinical and genetic data confirm the existence of a form of Noonan syndrome that is inherited in an autosomal recessive pattern and identify biallelic mutations in LZTR1.Genet Med advance online publication, 22 February 2018; doi:10.1038/gim.2017.249.

Bliss CM, Bowyer G, Anagnostou NA, Havelock T, Snudden CM, Davies H, de Cassan SC, Grobbelaar A, Lawrie AM, Venkatraman N et al. 2018. Assessment of novel vaccination regimens using viral vectored liver stage malaria vaccines encoding ME-TRAP. Sci Rep, 8 (1), pp. 3390. | Citations: 2 (European Pubmed Central) | Show Abstract | Read more

Heterologous prime-boost vaccination with viral vectors simian adenovirus 63 (ChAd63) and Modified Vaccinia Ankara (MVA) induces potent T cell and antibody responses in humans. The 8-week regimen demonstrates significant efficacy against malaria when expressing the pre-erythrocytic malaria antigen Thrombospondin-Related Adhesion Protein fused to a multiple epitope string (ME-TRAP). We tested these vaccines in 7 new 4- and 8- week interval schedules to evaluate safety and immunogenicity of multiple ChAd63 ME-TRAP priming vaccinations (denoted A), multiple MVA ME-TRAP boosts (denoted M) and alternating vectors. All regimens exhibited acceptable reactogenicity and CD8+ T cell immunogenicity was enhanced with a 4-week interval (AM) and with incorporation of additional ChAd63 ME-TRAP vaccination at 4- or 8-weeks (AAM or A_A_M). Induction of TRAP antibodies was comparable between schedules. T cell immunity against the ChAd63 hexon did not affect T cell responses to the vaccine insert, however pre-vaccination ChAd63-specific T cells correlated with reduced TRAP antibodies. Vaccine-induced antibodies against MVA did not affect TRAP antibody induction, and correlated positively with ME-TRAP-specific T cells. This study identifies potentially more effective immunisation regimens to assess in Phase IIa trials and demonstrates a degree of flexibility with the timing of vectored vaccine administration, aiding incorporation into existing vaccination programmes.

Agrawal S, Tapmeier T, Rahmioglu N, Kirtley S, Zondervan K, Becker C. 2018. The miRNA Mirage: How Close Are We to Finding a Non-Invasive Diagnostic Biomarker in Endometriosis? A Systematic Review. Int J Mol Sci, 19 (2), pp. 599-599. | Citations: 1 (European Pubmed Central) | Show Abstract | Read more

BACKGROUND: Endometriosis is a common disorder of the reproductive age group, characterised by the presence of ectopic endometrial tissue. The disease not only causes enormous suffering to the affected women, but also brings a tremendous medical and economic burden to bear on society. There is a long lag phase between the onset and diagnosis of the disease, mainly due to its non-specific symptoms and the lack of a non-invasive test. Endometriosis can only be diagnosed invasively by laparoscopy. A specific, non-invasive test to diagnose endometriosis is an unmet clinical need. The recent discovery of microRNAs (miRNAs) as modulators of gene expression, and their stability and specificity, make them an attractive candidate biomarker. Various studies on miRNAs in endometriosis have identified their cardinal role in the pathogenesis of the disease, and have proposed them as potential biomarkers in endometriosis. Rationale/Objectives: The aims of this review were to study the role of circulatory miRNAs in endometriosis, and bring to light whether circulatory miRNAs could be potential non-invasive biomarkers to diagnose the disease. SEARCH METHODS: Three databases, PubMed, EMBASE, and BIOSIS were searched, using a combination of Mesh or Emtree headings and free-text terms, to identify literature relating to circulating miRNAs in endometriosis published from 1996 to 31 December 2017. Only peer-reviewed, full-text original research articles in English were included in the current review. The studies meeting the inclusion criteria were critically assessed and checked using the QUADAS-2 (Quality Assessment of Diagnostic Accuracy Studies) tool. The dysregulated miRNAs were assessed regarding the concordance between the various studies and their role in the disease. OUTCOMES: Nine studies were critically analysed, and 42 different miRNAs were found to be dysregulated in them, with only one common miRNA (miR-20a) differentially expressed in more than one study. miR-17-5p/20a, miR-200, miR-199a, miR-143, and miR-145 were explored for their pivotal role in the aetiopathogenesis of endometriosis. Wider implications: It is emerging that miRNAs play a central role in the pathogenesis of endometriosis and have the potential of being promising biomarkers. Circulating miRNAs as a non-invasive diagnostic tool may shorten the delay in the diagnosis of the disease, thus alleviating the suffering of women and reducing the burden on health care systems. However, despite numerous studies on circulating miRNAs in endometriosis, no single miRNA or any panel of them seems to meet the criteria of a diagnostic biomarker. The disagreement between the various studies upholds the demand of larger, well-controlled systematic validation studies with uniformity in the research approaches and involving diverse populations.

Total publications on this page: 100

Total citations for publications on this page: 96