A deficiência de Glicose 6 Fosfato Desidrogenase (G6PD) é a enzimopatia mais comum, e sua base genética na população brasileira é reconhecidamente complexa. No entanto, a distribuição regional dos haplótipos do gene G6PD e o perfil de risco da sua deficiência enzimática ainda são pouco caracterizados, justificando a necessidade de estudos populacionais detalhados.

Caracterizar a heterogeneidade haplotípica do gene G6PD e a distribuição de haplótipos de relevância clínica nas cinco macrorregiões brasileiras, analisando a diversidade e a distribuição de haplótipos patogênicos.

Este estudo retrospectivo analisou a heterogeneidade haplotípica do gene G6PD a partir de um banco de dados laboratoriais anonimizados. Foram incluídas 915 crianças (456 meninas, 459 meninos) que realizaram triagem neonatal por NGS entre janeiro de 2021 e janeiro de 2024, a partir de amostras de mucosa oral, totalizando 1.372 haplótipos analisados. Do sistema de informação laboratorial, foram recuperados os dados demográficos (sexo, logradouro) e as variantes genéticas de cada indivíduo. A determinação dos haplótipos foi então realizada a partir das 20 variantes distintas contidas na região genômica do gene G6PD, utilizando o software BEAGLE (v5.5). Para cada macrorregião brasileira, foram calculadas a diversidade haplotípica (h), a frequência e o número de tipos de haplótipos patogênicos. O estudo recebeu aprovação ética (CAAE: 79180324.7.0000.0368).

A distribuição dos haplótipos por macrorregião foi: Centro-Oeste (CO) n = 475, Norte (N) n = 389, Sudeste (SE) n = 273, Nordeste (NE) n = 122 e Sul (S) n = 112. A diversidade haplotípica (h) e o número de haplótipos distintos variaram entre as regiões, sendo maior na região N (h = 0.450; 14-tipos), seguida por NE (h = 0.390; 6-tipos), CO (h = 0.386; 13-tipos), S (h = 0.357; 8-tipos), e menor na SE (h = 0.308; 7-tipos). Dos 22 haplótipos únicos identificados, 9 foram classificados como patogênicos, e sua prevalência combinada demonstrou um gradiente de risco, com as seguintes prevalências: NE (n = 5; 4.10%), N (n = 10; 2.57%), SE (n = 6; 2.20%), CO (n = 10; 2.11%) e S (n = 2; 1.79%). O perfil de risco também se mostrou distinto em sua natureza: as regiões N e CO apresentaram maior complexidade (5-tipos de haplótipos patogênicos cada), seguidas por SE e S (2-tipos cada), enquanto o risco na região NE foi concentrado em apenas 1 tipo. Um único haplótipo patogênico, que inclui a variante 202G>A, foi o mais prevalente em quatro das cinco regiões, com as seguintes frequências: NE (n = 5; 4.10%), SE (n = 5; 1.83%), N (n = 5; 1.29%), CO (n = 6; 1.26%) e S (n = 1; 0.89%).

A estrutura genética regionalizada do gene G6PD, evidenciada pela variação na diversidade haplotípica entre N (14-tipos) e S (7-tipos), pode refletir a distinta história demográfica do país. Essa heterogeneidade se estende ao perfil de risco clínico: a maior prevalência no NE (4.10%) contrasta com sua baixa complexidade (1-tipo patogênico), enquanto N e CO apresentam maior diversidade de risco (5-tipos cada). Este cenário sugere que a eficácia dos ensaios diagnósticos pode variar regionalmente. Isso reforça a pertinência de se avaliar estratégias de triagem adaptadas a cada localidade. Em conclusão, existe heterogeneidade na distribuição dos haplótipos do gene G6PD entre as macrorregiões brasileiras, que se manifesta tanto na variação da diversidade genética geral quanto no perfil de risco clínico, que difere em prevalência, complexidade e no tipo de haplótipo patogênico predominante em cada localidade.