Recent publications

Lenglet M, Robriquet F, Schwarz K, Camps C, Couturier A, Hoogewijs D, Buffet A, Knight SJ, Gad S, Couvé S et al. 2018. New lessons from an old gene: complex splicing and a novel cryptic exon in VHL gene cause erythrocytosis and VHL disease. Blood, pp. blood-2018-03-838235-blood-2018-03-838235. | 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 von Hippel-Lindau disease, which is characterized by the development of highly vascularized tumors. Here, we identified a new VHL cryptic-exon (termed E1') deep in intron 1 that is naturally expressed in many tissues. More importantly, we identified mutations in E1' in seven families with erythrocytosis (one homozygous case and six compound-heterozygous cases with a mutation in E1' in addition to a mutation in VHL coding sequences) and in one large family with typical VHL disease but without any alteration in the other VHL exons. In this study we have shown that the mutations induced a dysregulation of the VHL splicing with excessive retention of E1' and are associated with a downregulation of VHL protein expression. In addition, we have demonstrated a pathogenic role for synonymous mutations in VHL-Exon 2 that alter splicing through E2-skipping in five families with erythrocytosis or VHL disease. In all the studied cases, the mutations differentially impact splicing, correlating with phenotype severity. This study demonstrates that cryptic-exon-retention or 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 into the VHL-related-hypoxia-signaling pathway.

Batty EM, Chaemchuen S, Blacksell S, Richards A, 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.

Ternette N, Olde Nordkamp MJM, Muller 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, pp. e1700465. | 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 utilising 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 analysed the immunopeptidomes of 6 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 tumour and adjacent normal tissue and 130 of the peptides were found to have higher abundance in tumour than in normal tissues. To determine potential therapeutic target proteins, we calculated the average tumour-associated cohort coverage (aTaCC) that represents the percentage coverage of each protein in this cohort by peptides that had higher tumourial 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 the wide spectrum of target discovery. This article is protected by copyright. All rights reserved.

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

Todd JA. 2018. Evidence that UBASH3 is a causal gene for type 1 diabetes. Eur J Hum Genet, pp. 1-3. | Read more

Osgood JA, Knight JC. 2018. Translating GWAS in rheumatic disease: approaches to establishing mechanism and function for genetic associations with ankylosing spondylitis Briefings in Functional Genomics, | Read more

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.

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.

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.

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, | 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. | 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, pp. ASN.2017111179-ASN.2017111179. | 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), pp. e00368-18-e00368-18. | 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. | 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. | 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.

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 (Web of Science Lite) | 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. | 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, | 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 (Web of Science Lite) | 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, | 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: 1 (Scopus) | 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. | 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. | 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, | 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 MHC. 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 MHC 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 HLA genes using conventional mapping can be corrected using AltHapAlignR to allow more accurate quantification of gene expression for individual alleles and haplotypes. Availability: Source code freely available at https://github.com/jknightlab/AltHapAlignR. Contact: peter.humburg@gmail.com or julian@well.ox.ac.uk. Supplementary information: 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.

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

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), | 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.

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

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

Schwarze K, Buchanan J, Taylor JC, Wordsworth S. 2018. Are whole-exome and whole-genome sequencing approaches cost-effective? A systematic review of the literature. Genet Med, | Show Abstract | Read more

PurposeWe conducted a systematic literature review to summarize the current health economic evidence for whole-exome sequencing (WES) and whole-genome sequencing (WGS).MethodsRelevant studies were identified in the EMBASE, MEDLINE, Cochrane Library, EconLit and University of York Centre for Reviews and Dissemination databases from January 2005 to July 2016. Publications were included in the review if they were economic evaluations, cost studies, or outcome studies.ResultsThirty-six studies met our inclusion criteria. These publications investigated the use of WES and WGS in a variety of genetic conditions in clinical practice, the most common being neurological or neurodevelopmental disorders. Study sample size varied from a single child to 2,000 patients. Cost estimates for a single test ranged from $555 to $5,169 for WES and from $1,906 to $24,810 for WGS. Few cost analyses presented data transparently and many publications did not state which components were included in cost estimates.ConclusionThe current health economic evidence base to support the more widespread use of WES and WGS in clinical practice is very limited. Studies that carefully evaluate the costs, effectiveness, and cost-effectiveness of these tests are urgently needed to support their translation into clinical practice.GENETICS in MEDICINE advance online publication, 15 February 2018; doi:10.1038/gim.2017.247.

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

