RESUMO
Identification of candidate causal variants in regions associated with risk of common diseases is complicated by linkage disequilibrium (LD) and multiple association signals. Nonetheless, accurate maps of these variants are needed, both to fully exploit detailed cell specific chromatin annotation data to highlight disease causal mechanisms and cells, and for design of the functional studies that will ultimately be required to confirm causal mechanisms. We adapted a Bayesian evolutionary stochastic search algorithm to the fine mapping problem, and demonstrated its improved performance over conventional stepwise and regularised regression through simulation studies. We then applied it to fine map the established multiple sclerosis (MS) and type 1 diabetes (T1D) associations in the IL-2RA (CD25) gene region. For T1D, both stepwise and stochastic search approaches identified four T1D association signals, with the major effect tagged by the single nucleotide polymorphism, rs12722496. In contrast, for MS, the stochastic search found two distinct competing models: a single candidate causal variant, tagged by rs2104286 and reported previously using stepwise analysis; and a more complex model with two association signals, one of which was tagged by the major T1D associated rs12722496 and the other by rs56382813. There is low to moderate LD between rs2104286 and both rs12722496 and rs56382813 (r2 ≃ 0:3) and our two SNP model could not be recovered through a forward stepwise search after conditioning on rs2104286. Both signals in the two variant model for MS affect CD25 expression on distinct subpopulations of CD4+ T cells, which are key cells in the autoimmune process. The results support a shared causal variant for T1D and MS. Our study illustrates the benefit of using a purposely designed model search strategy for fine mapping and the advantage of combining disease and protein expression data.
Assuntos
Teorema de Bayes , Mapeamento Cromossômico/métodos , Diabetes Mellitus Tipo 1/genética , Predisposição Genética para Doença , Esclerose Múltipla/genética , Algoritmos , Mapeamento Cromossômico/estatística & dados numéricos , Haplótipos , Humanos , Subunidade alfa de Receptor de Interleucina-2/genética , Desequilíbrio de Ligação , Polimorfismo de Nucleotídeo Único , Processos EstocásticosRESUMO
The maintenance of peripheral naive T lymphocytes in humans is dependent on their homeostatic division, not continuing emigration from the thymus, which undergoes involution with age. However, postthymic maintenance of naive T cells is still poorly understood. Previously we reported that recent thymic emigrants (RTEs) are contained in CD31+CD25- naive T cells as defined by their levels of signal joint T cell receptor rearrangement excision circles (sjTRECs). Here, by differential gene expression analysis followed by protein expression and functional studies, we define that the naive T cells having divided the least since thymic emigration express complement receptors (CR1 and CR2) known to bind complement C3b- and C3d-decorated microbial products and, following activation, produce IL-8 (CXCL8), a major chemoattractant for neutrophils in bacterial defense. We also observed an IL-8-producing memory T cell subpopulation coexpressing CR1 and CR2 and with a gene expression signature resembling that of RTEs. The functions of CR1 and CR2 on T cells remain to be determined, but we note that CR2 is the receptor for Epstein-Barr virus, which is a cause of T cell lymphomas and a candidate environmental factor in autoimmune disease.