Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2018 Jul 25;9(1):2914.
doi: 10.1038/s41467-018-04960-9.

Unperturbed expression bias of imprinted genes in schizophrenia

Affiliations

Unperturbed expression bias of imprinted genes in schizophrenia

Attila Gulyás-Kovács et al. Nat Commun. .

Abstract

How gene expression correlates with schizophrenia across individuals is beginning to be examined through analyses of RNA-seq from postmortem brains of individuals with disease and control brains. Here we focus on variation in allele-specific expression, following up on the CommonMind Consortium (CMC) RNA-seq experiments of nearly 600 human dorsolateral prefrontal cortex (DLPFC) samples. Analyzing the extent of allelic expression bias-a hallmark of imprinting-we find that the number of imprinted human genes is consistent with lower estimates (≈0.5% of all genes), and thus contradicts much higher estimates. Moreover, the handful of putatively imprinted genes are all in close genomic proximity to known imprinted genes. Joint analysis of the imprinted genes across hundreds of individuals allowed us to establish how allelic bias depends on various factors. We find that age and genetic ancestry have gene-specific, differential effect on allelic bias. In contrast, allelic bias appears to be independent of schizophrenia.

PubMed Disclaimer

Conflict of interest statement

The authors declare no competing interests.

Figures

Fig. 1
Fig. 1
Quantifying allelic bias of expression in human individuals using the RNA-seq read count ratio statistic Sig. The strength of bias towards the expression of the maternal (red) or paternal (blue) allele of a given gene g in individual i is gauged with the count of RNA-seq reads carrying the reference allele (small closed circles), and the count of reads carrying an alternative allele (open squares) at all SNPs for which the individual is heterozygous. The allele with the higher read count tends to match the allele with the higher expression, but measurement errors may occasionally revert this tendency as seen for SNP s3 in gene g1 in individual i2
Fig. 2
Fig. 2
Across-individuals distribution of read count ratio S for each gene indicates substantial variation of allelic bias and that <1% of all genes are imprinted. The vertically arranged four main panels present the empirical distribution of Sg across all individuals given each gene g. The upper two panels are distinct representations (survival plot: 1 – ECDF, and “survival heatmap”) of the same three distributions corresponding to a: PEG10, b: ZNF331, and c: AFAP1. PEG10 and ZNF331, previously found to be imprinted in mice or in other human tissues, and one for AFAP1, a gene without prior evidence. The bottom two survival heatmaps present the distribution of Sg for the top 2% and 100% of the 5307 analyzed genes. These are ranked according to gene score defined as 1 – ECDF (0.9) in the bottom far right panels. The score of PEG10, ZNF331, and AFAP1 is marked by a,b,c, respectively, in green circles. As expected, PEG10 and ZNF331 both score high and rank within the top 30 of all genes suggesting they are also imprinted in the present context, the adult human DLFPC. The bottom panels also indicate that <1% of all genes might be imprinted
Fig. 3
Fig. 3
The top 50 genes ranked by the gene score defined, for gene g, as 1 – Fg(0.9), where Fg is the empirical cumulative distribution function (ECDF) for g. Thus 1 – Fg(0.9), is the fraction of individuals i for which Sig > 0.9. Note that the same ranking and score is shown in the bottom half of Fig. 2. The tan colored bars indicate the fraction of individuals with nearly unbiased expression (Eq. 2). Gene names (y axis) are colored according to prior imprinting status: known imprinted (blue), nearby candidate (green), and distant candidate (red). “X” characters next to gene names indicate mapping bias and/or cis-eQTL effects based on the reference/nonreference allele test (“Methods”), while “0” indicates that total allele count was insufficient for this test
Fig. 4
Fig. 4
Schizophrenia does not affect allelic bias of imprinted genes. Distribution of read count ratio for control, schizophrenic (SCZ), and affectic spectrum (AFF) individuals within each gene that has been considered as imprinted in the DLPFC brain area in this study
Fig. 5
Fig. 5
Allelic bias depends differentially on age across imprinted genes. The panels and colors are consistent with the imprinted genes and psychiatric diagnoses presented in Fig. 4. The differential dependence on age is apparent when comparing PEG3 or ZNF331 (negative dependence) to KCNK9 or RP13-487P22.1 (positive dependence) or to NDN or NLRP2 (no dependence)

References

    1. Noor A, et al. 5q11. 2 Duplication encompassing only the UBE3A gene is associated with developmental delay and neuropsychiatric phenotypes. Hum. Mutat. 2017;36:689–693. doi: 10.1002/humu.22800. - DOI - PubMed
    1. Rees E, et al. Analysis of copy number variations at 15 schizophrenia-associated loci. Br. J. Psychiatry, 2014;204:108–114. doi: 10.1192/bjp.bp.113.131052. - DOI - PMC - PubMed
    1. Gregg C, et al. High-resolution analysis of parent-of-origin allelic expression in the mouse brain. Science, 2010;329:643–648. doi: 10.1126/science.1190830. - DOI - PMC - PubMed
    1. Andergassen D, et al. Mapping the mouse allelome reveals tissue-specific regulation of allelic expression. eLife. 2017;6:e25125. doi: 10.7554/eLife.25125. - DOI - PMC - PubMed
    1. Baran Y, et al. The landscape of genomic imprinting across diverse adult human tissues. Genome Res. 2015;25:927–936. doi: 10.1101/gr.192278.115. - DOI - PMC - PubMed

Publication types