72
UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ-REITORIA DE PÓS-GRADUAÇÃO PROGRAMA DE PÓS-GRADUAÇÃO EM BIOQUÍMICA TATIANA KUMMER DA ROCHA PINHEIRO FATORES GENÉTICOS DE SUSCEPTIBILIDADE À HANSENÍASE NO RIO GRANDE DO NORTE NATAL 2017

UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

  • Upload
    others

  • View
    2

  • Download
    0

Embed Size (px)

Citation preview

Page 1: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE

PRÓ-REITORIA DE PÓS-GRADUAÇÃO PROGRAMA DE PÓS-GRADUAÇÃO EM BIOQUÍMICA

TATIANA KUMMER DA ROCHA PINHEIRO

FATORES GENÉTICOS DE SUSCEPTIBILIDADE À HANSENÍASE NO RIO GRANDE DO NORTE

NATAL 2017

Page 2: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

TATIANA KUMMER DA ROCHA PINHEIRO

FATORES GENÉTICOS DE SUSCEPTIBILIDADE À HANSENÍASE NO RIO GRANDE DO NORTE

Dissertação apresentada ao Programa de Pós-Graduação em Bioquímica da Universidade Federal do Rio Grande do Norte como requisito para obtenção do título de Mestre em Bioquímica. Orientador(a): Selma Maria Bezerra Jerônimo Co-orientadores: Carlos Gomes Maia Leonardo Capistrano Ferreira

NATAL 2017

Page 3: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

Universidade Federal do Rio Grande do Norte - UFRN Sistema de Bibliotecas - SISBI

Catalogação de Publicação na Fonte. UFRN - Biblioteca Setorial Prof. Leopoldo Nelson - -Centro de Biociências - CB

Pinheiro, Tatiana Kummer da Rocha.

Fatores genéticos de susceptibilidade à hanseníase no Rio Grande do Norte / Tatiana Kummer da Rocha Pinheiro. - Natal,

2017.

71 f.: il.

Dissertação (Mestrado) - Universidade Federal do Rio Grande do

Norte. Centro de Biociências. Programa de Pós-Graduação em

Bioquímica.

Orientadora: Profa. Dra. Selma Maria Bezerra Jerônimo. Coorientador: Prof. Dr. Carlos Gomes Maia.

Coorientador: Prof. Dr. Leonardo Capistrano Ferreira.

1. Reações Hansênicas - Dissertação. 2. M. leprae -

Dissertação. 3. Eritema nodoso hansênico - Dissertação. 4. Reação

reversa - Dissertação. 5. Immunochip - Dissertação. 6. SNP -

Dissertação. I. JERÔNIMO, Selma Maria Bezerra et al. II. Título.

Elaborado por KATIA REJANE DA SILVA - CRB-15/351

Page 4: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

TATIANA KUMMER DA ROCHA PINHEIRO

FATORES GENÉTICOS DE SUSCEPTIBILIDADE À HANSENÍASE NO RIO GRANDE DO NORTE

Dissertação apresentada ao Programa de Pós-Graduação em Bioquímica da Universidade Federal do Rio Grande do Norte como requisito para obtenção do título de Mestre em Bioquímica.

Aprovada em: 29/09/2017.

BANCA EXAMINADORA

Page 5: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

Aos meus pais,

Shirley Kummer da Rocha, e

Máspoli Ramassotti de Azevedo Pinheiro

Page 6: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

AGRADECIMENTOS

Aos meus familiares que sempre torcem pelo meu sucesso por sempre estarem comigo quando precisei. Principalmente aos meus pais, por serem exemplos de superação, de honestidade e dedicação. Aos amigos que participaram direta ou indiretamente para a conclusão deste trabalho, em especial a Raniery, por ter estado ao meu lado ao longo desta trajetória. À professora Selma Jerônimo, pela oportunidade de fazer parte do laboratório e da linha de pesquisa de hanseníase. Apesar de ter sido um grande desafio, foi um período de muito aprendizado. Sou grata por todas as oportunidades, por cada ensinamento e por ter sido tão compreensiva frente às minhas dificuldades no final do mestrado. Aos professores membros da banca de defesa, por terem aceite o convite. Aos professores Leonardo Capistrano e Carlos Maia, por todas as discussões e e ajuda com as análises estatísticas. Aos professores que fizeram parte da banca de qualificação, Dra. Adriana Uchoa, Dra. Adriana Rezende e Dra. Daniella Martins. Como também aos demais professores do programa de Pós-Graduação em Bioquímica da UFRN pelo aprendizado nas diversas áreas. A todos os colegas do Laboratório de Imunogenética: Francianne, Carol, Flávio, Glória, Freire, Paulo, Joanna, João Neto, Ilka e Washington, e aos alunos de iniciação científica da “Equipe Hansen”: Alesson, Larissa e Anny, pela ajuda na realização do trabalho nos mais diversos momentos.

À Margarita, sempre prestativa quando precisei. A equipe de técnicos de enfermagem de dermatologia do Hospital Giselda Trigueiro e do Hospital Universitário Onofre Lopes em Natal, do Centro Clínico Prof Vingt Um Rosado em Mossoró, e do Centro Clínico Integrado em Nova Cruz, que trabalhou diretamente conosco na coleta de amostras e auxiliando na coleta de informações para projeto de pesquisa da hanseníase. Em especial a Mércia e Eva, que me auxiliaram na coleta de informações em Mossoró. À Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), e ao Instituto Nacional de Ciência e Tecnologia - Doenças Tropicais (INCT-DT) pelo apoio financeiro.

Page 7: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

RESUMO

Hanseníase é uma doença infecciosa crônica, de evolução lenta, com um amplo espectro de manifestações clínicas, causada pela Mycobacterium leprae. O estado do Rio Grade do Norte tem baixo coeficiente de detecção de novos casos, mas algumas localidades apresentam áreas focais com altas taxas de detecção, como a cidade de Mossoró, que apresentou uma taxa de detecção de 39,73 por 100.000 habitantes em 2012. Após a exposição ao bacilo, cerca de 10% das pessoas evoluem para doença, com espectro de apresentações, variando do pólo tuberculóide (Tuberculóide-tuberculóide; borderline-tuberculóide), do pólo intermediário (Borderline-borderline) ao pólo lepromatoso (Borderline lepromatoso; Lepromatoso-lepromatoso). A Organização Mundial de Saúde classifica a doença de acordo com o número de lesões, para fins terapêuticos considerando a presença de até cinco lesões como paucibacilares (PB), e os casos com mais de cinco lesões como multibacilares (MB). Aproximadamente um terço das pessoas com hanseníase desenvolvem reações imunopatológicas, classificadas como tipo I, tipo II e neurite. A evolução pós-infecção com M. leprae é influenciada por fatores ambientais e pelo repertório genético do hospedeiro. O objetivo deste estudo foi analisar polimorfismos de nucleotídeo único (SNPs) utilizando um array contendo 196.524 SNPs (Immunochip) em pessoas expostas ao M. leprae. O desenho de estudo foi do tipo caso-controle com inclusão de casos de hanseníase e comunicantes. Os participantes foram caracterizados de acordo com a exposição ao M. leprae, considerando presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação operacional da hanseníase e com o tipo de reação, sendo mais elevada em casos MB e aqueles com reação do tipo II. Todos os participantes também foram genotipados com uso de Immunochip. Para análise dos dados de genotipagem foram considerados: hanseníase vs comunicantes, reação vs não-reação, e taxa de anticorpos anti-LID-NDO. Um total 55 SNPs evidenciaram associação à hanseníase e ao nível de anticorpos. No grupo hanseníase vs comunicantes, 10 SNPs demonstraram associação, sendo 3 relacionados a resposta imune (FASLG, TNFS18, EBF1 e ICOSLG), além de um SNP próximo a VDR, cuja proteína está relacionada com imunomodulação associado a vitamina D3 em monócitos, macrófagos e linfócitos ativados. No grupo reação vs não-reação foi observada a associação de 30 SNPs, 20 próximos a genes conhecidos como de susceptibilidade a psoríase, além de um SNP próximo ao gene UBD. Os níveis de anticorpos LID-NDO estavam associados a 15 SNPs, um deles relacionado ao gene THEMIS, que codifica uma proteína com papel regulador na seleção de células T. Dentre os 55 SNPs associados a um dos fenótipos, 4 se encontram em regiões codantes. Juntos os dados sugerem que a susceptibilidade à hanseníase e o desenvolvimento de reações hansênicas estão ligados à resposta imune celular e humoral do hospedeiro e potencialmente poderiam ser modulados. Palavras-chave: Reações hansênicas. M. leprae. Eritema nodoso hansênico. Reação reversa. Immunochip. SNP.

Page 8: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

ABSTRACT

Leprosy is a chronic, slowly progressing infectious disease with a broad spectrum of clinical manifestations and it is caused by Mycobacterium leprae infection The state of Rio Grande do Norte has a low coefficient detection for new cases, but some localities have focal areas with high detection rates, such as the city of Mossoró, which had a detection rate of 39.7 per 100,000 habitants in 2012. After exposure to the bacillus, about 10% of people develop disease, with a wide spectrum of presentations, ranging from the tuberculoid pole (tuberculoid-tuberculoid, borderline-tuberculoid), the borderline (borderline-borderline), to the lepromatous pole (lepromatous-lepromatous). The World Health Organization also classifies the disease according to the number of lesions, for therapeutic purposes, classifying up to five lesions as paucibacillary (PB), and and more than five lesions, as multibacillary (MB). Approximately one third of people with leprosy develop immunopathological reactions, classified as type I, type II and neuritis. The outcome of M. leprae infection is influenced by environmental factors and by the genetic repertoire of the host. The present study was used to analyze nucleotide polymorphisms (SNPs) using a set containing 196,524 SNPs (Immunochip) in people exposed to M. leprae. The study had a case-control design with recruitment of cases of leprosy and contacts. The participants were phenotypically characterized according to exposure to M. leprae, considering the presence of anti-LID-NDO antibodies and clinical presentation. The amount of antibodies varied according to the operational classification of the leprosy and the type of reaction, being higher in MB cases and those with type II reaction. All participants were also genotyped using Immunochip. For analysis of the genotyping data were considered: leprosy vs contacts, reaction vs non-reaction, and anti-LID-NDO antibody rate. A total of 55 SNPs showed association with leprosy and antibody levels. In the leprosy vs contacts group, 10 SNPs showed association, 3 related to the immune response (FASLG, TNFS18, EBF1 and ICOSLG), as well as one SNP close to VDR, whose protein is related to immunomodulation associated with vitamin D3 in monocytes, macrophages and activated lymphocytes. In the reaction vs. non-reaction group, the association of 30 SNPs, of which 20 were to genes known as susceptibility to psoriasis, and a SNP close to the UBD gene. Levels of LID-NDO antibodies were associated with 15 SNPs, one of them related to the THEMIS gene, which encodes a regulatory protein in T-cell selection. Of the 55 SNPs associated with one of the phenotypes, 4 are found in the coding regions. Together the data suggest that susceptibility to leprosy and the development of leprosy reactions are linked to the cellular and humoral immune response of the host and could potentially be modulated.

Key-words: Leprosy reactions. M. leprae. Erythema Nodosum Leprosum. Reversal reaction. Immunochip.SNP.

Page 9: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

LISTA DE FIGURAS

FIGURA 1 Distribuição global de novos casos de hanseníase em 2015.................................. 16

FIGURA 2 Taxa de detecção de hanseníase na cidade de Mossoró de 2003 a 2012 ............... 17

FIGURA 3 Variação clínico- imunopatológica da hanseníase ................................................. 19

FIGURA 4 População do Rio Grande do Norte incluída no estudo separadas em comunicantes e casos. .............................................................................................................................. 28

FIGURA 5 Esquema da metodologia de desenvolvimento do estudo, dividida em caracterização fenotípica e genotípica da população recrutada. ....................................... 29

FIGURA 6 Protocolo de genotipagem por Immunochip .......................................................... 32

FIGURA 7 Esquema do controle de qualidade por-indivíduo. ................................................ 34

FIGURA 8 Informação de sexo discordante.. .......................................................................... 35

FIGURA 9 Ausência de genótipos e taxa de heterozigosidade.. .............................................. 36

FIGURA 10 Indivíduos duplicados ou relacionados. ............................................................... 37

FIGURA 11 Avaliação da estratificação populacional.. .......................................................... 38

FIGURA 12 Esquema do controle de qualidade por-marcador. .............................................. 39

FIGURA 13 Análise comparativa entre os níveis de anticorpos específicos LID-NDO entre comunicantes e casos. ....................................................................................................... 42

FIGURA 14 Manhattan plot demonstrando p valores no genoma total em associação com Hanseníase. ....................................................................................................................... 44

Page 10: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

FIGURA 15 Manhattan plot demonstrando p valores no genoma total em associação com

reações hansênicas ............................................................................................................ 46

FIGURA 16 Parcelas de associação regional em 6p21.33 ....................................................... 47

FIGURA 17 Manhattan plot demonstrando p valores no genoma total em associação com taxa de anticorpos LID-NDO ................................................................................................... 49

Page 11: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

LISTA DE TABELAS

TABELA 1 Fatores de risco genéticos para a hanseníase identificados em estudos por associação de genes candidatos e análises de ligação de todo o genoma ......................... 24

TABELA 2 Características da população de estudo ................................................................ 41

TABELA 3 Associações com Hanseníase para 10 SNPs de susceptibilidade, de acordo com análise ............................................................................................................................... 45

TABELA 4 Associações com reações hansênicas para 30 SNPs de susceptibilidade, de acordo com análise ....................................................................................................................... 48

TABELA 5 Associações com taxas de anticorpos LID-NDO para 15 SNPs de susceptibilidade, de acordo com análise ........................................................................... 50

TABELA 6 55 SNPs associados à hanseníase, dentro de 36 regiões gênicas e intergênicas, de acordo com análise ........................................................................................................... 51

Page 12: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

LISTA DE ABREVIATURAS / SIGLAS

BAAR Bacilo álcool-ácido resistente

BB Borderline-borderline

BL Borderline-lepromatosa

BT Borderline-tuberculóide

ENH Eritema nodoso hansênico

Hemonorte Hemocentro do Rio Grande do Norte

HHC Comunicantes (Household Contacts)

HLA Antígeno leucocitário humano

HWE Equilíbrio de Hardy-Weinberg (Hardy–Weinberg equilibrium)

IB Índice baciloscópico

IFN Interferon

IL Interleucina

LD Desequilíbrio de ligação (Linkage disequilibrium)

LID-NDO Antígeno LID-1 ligado ao dissacarídio natural do PGL-I

LL Lepromatosa-lepromatosa

MB Multibacilares

MHC Complexo principal de histocompatibilidade

NF-κβ Fator nuclear κβ

OD Densidade óptica (Optical density)

OMS Organização Mundial de Saúde

PB Paucibacilares

PQT Tratamento poliquimioterápico

RR Reação reversa

SNPs Polimorfismos de nucleotídeo único (Single-nucleotide polymorphism)

Th Linfócito T auxiliar

Page 13: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

TT Tuberculóide

TCLE Termo de consentimento livre e esclarecido

TGF-β Fator de crescimento tumoral-β

TNF-α Fator de necrose tumoral alfa

TLR Receptor do tipo Toll

VDR Receptor de vitamina D (Vitamin D receptor)

Page 14: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

SUMÁRIO

1.1 Introdução ................................................................................................................... 14

1.2 Aspectos gerais ............................................................................................................ 14

1.3 Revisão da literatura .................................................................................................. 16

1.3.1 Epidemiologia da hanseníase .............................................................................. 16

1.3.2 Patogênese e classificação da hanseníase ........................................................... 17

1.3.3 Imunologia da hanseníase ................................................................................... 20

1.3.4 Genética da hanseníase ........................................................................................ 22

2 OBJETIVOS ............................................................................................................... 27

2.1 Objetivo geral .............................................................................................................. 27

2.2 Objetivos específicos ................................................................................................... 27

3 MATERIAIS E MÉTODOS ...................................................................................... 28

3.1 População e desenho de estudo .................................................................................. 28

3.2 Considerações éticas ................................................................................................... 30

3.3 Determinação de anticorpo anti-LID-NDO .............................................................. 30

3.4 Extração de DNA e genotipagem ............................................................................... 31

3.5 Controle de qualidade ................................................................................................ 33

3.5.1 Controle de qualidade por indivíduo ................................................................... 33

3.5.2 Controle de qualidade por marcador ................................................................... 38

3.6 Análise estatística ........................................................................................................ 40

4 RESULTADOS ........................................................................................................... 41

4.1 Características populacionais .................................................................................... 41

4.2 Determinação dos níveis de anticorpo anti-Mycobacterium leprae ........................ 42

Page 15: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

4.3 Análise genética ........................................................................................................... 43