BACKGROUND: Colorectal cancer (CRC) screening might be improved by using a measure of prior risk to modulate screening intensity or the faecal immunochemical test threshold. Intermediate molecular biomarkers could aid risk prediction by capturing both known and unknown risk factors. METHODS: We sampled normal bowel mucosa from the proximal colon, distal colon and rectum of 317 individuals undergoing colonoscopy. We defined cases as having a personal history of colorectal polyp(s)/cancer, and controls as having no history of colorectal neoplasia. Molecular analyses were performed for: telomere length (TL); global methylation; and the expression of genes in molecular pathways associated with colorectal tumourigenesis. We also calculated a polygenic risk score (PRS) based on CRC susceptibility polymorphisms. RESULTS: Bowel TL was significantly longer in cases than controls, but was not associated with blood TL. PRS was significantly and independently higher in cases. Hypermethylation showed a suggestive association with case:control status. No gene or pathway was differentially expressed between cases and controls. Gene expression often varied considerably between bowel locations. CONCLUSIONS: PRS and bowel TL (but not blood TL) may be clinically-useful predictors of CRC risk. Sample collection to assess these biomarkers is feasible in clinical practice, especially where population screening uses flexible sigmoidoscopy or colonoscopy.

Sung YJ, Winkler TW, de Las Fuentes L, Bentley AR, Brown MR, Kraja AT, Schwander K, Ntalla I, Guo X, Franceschini N et al. 2018. A Large-Scale Multi-ancestry Genome-wide Study Accounting for Smoking Behavior Identifies Multiple Significant Loci for Blood Pressure. Am J Hum Genet, 102 (3), pp. 375-400. | Citations: 1 (Web of Science Lite) | Show Abstract | Read more

Genome-wide association analysis advanced understanding of blood pressure (BP), a major risk factor for vascular conditions such as coronary heart disease and stroke. Accounting for smoking behavior may help identify BP loci and extend our knowledge of its genetic architecture. We performed genome-wide association meta-analyses of systolic and diastolic BP incorporating gene-smoking interactions in 610,091 individuals. Stage 1 analysis examined ∼18.8 million SNPs and small insertion/deletion variants in 129,913 individuals from four ancestries (European, African, Asian, and Hispanic) with follow-up analysis of promising variants in 480,178 additional individuals from five ancestries. We identified 15 loci that were genome-wide significant (p < 5 × 10-8) in stage 1 and formally replicated in stage 2. A combined stage 1 and 2 meta-analysis identified 66 additional genome-wide significant loci (13, 35, and 18 loci in European, African, and trans-ancestry, respectively). A total of 56 known BP loci were also identified by our results (p < 5 × 10-8). Of the newly identified loci, ten showed significant interaction with smoking status, but none of them were replicated in stage 2. Several loci were identified in African ancestry, highlighting the importance of genetic studies in diverse populations. The identified loci show strong evidence for regulatory features and support shared pathophysiology with cardiometabolic and addiction traits. They also highlight a role in BP regulation for biological candidates such as modulators of vascular structure and function (CDKN1B, BCAR1-CFDP1, PXDN, EEA1), ciliopathies (SDCCAG8, RPGRIP1L), telomere maintenance (TNKS, PINX1, AKTIP), and central dopaminergic signaling (MSRA, EBF2).

Smith D, Magri A, Bonsall D, Ip CLC, Trebes A, Brown A, Piazza P, Bowden R, Nguyen D, Ansari MA et al. 2018. Resistance analysis of genotype 3 hepatitis C virus indicates subtypes inherently resistant to nonstructural protein 5A inhibitors. Hepatology, | Show Abstract | Read more

