Schizophrenia Genes

Bioinformatics

A karyogram of the genes featured in GWAS of schizophrenia.

Eugene https://fizzics.netlify.app
04-26-2021

We’ve been working with the genetics of schizophrenia, looking at how the expression of genes that play a role in the condition vary across brain regions and neurodevelopment stages. Now, call me old-fashioned, but I like to know how things are laid out, so in this case where these genes are situated. Time for a plot.

The genes in question come from a GWAS study of the PGC. The images are build based on the ggbio and karyoploteR packages from bioconductor.

First, let’s load up the required packages:

Next, we get a list of the genes from the PGC. These are taken from one of the supplementary tables from the 2018 paper mentioned above. We only care about the ones that also feature in the Allen Brain Atlas, so we filter for those.

genes <- genes::pardinas() %>% 
  filter(genes %in% abagenes$hgnc_symbol) %>% 
  pull(genes)

From our gene list we build a data frame of gene information using Human build 38 (grch38 from the annotables package). And then turn his into an appropriate GRanges object (thank you GenomicRanges).

gene_table <- grch38 %>% 
  dplyr::filter(symbol %in% genes) %>% 
  dplyr::filter(chr %in% c(1:22, "X", "Y")) %>% 
  dplyr::mutate(strand = ifelse(strand == 1, "+", "-"))

gene_ranges <- makeGRangesFromDataFrame(gene_table, keep.extra.columns = T)
seqlevelsStyle(gene_ranges) <- "UCSC"

Now for the plot. The plotKaryotype() function does a lqyout of 24 chromosomes with their cytobands, the kpPlotMarkers() adds in the gene labels. The label colours depend on strand; red for + and green for -. I wasn’t too successful avoiding label overlaps, sorry about that.

kp <- plotKaryotype(genome="hg38")
kpPlotMarkers(kp, 
              data=gene_ranges, 
              labels=gene_ranges$symbol, 
              label.color = ifelse(gene_table$strand == "+", "darkred", "olivedrab4"),
              text.orientation = "horizontal",
              label.dist = 0.0001,
              r1=0.5, cex=0.6, adjust.label.position = T)

Our genes are pretty spread out, with no particular pattern of course. But I just feel I know them a little better now.

Corrections

If you see mistakes or want to suggest changes, please create an issue on the source repository.

Reuse

Text and figures are licensed under Creative Commons Attribution CC BY 4.0. Source code is available at https://github.com/eugene100hickey/fizzics, unless otherwise noted. The figures that have been reused from other sources don't fall under this license and can be recognized by a note in their caption: "Figure from ...".

Citation

For attribution, please cite this work as

Eugene (2021, April 26). Euge: Schizophrenia Genes. Retrieved from https://www.fizzics.ie/posts/2021-05-08-schizophrenia-genes/

BibTeX citation

@misc{eugene2021schizophrenia,
  author = {Eugene, },
  title = {Euge: Schizophrenia Genes},
  url = {https://www.fizzics.ie/posts/2021-05-08-schizophrenia-genes/},
  year = {2021}
}