4.3.1 Hanseníase vs. Comunicantes ............................................................................. 43

4.3.2 Reação vs. Não reação ......................................................................................... 46

4.3.3 Nível de anticorpos anti-LID-NDO ..................................................................... 49

5 DISCUSSÃO ............................................................................................................... 53

6 CONCLUSÕES ........................................................................................................... 60

REFERÊNCIAS ..................................................................................................................... 61

APÊNDICE 1: TERMO DE CONSENTIMENTO LIVRE E ESCLARECIDO (TCLE) 68

Page 16: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

14

1.1 Introdução

1.2 Aspectos gerais

A hanseníase é uma doença infecciosa crônica, complexa, de evolução lenta, causada

pela bactéria Mycobacterium leprae, apresentando um amplo espectro de manifestações

clínicas, principalmente sintomas dermatoneurológicos. Uma das principais características da

hanseníase é o comprometimento de nervos periféricos, o que lhe dá potencial para promover

deformidades e incapacidades temporárias ou permanentes, as quais levam à diminuição da

capacidade laborativa e da qualidade de vida do paciente (CARDOSO et al., 2011;

SCOLLARD et al., 2006). O M. leprae é uma bactéria intracelular com um longo tempo de

proliferação (COLE et al., 2001). A inexistência de meio de cultura e modelo experimental

para este microrganismo dificultam a realização de trabalhos experimentais para testar

hipóteses sobre a infecção e a doença.

A hanseníase ainda é um importante problema de saúde pública. Apesar da melhoria

do seu controle, após a implementação do tratamento poliquimioterápico, em 1981, novos

casos continuam a ocorrer em países endêmicos, sendo o Brasil o segundo país em número de

novos casos (WHO, 2016). O Rio Grande do Norte é um estado com número relativamente

baixo de casos. No entanto, existem regiões apresentando altas taxas de detecção, como

exemplo da cidade de Mossoró com uma taxa de detecção de 39,73 por 100.000 habitantes em

2012 (MINISTÉRIO DA SAÚDE, 2013), que somente recentemente foram reconhecidas

devido a detecção aprimorada de casos ativos, demonstrando um pico de novos casos em

2005 (QUEIROZ et al., 2010).

Após a exposição ao bacilo do M. leprae pelas vias áreas superiores, o indivíduo pode

ou não evoluir para doença. Após o desenvolvimento da hanseníase per se, o paciente poderá

obter cura espontânea ou desenvolver doença com formas variadas de apresentação. Ridley e

Jopling desenvolveram uma classificação de hanseníase em cinco formas clínicas (RIDLEY;

JOPLING, 1962), dependendo do nível de susceptibilidade e o tipo de resposta imune contra a

micobactéria (tuberculóide, borderline-tuberculóide, borderline-borderline, borderline-

lepromatosa e lepromatosa). Cerca de 30% das pessoas com hanseníase desenvolvem reações

imunopatológicas (reação reversa ou eritema nodoso hansênico), que são causas adicionais de

complicações e podem levar a sequelas permanentes (SCOLLARD et al., 2006).

Page 17: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

15

A evolução pós infecção com M. leprae é influenciada pelo tipo de resposta imune

elicitada, podendo esta ter um componente genético. A hipótese de que a composição genética

do hospedeiro influencia na susceptibilidade à hanseníase é a mais aceita na atualidade

(MARQUET; SCHURR, 2001; MONOT et al., 2005; REIBEL; CAMBAU; AUBRY, 2015).

O genoma de M. leprae não possui uma grande variabilidade, pois quatro cepas

geograficamente variadas de M. leprae compartilham 99,995% de identidade na sequência,

diferindo apenas em uma série de pseudogenes e sítios de polimorfismo, principalmente

polimorfismos de nucleotídeo único (SNPs). De forma característica, a M. leprae difere das

demais bactérias patogênicas e não-patogênicas pela grande proporção de pseudogenes, 1293

pseudogenes, com 50% do seu genoma isento de função (SINGH; COLE, 2011).

O risco de desenvolver formas sintomáticas e/ ou formas reacionais na infecção por M.

leprae é em parte geneticamente determinado (CARDOSO et al., 2011; GASCHIGNARD et

al., 2016; MILLER et al., 2004; SANTOS et al., 2002). Dessa forma, existem diversos

trabalhos tentando apontar genes associados à susceptibilidade à hanseníase e suas formas

reacionais (CHAITANYA et al., 2014; FAVA et al., 2017; LIU et al., 2015; NEELA et al.,

2015; SALES-MARQUES et al., 2017; XIANG et al., 2015).

Portanto, com a finalidade de contribuir para o conhecimento da associação de fatores

genéticos à susceptibilidade à hanseníase, o objetivo do presente estudo é analisar

polimorfismos associados com susceptibilidade à hanseníase e/ou formas reacionais em

pessoas oriundas do estado do Rio Grande do Norte.

Page 18: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

16 1.3 Revisão da literatura

1.3.1 Epidemiologia da hanseníase

A hanseníase afeta atualmente mais de 174 mil pessoas em todo o mundo. Três países

apresentam mais de 10 mil novos casos por ano (Figura 1): Índia com 127.326 novos casos

(60% dos novos casos globais), Brasil com 26.395 novos casos (13%) e Indonésia com 17.202

novos casos (8%) (WHO, 2016). Porém, o Brasil está em primeiro lugar em detecção de

novos casos de hanseníase, com 0,0128%, enquanto a Índia e Indonésia apresentam,

respectivamente, 0,0097% e 0,0067%, quando dividimos o número de novos casos pela

população.

Figura 1: Distribuição global de novos casos de hanseníase em 2015. (Adaptado de WHO, 2016).

A tendência de novos casos no Brasil é decrescente, porém ainda se nota uma alta

incidência nas regiões Norte, Nordeste e Centro-Oeste. O estado do Rio Grande do Norte

apresentou no ano de 2012 cerca de 318 novos casos da doença. Segundo parâmetros

mundiais, desde 2004 o estado atingiu a meta de eliminação da OMS (um novo caso a cada 10

mil habitantes), mas ainda existem regiões com altas taxas de detecção no estado, como

exemplo da cidade de Mossoró com uma taxa de detecção de 39,73 por 100.000 habitantes em

2012 (Figura 2) (MINISTÉRIO DA SAÚDE, 2013).

Page 19: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

17

Figura 2: Taxa de detecção de hanseníase na cidade de Mossoró de 2003 a 2012 (Nobre et al, 2015).

Mapeamento por georreferenciamento foi utilizado em Mossoró como ferramenta de

identificação de áreas de alto risco à hanseníase. Os estudos demonstram que em áreas

hiperendêmicas o risco de doença é alto em contatos sociais, além de uma exposição precoce

a M. leprae. Ademais, foi observado que não houve diferença significativa na taxa de

detecção entre contatos domiciliares e vizinhos (MOURA et al., 2013; QUEIROZ et al.,

2010).

1.3.2 Patogênese e classificação da hanseníase

Hanseníase é uma doença infecto-contagiosa crônica causada pelo agente etiológico

Mycobacterium leprae (ALTER et al., 2008). Essa bactéria é um bacilo álcool-ácido resistente

(BAAR) que não possui locomoção própria e não forma esporos, sendo um parasita

intracelular obrigatório, com proliferação muito lenta e não cultivável in vitro (SCOLLARD

et al., 2006; SUZUKI et al., 2012).

O reservatório principal de infecção por M. leprae é o homem, porém infecções

naturais já foram relatadas em animais silvestres como o tatu Dasypus novemcinctus e

diversas espécies de primatas (SUZUKI et al., 2012).

A infecção pelo M. leprae ocorre a partir da exposição do sistema respiratório a

gotículas transportadas pelo ar, maior em espaços onde existem indivíduos que apresentem

alta carga bacilar, pessoas com a forma multibacilar, (SUZUKI et al., 2012). Após o contato

com o M. leprae pelas vias aéreas superiores, em 10% dos indivíduos, o bacilo é difundido

pelos nervos periféricos e pele até as células de Schwann, no sistema nervoso periférico, onde

será fagocitado. O tropismo de M. leprae pelas células de Schwann é devido à ligação

Page 20: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

18 específica entre o M. leprae e a laminina-α2 da membrana basal anexa às células de Schwann.

O complexo M. leprae/laminina-α2 liga-se à α-distroglicana, uma proteína de membrana da

célula de Schwann, contribuindo para a invasão celular (RAMBUKKANA, 2000). A reação

imunológica ocasionada pela presença do patógeno na célula juntamente com a

desmielinização dos nervos periféricos leva ao dano neural atribuído à doença (SCOLLARD

et al., 2006). Portanto, o resultado da infecção e a manifestação clínica dependem da

imunidade celular do hospedeiro, a qual é a primeira linha de defesa contra a infecção por M.

leprae (SUZUKI et al., 2012). Após o período de incubação, o paciente pode evoluir para cura

espontânea ou desenvolver diferentes respostas clínicas (GOULART; PENNA; CUNHA,

2002).

A variação clínico-patológica da hanseníase pode ser resumida na Figura 3. As

manifestações clínicas características da hanseníase são distribuídas em dois pólos e três

formas intermediárias pelos critérios de Ridley-Jopling (RIDLEY; JOPLING, 1962): a forma

tuberculóide (TT), borderline-tuberculóide (BT), borderline-borderline (BB), borderline-

lepromatosa (BL) e a forma lepromatosa (LL).

No primeiro pólo temos a forma tuberculóide, com lesões definidas, hipo-pigmentadas

e anestésicas, vigorosa resposta celular e poucos bacilos, associada a um predomínio de

resposta imunológica do tipo Th1. Os pacientes com a forma TT podem apresentar danos

neurológicos de progressão rápida. No outro pólo, a forma lepromatosa apresenta lesões de

pele difusamente distribuídas e presença de bacilos na pele e nos nervos, com um predomínio

de resposta Th2. Os pacientes LL apresentam complicações a um maior prazo quando

comparados com pacientes TT (BRITTON; LOCKWOOD, 2004; RIDLEY; JOPLING,

1962).

Page 21: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

19

Figura 3: Variação clínico- imunopatológica da hanseníase (Adaptado de MISCH et al, 2010).

Entre estes dois pólos extremos, existem as formas intermediárias, ou borderlines, as

quais representam a maior parte dos casos: borderline-tuberculóide, borderline-borderline e

borderline-lepromatosa, com uma progressiva redução da resposta imune celular

acompanhada por numerosas lesões de pele e nervos, aumento da carga bacilar e dos níveis de

anticorpos (BRITTON; LOCKWOOD, 2004; GOULART; PENNA; CUNHA, 2002;

RIDLEY; JOPLING, 1962).

A Organização Mundial de Saúde (OMS) também classifica a doença de acordo com o

número de lesões, para fins terapêutico, e, por este método, as formas são classificadas como

paucibacilares (PB) os casos com até cinco lesões cutâneas, integrando as formas TT e BT, e

multibacilares os casos com mais de cinco lesões cutâneas (BT, BB, BL e LL) (WHO, 2014).

1.3.2.1 Tratamento da hanseníase

Em 1982 foi implementado o tratamento poliquimioterápico (PQT) formado pelos

medicamentos Rifampicina, Dapsona e Clofazimina pela Organização Mundial de Saúde

(OMS). Desde então, o controle de hanseníase tem melhorado devido a campanhas nos países

nos quais a doença é endêmica.

Page 22: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

20

Os regimes de tratamento são indicados de acordo com a classificação operacional dos

pacientes. Nos indivíduos paucibacilares, a recomendação é de seis doses mensais de

rifampicina 600mg e doses diárias de dapsona 100g. Para os pacientes multibacilares, são

recomendadas doze doses mensais de rifampicina 600mg e clafazimina 300 mg e doses

diárias de dapsona 100mg e clofazimina 50 mg (WHO, 1998).

1.3.3 Imunologia da hanseníase

As manifestações clínicas em resposta à infecção por M. leprae dependem de citocinas

no sítio de presença do parasita, determinando a forma clínica do paciente. Portanto, a

resposta imunológica do paciente contra o patógeno irá definir a severidade e evolução da

doença. Não se sabe integralmente a razão de uma pessoa evoluir para uma forma

tuberculóide ou para uma forma lepromatosa.

Na forma tuberculóide há uma intensa resposta de imunidade celular e uma grande

resposta imune Th1, limitando a progressão da doença a lesões de pele e dano nervoso

definido, sem resposta humoral, formando granulomas bem definidos compostos por

macrófagos epitelióides, cercados principalmente por grande número de linfócitos T e B e

outras células. Lesões tuberculóides exibem maior quantidade de IFN-γ, TNF, IL-2, IL-6 e IL-

12 (figura 3). A abundância de IL-2 e IFN-γ em lesões TT provavelmente contribui para o

estado de imunidade resistente. IFN-γ aumenta a produção de reativos intermediários de

oxigênio e nitrogênio por macrófagos, restrigindo a proliferação de micobactérias. IL-2

contribui para a defesa do hospedeiro por induzir a expansão clonal de células T imune

ativadas e por aumentar a produção de IFN-γ. A produção de TNF-α favorece uma ação

sinérgica autócrina para manter o macrófago ativado e formar o granuloma; enquanto IL-12

estimula células natural killer (NK) a liberar IFN-γ, que modula as células T em direção ao

padrão Th1 (BELONE et al., 2015; BRITTON; LOCKWOOD, 2004; MISCH et al., 2010;

REIBEL; CAMBAU; AUBRY, 2015; SAMPAIO et al., 2012).

Em contrapartida, a forma lepromatosa é caracterizada por uma baixa imunidade

mediada por células e uma resposta Th2 humoral predominante, levando a uma resposta

imune inadequada para o bacilo intracelular. Dessa forma, granulomas no pólo lepromatoso

são difusos, com macrófagos vacuolados e infiltrados linfocíticos escassos, com alta

concentração de IL-4 e IL-10 (figura 3). IL-4 tem um efeito imunorregulatório negativo sobre

a imunidade mediada por células, levando ao aumento da proliferação bacteriana ao bloquear

a proliferação dependente de IL-2 de células T humanas, por inibir receptores de IL-2, e

Page 23: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

21 bloquear a ativação de monócitos mediada pelo IFN-γ, e bloquear a geração de óxido nítrico.

IL-10 inibe a produção de citocinas por células CD4+ na ausência de IFN-γ ou IL-2

(BELONE et al., 2015; BRITTON; LOCKWOOD, 2004; MISCH et al., 2010; REIBEL;

CAMBAU; AUBRY, 2015; SAMPAIO et al., 2012).

Um estudo avaliou a expressão de genes em lesões de pacientes com hanseníase, na

qual no pólo tuberculóide foi observado o aumento da expressão de genes associados ao

processamento e apresentação de antígenos (UBD, PSMB4, CD1b, CD28 e CD79a),

antimicrobiano (cathepsina g e SLP1) e pró-inflamatório/Th1 (SLAM, receptor IL-1, IL-15,

IL-17, IL-12, receptor-beta e linfotoxina-α). Já em pacientes com hanseníase no pólo

lepromatoso, foi observado o aumento da expressão de genes associados à resposta anti-

inflamatória/Th2 (TGFβ1, IL-5 e proteína TGFB 2 latente) e a receptores de inibição (SIRP-

1α, LIR-7, LIR3 e LIR8) (BLEHARSKI et al., 2003).

O glicolípidio fenólico PGL-1, componente de membrana de M. leprae, é detectado na

maioria dos pacientes MB a partir de testes sorológicos, sendo sugerido como ferramenta

adjunta para uma categorização de hanseníase MB e PB, apesar de a maioria de pacientes PB

não possuirem anticorpos anti-PGL-1 detectáveis (BÜHRER-SÉKULA et al., 2007;

SPENCER et al., 2012). Um mimético sintético de PGL-I, o ND-O, comporta-se de forma

semelhante quando conjugado com proteínas transportadoras inertes, como a albumina

(SPENCER; BRENNAN, 2011). O antígeno LID-NDO foi desenvolvido pela conjugação de

dois antígenos, PGL-1 e LID-1, uma proteína fusionada composta pelos produtos dos genes

ML0405 e ML2331 de M. leprae (DUTHIE et al., 2014). Níveis de anticorpos anti-LID-1 se

mostraram elevados e persistentes em pacientes com hanseníase multibacilar que

desenvolveram uma reação, sugerindo, portanto, que esta é potencialmente uma ferramenta

preditiva sorológica (MIZOGUTI et al., 2015). Além disso, um estudo realizado pelo nosso

grupo demonstrou a eficácia de LID-NDO e LID-1 como método diagnóstico para hanseníase

e avaliação de comunicantes (AMORIM et al., 2016).

1.3.3.1 Reações hansênicas

Aproximadamente 70% dos pacientes com hanseníase TT ou BT desenvolvem um ou