Hepatitis C virus (HCV) genotype (gt) 3 is highly prevalent globally, with non-gt3a subtypes common in Southeast Asia. Resistance-associated substitutions (RASs) have been shown to play a role in treatment failure. However, the role of RASs in gt3 is not well understood. We report the prevalence of RASs in a cohort of direct-acting antiviral treatment-naive, gt3-infected patients, including those with rarer subtypes, and evaluate the effect of these RASs on direct-acting antivirals in vitro. Baseline samples from 496 gt3 patients enrolled in the BOSON clinical trial were analyzed by next-generation sequencing after probe-based enrichment for HCV. Whole viral genomes were analyzed for the presence of RASs to approved direct-acting antivirals. The resistance phenotype of RASs in combination with daclatasvir, velpatasvir, pibrentasvir, elbasvir, and sofosbuvir was measured using the S52 ΔN gt3a replicon model. The nonstructural protein 5A A30K and Y93H substitutions were the most common at 8.9% (n = 44) and 12.3% (n = 61), respectively, and showed a 10-fold and 11-fold increase in 50% effect concentration for daclatasvir compared to the unmodified replicon. Paired RASs (A30K + L31M and A30K + Y93H) were identified in 18 patients (9 of each pair); these combinations were shown to be highly resistant to daclatasvir, velpatasvir, elbasvir, and pibrentasvir. The A30K + L31M combination was found in all gt3b and gt3g samples. CONCLUSION: Our study reveals high frequencies of RASs to nonstructural protein 5A inhibitors in gt3 HCV; the paired A30K + L31M substitutions occur in all patients with gt3b and gt3g virus, and in vitro analysis suggests that these subtypes may be inherently resistant to all approved nonstructural protein 5A inhibitors for gt3 HCV. (Hepatology 2018).

Vecellio M, Cortes A, Roberts AR, Ellis J, Cohen CJ, Knight JC, Brown MA, Bowness P, Wordsworth BP. 2018. Evidence for a second ankylosing spondylitis-associated RUNX3 regulatory polymorphism. RMD Open, 4 (1), pp. e000628. | Show Abstract | Read more

Objectives: To explore the functions of RUNX3 single nucleotide polymorphisms (SNPs) associated with ankylosing spondylitis (AS). Methods: Individual SNP associations were evaluated in 4230 UK cases. Their effects on transcription factor (TF) binding, transcription regulation, chromatin modifications, gene expression and gene interactions were tested by database interrogation, luciferase reporter assays, electrophoretic mobility gel shifts, chromatin immunoprecipitation and real-time PCR. Results: We confirmed the independent association of AS with rs4265380, which was robust (P=4.7×10-6) to conditioning on another nearby AS-associated RUNX3 SNP (rs4648889). A RUNX3 haplotype incorporating both SNPs was strongly associated with AS (OR 6.2; 95% CI 3.1 to 13.2, P=1.4×10-8). In a large UK cohort, rs4265380 is associated with leucocyte counts (including monocytes). RUNX3 expression is lower in AS peripheral blood mononuclear cells than healthy controls (P<0.002), independent of rs4265380 genotype. Enhancer function for this RUNX3 region was suggested by increased luciferase activity (approximately tenfold; P=0.005) for reporter constructs containing rs4265380. In monocytes, there was differential allelic binding of nuclear protein extracts to a 50 bp DNA probe containing rs4265380 that was strongly augmented by lipopolysaccharide activation. TF binding also included the histone modifier p300. There was enrichment for histone modifications associated with active enhancer elements (H3K27Ac and H3K79Me2) that may be allele dependent. Hi-C database interrogation showed chromosome interactions of RUNX3 bait with the nearby RP4-799D16.1 lincRNA. Conclusions: The association of AS with this RUNX3 regulatory region involves at least two SNPs apparently operating in different cell types. Monocytes may be potential therapeutic targets in AS.

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

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

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

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

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

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

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

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

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

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

Tao C-L, Liu Y-T, Sun R, Zhang B, Qi L, Shivakoti S, Tian C-L, Zhang P, Lau P-M, Zhou ZH, Bi G-Q. 2018. Differentiation and Characterization of Excitatory and Inhibitory Synapses by Cryo-electron Tomography and Correlative Microscopy. J Neurosci, 38 (6), pp. 1493-1510. | Show Abstract | Read more

