[en] Genome-wide association studies have identified a number of single nucleotide polymorphisms associated with Schizophrenia (SZ). Moreover, increasing body of evidence suggests a complex connection of SNPs and epigenetic regulation of gene expression, which, up to now, is not fully understood. Integrative analyses that use genetic data from multi-omics studies to detect DNA methylation sites associated with gene expression and SZ phenotype might constitute a major approach able to elucidate how SZ-associated SNPs affect the disease traits throughout genetic regulation of transcriptional output. To test this hypothesis we performed an exploratory integrative quantitative association analysis to obtain summary statistics data for SZ severity associated SNPs and methylation sites of 10 drug-na¨¿ve SZ cases. We firstly identified a significant SZ-associated SNP (adjusted P < 8×10¿8 ), located on a non-coding region of chromosome 16 downstream a lncRNA gene (Long intergenic non-coding RNAs), a type of genes known for being dysregulated in SZ. At a less restrictive significant P (< 0.01), we found 341 common genes to significant SNP and CpG SZ-associated sites. We further investigate the association of these genes with a previous SZ-RNA sequencing study containing a set of 200 up and down regulated genes. 16 genes were identified as being differentially expressed in SZ. Remarkably, 3 of them (SHANK2, SGK1 and TCN2 ) had been previously described in the literature as being involved in SZ. The methodology presented here, might be novel and useful tool to further dilucidate SZ physiopathology.

[es] Los estudios de asociación de genoma completo han identificado varios polimorfismos de un solo nucleótido (SNP) asociados con esquizofrenia (SZ). La evidencia sugiere la existencia de una compleja conexión entre los SNPs y la regulación epigenética de la expresión génica, que no ha sido resuelta hasta el momento. Los análisis integrativos con datos genéticos de estudios multi-ómicos podrían constituir un enfoque importante para dilucidar cómo diferentes SNPs asociados con la SZ afectan al fenotipo de la enfermedad a través de la regulación transcripcional. Para comprobar esta hipótesis, se realizó un análisis cuantitativo de asociación en una cohorte de 10 sujetos con SZ con el objetivo de identificar SNPs y sitios de metilación asociados a la gravedad de la enfermedad. En primer lugar, se identificó un SNP significativo asociado con SZ (\textit{P}-valor ajustado $ < $ $8$$\times10^{-8}$), ubicado en una región no codificante del cromosoma 16 próximo a un gen lncRNA, un tipo de genes conocidos por estar disregulados en la enfermedad. Un análisis integrativo de los datos genéticos y de metilación obtenidos (\textit {P} $ <$ 0.01), permitió identificar 341 genes comunes a SNPs y sitios de metilación asociados a SZ. Finalmente, se realizó un análisis integrativo de genes diferencialmente expresados empleando datos de estudio previo de RNA-seq. De los 341 genes totales, 16 fueron identificados como diferencialmente expresados. Notablemente, 3 de ellos (\textit{SHANK2}, \textit{SGK1} y \textit{TCN2}) han sido descritos previamente en la literatura como genes involucrados en SZ. La metodología presentada aquí podría constituir una herramienta novedosa y útil para avanzar en el conocimiento de la fisiopatología de la SZ.