O reflexo de luz pupilar (PLR) é um fenótipo intermediário candidato
associado ao autismo. Neste estudo, objetivou-se expandir pesquisas
anteriores indicando a PLR como um fenótipo intermediário candidato para
o autismo, explorando a variabilidade epigenética associada a diferenças
individuais no desenvolvimento inicial da PLR.
Leia
o artigo
Os pesquisadores mediram um reflexo do cérebro que aparece muito
cedo, que se chama PLR: “Reflexo Pupilar à Luz” e em vez de
procurar mutações, o estudo investiga como os genes estão sendo
regulados desde o início da vida.
O que faz esse reflexo?
Quando a luz aumenta, a pupila fecha. Quanto mais rápido e
quanto mais ela fecha, mais eficiente está um circuito específico do
cérebro.
Esse circuito envolve:
- Retina;
- Nervo óptico;
- Tronco cerebral;
- Sistema Nervoso Autônomo (Simpático + Parassimpático)
Esse é um dos circuitos mais básicos do cérebro humano e por isso é perfeito para estudar desenvolvimento neurológico precoce.
Porque bebês autistas já mostram diferenças nesse reflexo, meses antes de qualquer comportamento aparecer.
Os pesquisadores já sabiam, de estudos anteriores, que:
| Trajetória da resposta pupilar em crianças que posteriormente são diagnosticadas com autismo | |
| Alterações no reflexo pupilar à luz ao longo do desenvolvimento | |
| Fase do desenvolvimento | O que acontece na pupila |
|---|---|
| Aos 9–10 meses | A pupila se fecha de forma mais intensa |
| De 9 → 14 meses | A resposta pupilar começa a enfraquecer |
| De 14 → 24 meses | A resposta pupilar torna-se mais lenta |
| Resumo baseado em estudos longitudinais sobre o reflexo pupilar à luz (PLR). | |
E isso significa que o circuito que controla a pupila está se
desenvolvendo de forma diferente.
A pergunta central para essa pesquisa, foi: “As diferenças
precoces no reflexo pupilar (PLR) podem ser explicadas por alterações
epigenéticas - especificamente metilação do DNA - e não por mutações na
sequência do DNA?”
(Se você não sabe o que é metilação do DNA, pense nela como uma
“etiqueta química” que pode aumentar ou reduzir a atividade de certos
genes, ajudando a regular quando e quanto um gene é expresso, sem
alterar a sequência do DNA.)
Para isso, eles coletaram:
- DNA da boca dos bebês aos 9 meses;
- Mediram a pupila aos 9, 14 e 24 meses.
A pergunta durante as coletas, foi:
“As marcas no DNA aos 9 meses preveem como a pupila vai se comportar
depois?”
Para velocidade da pupila (latência), encontraram
algo muito forte:
- 4 regiões do DNA extremamente significativas;
- 13 regiões maiores (DMRs) também significativas.
Isso foi especialmente forte entre 14 e 24 meses (que é justamente
quando as características autísticas começam a emergir).
Essas regiões regulam genes ligados a:
- Formação de neurônios;
- Crescimento de circuitos;
- Organização sináptica.
| Tabela 1. Resumo das análises EWAS por sonda e dos DMRs identificados para cada fenótipo | |||||
| Análise epigenômica ampla do reflexo pupilar à luz em lactentes | |||||
| N | λ |
Análise de DMPs associados ao PLR
|
Análise de DMRs associados ao PLR
|
||
|---|---|---|---|---|---|
| Sondas com significância estrita | Sondas significativas na fase de descoberta | DMRs significativos após ajuste de Bonferroni | |||
| Latência (ms) | |||||
| 9 meses | 48 | 1.16 | 0 | 28 | 1 |
| 14 meses | 47 | 1.14 | 2 | 27 | 1 |
| 24 meses | 40 | 1.02 | 1 | 21 | 1 |
| 9–14 meses | 44 | 0.94 | 0 | 6 | 1 |
| 14–24 meses | 37 | 1.20 | 1 | 60 | 9 |
| Amplitude (%) | |||||
| 9 meses | 46 | 1.07 | 0 | 11 | 1 |
| 14 meses | 47 | 1.03 | 0 | 21 | 1 |
| 24 meses | 37 | 1.05 | 0 | 28 | 5 |
| 9–14 meses | 42 | 1.06 | 0 | 28 | 3 |
| 14–24 meses | 34 | 0.96 | 0 | 15 | 7 |
| λ = fator de inflação genômica; sig. = significativo; Bonferroni = ajuste de Bonferroni. | |||||
Dois genes chaves apareceram, que são de alto risco para autismo no
banco SFARI:
- NR4A2: Responsável pelos neurônios dopaminérgicos (motivação,
foco, resposta sensorial);
- HNRNPU: Responsável pela formação de sinapses e redes
neurais.
Dessa forma, eles concluiram que o reflexo da pupila está
conectado a genes que constroem o cérebro autista.
Para quanto a pupila fecha (amplitude), o efeito foi
mais difuso, pois não acharam um único gene dominante, mas uma rede
enorme envolvendo:
Isso significa que a força da resposta da pupila é regulada por
muitos sistemas ao mesmo tempo, não só pelo circuito reflexo e por isso
os sinais epigenéticos são mais espalhados.
(Embora NR4A e HNRNPU não tenham sido identificados como alvos
epigenéticos diretos neste estudo, ambos integram os processos
biológicos enriquecidos identificados pelas análises funcionais,
particularmente aqueles relacionados ao neurodesenvolvimento, regulação
gênica e organização nuclear.
Eles aparecem no estudo de forma indireta, como parte das redes
biológicas representadas pelos termos de Ontologia Gênica.)
Esse estudo nos mostra algo muito profundo: que o autismo
começa em circuitos sensoriais básicos, não em
comportamento.
Ou seja, que antes da:
- Linguagem;
- Socialização e
- Cognição.
Já existe um cérebro que responde diferente à luz, regulado por
genes do neurodesenvolvimento que são modulados pela metilação do
DNA.
Isso significa que o autismo começa com uma diferença na
forma como o cérebro calibra o mundo sensorial.
Resumindo: essas diferenças iniciais na calibração sensorial podem
ser observadas precocemente por meio do reflexo pupilar à luz,
possivelmente associadas a alterações epigenéticas, como a metilação do
DNA.
É revolucionário pois eles mediram isso aos 9 meses, ou seja:
- Muito antes de qualquer diagnóstico;
- Usando um reflexo automático;
- Ligado a marcas químicas no DNA.
Isso significa que no futuro, pode-se:
- Identificar risco;
- Entender o tipo de cérebro;
- Oferecer suporte sensorial precoce sem rotular e/ou patologizar.
E isso não é “defeito” e sim um estilo de
neurodesenvolvimento mensurável no nível molecular.
A tabela 2 é central, porque ela mostra o nível mais fino do sinal
epigenético (CpG por CpG) antes de “subir” para DMRs e biologia.
Ela lista sondas CpG individuais cuja metilação está significativamente
associada à latência do reflexo pupilar à luz (PLR) em diferentes idades
da infância. Ou seja, em quais pontos específicos do genoma a
metilação do DNA se relaciona com o tempo que a pupila leva para
responder à luz?
| Tabela 2. Sondas diferencialmente metiladas associadas à latência do reflexo pupilar à luz | |||||
| p < 2.4 × 10⁻⁷ | |||||
| Sonda | Fenótipo de latência (meses) | Localização genômica (hg19) | Anotação gênica Illumina | Tamanho do efeito (β) | Valor de p |
|---|---|---|---|---|---|
| cg05148717 | 14 | chr6:28829171 | −0.00055 | 3.74 × 10−8 | |
| cg22367466 | 14 | chr12:6840532 | COPS7A | −0.00035 | 4.28 × 10−8 |
| cg09732535 | 24 | chr16:89331997 | 0.00066 | 1.85 × 10−8 | |
| cg15130433 | 14–24 | chr6:159240081 | EZR | 0.00025 | 2.21 × 10−7 |
| ¹ Um tamanho de efeito negativo (−β) indica hipermetilação associada a latência mais rápida (medidas transversais) ou a um aumento da latência ao longo do tempo (medidas de mudança). | |||||
Essa tabela identifica sondas CpG individuais com associações altamente significativas entre metilação do DNA e latência do reflexo pupilar à luz, evidenciando sinais epigenéticos pontuais que precedem e sustentam a análise regional por DMRs.
A Tabela 3 é o coração epigenético do artigo, porque sai do
nível “ponto a ponto” (Tabela 2) e mostra regulação coordenada em
regiões do genoma.
Ela apresenta Regiões Diferencialmente Metiladas
(DMRs) cuja metilação está associada à latência do
PLR em diferentes idades ou intervalos de
desenvolvimento.
Essa tabela sustenta que o sinal não é ruído estatístico.
| Tabela 3. DMRs significativamente associados a cada fenótipo de latência do PLR | ||||||
| Fenótipo de latência (meses) | Localização genômica (hg19) | Número de sondas | Sondas | Anotação gênica (Illumina) | Tamanho do efeito (β) | Valor de p ajustado por Bonferroni |
|---|---|---|---|---|---|---|
| 9 | chr10:74927623–74927863 | 8 | cg15213114; cg04167018; cg21416602; cg25138168; cg08571229; cg04749667; cg16124546; cg12276298 | FAM149B1; ECD | −0.007 | 7.25 × 10−3 |
| 14 | chr1:103573700–103573772 | 2 | cg26436330; cg20847625 | COL11A1 | 0.010 | 6.33 × 10−4 |
| 24 | chr10:45495981–45496216 | 3 | cg18382353; cg16512882; cg15078013 | C10orf25; ZNF22 | −0.020 | 8.49 × 10−4 |
| 9–14 | chr3:52569053–52569169 | 3 | cg18337363; cg08365687; cg13284614 | NT5DC2; LOC440957 | −0.009 | 7.35 × 10−3 |
| 14–24 | chr17:2699706–2699718 | 2 | cg05890550; cg25373595 | RAP1GAP2 | 0.020 | 6.53 × 10−5 |
| ¹ Um tamanho de efeito negativo (−β) indica hipermetilação associada a uma latência mais rápida (para medidas transversais) ou a um aumento da latência ao longo do tempo (para escores de mudança). | ||||||
| ² Esta linha detalha a DMR mais significativamente associada à mudança de latência entre 14–24 meses (ver Tabela Suplementar 2 para todas as DMRs significativas). | ||||||
Conceitualmente, a tabela 3 mostra que o PLR não é apenas um reflexo
periférico, ele reflete estados epigenéticos regionais ligados ao
neurodesenvolvimento.
Temporalmente, mostra diferentes regiões reguladas em idades distintas e
reforça a ideia de janelas epigenéticas sensíveis.
Metodologicamente, justifica análises funcionais (GO, SFARI).
As DMRs identificadas indicam que a latência do reflexo pupilar à luz está associada a padrões regionais coordenados de metilação do DNA, sugerindo regulação epigenética dinâmica de genes envolvidos no neurodesenvolvimento ao longo dos primeiros dois anos de vida.
A Tabela 4 apresenta as DMRs (Differentially Methylated Regions) mais
significativamente associadas à amplitude do reflexo pupilar à luz (PLR)
em diferentes idades ou intervalos de desenvolvimento.
Ela confirma e fortalece os achados CpG-por-CpG, mostrando que as
alterações não são isoladas, mas fazem parte de padrões regionais
biologicamente coerentes.
| Tabela 4. DMRs mais significativamente associados a cada fenótipo de amplitude do PLR | ||||||
| Fenótipo de amplitude (meses) | Localização genômica (hg19) | Número de sondas | Sondas | Anotação gênica Illumina | Tamanho do efeito (β) | Valor de p ajustado por Bonferroni |
|---|---|---|---|---|---|---|
| 9 | chr13:103452556–103453215 | 4 | cg06518779; cg21251000; cg15186648; cg15193473 | BIVM; KDELC1 | 0.04 | 0.000213 |
| 14 | chr11:111249659–111250201 | 5 | cg19126910; cg17390301; cg24049888; cg18316498; cg11362935 | POU2AF1 | -0.03 | 0.022200 |
| 24 | chr5:77253833–77253990 | 3 | cg25051331; cg09048251; cg07595776 | NA | -0.08 | 0.000823 |
| 9–14 | chr10:102295134–102295549 | 5 | cg07690778; cg26303175; cg07080220; cg08314679; cg07510080 | HIF1AN | 0.07 | 0.002670 |
| 14–24 | chr15:89786761–89787223 | 4 | cg22813622; cg01741397; cg15769724; cg06870609 | FANCI | -0.04 | 0.001220 |
Aa alterações regionais de metilação em genes e regiões regulatórias chave estão associadas à amplitude do reflexo pupilar à luz ao longo do desenvolvimento infantil, indicando que mecanismos epigenéticos coordenados participam da maturação das respostas sensoriais precoces: um fenótipo intermediário relevante para o autismo.
| Níveis de análise epigenômica apresentados nas tabelas do estudo | ||
| Tabela | Nível de análise | O que mostra |
|---|---|---|
| Tabela 2 | CpG individual (DMP) | Associações pontuais |
| Tabela 3 | Integração / anotação | Contexto funcional |
| Tabela 4 | Regiões (DMRs) | Efeitos coordenados e robustos |
Aqui, é onde o estudo “sobe um nível” e transforma CpGs/DMRs em
processos biológicos: ela mostra os termos de Ontologia Gênica (GO) que
estão significativamente enriquecidos entre os genes associados aos
fenótipos do reflexo pupilar à luz (PLR):
Latência do PLR: quão rápido a pupila responde.
Amplitude do PLR: quão forte é a contração.
Então, em vez de perguntar “qual CpG mudou?”, a aqui pergunta é:“Quais processos biológicos aparecem mais do que o esperado entre os genes associados ao PLR?”
| Tabela 5. Termos de Ontologia Gênica (GO) significativamente enriquecidos para cada fenótipo do PLR | ||||
| Número de termos GO significativos | Termo GO mais significativamente enriquecido (nº GO) | Enriquecimento (fold) | Valor de p (FDR) | |
|---|---|---|---|---|
| Latência do PLR | ||||
| 9 | 12 | Adesão celular homofílica via moléculas de adesão da membrana plasmática (GO:0007156) | 4.20 | 3.17 × 10−7 |
| 14 | 34 | Regulação negativa de processos do desenvolvimento (GO:0051093) | 1.84 | 6.95 × 10−3 |
| 24 | 0 | – | – | – |
| 9–14 | 16 | Adesão celular homofílica via moléculas de adesão da membrana plasmática (GO:0007156) | 6.62 | 8.27 × 10−22 |
| 14–24 | 16 | Regulação da organização de componentes celulares (GO:0051128) | 1.50 | 5.00 × 10−3 |
| Amplitude do PLR | ||||
| 9 | 0 | – | – | – |
| 14 | 73 | Regulação da atividade de glicuronosiltransferases (GO:1904223) | 24.41 | 3.86 × 10−8 |
| 24 | 0 | – | – | – |
| 9–14 | 22 | Desenvolvimento de sistemas (GO:0048731) | 1.41 | 1.30 × 10−3 |
| 14–24 | 0 | – | – | – |
| sig. = significativo após ajuste por FDR (False Discovery Rate). | ||||
A latência do PLR está fortemente associada a processos de adesão
celular, organização estrutural e regulação do desenvolvimento neural,
especialmente no primeiro ano de vida. Já a amplitude do PLR apresenta
associações mais pontuais, destacando processos metabólicos e de
desenvolvimento sistêmico em janelas específicas.
Ou seja, a Latência mostra a conectividade e arquitetura neural e a
Amplitude mostra a modulação funcional e metabólica.
A Tabela 6 é uma das mais importantes do trabalho, porque ela faz a
ponte direta entre os achados epigenéticos do PLR e genes já implicados
no autismo.
Ela lista Regiões Diferencialmente Metiladas (DMRs) que estão
significativamente associadas a fenótipos do PLR (latência ou
amplitude),coincidem com genes catalogados no banco SFARI, que reúne
evidências genéticas para autismo, permanecem significativas após
correções estatísticas rigorosas.
Não são CpGs isolados, mas sinais epigenéticos regionais que atingem
genes com relevância conhecida para o TEA.
| Tabela 6. DMRs significativamente associados anotados a genes relacionados ao autismo | |||||
| Localização genômica do DMR (hg19) | Sondas | Fenótipo associado | Gene | Pontuação SFARI do gene | Nível de evidência |
|---|---|---|---|---|---|
| chr11:19372012–19372234 | cg23330281; cg24137774; cg25909885 | Latência 14–24 meses | NAV2 | 2 | Candidato forte |
| chr7:2968559–2968595 | cg22989995; cg23770265 | Amplitude 24 meses | CARD11 | 2 | Candidato forte |
| chr13:38444227–38444490 | cg16409955; cg23021771 | Amplitude 9–14 meses | TRPC4 | 3 | Evidência sugestiva |
| ¹ A pontuação SFARI reflete a força da evidência da relevância do gene para o autismo. | |||||
Os genes NAV2, TRPC4 e CARD11 atuam em níveis complementares do neurodesenvolvimento - organização estrutural, excitabilidade neuronal e modulação neuroimune - oferecendo um arcabouço biológico plausível para a associação entre metilação do DNA, reflexo pupilar e risco para TEA.
| Integração funcional dos três genes analisados | ||
| Gene | Nível principal | O que regula |
|---|---|---|
| NAV2 | Estrutural | Montagem e orientação dos circuitos neurais |
| TRPC4 | Funcional | Excitabilidade neuronal e sinalização intracelular |
| CARD11 | Modulador | Estado inflamatório e modulação autonômica |
Sondas são uma sequência curta de DNA
sintético, fixada no chip, que se liga (hibridiza) a um local específico
do genoma, permitindo medir o nível de metilação de uma citosina
(CpG) naquele ponto.
As sondas capturam regulação gênica, não mutações.
A sonda gera um valor chamado β - beta value - cujo cálculo é feito da
seguinte forma:
β = metilação relativa (0 = não metilado | 1 = totalmente
metilado), cujo resultado significa:
β negativo: região mais metilada em um grupo específico.
β positivo: região menos metilada.
O que seriam várias sondas juntas num estudo?
Na tabela, aparecem cg23330281/cg24137774/cg25909885 e
isso significa que todas essas sondas CpG estão:
- Próximas no genoma;
- Dentro da mesma região diferencialmente metilada (DMR);
- E que juntas, apresentam alteração significativa de metilação,
reforçando que não é ruído, mas sim um sinal epigenético
regional (e é exatamente isso que define uma DMR Differentially
Methylated Region)
Como as sondas capturam regulação gênica e não as mutações, então,
muitas DMRs são encontradas em:
- Genes de adesão celular;
- Axon guidance;
- Desenvolvimento sináptico.
Isso significa: calibração sensorial do cérebro em
formação, que é exatamente o que esse artigo
demonstra.
“Reforçando que não é ruído, mas um sinal epigenético
regional”
Mas Jânice, o que isso significa?
Vamos por partes.
1. O que significa “ruído” nesse contexto?
Em dados de metilação, ruído significa:
- Variação aleatória do sinal;
- Erro técnico do array;
- Flutuação estatística em um único CpG;
- Efeito espúrio que não se replica nem tem coerência biológica.
Um exemplo de ruído: Uma única sonda CpG aparece “significativa”,
mas as sondas vizinhas não mostram nada parecido.
(Isso é comum em EWAS (Epigenome-Wide Association Study) e não
costuma ter valor biológico robusto.)
2. O que caracteriza um sinal epigenético real?
Um sinal epigenético biologicamente relevante apresenta:
- Múltiplas sondas CpG próximas entre si;
- Direção consistente do efeito (todas mais ou menos metiladas);
- Associação estatística após correção para múltiplos testes e
- Localização coerente com regiões regulatórias, promotors, enhancers e
corpos gênicos.
Isso indica regulação coordenada, não acaso.
3. Mas por que “regional”?
A alteração não ocorre em um único ponto mas sim em uma região do
genoma, que envolve vários CpGs contíguos e que costumam ser regulados
em conjunto.
A metilação não atua ponto a ponto, mas por blocos
regulatórios.
(Que é exatamente a definição de uma DMR.)
4.Por que várias sondas reforçam que não é ruído?
Porque estatisticamente a probabilidade de vários CpGs adjacentes
mudarem na mesma direção com significância no mesmo grupo por acaso é
extremamente baixa.
Biologicamente, isso reflete a ação de:
- DNMTs (DNA Methyltransferases);
- TETs(Ten-Eleven Translocation enzymes);
- Remodeladores de cromatina atuando sobre uma região funcional
inteira.
(Ou seja, um programa regulatório ativo, não flutuação
aleatória.)
Qual a Relação entre sondas, DMPs e
DMRs?
| Relação entre sondas, DMPs e DMRs | |
| Termo | Definição |
|---|---|
| Sonda (probe) | Mede o nível de metilação de um único sítio CpG |
| DMP (Differentially Methylated Position) | Sítio CpG individual com diferença de metilação estatisticamente significativa entre grupos |
| DMR (Differentially Methylated Region) | Região genômica contendo múltiplos sítios CpG com alterações coordenadas de metilação |
| CpG: região do DNA onde uma citosina é seguida por uma guanina; a metilação nesses sítios é um importante mecanismo de regulação epigenética. | |
EWAS (Epigenome-Wide Association Study) é
um estudo que busca associações entre marcas epigenéticas no DNA
(principalmente metilação) e fenótipos biológicos ou clínicos.
Um EWAS analisa, ponto a ponto do DNA, se diferenças na metilação estão
relacionadas a um determinado fenótipo, usando:
- Dados de arrays Illumina;
- Modelos estatísticos (geralmente regressão linear) e
- Ajustes para idade, sexo, composição celular, batch e múltiplos
testes.
Nesse estudo, o uso de EWAS permitiu detectar assinaturas epigenéticas
muito precoces antes do diagnóstico clínico de TEA, ligadas a uma função
sensorial básica (PLR) sem depender de mutações raras. Isso sustenta
fortemente a ideia de que o autismo pode envolver diferenças
regulatórias iniciais e não apenas alterações genéticas
fixas.
DMPs e DMRs dentro do EWAS
DMP (Differentially Methylated Position): um único CpG
associado ao fenótipo.
DMR (Differentially Methylated Region): conjunto de
CpGs próximos sinal mais robusto e biologicamente interpretável.
As DMRs dentro do EWAS são especialmente valorizadas porque reduzem
falso positivo e refletem regulação regional real.
EWAS é diferente de causalidade direta!
EWAS identifica associações e não prova causalidade por si só, mas
revela mecanismos plausíveis que podem ser testados em modelos
celulares, estudos longitudinais e integração com genética.
Ontologia Gênica é um sistema padronizado
de classificação usado em biologia molecular e bioinformática para
descrever, de forma consistente, o que os genes e seus produtos (RNAs,
proteínas) fazem, onde atuam e em quais processos biológicos participam
- independentemente da espécie.
Ela funciona como um vocabulário controlado + uma estrutura
hierárquica, permitindo comparar genes entre organismos diferentes e
interpretar grandes conjuntos de dados (RNA-seq, proteômica, metilação
etc.).