As key functional units in neural circuits, different types of neuronal synapses play distinct roles in brain information processing, learning, and memory. Synaptic abnormalities are believed to underlie various neurological and psychiatric disorders. Here, by combining cryo-electron tomography and cryo-correlative light and electron microscopy, we distinguished intact excitatory and inhibitory synapses of cultured hippocampal neurons, and visualized the in situ 3D organization of synaptic organelles and macromolecules in their native state. Quantitative analyses of >100 synaptic tomograms reveal that excitatory synapses contain a mesh-like postsynaptic density (PSD) with thickness ranging from 20 to 50 nm. In contrast, the PSD in inhibitory synapses assumes a thin sheet-like structure ∼12 nm from the postsynaptic membrane. On the presynaptic side, spherical synaptic vesicles (SVs) of 25-60 nm diameter and discus-shaped ellipsoidal SVs of various sizes coexist in both synaptic types, with more ellipsoidal ones in inhibitory synapses. High-resolution tomograms obtained using a Volta phase plate and electron filtering and counting reveal glutamate receptor-like and GABAA receptor-like structures that interact with putative scaffolding and adhesion molecules, reflecting details of receptor anchoring and PSD organization. These results provide an updated view of the ultrastructure of excitatory and inhibitory synapses, and demonstrate the potential of our approach to gain insight into the organizational principles of cellular architecture underlying distinct synaptic functions.SIGNIFICANCE STATEMENT To understand functional properties of neuronal synapses, it is desirable to analyze their structure at molecular resolution. We have developed an integrative approach combining cryo-electron tomography and correlative fluorescence microscopy to visualize 3D ultrastructural features of intact excitatory and inhibitory synapses in their native state. Our approach shows that inhibitory synapses contain uniform thin sheet-like postsynaptic densities (PSDs), while excitatory synapses contain previously known mesh-like PSDs. We discovered "discus-shaped" ellipsoidal synaptic vesicles, and their distributions along with regular spherical vesicles in synaptic types are characterized. High-resolution tomograms further allowed identification of putative neurotransmitter receptors and their heterogeneous interaction with synaptic scaffolding proteins. The specificity and resolution of our approach enables precise in situ analysis of ultrastructural organization underlying distinct synaptic functions.

Cahill TJ, Wijesurendra RS, Casadei B. 2018. Rivaroxaban in Stable Cardiovascular Disease. N Engl J Med, 378 (4), pp. 397. | Read more

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

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

Biswas S, Quante M, Leedham S, Jansen M. 2018. The metaplastic mosaic of Barrett's oesophagus. Virchows Arch, 472 (1), pp. 43-54. | Show Abstract | Read more

Barrett's oesophagus surveillance biopsies represent a significant share of the daily workload for a busy histopathology department. Given the emphasis on endoscopic detection and dysplasia grading, it is easy to forget that the benefits of these screening programs remain unproven. The majority of patients are at low risk of progression to oesophageal adenocarcinoma, and periodic surveillance of these patients is burdensome and costly. Here, we investigate the parallels in the development of Barrett's oesophagus and other scenarios of wound healing in the intestine. There is now increased recognition of the full range of glandular phenotypes that can be found in patients' surveillance biopsies, and emerging evidence suggests parallel pathways to oesophageal adenocarcinoma. Greater understanding of the conditions that favour progression to cancer in the distal oesophagus will allow us to focus resources on patients at increased risk.

Mitchell TJ, Turajlic S, Rowan A, Nicol D, Farmery JHR, O'Brien T, Martincorena I, Tarpey P, Angelopoulos N, Yates LR et al. 2018. Timing the Landmark Events in the Evolution of Clear Cell Renal Cell Cancer: TRACERx Renal. Cell, 173 (3), pp. 611-623.e17. | Citations: 3 (Web of Science Lite) | Show Abstract | Read more

Clear cell renal cell carcinoma (ccRCC) is characterized by near-universal loss of the short arm of chromosome 3, deleting several tumor suppressor genes. We analyzed whole genomes from 95 biopsies across 33 patients with clear cell renal cell carcinoma. We find hotspots of point mutations in the 5' UTR of TERT, targeting a MYC-MAX-MAD1 repressor associated with telomere lengthening. The most common structural abnormality generates simultaneous 3p loss and 5q gain (36% patients), typically through chromothripsis. This event occurs in childhood or adolescence, generally as the initiating event that precedes emergence of the tumor's most recent common ancestor by years to decades. Similar genomic changes drive inherited ccRCC. Modeling differences in age incidence between inherited and sporadic cancers suggests that the number of cells with 3p loss capable of initiating sporadic tumors is no more than a few hundred. Early development of ccRCC follows well-defined evolutionary trajectories, offering opportunity for early intervention.

