School of Mathematics and Statistics - Research Publications

Permanent URI for this collection

Search Results

Now showing 1 - 10 of 15
  • Item
    Thumbnail Image
    Identification and Re-consent of Existing Cord Blood Donors for Creation of Induced Pluripotent Stem Cell Lines for Potential Clinical Applications
    Abberton, KM ; McDonald, TL ; Diviney, M ; Holdsworth, R ; Leslie, S ; Delatycki, MB ; Liu, L ; Klamer, G ; Johnson, P ; Elwood, NJ (OXFORD UNIV PRESS, 2022-10-21)
    We aim to create a bank of clinical grade cord blood-derived induced pluripotent stem cell lines in order to facilitate clinical research leading to the development of new cellular therapies. Here we present a clear pathway toward the creation of such a resource, within a strong quality framework, and with the appropriate regulatory, government and ethics approvals, along with a dynamic follow-up and re-consent process of cord blood donors from the public BMDI Cord Blood Bank. Interrogation of the cord blood bank inventory and next generation sequencing was used to identify and confirm 18 donors with suitable HLA homozygous haplotypes. Regulatory challenges that may affect global acceptance of the cell lines, along with the quality standards required to operate as part of a global network, are being met by working in collaboration with bodies such as the International Stem Cell Banking Initiative (ISCBI) and the Global Alliance for iPSC Therapies (GAiT). Ethics approval was granted by an Institutional Human Research Ethics Committee, and government approval has been obtained to use banked cord blood for this purpose. New issues of whole-genome sequencing and the relevant donor safeguards and protections were considered with input from clinical genetics services, including the rights and information flow to donors, and commercialization aspects. The success of these processes has confirmed feasibility and utility of using banked cord blood to produce clinical-grade iPSC lines for potential cellular therapies.
  • Item
    Thumbnail Image
    Allele imputation for the killer cell immunoglobulin-like receptor KIR3DL1/S1
    Harrison, GR ; Leaton, LAA ; Harrison, E ; Kichula, K ; Viken, MJ ; Shortt, J ; Gignoux, C ; Lie, B ; Vukcevic, D ; Leslie, S ; Norman, P ; Kouyos, RD (PUBLIC LIBRARY SCIENCE, 2022-02)
    Highly polymorphic interaction of KIR3DL1 and KIR3DS1 with HLA class I ligands modulates the effector functions of natural killer (NK) cells and some T cells. This genetically determined diversity affects severity of infections, immune-mediated diseases, and some cancers, and impacts the course of immunotherapies, including transplantation. KIR3DL1 is an inhibitory receptor, and KIR3DS1 is an activating receptor encoded by the KIR3DL1/S1 gene that has more than 200 diverse and divergent alleles. Determination of KIR3DL1/S1 genotypes for medical application is hampered by complex sequence and structural variation, requiring targeted approaches to generate and analyze high-resolution allele data. To overcome these obstacles, we developed and optimized a model for imputing KIR3DL1/S1 alleles at high-resolution from whole-genome SNP data. We designed the model to represent a substantial component of human genetic diversity. Our Global imputation model is effective at genotyping KIR3DL1/S1 alleles with an accuracy ranging from 88% in Africans to 97% in East Asians, with mean specificity of 99% and sensitivity of 95% for alleles >1% frequency. We used the established algorithm of the HIBAG program, in a modification named Pulling Out Natural killer cell Genomics (PONG). Because HIBAG was designed to impute HLA alleles also from whole-genome SNP data, PONG allows combinatorial diversity of KIR3DL1/S1 with HLA-A and -B to be analyzed using complementary techniques on a single data source. The use of PONG thus negates the need for targeted sequencing data in very large-scale association studies where such methods might not be tractable.
  • Item
    Thumbnail Image
    Large-Scale Imputation of KIR Copy Number and HLA Alleles in North American and European Psoriasis Case-Control Cohorts Reveals Association of Inhibitory KIR2DL2 With Psoriasis
    Ahn, R ; Vukcevic, D ; Motyer, A ; Nititham, J ; Squire, DM ; Hollenbach, JA ; Norman, PJ ; Ellinghaus, E ; Nair, RP ; Tsoi, LC ; Oksenberg, J ; Foerster, J ; Lieb, W ; Weidinger, S ; Franke, A ; Elder, JT ; Jorgenson, E ; Leslie, S ; Liao, W (FRONTIERS MEDIA SA, 2021-06-11)
    Killer cell immunoglobulin-like receptors (KIR) regulate immune responses in NK and CD8+ T cells via interaction with HLA ligands. KIR genes, including KIR2DS1, KIR3DL1, and KIR3DS1 have previously been implicated in psoriasis susceptibility. However, these previous studies were constrained to small sample sizes, in part due to the time and expense required for direct genotyping of KIR genes. Here, we implemented KIR*IMP to impute KIR copy number from single-nucleotide polymorphisms (SNPs) on chromosome 19 in the discovery cohort (n=11,912) from the PAGE consortium, University of California San Francisco, and the University of Dundee, and in a replication cohort (n=66,357) from Kaiser Permanente Northern California. Stratified multivariate logistic regression that accounted for patient ancestry and high-risk HLA alleles revealed that KIR2DL2 copy number was significantly associated with psoriasis in the discovery cohort (p ≤ 0.05). The KIR2DL2 copy number association was replicated in the Kaiser Permanente replication cohort. This is the first reported association of KIR2DL2 copy number with psoriasis and highlights the importance of KIR genetics in the pathogenesis of psoriasis.
  • Item
    Thumbnail Image
    Transethnic analysis of the human leukocyte antigen region for ulcerative colitis reveals not only shared but also ethnicity-specific disease associations
    Degenhardt, F ; Mayr, G ; Wendorff, M ; Boucher, G ; Ellinghaus, E ; Ellinghaus, D ; ElAbd, H ; Rosati, E ; Huebenthal, M ; Juzenas, S ; Abedian, S ; Vahedi, H ; Thelma, BK ; Yang, S-K ; Ye, BD ; Cheon, JH ; Datta, LW ; Daryani, NE ; Ellul, P ; Esaki, M ; Fuyuno, Y ; McGovern, DPB ; Haritunians, T ; Hong, M ; Juyal, G ; Jung, ES ; Kubo, M ; Kugathasan, S ; Lenz, TL ; Leslie, S ; Malekzadeh, R ; Midha, V ; Motyer, A ; Ng, SC ; Okou, DT ; Raychaudhuri, S ; Schembri, J ; Schreiber, S ; Song, K ; Sood, A ; Takahashi, A ; Torres, EA ; Umeno, J ; Alizadeh, BZ ; Weersma, RK ; Wong, SH ; Yamazaki, K ; Karlsen, TH ; Rioux, JD ; Brant, SR ; Franke, A (OXFORD UNIV PRESS, 2021-03-01)
    Inflammatory bowel disease (IBD) is a chronic inflammatory disease of the gut. Genetic association studies have identified the highly variable human leukocyte antigen (HLA) region as the strongest susceptibility locus for IBD and specifically DRB1*01:03 as a determining factor for ulcerative colitis (UC). However, for most of the association signal such as delineation could not be made because of tight structures of linkage disequilibrium within the HLA. The aim of this study was therefore to further characterize the HLA signal using a transethnic approach. We performed a comprehensive fine mapping of single HLA alleles in UC in a cohort of 9272 individuals with African American, East Asian, Puerto Rican, Indian and Iranian descent and 40 691 previously analyzed Caucasians, additionally analyzing whole HLA haplotypes. We computationally characterized the binding of associated HLA alleles to human self-peptides and analyzed the physicochemical properties of the HLA proteins and predicted self-peptidomes. Highlighting alleles of the HLA-DRB1*15 group and their correlated HLA-DQ-DR haplotypes, we not only identified consistent associations (regarding effects directions/magnitudes) across different ethnicities but also identified population-specific signals (regarding differences in allele frequencies). We observed that DRB1*01:03 is mostly present in individuals of Western European descent and hardly present in non-Caucasian individuals. We found peptides predicted to bind to risk HLA alleles to be rich in positively charged amino acids. We conclude that the HLA plays an important role for UC susceptibility across different ethnicities. This research further implicates specific features of peptides that are predicted to bind risk and protective HLA proteins.
  • Item
    Thumbnail Image
    The UK Biobank resource with deep phenotyping and genomic data
    Bycroft, C ; Freeman, C ; Petkova, D ; Band, G ; Elliott, LT ; Sharp, K ; Motyer, A ; Vukcevic, D ; Delaneau, O ; O'Connell, J ; Cortes, A ; Welsh, S ; Young, A ; Effingham, M ; McVean, G ; Leslie, S ; Allen, N ; Donnelly, P ; Marchini, J (NATURE PORTFOLIO, 2018-10-11)
    The UK Biobank project is a prospective cohort study with deep genetic and phenotypic data collected on approximately 500,000 individuals from across the United Kingdom, aged between 40 and 69 at recruitment. The open resource is unique in its size and scope. A rich variety of phenotypic and health-related information is available on each participant, including biological measurements, lifestyle indicators, biomarkers in blood and urine, and imaging of the body and brain. Follow-up information is provided by linking health and medical records. Genome-wide genotype data have been collected on all participants, providing many opportunities for the discovery of new genetic associations and the genetic bases of complex traits. Here we describe the centralized analysis of the genetic data, including genotype quality, properties of population structure and relatedness of the genetic data, and efficient phasing and genotype imputation that increases the number of testable variants to around 96 million. Classical allelic variation at 11 human leukocyte antigen genes was imputed, resulting in the recovery of signals with known associations between human leukocyte antigen alleles and many diseases.
  • Item
    Thumbnail Image
    Clinical and Genetic Risk Factors Associated with Psoriatic Arthritis among Patients with Psoriasis
    Yan, D ; Ahn, R ; Leslie, S ; Liao, W (ADIS INT LTD, 2018-12)
    INTRODUCTION: Psoriatic arthritis (PsA) is a chronic, inflammatory arthritis that affects an estimated 30% of patients with psoriasis. PsA is underdiagnosed in primary care and dermatology clinics due to a variety of reasons, including failure of healthcare providers to ask about symptoms, overlap of symptoms and signs with other rheumatologic conditions, and lack of a specific diagnostic test. A delay in PsA diagnosis and treatment, even as short as 6 months, can lead to decreased quality of life, increased joint damage, and worse long-term physical function. In this study, we sought to identify the clinical and genetic factors that help discriminate patients with PsA from those with cutaneous psoriasis only. METHODS: We analyzed a cohort of 974 psoriasis patients at an academic medical center, of whom 175 had confirmed PsA, and performed univariate, multivariate, and predictive modeling to determine factors associated with PsA. RESULTS: The univariate analysis revealed significant positive associations of PsA with age, nail involvement, scalp involvement, skin fold involvement, elbow/knee involvement, psoriasis severity, plaque subtype, erythrodermic subtype, hypertension, type 2 diabetes, and coronary artery disease, and a significant negative association of PsA with the human leukocyte antigen (HLA)-C*06:02 allele. In the multivariate analysis, nail involvement, type 2 diabetes, and pustular psoriasis remained significantly associated with PsA, while HLA-C*06:02 positivity remained protective. There was a trend towards an association of PsA with older age, younger age of psoriasis onset, and skin fold involvement, while there was protective trend for smoking. A predictive model including both clinical and genetic factors showed reasonable discriminative ability between psoriasis and PsA, with an area under the curve of 0.87 for a receiver operating characteristic curve. CONCLUSION: This study identified a number of clinical and genetic features that could help stratify patients who are at higher risk for having PsA and for whom rheumatology referral may be beneficial.
  • Item
    No Preview Available
    Bayesian analysis of genetic association across tree-structured routine healthcare data in the UK Biobank
    Cortes, A ; Dendrou, CA ; Motyer, A ; Jostins, L ; Vukcevic, D ; Dilthey, A ; Donnelly, P ; Leslie, S ; Fugger, L ; McVean, G (NATURE PUBLISHING GROUP, 2017-09)
    Genetic discovery from the multitude of phenotypes extractable from routine healthcare data can transform understanding of the human phenome and accelerate progress toward precision medicine. However, a critical question when analyzing high-dimensional and heterogeneous data is how best to interrogate increasingly specific subphenotypes while retaining statistical power to detect genetic associations. Here we develop and employ a new Bayesian analysis framework that exploits the hierarchical structure of diagnosis classifications to analyze genetic variants against UK Biobank disease phenotypes derived from self-reporting and hospital episode statistics. Our method displays a more than 20% increase in power to detect genetic effects over other approaches and identifies new associations between classical human leukocyte antigen (HLA) alleles and common immune-mediated diseases (IMDs). By applying the approach to genetic risk scores (GRSs), we show the extent of genetic sharing among IMDs and expose differences in disease perception or diagnosis with potential clinical implications.
  • Item
    No Preview Available
    Genetics of gene expression in primary immune cells identifies cell type-specific master regulators and roles of HLA alleles
    Fairfax, BP ; Makino, S ; Radhakrishnan, J ; Plant, K ; Leslie, S ; Dilthey, A ; Ellis, P ; Langford, C ; Vannberg, FO ; Knight, JC (NATURE PORTFOLIO, 2012-05)
    Trans-acting genetic variants have a substantial, albeit poorly characterized, role in the heritable determination of gene expression. Using paired purified primary monocytes and B cells, we identify new predominantly cell type-specific cis and trans expression quantitative trait loci (eQTLs), including multi-locus trans associations to LYZ and KLF4 in monocytes and B cells, respectively. Additionally, we observe a B cell-specific trans association of rs11171739 at 12q13.2, a known autoimmune disease locus, with IP6K2 (P = 5.8 × 10(-15)), PRIC285 (P = 3.0 × 10(-10)) and an upstream region of CDKN1A (P = 2 × 10(-52)), suggesting roles for cell cycle regulation and peroxisome proliferator-activated receptor γ (PPARγ) signaling in autoimmune pathogenesis. We also find that specific human leukocyte antigen (HLA) alleles form trans associations with the expression of AOAH and ARHGAP24 in monocytes but not in B cells. In summary, we show that mapping gene expression in defined primary cell populations identifies new cell type-specific trans-regulated networks and provides insights into the genetic basis of disease susceptibility.
  • Item
    No Preview Available
    High-density mapping of the MHC identifies a shared role for HLA-DRB101:03 in inflammatory bowel diseases and heterozygous advantage in ulcerative colitis
    Goyette, P ; Boucher, G ; Mallon, D ; Ellinghaust, E ; Jostins, L ; Huang, H ; Ripke, S ; Gusareva, ES ; Annese, V ; Hauser, SL ; Oksenberg, JR ; Thomsent, I ; Leslie, S ; Daly, MJ ; Van Steen, K ; Duerr, RH ; Barrett, JC ; McGovern, DPB ; Schumm, LP ; Traherne, JA ; Carrington, MN ; Kosmoliaptsis, V ; Karsen, TH ; Franke, A ; Rioux, JD (NATURE PUBLISHING GROUP, 2015-02)
    Genome-wide association studies of the related chronic inflammatory bowel diseases (IBD) known as Crohn's disease and ulcerative colitis have shown strong evidence of association to the major histocompatibility complex (MHC). This region encodes a large number of immunological candidates, including the antigen-presenting classical human leukocyte antigen (HLA) molecules. Studies in IBD have indicated that multiple independent associations exist at HLA and non-HLA genes, but they have lacked the statistical power to define the architecture of association and causal alleles. To address this, we performed high-density SNP typing of the MHC in >32,000 individuals with IBD, implicating multiple HLA alleles, with a primary role for HLA-DRB1*01:03 in both Crohn's disease and ulcerative colitis. Noteworthy differences were observed between these diseases, including a predominant role for class II HLA variants and heterozygous advantage observed in ulcerative colitis, suggesting an important role of the adaptive immune response in the colonic environment in the pathogenesis of IBD.
  • Item
    Thumbnail Image
    Multi-Population Classical HLA Type Imputation
    Dilthey, A ; Leslie, S ; Moutsianas, L ; Shen, J ; Cox, C ; Nelson, MR ; McVean, G ; Browning, S (PUBLIC LIBRARY SCIENCE, 2013-02)
    Statistical imputation of classical HLA alleles in case-control studies has become established as a valuable tool for identifying and fine-mapping signals of disease association in the MHC. Imputation into diverse populations has, however, remained challenging, mainly because of the additional haplotypic heterogeneity introduced by combining reference panels of different sources. We present an HLA type imputation model, HLA*IMP:02, designed to operate on a multi-population reference panel. HLA*IMP:02 is based on a graphical representation of haplotype structure. We present a probabilistic algorithm to build such models for the HLA region, accommodating genotyping error, haplotypic heterogeneity and the need for maximum accuracy at the HLA loci, generalizing the work of Browning and Browning (2007) and Ron et al. (1998). HLA*IMP:02 achieves an average 4-digit imputation accuracy on diverse European panels of 97% (call rate 97%). On non-European samples, 2-digit performance is over 90% for most loci and ethnicities where data available. HLA*IMP:02 supports imputation of HLA-DPB1 and HLA-DRB3-5, is highly tolerant of missing data in the imputation panel and works on standard genotype data from popular genotyping chips. It is publicly available in source code and as a user-friendly web service framework.