mais tipos de reações imunopatológicas, conhecidas como estados reacionais hansênicos, que

ocorrem por alterações no sistema imunológico (BRITTON; LOCKWOOD, 2004). As

reações podem ser Reação Reversa (RR), ou reação tipo I, e Eritema Nodoso Hansênico

(ENH), ou reação tipo II (BRITTON; LOCKWOOD, 2004). Essas reações hansênicas são

Page 24: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

22 comuns durante os primeiros três meses de tratamento, mas podem ocorrer antes, durante ou

após o tratamento, e são tipicamente crônicas, requerendo até anos de tratamento com

corticosteróides ou outros medicamentos imunomoulatórios (DUPNIK et al., 2015).

A reação reversa, ou reação do tipo I, resulta em inchaço com inflamação em lesões da

pele com aumento do dano nervoso. Esse tipo de reação tem relação com um aumento de

resposta imune celular a antígenos M. leprae na pele e nervos, em pacientes borderline. Nas

lesões há um influxo de linfócitos e macrófagos jovens, resultando na desorganização de

granulomas preexistentes, edema e necrose. Células presentes nas lesões expressam as

citocinas pró-inflamatórias interferon gamma, interleucina 12 e óxido nítrico (BELONE et al.,

2015; MANDAL et al., 2015; WALKER; LOCKWOOD, 2008).

O eritema nodoso hansênico, ou reação do tipo II, ocorre em hanseníase lepromatosa

ou borderline-lepromatosa com a presença de nódulos subcutâneos e achados sistêmicos

como febre, artrite, nefrite e paniculite. Neutrófilos são um componente principal e

comumente formam microabcessos. Estudos indicam que esta patogênese pode ser

relacionada à formação e deposição de complexos imunes na pele e órgãos, altos níveis de

TNF-α foram detectados no plasma de pacientes com esta reação (BELONE et al., 2015;

KAHAWITA; LOCKWOOD, 2008; MANDAL et al., 2015).

Outro tipo de reação hansênica é a neurite, uma inflamação aguda dos nervos com dor,

edema local e perda rápida de função do nervo. Este tipo de reação pode ocorrer antes,

durante ou após o tratamento (WHO, 1998). A neurite é muitas vezes devastadora para a

saúde e o bem-estar do paciente, através do desenvolvimento de anestesia, paralisia e

deformidades de dedos das mãos e dos pés (SCOLLARD et al., 2006). Neste estudo, iremos

considerar a neurite como um tipo de reação, devido a diminuição da capacidade laborativa do

paciente.

1.3.4 Genética da hanseníase

Apesar de ser considerada uma doença infectocontagiosa crônica, a evolução da

hanseníase também é influenciada por genes, sendo, então, uma doença complexa na qual

características do hospedeiro e do ambiente são tão importantes quanto as características do

patógeno para determinar a evolução da doença (CARDOSO et al., 2011). Diversos estudos

realizados estabeleceram que entre os contatos próximos de pacientes, aqueles biologicamente

Page 25: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

23 relacionados ao caso estão mais em risco de desenvolver hanseníase (DOS SANTOS et al.,

2013; MOET et al., 2006).

Outros fatores de risco para o desenvolvimento de hanseníase incluem contato

próximo a pacientes, baixo nível educacional, idade (pessoas mais velhas estariam em contato

com o bacilo a mais tempo), falta de vacinação BCG e exposição a tatus (GASCHIGNARD et

al., 2016; MISCH et al., 2010; MOET et al., 2006).

O genoma de M. leprae não possui uma grande variabilidade. Além disso, quando

comparado a outras micobactérias, como a Mycobacterium tuberculosis, a M. leprae

apresenta redução evolutiva no seu genoma, via pseudogenes, como a perda da atividade de

mais de 2000 genes (COLE et al., 2001; MONOT et al., 2009).

Estudo de associação gênica realizado com vinte famílias vietnamitas identificou

ligação significante entre polarização da hanseníase e a região HLA/TNF (antígeno

leucocitário humano/fator de necrose tumoral) no cromossomo 6p21 (MIRA et al., 2003a). Já

o estudo conduzido por Mira et al em populações do Vietnã e do Brasil obtiveram como

resultado a identificação de variantes genéticas na região compartilhada pelos genes PARK2 e

PACRG, gene corregulado por PARK2 como fatores de risco para susceptibilidade à

hanseníase per se (MIRA et al., 2003b, 2004). O gene PARK2 codifica a parquina, uma

proteína E3-ubiquitina-ligase que está envolvida na via de ubiquitinação e degradação de

proteínas no proteossomo, estando também envolvida em apoptose, manutenção das funções

mitocondriais e ativação do fator de transcrição NF-κβ. A associação de que polimorfismos

nos genes PARK2/PACRG possam aumentar a susceptibilidade à hanseníase pode indicar que

o sistema ubiquitina-proteossomo talvez tenha um papel no controle da infecção (CARDOSO

et al., 2011; SCOLLARD et al., 2006).

Fatores de risco genéticos para a hanseníase foram identificados em diversos estudos

utilizando associação de genes candidatos e análises de ligação de todo o genoma (Tabela 1),

encontrando associações em SNPs focados em genes da região do complexo de

histocompatibilidade principal (MHC), receptores do tipo Toll (TLR) e a família das

citocinas, como CCDC122 e LACC1 no Vietnã (GRANT et al., 2012), LRRK2 e RIPK2, da

rede sinalização de NOD2, na Índia (MARCINEK et al., 2013), e NRAMP1, em Mali, que

pode influenciar respostas imunecelulares a patógenos intracelulares (MEISNER et al., 2001).

Page 26: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

24 Tabela 1: Fatores de risco genéticos para a hanseníase identificados em estudos por associação de genes candidatos e análises de ligação de todo o genoma

Gene País Referência

IL-10 Brasil MORAES et al., 2004; SANTOS et al., 2002 IFNγ Brasil REYNARD et al., 2003 MLB2 Brasil DE MESSIAS-REASON et al., 2007 KIRD2DS3 Brasil FRANCESCHI et al., 2008 MICA Brasil DO SACRAMENTO et al., 2012 HLA-G Brasil LUCENA-SILVA et al., 2013 VDR Brasil NEELA et al., 2015 MRC1 Brasil e Vietnã ALTER et al., 2010

CCDC122 China ZHANG et al., 2009 HLA-DR-DQ China ZHANG et al., 2009 LACC1 China ZHANG et al., 2009 LRRK2 China ZHANG et al., 2009 NOD2 China ZHANG et al., 2009 RIPK2 China ZHANG et al., 2009 TNFSF15 China ZHANG et al., 2009 BATF3 China LIU et al., 2015 CCDC88B China LIU et al., 2015 CDH18 China LIU et al., 2015 CIITA-SOCS1 China LIU et al., 2015 DEC1 China LIU et al., 2015 EGR2 China LIU et al., 2015 OPA1 China XIANG et al., 2015 TLR2 Etiópia BOCHUD et al., 2008 TLR4 Etiópia BOCHUD et al., 2009 LTA4H Etiópia TOBIN et al., 2010 HSPA1A India RAJALINGAM; MEHRA; SINGAL, 2000 IL-23R India ALI et al., 2013 LRRK2 India MARCINEK et al., 2013 RIPK2 India MARCINEK et al., 2013 IL-17F India CHAITANYA et al., 2014 LAMA2 Indonésia WIBAWA; SOEBONO; MATSUO, 2002 IL-12RBB2 Japão OHYAMA et al., 2008 LTA Malawi FITNESS et al., 2004 NRAMP1 Mali MEISNER et al., 2001 DEFB1 México PRADO-MONTES DE OCA et al., 2009 MLB2 Nepal SAPKOTA et al., 2010 TNF Tailândia VEJBAESYA et al., 2007 CCDC122 Vietnã GRANT et al., 2012 LACC1 Vietnã GRANT et al., 2012

Page 27: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

25

O primeiro estudo de varredura de associação do genoma (GWAS) para hanseníase foi

publicado em 2009, incluindo 706 casos e 1.225 controles de uma população chinesa

(ZHANG et al., 2009), no qual foram detectados 16 SNPs associados à hanseníase, atribuídos

a sete genes (Tabela 1): HLA-DR-DQ (complexo principal de histocompatibilidade classe II

DR), RIPK2 (receptor que interage com serina-treonina quinase 2), TNFSF15 (membro 15 da

superfamília de fator de necrose tumoral), LRRK2 (receptor quinase de repetição rico em

leucina 2), CCDC122 (domínio de bobina em espiral 122), LACC1 e NOD2 (domínio de

oligomerização de ligação de nucleótidos 2).

Um segundo GWAS realizado na população chinesa, utilizando 8.313 casos e 16.017

controles, identificou seis novos loci de susceptibilidade (Tabela 1): BATF3 (zíper de leucina

básico de tipo ATF), CCDC88B, CIITA-SOCS1 (transativador do complexo de

histocompatibilidade principal classe II - supressor de sinalização de citocinas 1), CDH18

(caderina 18), EGR2 (resposta inicial ao crescimento 2) e DEC1 (suprimido no câncer de

esôfago 1) (LIU et al., 2015).

SNPs associados a reações hansênicas foram encontrados próximos aos genes NINJ1

(CARDOSO et al., 2007), TLR1 (SCHURING et al., 2009), IL6 (SOUSA et al., 2012),

TNFSF8 (FAVA et al., 2015) e NOD2 (SALES-MARQUES et al., 2017). O primeiro GWAS

analisando reações hansênicas indentifou lncRNA como um fator de risco para

desenvolvimento de reação do tipo I (FAVA et al., 2017).

Estudos GWAS examinam um conjunto genômico de variantes genéticas em

diferentes indivíduos para determinar se qualquer variante está associada a uma característica,

baseado no principio de desequilíbrio de ligação (LD), que é a associação não aleatória entre

alelos em loci diferentes criado por forças evolutivas, como mutação, deriva e seleção

(VISSCHER et al., 2012). Porém existem limitações para este tipo de estudo, os chips de

genotipagem no GWAS não são adequados para escolher as variantes comuns restantes ou

identificar variantes raras, além de seu alto custo (em 2011 o valor médio por amostra era de

250 doláres) (CORTES; BROWN, 2011). Dessa forma, o Immunochip Consortium

desenvolveu um chip contendo aproximadamente 200.000 SNPs para genotipar loci de

doenças mediadas pela imunidade identificados por GWAS de variantes comuns utilizando

dados do 1000 Genomes Project e outros dados disponíveis (TRYNKA et al., 2011).

O Immunochip, chip de genotipagem da Illumina Infinium, possui um custo de

aproximadamente 20% do GWAS, e foi desenvolvido cobrindo as principais doenças

Page 28: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

26 autoimunes: artrite reumatóide, espondilite anquilosante, lúpus eritematoso sistêmico,

diabetes tipo 1, doença da tiróide autoimune, doença celíaca, esclerose múltipla, colite

ulcerativa, doença de Crohn e psoríase (CORTES; BROWN, 2011). Dessa forma, o

Immunochip é uma ferramenta valiosa na análise de SNPs associados à hanseníase por

representar aproximadamente 20% do genoma e incluir genes de resposta imune, tais genes

podem estar associados ao risco de desenvolver reações imunopatológicas, pois são um tipo

de resposta autoimune. Além disso, a hanseníase mimetiza algumas doenças autoimune, como

lúpus eritematoso sistêmico, psoaríase e artrite no aspecto clínico, portanto poderia apresentar

SNPs de susceptibilidade em comum (PRASAD et al., 2013).

A hipótese desse estudo é de que a composição genética do hospedeiro influencia na

susceptibilidade à hanseníase e no desenvolvimento de reações. Além disso, também tivemos

como hipótese que as reações hansênicas são, em parte, associadas à produção de anticorpos.

Page 29: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

27 2 OBJETIVOS

2.1 Objetivo geral

Avaliar polimorfismos genéticos (SNPs) associados à susceptibilidade à hanseníase em

indivíduos do estado do Rio Grande do Norte.

2.2 Objetivos específicos

1. Caracterizar fenotípicamente pacientes e comunicantes quanto a anticorpos anti-LID-

NDO;

2. Investigar loci adicionais associados a susceptibilidade à hanseníase e/ou reações

hansênicas utilizando Immunochip;

3. Avaliar a associação de SNPs à doença/infecção, reações hansênicas e quanto a

produção de anticorpos nos indivíduos estudados.

Page 30: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

28 3 MATERIAIS E MÉTODOS

3.1 População e desenho de estudo

O desenho de estudo foi observacional analítico tipo caso-controle. As análises

genéticas foram realizadas utilizando uma amostragem da população do RN constituída por

828 pessoas. Grupos foram separados pela cidade de moradia, forma clínica de acordo com a

classificação de Madrid, na qual as formas borderline são referidas como “dimorfa”, e estado

reacional. O DNA genômico foi obtido e as amostras foram genotipadas utilizando uma

matriz de SNPs associados a doenças auto-imunes (Immunochip). Dos 828 indivíduos

recrutados, amostras de DNA de 813 participantes foram genotipadas, das quais 13 eram

duplicados e 1 amostra foi excluída por falta de informação, obtendo então um total de 799

amostras (Figura 4). Os genótipos obtidos foram revisados em um controle de qualidade

(ANDERSON et al., 2011), e um total de 702 amostras (600 casos e 102 comunicantes)

selecionados para análises posteriores nas ferramentas PLINK v 1.7 (PURCELL et al., 2007)

e R (Figura 5).

Figura 4: População do Rio Grande do Norte incluída no estudo separadas em comunicantes e casos.

Page 31: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

29

Figura 5: Esquema da metodologia de desenvolvimento do estudo, dividida em caracterização fenotípica e

genotípica da população recrutada.

A amostra populacional do grupo caso foi constituída por moradores de áreas

endêmicas do Rio Grande do Norte. Um total de 714 pacientes foram recrutados em diferentes

centros de saúde entre 2003 e 2015. As amostras foram coletadas nas cidades de Natal

(Hospital Giselda Trigueiro – HGT e Hospital Universitário Onofre Lopes – HUOL),

Mossoró (Centro Clínico Prof Vingt Un Rosado) e Nova Cruz (Centro Clínico Integrado).

Foram incluídos nesse estudo aproximadamente 284 pacientes no momento do diagnóstico,

260 pacientes em tratamento para hanseníase (PQT), e 170 pacientes que haviam concluído o

PQT, mas estavam em tratamento para reações hansênicas, independente da forma clínica.

Todos os pacientes tiveram o diagnóstico confirmado por médico especialista seguindo os

critérios da OMS. Foram obtidas informações relativas à idade, sexo, exame físico, situação

Page 32: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

30 socioeconômica, forma clínica da hanseníase, data do diagnóstico, ocorrência de episódios

reacionais, índice baciloscópico e resultado de biópsia cutânea, todos resgatados dos

prontuários hospitalares localizados onde o paciente havia se tratado. Os pacientes foram

acompanhados quanto ao aparecimento de episódios reacionais durante o tratamento e até 1

ano após o fim da PQT. Todos os pacientes foram orientados quanto aos sintomas das

reações, e aconselhados a retornar ao local de tratamento.

Para a amostra populacional de comunicantes (grupo controle), foram incluídos

moradores saudáveis de domicílios com histórico de hanseníase e moradores de domicílios

vizinhos. As amostras foram coletadas nas cidades de Mossoró e Nova Cruz. Foram obtidas

informações relativas à idade, gênero, exame físico e socioeconômicas, seguindo protocolo

idêntico ao aplicado para o grupo de casos.

Após o recrutamento e caracterização clínica, o participante da pesquisa foi submetido a

uma coleta de aproximadamente 10 mL de sangue periférico, seguido os padrões de assepsia,

após a qual o material biológico foi acondicionado em tubos com EDTA, para a obtenção de

DNA, e tubos sem anticoagulante, para a obtenção de soro. O material biológico coletado foi

enviado ao Laboratório de Imunogenética, localizado no Instituto de Medicina Tropical (IMT)

na UFRN para processamento. O soro foi utilizado para análises sorológicas com ELISA, e o

DNA foi utilizado para genotipagem utilizando Immunochip.

3.2 Considerações éticas

Todos os indivíduos foram instruídos a respeito dos objetivos do estudo através de um

formulário de consentimento informado. O termo de consentimento e o protocolo de estudo

foram avaliados e aprovados pelo Comitê de Ética em Pesquisa da Universidade Federal do

Rio Grande do Norte, bem como pela Comissão Nacional de Ética em Pesquisa (CEP-UFRN

145/05; CONEP 12504, CAAE 0042.0.051.051-09).

3.3 Determinação de anticorpo anti-LID-NDO

O soro foi obtido após a centrifugação de sangue total, sendo estas amostras

amarzenados a temperatura inferior a -20ºC.

