Functional Genomic Analyses Reveal transfer RNA Modification Enzymes as Risk Genes for Bipolar I Disorder and Schizophrenia

Resumo Técnico: Análises Genômicas Funcionais de Enzimas de Modificação de RNA no Transtorno Bipolar Tipo I e Esquizofrenia

Problema
Embora os Estudos de Associação de Genoma Amplo (GWAS) tenham identificado milhares de variantes genéticas associadas a transtornos neuropsiquiátricos, como a Esquizofrenia (SCZ) e o Transtorno Bipolar Tipo I (BIPI), traduzir esses sinais genéticos amplos em mecanismos biológicos específicos e tratáveis permanece um desafio significativo. A regulação epigenética é um eixo de vulnerabilidade conhecido, contudo, o papel do "epitranscritoma" — especificamente as modificações de RNA mediadas por proteínas modificadoras de RNA (RMPs) — não foi sistematicamente avaliado no contexto do risco de doenças neuropsiquiátricas. Embora mais de 100 RMPs tenham sido identificadas, sua contribuição para a arquitetura genética de transtornos psiquiátricos permanece em grande parte inexplorada, em parte porque elas não foram categorizadas como um conjunto de genes funcionais para análises de vias padrão.

Metodologia
Os autores empregaram um pipeline genômico funcional de múltiplos estágios para priorizar sistematicamente genes RMP subjacentes ao risco para 16 transtornos neuropsiquiátricos. O estudo utilizou um conjunto curado de 123 RMPs humanas correspondentes a 31 tipos de modificação de RNA em seis biotipos de RNA.

O fluxo de trabalho analítico integrou as seguintes abordagens:

1. Teste de Associação Baseado em Genes: Utilizou MAGMA e hMAGMA (incorporando dados de interação de cromatina Hi-C) para testar associações de conjuntos de genes entre RMPs e estatísticas de resumo de GWAS do Psychiatric Genomics Consortium (PGC).
2. Estimativa de Hereditariedade: Aplicou GCTB–SBayesR para estimar a contribuição dos corpos gênicos de RMPs individuais para a hereditariedade da doença.
3. Mapeamento Fino (Fine-Mapping): Utilizou PolyFun e SuSiE para identificar variantes causais dentro de loci de risco, incorporando estimativas de hereditariedade por SNP.
4. Análises Integradas com QTL:
   • SMR (Summary-data-based Mendelian Randomization): Integrou dados de GWAS com dados de eQTL e sQTL de córtex cerebral (do BrainMeta) para testar associações conjuntas, seguido de testes HEIDI para distinguir pleiotropia de ligação (linkage).
   • TWAS (Transcriptome-Wide Association Studies): Utilizou FUSION com modelos de Córtex Cerebral do GTEx v8 para testar associações entre a expressão gênica geneticamente prevista e o risco da doença.
   • MR de Duas Amostras: Validou a consistência direcional usando IVW, MR-Egger e outros métodos.
   • Colocalização: Aplicou colocalização Bayesiana (coloc) para determinar se os sinais de GWAS e QTL compartilham variantes causais comuns (PP.H4 > 0,6).
5. Validação: Validou os candidatos priorizados usando dados de RNA-seq bulk da coorte BrainSeq Fase II (hipocampo, córtex pré-frontal dorsolateral e núcleo caudado) via análise de expressão diferencial e regressão logística ajustada para covariáveis.

Principais Contribuições e Resultados
O estudo identificou o NSUN2, uma citosina-5 metiltransferase de RNA, como um locus de risco compartilhado reprodutível tanto para BIPI quanto para SCZ.

• Associações do NSUN2: Embora SNPs individuais não tenham atingido significância em todo o genoma, o NSUN2 mostrou associações significativas ao nível de gene em testes MAGMA e hMAGMA para ambos os transtornos. As análises SMR e TWAS revelaram uma associação negativa entre a expressão do NSUN2 e o risco da doença, sugerindo um papel neuroprotetor. A análise de colocalização apoiou variantes causais compartilhadas entre os eQTLs do NSUN2 e o risco de BIPI.
• Mecanismos de Splicing: O estudo destacou que o risco é mediado pelo splicing alternativo. sQTLs específicos que levam ao salto de éxon (produzindo isoformas não codificantes ou truncadas) foram positivamente associados ao risco da doença, enquanto eventos que promovem isoformas codificantes de proteínas foram negativamente associados. A validação por RNA-seq bulk confirmou a redução da expressão de isoformas codificantes de NSUN2 em cérebros de pacientes com SCZ (especificamente no hipocampo masculino e no córtex pré-frontal dorsolateral feminino).
• Envolvimento de Enzimas de tRNA Mais Amplas: Além do NSUN2, o estudo identificou associações significativas para outras enzimas de modificação de tRNA:
  • NSUN6 e TYW5 atingiram significância em todo o genoma para SCZ.
  • TRMT61A e QTRT1 mostraram associações significativas com BIPI ou SCZ.
  • MRM1 foi associado ao BIPI.
  • Análises integrativas (incluindo dados do COLOCdb) confirmaram associações para esses genes através de múltiplas modalidades de QTL (eQTL, sQTL e mQTL).
• Priorização: Um sistema de pontuação ponderada classificou o NSUN2 como o gene de maior prioridade para BIPI e SCZ, seguido por NSUN6 e TYW5 para SCZ.

Significância e Alegações
O artigo alega que os riscos genéticos comuns convergem para variações regulatórias que afetam a maquinaria de modificação de RNA, nomeando especificamente a desregulação da modificação de tRNA como uma via molecular que contribui para a suscetibilidade à doença psiquiátrica.

As principais conclusões incluem:

• Mecanismo Convergente: Os achados sugerem que as enzimas de modificação de tRNA (escritores de modificações m5C, m1A, yW e Q) representam um mecanismo molecular convergente subjacente à patofisiologia do BIPI e SCZ, distinto das vias de m6A de mRNA mais comumente estudadas.
• NSUN2 como um Gene de Risco Compartilhado: O NSUN2 é identificado como um gene de risco compartilhado onde o risco genético atua através da desregulação da expressão e do splicing específicos de isoformas, em vez da abundância total do gene.
• Papel Neuroprotetor: Os dados apoiam um modelo onde a função intacta do NSUN2 é neuroprotetora, e sua desregulação (via redução de isoformas codificantes) contribui para a suscetibilidade à doença.
• Limitações: Os autores observam modestamente que a lista atual de RMPs não é abrangente e que a heterogeneidade da expressão das RMPs entre regiões cerebrais e tipos celulares continua sendo um desafio. Eles enfatizam que, embora o estudo priorize candidatos, é necessária validação adicional em modelos animais e organoides cerebrais para definir seus papéis específicos no neurodesenvolvimento.

O estudo não propõe novas aplicações clínicas ou propostas experimentais futuras além da necessidade de validar esses candidatos em sistemas de modelo para confirmar seus papéis nos processos de neurodesenvolvimento.