Gloyn AL, Drucker DJ. 2018. Precision medicine in the management of type 2 diabetes The Lancet Diabetes and Endocrinology, | Show Abstract | Read more

© 2018 Elsevier Ltd The study of type 2 diabetes has been driven by advances in human genetics, epigenetics, biomarkers, mechanistic studies, and large clinical trials, enabling new insights into disease susceptibility, pathophysiology, progression, and development of complications. Simultaneously, several new drug classes with different mechanisms of action have been introduced over the past two decades, accompanied by data about cardiovascular safety and non-glycaemic outcomes. In this Review, we critically examine the progress and integration of this new science into clinical practice, and review opportunities for enabling the use of precision medicine in the diagnosis and treatment of type 2 diabetes. We contrast the success in delivering personalised medicine for monogenic diabetes with the greater challenge of providing a precision medicine approach for type 2 diabetes, highlighting gaps, limitations, and areas requiring further study.

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

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

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

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

Muszkiewicz A, Liu X, Bueno-Orovio A, Lawson BAJ, Burrage K, Casadei B, Rodriguez B. 2018. From ionic to cellular variability in human atrial myocytes: an integrative computational and experimental study. Am J Physiol Heart Circ Physiol, 314 (5), pp. H895-H916. | Citations: 1 (Scopus) | Show Abstract | Read more

Variability refers to differences in physiological function between individuals, which may translate into different disease susceptibility and treatment efficacy. Experiments in human cardiomyocytes face wide variability and restricted tissue access; under these conditions, computational models are a useful complementary tool. We conducted a computational and experimental investigation in cardiomyocytes isolated from samples of the right atrial appendage of patients undergoing cardiac surgery to evaluate the impact of variability in action potentials (APs) and subcellular ionic densities on Ca2+ transient dynamics. Results showed that 1) variability in APs and ionic densities is large, even within an apparently homogenous patient cohort, and translates into ±100% variation in ionic conductances; 2) experimentally calibrated populations of models with wide variations in ionic densities yield APs overlapping with those obtained experimentally, even if AP characteristics of the original generic model differed significantly from experimental APs; 3) model calibration with AP recordings restricts the variability in ionic densities affecting upstroke and resting potential, but redundancy in repolarization currents admits substantial variability in ionic densities; and 4) model populations constrained with experimental APs and ionic densities exhibit three Ca2+ transient phenotypes, differing in intracellular Ca2+ handling and Na+/Ca2+ membrane extrusion. These findings advance our understanding of the impact of variability in human atrial electrophysiology. NEW & NOTEWORTHY Variability in human atrial electrophysiology is investigated by integrating for the first time cellular-level and ion channel recordings in computational electrophysiological models. Ion channel calibration restricts current densities but not cellular phenotypic variability. Reduced Na+/Ca2+ exchanger is identified as a primary mechanism underlying diastolic Ca2+ fluctuations in human atrial myocytes.

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

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

Faller KME, Atzler D, McAndrew DJ, Zervou S, Whittington HJ, Simon JN, Aksentijevic D, Ten Hove M, Choe C-U, Isbrandt D et al. 2018. Impaired cardiac contractile function in arginine:glycine amidinotransferase knockout mice devoid of creatine is rescued by homoarginine but not creatine. Cardiovasc Res, 114 (3), pp. 417-430. | Show Abstract | Read more