Do total de indivíduos selecionados para o estudo, 763 foram analisados quanto a

presença de anticorpos para LID-NDO a partir de Elisa, seguindo protocolo estabelecido em

Page 33: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

31 DUTHIE et al., 2014. Inicialmente, foi adicionado 1µg/mL do antígeno recombinante LID-

NDO ou controle negativo (apenas tampão) a placas de 96 poços Polysorp® Nunc (Thermo

Fisher Scientific, NY, USA), que foram incubadas a 4 ºC durante a noite. Após a incubação

foi realizado o bloqueio com PBS-Tween (PBST) contendo 1% de albumina bovina sérica

(BSA) (Sigma-Aldrich, St. Louis, MO, USA) à temperatura ambiente, por 1 hora. As

amostras de soro foram diluídas em PBST contendo 0,1% de BSA (1:200) e adicionadas às

placas em duplicata. Os poços foram lavados após incubação por 2 horas em temperatura

ambiente e foi adicionado o conjugado anticorpo-peroxidase (anti-IgG-HRP e anti-IgM-HRP)

(Rockland Immunochemicals, Gilbertsville, PA, USA). Os poços foram novamente lavados

após 1 hora de incubação à temperatura ambiente, aos quais foram adicionados

tetrametilbenzidina (TMB) (Kirkegaard and Perry Laboratories, Gaithersburg, MD, USA). A

reação foi parada após 15 minutos pela adição de ácido sulfúrico 1 N. A leitura da

absorbância das placas foi realizada utilizando o comprimento de onda de 450 nm. A

densidade optica (OD) de cada amostra foi expressa como a OD média de suas duplicadas

subtraída da placa controle negativo. Duas amostras controle foram adicionadas em cada

placa.

A quantificação dos anticorpos anti-LID-NDO foi realizada no ao Laboratório de

Imunogenética, localizado no IMT-UFRN.

3.4 Extração de DNA e genotipagem

O DNA genômico foi obtido a partir de células do sangue periférico por meio de

protocolo salting-out (MILLER, DYKES e POLESKY, 1988) e em seguida quantificado e

estocado a temperatura inferior a -20ºC.

Posteriormente, 200 ng de DNA de 828 indivíduos, sendo 714 pacientes e 114

comunicantes, foram adicionados a uma matriz de SNPs associados a doenças auto-imunes

(Infinium ImmunoArray BeadChip-Illumina). Dessas amostras, 813 foram genotipadas de

acordo com as instruções do fabricante no laboratório de referência CIDR (Center for

Inherited Disease Research), localizado na Universidade Johns Hopkins, situada em

Baltimore, Maryland, nos Estados Unidos.

A figura 6 mostra o protocolo estabelecido para a genotipagem por Immunochip. A

matriz é constituída por diversas esferas microscópicas de sílica codificadas com cópias

múltiplas de uma sonda de oligonucleotídeos ligadas e que têm como alvo um locus

Page 34: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

32 especifico no genoma. As amostras de DNA são amplificadas, incubadas, fragmentadas e

adicionadas à matriz (passos 1 a 5 na figura 6). Quando o DNA fragmentado passa sobre a

esfera microscópica, cada fragmento se ligará por complementaridade de bases a uma sonda,

parando uma base antes do locus de interesse (passo 6 na figura 6). Uma base é adicionada

enzimaticamente e os produtos são marcados com fluoróforos (passo 7 na figura 6), que

emitem sinais detectados pelo scanner, gerando uma imagem, na qual é possível analisar e

gerar relatórios referentes a relação alélica de cada locus.

Figura 6: Protocolo de genotipagem por Immunochip, figura adaptada de https://www.illumina.com/science/technology/beadarray-technology/infinium-assay.html.

Page 35: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

33 3.5 Controle de qualidade

Em estudos de associação caso-controle, com um alto número de marcadores, é

necessário realizar várias etapas de controle de qualidade para remover indivíduos e

marcadores com altas taxas de erro que introduzem viés ao estudo, pois poderiam levar a um

aumento de associações falso-negativas e falso-positivas.

Os passos realizados no controle de qualidade são primordiais para o sucesso do estudo

e são necessários antes de testar estatisticamente para associação. Porém, o impacto da

remoção de um marcador é maior do que a remoção de um indivíduo, pois cada marcador

removido de um estudo é uma associação em potencial que seria ignorada. Dessa forma, o

controle de qualidade por indivíduo é realizada antes do controle de qualidade por marcador,

evitando que marcadores sejam erroneamente removidos (ANDERSON et al., 2011).

3.5.1 Controle de qualidade por indivíduo

Após a genotipagem, o laboratório de referência forneceu um arquivo contendo um total

de 831 amostras, das quais 13 eram duplicatas cegas (o grupo que realizou a genotipagem não

conhecia) e 18 eram controles positivos. Para realizar o controle de qualidade, foram

removidas as amostras duplicadas, os controles positivos e uma amostra devido a

inconsistência de identificação. Dessa forma, realizamos o controle de qualidade em 799

amostras, obtendo um resultado de 702 amostras que passaram pelo controle de qualidade e

foram utilizadas nas análises de associação.

O controle de qualidade por indivíduo consiste em quatro passos realizados nos

programas PLINK e R, sumarizado na figura 7:

1. Identificação de indivíduos com informação de sexo discordante;

2. Identificação de indivíduos com genótipos faltantes e taxa de heterozigosidade;

3. Identificação de indivíduos duplicados ou relacionados;

4. Identificação de indivíduos com ancestralidade divergentes.

Page 36: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

34

Figura 7: Esquema do controle de qualidade por-indivíduo.

O primeiro passo no controle qualidade foi a identificação de indivíduos com

informação de sexo discordante, na qual a melhor maneira é calcular a taxa de

homozigosidade em todos os SNPs de cromossomo X para cada indivíduo na amostra e

compará-los com a taxa esperada, é esperado que amostras do sexo masculino tenham uma

taxa de homozigosidade maior que 0,8, e amostras do sexo feminino tenham uma taxa de

homozigosidade menor que 0,2. Ao utilizar a função –check-sex no PLINK identificamos 61

problemas ilustrados na figura 8.

Dos 61 problemas identificados, temos:

• 1 amostra codificada como sexo feminino foi genotipada como sexo masculino;

• 1 amostra codificada como desconhecida foi genotipada como sexo masculino;

• 2 amostras codificadas como sexo masculino foram genotipadas como sexo

feminino;

• 57 amostras codificadas como sexo feminino foram genotipadas como

inconclusivas, mas não foram excluídas da análise.

Page 37: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

35

Figura 8: Informação de sexo discordante. O eixo Y (F) é a medida de homozigosidade do cromossomo X, e o

eixo X é o sexo assinalado da amostra (0: desconhecido, 1: sexo masculino, 2: sexo feminino).

O segundo passo do controle de qualidade foi a identificação de indivíduos com

genótipos faltantes e taxa de heterozigosidade. Esta análise consistiu na exclusão de amostras

com DNA de baixa qualidade. 14 amostras foram excluídas nessa etapa, sendo 3 amostras

com proporção de ausência de genótipos ≥ 0,05, e 11 amostras com a taxa de

heterozigosidade menor que o limite determinado em três desvios padrões da média, como

ilustrados na figura 9.

Page 38: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

36

Figura 9: Ausência de genótipos e taxa de heterozigosidade. O eixo Y representa a taxa de heterozigosidade, e o eixo X a proporção de genótipos faltantes. As linhas pontilhadas representam os limites.

O terceiro passo do controle de qualidade consistiu na identificação de indíduos

duplicados ou relacionados, sendo de grande importância, porque se duplicados ou parentes

estão presentes na análise, um viés pode ser introduzido, já que os genótipos dentro das

famílias estarão sobre-representados. Para identificar indivíduos duplicados ou relacionados,

IBS (identidade por estado) e o IBD (identidade por descendência) em pares para todos os

pares de indivíduos foram calculados com base em um conjunto de 54.882 SNPs

independentes (r2 <0,2), usando a função --genome em PLINK.

O resultado desta análise está representado na figura 10. Dessa forma, 83 indivíduos

foram excluídos para que o estudo contenha somente pessoas geneticamente não-relacionada,

o que corresponde a mais de 10% das amostras. Porém o valor alto foi devido à forma com

que as informações acerca de parentesco foram obtidas. Além disso, principalmente em

Mossoró, muitas pessoas não sabiam que eram parentes.

Page 39: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

37

Figura 10: Indivíduos duplicados ou relacionados. O eixo Y representa a proporção de alelos IBD, e o eixo X o

Pr, os pares de alelos. Cada ponto indica um par de indivíduos e linhas tracejadas representam valores de IBD

esperados para diferentes pares relativos (Dup: duplicado; P-Off: descendência parental; Sib: irmãos; HS: meio-

irmãos; Avu: avuncular; FC: primos; Unr: não relacionado).

O quarto passo do controle de qualidade, representado na figura 11, consistiu na

identificação de indivíduos com ancestralidades divergentes, que devem ser removidos para

que os casos e controles incluídos no estudo possuam backgrounds étnicos os mais próximos

possíveis. Dessa forma, foi realizada uma análise de componentes principais com base em

distâncias pares do IBS. Cinco indivíduos que se encontram abaixo do limite foram

removidos.

Page 40: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

38

Figura 11: Avaliação da estratificação populacional, scatter plot dos dois primeiros componentes principais da

avaliação. A linha tracejada em vermelho indica o limite estabelecido.

3.5.2 Controle de qualidade por marcador

Após a realização do controle de qualidade por indivíduo, e da exclusão de 97

amostras, foi realizado o controle de qualidade por marcador no PLINK utilizando 702

amostras (600 casos e 102 controles) e 196.524 marcadores.

O controle de qualidade por marcador consistiu em quatro passos (figura 12):

1. Identificação de SNPs com excesso de genótipo faltante;

2. Identificação de SNPs com um desvio significativo no equilíbrio de Hardy-

Weinberg (HWE);

3. Identificação de SNPs com taxas de genótipo faltantes significativamente

diferentes entre casos e controles;

4. Remoção de marcadores com uma frequência de alelos menores muito baixa.

Page 41: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

39

Figura 12: Esquema do controle de qualidade por marcador.

O primeiro passo do controle de qualidade por marcador consistiu na identificação e

remoção de marcadores com excesso de genótipo faltante, pois esses poderiam representar

falso-positivos e reduzir a habilidade de identificar associações. Dessa forma, marcadores

com uma taxa de chamada inferior a 95% foram removidos do estudo.

O segundo passo compreendeu na exclusão de marcadores com desvio no equilíbrio

de Hardy-Weinberg, porque isso é indicativo de erro na genotipagem ou na chamada de

genótipos. Porém, desvios no HWE podem indicar seleção de um alelo em relação a outro,

então somente amostras de controles são utilizadas para testar desvios em HWE.

O terceiro passo foi a identificação e exclusão de SNPs com taxas de genótipo

faltantes significativamente diferentes entre casos e controles, sendo outro meio de reduzir

confudimento e remoção de SNP pouco genotipados.

O quarto e último passo do controle de qualidade por marcador consistiu na remoção

de marcadores com uma frequência de alelos menores muito baixa (MAF < 0,01), importante

pois são variantes difíceis de serem chamados com algoritmos de chamada utilizados

atualmente, apresentando-se como falso positivos.

Após a realização dos passos, 49.695 marcadores foram excluídos do estudo, obtendo,

então, um total de 146.829 marcadores.

Page 42: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

40 3.6 Análise estatística

Gráficos e valores médios foram gerados utilizando GraphPad Prism v.7.

Como o Immunochip apresenta centenas de milhares de marcadores, é importante

estabelecer padrões estatísticos para distinguir entre sinais causais e sinais falsos. Para os

estudos genéticos, as análises de associação caso-controle foram realizadas através do modelo

de regressão logística utilizando o programa PLINK v.107 para estimar a chance de o

indivíduo ser doente. A associação dos SNPs nos grupos foi testada utilizando o teste de

Cochran-Armitage corrigido para sexo e idade. Os testes compararam frequências dos alelos

dos SNPs com:

• Hanseníase vs. Comunicantes/ PB vs. Comunicantes/ MB vs. Comunicantes

• Reação vs Não-reação/ Reação Tipo 1 vs Não-reação/ Reação Tipo 2 vs Não-

reação

A associação de freqüências alélicas com taxa de anticorpos LID-NDO também foi

analisada, utilizando modelo de regressão linear.

O nível de significância (α) amplamente usado na ciência é de 5% de erro, esta taxa foi

corrigida pela correção de Bonferroni, na qual divide a taxa pelo número de hipóteses, que

seriam, nesse caso, o número de marcadores independentes, ou seja: p-valor = α/m =

0,05/54.945 = 0,00000091.

Dessa forma, consideramos um p-valor de 9,10 x 10-7 como limiar para associação

significativa.

Page 43: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

41 4 RESULTADOS

4.1 Características populacionais

A população estudada foi dividida entre casos de hanseníase e comunicantes, cujos

dados estão representados na tabela 2. A média de idade dos casos foi de 48,2 anos, maior

quando comparada à média dos comunicantes de 41,5 anos (p=0,0004). Houve um

predomínio de indivíduos do sexo masculino (55,6%) nos casos, ao contrário dos

comunicantes, no qual predomínou o sexo feminino (63,5%). Nos casos, a maioria dos

indivíduos residiam em Natal (75,3%), os comunicantes foram coletados em sua maioria em

Mossoró (92,8%). Entre os casos, houve um predomínio de formas multibacilares (68,9%),

forma clínica dimorfa (44%) e de pacientes que desenvolveram algum tipo de reação

hansênica (57,9%), sendo a maioria a reação do tipo I (41,3%). Alguns pacientes não

apresentaram algumas informações, como cidade, se desenvolveram reação e qual tipo de

reação desenvolveram (1, 16 e 36 pacientes sem estas informações, respectivamente), porém a

falta destas informações não irá influenciar muito nas análises genéticas.

Tabela 2: Características da população de estudo.

Variável Casos (n = 673) Comunicantes (n = 126) p-valor

Média de idade, anos Sexo, n (%)

Masculino

48,2 (±17) 41,5 (±21) 0,00041

374 (55,6) 46 (36,5) <0,00012

Feminino 299 (44,4) 80 (63,5)

Cidade, n (%)

Natal 507 (75,3) 6 (4,8) <0,00013

Mossoró 145 (21,5) 117 (92,8)

Nova Cruz 21 (3,2) 3 (2,4)

Classificação Operacional, n (%)

PB 208 (30,9) - -

MB 464 (68,9) -

Forma Clínica, n (%)

Indeterminada 13 (2) - -

Tuberculóide 177 (26,3) -

Dimorfa 296 (44) -

Lepromatosa 184 (27,3) -

Reação Hansênica, n (%)

Não teve reação 267 (39,7) - -

Teve reação 390 (57,9) -

Tipo I 161 (41,3) -

Tipo II 135 (34,6) -

Neurite 58 (14,9) -

A idade está representada com a media (± desvio padrão). 1 Teste t de Student; 2 Teste exato de Fisher; 3 Teste Chi Quadrado

Page 44: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

42 4.2 Determinação dos níveis de anticorpo anti-Mycobacterium leprae

Do total de indivíduos selecionados para o estudo, 763 foram analisados quanto a

presença de anticorpos LID-NDO para o M. leprae, sendo 637 casos e 126 controles (Figura

13). Os níveis de anticorpos variaram de acordo com a classificação operacional do caso (PB

ou MB) e o tipo de reação, sendo mais elevada em casos MB e reação do tipo II (ENH)

(Figuras 13.A e 13.B, respectivamente).

Figura 13: Análise comparativa entre os níveis de anticorpos específicos LID-NDO entre comunicantes e casos.

A densidade óptica (OD) para cada amostra foi obtida após subtrair o valor da OD obtida com o placebo. A)

Comparação entre comunicantes, casos PB e MB, B) Comparação entre diferentes tipos de reação. HHC:

comunicantes; PB: paucibacilares; MB: multibacilares; RR: Reação Hansênica; ENH: Eritema Nodoso

Hansênico.

Page 45: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

43 4.3 Análise genética

Do total de indivíduos selecionados para o estudo, 702 (600 casos e 102

comunicantes) foram analisados quanto à associação de SNPs.

As análises foram separadas em três grupos:

• Hanseníase vs comunicantes

• Reação vs não-reação

• Taxa de anticorpos LID-NDO

4.3.1 Hanseníase vs. Comunicantes

Primeiramente foram analisados os diferentes SNPs associados a doença ou infecção

subclínica utilizando amostras de casos de hanseníase per se e amostras de comunicantes,

respectivamente. Os resultados da análise estão demonstrados na figura 14, na qual os dez

SNPs mais fortemente associados à Hanseníase estão listados na tabela 3 e na tabela 6.

