Three general categories of genes with parent-of-origin expression bias emerged from our analysis in known, as well as newly identified loci, which we term consensus, complex, and single SNP loci (,
Table S6–
S12). Consensus loci have multiple SNPs, at least one SNP site above the
P value cutoff of 0.05 in each cross and 100% of SNPs agree with the direction of the parental expression bias in both the F1i and F1r cross. Complex loci have multiple SNPs, with one or more SNPs above the
P value cutoff and one or more SNP sites that differ (ie. biallelic expression, strain or opposite parental bias). Finally, a subset of genes had a single SNP site or multiple SNPs within 32 bp of each other (the size of a single read). Thirty-five of the 47 previously known imprinted genes displayed consensus imprinting in at least one brain region (
Table S6). However, several of these same genes also exhibited complex imprinting in other samples, such that 41 known imprinted genes were identified as complex in one or more brain regions. Seven known imprinted genes were identified based on a single SNP site (). Detailed analysis of the positions of SNPs with parental allelic bias within complex loci revealed genes in which monoallelic SNPs are confined to a specific exon (195 genes), to the 3’UTR (57 genes), or to the final exon and 3’UTR region (39 genes) (), suggesting that, in these genes, the parental allelic effect is restricted to only one or a few transcript isoforms. In a subset of these cases, the same parental bias is confirmed by multiple SNPs in the exon or 3’UTR (). A large proportion of the genes exhibited parental effects in the last exon (including 3’UTR region), but involved disagreements between SNP sites in the same region of the gene (classified as “other” (560 genes)). In some cases, as detailed below, these disagreements appear to be related to the fact that only a subset of the SNPs for a given complex gene are able to distinguish a specific imprinted isoform from other overlapping transcripts arising from the same locus.
Cadherin 15 (
cdh15), a gene prospectively linked to intellectual disability in humans (
31), emerges as a consensus imprinted locus, in which all three SNPs display preferential expression of the paternal allele in independent male and female samples () and by Sequenom on an independent cohort of animals (). Other notable consensus imprinted genes include
Bcl2l1, a major regulator of apoptosis linked to cancer (
32), and
Eif2c2 (also called
argonaute2), involved in microRNA and short-interfering RNA (siRNA)-mediated gene silencing (
33) (
Table S6).
Detailed analysis of complex loci revealed remarkable and so far unsuspected features of parent-of-origin transcription bias (
Table S7–
11). In the
Inpp5f locus, three isoforms have previously been described () with preferential paternal expression of
Inpp5f_v2 and
Inpp5f_v3, while
Inpp5f_v1 was reported biallelic (
34). In our analysis, SNP sites aligning to the
Inpp5f_v2 and
Inpp5f_v3 isoforms confirmed strict paternal expression. Four SNP sites located in exons shared by the
Inpp5f_v1 isoform and an overlapping UCSC annotated transcript (mKIAA0966) indicated a modest and nonsignificant paternal bias in expression in the adult mPFC and E15 brain. However, in the adult male and female POA 73% of transcription at these sites (
P<0.01, in F1i and F1r cross) originated from the maternal allele. A single SNP found in the first exon of
Inpp5f_v1, indicated a modest, nonsignificant maternal expression bias in POA. Thus, our approach resolved complex regional-, developmental stage- and isoform- specific parental bias in the transcriptome.
Recently, a highly complex form of imprinting has been described for the gene
H13, such that some
H13 isoforms are maternally expressed, whereas others are paternally expressed (
35). Our analysis confirmed these results (
fig. S9). Here we find that
Herc3, a host gene for the known PEG,
Nap1l5 showed features indicative of isoform specific imprinting in a manner very similar to
H13 (
fig. S10). Additional examples of complex parent-of-origin effects in the CNS transcriptome are presented in the SOM for
Lsm14a, Pafah1b3, and
Ndel1 (
fig. S11). Other notable complex loci include
cdh2 (neuronal-cadherin), which plays a central role in brain morphogenesis (
36), as well as
arnt2 (
aryl hydrocarbon receptor nuclear transclocator 2), a gene with multiple isoforms that regulates hypothalamic development in concert with other imprinted genes, such as
Ndn (
37). Many genes identified in our analysis exhibited complex patterns of parental allelic effects for which the underlying mechanism and functional significance are not yet clear.
Finally, several loci in our dataset did not display the classical pattern associated with parent-of-origin expression biases, but instead displayed significant differences in the relative expression of the maternal and paternal alleles in F1i versus F1r offspring, which we refer to as cross-effects. These effects were analyzed separately and the findings are detailed in the supplemental data (
fig. S12).