Aims: Creatine buffers cellular adenosine triphosphate (ATP) via the creatine kinase reaction. Creatine levels are reduced in heart failure, but their contribution to pathophysiology is unclear. Arginine:glycine amidinotransferase (AGAT) in the kidney catalyses both the first step in creatine biosynthesis as well as homoarginine (HA) synthesis. AGAT-/- mice fed a creatine-free diet have a whole body creatine-deficiency. We hypothesized that AGAT-/- mice would develop cardiac dysfunction and rescue by dietary creatine would imply causality. Methods and results: Withdrawal of dietary creatine in AGAT-/- mice provided an estimate of myocardial creatine efflux of ∼2.7%/day; however, in vivo cardiac function was maintained despite low levels of myocardial creatine. Using AGAT-/- mice naïve to dietary creatine we confirmed absence of phosphocreatine in the heart, but crucially, ATP levels were unchanged. Potential compensatory adaptations were absent, AMPK was not activated and respiration in isolated mitochondria was normal. AGAT-/- mice had rescuable changes in body water and organ weights suggesting a role for creatine as a compatible osmolyte. Creatine-naïve AGAT-/- mice had haemodynamic impairment with low LV systolic pressure and reduced inotropy, lusitropy, and contractile reserve. Creatine supplementation only corrected systolic pressure despite normalization of myocardial creatine. AGAT-/- mice had low plasma HA and supplementation completely rescued all other haemodynamic parameters. Contractile dysfunction in AGAT-/- was confirmed in Langendorff perfused hearts and in creatine-replete isolated cardiomyocytes, indicating that HA is necessary for normal cardiac function. Conclusions: Our findings argue against low myocardial creatine per se as a major contributor to cardiac dysfunction. Conversely, we show that HA deficiency can impair cardiac function, which may explain why low HA is an independent risk factor for multiple cardiovascular diseases.

Tyrrell HEJ, Church DN, Joseph J, Traill ZC, Sullivan ME, Tuthill MH, Verrill CL, Pintus EP, Dallas NL, Rogers PB et al. 2018. Changing Practice Evaluation-Stage 1 Seminoma: Outcomes With Adjuvant Treatment Versus Surveillance: Risk Factors for Recurrence and Optimizing Follow-up Protocols-Experience From a Supraregional Center. Clin Genitourin Cancer, 16 (3), pp. 240-244. | Show Abstract | Read more

BACKGROUND: Stage 1 seminoma is frequently cured by radical orchiectomy; however, the management strategies after this diagnosis vary in terms of the use of adjuvant treatment and the nature of the follow-up protocols. We analyzed stage 1 seminomas treated in the Thames Valley Cancer Network for outcomes to determine whether any factors are predictive of recurrence. We also studied relapses to determine the optimal follow-up schedule and protocol. MATERIALS AND METHODS: Data were obtained from centers within the Thames Valley Cancer Network for a 12-year period from 2004 to 2016. We identified 501 patients with stage 1 seminoma. RESULTS: Relapses occurred in 6.2% of the patients receiving adjuvant treatment and 6.1% of those who did not. The only statistically significant predictive factor identified for relapse was rete testis invasion, and the risk was greater when only stromal rete invasion was included, rather than pagetoid as well. A trend was seen toward an increased risk with increased tumor size, but the difference was not statistically significant. Recurrences developed within the first 2 years after surgery in nearly 75% of cases and were identified through surveillance computed tomography scans in 54.8% of the patients. All relapses were treated curatively. CONCLUSION: Active surveillance leads to excellent outcomes for stage 1 seminoma; however, adjuvant treatment should be reserved for those with high-risk disease. Follow-up schedules should include computed tomography imaging during the first 3 years, long-term measurement of tumor markers, and mechanisms for patients to be seen promptly should symptoms of tumor recurrence occur.

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

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

Kindt ASD, Fuerst RW, Knoop J, Laimighofer M, Telieps T, Hippich M, Woerheide MA, Wahl S, Wilson R, Sedlmeier E-M et al. 2018. Allele-specific methylation of type 1 diabetes susceptibility genes. J Autoimmun, 89 pp. 63-74. | Citations: 1 (Scopus) | Show Abstract | Read more