Na figura 14 foram destacados os SNPs mais fortemente associados com os subtipos

de hanseníase, MB e PB, comparados com comunicantes. Os três SNPs mais associados com

MB se encontram dentro das regiões intergênicas FASLG-TNFSF18 (OR=0,1359,

p=1,276x10-6), e LOC152663-MGC39584 (OR=0,2753, p=3,079x10-5), e do gene VDR

(OR=0,2355, p=9,472x10-6). Os três SNPs mais associados com PB se encontram dentro do

gene CTNNA2 (OR=2,4720, p=2,449x10-5) e das regiões intergênicas CLINT1-EBF1

(OR=0,4447, p=3,070x10-5) e STEAP4-ZNF804B (OR=0,3529, p=5,431x10-6).

Page 46: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

44

Figura 14: Manhattan plot demonstrando p valores no genoma total em associação com Hanseníase. O eixo y representa os p valores em –log10 de 146.829 SNPs e suas posições cromossômicas estão mostrados no eixo x. A linha horizontal em azul demonstra o limite de p< 1 x 10-4. Os genes em vermelho representam os três SNPs mais associados com Hanseníase MB, e os genes em verde representam os três SNPs mais associados com Hanseníase PB.

Page 47: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

45 Tabela 3: Associações com Hanseníase para 10 SNPs de susceptibilidade, de acordo com análise.

SNP Crm Posição Alelo menor Gene MAF_A MAF_U p-valor OR (L95-U95)

imm_1_171206127 1 171206127 G FASLG - TNFSF18 0,02667 0,07843 1,45x10-5 0,2229 (0,1131-0,4393)

1kg_5_102565099 5 102565099 C HISPPD1 0,3917 0,5294 9,18x10-5 0,5275 (0,3828-0,7268)

rs184946 5 158079394 A CLINT1 - EBF1 0,3158 0,4853 6,57x10-5 0,5190 (0,3761-0,7162)

rs727098 6 139911322 G LINC01625 - LOC100132735 0,2733 0,3824 7,22x10-5 0,4927 (0,3473-0,6988)

rs6465973 7 77671200 G MAGI2 0,2458 0,1324 6,62x10-5 2,741 (1,67-4,498)

rs2538909 7 88271503 A STEAP4 - ZNF804B 0,185 0,3088 6,32x10-6 0,4338 (0,3019-0,6233)

rs11251836 10 3335447 A PITRM1 - KLF6 0,09917 0,201 5,46x10-5 0,3965 (0,2530-0,6214)

rs10783218 12 48272743 A VDR 0,04333 0,1225 9,89x10-6 0,2825 (0,1613-0,4949)

rs3829948 14 92186579 C RIN3 0,145 0,2451 5,84x10-5 0,4326 (0,2875-0,6509)

imm_21_44431174 21 44431174 A LOC100129890 - ICOSLG 0,0425 0,1029 6,76x10-5 0,3033 (0,1686-0,5454)

Crm, cromossomo; MAF_A, alelo de menor frequência nos casos; MAF_U, alelo de menor frequência nos comunicantes; OR, odds ratio

Page 48: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

46 4.3.2 Reação vs. Não reação

Da mesma forma, foram avaliados os SNPs associados a reações hansênicas utilizando

amostras de casos com reações hansênicas e casos que não as desenvolveram.

Os resultados da análise estão demonstrados na figura 15, na qual os SNPs mais

fortemente associados a reações estão listados na tabela 4 e na tabela 6.

Na figura 15 foram destacados os SNPs associados mais fortemente com as reações

hansênicas de tipo I e tipo II em comparação com casos sem reação. Os cinco SNPs mais

fortemente associados com reação do tipo I estão dentro do gene UBD (OR=3,0380,

p=7,520x10-6) e da região intergênica FLJ45139-LOC391282 (OR=0,4953, p=4,331x10-6). Os

cinco SNPs mais fortemente associados com reação do tipo II estão dentro da região

intergênica MYADML-CRIM1 (OR=2,2750, p=3,287x10-5) e dos genes CCHCR1

(OR=2,3280, p=1,192x10-5) e SDSL (OR=4,4830, p=1,340x10-5).

Figura 15: Manhattan plot demonstrando p valores no genoma total em associação com reações hansênicas. O eixo y representa os p valores em –log10 de 146.829 SNPs e suas posições cromossômicas estão mostrados no eixo x. A linha horizontal em azul demonstra o limite de p<1 x 10-4. Os genes em vermelho representam os cinco SNPs mais associados com reação do tipo II, e os genes em verde representam os cinco SNPs mais associados com reação do tipo I.

Page 49: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

47

Na figura 15 observamos vários SNPs com significância na mesma região no

cromossomo 6, dessa forma, destacamos essa região (6p21.33) na figura 16. As cores dos

SNPs indicam o nível de LD (desequilíbrio de ligação) com os SNPs mais fortemente

associados de acordo com uma escala de 0 a 1 baseada em valores de r2 da população

europeia do projeto 1000 Genomes. No painel inferior foram adicionadas anotações de genes

do navegador de genoma da Universidade da Califórnia, Santa Cruz, mostrando os genes

próximos.

Figura 16: Parcelas de associação regional em 6p21.33. Os valores de p de SNPs genotipados são plotados como valores de -log10 contra suas posições físicas cromossômicas (hg19). Taxas de recombinação estimadas da população europeia de 1000 Genomes mostram a estrutura LD local. As cores dos SNPs indicam LD com os SNPs mais associados de acordo com uma escala de r2=0 a r2=1 baseada em valores de r2 da população europeia de 1000 Genomes. No painel inferior, anotações de genes do navegador de genoma da Universidade da Califórnia, Santa Cruz.

Page 50: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

48 Tabela 4: Associações com reações hansênicas para 30 SNPs de susceptibilidade, de acordo com análise. SNP Crm Posição Alelo menor Gene MAF_A MAF_U p-valor OR (L95-U95)

rs2193790 2 34299457 A MYADML - CRIM1 0,2845 0,1885 9,685x10-5 1,775 (1,33 - 2,368)

rs9501677 6 29497002 A LOC729653 - UBD 0,1569 0,07787 6,034x10-5 2,304 (1,533 - 3,465)

rs10947130 6 31013746 A MUC22 - HCG22 0,1452 0,2439 7,226x10-5 0,5489 (0,4082 - 0,7381)

rs10947131 6 31015936 A MUC22 - HCG22 0,1642 0,2766 1,311x10-5 0,5317 (0,4002 - 0,7064)

rs9295947 6 31019460 A MUC22 - HCG22 0,1628 0,2766 1,052x10-5 0,5277 (0,3971 - 0,7012)

rs2394427 6 31019562 A MUC22 - HCG22 0,1628 0,2766 1,052x10-5 0,5277 (0,3971 - 0,7012)

rs4713429 6 31021017 C MUC22 - HCG22 0,2082 0,3156 8,542x10-5 0,5803 (0,4423 - 0,7613)

rs9263726 6 31106499 A PSORS1C1 0,1261 0,05533 5,595x10-5 2,56(1,621 - 4,045)

rs2233945 6 31107361 A PSORS1C2 - PSORS1C1 0,1408 0,06762 5,777x10-5 2,405 (1,568 - 3,688)

rs12364 6 31110786 A CCHCR1 0,1569 0,075 3,584x10-5 2,359 (1,57 - 3,545)

rs9263739 6 31111356 A CCHCR1 0,1525 0,07377 2,945x10-5 2,42 (1,599 - 3,664)

rs9263745 6 31113326 A CCHCR1 0,1408 0,06557 3,341x10-5 2,497 (1,621 - 3,847)

rs2240064 6 31114573 G CCHCR1 0,456 0,3115 1,892x10-7 2,009 (1,545 - 2,612)

rs3131012 6 31115441 A CCHCR1 0,4604 0,3197 3,929x10-7 1,958 (1,51 - 2,538)

rs2073717 6 31122126 G CCHCR1 0,4032 0,2582 1,485x10-7 2,037 (1,562 - 2,657)

rs130077 6 31122330 A CCHCR1 0,1408 0,06762 5,777x10-5 2,405 (1,568 - 3,688)

rs2239524 6 31125569 A CCHCR1 0,3636 0,2418 2,496x10-6 1,921 (1,464 - 2,522)

rs9263785 6 31125823 C CCHCR1 0,1408 0,06762 5,777x10-5 2,405 (1,568 - 3,688)

rs9263800 6 31134599 A POU5F1 0,1408 0,06762 5,777x10-5 2,405 (1,568 - 3,688)

rs9263810 6 31137632 C POU5F1 0,1408 0,06762 5,777x10-5 2,405 (1,568 - 3,688)

rs3094193 6 31137774 C POU5F1 0,368 0,25 5,418x10-6 1,86 (1,424 - 2,43)

rs879882 6 31139452 A POU5F1 – PSORS1C3 0,4399 0,3238 5,995x10-6 1,82 (1,404 - 2,359)

rs4713447 6 31162963 G PSORS1C3 - HCG27 0,4487 0,5697 9,04x10-5 0,6118 (0,4784 - 7,824)

rs9501077 6 31167512 G HCG27 0,4487 0,5697 9,04x10-5 0,6118 (0,4784 - 0,7824)

rs12662501 6 31190850 A HCG27 - HLA-C 0,1569 0,2541 7,69x10-5 0,5525 (0,4117 - 0,7414)

rs1063646 6 31215627 A PSORS1C1 0,1529 0,07377 2,761x10-5 2,428 (1,604 - 3,676)

rs9366768 6 31276473 A HCG27 0,4487 0,5697 9,04x10-5 0,6118 (0,4784 - 0,7824)

rs4711658 6 41270869 G TREML2 0,5396 0,4016 3,589x10-6 1,773 (1,392 - 2,259)

imm_10_49794086 10 49794086 A LRRC18 - WDFY4 0,4457 0,3258 6,923x10-5 1,666 (1,296 - 2,142)

ccc-21-39331475-A-G 21 39331475 G FLJ45139 - LOC391282 0,1003 0,1803 8,378x10-5 0,4898 (0,3432 - 0,699)

Crm, cromossomo; MAF_A, alelo de menor frequência nos casos com reação; MAF_U, alelo de menor frequência nos casos sem reação; OR, odds ratio.

Page 51: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

49 4.3.3 Nível de anticorpos anti-LID-NDO

Também foram avaliados os SNPs associados a resultados contínuos de taxa de

anticorpos LID-NDO utilizando amostras de 675 indivíduos, sendo 574 casos e 101

comunicantes.

Os resultados da análise estão demonstrados na figura 17, na qual os SNPs mais

fortemente associados estão listados na tabela 5 e na tabela 6.

Figura 17: Manhattan plot demonstrando p valores no genoma total em associação com taxa de anticorpos LID-NDO. O eixo y representa os p valores em –log10 de 146.829 SNPs e suas posições cromossômicas estão mostrados no eixo x. A linha horizontal em azul demonstra o limite de p<1 x 10-5 e a linha horizontal em vermelho representa o limite de p<1x10-4.

Todos os SNPs associados nos três grupos analisados foram sumarizados na tabela 6, que lista esses SNPs, em qual cromossomos se encontram, o grupo resposta, seu p-valor, coeficientes de risco (OR e β), genes relacionados e região intergênica.

Page 52: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

50 Tabela 5: Associações com taxas de anticorpos LID-NDO para 15 SNPs de susceptibilidade, de acordo com análise.

SNP Crm Posição Alelo menor Gene p-valor β (L95-U95)

rs2252865 1 8422676 A RERE 5,3x10-5 -0,07127 (-0,1056 - -0,03694)

rs10779702 1 8423510 A RERE 3,826x10-5 -0,07162 (-0,1055 - -0,03776)

rs894875 1 8432136 G RERE - DKFZp566H0824 3,142x10-5 -0,07273 (-0,1067 - -0,03873)

rs6678140 1 8436802 G RERE - DKFZp566H0824 5,573x10-5 -0,06956 (-0,1032 - -0,03595)

rs2708633 1 8439625 C RERE - DKFZp566H0824 2,315x10-5 -0,07363 (-0,1075 - -0,03977)

rs13035389 2 199666111 A PLCL1 - SATB2 8,662x10-5 -0,05989 (-0,08961 - -0,03017)

1kg_5_173187315 5 173187315 G LOC729170 - CPEB4 5,692x10-5 0,1336 (0,06899 - 0,1983)

1kg_5_173187876 5 173187876 G LOC729170 - CPEB4 5,692x10-5 0,1336 (0,06899 - 0,1983)

1kg_5_173188023 5 173188023 G LOC729170 - CPEB4 5,692x10-5 0,1336 (0,06899 - 0,1983)

imm_6_127396776 6 127396776 G LOC728666 - RSPO3 8,454x10-5 -0,1553 (-0,2323 - -0,07837)

imm_6_128252771 6 128252771 A THEMIS 5,003x10-5 -0,07685 (-0,1137 - -0,03995)

rs1613984 10 44460433 A LINC00841 2,168x10-6 -0,1407 (-0,1984 - -0,08299)

imm_11_691201 11 691201 T TMEM80 4,38x10-5 0,09417 (0,04931 - 0,139)

1kg_11_113787520 11 113787520 A RBM7 - REXO2 9,089x10-5 -0,123 (-0,1842 - -0,06178)

rs2877400 18 56810775 A SEC11C 9,069x10-6 0,09168 (0,05151 - 0,1319)

Crm, cromossomo; β, coeficiente beta.

Page 53: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

51 Tabela 6: 55 SNPs associados à hanseníase, dentro de 36 regiões gênicas e intergênicas, de acordo com análise.

SNP Crm Grupo p-valor OR β Gene Região

imm_1_171206127 1 Casos vs comunicantes 1,45x10-5 0,2229 - FASLG | TNFSF18 Intergênica

1kg_5_102565099 5 Casos vs comunicantes 9,18x10-5 0,5275 - HISPPD1 Intron

rs184946 5 Casos vs comunicantes 6,57x10-5 0,519 - CLINT1 | EBF1 Intergênica

rs727098 6 Casos vs comunicantes 7,22x10-5 0,4927 - LOC645434 |LOC100132735 Intergênica

rs6465973 7 Casos vs comunicantes 6,62x10-5 2,741 - MAGI2 Intron

rs2538909 7 Casos vs comunicantes 6,32x10-6 0,4338 - STAP4 | ZNF804B Intergênica

rs11251836 10 Casos vs comunicantes 5,46x10-5 0,3965 - PITRM1 | KLF6 Intergênica

rs10783218 12 Casos vs comunicantes 9,89x10-6 0,2825 - VDR Intron

rs3829948 14 Casos vs comunicantes 5,84x10-5 0,4326 - RIN3 Intron

imm_21_44431174 21 Casos vs comunicantes 6,76x10-5 0,3033 - LOC100129890 | ICOSLG Intergênica

rs2193790 2 Reação vs não-reação 9,69x10-2 1,775 - MYADML | CRIM1 Intergênica

rs9501677 6 Reação vs não-reação 6,03x10-2 2,304 - LOC729653 | UBD Intergênica

rs10947130 6 Reação vs não-reação 7,23x10-2 0,5489 - LOC729792 | HCG22 Intergênica

rs10947131 6 Reação vs não-reação 1,31x10-2 0,5317 - LOC729792 | HCG22 Intergênica

rs9295947 6 Reação vs não-reação 1,05x10-2 0,5277 - LOC729792 | HCG22 Intergênica

rs2394427 6 Reação vs não-reação 1,05x10-2 0,5277 - LOC729792 | HCG22 Intergênica

rs4713429 6 Reação vs não-reação 8,54x10-2 0,5803 - LOC729792 | HCG22 Intergênica

rs9263726 6 Reação vs não-reação 5,60x10-2 2,560 - PSORS1C1 Codante

rs2233945 6 Reação vs não-reação 5,78x10-2 2,405 - PSORS1C2 | PSORS1C1 Intergênica

rs12364 6 Reação vs não-reação 3,58x10-2 2,359 - CCHCR1 Codante

rs9263739 6 Reação vs não-reação 2,95x10-2 2,42 - CCHCR1 Intron

rs9263745 6 Reação vs não-reação 3,34x10-2 2,497 - CCHCR1 Intron

rs2240064 6 Reação vs não-reação 1,89x10-4 2,009 - CCHCR1 Intron

rs3131012 6 Reação vs não-reação 3,93x10-4 1,958 - CCHCR1 Intron

rs2073717 6 Reação vs não-reação 1,49x10-4 2,037 - CCHCR1 Intron

rs130077 6 Reação vs não-reação 5,78x10-2 2,405 - CCHCR1 Codante

rs2239524 6 Reação vs não-reação 2,50x10-3 1,921 - CCHCR1 Intron

Page 54: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

52

rs9263785 6 Reação vs não-reação 5,78x10-2 2,405 - CCHCR1 Intron

