Ocn Tesis
Ocn Tesis
Ocn Tesis
Memoria presentada por Osmel Companioni Nápoles para optar al grado de Doctor
por la Universidad de Barcelona.
A Mayelín
AGRADECIMIENTOS
ÍNDICE
ÍNDICE DE FIGURAS ....................................................................................................................... 1
ÍNDICE DE TABLAS ......................................................................................................................... 3
HIPÓTESIS ................................................................................................................................... 61
OBJETIVOS .................................................................................................................................. 65
MATERIALES Y MÉTODOS ............................................................................................................ 69
1. ESTUDIO DE ASOCIACIÓN DE GENES DE SEÑALIZACIÓN DE HELICOBACTER PYLORI Y RIESGO DE CÁNCER GÁSTRICO ....... 71
1.1. Pacientes y muestras ............................................................................................................ 71
1.2. Selección de tagSNPs en genes de la vía de señalización de Helicobacter pylori.................. 71
1.3. Extracción de ADN y genotipado .......................................................................................... 72
1.4. Infección por Helicobacter pylori .......................................................................................... 73
1.5. Análisis estadísticos y bioinformáticos ................................................................................. 73
2. ESTUDIO DE ASOCIACIÓN GENÉTICA EN EL SEGUIMIENTO DE LESIONES PRECURSORAS GÁSTRICAS ............................ 74
2.1. Pacientes y muestras ............................................................................................................ 74
2.2. Histopatología ...................................................................................................................... 75
2.3. Infección por Helicobacter pylori .......................................................................................... 75
2.4. Prueba piloto de colecta, conservación y extracción de ADN ............................................... 76
2.4.1. Medición de cantidad y calidad del ADN ....................................................................................... 77
2.5. Diseño de solución estabilizadora de saliva y prueba de conservación ................................ 78
2.5.1. Protocolo de extracción de ADN a partir de saliva ........................................................................ 78
2.5.2. Prueba piloto de genotipado de ADN de saliva ............................................................................. 79
2.5.3. Colecta de saliva en hospitales ...................................................................................................... 79
2.6. Extracción de ADN de tejido gástrico parafinado ................................................................. 80
2.7. Extracción de ADN de sangre ................................................................................................ 81
2.8. Selección de genes candidatos.............................................................................................. 81
2.9. Selección de tagSNPs en los genes candidatos y análisis bioinformáticos............................ 86
2.10. Genotipado de SNPs............................................................................................................ 86
2.11. Análisis estadísticos de los resultados de genotipado de SNPs........................................... 86
3. ESTUDIO DE EXPRESIÓN GÉNICA POR MICROARRAY EN METAPLASIA INTESTINAL QUE PROGRESA A CÁNCER GÁSTRICO . 94
3.1. Pacientes y muestras ............................................................................................................ 94
3.2. Corte de bloques de parafina y evaluación histológica......................................................... 98
3.3. Extracción de ARN total ........................................................................................................ 98
3.4. Medición de cantidad y calidad del ARN............................................................................... 99
3.5. Prueba piloto de comparación de microarrays Human Gene 2.0 ST vs Almac Xcel ............ 100
3.6. Análisis de calidad y estadísticos ........................................................................................ 100
3.7. Análisis de enriquecimiento funcional por IPA y GSEA........................................................ 101
3.8. Validación de expresión génica por qRT-PCR ...................................................................... 103
3.8.1. Protocolo de qRT-PCR por sondas UPL y cebadores en plataforma Biomark .............................. 109
3.8.2. Procesamiento de los resultados de qRT-PCR por sondas UPL y cebadores en plataforma
Biomark ................................................................................................................................................. 110
4. ESTUDIO DE EXPRESIÓN GÉNICA E INMUNO-HISTOQUÍMICA DE SCHLAFEN5 EN LESIONES PRECURSORAS Y CÁNCER
GÁSTRICO ............................................................................................................................................... 112
4.1. Tissue arrays comerciales de tejido gástrico....................................................................... 112
4.2. Pacientes y muestras analizadas ........................................................................................ 112
4.3. Histopatología .................................................................................................................... 113
Índice
1. ESTUDIO DE ASOCIACIÓN DE GENES DE SEÑALIZACIÓN DE HELICOBACTER PYLORI Y RIESGO DE CÁNCER GÁSTRICO ..... 197
1.1. NOD2. Asociación genética y plausibilidad biológica .......................................................... 197
1.2. CD14. Asociación genética y plausibilidad biológica ........................................................... 200
1.3. TLR4 y NFKB1. Asociación genética y plausibilidad biológica ............................................. 202
2. ESTUDIO DE ASOCIACIÓN GENÉTICA EN EL SEGUIMIENTO DE LESIONES PRECURSORAS GÁSTRICAS.......................... 204
2.1. Asociación genética con la evolución de las lesiones precursoras gástricas ....................... 205
3. ESTUDIO DE EXPRESIÓN GÉNICA EN METAPLASIA INTESTINAL QUE PROGRESA A CÁNCER GÁSTRICO ........................ 212
3.1 Genes diferencialmente expresados..................................................................................... 212
3.1.1 Subtipos histológicos de metaplasia intestinal que progresan a cáncer gástrico ......................... 212
3.1.2 Metaplasia intestinal que no progresa a cáncer gástrico en comparación con mucosa gástrica sana
............................................................................................................................................................... 215
Índice
3.2 Vías de señalización celular y procesos biológicos implicados en la progresión de los subtipos
histológicos de metaplasia intestinal a cáncer gástrico ............................................................ 218
3.2.1 Gene Set Enrichment Analysis ...................................................................................................... 219
3.2.2 Ingenuity Pathway Analysis .......................................................................................................... 229
3.3. Validación de la expresión génica por qRT-PCR .................................................................. 232
4. ESTUDIO DE EXPRESIÓN GÉNICA E INMUNO-HISTOQUÍMICA DE SCHLAFEN5 EN LESIONES PRECURSORAS Y CÁNCER
GÁSTRICO ............................................................................................................................................... 232
ÍNDICE DE FIGURAS
FIGURA 1. DISTRIBUCIÓN MUNDIAL DE LA INCIDENCIA DE CÁNCER GÁSTRICO ESTANDARIZADA POR EDAD .................... 16
FIGURA 2. ESTRUCTURA ANATÓMICA DEL ESTÓMAGO. ...................................................................................... 16
FIGURA 3. DISPOSICIÓN DE LAS CAPAS CONSTITUYENTES DEL TRACTO DIGESTIVO. ................................................... 17
FIGURA 4. ESTRUCTURA HISTOLÓGICA DE LAS GLÁNDULAS OXÍNTICAS O FÚNDICAS Y PILÓRICAS ................................. 18
FIGURA 5. ADENOCARCINOMA GÁSTRICO DE TIPO INTESTINAL ............................................................................ 20
FIGURA 6. ADENOCARCINOMA GÁSTRICO DE TIPO DIFUSO ................................................................................. 20
FIGURA 7. MUCOSA GÁSTRICA COLONIZADA POR H. PYLORI ............................................................................... 21
FIGURA 8. MECANISMO DE ACCIÓN DE VACA EN EPITELIO GÁSTRICO .................................................................... 22
FIGURA 9. MECANISMOS DE INFECCIÓN DE HELICOBACTER PYLORI EN EL EPITELIO GÁSTRICO. ................................... 25
FIGURA 10. MECANISMO DE LA CARCINOGÉNESIS GÁSTRICA .............................................................................. 29
FIGURA 11. MUCOSA NORMAL TEÑIDA CON HEMATOXILINA-EOSINA ................................................................... 32
FIGURA 12. GASTRITIS CRÓNICA NO ATRÓFICA ................................................................................................ 32
FIGURA 13. GASTRITIS CRÓNICA ATRÓFICA ..................................................................................................... 33
FIGURA 14. METAPLASIA INTESTINAL COMPLETA ............................................................................................. 34
FIGURA 15. METAPLASIA INTESTINAL INCOMPLETA .......................................................................................... 34
FIGURA 16. DISPLASIA GÁSTRICA .................................................................................................................. 35
FIGURA 17. SCHLAFEN4 SE SOBRE-EXPRESA EN LA METAPLASIA GÁSTRICA. ............................................................ 45
FIGURA 18. MICROARRAYS DE TEJIDO USADOS PARA INMUNOHISTOQUÍMICA DE SCHLAFEN5. ................................ 112
FIGURA 19. ELECTROFORESIS DE MUESTRAS DE ADN EXTRAÍDAS DE HISOPOS BUCALES DRI‐CAPSULES. .................... 128
FIGURA 20. ELECTROFORESIS DE MUESTRAS DE ADN EXTRAÍDAS DE HISOPOS BUCALES BUCCALAMP™. ................... 129
FIGURA 21. ELECTROFORESIS DE MUESTRAS DE ADN EXTRAÍDAS DE SALIVA EMPLEANDO EL KIT ORAGENE................ 129
FIGURA 22. ELECTROFORESIS DEL ADN EXTRAÍDO DE SALIVA CONSERVADAS EN SOLUCIÓN ESTABILIZADORA # 4........ 131
FIGURA 23. ELECTROFORESIS DEL ADN EXTRAÍDO DE SALIVA COLECTADA EN LOS HOSPITALES EN SOLUCIÓN
ESTABILIZADORA # 4 ................................................................................................................................. 131
FIGURA 24. GENOTIPADO DEL SNP RS7578034 EN EL GEN IL1B EMPLEANDO LIGHTCYCLER II. ............................. 132
FIGURA 25. CLUSTER DE CALIDAD MEDIA EN EL GENOTIPADO SEQUENOM........................................................... 133
FIGURA 26. CLUSTERS DE CALIDAD ÓPTIMA EN EL GENOTIPADO SEQUENOM. ...................................................... 133
FIGURA 27. ELECTROFORESIS DE ARN TOTAL OBTENIDO DE MUESTRAS PARAFINADAS EN BIOANALYZER 2100 .......... 145
FIGURA 28. DENDROGRAMA DE TODOS LOS GRUPOS ANALIZADOS. ................................................................... 148
FIGURA 29. HEATMAP DE LA COMPARACIÓN MII-CG VS MII-NO CG. .............................................................. 149
FIGURA 30. HEATMAP DE LA COMPARACIÓN MIC-CG VS MIC-NO CG. ............................................................ 149
FIGURA 31. HEATMAP DE LA COMPARACIÓN MI-NO CG VS SANOS. ................................................................. 150
FIGURA 32. DIAGRAMA DE VENN DE LOS GENES SIGNIFICATIVOS EN LAS COMPARACIONES RELATIVAS A LA MII. ......... 151
FIGURA 33. DIAGRAMA DE VENN DE LOS GENES SIGNIFICATIVOS EN LAS COMPARACIONES RELATIVAS A LA MIC. ........ 152
FIGURA 34. NÚMERO DE GENES SOBRE Y SUB-EXPRESADOS EN LOS GRUPOS ANALIZADOS. ..................................... 152
FIGURA 35. ANÁLISIS DE GENES LEADING EDGE CON LOS CONJUNTOS GÉNICOS SOBRE-EXPRESADOS EN MII-CG......... 169
FIGURA 36. SUBANÁLISIS DE GENES LEADING EDGE DEL NÚCLEO 1 (DEGRADACIÓN PROTEASOMAL). ........................ 169
FIGURA 37. HEAT MAP DE LOS GENES LEADING EDGE DEL NÚCLEO 1. ................................................................. 170
FIGURA 38. SUBANÁLISIS DE GENES LEADING EDGE DEL NÚCLEO 2 (FOSFORILACIÓN OXIDATIVA).............................. 170
FIGURA 39. ESQUEMAS DE LAS VÍAS CANÓNICAS SOBRE-EXPRESADAS EN LOS GRUPOS MII-CG Y MIC-CG POR IPA... 173
FIGURA 40. REDES DE ALTO SCORE IPA DE LAS COMPARACIONES MII-CG VS MII-NO CG Y MI-NO CG VS SANOS. .. 175
FIGURA 41. HEATMAP DE LA QRT-PCR EN LA PLATAFORMA BIOMARK HD 96X96 .............................................. 180
FIGURA 42. GRÁFICO DE CAJAS DE LOS VALORES DE CT DEL GEN ACTB EN LOS GRUPOS MI-CG, MI-NO CG Y SANOS. 181
FIGURA 43. GRÁFICA DEL NIVEL DE EXPRESIÓN DEL MICROARRAY Y LA QRT-PCR DE LOS GENES SIGNIFICATIVOS. ........ 182
FIGURA 44. INDUCCIÓN DE SCHLAFEN5 POR EL IFNΑ EN LAS LÍNEAS CELULARES HL-60 Y JURKAT............................ 185
FIGURA 45. ANÁLISIS INMUNOHISTOQUÍMICO DE SCHLAFEN5 TRAS INDUCCIÓN POR IFNΑ EN LA LÍNEA JURKAT. ........ 185
FIGURA 46. INMUNO-HISTOQUÍMICA DE LA PROTEÍNA SCHLAFEN5.................................................................... 187
1
Índice de Figuras
FIGURA 47. INMUNO-HISTOQUÍMICA DE LA PROTEÍNA SCHLAFEN5 EN LAS MUESTRAS ESTUDIADAS ......................... 189
FIGURA 48. REPRESENTACIÓN GRÁFICA DE LA TINCIÓN ESTROMAL DE SCHLAFEN5................................................ 191
FIGURA 49. CURVAS DE SOBREVIVENCIA PARA LOS MODELOS MI-CG CON Y SIN ADICIÓN DE SCORE H ESTROMAL ...... 192
FIGURA 50. CURVAS DE SOBREVIVENCIA PARA LOS MODELOS MII-CG CON Y SIN ADICIÓN DE SCORE H ESTROMAL ..... 193
FIGURA 51. DOBLE MARCAJE DE SCHLAFEN5 CON MARCADORES DE SUPERFICIE. ................................................. 194
FIGURA 52. BLOQUES DE HAPLOTIPOS DE MUC2.......................................................................................... 208
FIGURA 53. DIAGRAMA DE VENN DE PROCESOS FUNCIONALES SOBRE-EXPRESADOS COMUNES Y ESPECÍFICOS ENTRE LOS
GRUPOS MII-CG Y MI-NO CG, SEGÚN LOS ANÁLISIS GSEA E IPA. ................................................................... 220
FIGURA 54. PROCESOS FUNCIONALES COMUNES A CONJUNTOS GÉNICOS CON DIFERENTES CRITERIOS DE SELECCIÓN EN LA
MII-CG................................................................................................................................................. 224
FIGURA 55. PROCESOS FUNCIONALES COMUNES A CONJUNTOS GÉNICOS CON DIFERENTES CRITERIOS DE SELECCIÓN EN LA
MI-NO CG. ........................................................................................................................................... 228
2
Índice de Tablas
ÍNDICE DE TABLAS
TABLA 1. TASAS DE TRANSICIÓN PARA LA PROGRESIÓN Y REGRESIÓN DE LAS LESIONES PRECURSORAS ......................... 28
TABLA 2. PRINCIPALES ESTUDIOS EPIDEMIOLÓGICOS DE SEGUIMIENTO DE LESIONES PRECURSORAS GÁSTRICAS. ............ 38
TABLA 3. RESUMEN DE ESTUDIOS DE ASOCIACIÓN DE GENOMA COMPLETO EN CÁNCER GÁSTRICO. ............................. 52
TABLA 4. ESTUDIOS DE ASOCIACIÓN GENÉTICA EN LESIONES PRECURSORAS DE CÁNCER GÁSTRICO. ............................. 56
TABLA 5. SOLUCIONES DE COLECTA Y ESTABILIZADORAS DE SALIVA....................................................................... 78
TABLA 6. GENES CANDIDATOS SELECCIONADOS A PARTIR DE ESTUDIOS DE MODELOS MURINOS DE CARCINOGÉNESIS
GÁSTRICA. ................................................................................................................................................ 83
TABLA 7. GENES CANDIDATOS SELECCIONADOS A PARTIR DEL ANÁLISIS DE LOS MECANISMOS DE INFECCIÓN DE
HELICOBACTER PYLORI. ............................................................................................................................... 84
TABLA 8. GENES CANDIDATOS SELECCIONADOS A PARTIR DE ESTUDIOS DE ASOCIACIÓN GENÉTICA CON CÁNCER GÁSTRICO
(GWAS) Y CON LESIONES PRECURSORAS GÁTRICAS (GENES CANDIDATOS). ............................................................ 85
TABLA 9. SNPS GENOTIPADOS CON SU INFORMACIÓN ASOCIADA. ....................................................................... 88
TABLA 10. MUESTRAS USADAS PARA PRUEBA PILOTO DE COMPARACIÓN DE MICROARRAYS DE EXPRESIÓN. ................. 95
TABLA 11. MUESTRAS USADAS EN EL ESTUDIO DE EXPRESIÓN DIFERENCIAL POR MICROARRAYS DE LOS SUBTIPOS
HISTOLÓGICOS DE MI QUE PROGRESAN A CG, RESPECTO A LOS QUE NO PROGRESAN Y A LA MUCOSA SANA. ................. 96
TABLA 12. GENES SIGNIFICATIVOS DEL MICROARRAY DE LAS COMPARACIONES MII-CG VS MII-NO CG, MIC-CG VS MIC-
NO CG Y MI-NO CG VS SANOS VALIDADOS POR QRT-PCR. ........................................................................... 105
TABLA 13. SERIE INDEPENDIENTE DE MUESTRAS USADAS EN LA VALIDACIÓN DE EXPRESIÓN POR QRT-PCR. ............... 107
TABLA 14. MUESTRAS DE LESIONES PRECURSORAS Y CÁNCER GÁSTRICO USADAS EN EL ESTUDIO INMUNO-HISTOQUÍMICO
DE SCHLAFEN5. ....................................................................................................................................... 113
TABLA 15. PRINCIPALES CARACTERÍSTICAS DE LA POBLACIÓN ANALIZADA EN EL ESTUDIO DE ASOCIACIÓN GENÉTICA CON
CÁNCER GÁSTRICO. ................................................................................................................................... 121
TABLA 16. CARACTERÍSTICAS DE LOS SNPS ANALIZADOS Y SU ASOCIACIÓN ESTADÍSTICA CON EL CG GLOBAL (CARDIAS + NO
CARDIAS/INTESTINAL + DIFUSO). ................................................................................................................. 124
TABLA 17. ASOCIACIÓN ESTADÍSTICA DE SNPS CON RESPECTO A LA LOCALIZACIÓN ANATÓMICA DEL CÁNCER GÁSTRICO.125
TABLA 18. ASOCIACIÓN ESTADÍSTICA DE SNPS PREVIAMENTE ASOCIADOS CON CÁNCER GÁSTRICO NO CARDIAS EN
INDIVIDUOS INFECTADOS POR HELICOBACTER PYLORI CON RESPECTO AL ESTATUS DE CAGA. ................................... 126
TABLA 19. ASOCIACIÓN ESTADÍSTICA DE SNPS CON RESPECTO AL SUBTIPO HISTOLÓGICO DE CÁNCER GÁSTRICO. ........ 126
TABLA 20. ASOCIACIÓN ESTADÍSTICA DE HAPLOTIPOS RESPECTO A LOS SUBTIPOS HISTOLÓGICOS Y ANATÓMICOS DE
CÁNCER GÁSTRICO. ................................................................................................................................... 127
TABLA 21. BÚSQUEDA DE POTENCIALES EQTLS ENTRE LOS SNPS ASOCIADOS Y NO ASOCIADOS A CG. ...................... 128
TABLA 22. PROMEDIO DE LA MASA DE ADN OBTENIDA Y PRESENCIA DE OXIDACIÓN EN LAS DIFERENTES SOLUCIONES
ESTABILIZADORAS DE SALIVA. ...................................................................................................................... 131
TABLA 23. RESULTADOS DEL GENOTIPADO DEL SNP RS7578034 EN GEN IL1B EN MUESTRAS DE ADN EXTRAÍDO DE
SALIVA. .................................................................................................................................................. 132
TABLA 24. CARACTERÍSTICAS AL RECLUTAMIENTO DE LA POBLACIÓN ESTUDIADA. ................................................. 135
TABLA 25. OR (IC 95%) DE LOS SNPS ASOCIADOS SIGNIFICATIVAMENTE CON PROGRESIÓN Y REGRESIÓN DE LAS LPGS.136
TABLA 26. ASOCIACIÓN CON LA EVOLUCIÓN DE LPGS EN EL PRESENTE Y EN EL ESTUDIO DE SEGUIMIENTO DE LPGS DE
SORIA. ................................................................................................................................................... 138
TABLA 27. ASOCIACIONES SIGNIFICATIVAS DE HAPLOTIPOS RESPECTO A LA EVOLUCIÓN DE LAS LPGS. ....................... 140
TABLA 28. ANÁLISIS DE ASOCIACIÓN EN PORTADORES DE LA COMBINACIÓN DE LOS FACTORES DE VIRULENCIA DE H. PYLORI
CAGA+VACAS1/M1 VS OTRAS COMBINACIONES, EN INDIVIDUOS INFECTADOS POR H. PYLORI. ............................... 141
TABLA 29. POTENCIALES EQTL ENTRE LOS SNPS ESTUDIADOS. ........................................................................ 143
TABLA 30. MEDICIÓN DE CANTIDAD Y CALIDAD DE LAS MUESTRAS DE ARN USADAS EN LA PRUEBA PILOTO DE
COMPARACIÓN DE MICROARRAYS HUMAN GENE 2.0 ST Y ALMAC XCEL. ............................................................ 145
TABLA 31. MEDICIÓN DE CANTIDAD Y CALIDAD DE LAS MUESTRAS USADAS EN EL MICROARRAY DE EXPRESIÓN. ........... 146
TABLA 32. GENES DIFERENCIALMENTE EXPRESADOS ENTRE LOS GRUPOS MII-CG VS MII-NO CG. .......................... 153
3
Índice de Tablas
TABLA 33. GENES DIFERENCIALMENTE EXPRESADOS ENTRE LOS GRUPOS MIC-CG VS MIC-NO CG. ........................ 157
TABLA 34. GENES DIFERENCIALMENTE EXPRESADOS ENTRE LOS GRUPOS MI-CG VS MI-NO CG. ............................ 158
TABLA 35. CONJUNTOS GÉNICOS SOBRE-EXPRESADOS EN MII-CG CON LOS VALORES MÁS EXTREMOS DEL PARÁMETRO
RANK AT MAX. ........................................................................................................................................ 162
TABLA 36. CONJUNTOS GÉNICOS SOBRE-EXPRESADOS EN MI-NO CG CON LOS VALORES MÁS EXTREMOS DEL PARÁMETRO
RANK AT MAX. ........................................................................................................................................ 164
TABLA 37. CONJUNTOS GÉNICOS COMPUESTOS POR AL MENOS 3 GENES LEADING EDGE SOBRE-EXPRESADOS EN EL GRUPO
MII-CG................................................................................................................................................. 166
TABLA 38. CONJUNTOS GÉNICOS COMPUESTOS POR AL MENOS 3 GENES LEADING EDGE SOBRE-EXPRESADOS EN EL GRUPO
MI-NO CG. ........................................................................................................................................... 167
TABLA 39. VÍAS CANÓNICAS DESREGULADAS E INFORMACIÓN FUNCIONAL OBTENIDA EN EL IPA.............................. 172
TABLA 40. REDES MOLECULARES OBTENIDAS EN EL IPA PARA LAS COMPARACIONES MII-CG VS MII-NO CG Y MI-NO CG
VS SANOS............................................................................................................................................... 177
TABLA 41. RESULTADOS DEL ANÁLISIS DE ONTOLOGÍA GÉNICA EN LAS COMPARACIONES DE ESTUDIO. ....................... 179
TABLA 42. RESULTADOS DE PRUEBA ANOVA DE LOS GENES DE REFERENCIA EN LOS GRUPOS MI-CG, MI-NO CG Y
CONTROLES SANOS................................................................................................................................... 181
TABLA 43. RESULTADOS DEL TEST DE PERMUTACIONES DEL PROGRAMA BOOTSRATIO EMPLEADO PARA ESTIMAR
DIFERENCIAS SIGNIFICATIVAS ENTRE GRUPOS MI-CG VS MI-NO CG. ................................................................. 183
TABLA 44. EXPRESIÓN GÉNICA DE SCHLAFEN5 EN LÍNEAS CELULARES EXPUESTAS A DIVERSAS CONDICIONES............... 184
TABLA 45. CUANTIFICACIÓN DE LA TINCIÓN ESTROMAL Y GLANDULAR DE SCHLAFEN5 MEDIANTE EL SCORE H EN LOS
GRUPOS ANALIZADOS. .............................................................................................................................. 190
4
Abreviaturas y Símbolos
ABREVIATURAS Y SÍMBOLOS
ADN: Ácido desoxirribonucleico.
DAB: Diaminobenzidina.
DE: Desviación estándar.
DISP: Displasia.
E: Extensión.
5
Abreviaturas y Símbolos
LPS: Lipopolisacárido.
6
Abreviaturas y Símbolos
V: Volumen.
VNTR: Variable number tandem repeat.
Vs: Volúmenes.
7
PRESENTACIÓN
Presentación
11
INTRODUCCIÓN
Introducción
15
Introducción
-Mucosa: Formada por un epitelio que forma las glándulas gástricas secretoras y
que está en contacto con el lumen gástrico. Está constituida, a su vez, por la lámina
propia o estroma, formada por vasos sanguíneos, tejido conectivo y células de
respuesta inmune como linfocitos, células plasmáticas, eosinófilos, macrófagos y
neutrófilos; y una fina capa de músculo liso denominada mucosa muscularis, que
separa la mucosa de la submucosa, y cuya contracción posibilita la secreción y
absorción de sustancias. Su epitelio superficial es formado por las foveolas gástricas
16
Introducción
hacia las cuales vierten su contenido las glándulas gástricas (entre 4 y 5 a cada
foveola) y estas a su vez vierten hacia el lumen del estómago. Además, presenta
células mucosas de superficie que secretan un mucus de pH básico, el cual protege
a la mucosa de las secreciones ácidas y sirve como una barrera física entre la
mucosa y los contenidos estomacales.
-Muscular (muscularis externa): Compuesta por dos capas de tejido muscular liso,
una capa interior orientada circularmente y otra exterior dispuesta
longitudinalmente. Entre ambas capas existe tejido conectivo que contiene al plexo
nervioso mientérico. Las contracciones de la muscularis externa bajo el control del
sistema nervioso entérico son responsables de la peristalsis, las contracciones que
producen acortamientos y alargamientos del tracto gastrointestinal para mezclar y
mover el alimento.
-Serosa: Estrato más alejado de la luz del estómago, compuesto por epitelio
escamoso simple denominado mesotelio y tejido conectivo subyacente. Grandes
vasos sanguíneos y linfáticos, así como troncos nerviosos, viajan a través de la
serosa hasta la pared del tracto digestivo. En la figura 3 se observa la disposición
histológica de las capas antes descritas.
17
Introducción
Por otra parte, cada área anatómica del estómago presenta una composición
celular característica (Figuras 2 y 4):
18
Introducción
19
Introducción
Figura 5. Adenocarcinoma gástrico de tipo intestinal. Tomado de Correa et al., 2012 14.
Figura 6. Adenocarcinoma gástrico de tipo difuso. Tomado de Correa et al., 2012 14.
20
Introducción
Figura 7. Mucosa gástrica colonizada por H. pylori (indicada por las flechas).
Tomado de Correa et al., 2012 14.
El hecho de que ciertas áreas geográficas (África, India, Malasia, India, China,
Colombia y Costa Rica) presenten elevada prevalencia de infección por H. pylori y
baja incidencia de CG se conoce como el enigma africano. Las posibles explicaciones
a este fenómeno han sido diferencias en la virulencia oncogénica de cepas de
H. pylori, co-infección con helmintos y desplazamiento de la respuesta inmune
desde Th1 a Th2 16. Existen evidencias en modelos animales y en humanos de que la
respuesta inmune ante la infección de Helicobacter puede ser cambiada de un típico
patrón Th1 a Th2 debido a la co-infección con helmintos. Este cambio ha sido
observado en áreas con una baja incidencia de CG como la costa pacífica
colombiana. La prevalencia del parasitismo intestinal en esta población es del 50%
con respecto al 25% en las áreas de alto riesgo de las montañas andinas 14.
La oncoproteína cagA está presente en alrededor del 60% de las cepas de H. pylori
en los países occidentales, y en casi todas las del este de Asia. cagA se encuentra
21
Introducción
presente en alrededor del 90% de las cepas de los Andes colombianos (alto riesgo) y
en aproximadamente el 70% de las cepas de la costa pacífica colombiana (bajo
riesgo) 18. Un reciente trabajo identificó que en las áreas de bajo riesgo colombianas
predominan las cepas de H. pylori de origen africano mientras que en las de alto
riesgo son más comunes las de origen europeo 19. Por lo que las cepas africanas
serían menos virulentas que las europeas condicionando el riesgo de CG en
dependencia del genotipo de las mismas.
22
Introducción
23
Introducción
E) Inhibición de la ATPsa H+- K+: El receptor del T4SS de H. pylori, denominado α5β1,
al interactuar directamente con H. pylori activa a ADAM17, la cual induce la
represión transcripcional de la ATPsa H+-K+, responsable de generar el gradiente de
protones H+ que garantiza el contenido ácido del estómago.
24
Introducción
Figura 9. Mecanismos de infección de Helicobacter pylori en el epitelio gástrico. Tomado de 17, 21, 23, 27, 28, 29, 30, 31, 32.
25
Introducción
En la infección por H. pylori existe una respuesta mayoritaria tipo Th1, sin embargo es
incapaz de eliminar la infección, contribuyendo a la inflamación crónica 27. En modelos
murinos y pacientes con una respuesta predominante Th2, existe una disminución del
riesgo de gastritis crónica atrófica multifocal (GCA) y CG. En ratones con respuesta Th2 la
infección con Helicobacter felis induce una leve gastritis a pesar de lograr una alta
colonización bacteriana. Pacientes infectados con una gastritis leve poseen una respuesta
tipo Th0 en la cual secretan tanto citoquinas Th1 como Th2, evidenciando que la
respuesta Th2 ejerce protección en comparación con la respuesta Th133.
5.1. Dieta
26
Introducción
vegetales en la reducción del CG, altos contenidos de los mismos reducen el riesgo de CG
en un 31% 36.
Sin embargo, los consumos elevados de sal, alcohol y carne roja procesada
(especialmente en infectados por H. pylori) 37 probablemente aumentan el riesgo de CG
no cardias. La mayoría de estas asociaciones provienen de estudios caso-control; los
pocos estudios de cohortes han sido más inconsistentes y se requieren más estudios de
cohorte para arribar a conclusiones definitivas. Las carnes procesadas aportan
nitrosaminas y precursores de su formación como los nitratos y nitritos. Las carnes rojas
aportan hierro orgánico que es un precursor de la formación endógena de nitrosaminas,
las cuales son cancerígenas en modelos de experimentación animal.
5.2. Tabaquismo
27
Introducción
ciudad de Cali que poseía un 65.8% de inmigrantes se demostró que la tasa de incidencia
de CG era 5 veces mayor en inmigrantes de los Andes respecto a los de la costa del
Pacífico. La prevalencia de la metaplasia intestinal siguió un comportamiento análogo a la
del CG en dichas poblaciones 14.
El proceso se inicia cuando H. pylori coloniza la mucosa gástrica de la región antral del
estómago, evitando la disminución del pH en el fundus y cuerpo del estómago. Dicha
infección por lo general provoca una gastritis crónica no atrófica que puede durar
28
Introducción
29
Introducción
ciclina que controla la transición G1-S del ciclo celular, y la amplificación de los proto-
oncogenes MET y K-sam contribuyen a la progresión y dispersión metastásica 43.
B) Evitar la muerte celular por apoptosis debido a mutaciones en oncogenes (Ej. MYC) o
alteraciones en los reguladores de la apoptosis, como el programa intrínseco, extrínseco,
BCL2, BAX, BAK, BIM, PUMA, IGF1/2 y citocromo C.
30
Introducción
31
Introducción
Figura 11. Mucosa normal teñida con hematoxilina-eosina. Tomado de Correa et al., 2012 14.
Figura 12. Gastritis crónica no atrófica. Tomado de Correa et al., 2012 14.
32
Introducción
Figura 13. Gastritis crónica atrófica. Tomado de Correa et al., 2012 14.
33
Introducción
Figura 14. Metaplasia intestinal completa. Tomado de Correa et al., 2012 14.
La MII está compuesta por células epiteliales del colón, específicamente células
columnares y caliciformes de tamaño variable secretoras de sialomucinas y sulfomucinas,
ausencia de enterocitos y rara presencia de células de Paneth. Presenta una importante
desestructuración glandular y desaparición parcial o total de las enzimas digestivas
(Figura 15). Existe expresión de mucinas gástricas o neutrales MUC1, MUC5AC, MUC6 y
de la mucina intestinal MUC2. Predomina en pacientes mayores de edad 24.
Figura 15. Metaplasia intestinal incompleta. Tomado de Correa et al., 2007 24.
La conexión causal entre la metaplasia intestinal con el CG fue propuesta por vez primera
en Java y Sumatra en 1938, comparando pacientes malayos con bajas frecuencias de CG,
atrofia y MI, respecto a inmigrantes chinos con altas frecuencias de tales patologías.
34
Introducción
Posteriormente, Jarvi y Lauren reportaron que los CGs de tipo intestinal estaban rodeados
por metaplasia intestinal más frecuentemente que los CGs de tipo difuso 13.
7.4. Displasia
Lesión neoplásica limitada al epitelio, sin invasión de la lámina propia, caracterizada por
atipia celular, esto es, células de formas irregulares, hipertróficas e hipercromáticas, y en
las que se pierde la polaridad apico-basal epitelial; además, hay presencia de frecuentes
mitosis y nucléolos prominentes (Figura 16). Esta des-diferenciación del epitelio es
progresiva y se convierte en carcinoma invasivo cuando la atipia celular alcanza la
membrana basal 14. Su diagnóstico presenta una baja concordancia inter-observador
(k=0.42) 48.
35
Introducción
36
Introducción
37
Introducción
Tabla 2. Principales estudios epidemiológicos de seguimiento de lesiones precursoras gástricas. Tomado y modificado de
de Vries et al., 2007 56 (con excepción de los estudios españoles).
Tamaño Tiempo de Progresión a CG
Diseño de seguimiento después Incidencia anual de CG Incidencia anual de
País de Progresión a CG (%) (% de todos los
estudio (n/100.000 personas años) CG (% de estudios)
muestra de endoscopia (años) estudios)
Gastritis Atrófica
Finlandia Prospectivo 116 15 (máximo) 9 (7.8%) 3.48 ≥517 208.6
23 (máximo) 10 (8.6%) ≥375
Suecia Prospectivo 61 2.7 (%) 0 (0%) 0
China Prospectivo 1240 5 (máximo) Ligera: 1 (0.1%) ≥ 16
Severa: 0 (0%) 0
Italia Prospectivo 42 4 (máximo) 0 (0%) 0
Italia Prospectivo 106 7 (media) 1 (0.9%) ≈135
Metaplasia Intestinal
Italia Prospectivo 261 9 (%) 12 (4.6%) 11.53 ≈ 511 572
Portugal Prospectivo 124 6 (máximo) Tipo I-III: 0 (0%) 0
China Prospectivo 842 5 (máximo) Superficial: 2 (0.8%) ≥ 156
Severa: 16 (2.7%) ≥ 546
UK Prospectivo 93 10 (máximo) 10 (11%) ≥ 1075
España (Soria) 57 Prospectivo 478 12.8 (MII-21%, MIC 19%) Global: 3.77 (2.51–5.67)/1000
MII: 16.5 (10.08-6.86)/1000,
MII vs MIC (HR 11.3 [3.8-33.9]
España (multi- Prospectivo 649 10 MIC: 3.23% Global: 3.09 (2.07-4.6)/1000
céntrico) MII: 6.85% MIC: 2.76/1000,MII: 5.76/1000
MII vs MIC: HR 2.75 (1.06-6.26)
Displasia
Italia Prospectivo 134 1.6 (%) Ligera: 5 (6.1%) 64.56 ≥ 3811 19 876
Moderada: 7 (22.6%) ≥ 14.113
Severa: 6 (46.2%) ≥ 28.846
Italia Prospectivo 125 10 (máximo) Bajo grado: 22(27.2%) ≥ 2716
Alto grado: 36 (81.8%) ≥ 8182
China Prospectivo 546 5 (máximo) Ligera: 14 (2.8%) ≥560
Moderada/Severa: 3 ≥ 1400
(7.0%)
38
Introducción
8. Metaplasia intestinal
39
Introducción
Liu Q propone que un nivel muy alto de CDX2 sería causal de la MIC, la cual es un
estadio completamente diferenciado con una baja tasa de progresión a CG,
mientras que un nivel menor sería responsable de la MII la cual es un estadio de
menor diferenciación y por ende más proclive de avanzar a CG 65.
40
Introducción
41
Introducción
42
Introducción
43
Introducción
Los genes Schlafen (dormir en alemán) constituyen una familia multigénica con un
alto grado de homología de secuencia, similitud funcional y proximidad
cromosomal. Su función es la regulación de la respuesta inmune mediante el
control de la diferenciación de timocitos a células T, así como la activación de las
mismas 82. En ratón esta familia está compuesta por 3 grupos basados, entre otras
características, en el tamaño de sus proteínas 82: grupo I (37-42 KDa, Schlafen1 y 2),
grupo II (58-68 KDa, Schlafen3 y 4) y grupo III (100-104 KDa, Schlafen5, 8, 9, 10, 14).
Todos sus miembros contienen un dominio Slfn de función desconocida, el cual no
ha sido encontrado en otras proteínas. Además, contienen un dominio N-terminal
AAA con una función hipotética en el metabolismo del ácido ribonucleico (ARN) y la
unión a GTP/ATP. Los grupos II y III poseen otro dominio llamado SWADL (Ser-Trp-
Ala-Asp-Leu). La localización subcelular de los miembros de los grupos I y II es el
citoplasma, mientras que los del grupo III lo hacen en el núcleo debido a una señal
de localización nuclear (RKRRR) 82.
Los genes de los grupos I y II fueron los primeros descubiertos mediante técnicas de
hibridación sustractiva para conocer genes reguladores de la diferenciación de
timocitos a células T 82. Posteriormente, usando técnicas similares aplicadas a la
infección por Lysteria monocitogenes fue descubierto el grupo III. En humanos
existen cinco isoformas (Slfn5, 11, 12, 13 y 14), a diferencia de en el ratón solo un
miembro tiene localización citoplasmática, Slfn12, el cual carece del dominio
helicasa y pertenece al grupo II, los restantes pertenecen al grupo III. La mayoría de
los Schlafen murinos no presentan homólogos en humanos y viceversa 82.
44
Introducción
45
Introducción
9. Susceptibilidad genética
Además de los factores de virulencia de H. pylori y los componentes ambientales, el
largo proceso inflamatorio de respuesta causado por la infección bacteriana es uno
de los factores desencadenantes del proceso neoplásico. La susceptibilidad genética
del huésped tiene un papel importante en la carcinogénesis gástrica como lo
evidencia el hecho de que la heredabilidad del CG se ha calculado en un 28% y que
la historia familiar es un factor de riesgo independiente de la infección por
H. pylori 90. Aunque el efecto de dicha infección es el principal factor de riesgo de
CG, los factores heredados que confieren susceptibilidad genética a la infección y a
la severidad de la reacción inflamatoria y otros mecanismos de tumorigénesis son
también una parte de este complejo proceso.
46
Introducción
47
Introducción
ambientales a que están expuestas dichas poblaciones 94. Con el objetivo de ganar
potencia estadística para generar resultados más confiables se han realizado meta-
análisis los cuales tienen en cuenta varios estudios incrementando el tamaño
muestral y por ende la potencia estadística para detectar una posible asociación 95.
48
Introducción
49
Introducción
Los estudios de asociación de genoma completo (GWAS, del inglés genome wide
association study) representan un avance sobre los estudios de genes candidatos
debido a que el genoma es inspeccionado en su totalidad mediante el genotipado
de cientos de miles de tagSNPs usando chips de ADN y poseen tamaños de muestra
en el orden de miles de casos y de controles que permiten tener un poder
estadístico adecuado para detectar el pequeño efecto de los SNPs en la causalidad
de las enfermedades complejas 117. En poblaciones asiáticas se han realizado
diversos estudios de asociación con marcadores de todo el genoma asociando
nuevos genes con el CG cardias y no cardias. En la tabla 3 se indican los SNPs
asociados junto con otras características de los estudios.
50
Introducción
51
Introducción
Variantes asociadas con las LPGs en genes de esta vía incluyen a CD14 260C>T 130,
TLR4 +896A/G 131, PTPN11 (rs2301756, rs12229892) 132, entre otros. El mecanismo
funcional por el que los SNPs de esta vía incrementan el riesgo de carcinogénesis
gástrica está relacionado con el proceso inflamatorio. De este modo TLR4 +896A/G
fue asociado con la hipoclorhidria, un fenotipo que es consecuencia de la GCA 131.
53
Introducción
Los alelos S de VNTR de MUC1 se asociaron con un mayor riesgo de gastritis atrófica
y metaplasia intestinal en población portuguesa, con una de las mayores incidencias
de CG en Europa 137. En la población colombiana los genotipos SS y LS de un
polimorfismo VNTR de MUC1 se asocian con la MI 138 mientras que en portugueses
los genotipos LL y SS son asociados con la metaplasia intestinal completa e
incompleta, respectivamente 139. Concluyendo que se observa una tendencia a la
asociación entre el alelo S de MUC1 con la metaplasia intestinal.
54
Introducción
55
Introducción
130 Asociación con lesiones -CD14 con MI, 1.45 (1.11-1.89) -CD14 260T Controles afectados por
Venezuela, CD14 C260T
precursoras: asociado con GCA
2007 (rs2569190)
expresión
Casos: GCA, MI y DISP -GCA, 289 incrementada de
TLR4 Asp299Gly -MI, 543 CD14
(rs4986790) -DISP, 118
-Asp299Gly
asociado con
NOD2 Controles: -Normales y respuesta reducida
del3020ensC, -Normales, gastritis superficial y gastritis al LPS
Gly908Arg gastritis crónica superficial, 87 -Gly908Arg
-Gastritis crónica, es un SNP no
(rs2066847)
996 sinónimo
NOD2
Gly908Arg
(rs2066845)
57
Introducción
140 Asociación con la evolución de las -SNPs en extremo 3’ de MUC2, (rs10794293, Ausencia de controles
España, 22 tagSNPs en
lesiones: rs3924453 y rs4077759) asociados con menor sanos
2010 MUC1, MUC2 y
riesgo de progresión
MUC6 Casos Casos
-GNA, GCA, MIC, MII, DISP que Progresión: 109 -SNPs en extremo 5’ de MUC2, (rs10902073,
progresan a LPGs de mayor rs10794281,
severidad. Controles rs2071174 and rs7944723) asociados con mayor
-GNA, GCA, MIC, MII, DISP que Controles riesgo de regresión
regresan a lesiones de menor Regresión:111
severidad o se mantienen Estables: 167
estables
59
HIPÓTESIS
Hipótesis
63
OBJETIVOS
Objetivos
67
MATERIALES Y MÉTODOS
Materiales y Métodos
Este estudio se basa en uno previo tipo caso control anidado en la cohorte de la
European Prospective Investigation into Cancer and Nutrition (EPIC), denominado
Eurgast. La cohorte EPIC se encuentra constituida por 521 457 sujetos sanos de
ambos sexos, y edades entre 35-70 años, reclutados entre los años 1992 y 2000 en
10 países europeos (incluyendo España), y que durante el seguimiento (media >15
años) han desarrollado diversas enfermedades 143. Los casos identificados de CG
(N=365) a la fecha del análisis, son adenocarcinomas definidos por el código C16 de
la 10ma revisión de la Clasificación Internacional de Enfermedades.
71
Materiales y Métodos
72
Materiales y Métodos
73
Materiales y Métodos
Se realizó una búsqueda de pacientes que entre 1995 y 2004 tenían un diagnóstico
histológico de GCA, MI o DISP. Del total de pacientes identificados fueron excluidos
12 073 pacientes por presentar algunos de los siguientes criterios de exclusión:
74
Materiales y Métodos
Se seleccionaron 1 538 pacientes para ser incluidos en el estudio que cumplían con
los siguientes criterios:
Los mismos fueron invitados a ser sometidos a una gastroscopia con toma de
muestras gástricas e independientemente de saliva, para la extracción del ADN. Así
como, a responder un cuestionario personal acerca de historia de tabaquismo, uso
de anti-inflamatorios no esteroideos, diagnóstico y tratamiento de infección por
Helicobacter pylori, historia familiar de CG, grupo sanguíneo, ovariotomía y uso de
hormonas anticonceptivas. Dicha información epidemiológica fue obtenida de las
historias clínicas o a partir de familiares en el caso de aquellos pacientes con un
punto final predefinido al seguimiento y los que no pudieron ser entrevistados.
Todos los participantes firmaron un consentimiento informado aprobado por el
comité ético de cada hospital participante. Para el presente estudio de asociación
genética el tamaño de muestra fue de 559 pacientes afectados con lesiones
precursoras gástricas.
2.2. Histopatología
75
Materiales y Métodos
Así, fueron colectadas muestras de mucosa bucal con hisopos bucales Dri‐Capsules
(SGC‐50, Isohelix, Inglaterra) (N=20) y BuccalAmp™ DNA Extraction Kit (Epicentre,
EUA) (N=20). También, muestras de saliva mediante el kit Oragene DNA (OG-500,
DNA Genotek, Canadá) (N=10) siguiendo las recomendaciones de los fabricantes.
Todas se incubaron a temperatura ambiente (TA) durante 30 días y,
posteriormente, se extrajo el ADN según el protocolo de los fabricantes.
76
Materiales y Métodos
77
Materiales y Métodos
Teniendo en cuenta que los mejores resultados de la anterior prueba piloto fueron
obtenidos con el ADN proveniente de saliva y que los kits eran de un coste elevado
(13.9 euros/muestra, Abyntek distribuidor en Europa de DNA Genotek) para
nuestro presupuesto, se planteó la posibilidad de confeccionar en el laboratorio un
similar de la solución estabilizadora de saliva que usan los kits Oragene DNA. Se
realizó una prueba piloto con las diferentes soluciones estabilizadoras de saliva
propuestas en la aplicación de patente asociada a dichos kits: “Compositions and
methods for obtaining nucleic acids from sputum. Appl. No 12/338.848” 159. Esta
prueba fue necesaria debido a que el lenguaje ambiguo de las patentes no permite
conocer cuál es la solución elegida por los autores como la de mejores resultados
(Tabla 5).
78
Materiales y Métodos
79
Materiales y Métodos
Se extrajo ADN del tejido gástrico de las biopsias del reclutamiento con el objetivo
de genotipar cagA y vacA de acuerdo al siguiente protocolo, basado en el método
del fenol: cloroformo: isoamilalcohol 158.
80
Materiales y Métodos
81
Materiales y Métodos
82
Materiales y Métodos
Tabla 6. Genes candidatos seleccionados a partir de estudios de modelos murinos de carcinogénesis gástrica.
84
Materiales y Métodos
85
Materiales y Métodos
Haciendo uso del navegador del sitio web del proyecto HapMap
(http://hapmap.ncbi.nlm.nih.gov/) en población de origen europeo (CEU) release
#28 (Phases 1, 2 & 3 - merged genotypes & frequencies, db SNP b 126) fueron
descargados los datos genotípicos de cada gen candidato seleccionado; así como 5 y
3 kilobases adicionales de sus regiones 5’ y 3’. Posteriormente, dichos datos
genotípicos fueron cargados en el programa informático Haploview versión 4.2 y
seleccionados los SNPs con MAF≥0.05. Para la definición de bloques en LD
(haplotipos) se aplicó el método de intervalos de confianza de Gabriel et al. 92,
utilizando los parámetros por defecto en Haploview v.4.2. Los tagSNPs fueron
seleccionados utilizando el método Aggressive Tagger con un umbral de r2≥0.8
entre cada par de tagSNPs y SNPs etiquetado; así como, haplotipos etiquetados con
MAF≥0.05 93. Para los polimorfismos significativos extraídos de GWAS y de estudios
de asociación con la evolución de las LPGs de Soria 140 no se realizó selección de
tagSNPs. Los análisis funcionales in silico de los SNPs genotipados fueron realizados
como se describió en el apartado 1.5 de esta sección. La cobertura de la variabilidad
de tagSNPs en las regiónes genéticas fue obtenida por Haploview eliminando los
polimorfismos que no eran tagSNPs del total de los existentes en los bloques de
haplotipos (Tabla 9).
86
Materiales y Métodos
87
Materiales y Métodos
Variant Effect
Selección a Gen SNP b Cromosoma c Posición d Alelos e MAF f PupaSuite h HWE i Cobertura j
Predictor g
A PTGS2 rs5277 1 184914820 G/C 0,192 Sinónimo, NMD Región conservada 0,373
A PTGS2 rs2745557 1 184915844 G/A 0,204 Intrónico, NMD, Región conservada 0,475
A PTGS2 rs689466 1 184917374 A/G 0,175 Región
RR 5’, RR Región conservada 0,201
A PTGS2 rs12042763 1 184918499 G/T 0,273 Región 5’, RR Triple hélice, 0,295
conservada
A RUNX3 rs6663310 1 25094966 T/C 0,334 Región 3’ ------------- 0,145 91 %
A RUNX3 rs9438858 1 25096492 C/A 0,391 Región 5’ Región conservada 0,098
A RUNX3 rs7553295 1 25096597 G/T 0,238 Región 3’ Región conservada 0,799
A RUNX3 rs9438859 1 25096908 T/C 0,232 Región 3’, RR Conservada, RR 0,297
A RUNX3 rs7514332 1 25098262 G/C 0,054 Region 3’ Región conservada 0,384
A RUNX3 rs2236850 1 25112928 T/C 0,408 Intrónico, RR Región conservada 0,503
A RUNX3 rs2282718 1 25113643 G/A 0,335 Intrónico, RR Región conservada 1,000
A RUNX3 rs9438876 1 25113703 A/G 0,489 Intrónico, RR Región conservada 0,116
A RUNX3 rs2236852 1 25116354 A/G 0,486 Intrónico ------------- 0,579
A RUNX3 rs7517302 1 25126904 T/C 0,362 Intrónico, RR Región conservada 0,616
A TFF1 rs225352 21 42650501 C/T 0,086 Region 3’ ------------- 1,000 100 %
A TFF1 rs178740 21 42650741 T/C 0,198 Region 3’ ------------- 0,012
A TFF1 rs1547374 21 42651964 A/G 0,312 Region 3’ ------------- 0,516
A TFF1 rs4919984 21 42653915 A/G 0,222 Region 3’ Región conservada 1,000
A TFF1 rs225355 21 42657551 T/C 0,492 Intrónico ------------- 0,114
A TFF1 rs225357 21 42659382 C/T 0,431 Intrónico Región conservada 0,638
A TFF1 rs4920094 21 42659402 G/A 0,428 Intrónico Región conservada 0,158
A TFF1 rs225358 21 42659581 T/C 0,303 Intrónico, Región conservada 0,319
región de
empalme
A TFF1 rs2156310 21 42659675 G/A 0,151 alternativo
UTR 5’ TFBS/Triple hélice/ 1,000
conservada
A TFF1 rs225359 21 42660505 A/G 0,327 Región 5’ ------------- 0,833
89
Materiales y Métodos
Variant Effect
Selección a Gen SNP b Cromosoma c Posición d Alelos e MAF f PupaSuite h HWE i Cobertura j
Predictor g
A TFF1 rs1788413 21 42661759 A/C 0,162 Region 3’ ------------- 0,864
A TFF1 rs9976977 21 42662543 G/A 0,369 Region 3’ ------------- 0,112
A TFF1 rs424694 21 42662841 C/T 0,413 Region 3’ ------------- 0,445
A TFF1 rs13047838 21 42665938 C/T 0,257 UTR 3’ Conservada, Triple 0,183
hélice
D TFF2 rs1079380 21 42643342 A/G 0,489 Intrónico Región conservada 0,309 No es tagSNP
Variant Effect
Selección a Gen SNP b Cromosoma c Posición d Alelos e MAF f PupaSuite h HWE i Cobertura j
Predictor g
B NFKBIA rs3138056 14 34938265 G/A 0,319 Región 3’, RR Región conservada 0,670 100 %
B NFKBIA rs4982269 14 34939618 T/C 0,338 Región 3’ Región conservada 0,536
B NFKBIA rs696 14 34940844 G/A 0,372 UTR 3’, NMD Región conservada 0,321
B NFKBIA rs3138054 14 34942058 G/A 0,154 Intrónico, RR Región conservada 0,859
B NFKBIA rs3138053 14 34944605 A/G 0,284 Región 3’ Región conservada 0,820
B NOD1 rs11536450 7 30434089 C/G 0,13 Intrónico, NMD ------------- 1,000 100 %
B NOD1 rs2907749 7 30452266 A/G 0,295 Intrónico, RR ------------- 0,074
B NOD1 rs7789045 7 30460547 T/A 0,463 Intrónico, NMD Región conservada 0,641
B NOD1 rs2906766 7 30466100 T/C 0,302 UTR 5’, RR Región conservada 0,661
B NOD1 rs3823773 7 30466486 A/G 0,087 Intrónico, RR ------------- 0,768
B NOD1 rs4272257 7 30470374 C/T 0,157 Intrónico, NMD, ------------- 1,000
B NOD1 rs2709803 7 30477059 C/T 0,217 Intrónico,
RR NMD Región conservada 1,000
B NOD1 rs4720004 7 30478744 T/C 0,144 Intrónico, NMD Región conservada 0,706
B NOD1 rs2256023 7 30481516 T/C 0,426 Intrónico, NMD Región conservada 0,924
B NOD1 rs17770244 7 30482337 A/G 0,061 Intrónico, NMD Región conservada 1,000
B PTPN1 rs17822304 12 111373075 C/A 0,049 Intrónico Región conservada 1,000 100 %
B 1
PTPN1 rs11614544 12 111389338 A/G 0,154 Intrónico, RR ------------- 0,285
B 1
PTPN1 rs11066320 12 111390798 G/A 0,437 Intrónico ------------- 0,778
B 1
PTPN1 rs11066322 12 111406912 G/A 0,189 Intrónico ------------- 0,128
B 1
PTPN1 rs11066323 12 111407744 G/A 0,103 Intrónico ------------- 0,613
B 1
PTPN1 rs7958372 12 111419936 C/T 0,085 Intrónico ------------- 0,561
B 1
PTPN1 rs7953150 12 111425954 G/A 0,286 Intrónico, RR ------------- 0,821
B 1
SRC rs6017901 20 35401507 T/C 0,043 Región 3’ ------------- 0,583 88 %
B SRC rs6094373 20 35402488 T/C 0,221 Región 3’ Región conservada 0,418
B SRC rs6017916 20 35404692 A/C 0,251 Región 3’, RR Región conservada 0,712
B SRC rs12106024 20 35412762 C/G 0,203 Intrónico, RR Región conservada 0,115
B SRC rs1547836 20 35414365 G/C 0,188 Intrónico, RR Región conservada 0,878
B SRC rs16986606 20 35417842 A/G 0,146 Intrónico, RR Región conservada 0,263
91
Materiales y Métodos
Variant Effect
Selección a Gen SNP b Cromosoma c Posición d Alelos e MAF f PupaSuite h HWE i Cobertura j
Predictor g
B SRC rs6017996 20 35420067 G/A 0,144 Región 3’, RR ------------- 0,851 88%
B SRC rs6018027 20 35424206 T/C 0,275 Intrónico ------------- 0,488
B SRC rs6063022 20 35426741 C/T 0,16 Intrónico ------------- 1,000
B SRC rs6018088 20 35432039 C/A 0,136 Intrónico ------------- 0,692
B SRC rs6018148 20 35439194 C/T 0,204 Intrónico ------------- 0,670
B SRC rs7269342 20 35441575 T/C 0,146 Intrónico, RR ------------- 0,710
B SRC rs3790150 20 35448513 A/G 0,163 Intrónico Región conservada 0,395
B SRC rs754626 20 35450754 A/C 0,223 Intrónico ------------- 0,503
B SRC rs754625 20 35450938 T/C 0,288 Intrónico ------------- 0,821
B SRC rs6018257 20 35455953 T/C 0,122 Intrónico, RR Región conservada 0,385
B SRC rs1570209 20 35462245 T/C 0,102 Intrónico ------------- 1,000
B SRC rs17785475 20 35466307 C/T 0,067 UTR 3’ Región conservada 0,711
C LRFN2 rs2494938 6 40644106 A/G 0,462 Intrónico ------------- 0,263 No es tagSNP
C MUC1 rs2070803 1 153424339 T/C 0,43 Región 3’ Región conservada 0,089 No es tagSNP
C PSCA rs2294008 8 143758933 C/T 0,437 Intrónico, UTR Región conservada 0,349 No es tagSNP
C PRKAA rs13361707 5 40827641 C/T 0,268 Intrónico,
5’ RR Región conservada 0,156 No es tagSNP
C 1
DNAH1 rs2285947 7 21550613 G/A 0,488 Intrónico Región conservada 0,194 No es tagSNP
C 1
ZBTB20 rs9841504 3 115845454 C/G 0,059 Intrónico, RR Región conservada 0,209 No es tagSNP
B CD14 rs11167532 5 139968081 G/A 0,45 Región 5’ ------------- 0,075 88 %
B CD14 rs778588 5 139987195 T/C 0,294 Región 3’ ------------- 0,181
B CD14 rs2569190 5 139993100 A/G 0,494 Intrónico, UTR Conservado/Triple 0,459
5’, RR hélice
B CD14 rs5744455 5 139993491 C/T 0,2 Región 5’, RR Región conservada 0,384
B CD14 rs1583005 5 140011721 C/T 0,418 Intrónico, NMD ------------- 0,390
B CDH1
(IK) rs16260 16 67328535 C/A 0,275 Región 5’ Región conservada 0,907 60 %
B CDH1 rs1125557 16 67367100 A/G 0,433 Intrónico, NMD Región conservada 0,345
B CDH1 rs12597188 16 67372327 G/A 0,317 Intrónico, NMD ------------- 1,000
92
Materiales y Métodos
Variant Effect
Selección a Gen SNP b Cromosoma c Posición d Alelos e MAF f PupaSuite h HWE i Cobertura j
Predictor g
B CDH1 rs9929218 16 67378447 G/A 0,292 Intrónico, NMD, ------------- 0,576
B CDH1 rs7186053 16 67396794 G/A 0,423 Intrónico,
RR NMD ------------- 0,705
B CDH1 rs4783573 16 67398089 A/G 0,356 Intrónico, NMD ------------- 0,419
D IL1RN rs2637988 2 113593250 A/G 0,389 Intrónico, ------------- 0,041 No es tagSNP
Región 3’
D MUC2 rs10902073 11 1050934 C/A 0,336 Intergénico, RR ------------- 0,097 34 %
D MUC2 rs10794281 11 1053149 T/C 0,397 Intergénico Región conservada 0,176
D MUC2 rs2071174 11 1063712 T/C 0,29 Región 5’, RR ------------- 1,000
D MUC2 rs7944723 11 1083710 C/G 0,218 Exónico ------------- 0,105
D MUC2 rs10794293 11 1088939 C/T 0,347 Intrónico, ------------- 0,358
D MUC2 rs3924453 11 1095806 G/A 0,263 Región 3’,
3’ ------------- 0,341
D MUC2 rs4077759 11 1095976 T/C 0,353 Región 3’,
intergénico ------------- 0,361
D MUC6 rs6597947 11 1027029 G/T 0,079 Región 3’
intergénico Región conservada 0,005 No es tagSNP
a
Se refiere a la procedencia de los genes candidatos. A: modelos murinos de carcinogénesis gástrica, B: Mecanismo de infección por H. pylori, C: Genes
asociados a GWAS en CG, D: Genes asociados a la evolución de LPGs en estudio de seguimiento de Soria 140. b Los SNPs se ordenan en sentido 5’-3’. c-f
Información de HapMap. e Alelo común/alelo variante. f MAF: Frecuencia del alelo menor en la muestra total estudiada. g h Efecto funcional predicho de
programas Variant Effect Predictor y Pupasuite respectivamente. NMD significa que el SNP se encuentra localizado en transcrito sujeto a mecanismo de
degradación de ARNm Nonsense mediated mRNA decay. Región conservada se refiere a una secuencia de ADN conservada evolutivamente entre ratón y
humano. TFBS se refiere a sitios de unión de factores de transcripción. Triple hélice se refiere a la formación de esa estructura en el ADN como
consecuencia del alelo variante. Región regulatoria se refiere a aquella que afecta la expresión génica. ----- Significa que la herramienta no predijo dicho
SNP. i p-valor de equilibrio de Hardy-Weinberg. j Cobertura de variabilidad de tagSNPs en la región genética.
93
Materiales y Métodos
94
Materiales y Métodos
G4 o MII: La MII afecta la casi totalidad de los fragmentos, aunque cabe la presencia
de áreas de MIC. Equivalente al estadio III de la clasificación de Filipe et al. 50.
Después de confirmar con este estudio piloto que las muestras eran válidas para su
análisis mediante microarrays de expresión, se procedió al análisis de los genes y
mecanismos funcionales desregulados en la progresión de MI a CG en la totalidad
de los casos a estudiar (Tabla 11).
95
Materiales y Métodos
Tabla 11. Muestras usadas en el estudio de expresión diferencial por microarrays de los subtipos histológicos de metaplasia intestinal
que progresan a cáncer gástrico, respecto a los que no progresan y a la mucosa sana.
Diagnóstico Extensión de Diagnóstico al Infección
Localización Años de
Código al MI al final del Estatus d por Sexo e Edad f Procedencia g
anatómica seguimiento c
reclutamiento a reclutamiento b seguimiento H. pylori
A7 G4 E4 Antro CG intestinal 10 MII-CG NO M 52 H. Soria
A8 G4 E4 Antro CG intestinal
stinalntestinaI 8.4 MII-CG SI H 69 H. Soria
A9 G4 E3 Antro CG Difuso
ntestinal 7.7 MII-CG NO H 45 H.Soria
A10 G4 E3 Antro CG intestinal 2 MII-CG NO H 73 H.Soria
A11 G4 E3 Antro CG intestinal 1 MII-CG NO H 60 H.Soria
A12 G4 E3 Antro CG intestinal 11.8 MII-CG SI M 66 H.Soria
A13 G3 E4 Antro CG 13.4 MII-No CG H 53 H.G Marañón
A14 G4 E4 Antro MII 15 MII-No CG SI H 64 H.Soria
A15 G3-G4 E3-E4 Antro MII 11 MII-No CG NO M 64 H.Soria
A16 G3-G4 E4 Antro MII ND MII-No CG NO H 64 H.Soria
A17 G3-G4 E3-E4 Antro MII 14 MII-No CG SI H 65 H.Soria
A18 G3-G4 E4 Antro MII ND MII-No CG NO M 74 H.Soria
A19 G3-G4 E4 Antro MII ND MII-No CG NO H 71 H.Soria
A20 G1 E4 Antro CG no cardias 2 MIC-CG ND H 57 Multicéntrico
A21 G1 E4 Cuerpo CG intestinal 8 MIC-CG SI H 57 Multicéntrico
A22 G1 E4 Antro CG intestinal 1 MIC-CG ND H 44 Multicéntrico
A23 G1 E4 Antro MIC 15.37 MIC-CG NO H 56 H.G Marañón
A24 G2 E4 Antro MIC 14 MIC-CG SI M 59.7 Multicéntrico
A25 G1 E4 Antro MIC 8.24 MIC-CG NO M 57
4 H.G
LPGsMarañón
A26 G1 E4 Antro MIC 5.95 MIC-CG NO M 78 H.G Marañón
A27 G1 E4 Antro CG no cardias 6.2 MIC-CG NO H 76 H.G Marañón
A28 G1 E4 Transiciona MIC 10.6 MIC-No CG NO H 67 H.G Marañón
A29 G2 E4 Antro
l MIC 10.5 MIC-No CG SI H 63 H.G Marañón
A30 G2 E4 Transiciona MIC 12.31 MIC-No CG NO H 74 H.G Marañón
A31 G2 E4 Antro
l MIC 13 MIC-No CG NO H 56.9 Multicéntrico
8 LPGs
96
Materiales y Métodos
La extracción de ARN total se realizó utilizando el kit Recover All Total Nucleic Acid
(AM1975, Ambion, EUA) de acuerdo al protocolo recomendado por los fabricantes,
con modificaciones menores:
98
Materiales y Métodos
99
Materiales y Métodos
Se realizó una prueba piloto con el objetivo de comparar el microarray Human Gene
2.0 ST (Affymetrix, EUA) respecto al Almac Xcel (Affymetrix, EUA), el cual es
específico de muestras parafinadas. El Human Gene 2.0 ST es un array en el que las
sondas hibridan a lo largo de todo el transcrito interrogando 40 716 de los mismos,
mediante 21 sondas/gen. En el Almac-Xcel las sondas hibridan en el extremo 3’ del
ARNm, interrogando alrededor de 97 000 transcritos y es equivalente en un 70% al
Human Genome U133 Plus 2.0 array, el más usado de Affymetrix 176. Previo a la
hibridación en el chip se amplificó el ADNc usando el WT-Ovation™ FFPE System V2
(Nugen, EUA), específico para ARN degradado de muestras parafinadas. En esta
prueba piloto se emplearon las muestras descritas en la tabla 10 y posteriormente
se procedió a hibridar muestras clínicas (Tabla 11), con el objetivo de conocer
procesos fisiológicos causales involucrados en la progresión desde MI a CG y que
caracterizan a la MI.
100
Materiales y Métodos
Los genes diferencialmente expresados entre los grupos MII-CG vs MII-No CG,
MIC-CG vs MIC-No CG, MI-No CG vs Sanos se obtuvieron aplicando la prueba t
moderada, específica para datos de microarrays (Linear Models for Microarray
Data, Limma), que tiene en cuenta la media, desviación estándar y N de los grupos a
comparar 179, y se ajustó para las comparaciones múltiples por la tasa de falsos
positivos (FDR) 175. Se consideraron las diferencias estadísticamente significativas si
p<0.05 y el Fold Change≥2 (genes sobre-expresados) ó ≤0.5 (genes sub-expresados).
Los análisis se realizaron en el programa estadístico R (v3.1.1) 180 con sus paquetes
estándar y, adicionalmente, los paquetes Biobase y Limma de Bioconductor 181. Se
comprobó si los genes diferencialmente expresados en las comparaciones MII-CG vs
MII-No CG y MIC-CG vs MIC-No CG también estaban diferencialmente expresados
en MI-No CG vs Sanos.
101
Materiales y Métodos
en cada grupo. Por otra parte se construyeron redes moleculares de alta puntuación
(score>15) que incluyen genes diferencialmente expresados interconectados con
otras moléculas. Un score ≥3 significa que es significativa (IC 99.9%) la probabilidad
de que dichos genes interactúen entre sí.
El programa bioinformático Gene Set Enrichment Analysis (GSEA) 183 identifica si los
transcritos de un experimento de microarrays están significativamente enriquecidos
en vías de señalización celular u otros procesos biológicos mediante la comparación
con conjuntos génicos diferencialmente expresados y representativos de los
anteriores, previamente depositados en la base de datos Molecular Signatures
Database (MSigDB). Para determinar los conjuntos génicos significativamente
enriquecidos calcula el Enrichment Score (ES) definido como la máxima desviación
respecto a 0 encontrada entre los conjuntos génicos de bases de datos y los
resultadosdel experimento, y su p-valor nominal asociado no ajustado, así como
ajustado para comparaciones múltiples mediante FDR.
102
Materiales y Métodos
En el análisis del GSEA no todos los componentes del conjunto génico contribuyen
de igual forma a la señal de enriquecimiento sino que existe un subconjunto que lo
hace de manera más significativa, denominados leading edge genes 185, siendo el
núcleo que contribuye a la señal de enriquecimiento ES. Teniendo esto en cuenta,
para las comparaciones con conjuntos génicos significativos MII-CG vs MII-No y
MI-No CG vs Sanos se determinaron los genes leading edge de los conjuntos génicos
significativos y se compararon con los genes sobre-expresados de dichas
comparaciones.
103
Materiales y Métodos
104
Materiales y Métodos
Tabla 12. Genes significativos del microarray de las comparaciones MII-CG vs MII-No CG, MIC-CG vs MIC-No CG y MI-No CG vs Sanos a ser
validados por qRT-PCR.
Sonda Tamaño
Gen Transcrito(s) Comparación Secuencia de los cebadores
UPL amplicon (pb)
ACTB NM_001101.3 Genes de referencia 64 97 ccaaccgcgagaagatga y ccagaggcgtacagggatag
B2M NM_004048.2 Genes de referencia 42 86 ttctggcctggaggctatc y tcaggaaatttgactttccattc
G6PD NM_000402.4, NM_001042351.2 Genes de referencia 22 91 ctgcagatgctgtgtctggt y tgcatttcaacaccttgacc
GAPDH NM_002046.3 Genes de referencia 24 70 cctcctcctaagatggtgtctg y gaccgatgcgtccaaatc
RPL29 NM_000992.2 Genes de referencia 53 81 ccagaggcgtacagggatag
caggctcccaaacgtacc y gcaccagtccttctgtcctc
AFAP1-AS1 NR_026892.1 MII-CG vs MII-No CG 14 60 tcaccgatttcaaagcttcc y ttcaggcgtgagtgattgtt
ATP5A1 NM_001001937.1, NM_004046.4 MII-CG vs MII-No CG 61 93 tgctattggtcaaaagagatcca y gtagccgacaccacaatgg
ATP6V0E1 NM_003945.3 MII-CG vs MII-No CG 45 112 cctcactgtgcctctcattg y caccaacatggtaatgataactcc
GNAS NM_000516.4,NM_001077488.2, MII-CG vs MII-No CG 3 72 aaggacaagcaggtctaccg y ggtgcttttaccagattctcca
NM_001077489.2, NM_001077490.1,
NM_016592.2,NM_080425.2, NM_080426.2
GSTA1 NM_145740.3 MII-CG vs MII-No CG 53 90 acggtgacagcgtttaacaa y ccgtgcattgaagtagtgga
HLA-C NM_001243042.1, NM_002117.5 MII-CG vs MII-No CG 32 116 acccggactcacattctcc y gcggtgtcgaaatacctcat
HLA-DQA1 NM_002122.3 MII-CG vs MII-No CG 68 101 accaagggccattgtgaat y aatcgggccagagaatagtg
HLA-DRB3 NM_022555.3 MII-CG vs MII-No CG, 41 113 gggctgttcatctacttcagga y caaagctggggcagaagat
MIC-No CG vs MIC-CG
HLA-DRB4 NM_021983.4 MII-CG vs MII-No CG, 41 98 gggacagggctgttcatcta y ccttgaatgtggtcatctgc
MIC-No CG vs MIC-CG
IGHG4 ENST00000390543 MII-CG vs MII-No CG 79 106 tcctccatcgagaaaaccat y ggctgacctggttcttggt
IK NM_006083.3 MII-CG vs MII-No CG 11 103 agctgacccagatcctttca y tcagcctcaggaggtttctt
IL1R2 NM_173343.1, NM_004633.3 MII-CG vs MII-No CG, 72 89 cacatagagagcgcctaccc y
MIC-No CG vs MIC-CG ggcacttcaatgtagttctcattatt
MGAM NM_004668.2 MII-CG vs MII-No CG 3 93 ccaggtgacatgggacataga y gcagtgacatttctggcattac
MME NM_000902.3, NM_007287.2, MII-CG vs MII- No CG 67 78 ggggaggctttatgtggaag y tctcggatctgtgcaatcaa
NM_007288.2, NM_007289.2
PGC NM_001166424.1, NM_002630.3 MII-CG vs MII- No CG 21 61 cttgtgggtgccctctgt y gttgaagcgggagtgactg
PI3 NM_002638.3 MII-CG vs MII- No CG 34 95 tgatcgtggtggtgttcct y acggcctttgacagtgtctt
PPIA NM_021130.3 MII-CG vs MII- No CG 48 97 atgctggacccaacacaaat y tctttcactttgccaaacacc
105
Materiales y Métodos
Sonda Tamaño
Gen Transcrito(s) Comparación Secuencia de los cebadores
UPL amplicon (pb)
RHOA NM_001664.2 MII-CG vs MII- No CG 8 83 gggagctagccaagatgaag y gtacccaaaagcgccaatc
TMEM25 NM_001144034.1, NM_001144035.1, MII-CG vs MII- No CG 79 109 cgcctctgtcatccttaatgt y caccagggcaaacaggac
NM_001144036.1, NM_001144037.1,
NM_001144038.1, NM_032780.3
ANPEP NM_001150.2 MI-No CG vs Sanos 18 75 catccatcagagatggcagac y tgctgaagagatcgttctgg
APOB NM_000384.2 MI-No CG vs Sanos 55 78 gacgacttttctaaatggaacttctac y
CDH17 NM_001144663.1, MI-No CG vs Sanos 2 72 cgtcaccagtaaccgtatttga y gtgcagtaagggtcccgata
ctcagttttgaatatggtgagttttt
CDX1 NM_001804.2 MI-No CG vs Sanos 70 76 acgccctacgagtggatg y tgtccttggtccgagtcttac
NM_004063.3
NM_001265.3 MI-No CG vs Sanos 34 82
CDX2 atcaccatccggaggaaag y tgcggttctgaaaccagatt
CLCA1 NM_001285.3 MI-No CG vs Sanos 62 78 ccaatggaagaatacaagcagtaa y
CPS1 NM_001122633.2, NM_001122634.2, MI-No CG vs Sanos 8 64 caagttttgcagtggaatcg
tgcctccctgacacttcttta y actgggtagccaatggtgtc
NM_001875.4
DMBT1 NM_004406.2,NM_007329.2, MI-No CG vs Sanos 21 77 accaacttaccggcattgac y cgacacctgtcacctccatt
FABP1 NM_017579.2
NM_001443.2 MI-No CG vs Sanos 78 90 tgatccaaaacgaattcacg y caccttccaactgaaccactg
MUC12 NM_001164462.1 MI-No CG vs Sanos 72 68 cctggaaaccttagcaccag y gacagacgcattgttttccat
MUC17 NM_001040105.1 MI-No CG vs Sanos 17 76 ggggtgaacatcacaaagcta y
MUC3A ENST00000319509, ENST00000422757, MI-No CG vs Sanos 18 102 gtggagatcctgtccctgag y cacctgctcatactcgctctc
gtgtgtacttggttcttaggaggac
ENST00000414964
OLFM4 NM_006418.4 MIC-CG vs MIC-NO CG, 24 74 atcaaaacacccctgtcgtc y gctgatgttcaccacaccac
MI-No CG vs Sanos
PRDM10 NM_020228.2, NM_199437.1, MI-No CG vs Sanos 36 88 gtgaggacacggatctgga y gcgttattgcactcctcaca
NM_199438.1, NM_199439.1
SI NM_001041.3 MI-No CG vs Sanos 86 65 ttttggcatccagattcgac y atccaggcagccaagaatc
SLC26A3 NM_000111.2 MI-No CG vs Sanos 2 74 ccatcatcgtgctgattgtc y agctgccaggacggactt
TFF3 NM_003226.3 MI-No CG vs Sanos 4 87 gctgctgctttgactccag y tggaggtgcctcagaaggt
CYP3A4 NM_001202855.2, NM_017460.5 MI-No CG vs Sanos 2 96 gatggctctcatcccagactt y agtccatgtgaatgggttcc
SLFN5 NM_144975.3 IHQ 84 86 agcaagcctgtgtgcattc y accactctgtctgaaaatactgga
SLFN12 NM_018042.3 IHQ 2 76 aacctggacactcttcactggt y atgcagtgtccaagcagaaa
SLFN12L NM_001195790.1 IHQ 67 96 ttgaccgagaaggaatggat y ctgggagacaggtgatgatgta
106
Materiales y Métodos
Tabla 13. Serie independiente de muestras usadas en la validación de expresión por qRT-PCR.
Volumen / Volumen 96 +
Componente
reacción (µl) 10% (µl)
PreAmp Master Mix (Fluidigm, 1005876B1) 1.0 105.6
Mezcla de cebadores STA 1.25 132.0
H2O 1.5 158.4
ADNc 1.25
Total 5.00
109
Materiales y Métodos
13. Mezclar en placa 4.4 µl de pre-mezcla con 3.6 µl de ADNc preamplificado (STA).
110
Materiales y Métodos
111
Materiales y Métodos
112
Materiales y Métodos
4.3. Histopatología
113
Materiales y Métodos
114
Materiales y Métodos
absoluto, EtOH 90%, EtOH 70%, e hidratadas en PBS 1X (NaCl 137 mM, KCl 2.7 mM,
Na2HPO4 10 mM, KH2PO4 1.8 mM, pH 7.4). La recuperación del antígeno se realizó
por ebullición en tampón citrato de sodio 10mM pH 6 durante 45 minutos. Después,
las muestras se incubaron durante 20’ a TA y se lavaron dos veces con tampón
TBS-T 1X (50 mM Tris, 150 mM NaCl, Triton X-100 0.001%, pH 7.5). A continuación
las secciones se sumergieron por 10’ en H2O2 para inactivar las peroxidasas
endógenas y se realizaron tres lavados en TBS-T 1X. Los sitios de unión no
específicos fueron bloqueados con solución Protein Block durante 10’.
115
Materiales y Métodos
116
Materiales y Métodos
10 mM, pH 6 durante 15’, se lavó con TBS 1X (50 mM Tris, 150 mM NaCl, pH 7.5) y
los sitios de unión no específicos fueron bloqueados con suero de cabra al 20%
disuelto en TBS-T 1X. El anticuerpo primario αCD20 (ab9475, Abcam, Reino Unido)
se incubó en dilución 1:50 a 4°C durante la noche. Después de lavar con TBS 1X, el
anticuerpo secundario biotinilado obtenido en cabra (BA-9200, Vector Laboratories,
EUA) fue incubado en una dilución 1:200 durante 30’ y se lavó con TBS 1X. A
continuación, las muestras fueron incubadas con reactivo ABC-AP (Vectastain,
Vector Laboratories, EUA) durante 30’ y lavadas con TBS 1X. El producto de la
reacción se visualizó por incubación con solución de sustrato de fosfatasa alcalina
roja durante 20’ (Red AP Substrate Kit SK-5100, Vector Laboratories, EUA) disuelta
en 100 mM Tris-HCl pH8.5. Posteriormente, las muestras fueron lavadas con agua
del grifo y TBS-T 1x. Después de este paso se llevó a cabo el mismo protocolo de
detección de Schlafen5 previamente descrito a partir del bloqueo de peroxidasas
con H2O2. Se utilizó TBS 1X para todos los lavados finales antes del secado al aire.
117
RESULTADOS
Resultados
121
Resultados
122
Resultados
123
Resultados
Tabla 16. Características de los SNPs analizados y su asociación estadística con el CG global (cardias + no cardias/intestinal + difuso).
Número de Número de
Análisis p-valor
Gen SNP a Localización b Genotipos d controles con casos con cada OR (IC 95%) g p-valor h
bioinformáticos c perm i
cada genotipo e genotipo f
CD14 b rs11167532 Región 5’ Región 5’ GG/GA/AA 382/582/198 112/158/54 0.97 (0.81-1.16) 0,720 0.188
rs7711117 No codificante RR, Conservado CC/CT/TT 770/439/67 203/142/19 1.14 (0.94-1.38) 0,202
rs778588 Región 3’ Región 3’ TT/TC/CC 667/526/89 189/143/32 1.05 (0.87-1.26) 0,614
rs4914* Sinónimo, L367L RR, Cons CC/CG/GG 957/261/13 80/19/1 0.91 (0.69-1.19) 0,479
rs2569190 UTR 5’ RR, Triple GG/GA/AA 307/621/214 51/79/42 0.92 (0.77-1.09) 0,341
hélice,Cons
rs1583005 IK, intrónico Transcripto NMD CC/CT/TT 244/249/90 53/71/27 1.20 (1.02-1.42) 0,032
NOD2 rs7202124 SNX20, intrónico NMD, Conservado AA/AG/GG 313/225/45 97/48/5 0.74 (0.61-0.89) 0,002 0.009
rs2111235 Intrónico NMD, RR CC/CT/TT 631/507/134 205/129/26 0.77 (0.64-0.93) 0,005
rs5743289* Intrónico NMD CC/CT/TT 478/190/19 15/11/1 1.05 (0.85-1.31) 0,634
rs3135500 UTR 3’ Conservado GG/GA/AA 472/598/214 42/66/20 1.13 (0.96-1.33) 0,150
rs3785142 CYLD, intrónico Conservado TT/TC/CC 343/614/317 39/56/33 0.92 (0.78-1.08) 0,324
TLR4 rs1329061 Intergénico Intergénico TT/TC/CC 572/580/129 49/64/15 1.16 (0.98-1.39) 0,089 0.24
rs1329060 Intergénico Intergénico GG/GA/AA 948/309/26 87/36/3 1.13 (0.90-1.42) 0,302
rs1329057* Intergénico RR, intergénico, AA/AG/GG 825/410/48 62/41/4 1.13 (0.92-1.39) 0,234
Conservado
rs10491851 Región 5’ Conservado GG/GT/TT 678/516/87 217/126/21 0.81 (0.67-0.99) 0,037
rs4986790* D299G Benigna, Cons AA/AG/GG 1134/133/3 316/45/0 1.22 (0.85-1.75) 0,273
rs4986791* T399I Deletérea, Cons CC/CT/TT 1124/134/5 309/45/0 1.23 (0.86-1.77) 0,156
rs11536889* UTR 3’ RR, Conservado GG/GC/CC 940/308/35 258/98/9 1.09 (0.87-1.36) 0,457
NFKB1 rs2085549* UTR 5’ Intergénico TT/TC/CC 754/450/79 70/30/6 1.00 (0.82-1.21) 0,965 0.577
rs2085548 Intergénico Intergénico GG/GA/AA 642/531/111 57/42/8 1.00 (0.83-1.20) 1,000
rs747559 Intergénico Intergénica, Cons AA/AG/GG 417/610/177 39/50/15 1.04 (0.87-1.24) 0,669
rs3774934* Intrónico RR, Cons GG/GA/AA 1045/225/14 145/35/1 0.99 (0.75-1.31) 0,972
rs1585214 Intrónico Triple hélice GG/GA/AA 387/635/250 119/182/61 0.89 (0.75-1.05) 0,178
rs4648022 Intrónico RR CC/CT/TT 1068/206/8 298/59/8 1.23 (0.93-1.62) 0,146
rs4648055 Intrónico Conservado GG/GA/AA 656/511/117 183/155/27 0.99 (0.82-1.18) 0,896
rs9790601 Intrónico Intrónico AA/AG/GG 624/532/125 170/161/32 1.03 (0.86-1.23) 0,747
rs1609798 Intrónico Intrónico, Cons CC/CT/TT 628/528/124 49/46/12 1.02 (0.85-1.22) 0,830
rs7677509 Región 3’ RR, Región 3’ CC/CT/TT 323/627/327 32/47/26 0.88 (0.74-1.03) 0,113
rs765789 Intrónico Cons, Intrónico TT/TA/AA 845/395/45 103/66/11 1.22 (0.99-1.49) 0,062
rs228611 Intrónico RR, Conservado GG/GA/AA 341/639/297 114/175/75 0.86 (0.73-1.01) 0,073
124
Resultados
a
Los SNPs están orientados en la dirección 5’-3’ de cadena codificante. SNPs marcados con * no son fiables debido a N<5 en el genotipo de riesgo.
b
El gen CD14 está orientado en la dirección 5’-3’ en la cadena no codificante. La localización de los tagSNPs fue buscada mediante el programa
Pupasuite v3.1. c Resultado de programas Variant Effect Predictor y Pupasuite v3.1, RR significa región regulatoria, NMD significa que el SNP reside en
transcripto sujeto a mecanismo Non sense-mediated mRNA decay, Cons significa conservado, que consiste en que el SNP se localiza en región
conservada evolutivamente entre ratón y humano. d Los genotipos son reportados como alelo común/alelo menor. e, f Número de controles y casos
con los genotipos indicados en d. g OR (IC 95%) calculados para el modelo log-aditivo. En el caso de rs4986790 and rs4986791, solo el modelo
codominante fue estimado. h p-valor asociado al OR del análisis de asociación por regresión logística. i p-valor del test de permutaciones por regiones
genéticas.
Tabla 17. Asociación estadística de SNPs con respecto a la localización anatómica del cáncer gástrico.
CG tipo cardias (N=107) CG tipo no cardias (N=181)
Número de
Número de p-valor de
b b p-valor de casos con p-valor de
Gen SNP casos con cada OR (IC 95%) p-valor c OR (IC 95%) b p-valor b permutación
permutación cada c heterogeneidad d
genotipo a
genotipo a
CD14 rs11167532 40/46/10 0.69 (0.49-0.95) 0.023 0.010* 51/74/33 1.13 (0.89-1.44) 0.328 0.907 0.007
rs778588 44/52/11 1.50 (1.09-2.05) 0.013 94/70/16 1.02 (0.79-1.31) 0.881 0.060
rs2569190 36/55/15 0.69 (0.51-0.93) 0.013 51/79/42 0.99 (0.78-1.25) 0.939 0.025
rs1583005 25/55/27 1.61 (1.20-2.14) 0.001 65/83/33 1.05 (0.84-1.32) 0.645 0.056
NOD2 rs7202124 58/40/9 0.92 (0.67-1.26) 0.595 0.786 117/58/5 0.60 (0.46-0.79) 0.00019 0.0003* 0.090
rs2111235 53/39/13 0.99 (0.74-1.34) 0.961 115/56/8 0.58 (0.44-0.76) 3.664e-05 0.021
NFKB1 rs765789 71/31/5 1.09 (0.76-1.56) 0.655 0.617 103/66/11 1.42 (1.09-1.85) 0.012 0.131 0.154
rs228611 34/47/26 0.93 (0.70-1.24) 0.635 61/85/35 0.79 (0.63-0.99) 0.041 0.243
a
Reportados para los homocigotos del alelo común/heterocigotos/homocigotos del alelo menor. b OR (IC95%) y p-valores son el resultado del
análisis de regresión logística con el modelo log-aditivo ajustado. c p-valores resultantes del test de permutaciones para las localizaciones cardias y
no cardias. La significación después de la corrección de Bonferroni fue fijada en 0.05/4 genes=0.0125 (*).d p-valor resultante del test de
heterogeneidad para saber si existen diferencias en el efecto del SNP entre los grupos.
125
Resultados
Tabla 18. Asociación estadística de SNPs previamente asociados con cáncer gástrico no cardias en individuos infectados por
Helicobacter pylori con respecto al estatus de cagA.
Infectados por H. pylori y cagA+ (controles N=584, CG no cardias N=151) a Infectados por H. pylori y cagA- (controles N=687, casos no
cardias, N=27)
Número de Número de Número de Número de
c c
Gen SNP controles con casos con cada OR (IC 95%) p-valor controles con casos con cada OR (IC 95%) c p-valor c
cada genotipo b genotipo b cada genotipo b genotipo b
NOD2 rs7202124 313/225/45 97/48/5 0.66 (0.47-0.91) 0.009 92/77/25 9/3//0 0.3 (0.08-1.08) 0.035
rs2111235 295/224/60 99/43/7 0.57 (0.42-0.79) 0.0003 86/83/23 7/4//1 0.68 (0.25-1.87) 0.446
rs5743289 409/54/21 99/44/8 1.23 (0.90–1.70) 0.204 148/43/4 4/7/1 3.76 (1.33-10.63) 0.012
a
El estatus de cagA fue calculado respecto a todos los casos de CG no cardias (163) y controles (779) infectados por Helicobacter pylori. b Reportados para los
homocigotos del alelo común/heterocigotos/homocigotos del alelo menor. c OR (IC 95%) y p-valores resultantes del análisis de regresión logística con el modelo log-
aditivo ajustado para país, sexo y edad.
Tabla 19. Asociación estadística de SNPs con respecto al subtipo histológico de cáncer gástrico.
CG Intestinal (N=126) CG Difuso (N=128)
Número de Número de
p-valor de p-valor de p-valor de
Gen SNP casos con cada OR (IC 95%) b p-valor b c casos con cada OR (IC 95%) b p-valor b
a permutación a permutación c heterogeneidad d
genotipo genotipo
CD1 rs778588 52/59/15 1.49 (1.12-1.97) 0.006 0.041* 76/41/10 0.80 (0.59-1.09) 0.151 0.652 0.010
4
NOD rs7202124 76/43/7 0.73 (0.54-0.99) 0.040 0.077 78/44/5 0.72 (0.53-0.99) 0.037 0.204 0.922
2 rs2111235 74/44/7 0.68 (0.50-0.93) 0.011 69/50/8 0.83 (0.62-1.10) 0.191 0.232
NFK rs1585214 38/66/21 0.94 (0.72-1.23) 0.665 0.727 46/65/16 0.76 (0.58-0.99) 0.044 0.213 0.129
B1 rs228611 40/53/32 0.93 (0.72-1.21) 0.605 39/74/15 0.74 (0.56-0.96) 0.024 0.144
rs765789 75/44/7 1.31 (0.96-1.80) 0.096 73/43/11 1.46 (1.07-1.99) 0.018 0.280
rs7677509 36/54/35 0.97 (0.75-1.26) 0.819 36/72/19 0.77 (0.59-1.00) 0.047 0.107
a b
Reportados para los homocigotos del alelo común/ heterocigotos/homocigotos del alelo menor. OR (IC 95%) y p-valores resultantes del análisis de regresión
logística con el modelo log-aditivo ajustado para país, sexo y edad. c p-valores resultantes del test de permutaciones para localizaciones cardias y no cardias. La
significación después de la corrección de Bonferroni fue fijada en 0.05/4 genes=0.0125, por lo que un p-valor<0.0125 fue considerado significativo (*). d p-valor
resultante del test de heterogeneidad realizado para saber si existen diferencias en el efecto del SNP entre los grupos.
126
Resultados
Tabla 20. Asociación estadística de haplotipos respecto a los subtipos histológicos y anatómicos de cáncer gástrico.
Frecuencia Frecuencia
Subtipo histológico o
Gen Haplotipo SNPs que componen el haplotipo a haplotípica haplotípica OR (IC 95%) b p-valor c
anatómico de CG
en controles en casos
CD14 ACTCAC rs11167532,rs7711117, 0.413 ----- 1.00 (Referencia) ----- -----
rs778588,rs4914,rs2569190,
rs1583005
CD14 GCCCGT Idem 0.128 0.212 2.15 (1.41-3.27) 0.0004 Cardias
CD14 GCCCGT Idem 0.128 0.193 1.71 (1.17-2.49) 0.006
* Intestinal
NOD2 ACCAT rs7202124, rs2111235, rs5743289, 0.355 ----- 1.00 ----- -----
rs3135500, rs3785142 (Referencia)
NOD2 GTCGC Idem 0.198 0.135 0.60 (0.43-0.86) 0.005 No cardias
NOD2 GTCGC Idem 0.198 0.141 0.64 (0.43-0.96) 0.031 Intestinal
TLR4 TGATACG rs1329061, rs1329060, rs1329057, ----- 1.00 ----- -----
rs10491851, rs4986790, 0.267 (Referencia)
rs4986791, rs11536889
TLR4 CGAGACG Idem 0.074 0.112 1.98 (1.14-3.47) 0.016 Cardias
NFKB1 TGAGACGACTTA rs2085549, rs2085548, rs747559, ----- 1.00 ----- -----
rs3774934, rs1585214, rs4648022, (Referencia)
rs4648055, rs9790601, rs1609798, 0.328
rs7677509, rs765789, rs228611
NFKB1 TGAGGTGACCAG Idem 0.066 0.095 1.78 (1.15-2.73) 0.009 No cardias
a b c
Orden de los SNPs que componen los haplotipos. Odds ratio (IC 95%) y p-valores, ajustados por país, sexo y edad para haplotipos asociados
significativamente (p<0.05) con los subtipos histológicos o anatómicos de CG, los cálculos fueron realizados para el modelo log-aditivo y el haplotipo
de referencia es el común.
127
Resultados
Tabla 21. Búsqueda de potenciales eQTLs entre los SNPs asociados y no asociados a CG.
SNPs SNPs no r2/p-valor o Transcripto Tejido y Bases
Gen
asociados a asociados a score a regulado referencia de
CG CG datos b
-8 153
rs11167532 ------ CD14 0.19/4.29 x 10 WDR55 Cerebro A
rs1583005 ------ CD14 0.36/3.17 x 10-15 WDR55 Cerebro 153 A
rs2569190 ------ CD14 0.23/1.54 x 10-9 WDR55 Cerebro 153 A
157
rs1583005 ------ CD14 30.50 WDR55 Diversos B
156
rs778588 ------ CD14 34.51 SRA1 Monocitos B
rs7202124 ------ NOD2 42.18 NOD2 Monocitos 156 B
156
rs7202124 ------ NOD2 16.77 SNX20 Monocitos B
rs5743289 ------ NOD2 3.903 CYLD Linfoblastoide 193 B
------ rs228611 NFKB1 7.7 MANBA Linfoblastoide 193 B
------ rs7677509 NFKB1 52.94 MANBA Diversos 157 B
156
rs747559 NFKB1 27.318 COL9A2 Monocitos B
rs4648055 NFKB1 23.828 COL9A2 Monocitos 156 B
155
rs1609798 NFKB1 0.012 MANBA Linfoblastoide B
a 2
Valor de r y p-valor asociado de la correlación entre la asociación del SNP con los niveles de
expresión génica para la base de datos A y score para la B. b A: Base de datos eQTL browser (NCBI), B:
Base de datos eQTL de Gilad/Pritchard.
El ADN obtenido por Isohelix Dri‐Capsules mostró una gran variación en la masa
obtenida (5.53 ± 5.17 μg), relaciones de absorbancia 260/280 indicativas de
contaminación proteica y con sales, y se observó degradación del mismo, como se
evidencia en la siguiente electroforesis (Figura 19).
128
Resultados
Mediante los hisopos bucales BuccalAmp™ DNA Extraction Kit el rendimiento fue
menor y también mostró una gran variabilidad (3.01 ± 3.62 μg), pero menor
degradación del ADN obtenido (Figura 20). Al igual que en el kit anterior las
relaciones de absorbancia 260/280 fueron indicativas de contaminación proteica y
con sales.
El ADN de saliva obtenido mediante este kit mostró una degradación insignificante
(Figura 21), relaciones de absorbancia 260/280 >1.8 demostrando poca
contaminación proteica y con sales, y una elevada masa de ADN (125 ± 55.9 µg),
suficiente para realizar el genotipado propuesto y que quedara material para
futuros proyectos.
129
Resultados
Como se obtuvó una mayor masa de ADN con la solución 4 respecto a la 3 se optó
por ella como la mejor. Su composición es CDTA 3.3 mM, Tris base 33 mM, SDS
0.6%, LiCl 0.67 M, Urea 0.67 M, EtOH 30%, pH 8, proteinasa K 150 ng/µl. El
conocimiento de la acción de sus componentes contribuyo a dicha elección:
La electroforesis del ADN extraído usando esta solución evidenció que la mayoría de
las muestras poseían un adecuado tamaño para ser utilizadas en el genotipado de
SNPs (Figura 22). Al ser aplicado este método a la colecta de mayores cantidades de
muestras de pacientes en condiciones no ideales respecto a las recomendadas, se
observó que la masa (23.5 ± 7.8 µg) y el estado de conservación del ADN (Figura 23)
seguía siendo adecuado para el genotipado de SNPs.
130
Resultados
Figura 22. Electroforesis en gel de agarosa 1% del ADN extraído a partir de muestras
de saliva conservadas durante 60 días en solución estabilizadora # 4 a
temperatura ambiente. Leyenda: 1, 2, 4-7 y 9-12 muestras, 3 y 8
marcador de peso de ADN 1 Kb ladder (Invitrogen, EUA).
Figura 23. Electroforesis en gel de agarosa 1% del ADN extraído a partir de muestras
de saliva colectadas por los hospitales en solución estabilizadora # 4.
Leyenda: 2-13 muestras, 1 y 14 marcador de peso de ADN 1 Kb ladder
(Invitrogen, EUA).
131
Resultados
también podría ser genotipado por otra plataforma diferente como la Sequenom, la
cual había sido propuesta por el CEGEN como la adecuada de acuerdo a las
características de nuestro estudio.
Tabla 23. Resultados del genotipado del SNP rs7578034 en gen IL1B
en muestras de ADN extraído de saliva.
Muestras Call Score
10001 Alelo X 0.94
10002 XY 0.99
10003 XY 0.98
10004 Alelo Y 0.9
10006 Alelo X 0.84
10007 XY 0.97
10008 Alelo X 0.87
10009 XY 0.93
10009 Alelo Y 0.9
10010 XY 0.95
Control + Alelo Y 0.98
Control + XY 0.98
Control + XY 0.91
Control + XY 0.98
Control + Alelo X 0.98
Control + Alelo Y 0.94
H2O Negativo 0.00
H2O Negativo 0.00
Figura 24. Genotipado del SNP rs7578034 en el gen IL1B empleando LightCycler II.
132
Resultados
133
Resultados
Los resultados del genotipado de los 141 SNPs analizados en cada uno de los grupos
de evolución de las lesiones se indican en la tabla anexa 1. La asociación más
significativa (OR=0.63, IC95% 0.48-0.81, p=0.0003, modelo aditivo) se observó en el
SNP rs10902073 de MUC2 asociado inversamente con la regresión de las LPGs; fue
el único polimorfismo significativo después de aplicar la corrección para
comparaciones múltiples mediante FDR. Los SNPs asociados con un incremento del
riesgo de progresión se encontraron en los genes CD14, DNAH11, MAPK3, PRKAA1,
RUNX3 y TFF1. Otros SNPs localizados en los genes CD14, CDH1, PTGS2, TFF1 y TFF2,
se asociaron inversamente con la progresión de las lesiones (Tabla 25A). Los
resultados del SNP rs10760634 no fueron considerados fiables debido al amplio
rango de los límites de confianza del 95% en el cálculo del OR.
El análisis de la regresión mostró que SNPs localizados en los genes CD14, IL1B,
IL1RN, NFKB1, TFF1 y TFF2 se asocian con una mayor regresión; mientras que SNPs
en CDH1, MUC1, MUC2, NFKB1 y TFF1 se asocian con una menor regresión.
También, existen algunos SNPs cuyos resultados no fueron considerados a causa del
amplio rango de los límites de confianza del 95% en el cálculo del OR (Tabla 25B). A
continuación, se exploró el sentido de la asociación de aquellos SNPs que se
asociaron tanto a progresión como a regresión de las lesiones y se eliminó del
análisis a aquellos en los que existía el mismo sentido en ambos análisis, rs1598861
y rs7674640 de NFKB1, por considerarlos no informativos respecto a la evolución de
las LPGs (Tabla 26).
134
Resultados
135
Resultados
Tabla 25. OR (IC 95%) de los SNPs asociados significativamente con progresión y regresión de las LPGs.
A) Progresión
Modelo Número de Número de
Gen SNP a Genotipos controles con casos con cada OR (IC 95%) p-valor Variant Effect PupaSuite f
genético b
cada genotipo c genotipo d Predictor e
CD14 rs11167532 log-aditivo GG/AG/AA 118/211/104 42/42/11 0.63 (0.46-0.87) 0.004 Región 5’ -------------
CD14 rs2569190 log-aditivo GG/AG/AA 113/224/123 34/47/15 0.65 (0.47-0.89) 0.006 Intrónico, UTR 5’, Conservada/
RR Triple hélice
CD14 rs5744455 log-aditivo CC/CT/TT 293/154/14 49/40/7 1.63 (1.13-2.36) 0.01 Región 5’, RR Conservada
CD14 rs1583005 log-aditivo CC/CT/TT 152/235/74 21/50/25 1.56 (1.13-2.16) 0.006 Intrónico, NMD -------------
CDH1 rs4783573 recesivo AA/AG/GG 196/205/60 42/48/5 0.37 (0.14-0.95) 0.02 Intrónico, NMD -------------
DNAH11 rs2285947 dominante GG/AG/AA 128/217/116 20/56/20 1.57 (1.01-2.46) 0.045 Intrónico Conservada
MAPK3 rs11865086 dominante AA/AC/CC 123/222/116 16/57/23 1.82 (1.02-3.23) 0.033 UTR 3’, NMD Conservada
NFKB1 rs1598861 * recesivo AA/AC/CC 201/184/77 46/42/8 0.45 (0.21-0.98) 0.028 Intrónico Conservada
NFKB1 rs7674640 * recesivo TT/CT/CC 115/217/130 30/51/14 0.46 (0.26-0.83) 0.004 Región 3’, RR Conservada
PRKAA1 rs13361707 log-aditivo CC/CT/TT 240/193/29 44/42/10 1.44 (1.02-2.02) 0.038 Intrónico, RR Conservada
PTGES rs10760634 * recesivo AA/AG/GG 392/69/2 79/14/3 7.44 (1.23-45.11) 0.031 Intrónico Conservada
PTGS2 rs5275 log-aditivo TT/CT/CC 216/194/50 58/32/6 0.64 (0.45-0.92) 0.013 UTR 3’ Conservada
PTGS2 rs4648276 dominante TT/CT/CC 315/131/14 77/19/0 0.52 (0.31-0.87) 0.03 Intrónico, NMD Conservada
RUNX3 rs6663310 dominante TT/CT/CC 215/191/56 35/51/10 1.61 (1.03-2.5) 0.035 Región 3’ -------------
TFF1 rs9976977 recesivo GG/AG/AA 176/231/54 29/47/20 1.46 (1.05-2.02) 0.023 Region 3’ -------------
TFF1 rs424694 log-aditivo CC/CT/TT 162/217/84 44/41/11 0.69 (0.5-0.96) 0.025 Region 3’ -------------
TFF2 rs1079380 log-aditivo AA/AG/GG 113/243/107 32/49/15 0.7 (0.51-0.98) 0.033 Intrónico Conservada
a
* Resultados para estos SNPs no considerados debido a que la dirección de asociación en el análisis de progresión y regresión tiene el mismo sentido (SNPs de NFKB1)
o a un excesivo rango en los límites de confianza del 95% del OR (rs10760634). b Modelo genético con mayor significación estadística. c d Número de controles
(regresión + estables) y de casos (progresión) de homocigotos para el alelo más común, heterocigotos y homocigotos para el alelo menos común,respectivamente.
eyf
Predicción de efecto funcional en los programas indicados. Conservada se refiere a región conservada evolutivamente entre humanos y ratón.
136
Resultados
B) Regresión
Gen SNP a Modelo Genotipos Número de Número de OR (IC 95%) p-valor PupaSuite f
genético b controles con casos con cada Variant Effect
cada genotipo c genotipo d Predictor e
CD14 rs778588 Sobre- TT/CT/CC 187/138/30 86/103/15 1.6 (1.13-2.27) 0.008 Región 3’ -------------
CDH1 rs16260 log-aditivo
dominante CC/AC/AA 173/146/33 113/80/10 0.76 (0.57-1) 0.048 Región 5’ Conservada
CDH1 rs12597188 log-aditivo GG/AG/A 152/159/43 101/84/15 0.75 (0.58-0.98) 0.036 Intrónico, NMD -------------
CDX2 rs4503658 * recesivo GG/GT/TT
A 312/41/2 169/30/5 1.54 (1.01-2.37) 0.047 UTR 5’, RR Conservada
IL1B rs1143627 dominante TT/CT/CC 165/152/36 77/103/24 1.45 (1.02-2.06) 0.038 Región 5’, RR Conservada
IL1RN rs2637988 recesivo AA/AG/G 139/165/51 80/80/44 1.64 (1.05-2.56) 0.031 Intrónico, Región -------------
MAP3K14 rs16939926 dominante TT/CT/CC
G 305/49/1 188/16/0 0.52 (0.29-0.94) 0.023 Intrónico
3’ Conservada
MUC1 rs2070803
* dominante TT/CT/CC 106/168/80 79/86/38 0.67 (0.47-0.96) 0.031 Conservada Conservada
MUC2 rs10902073 log-aditivo CC/AC/A 141/155/59 110/76/18 0.63 (0.48-0.81) 0.0003 * Intergénico, RR -------------
MUC2 rs10794281 log-aditivo TT/CT/CC
A 123/158/74 85/93/26 0.55 (0.34-0.9) 0.014 Intergénico Conservada
MUC6 rs6597947 * recesivo GG/GT/TT 307/46/2 172/25/6 5.38 (1.07-26.89) 0.025 Región 3’ Conservada
NFKB1 rs980455 log-aditivo AA/AG/G 144/153/55 99/97/41 1.29 (1.01-1.65) 0.038 Región 3’ -------------
NFKB1 rs1598861 recesivo AA/AC/CC
G 153/137/64 94/89/21 0.52 (0.31-0.88) 0.012 Intrónico Conservada
NFKB1 rs7674640 recesivo TT/CT/CC 87/164/103 58/104/41 0.62 (0.41-0.93) 0.019 Región 3’, RR Conservada
PTGS2 rs5277 recesivo GG/CG/C 225/109/21 139/61/4 0.32 (0.11-0.94) 0.021 Sinónimo, NMD Conservada
SRC rs6017901 * dominante TT/CT/CC
C 314/39/1 193/11/0 0.45 (0.22-0.89) 0.016 Región 3’ -------------
TFF1 rs4920094 dominante GG/AG/A 97/194/64 75/93/36 0.65 (0.45-0.93) 0.021 Intrónico Conservada
TFF1 rs225357 recesivo CC/CT/TT
A 121/179/54 61/97/46 1.62 (1.04-2.51) 0.032 Intrónico Conservada
TFF1 rs424694 recesivo CC/CT/TT 134/171/50 72/87/45 1.73 (1.11-2.7) 0.017 Region 3’ -------------
TFF1 rs13047838 dominante CC/CT/TT 178/152/25 121/74/9 0.72 (0.54-0.97) 0.026 UTR 3’ Conservada
TFF2 rs1079380 recesivo AA/AG/G 99/189/67 46/103/55 1.59 (1.06-2.39) 0.027 Intrónico Conservada
, Triple
a
* Resultados para estos SNPs no considerados
G debido a que el número de individuos con el genotipo de riesgo es muy bajo o el intervalo de confianza
hélice
b c d
muy amplio. Modelo genético con mayor significación estadística. Número de controles (regresión + estables) y de casos (progresión) de
homocigotos para el alelo más común, heterocigotos y homocigotos para el alelo menos común,respectivamente. e y f Localización del SNP y predicción
de efecto funcional en los programas indicados. Conservada se refiere a región conservada evolutivamente entre humanos y ratón.
137
Resultados
138
Resultados
Cuando comparamos los resultados obtenidos en este estudio con los obtenidos
con SNPs genotipados en un estudio previo de seguimiento de lesiones precursoras,
realizado en la provincia de Soria por nuestro grupo, se confirmaron las
asociaciones inversas de CD14 (rs2569190 y rs5744455), TFF2 (rs1079380) con la
progresión. También que se mantuvo la asociación de MUC2 (rs10902073 y
rs10794281) y TFF1 (rs424694) con la regresión (Tabla 26). Por el contrario, no se
replicaron en el presente estudio las asociaciones de CDH1 (rs1125557), TFF1
(rs9976977, rs13047838 y rs225358), MUC2 (rs7944723, rs10794293, rs3924453 y
rs4077759), IL1RN (rs2637988) e IL1A (rs17561).
Al analizar genes que mostraron asociación a nivel de SNP y haplotipos, este análisis
confirmó las asociaciones a nivel de SNP (Tabla 27), es decir, correspondiéndose con
lo esperado atendiendo al sentido de asociación de sus SNPs componentes. Se
encontró asociación con mayor riesgo de progresión de los haplotipos TGTG de
CD14 (OR=1.90, IC95% 1.26-2.87), así como CAAC (OR=1.82, IC95% 1.10-2.98) y
CAGC (OR=1.89, IC95% 1.03-3.46) de TFF1. La combinación de CD14 contiene al
alelo de riesgo T de rs57444455 que se asocia positivamente con la progresión. El
haplotipo GCCGAG de PTGS2 (OR=0.55, IC95% 0.31-0.99) se asocia inversamente
con la progresión, y el alelo C de rs5275 se asocia inversamente con la progresión.
139
Resultados
Tabla 28. Análisis de asociación en portadores de la combinación de los factores de virulencia de H. pylori cagA+vacAs1/m1 vs otras
combinaciones, en individuos infectados por H. pylori.
Progresión en cagA+vacAs1/m1 Progresión en el resto de combinaciones
a Modelo # de controles # de casos c d # de controles # de casos
Gen SNP OR (IC 95%) p-valor OR (IC 95%) c p-valor d
genetico b con cada con cada con cada con cada
genotipo b genotipo c genotipo b genotipo c
CD14 rs5744455* recesivo 72/35/1 13/7/3 16.05 (1.59-162.2) 0.010 129/64/4 17/14/1 1.56 (0.17-14.39) 0.708
CDX1 rs6894617 recesivo 52/49/7 8/10/5 4.01 (1.15-14.02) 0.038 110/75/13 21/9/2 0.95 (0.2-4.42) 0.946
GAST rs12453761 log-aditivo 31/51/26 11/10/2 0.49 (0.25-0.97) 0.034 50/106/42 9/16/7 0.95 (0.55-1.65) 0.865
NFKB1 rs7674640* recesivo 26/46/36 7/15/1 0.09 (0.01-0.7) 0.001 51/69/51 11/15/6 0.67 (0.26-1.71) 0.382
RUNX3 rs2282718 dominante 54/44/10 6/13/4 2.83 (1.04-7.73) 0.033 84/93/21 11/13/8 1.41 (0.64-3.07) 0.387
SRC rs6017916 recesivo 58/43/7 11/6/6 5.09 (1.53-17) 0.011 112/68/18 16/12/4 1.43 (0.45-4.53) 0.557
SRC rs6094373 recesivo 64/38/6 13/5/5 4.72 (1.3-17.12) 0.024 122/62/14 18/10/4 1.88 (0.58-6.11) 0.319
TFF1 rs225355 sobredominant 26/54/28 10/6/7 0.35 (0.13-0.96) 0.033 63/88/45 10/15/6 1.15 (0.54-2.46) 0.717
e
CD14 rs11167532* recesivo 37/47/24 9/11/3 0.53 (0.14-1.92) 0.302 61/90/47 11/19/2 0.21 (0.05-0.93) 0.012
MAP3K14 rs16939926* log-aditivo 96/12 18/5 2.22 (0.7-7.08) 0.194 177/21/0 24/8/0 2.81 (1.12-7.04) 0.037
MAP3K14 rs2074292 log-aditivo 24/63/21 5/15/3 0.86 (0.42-1.76) 0.684 42/108/46 11/17/3 0.53 (0.29-0.95) 0.030
MAP3K14 rs4247364 dominante 43/51/14 5/15/3 2.38 (0.82-6.89) 0.091 84/93/20 6/20/6 3.22 (1.27-8.18) 0.007
PTGS2 rs4648276* dominante 73/34/1 17/6/0 0.74 (0.27-2.03) 0.547 129/61/6 27/5/0 0.36 (0.13-0.97) 0.027
PTGS2 rs5275* dominante 53/46/9 14/7/2 0.62 (0.25-1.55) 0.302 85/88/25 24/7/1 0.25 (0.11-0.59) 0.001
PTGES rs10760634 log-aditivo 93/14/1 22/1/0 0.29 (0.04-2.21) 0.562 164/33/1 22/8/2 2.34 (1.14-4.83) 0.027
RUNX3 rs2282718 recesivo 54/44/10 6/13/4 2.06 (0.59-7.27) 0.279 84/93/21 11/13/8 2.81 (1.12-7.04) 0.037
RUNX3 rs6663310* dominante 56/41/11 12/9/2 0.99 (0.4-2.43) 0.978 88/80/29 8/19/5 2.42 (1.04-5.66) 0.032
NFKBIA rs696 dominante 38/60/10 9/13/1 0.84 (0.33-2.13) 0.72 74/93/31 5/21/5 3.1 (1.14-8.43) 0.015
WNT1 rs4760663 log-aditivo 38/49/21 11/10/2 0.62 (0.31-1.21) 0.147 69/95/33 7/15/10 1.74 (1.02-2.96) 0.041
a
El asterisco * indica SNPs que fueron asociados con la evolución de las LPGs en el análisis general y los no marcados con * indica que son nuevas
asociaciones. b Modelo genético donde fue obtenida la asociación más significativa. c, d Controles (regresión + estables) y casos (progresión).
d
OR (IC 95%) y p-valor asociado del modelo genético con la asociación más significativa.
141
Resultados
Como en el caso de la asociación con CG los SNPs de CD14 asociados con la evolución de
las LPGs, rs11167532, rs1583005, rs778588 y rs2569190 se asocian con niveles de
expresión génica de ARNm y los mismos se encuentran en desequilibrio de ligamiento.
Otros SNPs asociados con la evolución de las LPGs pertenecientes a genes NFKB1, MAPK3
e IL1B fueron encontrados como eQTLs reguladores de la expresión génica. Mientras que
SNPs no asociados con la evolución de LPGs también demostraron ser eQTLs (Tabla 29).
143
Resultados
3.2. Prueba piloto de comparación de microarrays Human Gene 2.0 ST versus Almac
Xcel
Los análisis de calidad demostraron que ambos arrays eran adecuados para trabajar con
tejido parafinado. Sin embargo, con el Almac-Xcel se obtuvo una mayor cantidad de
genes diferencialmente expresados (N=266) respecto al Human Gene 2.0 ST (N=11). Entre
los genes diferencialmente expresados obtenidos con el Almac-Xcel se encontraron
múltiples genes codificantes para proteínas con funciones en la fisiología gastrointestinal
que se conoce están desreguladas en la MI, Ej. OLFM4, LYZ, CDH17, etc. 81. Los transcritos
de relevancia en la MI obtenidos con el Human Gene 2.0 ST se restringían a ATP4A/B
(ATPase, H+/K+ exchanging), GIF (gastric intrinsic factor) y PGA3 (pepsinogen 3), los
cuales a su vez estaban comprendidos en los resultados del Almac-Xcel. Con base en su
amplia cobertura y que el Almac Xcel es un array optimizado para tejido se decició por el
mismo para el análisis de muestras clínicas.
144
Resultados
Tabla 30. Medición de cantidad y calidad de las muestras de ARN usadas en la prueba piloto de comparación de microarrays
Human Gene 2.0 ST y Almac Xcel.
Tipo de MI al Extensión de MI Localización Fin de Infección c(ARN) A c(ADN)
Código Estatus c Sexo d RIN
reclutamiento a al reclutamiento b anatómica seguimiento por H. pylori (ng/µl) 260/280 (ng/µl) e
A1 No presenta No presenta Cuerpo Sanos Sanos SI H 593.35 2.1 2.1 7.5
A2 No presenta No presenta Cuerpo Sanos Sanos NO H 465.85 1.97 2.2 12
A3 No presenta No presenta Cuerpo Sanos Sanos NO M 487.25 2 2.2 9.8
A4 G1 E4 Antro CG intestinal MII-CG SI H 148.56 2.03 2.4 9.6
A5 G3 E4 Antro CG intestinal MII-CG SI M 304.0 2.02 2.4 7
A6 G3 E4 Antro CG mixto MII-CG NO H 918.74 2.06 2.4 10.3
a,b
Según se indica en apartado 3.1 de Materiales y Métodos.c Sanos: Mucosa gástrica sana, MII-CG: Metaplasia intestinal incompleta que
progresa a CG d H: Hombre, M: Mujer. e Concentración de ADN residual en la muestra de ARN.
Figura 27. Electroforesis de ARN total obtenido de muestras parafinadas en Bioanalyzer 2100 (Agilent, EUA).
145
Resultados
Tabla 31. Medición de cantidad y calidad de las muestras usadas en el microarray de expresión.
Extensión de
Diagnóstico al Localización Diagnóstico al final Años de c(ARN) c(ADN)
Código MI al Estatus d Sexo e Edad f RIN
reclutamiento a anatómica de seguimiento seguimiento c (ng/µl) (ng/µl) g
reclutamientob
A7 G4 E4 Antro CG intestinal 10 MII-CG M 52 198.7 2.4 3.5
A8 G4 E4 Antro stinalntestinaIntesti
CG intestinal 8.4 MII-CG H 69 199.5 2.5 6.8
A9 G4 E3 Antro nal Difuso
CG 7.7 MII-CG H 45 254.8 2.4 7.9
A10 G4 E3 Antro CG intestinal 2 MII-CG H 73 188 2.4 5.6
A11 G4 E3 Antro CG intestinal 1 MII-CG H 60 303.9 2.1 10
A12 G4 E3 Antro CG intestinal 11.8 MII-CG M 66 157.8 2.6 8.6
A13 G3 E4 Antro CG 13.4 MII-No CG H 53 38.7 2.1 1.3
A14 G4 E4 Antro MII 15 MII-No CG H 64 53.6 2.3 4.8
A15 G3-G4 E3-E4 Antro MII 11 MII-No CG M 64 144 2.6 <2
A16 G3-G4 E4 Antro MII ND MII-No CG H 64 29.7 NA 1.8
A17 G3-G4 E3-E4 Antro MII 14 MII-No CG H 65 114.8 1.3 14
A18 G3-G4 E4 Antro MII ND MII-No CG M 74 114.6 1.5 <2
A19 G3-G4 E4 Antro MII ND MII-No CG H 71 136.1 2.1 8.7
A20 G1 E4 Antro CG no cardias 2 MIC-CG H 57 378 2.3 3.26
A21 G1 E4 Cuerpo CG intestinal 8 MIC-CG H 57 183 2.4 <2
A22 G1 E4 Antro CG intestinal 1 MIC-CG H 44 286 2.3 <2
A23 G1 E4 Antro MIC 15.37 MIC-CG H 56 69.3 2.4 6.4
A24 G2 E4 Antro MIC 14 MIC-CG M 59.74 62.9 2.2 8.14
A25 G1 E4 Antro MIC 8.24 MIC-CG M 57 215.5 2.1 9.65
A26 G1 E4 Antro MIC 5.95 MIC-CG M 78 119.6 2 12.3
A27 G1 E4 Antro CG no cardias 6.2 MIC-CG H 76 356.63 2.2 11.3
A28 G1 E4 Transicional MIC 10.6 MIC-No CG H 67 110.9 2.6 <2
A29 G2 E4 Antro MIC 10.5 MIC-No CG H 63 38.2 2.4 1.5
A30 G2 E4 Transicional MIC 12.31 MIC-No CG H 74 50.7 2.4 <2
146
Resultados
Extensión de Diagnóstico al
Diagnóstico al Localización Años de c(ARN) c(ADN)
Código MI al final de Estatus d Sexo e Edad f RIN
reclutamiento a anatómica seguimiento c (ng/µl) (ng/µl) g
reclutamientob seguimiento
A31 G2 E4 Antro MIC 13 MIC-No CG H 56.98 104.9 2.4 <2
A32 G1 E4 Antro MII 14 MIC-No CG M 54 45.2 2.3 <2
A33 E4 G1 Antro MII 15 MIC-No CG H 57 98.6 2.5 12.5
A34 E4 G1 Antro MIC 18 MIC-No CG H 57 57.4 2.4 5.9
A35 E4 G1 Antro MIC 14 MIC-No CG H 62 67.7 2.2 6.3
A36 E4 G1 Antro MIC 15 MIC-No CG H 62 181 2.4 10.2
A37 -------- -------- Antro -------- -------- Sanos M 16 58 2.4 2.9
A38 -------- -------- Antro -------- -------- Sanos H 42 99.7 2.4 3.3
A39 -------- -------- Antro -------- -------- Sanos H 39 78.5 2 35.6
A40 -------- -------- Antro -------- -------- Sanos H 41 274 2.3 9.34
A41 -------- -------- Antro -------- -------- Sanos M 72 104 1.1 9.56
A42 -------- -------- Antro -------- -------- Sanos M 86 108 2.1 <2
A43 -------- -------- Antro -------- -------- Sanos M 18 77.4 2.2 <2
A44 -------- -------- Cuerpo -------- -------- Sanos H 27 139 2.3 3.79
A45 -------- -------- Antro -------- -------- Sanos H 68 111 1.6 3.25
A46 -------- -------- Antro -------- -------- Sanos M 31 244 2.2 12.3
A47 -------- -------- Antro -------- -------- Sanos M 19 476 2.4 5.47
A48 -------- -------- Cuerpo -------- -------- Sanos H 73 472 2.5 8.25
A49 -------- -------- Cuerpo -------- -------- Sanos M 58 409 2.4 3.96
A50 -------- -------- Antro -------- -------- Sanos M 11 286 2.3 2.9
A51 -------- -------- Cuerpo -------- -------- Sanos H 69 153 2.3 3.13
a
Tipos histológicos de MI (G1-G4) descritos en apartado 3.1 de Materiales y Métodos. b E3>50%, E4>75%: Extensión de la MI en la muestra. En el caso de que
existan dos diferentes en a ó b en la misma muestra se corresponden a los cortes inicial y final. c Se refiere al intervalo de tiempo entre el reclutamiento y el final de
seguimiento. d MII/MIC-CG: Muestra de MII o de MIC al reclutamiento que progresa a CG al final del seguimiento. MII/MIC-No CG son controles de MII y MIC que no
progresan a CG desde el reclutamiento hasta el final de seguimiento. e M: Mujer, H: Hombre. ND: No disponible. f Edad en el reclutamiento. g Concentración de ADN
residual en la muestra de ARN.
147
Resultados
Figura 28. Dendrograma de todos los grupos analizados usando la distancia Average
y método Complete en el programa informático Combat. Leyenda: Indica
si la muestra progresa (MII-CG/MIC-CG) o no (MII-No CG/MIC-No CG) a
CG, o si proviene de mucosa gástrica sana (Sanos).
148
Resultados
149
Resultados
Se analizaron los datos de expresión de 57 263 sondas con desviación estándar (DE)
<0.34 (30% de la DE global) de un total de 81 804, por lo que fueron eliminadas
24 541 sondas del análisis. El motivo de filtrar es eliminar las sondas que no se
expresan suficientemente, produciendo “ruido” al hacer el análisis estadístico. El
punto de corte escogido corresponde a las sondas que tienen una DE>30% del
resto, y esto corresponde a una DE de 0.34.
150
Resultados
Figura 32. Diagrama de Venn que superpone los genes significativos en las
comparaciones relativas a la MII.
151
Resultados
Figura 33. Diagrama de Venn que superpone los genes significativos en las
comparaciones relativas a la MIC.
600
Sobre-expresados
400 Sub-expresados
200
0
G
G
C
C
-C
II-
o
IC
I-N
M
152
Resultados
Tabla 32. Genes diferencialmente expresados entre los grupos MII-CG vs MII-No CG.
Sign en Sign en MII- Candidatos
Fold Asociados con MI o
Símbolo del gen a Nombre del gen b p-valor c MII-CG vs No CG vs de mayor
Change CG d
Sanos e Sanos f interés g
HLA-DRB1 major histocompatibility complex, class II, DR beta 1 0.260 0.014 SI 195 SI NO SI
RNU2-2 RNA, U2 small nuclear 2 0.297 0.008 Nueva molécula SI NO SI
MME membrane metallo-endopeptidase 0.363 0.007 Nueva molécula NO SI NO
MGAM maltase-glucoamylase (alpha-glucosidase) 0.374 0.000 SI 196 NO SI NO
IK IK cytokine, down-regulator of HLA II 0.384 0.002 SI 192 SI NO SI
TMEM25 transmembrane protein 25 0.386 0.005 Nueva molécula NO SI NO
GSTA1 glutathione S-transferase alpha 1 0.386 0.029 SI 197 NO SI NO
PPIA Peptidylprolyl isomerase A (cyclophilin A) 0.393 0.027 SI 198 SI NO SI
KARS lysyl-tRNA synthetase 0.396 0.032 Nueva molécula NO NO SI
RHOA ras homolog family member A 0.398 0.000 SI 199 NO SI NO
CDKL3 cyclin-dependent kinase-like 3 0.403 0.000 Nueva molécula SI NO SI
TMEM236 transmembrane protein 236 0.409 0.001 Nueva molécula NO SI NO
RBP2 retinol binding protein 2, cellular 0.411 0.021 SI 200 SI SI NO
CP ceruloplasmin (ferroxidase) 0.412 0.040 SI 201 SI NO SI
GIP gastric inhibitory polypeptide 0.433 0.017 Nueva molécula SI SI NO
XPNPEP2 X-prolyl aminopeptidase (aminopeptidase P) 2, 0.442 0.020 Nueva molécula SI SI NO
GNAS membrane-bound
GNAS complex locus 0.445 0.009 SI 202 NO NO SI
SART3 squamous cell carcinoma antigen recognized by T 0.448 0.000 SI 203 NO NO SI
RBBP7 retinoblastoma
cells 3 binding protein 7 0,449 0.001 Nueva molécula SI NO SI
AQP10 aquaporin 10 0.465 0.014 Nueva molécula NO SI NO
SLC13A2 solute carrier family 13 (sodium-dependent 0.470 0.007 Nueva molécula SI SI NO
DNAJB7 dicarboxylate
DnaJ (Hsp40) homolog, subfamily transporter),
B, member 7 0.482 0.001 Nueva molécula NO SI NO
XDH member 2
xanthine dehydrogenase 0.486 0.006 Nueva molécula SI SI NO
DTNB dystrobrevin, beta 0.487 0.001 Nueva molécula NO NO SI
LCT lactase 0.493 0.002 Nueva molécula NO NO SI
153
Resultados
ATP1B3 ATPase, Na+/K+ transporting, beta 3 polypeptide 0.499 0.002 Nueva molécula SI SI NO
FAM133B/ family with sequence similarity 133, member B 2.001 0.004 Nueva molécula NO NO SI
FAM133CP/FAM133
DP
MORF4L1 mortality factor 4 like 1 2.006 0.016 Nueva molécula NO NO SI
CXCL14 chemokine (C-X-C motif) ligand 14 2.007 0.028 SI 204 NO NO SI
LAPTM5 lysosomal protein transmembrane 5 2.011 0.010 Nueva molécula NO NO SI
NUCB2 nucleobindin 2 2.013 0.003 SI 205 NO SI NO
RGS16 regulator of G-protein signaling 16 2.022 0.000 Nueva molécula NO NO SI
HIST2H3A / HIST2H3C / histone cluster 2, H3a, H3c, H3d 2.024 0.048 Nueva molécula NO NO SI
HIST2H3D
MYOF myoferlin 2.035 0.006 Nueva molécula SI NO SI
C1R complement component 1, r subcomponent 2.048 0.005 Nueva molécula NO NO SI
CCT6A chaperonin containing TCP1, subunit 6A (zeta 1) 2.053 0.002 Nueva molécula NO NO SI
ABCC5 ATP-binding cassette, sub-family C (CFTR/MRP), 2.056 0.013 SI 206 NO SI NO
HLA-A major
member histocompatibility
5 complex, class I, A 2.059 0.029 SI 207 SI NO SI
RPS6KA2 ribosomal protein S6 kinase, 90kDa, polypeptide 2 2.059 0.001 Nueva molécula NO NO SI
EIF5B eukaryotic translation initiation factor 5B 2.059 0.004 Nueva molécula NO NO SI
C1QBP complement component 1, q subcomponent binding 2.076 0.002 Nueva molécula SI NO SI
CTSB cathepsin
protein B 2.078 0.001 SI 208 NO NO SI
IGKC immunoglobulin kappa constant 2.089 0.014 SI 209 SI SI NO
UHMK1 U2AF homology motif (UHM) kinase 1 2.092 0.012 Nueva molécula NO NO SI
NOL7 nucleolar protein 7, 27kDa 2.095 0.001 Nueva molécula NO NO SI
MEGF6 multiple EGF-like-domains 6 2.104 0.000 Nueva molécula NO SI NO
SMNDC1 survival motor neuron domain containing 1 2.104 0.013 Nueva molécula SI NO SI
ZNF259 zinc finger protein 259 2.110 0.022 Nueva molécula SI NO SI
CANX calnexin 2.117 0.027 Nueva molécula SI NO SI
154
Resultados
BASP1 brain abundant, membrane attached signal protein 2.704 0.005 Nueva molécula NO SI NO
BPIFB1 1
BPI fold containing family B, member 1 2.866 0.009 Nueva molécula NO NO SI
DPCR1 diffuse panbronchiolitis critical region 1 1.519 0.024 SI NO SI NO
HLA-DQA1 major histocompatibility complex, class II, DQalpha1 3.182 0.024 Gastroenterology.
SI 219 SI NO SI
HLA-DRB1 /HLA- major histocompatibility complex, class II, DR beta 3.368 0.048
2010
SI 195 SI NO SI
DRB3/ HLA-DRB5 complex Jul;139(1):213-
220
LTF lactotransferrin 3.432 0.019 25.e3.
SI NO SI NO
221
GKN2 gastrokine 2 3.499 0.025 SI SI SI NO
CXCL17 chemokine (C-X-C motif) ligand 17 3.506 0.002 Nueva molécula NO SI NO
GKN1 gastrokine 1 3.834 0.038 SI 222 SI SI NO
HLA-C major histocompatibility complex, class I, C 4.417 0.011 SI 223 NO NO SI
PGC progastricsin (pepsinogen C) 5.004 0.002 SI 224 SI SI NO
a
En el caso de la existencia de más de 1 gen significa que la sonda reconoce varios transcritos. b Fold change es el promedio de expresión en MII-CG dividido
entre MII-No CG. Los genes están ordenados de menor a mayor por dicho valor. c p-valor resultante del t test moderado sin ser ajustado por comparaciones
múltiples. d Genes asociados con MI o CG en estudios de expresión génica, proteíca, de asociación genética, genómicos o funcionales. e Significación estadística
en la comparación MII-CG vs Sanos. f Significación estadística en la comparación MII-No CG vs Sanos. g Se refiere a aquellos genes de interés del diagrama de
Venn de la Figura 32.
156
Resultados
Tabla 33. Genes diferencialmente expresados entre los grupos MIC-CG vs MIC-No CG.
Sign en
Sign en Candidatos
a Fold c Asociados con MIC-No
Símbolo del gen Nombre del gen b p-valor MIC-CG vs de mayor
Change MI o CG d CG vs
Sanos e interés g
Sanos f
REG1B regenerating islet-derived 1 beta 0.421 0.011 SI 214 NO SI NO
GPR110 G protein-coupled receptor 110 0.440 0.015 Nueva molécula NO SI NO
CEACAM20 Carcinoembryonic antigen-related cell adhesion 0.461 0.046 Nueva molécula SI SI NO
molecule 20 similar 225
ANAPC5 anaphase promoting complex subunit 5 0.461 0.031 Nueva molécula NO NO SI
EPB41L3 erythrocyte membrane protein band 4.1-like 3 0.474 0.014 SI 226 NO SI NO
IL1R2 interleukin 1 receptor, type II 0.475 0.017 SI 102 SI NO SI
CXCL17 chemokine (C-X-C motif) ligand 17 0.484 0.038 Nueva molécula SI SI NO
IGHG1 / IGHM / IGHV4-31 immunoglobulin heavy constant gamma 1-4-31 0.493 0.014 SI 209 SI NO SI
PRSS1 protease, serine, 1 (trypsin 1) 0.495 0.033 Nueva molécula NO NO SI
HLA-DRB4 major histocompatibility complex, class II, DR 2.049 0.022 SI 195 NO NO SI
GP2 beta 4
glycoprotein 2 (zymogen granule membrane) 2.063 0.018 Nueva molécula SI NO SI
264
HOXA13 homeobox A13 2.076 0.028 SI SI NO SI
IGFBP5 insulin-like growth factor binding protein 5 2.127 0.004 Nueva molécula SI NO SI
OLFM4 olfactomedin 4 3.332 0.0002 SI 227 SI NO SI
HLA-DRB1 / HLA-DRB3 / major histocompatibility complex, class II, DR 3.475 0.021 SI 195 SI NO SI
HLA-DRB5 beta 1-5
a
En el caso de la existencia de más de 1 gen significa que la sonda reconoce varios transcritos. b Fold change es el promedio de expresión en MIC-CG
dividido entre MIC-No CG. Los genes están ordenados de menor a mayor por dicho valor. c p-valor resultante del t test moderado sin ser ajustado por
comparaciones múltiples. d Genes asociados con MI o CG en estudios de expresión génica, proteíca, de asociación genética, genómicos o funcionales
e
Significación estadística en la comparación MIC-CG vs Sanos. f Significación estadística en la comparación MIC-No CG vs Sanos. g Se refiere a aquellos
genes de interés del diagrama de Venn de la Figura 33.
157
Resultados
Tabla 34. Genes diferencialmente expresados entre los grupos MI-CG vs MI-No CG.
Fold d
Símbolo del gen a Nombre del gen b p-valor c Asociados con MI o CG
Change
HLA-DRB1 major histocompatibility complex, class II, DR beta 1 0.369 0.0052 SI 195
CEACAM20 Carcinoembryonic antigen-related cell adhesion molecule 20 0.457 0.0071 Nueva molécula
similar 225
REG1B regenerating islet-derived 1 beta 0.463 0.0026 SI 214
RBP2 retinol binding protein 2, cellular 0.480 0.0038 SI 200
HSP90AB1 heat shock protein 90kDa alpha (cytosolic), class B member 1 2.028 0.0005 SI 216
HLA-DRB4 major histocompatibility complex, class II, DR beta 4 1.09 0.0014 SI 195
OLFM4 olfactomedin 4 1.10 0.002 SI 227
HLA-DRB1/HLA-DRB3/HLA-DRB5 major histocompatibility complex, class II, DR beta 1-5 1.77 0.0025 SI 195
a b
///
EnLOC100507709 ///
el caso de la existencia de más de 1 gen significa que la sonda reconoce varios transcritos. Fold change es el promedio de expresión
en MI-CG dividido entre MI-No CG. Los genes están ordenados de menor a mayor por dicho valor. c p-valor resultante del t test
LOC100507714
moderado sin ser ajustado por comparaciones múltiples. d Genes asociados con MI o CG en estudios de expresión génica, proteíca, de
asociación genética, genómicos o funcionales.
158
Resultados
159
Resultados
Con el objetivo de validar nuestros resultados con el GSEA respecto a otros antes
publicados, datos de expresión génica por microarray depositados en GEO
(http://www.ncbi.nlm.nih.gov/gds, GSE47797), de metaplasia intestinal y mucosa
gástrica sana fueron comparados 186. En concordancia con nuestros resultados los
principales procesos biológicos sobre-expresados en la metaplasia intestinal
usando el catálogo c2All.v.5 fueron el efecto Warburg no tumoral (N=17),
metabolismo lipídico (N=15), diferenciación intestinal (N=6), ciclo y proliferación
celular (N=5), inflamación (N=5), supresores tumorales (N=4), cáncer gástrico (N=3),
glicosilación proteica aberrante (N=2), apoptosis (N=2), metabolismo de
xenobióticos (N=2) y respuesta al daño genómico (N=2). En el caso de C3All.v5.0
encontramos 6 conjuntos génicos sobre-expresados compuestos por factores
160
Resultados
161
Resultados
Tabla 35. Conjuntos génicos sobre-expresados en MII-CG con los valores más extremos del parámetro Rank at max.
p-valor q-valor Rank at
Conjuntos génicos Proceso funcional ES
nominal a FDR b max c
REACTOME_G_BETA_GAMMA_SIGNALLING_THROUGH_PI3KGAMMA Ciclo y Proliferación celular 0,593 0,007 0,024 895
BERENJENO_TRANSFORMED_BY_RHOA_FOREVER_DN Oncogenes 0,627 0,000 0,002 1908
YU_MYC_TARGETS_UP Oncogenes 0,512 0,003 0,016 2259
SCHOEN_NFKB_SIGNALING Inflamación 0,517 0,002 0,027 2552
PEDERSEN_METASTASIS_BY_ERBB2(her2)_ISOFORM_1 Invasión y metástasis 0,489 0,005 0,027 2804
CROONQUIST_NRAS_SIGNALING_UP Oncogenes 0,589 0,000 0,003 2872
SCIAN_CELL_CYCLE_TARGETS_OF_TP53_AND_TP73_DN Supresores tumorales 0,598 0,004 0,018 2879
SARTIPY_NORMAL_AT_INSULIN_RESISTANCE_UP Genes regulados por insulina 0,646 0,000 0,001 2898
ASTIER_INTEGRIN_SIGNALING Ciclo y Proliferación celular 0,449 0,003 0,031 3343
PEDERSEN_METASTASIS_BY_ERBB2_ISOFORM_4 Invasión y metástasis 0,419 0,000 0,023 3368
REN_BOUND_BY_E2F Ciclo y Proliferación celular 0,571 0,000 0,001 3374
ZHOU_CELL_CYCLE_GENES_IN_IR_RESPONSE_6HR Respuesta a daño genómico 0,423 0,003 0,032 3429
MARKEY_RB1_ACUTE_LOF_UP Supresores tumorales 0,369 0,000 0,037 3592
MARZEC_IL2_SIGNALING_UP Inflamación 0,437 0,000 0,016 3796
PID_INTEGRIN5_PATHWAY Adhesión celular 0,631 0,003 0,018 3850
LU_TUMOR_ANGIOGENESIS_UP Angiogénesis 0,564 0,002 0,034 3855
PHONG_TNF_RESPONSE_VIA_P38_PARTIAL Inflamación 0,366 0,002 0,050 3878
ZHOU_CELL_CYCLE_GENES_IN_IR_RESPONSE_24HR Respuesta a daño genómico 0,429 0,000 0,014 3976
WANG_TUMOR_INVASIVENESS_UP Invasión y metástasis 0,420 0,000 0,002 6767
ST_INTEGRIN_SIGNALING_PATHWAY Adhesión celular 0,398 0,006 0,048 6779
GOTZMANN_EPITHELIAL_TO_MESENCHYMAL_TRANSITION_DN Invasión y metástasis 0,375 0,000 0,026 6789
REACTOME_CELL_CYCLE_MITOTIC Ciclo y Proliferación celular 0,382 0,000 0,011 6864
KEGG_CELL_CYCLE Ciclo y Proliferación celular 0,419 0,000 0,018 6872
162
Resultados
Tabla 36. Conjuntos génicos sobre-expresados en MI-No CG con los valores más extremos del parámetro Rank at max.
p-valor q-valor Rank at
Conjuntos génicos Proceso funcional ES
nominal a FDR b max c
REACTOME_CHYLOMICRON_MEDIATED_LIPID_TRANSPORT Metabolismo lipídico 0,843 0,000 0,157 86
REACTOME_LIPOPROTEIN_METABOLISM Metabolismo lipídico 0,748 0,000 0,123 86
REACTOME_ABCA_TRANSPORTERS_IN_LIPID_HOMEOSTASIS Metabolismo lipídico 0,683 0,009 0,202 167
REACTOME_TERMINATION_OF_O_GLYCAN_BIOSYNTHESIS Glicosilación proteica aberrante 0,716 0,002 0,144 210
WANG_NEOPLASTIC_TRANSFORMATION_BY_CCND1_MYC Oncogenes 0,570 0,002 0,191 284
DAZARD_UV_RESPONSE_CLUSTER_G24 Respuesta a daño genómico 0,671 0,000 0,152 312
LIU_CDX2_TARGETS_UP Diferenciación intestinal 0,738 0,000 0,193 385
REACTOME_LIPID_DIGESTION_MOBILIZATION_AND_TRANSPORT Metabolismo lipídico 0,665 0,000 0,157 396
WOTTON_RUNX_TARGETS_DN Supresores tumorales 0,624 0,008 0,201 437
MOOTHA_GLUCONEOGENESIS Efecto Warburg no tumoral 0,651 0,000 0,168 527
MCMURRAY_TP53_HRAS_COOPERATION_RESPONSE_UP Supresores tumorales 0,533 0,022 0,241 633
GENTILE_UV_HIGH_DOSE_UP Respuesta a daño genómico 0,605 0,021 0,170 661
WANG_BARRETTS_ESOPHAGUS_UP Diferenciación intestinal 0,744 0,000 0,144 829
PID_IL8CXCR1_PATHWAY Inflamación 0,662 0,000 0,140 860
BAUS_TFF2_TARGETS_UP Procesamiento antigénico 0,789 0,000 0,167 952
PID_AMB2_NEUTROPHILS_PATHWAY Inflamación 0,498 0,014 0,241 1044
MURATA_VIRULENCE_OF_H_PILORI H.pylori 0,669 0,002 0,137 1081
KANNAN_TP53_TARGETS_UP Supresores tumorales 0,439 0,019 0,242 1165
KEGG_ETHER_LIPID_METABOLISM Metabolismo lipídico 0,532 0,029 0,246 1240
BROWN_MYELOID_CELL_DEVELOPMENT_UP Inflamación 0,422 0,014 0,240 4270
REACTOME_PHOSPHOLIPID_METABOLISM Metabolismo lipídico 0,425 0,002 0,161 4273
GENTILE_UV_LOW_DOSE_UP Respuesta a daño genómico 0,682 0,000 0,138 4283
AMIT_EGF_RESPONSE_120_HELA Proliferación celular 0,425 0,023 0,241 4300
ZUCCHI_METASTASIS_DN Invasión y metastasis 0,505 0,039 0,236 4466
REACTOME_RIP_MEDIATED_NFKB_ACTIVATION_VIA_DAI Inflamación 0,621 0,032 0,242 4586
REACTOME_TRAF6_MEDIATED_NFKB_ACTIVATION Inflamación 0,569 0,016 0,241 4586
164
Resultados
Tabla 37. Conjuntos génicos compuestos por al menos 3 genes leading edge sobre-expresados en el grupo MII-CG.
Genes leading edge sobre-expresados en Rank at
Conjuntos génicos a Proceso funcional b ES c Nd
la MII-CG e max f
IGLESIAS_E2F_TARGETS_UP Ciclo y proliferación celular 0,418 76 CD53, HSP90AA1, MORF4L1, C3, RBBP7, 5813
LAPTM5
CHICAS_RB1_TARGETS_CONFLUENT Supresores tumorales 0,361 183 CLU, ATP6V0E1, CAV1, EIF5B, MYOF, C1R 4907
GOLDRATH_ANTIGEN_RESPONSE Presentación y procesamiento 0,42 171 HLA-A, C3, CD24, RAN 6476
antigénico
LABBE_WNT3A_TARGETS_UP Ciclo y proliferación celular 0,39 53 RBBP7, CAV1, FABP5, CD24 6333
CHICAS_RB1_TARGETS_GROWING Supresores tumorales 95 CLU, C1R, CAV1, FABP5 4697
ROME_INSULIN_TARGETS_INMUSCLE_UP Genes regulados por insulina 0,417 208 HLA-C, KARS, CCT6A 6051
BERENJENO_TRANSFORMED_BY_ Oncogenes 0,428 157 CAV, CD24, C1R 4980
RHOA_DN
REACTOME_CELL_CYCLE Ciclo y proliferación celular 0,355 149 NHP2, HSP90AA1, RBBP7 5615
SANA_TNF_SIGNALING_DN Inflamación 0,417 32 ATP5A1, CLU, MYOF 5458
MENSSEN_MYC_TARGETS Oncogenes 0,615 30 HSP90AB1, C1QBP, HSP90AA1 5289
HEDENFALK_BREAST_CANCERBRCA1_VS_ Supresores tumorales 0,424 67 PDCD5, CCT6A, FABP5 5652
BRCA2
ALONSO_METASTASIS_UP Invasión y metástasis 0,477 88 NHP2, CAV1, RAN 5250
,
WANG_ESOPHAGUS_CANCER_VS_NORM Cáncer de esófago 0,477 56 HLA-A, C1R, LAPTM5 5578
a
AL_UP
Conjuntos génicos sobre-expresados en la MII-CG ordenados decrecientemente por el número de genes leading edge sobre-expresados en la
MII-CG. b Proceso funcional al que pertenecen los conjuntos génicos. c Enrichment score. d Número de genes leading edge que componen el
conjunto génico. e Genes leading edge que a su ves están sobre-expresados en la MII-CG. f Posición en la lista ranqueada en la cual se obtiene el
mayor valor de ES, los conjuntos génicos más interesantes poseén Rank at max muy pequeño o muy grande.
166
Resultados
Tabla 38. Conjuntos génicos compuestos por al menos 3 genes leading edge sobre-expresados en el grupo MI-No CG.
Genes leading edge sobre-expresados en la Rank at
Conjuntos génicos a Proceso funcional b ES c Nd
MII-CG e max f
BILD_HRAS_ONCOGENIC_SIGNATURE Oncogenes 0,392 54 TRIB1, ELOVL7, AGPAT9, ACOT7, FGFBP1, PI3, 2504
CXCL1, SLC25A37, SERPINB5, TBX3, LYN, CD55,
NT5E, TMEM45B, GK
VECCHI_GASTRIC_CANCER_EARLY_UP Cáncer Gástrico 0,508 133 ACAA2, LAPTM4B, RAD18, CDH3, EPPK1, CASC5, 4205
MACC1, OVOL1, TRIM31, CLDN3, CAMK2N1,
VECCHI_GASTRIC_CANCER_ADVANCED_VS Cáncer Gástrico 0,614 59 CDH17
ATP10B, FAM3D, CLDN15, BTNL8, CFTR, REEP6, 2735
EARLY_DN CLRN3, DHRS11, TM4SF20, GPA33, REG3A
KEGG_DRUG_METABOLISM_OTHER_ Metabolismo de 0,673 12 UGT2B7, NAT2, TPMT, UGT2A3, UGT1A6, UGT1A7, 1348
ENZYMES xenobióticos CDA, CES2, CYP3A4
OHGUCHI_LIVER_HNF4A_TARGETS_DN Diferenciación intestinal 0,605 36 CFI, NLRP6, SULT1B1, SLC30A10, SLC3A1, S100A10, 1356
ETHE1, HSD17B2, MTTP
BAUS_TFF2_TARGETS_UP Procesamiento antigénico 0,789 11 IGHM, RBP2, DGKA, TFF3, ITLN1, ABCG2, ZG16, 952
APOA4
BROWN_MYELOID_CELL_ Inflamación 0,422 52 UGT2B17, HEBP1, CLDN1, CCND2, CEACAM1, XDH, 4270
DEVELOPMENT_UP GLRX, ANPEP
LABBE_TARGETS_OF_TGFB1_AND_WNT3A Proliferación celular 0,450 30 TPK1, TFRC, LGALS4, SLC30A4, HNF4A, VIL1, VNN1, 2919
_DN SERPINA1
ALCALA_APOPTOSIS Apoptosis 0,524 30 GPI, MTCH2,CASP1, HADHA, QPRT, ATP1A1 3103
KEGG_ARGININE_AND_PROLINE_ Efecto Warburg no 0,605 16 NAGS, MAOB, PRODH, AGMAT, OAT, ABP1 1415
METABOLISM tumoral
SERVITJA_ISLET_HNF1A_TARGETS_DN Diferenciación intestinal 0,515 35 NR1H4, CDHR2, RNF186, MTMR11, TM4SF4, ACE2 1707
REACTOME_AMINO_ACID_AND_ Efecto Warburg no 0,618 13 SLC15A2, SLC6A20, SLC1A1, SLC15A1, SLC7A9, 1480
OLIGOPEPTIDE_SLC_TRANSPORTERS tumoral SLC6A19
167
Resultados
Figura 35. Análisis de genes leading edge realizado con todos los conjuntos génicos
sobre-expresados en MII-CG. La acumulación de color verde indica que
los conjuntos génicos comparten genes leading edge.
169
Resultados
Figura 37. Heat map de los genes leading edge del núcleo 1 de la figura 35. Los
colores rosado y rojo indican sobre-expresión génica.
Figura 38. Subanálisis de genes leading edge del núcleo 2 (fosforilación oxidativa)
de la figura 35.
170
Resultados
171
Resultados
Figura 39. Esquemas de las vías canónicas presentación y procesamiento antigénico sobre-expresada en los grupos MII-CG (A) y MIC-CG
(B) obtenidas con el programa Ingenuity Pathway Analysis. Las moléculas en colores rojo y verde se encuentran sobre y sub-
expresadas, respectivamente.
173
Resultados 174
Resultados
Figura 40. Redes de alto score IPA de las comparaciones MII-CG vs MII-No CG (A-B) y MI-No CG vs Sanos (C). La intensidad del color en los
nodos indica el grado de expresión: verde sub-expresado, rojo sobre-expresado y sin color no tiene diferencia de expresión. La
longitud de las líneas discontinuas refleja la evidencia en la literatura en la relación entre moléculas. Las líneas no discontinuas
indican interacción proteína-proteína. Las flechas indican la dirección de asociación.
175
Resultados
La red molecular con mayor score (44) en el grupo MI-No CG contiene a varias
moléculas del metabolismo lipídico como APOA4, APOB, APOC3, APOA1, FABP1/2
MTTP y HNF4A, sobre-expresadas en la MI-No CG y que son componentes de los
conjuntos génicos significativos del GSEA (Tabla 40, Figura 40C). El factor de
transcripción de hepatocitos HNF4A, con un rol en el desarrollo intestinal, esta
sobre-expresado en los datos e interactúa con varias lipoproteínas, siendo un nodo
de esta red molecular y quizás responsable del incremento del metabolismo lipídico
anteriormente mencionado.
176
Resultados
Tabla 40. Redes moleculares obtenidas en el IPA para las comparaciones MII-CG vs MII-No CG y MI-No CG vs Sanos.
Comúnes
Redes moleculares Comparación Score Moléculas centrales (Fold Change) a
con GSEA b
Enfermedades del tejido MII-CG vs MII- 39 Sobre-expresadas: HLA-DQB1, IL1R2, HLA-C, HLA-A, CANX, IGHM, CD24, HLA-A, CD24,
conectivo, enfermedades No CG CAV1, HSP90AA1, C3, CTSB, UHMK1, CLU, PDCD5, LTF, FABP, HLA-DRB1, C3, LTF
inflamatorias, esqueléticas y RHOA
musculares Sub-expresadas: MME, HLA-DRB1, RHOA, GNAS, ATP1B3
Señalización e interacción MII-CG vs MII- 16 Sobre-expresadas: CXCL14, MORF4L1, PI3, RGS16, HLA-DQA1, HLA-DRB4, NUCB2
celular, desarrollo del No CG CRISPLD2, NUCB2, RPS6KA2
sistema hematológico, Sub-expresadas: RBBP7, GSTA1
respuesta inmune humoral
Muerte y sobrevivencia MII-CG vs MII- 16 Sobre-expresadas: CCDC144B, C1R, HLA-DRB3, IGHG1, ANP32B, GNL3, C1R, ANP32B
celular, enfermedad No CG LYZ
inmunológica y metabólica Sub-expresadas: RBP2, DNAJB7, CP, GIP
Cáncer, daño orgánico, MIC-CG vs 5 No son redes de alto score (>15)
muerte y sobrevivencia MIC-No CG
Respuesta
celular anti-microbiana, MIC-CG vs 3 No son redes de alto score (>15)
inflamatoria, señalización e MIC-No CG
interacción celular
Respuesta inflamatoria, MIC-CG vs 3 No son redes de alto score (>15)
desarrollo tisular y MIC-No CG
morfología celular
177
Resultados
Comúnes
Redes Comparación Score Moléculas centrales a
con GSEA b
Metabolismo lipídico, MI-No CG vs 44 Sobre-expresadas: ABCC2, ABCG5, ABCG8, APOA1, APOA4, APOB, APOC3, ABCC2,
Transporte molecular, Sanos CDX1, CREB3L3, CYP3A4, DGKA, DMBT1, F3, FABP1, FABP2, FUT6, HNF4A, APOA4,
Bioquímica de moléculas
HNF4G, MME, MTTP, MUC4, NMUR2, NR1H4, NR1I2, RBP2, CDX1,
pequeñas
SLC2A5,SLC52A1,TNFRSF11A,TPMT FABP1,
FABP2,
Sub-expresadas: FOXA2, TIMP3 HNF4A,
MUC4,
Muerte y sobrevivencia MI-No CG vs 41 Sobre-expresadas: ACHE, ACSL5, ANXA2, CASP1, CD55, CEACAM1, CIDEB, CASP1,
TPMT, CD55,
cellular, Desarrollo y function Sanos CIDEC, CYCS, CYP4F3, EPHX2, GPRC5A, HMOX1, IL32, IL2RG,PRODH, MDK, CEACAM1,
CYP3A4,
del sistema digestivo,
PFKFB2, PIGR, S100A10, SERPINB5, TFF3, TRIB1, VDR TRIB1,
NR1H4,
Morfología de organos
Sub-expresadas: CCKAR, CSTA, MAL, CXCL1, DGKD, GAST CXCL1,
DGKA, RBP2,
EPHX2,
SLC2A5
PRODH,
Cáncer, daño y MI-No CG vs 37 Sobre-expresadas: PSMD1, ARL4A, ATP1B3, ATP7B, CCND2, CDKN1A, SERPINB5,
----------
anormalidades a nivel de Sanos CLDN1, COL17A1, DEFA5, EPCAM, GCG, HSD17B2, IgG e IgM, LGALS3, S100A10,
organism, Ciclo celular
LYN, MSLN, NT5E, SERPINA1, SLAMF7, SLC39A4, TFRC, TOP1 TFF3
Tabla 41. Resultados del análisis de ontología génica en las comparaciones de estudio.
MI-No CG vs Sanos Metabolismo de: flavonoides (p=4.82E-07), UDP-glucoronosil transferasa (p=1.32E- Complejo MHC (p=2.35E-08), región
glucuronidación (p=4.35E-09), ácido 07),unión a ácidos nucleicos (p=4.64E-05), extracelular (p=7.77E-04), complejos proteicos
urónico (p=4.35E-09), glucuronato UDP-glicosil transferasa (p=6.97E-04), (p=6.38E-03) y macromoleculares (p=6.38E-03),
(p=4.35E-09), xenobióticos (p=5.17E-08), unión de iones(p=1.24E-06), membranas del retículo endoplásmático y
monosacáridos (p=2.88E-07), carbohidratos transferencia de grupos glicosilo nuclear externa (p=1.82E-03), complejos
(p=2.98E-06), moléculas pequeñas (p=7.41E-03), actividad sulfo transferasa proteicos (p=6.38E-03) y macromoleculares
(p=1.55E-04). Homeostasis del colesterol (p=1.17E-02), unión de ácidos carboxílicos (p=3.82E-03), matriz extracelular (p=4.88E-02)
(p=2.03E-02), esterol (p=2.03E-02) y lípidos (p=2.01E-03), transportadores de
(p=2.59E-02). Metabolismo de lípidos colesterol y esterol (p=2.01E-03), unión a
(p=6.83E-08) y ácidos grasos (p=1.55E-02) lípidos (p=1.13E-02)
179
Resultados
180
Resultados
Tabla 42. Resultados de prueba ANOVA de los genes de referencia en los grupos
MI-CG, MI-No CG y controles sanos.
Gen de
Comparación ANOVA Observaciones
referencia
RPL29 MI-CG vs MI-No CG < 0.0001
RPL29 MI-CG vs MI-No CG vs Sanos < 0.0001
B2M MI-CG vs MI-No CG < 0.0001
B2M MI-CG vs MI-No CG vs Sanos < 0.0001
G6PD Casos MI-CG vs Sanos 0.2023 Falló en MI-No CG
GAPDH Casos MI-CG vs MI-No CG 0.6617 Falló en controles sanos
ACTB MI-CG vs MI-No CG 0.4478 Gen escogido
ACTB MI-CG vs MI-No CG vs Sanos 0.3946 Gen escogido
Ct de ACTB Casos y Controles
40
30
20
10
0
M oC 1
M o 2
n 4
I- 1
I- 2
I-C 3
I- 4
I-C 5
M MI- G6
M No G7
o G3
no 1
s2
G
I-N G
Sa CG
M G
M CG
M CG
M G
M CG
Sa os
I-N C
I-N C
I-C
I- C
M
Figura 42. Gráfico de cajas de los valores de Ct del gen ACTB en las muestras
analizadas de los grupos MI-CG, MI-No CG y Sanos.
181
Resultados
Sanos en los resultados de qRT-PCR (Tabla 43 y Figura 43). Sin embargo en solo 21
de ellos fue validada el mismo sentido de la diferencia de expresión, para un 56.75%
de efectividad. Los genes GSTA1, IK, PPIA, GNAS del grupo MII-CG se encuentran
sub-expresados en el microarray mientras que están sobre-expresados en la qRT-
PCR. Los genes IGHG4, AFAP1-AS1, ATP5A1, RHOA, ATP6V0E1, HLA-DRB4, PGC, HLA-
DQA1, OLFM4, HLA-DRB4 estan sobre-expresados en MII-CG y MIC-CG mediante
ambas tecnologías; pero tuvieron un mayor Fold Change en la qRT-PCR que en el
microarray. Respecto al grupo MI-No CG la mayoría de los genes (TMEM25, GSTA1,
AFAP-AS1, MUC3A, CLCA1, OLFM4, ANPEP, FABP1) se encuentran sobre-expresados
mediante ambas tecnologías. Las excepciones son CDH17 y CPS1, sobre-expresados
en el microarray y sub-expresados en la qRT-PCR.
Figura 43. Gráfica del nivel de expresión del microarray y la qRT-PCR de los genes
significativos por ambas metodologías.
182
Resultados
Tabla 43. Resultados del test de permutaciones del programa Bootsratio empleado
para estimar diferencias significativas en grupos MII-CG, MIC-CG y MI-No CG.
Genes Tamaño de Fold Change en Fold Change en Grupo donde es p-valor en qRT- Validado
a b c d e
significativos muestra qRT-PCR microarray significativo PCR
del microarray
IGHG4 11/7 12.82 2.67 MII-CG p<0.0005 SI
GSTA1 13/7 20.18 0.39 MII-CG 0.0015 NO
IK 13/6 28.12 0.38 MII-CG 0.006 NO
PPIA 12/7 30.54 0.39 MII-CG p<0.0005 NO
GNAS 11/7 31.56 0.45 MII-CG p<0.0005 NO
SLFN5 6/4 33.46 No sign MII-CG 5.00E-04 SI
AFAP1-AS1 12/7 41.53 2.16 MII-CG p<0.0005 SI
ATP5A1 12/7 46.39 1.3 MII-CG 5.00E-04 SI
RHOA 13/7 55.91 2.16/0.4 MII-CG 0.001 SI
ATP6V0E1 10/6 57.18 2.49 MII-CG 0.0075 SI
HLA-DRB4 10/5 67.13 2.25 MII-CG 0.031 SI
PGC 13/8 67.59 5.01 MII-CG 0.003 SI
HLA-DQA1 50/3 288.90 3.18 MII-CG 0.0345 SI
OLFM4 22/29 33,30 3.33 MIC-CG 5.00E-04 SI
HLA-DRB4 14/14 68,81 2.05 MIC-CG 0.0335 SI
PGC 39/9 0.65 0.13 MI-No CG 0.0075 SI
TMEM25 6/2 15.8 2.78 MI-No CG p<0.0005 SI
GSTA1 26/4 17.35 3.81 MI-No CG p<0.0005 SI
AFAP1-AS1 33/7 35.67 4.52 MI-No CG p<0.0005 SI
TFF3 38/8 7.18 9.08 MI-No CG 0.009 SI
MUC3A 27/6 30.68 15.9 MI-No CG p<0.0005 SI
CDH17 30/3 0.25 18.02 MI-No CG p<0.0005 NO
CPS1 32/6 0.37 26.25 MI-No CG p<0.0005 NO
CLCA1 34/2 55.26 33.77 MI-No CG p<0.0005 SI
OLFM4 36/3 56.36 37.04 MI-No CG p<0.0005 SI
ANPEP 31/4 62.3 49.98 MI-No CG p<0.0005 SI
FABP1 25/2 70.85 80.58 MI-No CG p<0.0005 SI
a
Número de muestras analizadas en los grupos de casos (MI-CG)/controles (MI-No CG).
b
Media de valores de expresión de casos (MI-No CG) dividido entre controles (Sanos) en la
qRT-PCR. c Media de expresión de casos dividido entre controles en el microarray. d Se refiere a
el microarray. e p-valor del test de permutaciones que evalua la significación de la diferencia de
expresión entre casos y controles en la qRT-PCR.
183
Resultados
Cuando fueron incubadas estas líneas con H. pylori solo se observó un incremento
en la expresión de Schlafen5 en la HL-60, mientras que en las Jurkat se observó un
decremento de la expresión, aunque ambos p-valores son nominales.
184
Resultados
Figura 44. Inducción de Schlafen5 por el IFNα en las líneas celulares HL-60 y Jurkat.
185
Resultados
186
Resultados
Figura 46. Inmuno-histoquímica de Schlafen5. A) Mucosa gástrica normal. B) Control negativo. C) Control positivo (nasofaringe).
D) Control positivo (íleon) con enfermedad de Crohn. E) Validación usando el anticuerpo HPA 017760 (Sigma-Aldrich, EUA).
Imágenes a 200X.
187
Resultados 188
Resultados
Figura 47. Inmuno-histoquímica de la proteína Schlafen5. A) Gastritis no atrófica. B) Gastritis crónica atrófica. C) MIC que no progresa a
CG. D) MIC que progresa CG. E) MII que no progresa CG. F) MII que progresa CG. G) CG tipo intestinal. H) CG tipo difuso.
Imágenes a 200X.
189
Resultados
Tabla 45. Cuantificación de la tinción estromal y glandular de Schlafen5 mediante el score H en los grupos analizados.
Score H
Score H
Infección por estromal de
Diagnóstico b Edad c Sexo (% Score H estromal de Schlafen5 categorizado glandular de
N H. pylori Schlafen5 d
Histológico a (% ± DE) hombres) N (media) e Schlafen5
(% +) (media ±
(media ± DE)
mediana)
Negativo Débil + Moderado +
Normal 41 42.5 (18.8) 51.2 0 3.87 (1.25) 37 (90.24) 4 (9.76) 0 (0) 0.89 (0.00)
GNA 26 44.2 (12.3) 46.1 88.5 5.84 (2.53) 21 (80.77) 5 (19.23) 0 (0) 8.54 (0.42)
GCA 49 48.5 (12.6) 38.8 49.0 4.56 (0.84) 40 (81.63) 9 (18.37) 0 (0) 5.15 (0.00)
MIC 44 53.1 (10.5) 45.4 63.6 19.73 (5.73) 25 (56.82) 18 (40.91) 1 (2.27) 14.78 (5.82)
MIC-CG 6 62.0 (6.6) 50 33.3 49.33 (50.68) 2 (33.33) 4 (66.67) 0 (0) 23.23 (14.82)
MIC- No CG 38 51.7 (10.4) 44.7 68.4 15.05 (2.53) 23 (60.53) 14 (36.84) 1 (2.63) 13.41 (5.82)
MII 54 58.3 (12.8) 50 51.8 20.95 (5.11) 32 (59.26) 20 (37.04) 2 (3.7) 7.69 (0.00)
MII-CG 19 66.5 (13.6) 63.2 42.1 45.94 (40.00) 3 (15.79) 14 (73.68) 2 (10.53) 7.99 (1.00)
MII-No CG 35 53.8 (9.9) 42.9 57.1 7.39 (1.14) 29 (82.86) 6 (17.14) 0 (0) 7.52 (0.00)
MI 98 55.9 (12.1) 48 57.1 20.34 (31.8) 57 (58.04) 38 (38.97) 3 (2.98) 10.7 (1.7)
MI-CG 25 65.4 (12.4) 60 40 47.63 (40.18) 5 (24.56) 18 (70.17) 2 (5.26) 11.6 (1.48)
MI-No CG 73 52.7 (10.1) 43.8 63.0 11.22 (2.00) 52 (71.69) 20 (26.99) 1 (1.31) 10.4 (1.00)
CG 67 72.6 (11.0) 61.2 25.4 7.9 (1.38) 51 (76.12) 16 (23.88) 0 (0) 3.22 (0.00)
a
Diagnóstico histológico al reclutamiento. b Tamaño de muestra. c Edad del paciente en el diagnóstico al reclutamiento. d Score H calculado según
%Células+ con intensidad Baja x 1 + %Células+ con intensidad Media x 2 + %Células+ con intensidad Alta x 3 ± mediana. e Categorización del score
H según 10>scoreH≥0 (Negativo), 100≥scoreH≥10 (Débil Positivo), 200≥scoreH≥100 (Moderado Positivo) y 300≥scoreH≥200 (Fuerte Positivo) 191.
190
Tinción estromal de Schlafen5
N=25 N=19
150
N=73
N=6
100
N=67
50
N=41 N=26 N=49
r G CG
I-N G
A
A
S
M G
o
O
-C
C
I-C
ic
G
G
N
II-
tr
o
IC
M
SA
as
M
ce
án
C
Figura 48. Representación gráfica de la tinción estromal de Schlafen5. Corresponde
a los resultados de la tabla 45.
191
significativa la diferencia entre ambos modelos logísticos por el test de razón de
verosimilitud (p<0.0001) (Figura 49).
Figura 49. Curvas de sobrevivencia ROC para los modelos MI-CG con adición de
Score H estromal de Schlafen5, AUC=91.2% (--------) y MI-CG sin adición de
Score H estromal de Schlafen5, AUC=78.3% ( ).
192
Figura 50. Curvas de sobrevivencia ROC para los modelos MII-CG con adición de
Score H estromal de Schlafen5, AUC=95.3% (--------) y MII-CG sin adición de
Score H estromal de Schlafen5, AUC=75.6% ( ).
193
Figura 51. Doble marcaje de Schlafen5 con marcadores de superficie. A-B) CD20
células B. C-D) CD2 células T. E-F) MAC2 células mieloides. Schlafen5 marca
los núcleos con color pardo, CD20 y CD2 con rosado y Mac2 con azul. A la
izquierda magnificación 200X y a la derecha 400X.
194
DISCUSIÓN
Discusión
Las asociaciones más significativas del estudio se observan con los SNPs de NOD2
rs7202124 y rs2111235, los cuales manifiestan asociaciones inversas con el riesgo
de CG no cardias aún después de aplicar la corrección para comparaciones
múltiples. La variante rs7202124 es un eQTL en cis que controla el nivel de
expresión de NOD2 en monocitos por lo que es posible que dicha asociación
genética refleje este hecho. El mapeo de eQTLs es una técnica estadística que
correlaciona los genotipos de polimorfismos genéticos con la expresión de ARNm
con el objetivo de localizar variantes genéticas que regulan la expresión de cada
transcrito 228.
197
Discusión
que permiten confiar en la validez de dicha asociación. A pesar de que estos mismos
SNPs de NOD2 no fueron asociados con el subtipo intestinal despúes de aplicar el
test de corrección para comparaciones múltiples y el de heterogeneidad, el mismo
haplotipo GTCGC (frecuencia 0.141) también se asocia negativa y nominalmente
(p=0.031) con el CG intestinal. Lo cual indica la tendencia a la asociación de estos
SNPs y el análisis haplotípico con mayor potencia es capaz de detectar la asociación.
198
Discusión
Los SNPs genotipados fueron seleccionados por ser tagSNPs, es decir SNPs
marcadores de bloques de haplotipos en desequilibrio de ligamiento con otros SNPs
que pudieran tener un efecto funcional, por lo que es posible que dichos tagSNPs
no tengan un efecto funcional en sí mismo. Algunas maneras de saberlo son
mediante estudios de mapeo fino o funcionales, que no fueron realizados en el
presente trabajo.
199
Discusión
mecanismo NMD 235. El programa Pupasuite no reconoció los alelos de dichos SNPs
como que muestren un reconocimiento diferencial por enhancers (ESE) o
silenciadores (ESS) del splicing alternativo, fenómeno que puede crear un exón en el
interior de un intrón propiciando la degradación NMD, como se ha reportado en
diabetes mellitus 236. Basado en lo cual, no se puede concluir una función de estos
SNPs relacionada con el mecanismo NMD.
Los estudios publicados de CD14 se han centrado en el SNP rs2569190 (-260 C/T),
que se encuentra localizado colindante al sitio de unión al factor de transcripción
Sp1, el cual influye en la expresión de CD14 240. A diferencia de estos resultados,
estudios en poblaciones de Polonia, EUA 105 y Taiwán 241 no encontraron asociación
de esta variante con el CG no cardias, un estudio en Japón encontró una asociación
inversa con el subtipo intestinal 242 y un meta-análisis solo encuentra asociación
para el SNP rs2569190 (-260C/T) en pacientes infectados por H. pylori 106. Por lo que
sabemos, no se ha realizado ningún estudio con respecto a la asociación de CD14
con GC cardias.
200
Discusión
En este sentido, la búsqueda de eQTLs entre todos los SNPs encontró que las
variantes asociadas al CG de la región de CD14 rs11167532, rs1583005 , rs2569190,
y rs778588 han sido encontradas previamente como reguladores de la expresión
génica de los genes SRA1 156 y WDR55 157 en monocitos. Este tipo celular se
acumula en el estroma gástrico a consecuencia de la inflamación generada por
Helicobacter pylori, transformándose a macrófagos que incrementan la inflamación
por un aumento de la generación de radicales libres 27. Respecto a la variante
rs2569190 (-260T) se ha descrito experimentalmente que regula los niveles de CD14
presentes en la membrana de monocitos sanguíneos, estando asociado los
homocigotos del alelo de riesgo con el incremento de la expresión de CD14 245, lo
cual está en concordancia con que sea identificado como un eQTL en estudios de
este tipo. Existen múltiples casos en los que tagSNPs de estudios GWAS asociados
con cáncer son eQTLs como en cáncer de próstata 244 y colorectal 245. Por lo que se
puede sugerir que la asociación identificada en CD14 (rs2569190) sea mediada por
un control de la expresión de CD14 en monocitos ejercido por este SNP, para el
resto de las variantes es necesario hacer análisis o experimentos in vitro adicionales
que permitan aclarar dicho efecto.
Debido a que la infección por H. pylori no se asocia con CG del cardias 230 y, sin
embargo, CD14 es receptor del lipopolisacárido de bacterias gram-negativas como
H. pylori, es posible que otros mecanismos puedan explicar la asociación de CD14
con el CG del cardias. En este sentido, CD14 podría actuar como receptor para el
virus de Epstein-Barr, una infección que sí está asociada con el CG del cardias 246, ya
que se ha planteado la posible existencia de un nuevo receptor diferente a CR2 , el
receptor del virus de Epstein-Barr 247. CD14 es un candidato por ser un receptor
poli-específico que interactúa con el dipéptido muramilo del LPS, peptidoglicano,
lipoarabinomannan y fosfolípidos 248. Lamentablemente, no se dispone de datos
sobre la infección por el virus de Epstein-Barr en la población estudiada. Una
evidencia que pudiera apoyar la anterior aseveración es el hecho que CD14 no se
asocia en nuestra población con CG no cardias, localización anatómica fuertemente
asociada con la infección por H. pylori.
201
Discusión
TLR4 también pertenece a la familia de los PRRs y junto a CD14 y MD-2 participa en
el reconocimiento del LPS de H. pylori. Después de su activación las proteínas
adaptadoras MyD88, TIRAP, TRAM y TRIF se unen a los dominios TIR (toll-
interleukin (IL)-1 receptor) del TLR4 y activan las quinasas IRAK, MAPK y TRAF
induciendo la activación de los factores transcripcionales NFKB1, AP-1 e IRF. En
relación con la falta de asociación a nivel de SNPs individuales de TLR4 obtenida en
este estudio es interesante apuntar que la infección de la línea celular HEK293 con
H. pylori es capaz de inducir NFKB en células transfectadas con TLR2 o TLR5 pero no
por TLR4 251 y que un estudio afirma que el reconocimiento de H. pylori por el
epitelio gástrico es independiente de TLR4 252.
Respecto a NFKB1 existió una tendencia a la asociación ya que los SNPs rs765789,
así como, rs228611 y rs1585214 se encontraron positiva y negativamente asociados
con CG no cardias y con el subtipo difuso; aunque ninguna de estas asociaciones fue
significativa después del test de comparaciones múltiples. Sin embargo, el haplotipo
TGAGGTGACCAG se encuentra positivamente asociado con el CG no cardias. Esta
combinación está formada por el alelo de riesgo A de rs765789 y los alelos salvajes
202
Discusión
203
Discusión
Solo el 17.2% de los pacientes progresan a una LPG de mayor severidad entre el
reclutamiento y el final de seguimiento, lo cual explicaría que la mayoría de las
asociaciones genéticas significativas fueran obtenidas al comparar el grupo de la
regresión con el resto.
204
Discusión
misma lesión (46.3%), seguidos de los que regresaron (36.5%) en comparación con
sólo un 17.2% de progresión.
Las principales ventajas de este estudio son que los pacientes pertenecen a distintas
regiones españolas, su largo tiempo de seguimiento (media 12 años) y el relativo
elevado tamaño de muestra. Su principal limitante es que la toma de biopsias fue
más rigurosa al final del seguimiento que al reclutamiento, debido al cumplimiento
del protocolo de Sydney de muestreo sólo al final del seguimiento. Sin embargo,
esto fue compensado solicitando que cada biopsia al reclutamiento tuviera al
menos 3 fragmentos con representación de cuerpo y antro. Adicionalmente,
algunos tagSNPs no pudieron ser genotipados reduciendose la cobertura de la
variabilidad de los genes donde fueron seleccionados.
Mediante este estudio se ha pretendido por una parte replicar SNPs previamente
asociados con la evolución de las LPGs en el estudio de Soria y, por otra, estudiar si
la variabilidad en nuevos genes candidatos se asocia estadísticamente con la
evolución en el tiempo de las lesiones precursoras gástricas. Los genes estudiados
realizan funciones claves en la fisiología gastrointestinal como se evidencia en los
modelos murinos de carcinogénesis gástrica, que se caracterizan por el hecho de
que al cabo de un tiempo aparece alguna o varias de las LPGs seguidas del
desarrollo de CG.
Como resultado de este estudio se asociaron SNPs con un incremento del riesgo de
progresión en los genes CD14, DNAH11, MAPK3, PRKAA1, RUNX3 y TFF1. Otras
variantes localizadas en los genes CD14, CDH1, PTGS2, TFF1 y TFF2 se asociaron
inversamente con la progresión de las lesiones (Tabla 25 A). El análisis de la
regresión mostró que SNPs en los genes CD14, IL1B, IL1RN, NFKB1, TFF1 y TFF2 se
asociaban con una mayor regresión, mientras que aquellos en CDH1, MUC1, MUC2
y NFKB1 se asociaron con una menor regresión (Tabla 25 B). El SNP más significativo
y el único que se mantiene así despúes de aplicar la corrección para comparaciones
múltiples es rs10902073 de MUC2, el cual es un tagSNP localizado en una región
205
Discusión
El análisis bioinformático mostró que SNPs asociados con la evolución de las LPGs se
encuentran localizados en zonas regulatorias, UTR 5’, UTR 3’, o son eQTLs
reguladores de la expresión génica. Aunque hay que tener presente que son
predicciones funcionales y que la mayoría de los SNPs se seleccionaron por ser
tagSNPs, por lo que su funcionalidad sería, en todo caso, casual. Sin embargo, es
relevante destacar que CD14 es uno de los genes con mayor número de SNPs
asociados con la evolución de LPGs, los cuales a su vez fueron encontrados como
eQTLs reguladores de la expresión génica, pero los mismos controlan los niveles de
expresión de WDR55 en lugar de CD14. Además, no existen diferencias significativas
206
Discusión
Uno de los trabajos antes publicados se realizó con muestras del estudio de
seguimiento de Soria y se evaluaron SNPs en los genes MUC1, MUC2 y MUC6 140.
Los SNPs (rs10794293, rs3924453, rs4077759) en el extremo 3’ de MUC2 se
asociaron con un menor riesgo de progresión de las lesiones. Mientras que los
localizados en el extremo 5’ se asociaron con el aumento (rs2071174) o disminución
(rs10902073, rs10794281, rs7944723) de la probabilidad de regresión, indicando un
rol de este gen en evitar la progresión de las LPGs a estadios de mayor severidad. El
análisis haplotípico fue el más relevante, ya que ambas combinaciones CTCCCG y
CATAGAAC de la región 5’ y 3’ se asociaron positivamente con la regresión e
inversamente con la progresión, respectivamente 140.
207
Discusión
Figura 52. Bloques de haplotipos de MUC2, el bloque que comprende a los SNPs
rs10902073 y rs10794281 (marcados en negro) más significativos del
estudio que se encuentran en el extremo 5’ del gen.
Otros trabajos han encontrado que los genes COX2 -765C 124, las citoquinas pro-
inflamatorias IL8 -251, MIF -173 128 y XRCC1 y OGG1 134, del mecanismo de
reparación del ADN, se asocian con la evolución de las LPGs. La mayoría de estos
208
Discusión
También, los SNPs rs1583005 (OR 1.56, IC95% 1.13-2.16) y rs574445 (OR 1.63,
IC95% 1.13-2.36) de CD14 se asociaron con una mayor progresión de las LPGs. Lo
cual fue corroborado en el haplotipo TGTG que se asocia igualmente con una mayor
209
Discusión
El alelo T del SNP rs4503658 de CDX2 presenta una asociación nominal (OR 1.54,
IC95% 1.01-2.37, p=0.047) y positiva con la regresión, sin embargo el haplotipo CAT
que contiene a dicho alelo presenta un OR mayor (1.78, IC95% 1.11-2.85, p=0.0168)
210
Discusión
que la variante asociada. CDX2 es el agente causal más importante o master gene
de la metaplasia intestinal 52, ya que induce la expresión de los factores de
reprogramación de desdiferenciación celular SALL4 y KLF5 63 y genes de
diferenciación intestinal como la sacarasa isomaltasa, cadherina LI, MUC2, FURIN,
entre otros 64. Además, es un inhibidor de la reparación del ADN no homóloga en
cáncer de colón 259, por lo que pudiera tener un rol similar en estómago.
MAP3K14 es una quinasa de tipo serina/treonina que forma parte de las vías de
señalización ERK y JNK que activan a NFKB 170. El alelo de riesgo C del SNP
rs16939926 de este gen tiene bajo tamaño de muestra, el haplotipo CCTC muestra
un efecto protector de la regresión (OR 0.52, IC95% 0.28-0.95) pero con baja
frecuencia (0.039) y una significación no robusta por lo que no reviste importancia.
211
Discusión
Mediante el análisis del perfil de expresión génica por microarrays usando el array
Almac-Xcel (Affymetrix, EUA) específico para muestras parafinadas, hemos
encontrado 106 genes diferencialmente expresados en la MII que progresa a CG
respecto a la que no progresa. Aunque las diferencias de expresión no son muy
grandes, algunos genes poseen funciones en la carcinogénesis que pudieran influir
en la progresión a CG desde la MII. Los candidatos de mayor interés son transcritos
significativos en MII-CG vs MII-No CG (N=36) y en la comparación MII-CG vs Sanos
(N=22), excluyendo transcritos no específicos compartidos con la comparación MII-
No CG vs Sanos (Figura 32, Tabla 32).
212
Discusión
213
Discusión
214
Discusión
215
Discusión
sección Revisión Bibliográfica) como CDX2, ACE2, KRT20, MUC13, OLFM4, REG4 81,
FABP1, MEP1B, SI, SLC6A19 214, CDX1, MTTP, CEACAM6 79, DGKQ 78, APOA1, APOA4,
APOB, ALDOB, CDH17, CLDN3, CLDN4, GGT1, HNF4A, VIL1 80, MUC2 266, MUC4 267,
TFF3 268 y CPS1 269 (Tabla anexa 2). La amplia similitud encontrada con genes
previamente hallados desregulados en la MI podemos considerarla como una
validación de la veracidad de estos resultados.
CDX1 y CDX2 presentan Fold Changes no tan altos, de 4.5 y 3.28 respectivamente,
por lo que es posible que su nivel de expresión no necesite ser tan elevado para
seguir propiciando la diferenciación intestinal en el estómago una vez establecida la
misma. De manera muy interesante encontramos la sobre-expresión de otro
miembro de la familia génica homeobox, HOXA13 (Fold Change 3.62, p-valor
ajustado 9,0725E-07), el cual no ha sido reportado previamente desregulado en la
MI, participa en el desarrollo embrionario del eje antero-posterior. Así como, la sub-
expresión de HOTAIRM1 (Fold Change 0.45, p-valor ajustado 0.014), un ARN
antisentido que regula negativamente a esta familia génica y que ha sido asociado
con mala sobrevivencia en CG 270. Por lo que es posible que HOXA13 contribuya
sinergísticamente con CDX1/2 a la diferenciación intestinal. La sobre-expresión de
este gen también fue observada en la MIC-CG, resaltando su posible funcionalidad.
216
Discusión
Otros grupos de transcritos con alta sobre-expresión (Fold Change 2.12-4.73) son
seis miembros de la familia de las sulfotransferasas (SULT1A1-4, SULT1E1,
SULT1B1). Enzimas que participan en la conjugación de grupos sulfato a drogas y
compuestos xenobióticos y carcinógenos. Además, existen varios genes de la familia
de las UDP-glucoronosiltransferasas (UGTs) sobre-expresados (Fold Change 2.01-
3.63) los cuales participan en la modificación química y eliminación de compuestos
xenobióticos potencialmente tóxicos. Previamente se ha encontrado un aumento
de la actividad sulfotransferasa en el cáncer gástrico 272. También, se encuentran
sobre-expresadas varias proteínas que contienen motivos tripartitas
(TRIM14/15/31/36/40, Fold Change 2.01-3.05), las cuales activan a NFKB
incrementando la respuesta inflamatoria 273.
217
Discusión
Diversos miembros (N=17) de los ARNs nucleolares pequeños C/D box (SNORDs116)
se encuentran sub-expresados. Los mismos son reguladores de la metilación de
ARNs ribosomales influyendo esta modificación química en la estructura
tridimensional de los mismos y su interacción con proteínas ribosomales.
Microdeleciones en el cluster de genes SNORDs116 contribuyen al síndrome de
Prader-Willi y este cluster controla la expresión génica de alrededor de 200 genes
274
. Otras proteínas gástricas reguladas negativamente son las mucinas MUC1,
MUC5AC, MUC5B y MUC6, cuya disminución de la expresión ha sido reportada
previamente en la MIC 14.
218
Discusión
Existen diversos métodos para este análisis pero no brindan soluciones estadísticas
exactas, sino que son algoritmos exploratorios que permiten dar una idea de los
procesos implicados, los cuales deben ser evaluados sobre la base de su
plausibilidad biológica. La diversidad de variables que influencian estos métodos
hacen que los tests de corrección para comparaciones múltiples sean solo una
solución parcial a su baja especificidad 275. Hemos utilizado dos programas
bioinformáticos (Gene Set Enrichment Análisis [GSEA] e Ingenuity Pathway Analysis
[IPA]) que hacen uso del enriquecimiento funcional con el objetivo de identificar
vías de señalización celular y procesos biológicos desregulados en la progresión
desde los subtipos de MI a CG (apartado 3.4.1 de la sección Resultados) y para
comparar la MI que no progresa a CG con la mucosa gástrica sana.
La principal ventaja del análisis del programa GSEA es que tiene en cuenta todos los
genes analizados independientemente de su significación estadística o nivel de
expresión. Permitiendo que genes con un efecto pequeño, que no pasarían los
criterios de selección habituales (p<0.05 y Fold Change>1.5 o 2), contribuyan al
enriquecimiento de un conjunto génico. Su principal limitante es que asume que los
genes con mayores diferencias de expresión son los más importantes 275.
Se observa una drástica diferencia entre la MII respecto a la MIC que progresa a CG;
ya que en la primera se obtienen 124 conjuntos génicos sobre-expresados (Tabla
anexa 3) y relacionados con diferentes procesos carcinogénicos, mientras que solo 2
sin plausibilidad biológica en el caso de la MIC. Esto concuerda con el hecho de un
mayor número de genes desregulados implicados en la carcinogénesis y específicos
de la MII-CG (N=58) en comparación con la MIC-CG (N=10). El mecanismo por el
cual las MIC progresan a CG es desconocido, sin embargo es probable que en este
caso sea relevante el papel de unos pocos genes diferencialmente expresados en
lugar de vías de señalización celular, que sea necesario un mayor tiempo para que
en la MIC comience un proceso de desregulación de la expresión génica como en el
grupo MII-CG o que al pasar del tiempo la MIC se transforme en MII y exhiba una
desregulación en la expresión de genes responsables de la progresión tumoral como
la observada en esta última. Pudiera ser de importancia en este aspecto la molécula
de adhesión celular que promueve la proliferación celular en el CG, OLFM4 265, la
cual se encuentra sobre-expresada en la MIC-CG.
219
Discusión
220
Discusión
221
Discusión
222
Discusión
Otros estudios moleculares han mostrado diferencias entre la MII y la MIC que
apoyan el hecho epidemiológico de un mayor riesgo de progresión para la MII
respecto a la MIC. Por ejemplo, la presencia de H. pylori intraepitelial o intercelular
se ha observado con mayor frecuencia en la MII en comparación con la MIC,
induciendo una mayor respuesta inflamatoria debido a la existencia de una
estructura formada por cagA, vacA, proteínas poliubiquitinadas, componentes del
proteasoma y proteínas oncogénicas SHP2 y ERK 75 286. Otro proceso es una menor
expresión de CDX2 en la MII en comparación con la MIC, induciendo un estadio
indiferenciado susceptible de transformación tumoral. En el mismo la acción
supresora tumoral de CDX2 sería menor en la MII permitiendo la proliferación de
células metaplásicas susceptibles de convertirse en tumorales 65 287. SOX2, una
proteína esencial para mantener la pluripotencia de las células madres
embrionarias, muestra una alta expresión en la mayoría de las MII (85%) mientras
que solo en el 7% de las MIC 288. También, la vía de señalización Shh tiene un
menor nivel de activación en la MII 70. Esta vía regula la proliferación y
diferenciación gástrica, y disminuye su activación a partir de la atrofia gástrica. La
inhibición de esta vía conlleva a un programa de diferenciación alterado y pérdida
223
Discusión
En el análisis del GSEA no todos los genes contribuyen de igual manera a un proceso
identificado, sino que existe un subconjunto de los mismos (leading edge) que
contribuyen más a la señal de enriquecimiento (enrichment score). Debido a lo
anterior, a la baja especificidad del análisis de enriquecimiento funcional y la alta
cantidad de conjuntos génicos sobre-expresados (N=124) en la MII-CG, y en un
intento por restringir los mismos a los más relevantes se identificaron los genes
leading edge de los mismos compartidos con los genes sobre-expresados en este
grupo 260.
224
Discusión
225
Discusión
226
Discusión
227
Discusión
228
Discusión
específicos y otros comunes a los grupos MII-CG y MI-No CG (Figura 53), siendo
estos últimos característicos de la MI independientemente de su progresión
tumoral. Los procesos funcionales sobre-expresados en la MI-No CG poseén una
alta similitud con resultados previos de la literatura. Mediante la superposición de
diferentes criterios de selección hemos reducido la gran cantidad de procesos
funcionales existentes en MII-CG y MI-No CG a los más relevantes.
La mayoría de las vías obtenidas son más significativas y con un mayor número de
moléculas en la MII-CG que en la MIC-CG (Figura 40); hecho que coincide con el
mayor riesgo de progresión a CG de la MII en comparación con la MIC obtenido en
estudios epidemiológicos 57 60. Todas las vías sobre-expresadas en MII-CG y MIC-CG
(Tabla 40) contienen a moléculas sobre-expresadas de HLA clase I/II y componentes
de anticuerpos. De manera coincidente, en el análisis del GSEA también se obtienen
vías implicadas en la presentación antigénica (MII-CG y MI-No CG), maduración del
fagosoma (fagocitosis, MII-CG) y maduración de células dendríticas (MII-CG).
229
Discusión
La red molecular construida por el IPA con mayor score en la MII-CG fue asociada
con enfermedades del tejido conectivo e inflamatorias, estando compuesta por 21
moléculas centrales diferencialmente expresadas en la MII-CG y con una extensa
inter-relación (Figura 40A). Por ejemplo, HLA-A, C, DRB1, DQB1, IL1R2, FABP5
interactúan con inmunoglobulinas. También varias moléculas interactúan con el
oncogén RHOA, mutado frecuentemente en CG 46 y sub-expresado en la MII-CG.
Otras moléculas sobre-expresadas interaccionan con el factor transcripcional NFKB,
inductor de la expresión de citoquinas pro-inflamatorias. El IPA no construyó
ninguna red molecular de alto score (>15) en la MIC-CG, debido a la pequeña
cantidad y baja inter-relación de los genes diferencialmente expresados en este
grupo, lo cual apoya el fenómeno de un mayor riesgo de progresión a CG de la MII
respecto a la MIC.
230
Discusión
Al igual que en el análisis con el GSEA las muestras de MII y MIC que no progresan a
CG fueron unidas y comparadas con los controles sanos, con el objetivo de
identificar mecanismos funcionales característicos del grupo MI-No CG. Observamos
que en esta comparación aparecen sobre-expresadas 3 vías pertenecientes a la
familia de receptores nucleares pregnano X (PXR/NR1I2), farnesoide X (FXR/NR1H4),
hepático X (LXR) y retinoide X (RXRG). PXR es un factor transcripcional activado por
compuestos xenobióticos, tras lo cual actúa como un regulador transcripcional de
los citocromos P450, encargados de metabolizar y eliminar los mismos.
Coincidentemente está sobre-expresada la vía metabolismo de xenobióticos,
obtenida como desregulada empleando el programa GSEA y a nivel de genes
individuales. El receptor FXR es miembro de la subfamilia de receptores nucleares
activados por esteroides, siendo un receptor de ácidos biliares y activador
transcripcional de la síntesis y transporte de los mismos.
Las funciones moleculares del IPA confirman, al igual que el GSEA y los genes
diferencialmente expresados, que el metabolismo lipídico está incrementado en la
MI-No CG, siendo dos de los reguladores upstream en este grupo los factores
transcripcionales nucleares hepáticos HNF1A/4A, inductores de este proceso. En la
red de mayor score (29) se observa a dicho proceso como el más relevante,
existiendo un nodo central compuesto por HNF4A sobre-expresado e interactuando
directamente con otras moléculas sobre-expresadas del metabolismo lipídico como
APOA, APOB, MTTP, APOC3 y ABCG8; así como, con el citocromo P450 CYP3A4 y el
receptor pregnano NR1I2, evidenciando también la relevancia del metabolismo de
xenobióticos en dicha red. Además, la función transporte molecular y bioquímica de
moléculas pequeñas pudiera estar reflejando el hecho de que en este grupo fue
obtenido el efecto Warburg como el proceso molecular más frecuentemente sobre-
expresado en el análisis de GSEA. De esta manera, en la MI existe un incremento del
metabolismo lipídico corroborado mediante el análisis de genes diferencialmente
expresados, el de funciones y redes moleculares del IPA y el GSEA.
231
Discusión
Se intento validar por qRT-PCR las diferencias de expresión entre casos y controles
de 37 genes (MII-CG N=19, MIC-CG N=3, MI-No CG N=18) significativos del
microarray de expresión, así como de 3 integrantes de la familia Schlafen entre los
que se encuentra Schlafen5. El mismo no estaba diferencialmente expresado en el
microarray; sin embargo se encuentra sobre-expresado en la qRT-PCR; lo cual
concuerda con su sobre-expresión inmuno-histoquímica en ambos subtipos
histológicos de MI-CG en comparación con MI-No CG descrita en la sección 4 de
Discusión. Se han validado las diferencias de expresión en un 56.75% de los genes
elegidos, la baja efectividad de la qRT-PCR la atribuimos a la mala calidad de las
muestras de ARN procedentes de tejido parafinado con una media de 16.4 años de
antigüedad y al hecho de que la plataforma utilizada no ha sido validada con dicho
ARN de tan mala calidad sino que aconsejan usar muestras de RIN alto (7-8)
(comunicación personal).
232
Discusión
233
Discusión
234
Discusión
Es importante destacar que el grupo MII-CG posee más muestras y mayores valores
de score H en las categorías Débil+ (73.68) y Moderado+ (10.53) en comparación
con el grupo MIC-CG (Débil+ (66.67) y Moderado+ (0)), y que el incremento del área
bajo la curva para el modelo de MII-CG es mayor que el que se obtiene para el
modelo MI-CG que agrupa los dos subtipos de metaplasia (MII y MIC). Todo ello
evidencia una mayor expresión de Schlafen5 en las MII con respecto a las MIC que
progresan a CG; así como, que la expresión de SLFN5 contribuye a mejorar el
modelo de predicción del CG en comparación a cuando se analizan sólo ambos tipos
histológicos de MI. El modelo de predicción basado solamente en la MIC no se
realizó debido al bajo tamaño de muestra de este grupo (N=6).
235
Discusión
236
Discusión
237
Discusión
Concluimos que Schlafen5 pudiera ser inducido por la acción del IFNα en células T
de la mucosa gástrica, e influir en la progresión desde MI a CG debido a la
disminución de la respuesta anti-tumoral de células T. Su tinción
inmunohistoquímica e interpretación cuantitativa en el estroma gástrico de
pacientes con metaplasia intestinal, especialmente en los que tienen una MII,
podría ser un biomarcador útil para detectar aquellos con mayor riesgo de
desarrollar un cáncer gástrico en el futuro. Estudios de validación en series
independientes de pacientes con metaplasia intestinal que progresan y no
progresan a CG; asi como tinción inmunohisto-química de marcador CD30 para
chequear si las células T se encuentran activas funcionalmente 320 son necesarios
para apoyar esta afirmación.
238
CONCLUSIONES
Conclusiones
3. Otros genes asociados con la evolución de lesiones precursoras del CG pero que
no superan la corrección para comparaciones múltiples, son CD14, PTGS2 y TFF1
que se asocian con la progresión, mientras que CDH1, TFF1 y MUC2 se asocian con
la regresión.
4. Otros genes que se observan asociados por primera vez con la evolución de las
lesiones precursoras del CG aunque no superan la corrección para comparaciones
múltiples, son DNAH11, PRKAA1 y PTGS2, que se asocian con la progresión; y
PTGES, CDX2, CDH1, MAP3K14 y PTPN11, que se asocian con la regresión.
241
Conclusiones
10. Es posible que otros genes homeobox diferentes de CDX1/CDX2 como HOXA13
(MI-No CG) y HOXC11 (MII-CG) contribuyan a la diferenciación intestinal en la MI.
242
REFERENCIAS BIBLIOGRÁFICAS
Referencias Bibliográficas
1. Piazuelo M. B. & Correa P. Gastric cancer: Overview. Colomb. Med. 44, 192–201
(2013).
5. Correa, P. Gastric cancer: Overview. Gastroenterol. Clin. North Am. 42, 211–217
(2014).
8. Hall, J. E. (2010). Guyton and Hall Textbook of Medical Physiology. 12ma edición.
St. Louis: W. B. Saunders.
10. Hansson, L. E., Sparén, P. & Nyrén, O. Increasing incidence of carcinoma of the
gastric cardia in Sweden from 1970 to 1985. Br. J. Surg. 80, 374–377 (1993).
13. Lauren, P. The two histological main types of gastric carcinoma: Diffuse and so
called intestinal-type carcinoma. Acta Pathol. Microbiol. Scand. 64, 31–49 (1965).
14. Correa, P. & Piazuelo, M. B. The gastric precancerous cascade. J. Dig. Dis. 13, 2–9
(2012).
245
Referencias Bibliográficas
15. Parkin, D. M. The global health burden of infection-associated cancers in the year
2002. Int. J. Cancer 118, 3030–3044 (2006).
16. Ghoshal, U. C., Chaturvedi, R. & Correa, P. The enigma of Helicobacter pylori
infection and gastric cancer. Indian J. Gastroenterol. 29, 95–100 (2011).
17. Stein, M., Ruggiero, P., Rappuoli, R. & Bagnoli, F. Helicobacter pylori CagA: From
pathogenic mechanisms to its use as an anti-cancer vaccine. Front. Immunol. 4, 328
(2013).
18. Bravo, L. E., van Doom, L. J., Realpe, J. L. & Correa, P. Virulence-associated
genotypes of Helicobacter pylori: Do they explain the African enigma? Am. J.
Gastroenterol. 97, 2839–2842 (2002).
19. Kodaman, N. et al. Human and Helicobacter pylori coevolution shapes the risk of
gastric disease. Proc. Natl. Acad. Sci. U. S. A. 111, 1455–1460 (2014).
20. Palframan, S. L., Kwok, T. & Gabriel, K. Vacuolating cytotoxin A (VacA), a key toxin
for Helicobacter pylori pathogenesis. Front. Cell. Infect. Microbiol. 2, (2012).
22. Cover, T. L. & Blanke, S. R. Helicobacter pylori VacA, a paradigm for toxin
multifunctionality. Nat. Rev. Microbiol. 3, 320–332 (2005).
25. Donnelly, J. M. et al. Mesenchymal stem cells induce epithelial proliferation within
the inflamed stomach. Am. J. Physiol. Gastrointest. Liver Physiol. 306,
G1075– G1088 (2014).
26. Camilo, V. et al. Helicobacter pylori and the BMP pathway regulate CDX2 and SOX2
expression in gastric cells. Carcinogenesis 33, 1985–1992 (2012).
27. Mobley, H. L. T., Mendz, G. L. & Hazell, S. L. (2001). Helicobacter pylori. Physiology
and Genetics. Washington DC: ASM.
28. Higashi, H. et al. Biological activity of the Helicobacter pylori virulence factor CagA
is determined by variation in the tyrosine phosphorylation sites. Proc. Natl. Acad.
Sci. U. S. A. 99, 14428–14433 (2002).
246
Referencias Bibliográficas
29. Kwok, T. et al. Helicobacter exploits integrin for type IV secretion and kinase
activation. Nature 449, 862–866 (2007).
30. Selbach, M., Moese, S., Hauck, C. R., Meyer, T. F. & Backert, S. Src is the kinase of
the Helicobacter pylori CagA protein in vitro and in vivo. J. Biol. Chem. 277,
6775–6778 (2002).
31. Brandt, S., Kwok, T., Hartig, R., König, W. & Backert, S. NF-κB activation and
potentiation of proinflammatory responses by the Helicobacter pylori CagA
protein. Proc. Natl. Acad. Sci. U. S. A. 102, 9300–9305 (2005).
34. Fox, J. G. et al. Concurrent enteric helminth infection modulates inflammation and
gastric immune responses and reduces Helicobacter-induced gastric atrophy. Nat.
Med. 6, 536–542 (2000).
35. Wang, Q. et al. Consumption of fruit, but not vegetables, may reduce risk of gastric
cancer: Results from a meta-analysis of cohort studies. Eur. J. Cancer 50, 1498–509
(2014).
36. Mendez, M. A. et al. Cereal fiber intake may reduce risk of gastric
adenocarcinomas: The EPIC-EURGAST study. Int. J. Cancer 121, 1618–1623 (2007).
38. Agudo, A. et al. Impact of cigarette smoking on cancer risk in the European
prospective investigation into cancer and nutrition study. J. Clin. Oncol. 30,
4550–4557 (2012).
39. Abnet, C. C. et al. Non-steroidal anti-inflammatory drugs and risk of gastric and
oesophageal adenocarcinomas: Results from a cohort study and a meta-analysis.
Br. J. Cancer 100, 551–557 (2009).
41. Hayakawa, Y. et al. Mouse models of gastric cancer. Cancers (Basel). 5, 92–130
(2013).
247
Referencias Bibliográficas
43. Nobili, S. et al. Genomic and genetic alterations influence the progression of gastric
cancer. World J. Gastroenterol. 17, 290–299 (2011).
44. Oliveira, C., Pinheiro, H., Figueiredo, J., Seruca, R. & Carneiro, F. E-cadherin
alterations in hereditary disorders with emphasis on hereditary diffuse gastric
cancer. Prog. Mol. Biol. Transl. Sci. 116, 337–359 (2013).
45. Hanahan, D. & Weinberg, R. A. Hallmarks of cancer: The next generation. Cell 144,
646–674 (2011).
47. Dixon, M. F., Genta, R. M., Yardley, J. H. & Correa, P. Classification and grading of
gastritis. The updated Sydney System. Am. J. Surg. Pathol. 20, 1161–1181 (1996).
48. Den Hoed, C. M. et al. Follow-up of premalignant lesions in patients at risk for
progression to gastric cancer. Endoscopy 45, 249–256 (2013).
49. Vannella, L. et al. Risk factors for progression to gastric neoplastic lesions in
patients with atrophic gastritis. Aliment. Pharmacol. Ther. 31, 1042–1050 (2010).
50. Filipe, M. I. et al. Intestinal metaplasia types and the risk of gastric cancer: A cohort
study in Slovenia. Int. J.Cancer 57, 324–329 (1994).
52. Barros, R., Freund, J. N., David, L. & Almeida, R. Gastric intestinal metaplasia
revisited: Function and regulation of CDX2. Trends Mol. Med. 18, 555–563 (2012).
54. De Vries, A. C. et al. Gastric cancer risk in patients with premalignant gastric
lesions: A nationwide cohort study in the Netherlands. Gastroenterology 134,
945–952 (2008).
248
Referencias Bibliográficas
56. De Vries, A. C., Haringsma, J. & Kuipers, E. J. The detection, surveillance and
treatment of premalignant gastric lesions related to Helicobacter pylori infection.
Helicobacter 12, 1–15 (2007).
59. Kang, K. P. et al. Role of intestinal metaplasia subtyping in the risk of gastric cancer
in Korea. J. Gastroenterol. Hepatol. 24, 140–148 (2009).
60. González, C. A., Sanz-Anquela, J. M., Gisbert, J. P. & Correa, P. Utility of subtyping
intestinal metaplasia as matker of gastric cancer risk. A review of evidence. Int. J.
Cancer 133, 1023–1032 (2013).
62. Correa, P., Piazuelo, M. B. & Wilson, K. T. Pathology of gastric intestinal metaplasia:
Clinical implications. Am. J. Gastroenterol. 105, 493–498 (2010).
63. Fujii, Y. et al. CDX1 confers intestinal phenotype on gastric epithelial cells via
induction of stemness-associated reprogramming factors SALL4 and KLF5. Proc.
Natl. Acad. Sci. U. S. A. 109, 20584–20589 (2012).
65. Liu, Q. et al. CDX2 expression is progressively decreased in human gastric intestinal
metaplasia, dysplasia and cancer. Mod. Pathol. 20, 1286–1297 (2007).
66. Farinati, F. et al. Helicobacter pylori, inflammation, oxidative damage and gastric
cancer: A morphological, biological and molecular pathway. Eur. J. Cancer Prev. 17,
195–200 (2008).
67. Verhulst, M. L., Van Oijen, A. H. A. M., Roelofs, H. M. J., Peters, W. H. M. & Jansen,
J. B. M. J. Antral glutathione concentration and glutathione S-transferase activity in
patients with and without Helicobacter pylori. Dig. Dis. Sci. 45, 629–632 (2000).
68. Osaki, M. et al. Expression of RUNX3 protein in human gastric mucosa, intestinal
metaplasia and carcinoma. Eur. J. Clin. Invest. 34, 605–612 (2004).
249
Referencias Bibliográficas
69. Shiao, Y. H., Rugge, M., Correa, P., Lehmann, H. P. & Scheer, W. D. P53 alteration in
gastric precancerous lesions. Am. J. Pathol. 144, 511–517 (1994).
70. Shiotani, A. et al. Evidence that loss of sonic hedgehog is an indicator of Helicobater
pylori-induced atrophic gastritis progressing to gastric cancer. Am. J. Gastroenterol.
100, 581–587 (2005).
71. Lee, J. H. et al. Frequent CpG island methylation in precursor lesions and early
gastric adenocarcinomas. Oncogene 23, 4646–4654 (2004).
73. Jeong, C. W., Lee, J. H., Sohn, S. S., Ryu, S. W. & Kim, D. K. Mitochondrial
microsatellite instability in gastric cancer and gastric epithelial dysplasia as a
precancerous lesion. Cancer Epidemiol. 34, 323–327 (2010).
75. Necchi, V. et al. Intracellular, intercellular, and stromal invasion of gastric mucosa,
preneoplastic lesions, and cancer by Helicobacter pylori. Gastroenterology 132,
1009–1023 (2007).
76. Meireles, S. I. et al. Molecular classifiers for gastric cancer and nonmalignant
diseases of the gastric mucosa. Cancer Res. 64, 1255–1265 (2004).
80. Chen, X. et al. Variation in gene expression patterns in human gastric cancers. Mol.
Biol. Cell 14, 3208–3215 (2003).
81. Lee, H. J. et al. Gene expression profiling of metaplastic lineages identifies CDH17
as a prognostic marker in early stage gastric cancer. Gastroenterology 139,
358–366 (2010).
250
Referencias Bibliográficas
82. Mavrommatis, E., Fish, E. N. & Platanias, L. C. The Schlafen family of proteins and
their regulation by interferons. J. Interf. Cytokine Res. 33, 206–210 (2013).
86. Van Zuylen, W. J. et al. Macrophage activation and differentiation signals regulate
Schlafen-4 gene expression: Evidence for Schlafen-4 as a modulator of
myelopoiesis. PLoS One 6, (2011).
89. Geserick, P., Kaiser, F., Klemm, U., Kaufmann, S. H. E. & Zerrahn, J. Modulation of T
cell development and activation by novel members of the Schlafen (slfn) gene
family harbouring an RNA helicase-like motif. Int. Immunol. 16, 1535–1548 (2004).
91. Hofker, M. H., Fu, J. & Wijmenga, C. The genome revolution and its role in
understanding complex diseases. B. B. A. Mol. Basis Dis. 1842, 1889–1895 (2014).
92. Gabriel, S. B. et al. The structure of haplotype blocks in the human genome. Science
296, 2225–2229 (2002).
93. Smith, A. V. Manipulating HapMap data using HaploView. Cold Spring Harb. Protoc.
3, (2008).
94. Pasche, B. & Yi, N. Candidate gene association studies: Successes and failures. Curr.
Opin. Genet. Dev. 20, 257–261 (2010).
95. Gao, L., Nieters, A. & Brenner, H. Meta-analysis: Tumour invasion-related genetic
polymorphisms and gastric cancer susceptibility. Aliment. Pharmacol. Ther. 28,
565–573 (2008).
251
Referencias Bibliográficas
96. Loh, M. et al. Meta-analysis of genetic polymorphisms and gastric cancer risk:
Variability in associations according to race. Eur. J. Cancer 45, 2562–2568 (2009).
97. González, C. A., Sala, N. & Capellá, G. Genetic susceptibility and gastric cancer risk.
Int. J. Cancer 100, 249–260 (2002).
100. Vincenzi, B. et al. Interleukin 1β-511T gene (IL1β) polymorphism is correlated with
gastric cancer in the Caucasian population: Results from a meta-analysis. Oncol.
Rep. 20, 1213–1220 (2008).
102. Persson, C., Canedo, P., Machado, J. C., El-Omar, E. M. & Forman, D.
Polymorphisms in inflammatory response genes and their association with gastric
cancer: A HuGE systematic review and meta-analyses. Am. J. Epidemiol. 173,
259–270 (2011).
104. Tahara, T., Shibata, T., Hirata, I., Nakano, H. & Arisawa, T. CD14 promoter-159
polymorphism is associated with reduced risk of intestinal-type gastric cancer in a
Japanese population. Dig. Dis. Sci. 54, 1508–1512 (2009).
105. Hold, G. L. et al. CD14-159C/T and TLR9-1237T/C polymorphisms are not associated
with gastric cancer risk in Caucasian populations. Eur. J. Cancer Prev. 18, 117–119
(2009).
106. Wang, J. J. et al. Association between CD14 gene polymorphisms and cancer risk: A
meta-analysis. PLoS One 9, (2014).
107. Zhou, W. et al. The -159C/T polymorphism in the CD14 gene and cancer risk: A
meta-analysis. Onco. Targets. Ther. 7, 5–12 (2013).
108. Liu, J., He, C., Xu, Q., Xing, C. & Yuan, Y. NOD2 polymorphisms associated with
cancer risk: A meta-analysis. PLoS One 9, e89340 (2014).
252
Referencias Bibliográficas
111. Zhang, K., Zhou, B., Wang, Y., Rao, L. & Zhang, L. The TLR4 gene polymorphisms and
susceptibility to cancer: A systematic review and meta-analysis. Eur. J. Cancer 49,
946–954 (2013).
112. Lo, S. S., Chen, J. H., Wu, C. W. & Lui, W. Y. Functional polymorphism of NFKB1
promoter may correlate to the susceptibility of gastric cancer in aged patients.
Surgery 145, 280–285 (2009).
113. Wang, G. Y., Lu, C. Q., Zhang, R. M., Hu, X. H. & Luo, Z. W. The E-cadherin gene
polymorphism -160C->A and cancer risk: A HuGE review and meta-analysis of 26
case-control studies. Am. J. Epidemiol. 167, 7–14 (2008).
114. Gao, L., Nieters, A. & Brenner, H. Cell proliferation-related genetic polymorphisms
and gastric cancer risk: Systematic review and meta-analysis. Eur. J. Hum. Genet.
17, 1658–1667 (2009).
116. Chen, B., Zhou, Y., Yang, P. & Wu, X. T. ERCC2 Lys751Gln and Asp312Asn
polymorphisms and gastric cancer risk: A meta-analysis. J. Cancer Res. Clin. Oncol.
137, 939–946 (2011).
117. Visscher, P. M., Brown, M. A., McCarthy, M. I. & Yang, J. Five years of GWAS
discovery. Am. J. Hum. Genet. 90, 7–24 (2012).
119. Abnet, C. C. et al. A shared susceptibility locus in PLCE1 at 10q23 for gastric
adenocarcinoma and esophageal squamous cell carcinoma. Nat. Genet. 42, 764–
767 (2010).
121. Shi, Y. et al. A genome-wide association study identifies new susceptibility loci for
non-cardia gastric cancer at 3q13.31 and 5p13.1. Nat. Genet. 43, 1215–1218
(2011).
253
Referencias Bibliográficas
123. Peleteiro, B. et al. Association between cytokine gene polymorphisms and gastric
precancerous lesions: Systematic review and meta-analysis. Cancer Epidemiol.
Biomarkers Prev. 19, 762–776 (2010).
124. Pereira, C. et al. -765G > C COX-2 polymorphism may be a susceptibility marker for
gastric adenocarcinoma in patients with atrophy or intestinal metaplasia. World J.
Gastroenterol. 12, 5473–5478 (2006).
127. Togawa, S. et al. Interleukin-2 gene polymorphisms associated with increased risk
of gastric atrophy from Helicobacter pylori infection. Helicobacter 10, 172–178
(2005).
128. Li, Z. W. et al. Inflammatory cytokine gene polymorphisms increase the risk of
atrophic gastritis and intestinal metaplasia. World J. Gastroenterol. 16, 1788–1794
(2010).
132. Pabalan, N., Singh, N., Pineda, M. R. & Jarjanazi, H. Meta-analysis of the association
between PTPN11 G/A polymorphism at intron 3 with risk of gastric atrophy among
East Asians. J. Gastrointest. Cancer 45, 319–324 (2014).
254
Referencias Bibliográficas
134. Li, W. Q. et al. Association between genetic polymorphisms of DNA base excision
repair genes and evolution of precancerous gastric lesions in a Chinese population.
Carcinogenesis 30, 500–505 (2009).
135. You, W. C. et al. Genetic polymorphisms of CYP2E1, GSTT1, GSTP1, GSTM1, ALDH2,
and ODC and the risk of advanced precancerous gastric lesions in a Chinese
population. Cancer Epidemiol. Biomarkers Prev. 14, 451–458 (2005).
136. Chen, S. Y. et al. Modification effects of GSTM1, GSTT1 and CYP2E1 polymorphisms
on associations between raw salted food and incomplete intestinal metaplasia in a
high-risk area of stomach cancer. Int. J. Cancer 108, 606–612 (2004).
137. Carvalho, F. et al. MUC1 gene polymorphism and gastric cancer an epidemiological
study. Glycoconj. J. 14, 107–111 (1997).
138. Silva, F. et al. MUC1 polymorphism confers increased risk for intestinal metaplasia
in a Colombian population with chronic gastritis. Eur. J. Hum. Genet. 11, 380–384
(2003).
139. Silva, F. et al. MUC1 gene polymorphism in the gastric carcinogenesis pathway. Eur.
J. Hum. Genet. 9, 548–552 (2001).
140. Marín, F. et al. Genetic variation in MUC1, MUC2 and MUC6 genes and evolution of
gastric cancer precursor lesions in a long-term follow-up in a high-risk area in
Spain. Carcinogenesis 33, 1072–1080 (2012).
142. Hishida, A. et al. Toll-like receptor 4 +3725 G/C polymorphism, Helicobacter pylori
seropositivity, and the risk of gastric atrophy and gastric cancer in Japanese.
Helicobacter 14, 47–53 (2009).
143. Riboli, E. et al. European Prospective Investigation into Cancer and Nutrition (EPIC):
Study populations and data collection. Public Health Nutr. 5, 1113–1124 (2002).
144. Carneiro, F. et al. Pathology findings and validation of gastric and esophageal
cancer cases in a European cohort (EPIC/EURGAST). Scand. J. Gastroenterol. 42,
618–627 (2007).
145. Sala, N. et al. Prostate stem-cell antigen gene is associated with diffuse and
intestinal gastric cancer in Caucasians: Results from the EPIC-EURGAST study. Int. J.
Cancer 130, 2417–2427 (2012).
255
Referencias Bibliográficas
146. Palli, D. et al. CagA+ Helicobacter pylori infection and gastric cancer risk in the EPIC-
EURGAST study. Int. J. Cancer 120, 859–867 (2007).
147. So, H. C. & Sham, P. C. Multiple testing and power calculations in genetic
association studies. Cold Spring Harb. Protoc. 2011, pdb.top95 (2011).
149. McLaren, W. et al. Deriving the consequences of genomic variants with the
Ensembl API and SNP Effect Predictor. Bioinformatics 26, 2069–2070 (2010).
150. Reumers, J. et al. Joint annotation of coding and non-coding single nucleotide
polymorphisms and mutations in the SNPeffect and PupaSuite databases. Nucleic
Acids Res. 36, 825–829 (2008).
152. Schadt, E. E. et al. Mapping the genetic architecture of gene expression in human
liver. PLoS Biol. 6, e107 (2008).
153. Gibbs, J. R. et al. Abundant quantitative trait loci exist for DNA methylation and
gene expression in human brain. PLoS Genet. 6, e1000952 (2010).
154. Myers, A. J. et al. A survey of genetic human cortical gene expression. Nat. Genet.
39, 1494–1499 (2007).
156. Zeller, T. et al. Genetics and beyond - the transcriptome of human monocytes and
disease susceptibility. PLoS One 5, (2010).
158. Ausubel, F. M. et al. (2003). Short Protocols in Molecular Biology. 3era edición.
New York: John Wiley & Sons.
159. Birnboim, C. H. (2004). Composition and methods for obtaining nucleic acids from
sputum. <http://www.google.com.na/patents/WO2003104251A9?cl=en>
160. Miller, S. A., Dykes, D. D. & Polesky, H. F. A simple salting out procedure for
extracting DNA from human nucleated cells. Nucleic Acids Res. 16, 1215 (1988).
256
Referencias Bibliográficas
161. Fukamachi, H., Ito, K. & Ito, Y. Runx3 -/- gastric epithelial cells differentiate into
intestinal type cells. Biochem. Biophys. Res. Commun. 321, 58–64 (2004).
163. Lefebvre, O. et al. Gastric mucosa abnormalities and tumorigenesis in mice lacking
the pS2 trefoil protein. Science 274, 259–262 (1996).
165. Oshima, H., Oshima, M., Inaba, K. & Taketo, M. M. Hyperplastic gastric tumors
induced by activated macrophages in COX-2/mPGES-1 transgenic mice. EMBO J. 23,
1669–1678 (2004).
166. Mutoh, H. et al. Cdx1 induced intestinal metaplasia in the transgenic mouse
stomach: Comparative study with Cdx2 transgenic mice. Gut 53, 1416–1423 (2004).
170. Maeda, S. et al. H. pylori activates NF-κB through a signaling pathway involving IκB
kinases, NF-κB-inducing kinase, TRAF2, and TRAF6 in gastric cancer cells.
Gastroenterology 119, 97–108 (2000).
171. Meyer-ter-Vehn, T., Covacci, A., Kist, M. & Pahl, H. L. Helicobacter pylori activates
mitogen-activated protein kinase cascades and induces expression of the proto-
oncogenes c-fos and c-jun. J. Biol. Chem. 275, 16064–16072 (2000).
173. Jin, G. et al. Genetic variants at 6p21.1 and 7p15.3 are associated with risk of
multiple cancers in Han Chinese. Am. J. Hum. Genet. 91, 928–934 (2012).
257
Referencias Bibliográficas
175. Benjamini, Y. & Hochberg, Y. Controlling the False Discovery Rate: A practical and
powerful approach to multiple testing. J. Roy. Stat. Soc. B Met. 57, 289–300 (1995).
176. Affymetrix, Inc. Almac Xcel™ Array for FFPE Profiling Data Sheet.
<http://www.carrerasresearch.org/almac-xcel-array-for-ffpe-profiling_38671. pdf>
178. Johnson, W. E., Li, C. & Rabinovic, A. Adjusting batch effects in microarray
expression data using empirical Bayes methods. Biostatistics 8, 118–127 (2007).
179. Smyth, G. K., Michaud, J. & Scott, H. S. Use of within-array replicate spots for
assessing differential expression in microarray experiments. Bioinformatics 21,
2067–2075 (2005).
182. Krämer, A., Green, J., Pollard, J. & Tugendreich, S. Causal analysis approaches in
Ingenuity Pathway Analysis. Bioinformatics 30, 523–530 (2014).
184. Sonachalam, M., Shen, J., Huang, H. & Wu, X. Systems biology approach to identify
gene network signatures for colorectal cancer. Front. Genet. 3, (2012).
185. The Broad Institute. Gene Set Enrichment Analysis GSEA User Guide. Broad Inst.
1–59 (2009).
186. Hanada, K. et al. Helicobacter pylori infection introduces DNA double-strand breaks
in host cells. Infect. Immun. 82, 4182–4189 (2014).
187. Wisnieski, F. Reference genes for quantitative RT-PCR data in gastric tissues and
cell lines. World J. Gastroenterol. 19, 7121 (2013).
258
Referencias Bibliográficas
190. Neumann, B., Zhao, L., Murphy, K. & Gonda, T. J. Subcellular localization of the
Schlafen protein family. Biochem. Biophys. Res. Commun. 370, 62–66 (2008).
191. McCarty, K. S., Miller, L. S., Cox, E. B. & Konrath, J. Estrogen receptor analyses.
Correlation of biochemical and immunohistochemical methods using monoclonal
antireceptor antibodies. Arch. Pathol. Lab. Med. 109, 716–721 (1985).
194. Durães, C. et al. Genetic variants in the IL1A gene region contribute to intestinal-
type gastric carcinoma susceptibility in European populations. Int. J. Cancer 135,
1343–1355 (2014).
195. Wee, A., Teh, M. & Kang, J. Y. Association of Helicobacter pylori with HLA-DR
antigen expression in gastritis. J. Clin. Pathol. 45, 30–33 (1992).
197. Nguyen, T. V et al. Genetic polymorphisms in GSTA1, GSTP1, GSTT1, and GSTM1
and gastric cancer risk in a Vietnamese population. Oncol. Res. 18, 349–355 (2010).
200. Li, L. et al. Critical role of histone demethylase RBP2 in human gastric cancer
angiogenesis. Mol. Cancer 13, 81 (2014).
201. Aquino, P. F. et al. Are gastric cancer resection margin proteomic profiles more
similar to those from controls or tumors? J. Proteome Res. 11, 5836–5842 (2012).
259
Referencias Bibliográficas
202. Ikuta, K., Seno, H. & Chiba, T. Molecular changes leading to gastric cancer: A
suggestion from rare-type gastric tumors with GNAS mutations. Gastroenterology
146, 1417–1418 (2014).
203. Niiya, F. et al. Expression of SART3 tumor-rejection antigen in gastric cancers. Jpn.
J. Cancer Res. 91, 337–342 (2000).
208. Qian, Z. et al. Whole genome gene copy number profiling of gastric cancer
identifies PAK1 and KRAS gene amplification as therapy targets. Gene.
Chromosome. Canc. 53, 883–894 (2014).
209. Li, S., Lu, A. P., Zhang, L. & Li, Y. D. Anti-Helicobacter pylori immunoglobulin G (IgG)
and IgA antibody responses and the value of clinical presentations in diagnosis of
H. pylori infection in patients with precancerous lesions. World J. Gastroenterol. 9,
755–758 (2003).
210. Nam, K. H. et al. Caveolin 1 expression correlates with poor prognosis and focal
adhesion kinase expression in gastric cancer. Pathobiology 80, 87–94 (2013).
211. Watanabe, Y. et al. HLA-DQB1 locus and gastric cancer in Helicobacter pylori
infection. J. Gastroenterol. Hepatol. 21, 420–424 (2006).
212. Wang, Y. C. et al. CD24 mediates gastric carcinogenesis and promotes gastric
cancer progression via STAT3 activation. Apoptosis 19, 643–656 (2014).
213. Xu, H. Y. et al. Transfection of PDCD5 effect on the biological behavior of tumor
cells and sensitized gastric cancer cells to cisplatin-induced apoptosis. Dig. Dis. Sci.
57, 1847–1856 (2012).
260
Referencias Bibliográficas
215. Xie, Ying, Guo, Z., Wang, Y. & Zhao, Y. Single-nucleotide polymorphisms of
microRNA processing machinery genes are associated with risk for gastric cancer.
Onco. Targets. Ther. 8, 567–571 (2015).
216. Wang, J., Cui, S., Zhang, X., Wu, Y. & Tang, H. High expression of heat shock protein
90 is associated with tumor aggressiveness and poor prognosis in patients with
advanced gastric cancer. PLoS One 8, e62876 (2013).
218. Wu, W. et al. S100A9, GIF and AAT as potential combinatorial biomarkers in gastric
cancer diagnosis and prognosis. Proteomics. Clin. Appl. 6, 152–162 (2012).
219. Magnusson P. K. E. et al. Gastric cancer and human leukocyte antigen: Distinct DQ
and DR alleles are associated with development of gastric cancer and infection by
Helicobacter pylori. Cancer Res. 61, 2684–2689 (2001).
220. Lee, J. Y., Eom, E. M., Kim, D. S., Ha-Lee, Y. M. & Lee, D. H. Analysis of gene
expression profiles of gastric normal and cancer tissues by SAGE. Genomics 82,
78–85 (2003).
221. Dai, J., Zhang, N., Wang, J., Chen, M. & Chen, J. Gastrokine-2 is downregulated in
gastric cancer and its restoration suppresses gastric tumorigenesis and cancer
metastasis. Tumour Biol. 35, 4199–4207 (2014).
223. Ishigami, S. et al. Cancerous HLA class I expression and regulatory T cell infiltration
in gastric cancer. Cancer Immunol. Immunother. 61, 1663–1669 (2012).
224. Ning, P. F., Liu, H. J. & Yuan, Y. Dynamic expression of pepsinogen C in gastric
cancer, precancerous lesions and Helicobacter pylori associated gastric diseases.
World J. Gastroenterol. 11, 2545–2548 (2005).
225. Shi, J., Xu, S., He, P. & Xi, Z. Expression of carcinoembryonic antigen-related cell
adhesion molecule 1(CEACAM1) and its correlation with angiogenesis in gastric
cancer. Pathol. Res. Pract. 210, 473–476 (2014).
261
Referencias Bibliográficas
226. Li, X. et al. miRNA-223 promotes gastric cancer invasion and metastasis by
targeting tumor suppressor EPB41L3. Mol. Cancer Res. 9, 824–833 (2011).
228. Zeller, T. et al. Genetics and beyond--the transcriptome of human monocytes and
disease susceptibility. PLoS One 18, e10693 (2010).
229. Corridoni, D., Arseneau, K. O., Cifone, M. G. & Cominelli, F. The dual role of nod-like
receptors in mucosal innate immunity and chronic intestinal inflammation. Front.
Immunol. 5, (2014).
230. Helicobacter and Cancer Collaborative Group. Gastric cancer and Helicobacter
pylori: a combined analysis of 12 case control studies nested within prospective
cohorts. Gut 49, 347–353 (2001).
231. Hnatyszyn, A., Szalata, M., Stanczyk, J., Cichy, W. & Slomski, R. Association of
c.802C>T polymorphism of NOD2/CARD15 gene with the chronic gastritis and
predisposition to cancer in H. pylori infected patients. Exp. Mol. Pathol. 88,
388–393 (2010).
232. Wang, P. et al. Association of NOD1 and NOD2 genes polymorphisms with
Helicobacter pylori related gastric cancer in a Chinese population. World J.
Gastroenterol. 18, 2112–2120 (2012).
234. Lewis, B. P., Green, R. E. & Brenner, S. E. Evidence for the widespread coupling of
alternative splicing and nonsense-mediated mRNA decay in humans. Proc. Natl.
Acad. Sci. U. S. A. 100, 189–192 (2003).
236. Nicod, N., Pradas-Juni, M. & Gomis, R. Role of the single nucleotide polymorphism
rs7903146 of TCF7L2 in inducing nonsense-mediated decay. Springerplus 3, 41
(2014).
237. Karhukorpi, J. et al. Effect of CD14 promotor polymorphism and H. pylori infection
and its clinical outcomes on circulating CD14. Clin Exp Immunol 128, 326–332
(2002).
262
Referencias Bibliográficas
238. Smith, M. G., Hold, G. L., Tahara, E. & El-Omar, E. M. Cellular and molecular aspects
of gastric cancer. World J. Gastroenterol. 12, 2979–2990 (2006).
239. Li, K. et al. CD14 regulates gastric cancer cell epithelial-mesenchymal transition and
invasion in vitro. Oncol. Rep. 30, 2725–2732 (2013).
240. Baldini, M. et al. A Polymorphism* in the 5’ flanking region of the CD14 gene is
associated with circulating soluble CD14 levels and with total serum
immunoglobulin E. Am. J. Respir. Cell Mol. Biol. 20, 976–983 (1999).
242. Tahara, T., Arisawa, T., Shibata, T., Hirata, I. & Nakano, H. Association of
polymorphism of TLR4 and CD14 genes with gastroduodenal diseases in Japan.
Inflammopharmacology 15, 124–128 (2007).
243. De Aguiar, B. B. et al. CD14 expression in the first 24h of sepsis: Effect of -260C>T
CD14 SNP. Immunol. Invest. 37, 752–769 (2008).
245. Loo, L. W. M. et al. Cis-expression QTL analysis of established colorectal cancer risk
variants in colon tumors and adjacent normal tissue. PLoS One 7, (2012).
246. Corvalan, A. et al. Epstein-Barr virus in gastric carcinoma is associated with location
in the cardia and with a diffuse histology: A study in one area of Chile. Int. J. Cancer
94, 527–530 (2001).
247. Yoshiyama, H., Imai, S., Shimizu, N. & Takada, K. Epstein-Barr virus infection of
human gastric carcinoma cells: Implication of the existence of a new virus receptor
different from CD21. J. Virol. 71, 5688–5691 (1997).
248. Pugin, J. et al. CD14 is a pattern recognition receptor. Immunity 1, 509–516 (1994).
263
Referencias Bibliográficas
251. Smith, M. F. et al. Toll-like receptor (TLR) 2 and TLR5, but not TLR4, are required for
Helicobacter pylori-induced NF-κB activation and chemokine expression by
epithelial cells. J. Biol. Chem. 278, 32552–32560 (2003).
253. Sun, X. F. & Zhang, H. NFKB and NFKBI polymorphisms in relation to susceptibility
of tumour and other diseases. Histology and Histopathology 22, 1387–1398 (2007).
254. Kufe, D. W. Mucins in cancer: Function, prognosis and therapy. Nat. Rev. Cancer 9,
874–885 (2009).
255. Conze, T. et al. MUC2 mucin is a major carrier of the cancer-associated sialyl-Tn
antigen in intestinal metaplasia and gastric carcinomas. Glycobiology 20, 199–206
(2010).
256. Boltin, D. & Niv, Y. Mucins in gastric cancer - An update. J. Gastrointest. Dig. Syst. 3,
15519 (2013).
258. Kim, E. J. et al. Helicobacter pylori infection enhances gastric mucosal inflammation
in individuals carrying the 260-T allele of the CD14 gene. Gut Liver 7, 317–322
(2013).
259. Renouf, B. et al. Cdx2 homeoprotein inhibits non-homologous end joining in colon
cancer but not in leukemia cells. Nucleic Acids Res. 40, 3456–3469 (2012).
261. Hensel, F. Ten-year follow-up of a prospective trial for the targeted therapy of
gastric cancer with the human monoclonal antibody PAT-SC1. Oncol. Rep. 31,
1059–1066 (2014).
262. Sasaki, T. et al. AKT activation and telomerase reverse transcriptase expression are
concurrently associated with prognosis of gastric cancer. Pathobiology 81, 36–41
(2013).
264
Referencias Bibliográficas
263. Cook, M. B. et al. Iron in relation to gastric cancer in the α-tocopherol, β-carotene
Cancer Prevention Study. Cancer Epidemiol. Biomarkers Prev. 21, 2033–2042
(2012).
264. Han, Y. et al. Identification and validation that up-expression of HOXA13 is a novel
independent prognostic marker of a worse outcome in gastric cancer based on
immunohistochemistry. Med. Oncol. 30, 564 (2013).
265. Liu, R. et al. Depletion of OLFM4 gene inhibits cell growth and increases
sensitization to hydrogen peroxide and tumor necrosis factor-alpha induced-
apoptosis in gastric cancer cells. J. Biomed. Sci. 19, 38 (2012).
267. De Bolos, C., Real, F. X. & Lopez-Ferrer, A. Regulation of mucin and glycoconjugate
expression: From normal epithelium to gastric tumors. Front. Biosci. 6,
D1256–D1263 (2001).
268. Wong, W. M., Poulsom, R. & Wright, N. A. Trefoil peptides. Gut 44, 890–895 (1999).
269. Liu, T. H., Li, D. C., Gu, C. F. & Ye, S. F. Carbamyl phosphate synthetase I. A novel
marker for gastric carcinoma. Chin. Med. J. (Engl). 102, 630–638 (1989).
270. Deng, Q. et al. Prognostic value of long non-coding RNA HOTAIR in various cancers.
PLoS One 10, e110059 (2014).
272. Chandrasekaran, E. V et al. Potential tumor markers for human gastric cancer: An
elevation of glycan:sulfotransferases and a concomitant loss of α1,2-
fucosyltransferase activities. J. Cancer Res. Clin. Oncol. 133, 599–611 (2007).
273. Pradeep, D. U. et al. TRIM protein mediated regulation of inflammatory and innate
immune signaling and its association with antiretroviral activity. J. Virol. 87,
257–272 (2013).
274. Falaleeva, M., Surface, J., Shen, M., de la Grange, P., Stamm, S. SNORD116 and
SNORD115 change expression of multiple genes and modify each other's activity.
Gene. 572, 266-73 (2015)
265
Referencias Bibliográficas
277. Cito, L., Pentimalli, F., Forte, I., Mattioli, E. & Giordano, A. Rb family proteins in
gastric cancer (review). Oncol. Rep. 24, 1411–1418 (2010).
278. Byun, E., Lim, J. W., Kim, J. M. & Kim, H. α-lipoic acid inhibits Helicobacter pylori
induced oncogene expression and hyperproliferation by suppressing the activation
of NADPH oxidase in gastric epithelial cells. Mediators Inflamm. 2014, 1–12 (2014).
280. Baenke, F., Peck, B., Miess, H. & Schulze, A. Hooked on fat: The role of lipid
synthesis in cancer metabolism and tumour development. Dis. Model. Mech. 6,
1353–1363 (2013).
283. Tsang, Y. H. et al. Helicobacter pylori CagA targets gastric tumor suppressor RUNX3
for proteasome-mediated degradation. Oncogene 29, 5643–5650 (2010).
284. Eguchi, H., Herschenhous, N., Kuzushita, N. & Moss, S. F. Helicobacter pylori
increases proteasome-mediated degradation of p27kip1 in gastric epithelial cells.
Cancer Res. 63, 4739–4746 (2003).
285. Valenzuela, M. et al. Helicobacter pylori-induced loss of survivin and gastric cell
viability is attributable to secreted bacterial gamma-glutamyl transpeptidase
activity. J. Infect. Dis. 208, 1131–1141 (2013).
286. Necchi, V., Sommi, P., Ricci, V. & Solcia, E. In vivo accumulation of Helicobacter
pylori products, NOD1, ubiquitinated proteins and proteasome in a novel
cytoplasmic structure. PLoS One 5, (2010).
266
Referencias Bibliográficas
287. Shiotani, A. et al. Re-expression of sonic hedgehog and reduction of CDX2 after
Helicobacter pylori eradication prior to incomplete intestinal metaplasia. Int. J.
Cancer 121, 1182–1189 (2007).
290. Baik, S., Youn, H. & Chung, M. Increased oxidative DNA damage in Helicobacter
pylori-infected human gastric mucosa. Cancer Res. 56, 1279–1282 (1996).
291. Yuan, X. et al. Activation of TLR4 signaling promotes gastric cancer progression by
inducing mitochondrial ROS production. Cell Death Dis. 4, e794 (2013).
292. Zhang, S. et al. Homeostasis of redox status derived from glucose metabolic
pathway could be the key to understanding the Warburg effect. Am. J. Cancer Res.
5, 928–944 (2015).
293. Barrett, M. T. et al. Evolution of neoplastic cell lineages in Barrett oesophagus. Nat.
Genet. 22, 106–109 (1999).
294. Tamási, V., Monostory, K., Prough, R. A. & Falus, A. Role of xenobiotic metabolism
in cancer: Involvement of transcriptional and miRNA regulation of P450s. Cell. Mol.
Life Sci. 68, 1131–1146 (2011).
296. Tahara, T. et al. COMT gene val158met polymorphism influences the severity of
intestinal metaplasia in Helicobacter pylori infected older subjects.
Hepatogastroenterology. 56, 411–415 (2009).
297. Joncquel Chevalier Curt, M. et al. Alteration or adaptation, the two roads for human
gastric mucin glycosylation infected by Helicobacter pylori. Glycobiology 25, 617–
631 (2015).
298. Karasawa, F. et al. Essential role of gastric gland mucin in preventing gastric cancer
in mice. J Clin Invest, 122, 923-34 (2012).
299. Gomes, C. et al. Glycoproteomic analysis of serum from patients with gastric
precancerous lesions. J. Proteome Res. 12, 1454–1466 (2013).
267
Referencias Bibliográficas
301. Pinho, SS., Reis CA. Glycosylation in cancer: mechanisms and clinical implications.
Nat Rev Cancer, 15, 540-55 (2015).
302. Kaebisch, R., Mejías-Luque, R., Prinz, C. & Gerhard, M. Helicobacter pylori
cytotoxin-associated gene A impairs human dendritic cell maturation and function
through IL-10-mediated activation of STAT3. J. Immunol. 192, 316–323 (2014).
303. Borlace, G. N. et al. A role for altered phagosome maturation in the long-term
persistence of Helicobacter pylori infection. Am. J. Physiol. Gastrointest. Liver
Physiol. 303, G169–G179 (2012).
305. Liu, R. et al. Mechanism of cancer cell adaptation to metabolic stress: Proteomics
identification of a novel thyroid hormone-mediated gastric carcinogenic signaling
pathway. Mol. Cell. Proteomics 8, 70–85 (2009).
306. Bubenik, G. A. Thirty four years since the discovery of gastrointestinal melatonin. J.
Physiol. Pharmacol. 59, 33–51 (2008).
307. Goto, T. et al. Enhanced expression of inducible nitric oxide synthase and
nitrotyrosine in gastric mucosa of gastric cancer patients. Clin. Cancer Res. 5,
1411–1415 (1999).
308. Hong, L. et al. The value of MG7-Ag and COX-2 for predicting malignancy in gastric
precancerous lesions. Cell Biol. Int. 34, 873–876 (2010).
309. Sun, Y., Tian, M. M., Zhou, L. X., You, W. C. & Li, J. Y. Value of c-Met for predicting
progression of precancerous gastric lesions in rural Chinese population. Chin. J.
Cancer Res. 24, 18–22 (2012).
310. Lindholm, C., Quiding-Järbrink, M., Lönroth, H., Hamlet, A. & Svennerholm, A. M.
Local cytokine response in Helicobacter pylori-infected subjects. Infect. Immun. 66,
5964–5971 (1998).
268
Referencias Bibliográficas
313. Mittal, D., Gubin, M. M., Schreiber, R. D. & Smyth, M. J. New insights into cancer
immunoediting and its three component phases-elimination, equilibrium and
escape. Curr. Opin. Immunol. 27, 16–25 (2014).
314. Shafaghi, A. et al. Serum gastrin and the pepsinogen I/II ratio as markers for
diagnosis of premalignant gastric lesions. Asian Pacific J. Cancer Prev. 14,
3931–3936 (2013).
316. Joo, M., Kim, H., Kim, M. K., Yu, H. J. & Kim, J. P. Expression of Ep-CAM in intestinal
metaplasia, gastric epithelial dysplasia and gastric adenocarcinoma.
J. Gastroenterol. Hepatol. 20, 1039–1045 (2005).
319. Junnila, S., Kokkola, A., Lindsberg, M. L. K., Puolakkainen, P. & Monni, O. Genome-
wide gene copy number and expression analysis of primary gastric tumors and
gastric cancer cell lines. BMC Cancer 10, 73 (2010).
320. Muta, H., Podack ER. CD30: from basic research to cancer therapy. Immunol Res.
57,151-8 (2013).
269
Tablas Anexas
TABLAS ANEXAS
Tabla Anexa 1. Resultados globales del genotipado de SNPs en el estudio de asociación en el seguimiento de LPGs.
Posición Progresión Progresión Regresión Regresión Bloque de
Gen SNP a Cromosoma b Alelos d MAF e EHW f
cromosomal c Controles g Casos h Controles i Casos j haplotipos k
CDX1 rs6894617 5 149521614 C/T 0,259 0,397 249/187/27 56/29/10 199/128/27 106/88/10 Recomb
CDX1 rs2302275 5 149527021 C/G 0,47 0,852 128/232/100 27/43/24 102/168/82 53/107/42 1
CDX1 rs10040224 5 149527350 C/G 0,062 1,000 406/56/1 81/15/0 309/45/1 178/26/0 1
CDX1 rs887343 5 149529016 G/A 0,366 0,550 182/224/57 41/42/13 147/163/45 76/103/25 1
CDX1 rs2282812 5 149533940 G/A 0,366 1,000 184/215/61 41/39/15 142/161/50 83/93/26 2
CDX1 rs10063514 5 149546502 C/T 0,068 0,712 400/61/1 80/15/0 306/47/1 174/29/0 Recomb
CDX2 rs2481952 13 27434901 C/T 0,498 0,309 121/219/123 27/48/20 92/173/89 56/94/54 1
CDX2 rs2504212 13 27435537 C/A 0,211 0,008 297/135/30 64/30/2 231/104/19 130/61/13 1
CDX2 rs4503658 13 27442919 G/T 0,08 0,006 397/59/7 84/12/0 312/41/2 169/30/5 1
CDX2 rs3812863 13 27443268 A/G 0,398 0,441 174/213/76 40/41/15 141/164/50 73/90/41 Recomb
GAST rs12453761 17 37120839 A/G 0,48 0,405 129/221/113 32/45/19 109/164/82 52/102/50 No es tagSNP
IL1A rs17561 2 113253684 G/T 0,289 0,178 239/179/45 53/39/4 178/146/31 114/72/18 No es tagSNP
IL1B rs3917366 2 113300351 G/T 0,246 0,381 258/180/24 60/33/3 203/132/19 115/81/8 1
IL1B rs1143633 2 113306938 G/A 0,382 0,168 167/232/61 30/50/15 120/177/55 77/105/21 1
IL1B rs1143627 2 113310858 T/C 0,338 0,409 199/215/48 43/40/12 165/152/36 77/103/24 2
PTGES rs12001450 9 131536460 T/C 0,192 0,073 296/155/11 59/35/2 226/119/9 129/71/4 1
PTGES rs4636306 9 131536838 G/A 0,209 0,567 275/150/18 54/33/4 206/113/15 123/70/7 1
PTGES rs2302821 9 131541702 A/C 0,073 0,294 396/60/4 80/15/1 310/41/3 166/34/2 Recomb
PTGES rs10760634 9 131543692 A/G 0,078 1,000 392/69/2 79/14/3 290/61/4 181/22/1 2
PTGES rs4837404 9 131545480 A/G 0,369 1,000 183/215/63 39/45/12 138/166/50 84/94/25 2
PTGES rs10739757 9 131547731 T/C 0,129 1,000 349/107/7 71/21/3 260/87/7 160/41/3 2
PTGES rs7872802 9 131555551 A/G 0,134 0,841 346/107/9 72/22/2 257/90/8 161/39/3 3
PTGS2 rs2206593 1 184909052 G/A 0,052 1,000 415/47/1 82/12/1 317/35/2 180/24/0 Recomb
271
Tablas Anexas
Tabla Anexa 2. Genes sobre (Fold Change≥3) o sub-expresados (Fold Change≤0.5) en la comparación MI-No CG vs Sanos.
Fold p-valor
Símbolo del gen a Nombre del gen Asociados con MI d
Change b ajustado c
GAST gastrin 0.110 1.85E-06 SI 41
SST somatostatin 0.127 7.18E-09 SI 41
PGC progastricsin (pepsinogen C) 0.130 2.20E-05 SI 314
GKN1 gastrokine 1 0.141 1.02E-03 SI 222
GKN2 gastrokine 2 0.195 1.87E-03 SI 221
SLC5A5 solute carrier family 5 (sodium iodide symporter), member 5 0.228 5.85E-04 Nueva molécula similar 315
CLU clusterin 0.232 6.51E-06 Nueva molécula
GLUL glutamate-ammonia ligase 0.247 9.50E-09 Nueva molécula
ASS1 argininosuccinate synthase 1 3,010 1,95E-06 Nueva molécula
GLS glutaminase 3,013 3,02E-07 Nueva molécula
SLC22A18AS solute carrier family 22 (organic cation transporter), member 3,025 1,97E-08 Nueva molécula similar 315
TRIM40 tripartite motif containing 40
18 antisense 3,043 2,09E-04 Nueva molécula
ATP1A1 ATPase, Na+/K+ transporting, alpha 1 polypeptide 3,049 9,28E-08 Nueva molécula
IRF4 interferon regulatory factor 4 3,065 2,20E-05 Nueva molécula
USP2 ubiquitin specific peptidase 2 3,067 7,83E-09 Nueva molécula
CTSZ cathepsin Z 3,090 9,47E-09 Nueva molécula
HOXB6 homeobox B6 3,118 4,74E-07 Nueva molécula
S100A10 S100 calcium binding protein A10 3,133 1,58E-06 Nueva molécula
CDHR5 cadherin-related family member 5 3,145 9,94E-09 Nueva molécula
ADH6 alcohol dehydrogenase 6 (class V) 3,147 5,83E-11 Nueva molécula
EPCAM epithelial cell adhesion molecule 3,166 8,06E-08 SI 316
MPP1 membrane protein, palmitoylated 1, 55kDa 3,166 2,05E-05 Nueva molécula
HNF4G hepatocyte nuclear factor 4, gamma 3,169 1,37E-07 Nueva molécula similar 80
SLC39A5 solute carrier family 39 (metal ion transporter), member 5 3,194 2,24E-09 Nueva molécula similar 315
IL2RG interleukin 2 receptor, gamma 3,217 1,11E-06 Nueva molécula
277
Tablas Anexas
Fold p-valor
Símbolo del gen a Nombre del gen Asociados con MI d
Change b ajustado c
UGT2A3 UDP glucuronosyltransferase 2 family, polypeptide A3 3,238 8,31E-06 Nueva molécula
GLRX glutaredoxin (thioltransferase) 3,245 3,60E-07 Nueva molécula
CDX2 caudal type homeobox 2 3,284 1,13E-10 SI 64
PDZD3 PDZ domain containing 3 3,300 2,47E-08 Nueva molécula
MS4A10 membrane-spanning 4-domains, subfamily A, member 10 3,307 7,97E-05 Nueva molécula
FAM211A family with sequence similarity 211, member A 3,350 7,52E-08 Nueva molécula
ETHE1 ethylmalonic encephalopathy 1 3,398 4,78E-07 Nueva molécula
TM6SF2 transmembrane 6 superfamily member 2 3,411 3,54E-08 Nueva molécula
CALML4 calmodulin-like 4 3,432 9,50E-07 Nueva molécula
CEACAM18 carcinoembryonic antigen-related cell adhesion molecule 18 3,439 7,72E-07 Nueva molécula similar 225
SLC1A1 solute carrier family 1 (neuronal/epithelial high affinity 3,471 2,07E-06 Nueva molécula
CD55 CD55 molecule,
glutamate decay accelerating
transporter, system Xag),factor
memberfor complement
1 3,484 1,08E-06 SI 261
TMEM150B transmembrane protein
(Cromer blood group) 150B 3,497 5,63E-11 Nueva molécula
LOC100124692 maltase-glucoamylase (alpha-glucosidase) pseudogene 3,508 9,91E-06 SI 196
BTNL8 butyrophilin-like 8 3,511 2,84E-04 Nueva molécula
MSLN mesothelin 3,519 1,17E-06 Nueva molécula
KLK1 kallikrein 1 3,524 1,64E-09 Nueva molécula
ACSL5 acyl-CoA synthetase long-chain family member 5 3,524 1,71E-07 Nueva molécula
CFTR cystic fibrosis transmembrane conductance regulator (ATP-binding 3,530 1,12E-10 SI 206
cassette sub-family C, member 7)
FOLH1 /// FOLH1B folate hydrolase (prostate-specific membrane antigen) 1 /// 3,540 2,67E-05 Nueva molécula
ABCG8 ATP-binding cassette,
folate hydrolase 1B sub-family G (WHITE), member 8 3,550 6,51E-07 Nueva molécula
IGHA1/2/IGHG1/2/3/ immunoglobulin heavy locus /// immunoglobulin heavy constant 3,551 1,96E-06 SI 209
IGHM/IGHV4-31
ACHE acetylcholinesterase 3,570 2,85E-10 Nueva molécula
APOBEC1 apolipoprotein B mRNA editing enzyme, catalytic 3,579 3,26E-06 Nueva molécula
VIL1 villin 1
polypeptide 1 3,602 1,58E-05 SI 80
278
Tablas Anexas
Fold p-valor
Símbolo del gen a Nombre del gen Asociados con MI d
Change b ajustado c
HOXA13 homeobox A13 3,624 1,98E-04 SI 264
FLVCR2 feline leukemia virus subgroup C cellular receptor family, 3,628 6,78E-10 Nueva molécula
TTLL6 tubulin
membertyrosine
2 ligase-like family, member 6 3,628 2,11E-07 Nueva molécula
UGT1A1 /3-10 UDP glucuronosyltransferase 1 family, polypeptide A1 /// 3,637 6,34E-06 Nueva molécula
FAM84A family with sequence similarity
UDP glucuronosyltransferase 84, member
1 family, A
polypeptide A10 /// 3,657 2,33E-07 Nueva molécula
CDA cytidine deaminase
UDP glucuronosyltransferase 1 family, polypeptide A3 /// 3,695 7,51E-06 Nueva molécula
ABP1 amiloride binding protein 1 (amine
UDP glucuronosyltransferase 1 family, oxidase
polypeptide
(copper-A4 /// 3,703 2,90E-08 Nueva molécula
CCL14 /15 chemokine
containing)) /// CCL14-CCL15 3,704 3,73E-06 Nueva molécula
UDP glucuronosyltransferase
(C-C motif) ligand 114family, polypeptide A5 ///
ANXA13 annexin A13
readthrough /// chemokine (C-C motif) 3,710 1,47E-06 Nueva molécula
UDP glucuronosyltransferase 1 family, polypeptide
ligand 15 A6 ///
HHLA2 HERV-H LTR-associating 2 3,739 1,40E-06 Nueva molécula
UDP glucuronosyltransferase 1 family, polypeptide A7 ///
TMEM139 transmembrane protein 139 3,762 2,65E-10 Nueva molécula
LDHA UDP glucuronosyltransferase
lactate dehydrogenase A 1 family, polypeptide A8 /// 3,769 1,37E-09 Nueva molécula
BCMO1 UDP glucuronosyltransferase
beta-carotene 15,15'-monooxygenase1 family, 1polypeptide A9 3,785 5,96E-06 Nueva molécula
ALPI alkaline phosphatase, intestinal 3,793 2,64E-07 Nueva molécula
SLC35G1 solute carrier family 35, member G1 3,804 1,12E-07 Nueva molécula similar 315
PIGR polymeric immunoglobulin receptor 3,809 2,98E-08 SI 317
GSTA1 glutathione S-transferase alpha 1 3,811 5,12E-04 Nueva molécula
NT5E 5'-nucleotidase, ecto (CD73) 3,822 5,04E-06 Nueva molécula
MME membrane metallo-endopeptidase 3,830 2,93E-05 Nueva molécula
FABP1 fatty acid binding protein 1, liver 3,973 1,93E-09 SI 80
CLDN3 claudin 3 3,978 4,95E-13 SI 80
NR1I2 nuclear receptor subfamily 1, group I, member 2 4.045 1.88E-09 Nueva molécula
PEPD peptidase D 4.067 1.61E-08 Nueva molécula
TMEM45B transmembrane protein 45B 4.095 2.34E-09 Nueva molécula
GBA3 glucosidase, beta, acid 3 (cytosolic) 4.121 1.58E-05 Nueva molécula
CREB3L3 cAMP responsive element binding protein 3-like 3 4.132 4.37E-08 Nueva molécula
279
Tablas Anexas
Fold p-valor
Símbolo del gen a Nombre del gen Asociados con MI d
Change b ajustado c
KRT20 keratin 20 4.252 3.41E-06 SI 81
MYO7B myosin VIIB 4.341 1.61E-12 Nueva molécula
CCL25 chemokine (C-C motif) ligand 25 4.402 8.47E-08 Nueva molécula
MOGAT3 monoacylglycerol O-acyltransferase 3 4.426 6.19E-14 Nueva molécula
GCNT3 glucosaminyl (N-acetyl) transferase 3, mucin type 4.432 7.16E-11 Nueva molécula
ABCG2 ATP-binding cassette, sub-family G (WHITE), member 2 4.457 7.18E-09 Nueva molécula
CDX1 caudal type homeobox 1 4.497 7.83E-09 SI 79
OSTalpha organic solute transporter alpha 4.500 4.04E-08 Nueva molécula
AFAP1-AS1 AFAP1 antisense RNA 1 (non-protein coding) 4.522 4.13E-08 Nueva molécula
SH3D21 SH3 domain containing 21 4.535 3.41E-13 Nueva molécula
ATP1B3 ATPase, Na+/K+ transporting, beta 3 polypeptide 4.588 4.66E-10 Nueva molécula
REEP6 receptor accessory protein 6 4.669 9.79E-11 Nueva molécula
SULT1E1 sulfotransferase family 1E, estrogen-preferring, member 1 4.727 6.74E-09 SI 272
SLC5A9 solute carrier family 5 (sodium/glucose cotransporter), 4.787 1.32E-08 Nueva molécula similar 315
CAMK2N1 calcium/calmodulin-dependent
member 9 protein kinase II inhibitor 1 4.827 6.08E-09 Nueva molécula
A1CF APOBEC1 complementation factor 4.837 6.82E-12 Nueva molécula
CLRN3 clarin 3 4.887 4.04E-10 Nueva molécula
VNN1 vanin 1 4.894 2.14E-11 Nueva molécula
ACY3 aspartoacylase (aminocyclase) 3 4.959 2.45E-10 Nueva molécula
PCK2 phosphoenolpyruvate carboxykinase 2 (mitochondrial) 4.969 2.65E-08 Nueva molécula
CEACAM20 Carcinoembryonic antigen-related cell adhesion molecule 20 4.990 2.65E-05 Nueva molécula similar 225
HSD17B2 hydroxysteroid (17-beta) dehydrogenase 2 5.049 6.39E-09 SI 318
DGKA diacylglycerol kinase, alpha 80kDa 5.127 1.18E-11 Nueva molécula
SLC15A1 solute carrier family 15 (oligopeptide transporter), member 5.141 2.36E-05 Nueva molécula similar 315
RBP2 retinol
1 binding protein 2, cellular 5.209 7.58E-07 Nueva molécula
SLC4A7 solute carrier family4, sodium bicarbonate cotransporter, member7 5.213 2.97E-10 Nueva molécula similar 315
280
Tablas Anexas
Fold p-valor
Símbolo del gen a Nombre del gen Asociados con MI d
Change b ajustado c
SERPINA1 serpin peptidase inhibitor, clade A (alpha-1 antiproteinase, 5.213 1.34E-08 Nueva molécula
CIDEB cell death-inducing
antitrypsin), member DFFA-like
1 effector b 5.285 4.04E-11 Nueva molécula
DHRS11 dehydrogenase/reductase (SDR family) member 11 5.307 5.99E-13 Nueva molécula
EFNA2 ephrin-A2 5.318 4.72E-15 Nueva molécula
GK glycerol kinase 5.408 2.45E-10 Nueva molécula
ITLN1 intelectin 1 (galactofuranose binding) 5.583 1.65E-09 Nueva molécula
GIP gastric inhibitory polypeptide 5.606 2.65E-08 Nueva molécula
TM4SF20 transmembrane 4 L six family member 20 5.614 1.64E-09 Nueva molécula
DPP4 dipeptidyl-peptidase 4 5.645 2.58E-08 Nueva molécula
XPNPEP2 X-prolyl aminopeptidase (aminopeptidase P) 2, membrane- 5.756 1.99E-08 Nueva molécula
SPINK4 serine
bound peptidase inhibitor, Kazal type 4 5.780 1.51E-12 Nueva molécula
PLA2G2A phospholipase A2, group IIA (platelets, synovial fluid) 5.804 3.90E-05 Nueva molécula
HEPH hephaestin 5.877 9.01E-11 Nueva molécula
GDA guanine deaminase 5.901 1.15E-09 Nueva molécula
SLC6A8 solute carrier family 6 (neurotransmitter transporter, 5.926 1.37E-09 Nueva molécula similar 315
PRAP1 proline-rich acidic protein
creatine), member 8 1 6.050 9.10E-14 Nueva molécula
SLC46A3 solute carrier family 46, member 3 6.143 2.76E-10 Nueva molécula similar 315
GPA33 glycoprotein A33 (transmembrane) 6.199 4.49E-15 SI, Marcador de Barret
SLC7A9 solute carrier family 7 (glycoprotein-associated amino acid 6.238 3.32E-10 Nueva molécula similar 315
ABCC13 ATP-binding cassette,
transporter light chain,sub-family C (CFTR/MRP),
bo,+ system), member 9 member 13, 6.264 3.61E-09 Nueva molécula
CES2 carboxylesterase
pseudogene 2 6.404 3.80E-12 Nueva molécula
CLDN4 claudin 4 6.422 2.57E-12 SI 80
ZG16 zymogen granule protein 16 homolog (rat) 6.658 3.98E-10 Nueva molécula
MUC4 mucin 4, cell surface associated 6.676 1.37E-06 SI 267
APOC3 apolipoprotein C-III 6.797 4.99E-06 Nueva molécula
MUC13 mucin 13, cell surface associated 6.902 5.11E-13 SI 81
281
Tablas Anexas
Fold p-valor
Símbolo del gen a Nombre del gen Asociados con MI d
Change b ajustado c
CEACAM6 carcinoembryonic antigen-related cell adhesion molecule 6 6.931 4.02E-10 SI 225
APOA1 apolipoprotein
(non-specific cross
A-I reacting antigen) 7.126 3.72E-06 SI 80
APOA4 apolipoprotein A-IV 7.165 2.40E-07 SI 80
ABCC2 ATP-binding cassette, sub-family C (CFTR/MRP), member 2 7.356 5.15E-07 Nueva molécula
SLC17A4 solute carrier family 17 (sodium phosphate), member 4 7.418 2.41E-11 Nueva molécula similar 315
LINC00483 long intergenic non-protein coding RNA 483 7.537 1.19E-10 Nueva molécula
RGS2 regulator of G-protein signaling 2, 24kDa 7.558 1.61E-10 Nueva molécula
TM4SF4 transmembrane 4 L six family member 4 7.934 1.45E-11 Nueva molécula
CYP3A4 cytochrome P450, family 3, subfamily A, polypeptide 4 8.000 1.32E-07 SI 319
DEFA5 defensin, alpha 5, Paneth cell-specific 8.219 1.67E-05 Nueva molécula
SLC5A1 solute carrier family 5 (sodium/glucose cotransporter), 8.427 1.26E-08 Nueva molécula similar 315
SLC13A2 solute
member carrier
1 family 13 (sodium-dependent dicarboxylate 8.586 1.44E-13 Nueva molécula similar 315
ACE2 angiotensin
transporter),I converting
member 2 enzyme (peptidyl-dipeptidase A) 2 8.700 7.16E-11 SI 81
HKDC1 hexokinase domain containing 1 8.951 2.26E-15 Nueva molécula
TFF3 trefoil factor 3 (intestinal) 9.082 7.66E-12 SI 268
CHP2 calcineurin-like EF hand protein 2 9.266 4.13E-14 Nueva molécula
REG3A regenerating islet-derived 3 alpha 10.098 2.93E-05 Nueva molécula similar 214
MALL MARVEL domain-containing protein 1-like /// mal, T-cell 10.793 1.02E-13 Nueva molécula
MOGAT2 monoacylglycerol O-acyltransferase 2
differentiation protein-like 11.345 5.32E-15 Nueva molécula
ONECUT2 one cut homeobox 2 11.448 1.88E-15 Nueva molécula
BTNL3 butyrophilin-like 3 12.711 1.21E-09 Nueva molécula
SLC6A19 solute carrier family 6 (neutral amino acid transporter), member 19 13.567 9.98E-15 SI 315
MEP1B meprin A, beta 14.867 1.44E-13 SI 214
MUC3A /// MUC3B mucin-3A-like ///mucin 3A, cell surface associated /// mucin 3B, cell surface 15.900 5.32E-16 Nueva molécula similar 14
associated
CDH17 cadherin 17, LI cadherin (liver-intestine) 18.026 1.34E-17 SI 80
TMPRSS15 transmembrane protease, serine 15 18.481 2.29E-12 Nueva molécula
282
Tablas Anexas
Fold p-valor
Símbolo del gen a Nombre del gen Asociados con MI d
Change b ajustado c
MUC12 mucin 12, cell surface associated 18.713 1.17E-08 Nueva molécula similar 14
SLC26A3 solute carrier family 26, member 3 20.210 6.36E-11 Nueva molécula similar 315
MUC17 mucin 17, cell surface associated 20.521 3.26E-15 Nueva molécula similar 14
MTTP microsomal triglyceride transfer protein 25.142 5.42E-14 SI 79
CPS1 carbamoyl-phosphate synthase 1, mitochondrial 26.246 7.02E-19 SI 269
MUC2 mucin 2, oligomeric mucus/gel-forming 26.300 1.26E-23 SI 266
SI sucrase-isomaltase (alpha-glucosidase) 30.086 1.05E-15 SI 214
REG4 regenerating islet-derived family, member 4 32.200 1.16E-13 SI 81
CLCA1 chloride channel accessory 1 33.778 6.18E-16 Nueva molécula
OLFM4 olfactomedin 4 37.040 1.10E-18 SI 81, 227
APOB apolipoprotein B (including Ag(x) antigen) 43.231 1.16E-13 SI 80 81
ANPEP alanyl (membrane) aminopeptidase 49.970 2.86E-23 Nueva molécula
ALDOB aldolase B, fructose-bisphosphate 58.404 7.64E-18 SI 80
FABP1 /// PRDM10 fatty acid binding protein 1, liver /// PR domain containing 10 80.560 1.26E-23 SI 80
DMBT1 deleted in malignant brain tumors 1 98.224 2.20E-23 SI 271
a
En el caso de la existencia de más de 1 gen significa que la sonda reconoce varios genes. b Fold Change es el promedio de la expresión de MI-
No CG menos mucosa gástrica sana; los genes están ordenados de menor a mayor por dicho valor. c p-valor resultante del test de
permutaciones corregido mediante FDR (tasa de falsos positivos). d Genes asociados con MI en estudios de expresión, de asociación genética,
proteómicos y otros. “Nueva molécula” significa que es la primera vez que se identifica este gen como diferencialmente expresado en la MI.
“Nueva molécula similar” significa que otro miembro de su familia génica ha sido previamente asociado con la MI.
283
Tablas Anexas
Tabla Anexa 3. Conjuntos génicos sobre-expresados resultado del GSEA catálogo c2All.v.5 para la comparación MII-CG vs MII-No CG.
p-valor q-valor Rank at
Nombre Procesos funcionales a Nb ES c
nominal d FDR e max f
REACTOME_G_BETA_GAMMA_SIGNALLING_THROUGH_PI3KGAMMA Ciclo y Proliferación celular 20 0,5932 0,0070 0,0235 895
REN_BOUND_BY_E2F Ciclo y Proliferación celular 52 0,5709 0,0000 0,0009 3374
REACTOME_AUTODEGRADATION_OF_CDH1_BY_CDH1_APC_C Ciclo y Proliferación celular 52 0,5589 0,0000 0,0013 5579
REACTOME_APC_C_CDH1_MEDIATED_DEGRADATION_OF_CDC20_AND_O Ciclo y Proliferación celular 58 0,5548 0,0000 0,0010 5579
THER_APC_C_CDH1_TARGETED_PROTEINS_IN_LATE_MITOSIS_EARLY_G1
REACTOME_APC_C_CDC20_MEDIATED_DEGRADATION_OF_MITOTIC_PRO Ciclo y Proliferación celular 59 0,5440 0,0000 0,0014 5579
TEINS
REACTOME_G2_M_CHECKPOINTS Ciclo y Proliferación celular 36 0,5457 0,0017 0,0090 4518
OLSSON_E2F3_TARGETS_DN Ciclo y Proliferación celular 40 0,5378 0,0000 0,0080 4553
REACTOME_MAPK_TARGETS_NUCLEAR_EVENTS_MEDIATED_BY_MAP_KINASES Ciclo y Proliferación celular 28 0,5265 0,0068 0,0239 4787
REACTOME_CELL_CYCLE_CHECKPOINTS Ciclo y Proliferación celular 101 0,5191 0,0000 0,0007 5579
KARLSSON_TGFB1_TARGETS_UP Ciclo y Proliferación celular 110 0,5140 0,0000 0,0006 4733
ISHIDA_E2F_TARGETS Ciclo y Proliferación celular 45 0,5046 0,0000 0,0112 6333
REACTOME_REGULATION_OF_MITOTIC_CELL_CYCLE Ciclo y Proliferación celular 70 0,5014 0,0000 0,0034 6112
BIOCARTA_P38MAPK_PATHWAY Ciclo y Proliferación celular 35 0,5000 0,0050 0,0253 6659
PID_PI3KCIPATHWAY Ciclo y Proliferación celular 40 0,4723 0,0032 0,0361 6697
ASTIER_INTEGRIN_SIGNALING Ciclo y proliferación celular 54 0,4491 0,0034 0,0312 3343
KEGG_CELL_CYCLE Ciclo y Proliferación celular 109 0,4192 0,0000 0,0177 6872
IGLESIAS_E2F_TARGETS_UP Ciclo y Proliferación celular 142 0,4177 0,0000 0,0117 5813
WHITFIELD_CELL_CYCLE_G2_M Ciclo y Proliferación celular 190 0,3930 0,0000 0,0138 5781
LABBE_WNT3A_TARGETS_UP Ciclo y Proliferación celular 99 0,3905 0,0063 0,0450 6333
CHANG_CYCLING_GENES Ciclo y Proliferación celular 123 0,3883 0,0000 0,0324 6140
REACTOME_CELL_CYCLE_MITOTIC Ciclo y Proliferación celular 273 0,3824 0,0000 0,0113 6864
WHITFIELD_CELL_CYCLE_M_G1 Ciclo y Proliferación celular 128 0,3820 0,0016 0,0341 4529
WHITFIELD_CELL_CYCLE_G2 Ciclo y Proliferación celular 152 0,3637 0,0000 0,0428 6121
REACTOME_CELL_CYCLE Ciclo y Proliferación celular 329 0,3550 0,0000 0,0246 5615
284
Tablas Anexas
Tabla Anexa 4. Conjuntos génicos sobre-expresados resultado del GSEA catálogo c2All.v.5 para la comparación MI-No CG vs Sanos.
a p-valor q-valor Rank at
Nombre Procesos funcionales Nb ES c
nominal d FDR e max f
KEGG_STARCH_AND_SUCROSE_METABOLISM Efecto Warburg no tumoral 33 0,724 0,002 0,127 1852
REACTOME_AMINO_ACID_TRANSPORT_ACROSS_THE_PLASMA_MEMB Efecto Warburg no tumoral 28 0,666 0,000 0,125 1480
RANE
MOOTHA_GLUCONEOGENESIS Efecto Warburg no tumoral 29 0,651 0,000 0,168 527
REACTOME_GLYCOLYSIS Efecto Warburg no tumoral 24 0,640 0,002 0,124 3435
REACTOME_AMINO_ACID_AND_OLIGOPEPTIDE_SLC_TRANSPORTERS Efecto Warburg no tumoral 44 0,618 0,000 0,142 1480
KEGG_ARGININE_AND_PROLINE_METABOLISM Efecto Warburg no tumoral 48 0,605 0,004 0,153 1415
REACTOME_GLUCOSE_METABOLISM Efecto Warburg no tumoral 56 0,593 0,006 0,126 3503
KEGG_ALANINE_ASPARTATE_AND_GLUTAMATE_METABOLISM Efecto Warburg no tumoral 25 0,563 0,018 0,242 2144
REACTOME_GLUCONEOGENESIS Efecto Warburg no tumoral 28 0,553 0,037 0,245 3118
KEGG_GALACTOSE_METABOLISM Efecto Warburg no tumoral 23 0,546 0,030 0,222 1852
REACTOME_PURINE_METABOLISM Efecto Warburg no tumoral 31 0,538 0,039 0,239 1793
KEGG_GLYCOLYSIS_GLUCONEOGENESIS Efecto Warburg no tumoral 55 0,511 0,016 0,201 1249
PID_HIF1_TFPATHWAY Efecto Warburg no tumoral 62 0,510 0,006 0,168 1741
REACTOME_TRANSPORT_OF_INORGANIC_CATIONS_ANIONS_AND_AM Efecto Warburg no tumoral 82 0,509 0,000 0,130 1480
INO_ACIDS_OLIGOPEPTIDES
REACTOME_TRANSPORT_OF_GLUCOSE_AND_OTHER_SUGARS_BILE_SA Efecto Warburg no tumoral 80 0,492 0,002 0,158 1446
LTS_AND_ORGANIC_ACIDS_METAL_IONS_AND_AMINE_COMPOUNDS
MOOTHA_GLYCOGEN_METABOLISM Efecto Warburg no tumoral 19 0,473 0,039 0,244 3304
REACTOME_TRANSPORT_OF_VITAMINS_NUCLEOSIDES_AND_RELATED Efecto Warburg no tumoral 23 0,435 0,029 0,241 3282
_MOLECULES
REACTOME_METABOLISM_OF_AMINO_ACIDS_AND_DERIVATIVES Efecto Warburg no tumoral 165 0,418 0,016 0,239 3285
GROSS_HYPOXIA_VIA_ELK3_AND_HIF1A_UP Efecto Warburg no tumoral 129 0,407 0,026 0,241 3796
REACTOME_METABOLISM_OF_CARBOHYDRATES Efecto Warburg no tumoral 200 0,379 0,002 0,234 3526
REACTOME_CHYLOMICRON_MEDIATED_LIPID_TRANSPORT Metabolismo lipídico 15 0,843 0,000 0,157 86
REACTOME_LIPOPROTEIN_METABOLISM Metabolismo lipídico 24 0,748 0,000 0,123 86
REACTOME_PEROXISOMAL_LIPID_METABOLISM Metabolismo lipídico 19 0,737 0,002 0,162 2643
REACTOME_ABCA_TRANSPORTERS_IN_LIPID_HOMEOSTASIS Metabolismo lipídico 16 0,683 0,009 0,202 167
290
Tablas Anexas
PUBLICACIÓN
295
Publicación
296
Publicación
297
Publicación
298
Publicación
299
Publicación
300
Publicación
301
Publicación
302
Publicación
303
Publicación
304
Publicación
305