The susceptibility to autoimmune diseases is influenced by genes encoding major histocompatibility complex (MHC) proteins. By examining the epigenetic methylation maps of cord blood samples, we found marked differences in the methylation status of CpG sites within the MHC genes (cis-metQTLs) between carriers of the type 1 diabetes risk haplotypes HLA-DRB1*03-DQA1*0501-DQB1*0201 (DR3-DQ2) and HLA-DRB1*04-DQA1*0301-DQB1*0302 (DR4-DQ8). These differences were found in children and adults, and were accompanied by reduced HLA-DR protein expression in immune cells with the HLA-DR3-DQ2 haplotype. Extensive cis-metQTLs were identified in all 45 immune and non-immune type 1 diabetes susceptibility genes analyzed in this study. We observed and validated a novel association between the methylation status of CpG sites within the LDHC gene and the development of insulin autoantibodies in early childhood in children who are carriers of the highest type 1 diabetes risk genotype. Functionally relevant epigenetic changes in susceptibility genes may represent therapeutic targets for type 1 diabetes.

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

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

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

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

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

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

Delgado V, Di Biase L, Leung M, Romero J, Tops LF, Casadei B, Marrouche N, Bax JJ. 2017. Structure and Function of the Left Atrium and Left Atrial Appendage: AF and Stroke Implications. J Am Coll Cardiol, 70 (25), pp. 3157-3172. | Citations: 2 (Scopus) | Show Abstract | Read more

Atrial fibrillation (AF) and stroke are important major health problems that share common risk factors and frequently coexist. Left atrial (LA) remodeling is an important underlying substrate for AF and stroke. LA dilation and dysfunction form a prothrombotic milieu characterized by blood stasis and endothelial dysfunction. In addition, alterations of the atrial cardiomyocytes, increase of noncollagen deposits in the interstitial space and fibrosis, favor the occurrence of re-entry that predisposes to AF. Eventually, AF further impairs LA function and promotes LA remodeling, closing a self-perpetuating vicious circle. Multimodality imaging provides a comprehensive evaluation of several aspects of LA remodeling and offers several parameters to identify patients at risk of AF and stroke. How multimodality imaging can be integrated in clinical management of patients at risk of AF and stroke is the focus of the present review paper.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Schwerd T, Bryant RV, Pandey S, Capitani M, Meran L, Cazier J-B, Jung J, Mondal K, Parkes M, Mathew CG et al. 2018. NOX1 loss-of-function genetic variants in patients with inflammatory bowel disease. Mucosal Immunol, 11 (2), pp. 562-574. | Citations: 2 (Scopus) | Show Abstract | Read more

Genetic defects that affect intestinal epithelial barrier function can present with very early-onset inflammatory bowel disease (VEOIBD). Using whole-genome sequencing, a novel hemizygous defect in NOX1 encoding NAPDH oxidase 1 was identified in a patient with ulcerative colitis-like VEOIBD. Exome screening of 1,878 pediatric patients identified further seven male inflammatory bowel disease (IBD) patients with rare NOX1 mutations. Loss-of-function was validated in p.N122H and p.T497A, and to a lesser degree in p.Y470H, p.R287Q, p.I67M, p.Q293R as well as the previously described p.P330S, and the common NOX1 SNP p.D360N (rs34688635) variant. The missense mutation p.N122H abrogated reactive oxygen species (ROS) production in cell lines, ex vivo colonic explants, and patient-derived colonic organoid cultures. Within colonic crypts, NOX1 constitutively generates a high level of ROS in the crypt lumen. Analysis of 9,513 controls and 11,140 IBD patients of non-Jewish European ancestry did not reveal an association between p.D360N and IBD. Our data suggest that loss-of-function variants in NOX1 do not cause a Mendelian disorder of high penetrance but are a context-specific modifier. Our results implicate that variants in NOX1 change brush border ROS within colonic crypts at the interface between the epithelium and luminal microbes.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Lashin HMS, Nadkarni S, Oggero S, Jones HR, Knight JC, Hinds CJ, Perretti M. 2018. Microvesicle Subsets in Sepsis Due to Community Acquired Pneumonia Compared to Faecal Peritonitis. Shock, 49 (4), pp. 393-401. | Citations: 2 (Scopus) | Show Abstract | Read more

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

Dendrou CA, Fugger L. 2017. Immunomodulation in multiple sclerosis: promises and pitfalls. Curr Opin Immunol, 49 pp. 37-43. | Citations: 4 (Web of Science Lite) | Show Abstract | Read more

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

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

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

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

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

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

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

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

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

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

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

Total publications on this page: 100

Total citations for publications on this page: 139