rs9263800 6 Reação vs não-reação 5,78x10-2 2,405 - POU5F1 Intron

rs9263810 6 Reação vs não-reação 5,78x10-2 2,405 - POU5F1 Intron

rs3094193 6 Reação vs não-reação 5,42x10-3 1,860 - POU5F1 Intron

rs879882 6 Reação vs não-reação 6,00x10-3 1,82 - POU5F1 | LOC100130889 Intergênica

rs4713447 6 Reação vs não-reação 9,04x10-5 0,6118 - LOC100130889 | HCG27 Intergênica

rs9501077 6 Reação vs não-reação 9,04x10-5 0,6118 - HCG27 Intron

rs12662501 6 Reação vs não-reação 7,69x10-5 0,5525 - HCG27 | HLA-C Intergênica

rs1063646 6 Reação vs não-reação 2,76x10-2 2,428 - PSORS1C1 Codante

rs9366768 6 Reação vs não-reação 9,04x10-5 0,6118 - HCG27 Intron

rs4711658 6 Reação vs não-reação 3,59x10-3 1,773 - TREML2 Intron

imm_10_49794086 10 Reação vs não-reação 6,92x10-2 1,666 - LRRC18 | WDFY4 Intergênica

ccc-21-39331475-A-G 21 Reação vs não-reação 8,38x10-2 0,4898 - FLJ45139 | LOC391282 Intergênica

rs2252865 1 Taxa de anticorpos 5,3x10-5 - -0,07127 RERE Intron

rs10779702 1 Taxa de anticorpos 3,826x10-5 - -0,07162 RERE Intron

rs894875 1 Taxa de anticorpos 3,142x10-5 - -0,07273 RERE - DKFZp566H0824 Intergênica

rs6678140 1 Taxa de anticorpos 5,573x10-5 - -0,06956 RERE - DKFZp566H0824 Intergênica

rs2708633 1 Taxa de anticorpos 2,315x10-5 - -0,07363 RERE - DKFZp566H0824 Intergênica

rs13035389 2 Taxa de anticorpos 8,662x10-5 - -0,05989 PLCL1 - SATB2 Intergênica

1kg_5_173187315 5 Taxa de anticorpos 5,692x10-5 - 0,1336 LOC729170 - CPEB4 Intergênica

1kg_5_173187876 5 Taxa de anticorpos 5,692x10-5 - 0,1336 LOC729170 - CPEB4 Intergênica

1kg_5_173188023 5 Taxa de anticorpos 5,692x10-5 - 0,1336 LOC729170 - CPEB4 Intergênica

imm_6_127396776 6 Taxa de anticorpos 8,454x10-5 - -0,1553 LOC728666 - RSPO3 Intergênica

imm_6_128252771 6 Taxa de anticorpos 5,003x10-5 - -0,07685 THEMIS Intron

rs1613984 10 Taxa de anticorpos 2,168x10-6 - -0,1407 LINC00841 Intergênica

imm_11_691201 11 Taxa de anticorpos 4,38x10-5 - 0,09417 TMEM80 Intron

1kg_11_113787520 11 Taxa de anticorpos 9,089x10-5 - -0,123 RBM7 - REXO2 Intergênica

rs2877400 18 Taxa de anticorpos 9,069x10-6 - 0,09168 SEC11C Intron

Crm, cromossomo; OR, odds ratio; β, coeficiente beta.

Page 55: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

53

5 DISCUSSÃO

O presente estudo foi dividido em duas partes: caracterização fenotípica e genotípica

da população exposta a M. leprae do Rio Grande do Norte, quanto a susceptibilidade à

hanseníase. A caracterização fenotípica foi realizada a partir da coleta e montagem de um

banco de dados referentes a pacientes e comunicantes, e pela detecção de anticorpos anti-LID-

NDO. Os resultados sorológicos também foram utilizados na análise genotípica, como um

grupo resposta para a detecção de SNPs associados a susceptibilidade à formação de

anticorpos anti-M. leprae.

Foi observada uma elevada produção de anticorpos em pacientes MB quando

comparados a PB, estando de acordo com o observado em diferentes estudos sorológicos

(AMORIM et al., 2016; BÜHRER-SÉKULA et al., 2003; FABRI et al., 2015).

A maioria dos indivíduos expostos a M. leprae desenvolvem resposta imune que

controla a proliferação do microorganismo e não desenvolvem hanseníase (SCOLLARD et

al., 2006), o que pode ser observado nos comunicantes, os quais o valor médio de anticorpos

não teve uma diferença significativa com relação aos pacientes PB.

A reatividade imune excessiva está envolvida com a patogênese (BARREIRO;

QUINTANA-MURCI, 2010). No nosso estudo observamos uma quantidade de anticorpos

mais elevada em pacientes MB, quando comparados com comunicantes e pacientes PB, e

naqueles que desenvolveram reação do tipo II, quando comparados com pacientes que não

desenvolveram reações, reação tipo I e neurite. Isso dá suporte a hipótese que excesso de

anticorpos pode ter um papel na patogênese da hanseníase e suas reações. Estudo de

expressão gênica com população oriunda do mesmo grupo etnico demonstrou que vias

relacionadas à imunidade estavam representadas nos perfis transcricionais de pacientes que

desenvolveram reações hansênicas tipo I e tipo II, com a via de complemento e coagulação

comum a ambos (DUPNIK et al., 2014).

A caracterização genotípica da população de estudo foi realizada a partir de

genotipagem por Immunochip, na qual foi feita a análise de SNPs associados a três grupos:

hanseníase vs comunicantes, reação vs não-reação, e taxa de anticorpos LID-NDO.

Primeiramente, analisamos os SNPs mais fortemente relacionados com doença

(casos) e infecção subclínica (comunicantes), utilizando uma abordagem caso/controle com

Page 56: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

54 regressão logística. Os dez SNPs mais fortemente associados foram listados, indicando o

cromossomo, sua posição, gene mais próximo, p-valor e Odds Ratio (OR). OR é uma

estimativa de risco relativo, na qual quando maior que 1 indica risco, e quando menor que 1

indica proteção à doença (STATISTICS SOLUTIONS, 2016).

Nesta análise, alguns SNPs foram associados a genes relacionados a resposta imune,

como FASLG (Fas ligante), TNFSF18 (Membro 18 da superfamília TNF), EBF1 (Fator

inicial 1 de células B) e ICOSLG (fator coestimulador indutivo de células T). Fas ligante é

um gene relacionado a regulação do sistema imune que codifica uma proteína

transmembrana indutora de apoptose pela ligação a Fas. Sua expressão aumentada em

macrófagos infectados por micobactérias indica um mecanismo de evasão que os protege de

linfócitos expressando Fas (MUSTAFA et al., 2001).

Também é interessante observar nessa análise que um SNP associado à hanseníase

está próximo do gene VDR (receptor de vitamina D). A vitamina D3 desempenha um papel

imunomodulador mediado através da ligação a VDR em monócitos, macrófagos e linfócitos

ativados, sua deficiência foi vista correlacionada negativamente com o aumento de citocinas

pró-inflamatórias encontradas em lesões de reação hansênica tipo II (LU’O’NG; HOÀNG

NGUYỄN, 2012). Um estudo de 2015 demonstrou que indivíduos com reações hansênicas

apresentaram níveis muito baixo de expressão de VDR associados a níveis baixos de

vitamina D3 séricos, tais indivíduos possuíam altos índices baciloscópicos e todos

apresentaram reação tipo II (MANDAL et al., 2015).

O receptor de vitamina D (VDR) possui vários SNPs localizados em locais de

restrição próximos a região 3’UTR (ARAÚJO et al., 2017), dessa forma a ocorrência de

hanseníase associada a polimorfismos em VDR tem sido investigada em várias populações,

tendo como justificativa o seu papel imunomodulador, além de regular a expressão de vários

genes envolvidos na homeostase do cálcio e fosfato, na proliferação e na diferenciação

celular (LU’O’NG; HOÀNG NGUYỄN, 2012; OLIVEIRA et al., 2017). Além disso, a

vitamina D é um estimulante essencial para a imunidade inata contra micobactérias, uma vez

que a ligação da vitamina D a VDR induz a síntese de peptídeos antimicrobianos, como a

catelicidina e β-defensinas, que destroem a membrana celular bacteriana (ARANOW, 2011).

Quando observamos os valores de OR dos SNPs associados à hanseníase, todos

apresentaram valores de OR menor que 1, indicando proteção contra o desenvolvimento de

doença, com exceção de um SNP relacionado ao gene MAGI2, que obteve um OR maior que

1, indicando risco à hanseníase. MAGI2 codifica uma proteína que age como molécula de

Page 57: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

55 andaime em junções sinápticas, montando receptores de neurotransmissores e proteínas de

adesão celular, podendo desempenhar um papel na regulação da sinalização mediada por

activina em células neuronais, não tendo um papel direto na hanseníase.

Além da análise hanseníase vs comunicantes, também foram avaliados os SNPs

associados à hanseníanse paucibacilar e multibacilar vs comunicantes. Todos SNPs

associados aos subgrupos também apresentaram um OR abaixo de 1, com exceção de

CTNNA2, associado à hanseníase PB. CTNNA pode funcionar como um agente de ligação

entre os receptores de adesão de caderina e o citoesqueleto para regular a adesão célula-

célula e diferenciação do sistema nervoso. Os SNPs associados a esse gene e o associado à

região intergênica LOC152663/MGC39584 em hanseníase MB não atigiram o p-valor

adotado na análise hanseníse vs comunicantes, mas quando cada subgrupo foi avaliado

separadamente o p-valor desses SNPs aumentaram.

Com relação aos SNPs com diferenças de frequência em reações hansênicas, foram

utilizados casos com reações hansênicas e casos que não as desenvolveram, utilizando

abordagem caso/controle com regressão logística. Foram encontrados trinta SNPs mais

fortemente associados, vinte dos quais apresentaram uma OR maior que 1, indicando risco

ao desenvolvimento de reações hansênicas, estando relacionados a MYADML/CRIM1,

LOC729653/UBD, CCHCR1, PSORS1C1, TREML2 e LRRC18/WDFY4. Um dos SNPs

listados se encontra próximo a uma região intergênica com o gene UBD (ubiquitina D).

Ubiquitinação é um importante processo bioquímico que controla diversos aspectos

da função proteica, e também está envolvido na indução de apoptose, controle do ciclo

celular e ativação de NF-kB, este último por sua vez controla inúmeros processos, como

imunidade, inflamação e apoptose. A ubiquitinação foi vista como um processo comum ao

espectro da hanseníase ao observar que o gene UBD é altamente expresso em hanseníase

per se (BELONE et al., 2015).

Outro SNP associado à reação hansênica se encontra próxima à região intergênica

LRRC18/WDFY4, região essa associada à susceptibilidade ao lúpus eritematoso sistêmico

(VISSCHER et al., 2012).

Diversos SNPs no cromossomo 6 associados a reações hansênicas foram detectados

próximos a genes conhecidos como de susceptibilidade a psoríase (PSORS1C1, PSORS1C2,

PSORS1C3, CCHCR1, POU5F1, HCG27 e HLA-C), estando em desequilíbrio de ligação

Page 58: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

56 (figura 16). LD indica a associação não-aleatória de alelos em dois ou mais loci, não

necessariamente no mesmo cromossomo (BALDING, 2006).

Diversos estudos genéticos foram realizados para entender a susceptibilidade

genética à psoríase, levando a identificação de pelo menos 19 loci de susceptibilidade,

envolvendo inflamação, apresentação de antígenos, sinalização celular e regulação

transcricional (LIANG et al., 2017). Um desses loci de susceptibilidade se encontra no

cromossomo 6, contendo, entre outros, os genes PSORS1C1, PSORS1C2, PSORS1C3,

CCHCR1, POU5F1, HCG27 e HLA-C (BOWCOCK; KRUEGER, 2005). Dentre eles,

CCHCR1 foi bem estudado quanto a seu papel na psoríase. CCHCR1 é superexpresso em

queratinócitos em lesões psoriáticas, e sua expressão é regulada por IFN-γ. Essa proteína

regula negativamente a diferenciação e proliferação de queratinócitos (ASUMALAHTI et

al., 2000).

Psoríase é uma doença inflamatória crônica mediada por células-T caracterizada por

placas escamosas espessadas na pele (BOWCOCK; KRUEGER, 2005). Essa doença

complexa é caracterizada por respostas autoimunes e autoinflamatórias, e é densamente

estudada quanto a susceptibilidade, de tal forma que SNPs associados à susceptibilidade à

psoríase foram adicionados no Immunochip (CORTES; BROWN, 2011; LIANG et al.,

2017).

Hanseníase e psoríase são doenças estigmatizadas desde tempos bíblicos, de tal

forma que psoríase já foi considerada uma forma de hanseníase (DOGRA; KAUR;

KUMAR, 2003). Ambas são doenças com perfis autoimune, desse modo, genes de

susceptibilidade seriam compartilhados, o que foi demonstrado em uma análise integrativa

comparando genes de susceptibilidade à hanseníase e a outras doenças autoimunes e

infecciosas. Foi visto que 19% destes genes são compartilhados com psoríase (ZHANG et

al., 2016). Alguns autores propõem a hipótese de que a que psoríase protegeria a progressão

clínica de hanseníase, devido a imunidade inata (BASSUKAS; GAITANIS; HUNDEIKER,

2012). Ao contrário, outras teorias sugerem que pacientes com hanseníase estariam

protegidos contra psoríase devido a neuropatia causada pelo M. leprae na pele (DOGRA;

KAUR; KUMAR, 2003). Dessa forma, hanseníase e psoríase raramente estariam presentes

no mesmo paciente, apesar de existirem relatos de casos nos quais foram observados a

coexistência das duas doenças (KUMAR et al., 1992; RAIOL et al., 2015).

Page 59: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

57

Interessante observar que, dentre os SNPs encontrados associados a genes de risco a

psoríase, todos apresentaram OR maior que 1, com exceção dos SNPs relacionados a

PSORS1C3/HCG27, HCG27 e PSORS1C1. Dessa forma os SNPs de risco à psoríase

também seriam de risco à hanseníase, em especial os SNPs próximos ao gene CCHCR1,

demonstrando mais uma vez que as duas doenças possuem perfis de resposta similares.

Também foram avaliados os SNPs com diferenças de frequência em variável contínua

de nível de anticorpos LID-NDO, utilizando regressão linear. Este tipo de análise é conduzido

assumindo uma relação linear entre as variáveis, ao invés de caso/controle. Além disso, é

utilizado o coeficiente beta como coeficiente de significância, sendo positivo ou negativo.

Quando o coeficiente beta é positivo, significa que para cada aumento de uma unidade na

variável preditora, a variável dependente aumentará pelo valor de beta, e o contrário quando o

coeficiente beta é negativo (STATISTICS SOLUTIONS, 2013). Dessa forma, todos os SNPs

associados a taxas de anticorpos são de proteção, com a exceção de quatro SNPs localizados

em LOC729170/CPEB4 e TMEM80.

Dentre os 15 SNPs mais fortemente associados nesta análise, é importante destacar o

SNP relacionado com o gene THEMIS, gene associado à susceptibilidade à esclerose múltipla,

que codifica uma proteína com papel regulador na seleção de células T (VISSCHER et al.,

2012). Cinco SNPs se encontram próximos ao gene RERE, o qual codifica uma proteína

contendo repetição de dipeptídeo de ácido glutâmico-arginina que controla a sinalização do

ácido retinóico, e sua expressão aumentada desencadeia a apopotose. Esse gene também foi

associado à susceptibilidade a psoríase (TSOI et al., 2012).

Observando todos os 55 SNPs associados no nosso estudo, podemos ver que somente

quatro se encontram dentro de regiões codantes (rs9263726 e rs1063646 em PSORS1C1, e

rs12364 e rs130077 em CCHCR1), todos associados a reações hansênicas. Um estudo

realizado por Liu et al., em 2017, demonstrou a importância de estudar variantes codantes em

doenças complexas, utilizando GWAS, com a descoberta de sete novas variantes codantes

(NCKIPSD, CARD9, IL23R, FLG, TYK2, SLC29A3 e IL27). Além disso, foi visto que os

efeitos genéticos individuais (ORs) de variantes de risco de codificação de baixa freqüência e

raras são mais fortes do que as variantes comuns, embora suas contribuições para o risco

sejam, individual ou cumulativamente, mais fracas (devido à sua baixa freqüência) do que as

comuns (LIU et al., 2017). Dessa forma, se torna importante um estudo mais detalhado desses

quatro SNPs.

Page 60: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

58

É interessante observar que quando analisamos os subgrupos de hanseníase (MB e PB)

e reações hansênicas (tipo I e tipo II), na primeira e segunda análise, percebemos que alguns

SNPs obtiveram valores de OR e p-valor diferentes daqueles relacionados ao grupo maior.

Isso demonstra a importância em analisar os subgrupos separadamente em hanseníase, já que

a doença apresenta um espectro patológico entre os pólos, além dos diferentes padrões de

resposta imune nas reações hansênicas.

O objetivo de estudos de associação populacionais é identificar padrões de

polimorfismos que variam sistematicamente entre indivíduos com diferentes estados de

doença e podem representar os efeitos de alelos de risco ou de proteção (BALDING, 2006).

Porém, foi visto que desde o advento de estudos GWAS em 2006, loci importantes do genoma

têm efeitos de tamanho pequenos e que, juntos, as associações significativas explicam apenas

uma pequena fração da variância genética prevista. Essa herdabilidade faltante de muitos

traços é representada por SNPs comuns com efeitos bem abaixo da significância estatística do

genoma. Dessa forma, Boyle et al, 2017, propôs que o risco de doença seria em grande parte

conduzido por genes sem relevância direta para a doença e é propagado através de redes

reguladoras para um número muito menor de genes principais com efeitos diretos. Isso

explicaria o porquê de termos visto tantos SNPs com p-valores abaixo do p-valor selecionado

para o tipo de estudo, e SNPs relacionados a genes que teoricamente não possuem uma

relevância direta para a hanseníase, como por exemplo HISPPD, que codifica um membro da

família de proteínas da fosfatase ácida da histidina, RIN3, cuja proteína funciona como uma

troca de nucleótidos de guanina para RAB5B e RAB31, e SEC11C.

O fato de um SNP estar dentro ou próximo a um gene não significa que este é um gene

causal para a doença, sendo necessária a validação dos SNPs em outra população. Dessa

forma, como perspectivas do nosso trabalho temos a escolha de SNPs mais associados e sua

validação em outras populações brasileiras.

Estudos funcionais são importantes para entender como produtos dos genes candidatos

estariam atuando na susceptibilidade à doença. Como no estudo de Manry et al, 2017, no qual

foi estudada a modulação da resposta transcricional à infecção e mecanismos de controle da

doença a partir de uma pesquisa genômica de loci de expressão quantitativa (eQTL) que estão

associados à variação da transcrição antes e após a estimulação com antígenos M. leprae em

células de sangue total (MANRY et al., 2017). Assim, também se torna necessário um estudo

funcional após a validação dos SNPs.

Page 61: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

59

Em síntese, nossos dados demonstram a importância do Immunochip como ferramenta

para a análise de polimorfismos associados a hanseníase, sendo o primeiro estudo utilizando

esta estratégia em hanseníase. Observamos no total 55 SNPs associados à hanseníase (10 em

hanseníase, 30 em reações hansênicas e 15 em taxas de produção de anticorpos), dentro de 36

regiões gênicas e intergênicas, sendo 4 SNPs em regiões codantes. Os dados genéticos e

sorológicos sugerem que o aumento da susceptibilidade para hanseníase e reações hansênicas

estão ligados à resposta imune celular e humoral do hospedeiro.

Page 62: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

60 6 CONCLUSÕES

1. Os níveis de produção de anticorpos são maiores em casos MB e em reações do tipo

II, quando comparados com casos PB e as outras reações hansênicas,

respectivamente;

2. Foram observados 10 SNPs associados a susceptibilidade à hanseníase, 30 SNPs

associados a susceptibilidade à reação hansênica e 15 SNPs associados a taxas de

produção de anticorpos;

3. A hanseníase possui genes de susceptibilidade comum a outras doenças autoimunes

(lúpus eritematoso sistêmico, Esclerose múltipla e psoríase).

Page 63: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

61 REFERÊNCIAS

ALI, S. et al. IL12B SNPs and copy number variation in IL23R gene associated with susceptibility to leprosy. Journal of Medical Genetics, v. 50, n. 1, p. 34–42, 13 jan. 2013.

ALTER, A. et al. Leprosy as a genetic model for susceptibility to common infectious diseases. Human Genetics, v. 123, n. 3, p. 227–235, 2008.

ALTER, A. et al. Genetic and functional analysis of common MRC1 exon 7 polymorphisms in leprosy susceptibility. Human Genetics, v. 127, n. 3, p. 337–348, 25 mar. 2010.

AMORIM, F. M. et al. Identifying Leprosy and Those at Risk of Developing Leprosy by Detection of Antibodies against LID-1 and LID-NDO. PLoS Neglected Tropical Diseases, v. 10, n. 9, p. 1–17, 2016.

ANDERSON, C. A et al. Europe PMC Funders Group Data quality control in genetic case-control association studies. Nature protocols, v. 5, n. 9, p. 1564–1573, 2011.

ARANOW, C. Vitamin D and the immune system. Journal of investigative medicine : the official publication of the American Federation for Clinical Research, v. 59, n. 6, p. 881–6, 1 ago. 2011.

ARAÚJO, T. G. et al. A novel vitamin D receptor polymorphism associated with leprosy. Journal of dermatological science, v. 89, n. 3, p. 304–307, 2017.

ASUMALAHTI, K. et al. A candidate gene for psoriasis near HLA-C, HCR (Pg8), is highly polymorphic with a disease-associated susceptibility allele. Human Molecular Genetics, v. 9, n. 10, p. 1533–1542, 2000.

BALDING, D. J. A tutorial on statistical methods for population association studies. Nat Rev Genet, v. 7, n. 10, p. 781–791, 2006.

BARREIRO, L. B.; QUINTANA-MURCI, L. From evolutionary genetics to human immunology: how selection shapes host defence genes. Nature reviews. Genetics, v. 11, n. 1, p. 17–30, jan. 2010.

BASSUKAS, I. D.; GAITANIS, G.; HUNDEIKER, M. Leprosy and the natural selection for psoriasis. Medical Hypotheses, v. 78, n. 1, p. 183–190, jan. 2012.

BELONE, A. DE F. F. et al. Genome-wide screening of mRNA expression in leprosy patients. Frontiers in Genetics, v. 6, n. NOV, p. 1–12, 2015.

BLEHARSKI, J. R. et al. Use of Genetic Profiling in Leprosy to Discriminate Clinical Forms of the Disease. Science, v. 301, n. 5639, p. 1527–1530, 12 set. 2003.

BOCHUD, P. et al. Toll‐Like Receptor 2 (TLR2) Polymorphisms Are Associated with Reversal Reaction in Leprosy. The Journal of Infectious Diseases, v. 197, n. 2, p. 253–261, 2008.

BOCHUD, P.-Y. et al. Polymorphisms in Toll-like receptor 4 (TLR4) are associated with protection against leprosy. European Journal of Clinical Microbiology & Infectious Diseases, v. 28, n. 9, p. 1055–1065, 9 set. 2009.

Page 64: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

62 BOWCOCK, A. M.; KRUEGER, J. G. Getting under the skin: the immunogenetics of psoriasis. Nature reviews. Immunology, v. 5, n. 9, p. 699–711, 2005.

BOYLE, E. A.; LI, Y. I.; PRITCHARD, J. K. An Expanded View of Complex Traits: From Polygenic to Omnigenic. Cell, v. 169, n. 7, p. 1177–1186, 2017.

BRITTON, W. J.; LOCKWOOD, D. N. J. Leprosy. Lancet, v. 363, n. 9416, p. 1209–1219, 2004.

BÜHRER-SÉKULA, S. et al. Simple and fast lateral flow test for classification of leprosy patients and identification of contacts with high risk of developing leprosy. Journal of clinical microbiology, v. 41, n. 5, p. 1991–5, maio 2003.

BÜHRER-SÉKULA, S. et al. The ML flow test as a point of care test for leprosy control programmes: potential effects on classification of leprosy patients. Leprosy review, v. 78, n. 1, p. 70–9, mar. 2007.

CARDOSO, C. C. et al. Ninjurin 1 asp110ala single nucleotide polymorphism is associated with protection in leprosy nerve damage. Journal of neuroimmunology, v. 190, n. 1–2, p. 131–8, 1 out. 2007.

CARDOSO, C. C. et al. Leprosy susceptibility: genetic variations regulate innate and adaptive immunity, and disease outcome. Future Microbiology, v. 6, n. 5, p. 533–549, maio 2011.

CHAITANYA, V. S. et al. Interleukin-17F single-nucleotide polymorphism (7488T&gt;C) and its association with susceptibility to leprosy. International Journal of Immunogenetics, v. 41, n. 2, p. 131–137, abr. 2014.

COLE, S. T. et al. Massive gene decay in the leprosy bacillus. Nature, v. 409, n. 6823, p. 1007–11, 2001.

CORTES, A.; BROWN, M. A. Promise and pitfalls of the Immunochip. Arthritis Res Ther, v. 13, n. 1, p. 101, 2011.

DE MESSIAS-REASON, I. J. et al. The association between mannan-binding lectin gene polymorphism and clinical leprosy: new insight into an old paradigm. The Journal of infectious diseases, v. 196, n. 9, p. 1379–1385, 2007.

DOGRA, S.; KAUR, I.; KUMAR, B. Leprosy and Psoriasis: An Enigmatic Relationship. International Journal of Leprosy and Other Mycobacterial Diseases, v. 71, n. 4, p. 341, dez. 2003.

DOS SANTOS, D. S. et al. Kinship and leprosy in the contacts of leprosy patients: Cohort at the Souza Araújo outpatient clinic, Rio de Janeiro, RJ, 1987-2010. Journal of Tropical Medicine, v. 2013, 2013.

DUPNIK, K. M. et al. Transcriptional changes that characterize the immune reactions of leprosy. Journal of Infectious Diseases, v. 211, n. 10, p. 1658–1676, 2015.

DUTHIE, M. S. et al. A rapid ELISA for the diagnosis of MB leprosy based on complementary detection of antibodies against a novel protein-glycolipid conjugate. Diagnostic Microbiology and Infectious Disease, v. 79, n. 2, p. 233–239, 2014.

FABRI, A. DA C. O. C. et al. Antigen-specific assessment of the immunological status of

Page 65: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

63 various groups in a leprosy endemic region. BMC infectious diseases, v. 15, p. 218, 2015.

FAVA, V. M. et al. Association of TNFSF8 Regulatory Variants With Excessive Inflammatory Responses but not Leprosy Per Se. The Journal of Infectious Diseases, v. 211, n. 6, p. 968–977, 15 mar. 2015.

FAVA, V. M. et al. A genome wide association study identifies a lncRna as risk factor for pathological inflammatory responses in leprosy. PLOS Genetics, v. 13, n. 2, p. e1006637, 2017.

FITNESS, J. et al. Large-scale candidate gene study of leprosy susceptibility in the Karonga district of northern Malawi. The American journal of tropical medicine and hygiene, v. 71, n. 3, p. 330–40, set. 2004.

FRANCESCHI, D. S. A. et al. Association between killer-cell immunoglobulin-like receptor genotypes and leprosy in Brazil. Tissue Antigens, v. 72, n. 5, p. 478–482, nov. 2008.

GASCHIGNARD, J. et al. Pauci- and Multibacillary Leprosy: Two Distinct, Genetically Neglected Diseases. PLoS Neglected Tropical Diseases, v. 10, n. 5, p. 1–20, 2016.

GOULART, I. M. B.; PENNA, G. O.; CUNHA, G. Immunopathology of leprosy: The complexity of the mechanisms of host immune response to Mycobacterium leprae. Revista da Sociedade Brasileira de Medicina Tropical, v. 35, n. 4, p. 365–375, 2002.

GRANT, A. V. et al. Crohn’s Disease Susceptibility Genes are Associated With Leprosy in the Vietnamese Population. Journal of Infectious Diseases, v. 206, n. 11, p. 1763–1767, 1 dez. 2012.

KAHAWITA, I. P.; LOCKWOOD, D. N. J. Towards understanding the pathology of erythema nodosum leprosum. Transactions of the Royal Society of Tropical Medicine and Hygiene, v. 102, n. 4, p. 329–337, 2008.

KUMAR, B. et al. THE RARE COEXISTENCE OF LEPROSY AND PSORIASIS. International Journal of Dermatology, v. 31, n. 8, p. 551–554, 1 ago. 1992.

LIANG, Y. et al. Psoriasis: a mixed autoimmune and autoinflammatory disease. Current Opinion in Immunology, v. 49, p. 1–8, 2017.

LIU, H. et al. Discovery of six new susceptibility loci and analysis of pleiotropic effects in leprosy. Nature Genetics, v. 47, n. 3, p. 267–271, 2015.

LIU, H. et al. Genome-wide Analysis of Protein-Coding Variants in Leprosy. The Journal of investigative dermatology, 22 ago. 2017.

LU’O’NG, K. VINH QUỐC; HOÀNG NGUYỄN, L. T. Role of the Vitamin D in Leprosy. The American Journal of the Medical Sciences, v. 343, n. 6, p. 471–482, jun. 2012.

LUCENA-SILVA, N. et al. The +3187A/G HLA-G polymorphic site is associated with polar forms and reactive reaction in leprosy. Molecular Genetics & Genomic Medicine, v. 1, n. 3, p. 123–130, 1 set. 2013.

MANDAL, D. et al. Vitamin D receptor expression levels determine the severity and complexity of disease progression among leprosy reaction patients. New microbes and new infections, v. 6, p. 35–9, 2015.

Page 66: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

64 MANRY, J. et al. Deciphering the genetic control of gene expression following Mycobacterium leprae antigen stimulation. PLoS genetics, v. 13, n. 8, p. e1006952, ago. 2017.

MARCINEK, P. et al. LRRK2 and RIPK2 Variants in the NOD 2-Mediated Signaling Pathway Are Associated with Susceptibility to Mycobacterium leprae in Indian Populations. PLoS ONE, v. 8, n. 8, p. e73103, 28 ago. 2013.

MARQUET, S.; SCHURR, E. Genetics of susceptibility to infectious diseases: Tuberculosis and leprosy as examples. Drug Metabolism and Disposition, v. 29, n. 4 II, p. 479–483, 2001.

MEISNER, S. J. et al. Association of NRAMP1 polymorphism with leprosy type but not susceptibility to leprosy per se in west Africans. The American journal of tropical medicine and hygiene, v. 65, n. 6, p. 733–5, dez. 2001.

MILLER, E. N. et al. Genome-wide scans for leprosy and tuberculosis susceptibility genes in Brazilians. Genes and immunity, v. 5, n. 1, p. 63–67, 2004.

MIRA, M. T. et al. Segregation of HLA/TNF region is linked to leprosy clinical spectrum in families displaying mixed leprosy subtypes. Genes and immunity, v. 4, p. 67–73, 2003a.

MIRA, M. T. et al. Chromosome 6q25 is linked to susceptibility to leprosy in a Vietnamese population. Nature genetics, v. 33, n. 3, p. 412–415, 2003b.

MIRA, M. T. et al. Susceptibility to leprosy is associated with PARK2 and PACRG. Nature, v. 427, n. 6975, p. 636–640, 2004.

MISCH, E. A et al. Leprosy and the human genome. Microbiology and molecular biology reviews : MMBR, v. 74, n. 4, p. 589–620, 2010.

MIZOGUTI, D. DE F. et al. Multibacillary leprosy patients with high and persistent serum antibodies to leprosy IDRI diagnostic-1/LID-1: higher susceptibility to develop type 2 reactions. Memórias do Instituto Oswaldo Cruz, v. 110, n. 7, p. 914–920, 2015.

MOET, F. J. et al. Physical distance, genetic relationship, age, and leprosy classification are independent risk factors for leprosy in contacts of patients with leprosy. The Journal of infectious diseases, v. 193, n. 3, p. 346–353, 2006.

MONOT, M. et al. On the Origin of Leprosy. Science, v. 308, n. 5724, p. 1040–1042, 13 maio 2005.

MONOT, M. et al. Comparative genomic and phylogeographic analysis of Mycobacterium leprae. Nature genetics, v. 41, n. 12, p. 1282–1289, 2009.

MORAES, M. O. et al. Interleukin-10 promoter single-nucleotide polymorphisms as markers for disease susceptibility and disease severity in leprosy. Genes and Immunity, v. 5, n. 7, p. 592–595, 12 nov. 2004.

MOURA, M. L. N. et al. Active surveillance of Hansen’s Disease (leprosy): importance for case finding among extra-domiciliary contacts. PLoS neglected tropical diseases, v. 7, n. 3, p. e2093, jan. 2013.

MUSTAFA, T. et al. Increased expression of fas ligand in human tuberculosis and leprosy

Page 67: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

65 lesions: a potential novel mechanism of immune evasion in mycobacterial infection . Scand.J.Immunol. , v. 54, n. 0300–9475, p. 630–639, 2001.

NEELA, V. S. K. et al. Association of Taq I, Fok I and Apa I polymorphisms in Vitamin D Receptor (VDR) gene with leprosy. Human Immunology, v. 76, n. 6, p. 402–405, jun. 2015.

OHYAMA, H. et al. Differential Effects of Polymorphisms in the 5′ Flanking Region of IL12RB2 on NK- and T-Cell Activity. Journal of Interferon & Cytokine Research, v. 28, n. 9, p. 563–570, set. 2008.

OLIVEIRA, A. L. G. DE et al. Vitamin D receptor expression and hepcidin levels in the protection or severity of leprosy: a systematic review. Microbes and Infection, v. 19, n. 6, p. 311–322, 2017.

PRADO-MONTES DE OCA, E. et al. SNP 668C (−44) alters a NF-κB1 putative binding site in non-coding strand of human β-defensin 1 (DEFB1) and is associated with lepromatous leprosy. Infection, Genetics and Evolution, v. 9, n. 4, p. 617–625, jul. 2009.

PRASAD, S. et al. Leprosy revealed in a rheumatology clinic: A case series. International Journal of Rheumatic Diseases, v. 16, n. 2, p. 129–133, abr. 2013.

PURCELL, S. et al. PLINK: A tool set for whole-genome association and population-based linkage analyses. American Journal of Human Genetics, v. 81, n. 3, p. 559–575, 2007.

QUEIROZ, J. W. et al. Geographic information systems and applied spatial statistics are efficient tools to study Hansen’s disease (leprosy) and to determine areas of greater risk of disease. American Journal of Tropical Medicine and Hygiene, v. 82, n. 2, p. 306–314, 2010.

RAIOL, T. K. A. et al. Leprosy associated with psoriasis. Lepr Rev, v. 86, p. 368–373, 2015.

RAJALINGAM, R.; MEHRA, N. K.; SINGAL, D. P. Polymorphism in heat-shock protein 70-1 (HSP70-1) gene promoter region and susceptibility to tuberculoid leprosy and pulmonary tuberculosis in Asian Indians. Indian journal of experimental biology, v. 38, n. 7, p. 658–62, jul. 2000.

RAMBUKKANA, A. How does Mycobacterium leprae target the peripheral nervous system? Trends in Microbiology, v. 8, n. 1, p. 23–28, 2000.

REIBEL, F.; CAMBAU, E.; AUBRY, A. Update on the epidemiology, diagnosis, and treatment of leprosy. Medecine et Maladies Infectieuses, 2015.

REYNARD, M. P. et al. Allele frequencies for an interferon-gamma microsatellite in a population of Brazilian leprosy patients. European Journal of Immunogenetics, v. 30, n. 2, p. 149–151, 1 abr. 2003.

RIDLEY, D. S.; JOPLING, W. H. A classification of leprosy for research purposes. Leprosy review, v. 33, p. 119–28, abr. 1962.

SALES-MARQUES, C. et al. Genetic polymorphisms of the IL6 and NOD2 genes are risk factors for inflammatory reactions in leprosy. PLoS neglected tropical diseases, v. 11, n. 7, p. e0005754, jul. 2017.

SAMPAIO, P. B. et al. Spatial analysis of new cases of leprosy in the State of Espírito Santo,

Page 68: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

66 Brazil, between 2004 and 2009. Revista da Sociedade Brasileira de Medicina Tropical, v. 45, n. 3, p. 380–384, jun. 2012.

SANTOS, A. R. et al. Role of Tumor Necrosis Factor–α and Interleukin‐10 Promoter Gene Polymorphisms in Leprosy. The Journal of Infectious Diseases, v. 186, n. 11, p. 1687–1691, 1 dez. 2002.

SAPKOTA, B. R. et al. Association of TNF, MBL, and VDR polymorphisms with leprosy phenotypes. Human Immunology, v. 71, n. 10, p. 992–998, out. 2010.

SCHURING, R. P. et al. Polymorphism N248S in the Human Toll‐Like Receptor 1 Gene Is Related to Leprosy and Leprosy Reactions. The Journal of Infectious Diseases, v. 199, n. 12, p. 1816–1819, 15 jun. 2009.

SCOLLARD, D. M. et al. The continuing challenges of leprosy. Clinical Microbiology Reviews, v. 19, n. 2, p. 338–381, 2006.

SERGIO DO SACRAMENTO, W. et al. Frequencies of MICA alleles in patients from southern Brazil with multibacillary and paucibacillary leprosy. International Journal of Immunogenetics, v. 39, n. 3, p. 210–215, 1 jun. 2012.

SINGH, P.; COLE, S. T. Mycobacterium leprae : genes, pseudogenes and genetic diversity. Future Microbiology, v. 6, n. 1, p. 57–71, 16 jan. 2011.

SOUSA, A. L. M. et al. Genetic and Immunological Evidence Implicates Interleukin 6 as a Susceptibility Gene for Leprosy Type 2 Reaction. The Journal of Infectious Diseases, v. 205, n. 9, p. 1417–1424, 1 maio 2012.

SPENCER, J. S. et al. Identification of serological biomarkers of infection, disease progression and treatment efficacy for leprosy. Memórias do Instituto Oswaldo Cruz, v. 107, n. suppl 1, p. 79–89, dez. 2012.

SPENCER, J. S.; BRENNAN, P. J. The role of Mycobacterium leprae phenolic glycolipid I (PGL-I) in serodiagnosis and in the pathogenesis of leprosy. Leprosy review, v. 82, n. 4, p. 344–57, dez. 2011.

SUZUKI, K. et al. Current status of leprosy: epidemiology, basic science and clinical perspectives. The Journal of dermatology, v. 39, n. 2, p. 121–9, 2012.

TOBIN, D. M. et al. The lta4h locus modulates susceptibility to mycobacterial infection in zebrafish and humans. Cell, v. 140, n. 5, p. 717–30, 5 mar. 2010.

TRYNKA, G. et al. Dense genotyping identifies and localizes multiple common and rare variant association signals in celiac disease. Nature genetics, v. 43, n. 12, p. 1193–1201, 2011.

TSOI, L. C. et al. Identification of 15 new psoriasis susceptibility loci highlights the role of innate immunity. Nature genetics, v. 44, n. 12, p. 1341–8, dez. 2012.

VEJBAESYA, S. et al. NRAMP1 and TNF-? polymorphisms and susceptibility to tuberculosis in Thais. Respirology, v. 12, n. 2, p. 202–206, 1 mar. 2007.

VISSCHER, P. M. et al. Five years of GWAS discovery. American journal of human genetics, v. 90, n. 1, p. 7–24, 13 jan. 2012.

Page 69: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

67 WALKER, S. L.; LOCKWOOD, D. N. J. Leprosy type 1 (reversal) reactions and their management. Leprosy review, v. 79, n. 4, p. 372–386, 2008.

WHO. WHOTreatmentGuide.pdf. p. 28, 1998.

WHO. Global leprosy update 2013. Weekly epidemiological record, v. 89, n. 36, p. 389–400, 2014.

WHO. Leprosy: weekly epidemiological record, Septembre 2016. World Health Organisation Weekly epidemiological record, v. 91, n. 35, p. 405–420, 2016.

WIBAWA, T.; SOEBONO, H.; MATSUO, M. Association of a missense mutation of the laminin alpha2 gene with tuberculoid type of leprosy in Indonesian patients. Tropical medicine & international health : TM & IH, v. 7, n. 7, p. 631–6, jul. 2002.

XIANG, Y.-L. et al. Common variants of OPA1 conferring genetic susceptibility to leprosy in Han Chinese from Southwest China. Journal of Dermatological Science, v. 80, n. 2, p. 133–141, nov. 2015.

ZHANG, D. F. et al. Integrative analyses of leprosy susceptibility genes indicate a common autoimmune profile. Journal of Dermatological Science, v. 82, n. 1, p. 18–27, 2016.

ZHANG, F.-R. R. et al. Genomewide association study of leprosy. New England Journal of Medicine, v. 361, n. 27, p. 2609–2618, 2009.

Page 70: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

68 APÊNDICE 1: TERMO DE CONSENTIMENTO LIVRE E ESCLARECIDO (TCLE)

Page 71: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

69

Universidade Federal do Rio Grande do Norte Centro de Biociências

Universidade Federal da Bahia Serviço de Imunologia

TERMO DE CONSENTIMENTO LIVRE E ESCLARECIDO (2)

Título do Projeto: INSTITUTO NACIONAL DE ESTUDOS EM DOENÇAS TROPICAIS: MARCADORES GENÉTICOS E IMUNOLÓGICOS NA HANSENÍASE

Pesquisadores: Selma M.B. Jerônimo, Paulo Roberto Machado, Edgar M. Carvalho, Pedro Bezerra da Trindade Neto, Shirley Moreira, Márcia C.S. Dias, Goretti Fernandes de Oliveira, Lúcio Leopoldino de Medeiros, Kaliane Fernandes Vale, Sérgio R. F. Araújo, Kathryn Dupnik, Mauricio L. Nobre

.I. Esclarecimentos

Estamos convidando você para participar, como voluntário, de uma pesquisa sobre hanseníase. Com essa pesquisa nós queremos aprender porque somente algumas pessoas quando contaminadas com o germe que causa esta doença a desenvolve.. Para isto, vamos realizar um estudo que usa DNA e outros componentes do sangue. O DNA é uma substância encontrada dentro das nossas células, herdada dos nossos pais e transmitidas por nós aos nossos filhos. Iremos estudar “pedacinhos” do DNA que na linguagem científica se chama de genes, que podem estar relacionados com a capacidade (sensibilidade) ou não (resistência) de desenvolver doenças.

Esta pesquisa será feita em colaboração com a Universidade Federal da Bahia. Nesta pesquisa nós iremos realizar um estudo comparativo entre famílias oriundas da Bahia e do Rio Grande do Norte, com o objetivo de identificar formas de melhor prevenir esta doença. Caso você aceite participar do estudo você será submetido a uma consulta médica e deverá permitir a realização de coleta de sangue. Solicitamos permissão para revisar o seu prontuário que está no posto de saúde ou hospital onde você está ou foi tratado para hanseníase. Para alguns participantes que têm formas clínicas da hanseníase com maior risco de desenvolver complicações, denominadas cientificamente reações reversa ou eritema nodoso, solicitamos permissão para contatar mensalmente, por telefone, para saber sobre a resposta ao tratamento e a presença de alguma dessas complicações.

Você também deverá responder um questionário a respeito da sua saúde. Quando você estiver respondendo esse questionário se alguma pergunta lhe causar constrangimento você poderá deixar de respondê-la.

Vamos, a seguir, explicar a você o motivo da coleta do material e os riscos e desconfortos dos procedimentos. Coleta de sangue serão coletados 10ml de sangue (mais ou menos uma colher de sopa) para a realização dos estudos genéticos e um adicional de 15 ml para aquelas pessoas que participarão dos estudos referentes a resposta de defesa. O sangue será coletado com material estéril e descartável. Uma segunda amostra de sangue (aproximadamente 15 ml) será coletada para aqueles que tiverem complicações como reações reversa ou eritema nodoso. Os riscos e desconfortos deste procedimento são pequenos e poderá ser desmaio, sangramento, manchas arroxeadas ou infecção no local da coleta. Iremos minimizar estes riscos com os cuidados: limpeza no local da coleta usando-se álcool e pressionando este local, por alguns minutos, após a retirada da agulha. Com a parte líquida do sangue (soro ou plasma), iremos determinar se voce foi contaminado com o germe que causa hanseníase; exames laboratoriais, como hemograma, avaliação da função do fígado e exposição a germes que causam hepatites e a e Síndrome da Imunodeficiência Adquirida (SIDA) serão realizados. O exame para detectar o vírus HIV só será feito se você autorizar. Das células, que é a parte figurada do sangue, iremos isolar o DNA e RNA e estudar seu sistema de defesa, a sua imunidade. Se for preciso, você poderá ser submetido a uma biópsia de pele. Este exame pode ser doloroso. Esta coleta será feita por um médico especializado e ele deverá usar um anestésico local. Parte do material coletado será usado para avaliação do estágio atual da sua doença. Iremos também registrar o local da sua residência. Para isto iremos utilizar um aparelho chamado GPS. Este aparelho captura sinais de satélites presentes na atmosfera e permite com precisão localizar um ponto na superfície da Terra. O objetivo deste procedimento é identificar na sua cidade áreas onde ocorrem com mais freqüência casos de doença. Isto permitirá direcionar medidas de intervenção que podem facilitar o controle da hanseníase. Sua participação é importante para a nossa pesquisa. Ao responder os questionários você estará contribuindo com informações que nos permitirão mapear áreas onde existe hanseníase no Rio Grande do Norte. Enfim,

Versão 3_CAAE - 0042.0.051.051-09 – INCT_ DT: MARCADORES GENÉTICOS E IMUNOLÓGICOS NA HANSENÍASE_CASO

1 DE 2

Page 72: UNIVERSIDADE FEDERAL DO RIO GRANDE DO NORTE PRÓ … · presença de anticorpos anti LID-NDO e apresentação clínica. A quantidade de anticorpos variou de acordo com a classificação

70

sua participação nesta pesquisa contribuirá para que uma melhor compreensão sobre a hanseníase. Mas, se você decidir não participar nesta pesquisa o tratamento para hanseníase ou outra doença que for diagnosticada em você não será diferente, além do que, você poderá retirar seu consentimento a qualquer momento da pesquisa, sem que essa desistência traga qualquer prejuízo a você ou sua família.

Garantimos que os resultados dos seus exames serão mantidos em sigilo e privacidade onde somente pesquisadores deste estudo, terão acesso. Cada indivíduo e todas as amostras coletadas serão identificados por código de barra. Não haverá disponibilização de seus resultados para terceiros (seguradora, empregadores, etc). Se qualquer relatório ou publicação resultar deste estudo, seu nome não será revelado. Porém, caso você deseje poderá receber os resultados de todos os seus exames. Assumimos o compromisso de que o material coletado será usado apenas para o que foi esclarecido aqui. Pedimos sua permissão para guardar (armazenar) o restante do DNA, do soro e do plasma coletados nessa pesquisa, no Laboratório de Imunogenética, na Universidade Federal do Rio Grande do Norte. Este material será guardado em tubos identificados por um código durante um período de 5 anos. Qualquer outro estudo que venhamos a desenvolver e que usará informação sobre sua doença ou material do seu sangue será primeiramente submetido ao Comitê de Ética em Pesquisa da UFRN, à Comissão Nacional de Ética em Pesquisa (CONEP), e solicitada sua autorização. Lembramos que você poderá remover sua permissão para o armazenamento do material em qualquer época, sem qualquer prejuízo para você. Sua participação nesta pesquisa é voluntária e por quatro anos. Você não será pago por participar desta pesquisa, mas será ressarcido (reembolsado) por qualquer gasto decorrente da mesma e será indenizado por qualquer dano também comprovadamente decorrente de sua participação nesta pesquisa. Se você tem, ou tiver, alguma dúvida sobre os esclarecimentos feitos aqui, riscos, benefícios, seus direitos ou qualquer outra dúvida a respeito desta pesquisa entrar em contato com a Dra Selma Jeronimo, ou outro membro da equipe, pelo telefone 84-3215-3428, no Departamento de Bioquímica da UFRN ou por correio eletrônico ([email protected]). Este projeto foi avaliado pelo Comitê de Ética em Pesquisa da Universidade Federal do Rio Grande do Norte e pela Comissão Nacional de Ética em Pesquisa (CONEP/CNS/Ministério da Saúde, Brasília). O Comitê de Ética da Universidade Federal do Rio Grande do Norte fica localizado no Campus Universitário (84-3215- 3135). O número do protocolo de pesquisa é CAAE 0042.0.051.051-09. 2. Consentimento Livre

Concordo em participar voluntariamente do estudo “INCT-DOENÇAS TROPICAIS: MARCADORES GENÉTICOS E IMUNOLÓGICOS NA HANSENÍASE”. Fui devidamente esclarecido (a) quanto ao material que será coletado, os testes que farei e o questionário que responderei. Entendi todos os riscos e desconfortos que poderei ter, assim como, os benefícios para mim e para a comunidade que este estudo poderá trazer.

Autorizo ou não os pesquisadores o seguinte: Sim Não • pesquisar o vírus HIV no meu sangue • usar parte do material da biópsia de pele • armazenar o restante do meu DNA, soro e plasma no Laboratório de

Imunogenética da Universidade Federal do Rio Grande do Norte.

Por fim, quanto aos resultados dos meus exames manifesto o desejo de:

• Receber todos os resultados Não receber os resultado DATA / / Assinatura do Participante DATA / / Assinatura do Pesquisador Responsável

Versão 3_CAAE - 0042.0.051.051-09 – INCT_ DT: MARCADORES GENÉTICOS E IMUNOLÓGICOS NA HANSENÍASE_CASO

2 D