Transcriptome-wide identification of novel imprinted genes in neonatal mouse brain

Dissecting the Immune System through Gene Regulation

The immune system plays a dual role in human health, functioning both as a protector against pathogens and, at times, as a contributor to disease. This feature emphasizes the importance to uncover the underlying causes of its malfunctions, necessitating an in-depth analysis in both pathological and physiological conditions to better understand the immune system and immune disorders. Recent advances in scientific technology have enabled extensive investigations into gene regulation, a crucial mechanism governing cellular functionality. Studying gene regulatory mechanisms within the immune system is a promising avenue for enhancing our understanding of immune cells and the immune system as a whole. The gene regulatory mechanisms, revealed through various methodologies, and their implications in the field of immunology are discussed in this chapter.

Nerve transcriptomes in autoimmune and genetic demyelinating neuropathies: Pathogenic pathway assessment of nerve demyelination

The pathogenesis of autoimmune demyelinating neuropathies is poorly understood compared to inherited demyelinating forms. We performed whole transcriptome (RNA-Seq) using nerve biopsy tissues of patients with different autoimmune and inherited demyelinating neuropathies (CIDP n = 10, POEMS n = 18, DADS n = 3, CMT1 n = 3) versus healthy controls (n = 6). A limited number of differentially expressed genes compared to healthy controls were identified (POEMS = 125, DADS = 15, CMT = 14, CIDP = 5). Divergent pathogenic pathways including inflammatory, demyelinating and neurite regeneration such as with the triggering receptor expressed on myeloid cells (TREM1) part of the immunoglobulin superfamily and RhoGD1 are found. Shared and discordant pathogenic injury are discovered between autoimmune and inherited forms.

Origin and Evolution of Marsupial-specific Imprinting Clusters Through Lineage-specific Gene Duplications and Acquisition of Promoter Differential Methylation

Genomic imprinting is a parent-of-origin-specific expression phenomenon that plays fundamental roles in many biological processes. In animals, imprinting is only observed in therian mammals, with ∼200 imprinted genes known in humans and mice. The imprinting pattern in marsupials has been minimally investigated by examining orthologs to known eutherian imprinted genes. To identify marsupial-specific imprinting in an unbiased way, we performed RNA-seq studies on samples of fetal brain and placenta from the reciprocal cross progeny of two laboratory opossum stocks. We inferred allele-specific expression for >3,000 expressed genes and discovered/validated 13 imprinted genes, including three previously known imprinted genes, Igf2r, Peg10, and H19. We estimate that marsupials imprint ∼60 autosomal genes, which is a much smaller set compared with eutherians. Among the nine novel imprinted genes, three noncoding RNAs have no known homologs in eutherian mammals, while the remaining genes have important functions in pluripotency, transcription regulation, nucleolar homeostasis, and neural differentiation. Methylation analyses at promoter CpG islands revealed differentially methylated regions in five of these marsupial-specific imprinted genes, suggesting that differential methylation is a common mechanism in the epigenetic regulation of marsupial imprinting. Clustering and co-regulation were observed at marsupial imprinting loci Pou5f3-Npdc1 and Nkrfl-Ipncr2, but eutherian-type multi-gene imprinting clusters were not detected. Also differing from eutherian mammals, the brain and placenta imprinting profiles are remarkably similar in opossums, presumably due to the shared origin of these organs from the trophectoderm. Our results contribute to a fuller understanding of the origin, evolution, and mechanisms of genomic imprinting in therian mammals.

Inference of putative cell-type-specific imprinted regulatory elements and genes during human neuronal differentiation

Genomic imprinting results in gene expression bias caused by parental chromosome of origin and occurs in genes with important roles during human brain development. However, the cell-type and temporal specificity of imprinting during human neurogenesis is generally unknown. By detecting within-donor allelic biases in chromatin accessibility and gene expression that are unrelated to cross-donor genotype, we inferred imprinting in both primary human neural progenitor cells and their differentiated neuronal progeny from up to 85 donors. We identified 43/20 putatively imprinted regulatory elements (IREs) in neurons/progenitors, and 133/79 putatively imprinted genes in neurons/progenitors. Although 10 IREs and 42 genes were shared between neurons and progenitors, most putative imprinting was only detected within specific cell types. In addition to well-known imprinted genes and their promoters, we inferred novel putative IREs and imprinted genes. Consistent with both DNA methylation-based and H3K27me3-based regulation of imprinted expression, some putative IREs also overlapped with differentially methylated or histone-marked regions. Finally, we identified a progenitor-specific putatively imprinted gene overlapping with copy number variation that is associated with uniparental disomy-like phenotypes. Our results can therefore be useful in interpreting the function of variants identified in future parent-of-origin association studies.

Variable allelic expression of imprinted genes at the Peg13, Trappc9, Ago2 cluster in single neural cells

Genomic imprinting is an epigenetic process through which genes are expressed in a parent-of-origin specific manner resulting in mono-allelic or strongly biased expression of one allele. For some genes, imprinted expression may be tissue-specific and reliant on CTCF-influenced enhancer-promoter interactions. The Peg13 imprinting cluster is associated with neurodevelopmental disorders and comprises canonical imprinted genes, which are conserved between mouse and human, as well as brain-specific imprinted genes in mouse. The latter consist of Trappc9, Chrac1 and Ago2, which have a maternal allelic expression bias of ∼75% in brain. Findings of such allelic expression biases on the tissue level raise the question of how they are reflected in individual cells and whether there is variability and mosaicism in allelic expression between individual cells of the tissue. Here we show that Trappc9 and Ago2 are not imprinted in hippocampus-derived neural stem cells (neurospheres), while Peg13 retains its strong bias of paternal allele expression. Upon analysis of single neural stem cells and in vitro differentiated neurons, we find not uniform, but variable states of allelic expression, especially for Trappc9 and Ago2. These ranged from mono-allelic paternal to equal bi-allelic to mono-allelic maternal, including biased bi-allelic transcriptional states. Even Peg13 expression deviated from its expected paternal allele bias in a small number of cells. Although the cell populations consisted of a mosaic of cells with different allelic expression states, as a whole they reflected bulk tissue data. Furthermore, in an attempt to identify potential brain-specific regulatory elements across the Trappc9 locus, we demonstrate tissue-specific and general silencer activities, which might contribute to the regulation of its imprinted expression bias.

DNA Methylation Analysis of Imprinted Genes in the Cortex and Hippocampus of Cross-Fostered Mice Selectively Bred for Increased Voluntary Wheel-Running

We have previously shown that high runner (HR) mice (from a line genetically selected for increased wheel-running behavior) have distinct, genetically based, neurobiological phenotypes as compared with non-selected control (C) mice. However, developmental programming effects during early life, including maternal care and parent-of-origin-dependent expression of imprinted genes, can also contribute to variation in physical activity. Here, we used cross-fostering to address two questions. First, do HR mice have altered DNA methylation profiles of imprinted genes in the brain compared to C mice? Second, does maternal upbringing further modify the DNA methylation status of these imprinted genes? To address these questions, we cross-fostered all offspring at birth to create four experimental groups: C pups to other C dams, HR pups to other HR dams, C pups to HR dams, and HR pups to C dams. Bisulfite sequencing of 16 imprinted genes in the cortex and hippocampus revealed that the HR line had altered DNA methylation patterns of the paternally imprinted genes, Rasgrf1 and Zdbf2, as compared with the C line. Both fostering between the HR and C lines and sex modified the DNA methylation profiles for the paternally expressed genes Mest, Peg3, Igf2, Snrpn, and Impact. Ig-DMR, a gene with multiple paternal and maternal imprinted clusters, was also affected by maternal upbringing and sex. Our results suggest that differential methylation patterns of imprinted genes in the brain could contribute to evolutionary increases in wheel-running behavior and are also dependent on maternal upbringing and sex.

The evolution of widespread recombination suppression on the dwarf hamster (Phodopus) X chromosome

The X chromosome of therian mammals shows strong conservation among distantly related species, limiting insights into the distinct selective processes that have shaped sex chromosome evolution. We constructed a chromosome-scale de novo genome assembly for the Siberian dwarf hamster (Phodopus sungorus), a species reported to show extensive recombination suppression across an entire arm of the X chromosome. Combining a physical genome assembly based on shotgun and long-range proximity ligation sequencing with a dense genetic map, we detected widespread suppression of female recombination across ∼65% of the Phodopus X chromosome. This region of suppressed recombination likely corresponds to the Xp arm, which has previously been shown to be highly heterochromatic. Using additional sequencing data from two closely related species (P. campbelli and P. roborovskii), we show that recombination suppression on Xp appears to be independent of major structural rearrangements. The suppressed Xp arm was enriched for several transposable element families and de-enriched for genes primarily expressed in placenta, but otherwise showed similar gene densities, expression patterns, and rates of molecular evolution when compared to the recombinant Xq arm. Phodopus Xp gene content and order was also broadly conserved relative to the more distantly related rat X chromosome. These data suggest that widespread suppression of recombination has likely evolved through the transient induction of facultative heterochromatin on the Phodopus Xp arm without major changes in chromosome structure or genetic content. Thus, substantial changes in the recombination landscape have so far had relatively subtle influences on patterns of X-linked molecular evolution in these species.

Identification of imprinted genes in the skeletal muscle of newborn piglets by high-throughput sequencing

Imprinted genes - exhibiting parent-specific transcription - play essential roles in the process of mammalian development and growth. Skeletal muscle growth is crucial for meat production. To further understand the role of imprinted genes during the porcine skeletal muscle growth, DNA-seq and RNA-seq were used to explore the characteristics of imprinted genes from porcine reciprocal crosses. A total of 584 545 single-nucleotide variations were discovered in the DNA-seq data of F0 parents, heterozygous in two pig breeds (Yorkshire and Min pigs) but homozygous in each breed. These single-nucleotide variations were used to determine the allelic-specific expression in F1 individuals. Finally, eight paternal expression sites and three maternal expression sites were detected, whereas two paternally expressed imprinted genes (NDN and IGF2) and one maternally expressed imprinted gene (H1-3) were validated by Sanger sequencing. DNA methylation regulates the expression of imprinted genes, and all of the identified imprinted genes in this study were predicted to possess CpG islands. PBX1 and YY1 binding motifs were discovered in the promoter regions of all three imprinted genes, which were candidate elements regulating the transcription of imprinted genes. For these identified imprinted genes, IGF2 and NDN promoted muscle growth whereas H1-3 inhibited cell proliferation, corroborating the 'parental conflict' theory that paternally expressed imprinted genes assisted descendants' growth whereas maternally expressed imprinted genes inhibited it. This study discovered porcine imprinted genes in skeletal muscle and was the first to reveal that H1-3 was expressed by the maternal allele to our knowledge. Our findings provided valuable resources for the potential utilization of imprinted genes in pig breeding.

The Role of Long Non-coding RNAs in Human Imprinting Disorders: Prospective Therapeutic Targets

Genomic imprinting is a term used for an intergenerational epigenetic inheritance and involves a subset of genes expressed in a parent-of-origin-dependent way. Imprinted genes are expressed preferentially from either the paternally or maternally inherited allele. Long non-coding RNAs play essential roles in regulating this allele-specific expression. In several well-studied imprinting clusters, long non-coding RNAs have been found to be essential in regulating temporal- and spatial-specific establishment and maintenance of imprinting patterns. Furthermore, recent insights into the epigenetic pathological mechanisms underlying human genomic imprinting disorders suggest that allele-specific expressed imprinted long non-coding RNAs serve as an upstream regulator of the expression of other protein-coding or non-coding imprinted genes in the same cluster. Aberrantly expressed long non-coding RNAs result in bi-allelic expression or silencing of neighboring imprinted genes. Here, we review the emerging roles of long non-coding RNAs in regulating the expression of imprinted genes, especially in human imprinting disorders, and discuss three strategies targeting the central long non-coding RNA UBE3A-ATS for the purpose of developing therapies for the imprinting disorders Prader-Willi syndrome and Angelman syndrome. In summary, a better understanding of long non-coding RNA-related mechanisms is key to the development of potential therapeutic targets for human imprinting disorders.

An integrated approach in gene-expression landscape profiling to identify housekeeping and tissue-specific genes in cattle

Context High-throughput transcriptome sequencing (RNA-Seq) has been widely applied in cattle studies. Public databases such as the National Center for Biotechnology Information (NCBI) contain large collections of gene expression data from various cattle tissues that can be used in gene expression analysis research Aims This study was conducted to investigate patterns of transcriptome variation across tissues of cattle through large-scale identification of housekeeping genes (i.e. those crucial to maintaining basic cellular activity) and tissue-specific genes in cattle tissues. Methods Using data available in the NCBI Sequence Read Archive database, we analysed 1377 transcriptome data sequences from 60 bovine tissue types, identified tissue-specific and housekeeping genes, and set up a web-based bovine gene expression analysis tool. Key results We found 101 genes widely expressed in almost all tissue and screened out five housekeeping genes: RPL35A, eIF4A2, GAPDH, IPO5 and PAK2. Focusing on 12 major organs, we found 861 genes specifically expressing in these tissues. Furthermore, 187 significantly differentially expressed genes were found among six types of muscle tissues. All expression data were made available at our new website http://cattleExp.org, which can be freely accessed for future gene expression analyses. Conclusions The housekeeping genes and tissue-specific genes identified will provide more information for researchers studying gene expression in cattle. Implications The web-based cattle gene expression analysis tool will make it easy for researchers to access large public datasets. Users can easily access all publicly available RNA data and upload their own RNA-Seq data.

Parental Bias Has Benefits

In this issue, Laukoter et al., 2020 report that parent-of-origin-dependent expression is homogeneous across distinct cortical cell types and within individual populations. Conversely, they observe preferential sensitivity of astrocytes to altered doses of imprinted loci.

A novel method for using RNA-seq data to identify imprinted genes in social Hymenoptera with multiply mated queens

Genomic imprinting results in parent-of-origin-dependent gene expression biased towards either the maternally or paternally derived allele at the imprinted locus. The kinship theory of genomic imprinting argues that this unusual expression pattern can be a manifestation of intra-genomic conflict between the maternally and paternally derived halves of the genome that arises because they are not equally related to the genomes of social partners. The theory thus predicts that imprinting may evolve wherever there are close interactions among asymmetrically related kin. The social Hymenoptera with permanent caste differentiation are suitable candidates for testing the kinship theory because haplodiploid sex determination creates strong relatedness asymmetries and nursing workers interact closely with kin. However, progress in the search for imprinted genes in the social Hymenoptera has been slow, in part because tests for imprinting rely on reciprocal crosses that are impossible in most species. Here, we develop a method to systematically search for imprinting in haplodiploid social insects without crosses, using instead samples of pooled individuals collected from natural colonies. We tested this protocol using data available for the leaf-cutting ant Acromyrmex echinatior, providing the first genome-wide search for imprinting in any ant. Although we identified several genes as potentially imprinted, none of the four genes tested could be verified as imprinted using digital droplet PCR, highlighting the need for higher quality genomic assemblies that accurately map duplicated genes.

Quantifying Genomic Imprinting at Tissue and Cell Resolution in the Brain

Imprinted genes are a group of ~150 genes that are preferentially expressed from one parental allele owing to epigenetic marks asymmetrically distributed on inherited maternal and paternal chromosomes. Altered imprinted gene expression causes human brain disorders such as Prader-Willi and Angelman syndromes and additional rare brain diseases. Research data principally obtained from the mouse model revealed how imprinted genes act in the normal and pathological brain. However, a better understanding of imprinted gene functions calls for building detailed maps of their parent-of-origin-dependent expression and of associated epigenetic signatures. Here we review current methods for quantifying genomic imprinting at tissue and cell resolutions, with a special emphasis on methods to detect parent-of-origin dependent expression and their applications to the brain. We first focus on bulk RNA-sequencing, the main method to detect parent-of-origin-dependent expression transcriptome-wide. We discuss the benefits and caveats of bulk RNA-sequencing and provide a guideline to use it on F1 hybrid mice. We then review methods for detecting parent-of-origin-dependent expression at cell resolution, including single-cell RNA-seq, genetic reporters, and molecular probes. Finally, we provide an overview of single-cell epigenomics technologies that profile additional features of genomic imprinting, including DNA methylation, histone modifications and chromatin conformation and their combination into sc-multimodal omics approaches, which are expected to yield important insights into genomic imprinting in individual brain cells.

Trappc9 deficiency causes parent-of-origin dependent microcephaly and obesity

Some imprinted genes exhibit parental origin specific expression bias rather than being transcribed exclusively from one copy. The physiological relevance of this remains poorly understood. In an analysis of brain-specific allele-biased expression, we identified that Trappc9, a cellular trafficking factor, was expressed predominantly (~70%) from the maternally inherited allele. Loss-of-function mutations in human TRAPPC9 cause a rare neurodevelopmental syndrome characterized by microcephaly and obesity. By studying Trappc9 null mice we discovered that homozygous mutant mice showed a reduction in brain size, exploratory activity and social memory, as well as a marked increase in body weight. A role for Trappc9 in energy balance was further supported by increased ad libitum food intake in a child with TRAPPC9 deficiency. Strikingly, heterozygous mice lacking the maternal allele (70% reduced expression) had pathology similar to homozygous mutants, whereas mice lacking the paternal allele (30% reduction) were phenotypically normal. Taken together, we conclude that Trappc9 deficient mice recapitulate key pathological features of TRAPPC9 mutations in humans and identify a role for Trappc9 and its imprinting in controlling brain development and metabolism.

Parental allele-specific protein expression in single cells In vivo

Allelic expression from each parent-of-origin is important as a backup and to ensure that enough protein products of a gene are produced. Thus far, it is not known how each cell throughout a tissue differs in parental allele expression at the level of protein synthesis. Here, we measure the expression of the Ribosomal protein L13a (Rpl13a) from both parental alleles simultaneously in single cells in the living animal. We use genome-edited Drosophila that have a quantitative reporter of protein synthesis inserted into the endogenous Rpl13a locus. We find that individual cells can have large (>10-fold) differences in protein expression between the two parental alleles. Cells can produce protein from only one allele oftentimes, and time-lapse imaging of protein production from each parental allele in each cell showed that the imbalance in expression from one parental allele over the other can invert over time. We also identify the histone methyltransferase EHMT to be involved in the protein synthesis dynamics within cells.

RNA-Seq in 296 phased trios provides a high-resolution map of genomic imprinting

Background

Identification of imprinted genes, demonstrating a consistent preference towards the paternal or maternal allelic expression, is important for the understanding of gene expression regulation during embryonic development and of the molecular basis of developmental disorders with a parent-of-origin effect. Combining allelic analysis of RNA-Seq data with phased genotypes in family trios provides a powerful method to detect parent-of-origin biases in gene expression.

Results

We report findings in 296 family trios from two large studies: 165 lymphoblastoid cell lines from the 1000 Genomes Project and 131 blood samples from the Genome of the Netherlands (GoNL) participants. Based on parental haplotypes, we identified > 2.8 million transcribed heterozygous SNVs phased for parental origin and developed a robust statistical framework for measuring allelic expression. We identified a total of 45 imprinted genes and one imprinted unannotated transcript, including multiple imprinted transcripts showing incomplete parental expression bias that was located adjacent to strongly imprinted genes. For example, PXDC1, a gene which lies adjacent to the paternally expressed gene FAM50B, shows a 2:1 paternal expression bias. Other imprinted genes had promoter regions that coincide with sites of parentally biased DNA methylation identified in the blood from uniparental disomy (UPD) samples, thus providing independent validation of our results. Using the stranded nature of the RNA-Seq data in lymphoblastoid cell lines, we identified multiple loci with overlapping sense/antisense transcripts, of which one is expressed paternally and the other maternally. Using a sliding window approach, we searched for imprinted expression across the entire genome, identifying a novel imprinted putative lncRNA in 13q21.2. Overall, we identified 7 transcripts showing parental bias in gene expression which were not reported in 4 other recent RNA-Seq studies of imprinting.

Conclusions

Our methods and data provide a robust and high-resolution map of imprinted gene expression in the human genome.

Electronic supplementary material

The online version of this article (10.1186/s12915-019-0674-0) contains supplementary material, which is available to authorized users.

Technical advances contribute to the study of genomic imprinting

Genomic imprinting in mammals was discovered over 30 years ago through elegant embryological and genetic experiments in mice. Imprinted genes show a monoallelic and parent of origin-specific expression pattern; the development of techniques that can distinguish between expression from maternal and paternal chromosomes in mice, combined with high-throughput strategies, has allowed for identification of many more imprinted genes, most of which are conserved in humans. Undoubtedly, technical progress has greatly promoted progress in the field of genomic imprinting. Here, we summarize the techniques used to discover imprinted genes, identify new imprinted genes, define imprinting regulation mechanisms, and study imprinting functions.

A statistical normalization method and differential expression analysis for RNA-seq data between different species

Background

High-throughput techniques bring novel tools and also statistical challenges to genomic research. Identifying genes with differential expression between different species is an effective way to discover evolutionarily conserved transcriptional responses. To remove systematic variation between different species for a fair comparison, normalization serves as a crucial pre-processing step that adjusts for the varying sample sequencing depths and other confounding technical effects.

Results

In this paper, we propose a scale based normalization (SCBN) method by taking into account the available knowledge of conserved orthologous genes and by using the hypothesis testing framework. Considering the different gene lengths and unmapped genes between different species, we formulate the problem from the perspective of hypothesis testing and search for the optimal scaling factor that minimizes the deviation between the empirical and nominal type I errors.

Conclusions

Simulation studies show that the proposed method performs significantly better than the existing competitor in a wide range of settings. An RNA-seq dataset of different species is also analyzed and it coincides with the conclusion that the proposed method outperforms the existing method. For practical applications, we have also developed an R package named “SCBN”, which is freely available at http://www.bioconductor.org/packages/devel/bioc/html/SCBN.html.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2745-1) contains supplementary material, which is available to authorized users.

Genomic Imprinting and Physiological Processes in Mammals

Complex multicellular organisms, such as mammals, express two complete sets of chromosomes per nucleus, combining the genetic material of both parents. However, epigenetic studies have demonstrated violations to this rule that are necessary for mammalian physiology; the most notable parental allele expression phenomenon is genomic imprinting. With the identification of endogenous imprinted genes, genomic imprinting became well-established as an epigenetic mechanism in which the expression pattern of a parental allele influences phenotypic expression. The expanding study of genomic imprinting is revealing a significant impact on brain functions and associated diseases. Here, we review key milestones in the field of imprinting and discuss mechanisms and systems in which imprinted genes exert a significant role.

The emerging landscape of in vitro and in vivo epigenetic allelic effects

Epigenetic mechanisms that cause maternally and paternally inherited alleles to be expressed differently in offspring have the potential to radically change our understanding of the mechanisms that shape disease susceptibility, phenotypic variation, cell fate, and gene expression. However, the nature and prevalence of these effects in vivo have been unclear and are debated. Here, I consider major new studies of epigenetic allelic effects in cell lines and primary cells and in vivo. The emerging picture is that these effects take on diverse forms, and this review attempts to clarify the nature of the different forms that have been uncovered for genomic imprinting and random monoallelic expression (RME). I also discuss apparent discrepancies between in vitro and in vivo studies. Importantly, multiple studies suggest that allelic effects are prevalent and can be developmental stage- and cell type-specific. I propose some possible functions and consider roles for allelic effects within the broader context of gene regulatory networks, cellular diversity, and plasticity. Overall, the field is ripe for discovery and is in need of mechanistic and functional studies.

Mapping the mouse Allelome reveals tissue-specific regulation of allelic expression

To determine the dynamics of allelic-specific expression during mouse development, we analyzed RNA-seq data from 23 F1 tissues from different developmental stages, including 19 female tissues allowing X chromosome inactivation (XCI) escapers to also be detected. We demonstrate that allelic expression arising from genetic or epigenetic differences is highly tissue-specific. We find that tissue-specific strain-biased gene expression may be regulated by tissue-specific enhancers or by post-transcriptional differences in stability between the alleles. We also find that escape from X-inactivation is tissue-specific, with leg muscle showing an unexpectedly high rate of XCI escapers. By surveying a range of tissues during development, and performing extensive validation, we are able to provide a high confidence list of mouse imprinted genes including 18 novel genes. This shows that cluster size varies dynamically during development and can be substantially larger than previously thought, with the Igf2r cluster extending over 10 Mb in placenta.

DOI:http://dx.doi.org/10.7554/eLife.25125.001

Assessment of imprinting- and genetic variation-dependent monoallelic expression using reciprocal allele descendants between human family trios

Genomic imprinting is an important epigenetic process that silences one of the parentally-inherited alleles of a gene and thereby exhibits allelic-specific expression (ASE). Detection of human imprinting events is hampered by the infeasibility of the reciprocal mating system in humans and the removal of ASE events arising from non-imprinting factors. Here, we describe a pipeline with the pattern of reciprocal allele descendants (RADs) through genotyping and transcriptome sequencing data across independent parent-offspring trios to discriminate between varied types of ASE (e.g., imprinting, genetic variation-dependent ASE, and random monoallelic expression (RME)). We show that the vast majority of ASE events are due to sequence-dependent genetic variant, which are evolutionarily conserved and may themselves play a cis-regulatory role. Particularly, 74% of non-RAD ASE events, even though they exhibit ASE biases toward the same parentally-inherited allele across different individuals, are derived from genetic variation but not imprinting. We further show that the RME effect may affect the effectiveness of the population-based method for detecting imprinting events and our pipeline can help to distinguish between these two ASE types. Taken together, this study provides a good indicator for categorization of different types of ASE, opening up this widespread and complex mechanism for comprehensive characterization.

Conservation of Repeats at the Mammalian KCNQ1OT1-CDKN1C Region Suggests a Role in Genomic Imprinting

KCNQ1OT1 is located in the region with the highest number of genes showing genomic imprinting, but the mechanisms controlling the genes under its influence have not been fully elucidated. Therefore, we conducted a comparative analysis of the KCNQ1/KCNQ1OT1-CDKN1C region to study its conservation across the best assembled eutherian mammalian genomes sequenced to date and analyzed potential elements that may be implicated in the control of genomic imprinting in this region. The genomic features in these regions from human, mouse, cattle, and dog show a higher number of genes and CpG islands (detected using cpgplot from EMBOSS), but lower number of repetitive elements (including short interspersed nuclear elements and long interspersed nuclear elements), compared with their whole chromosomes (detected by RepeatMasker). The KCNQ1OT1-CDKN1C region contains the highest number of conserved noncoding sequences (CNS) among mammals, where we found 16 regions containing about 38 different highly conserved repetitive elements (using mVista), such as LINE1 elements: L1M4, L1MB7, HAL1, L1M4a, L1Med, and an LTR element: MLT1H. From these elements, we found 74 CNS showing high sequence identity (>70%) between human, cattle, and mouse, from which we identified 13 motifs (using Multiple Em for Motif Elicitation/Motif Alignment and Search Tool) with a significant probability of occurrence, 3 of which were the most frequent and were used to find transcription factor-binding sites. We detected several transcription factors (using JASPAR suite) from the families SOX, FOX, and GATA. A phylogenetic analysis of these CNS from human, marmoset, mouse, rat, cattle, dog, horse, and elephant shows branches with high levels of support and very similar phylogenetic relationships among these groups, confirming previous reports. Our results suggest that functional DNA elements identified by comparative genomics in a region densely populated with imprinted mammalian genes may be related to the regulation of imprinted gene expression.

In Vitro Corticogenesis from Embryonic Stem Cells Recapitulates the In Vivo Epigenetic Control of Imprinted Gene Expression

In vitro corticogenesis from embryonic stem cells (ESCs) is an attractive model of cortical development and a promising tool for cortical therapy. It is unknown to which extent epigenetic mechanisms crucial for cortex development and function, such as parental genomic imprinting, are recapitulated by in vitro corticogenesis. Here, using genome-wide transcriptomic and methylation analyses on hybrid mouse tissues and cells, we find a high concordance of imprinting status between in vivo and ESC-derived cortices. Notably, in vitro corticogenesis strictly reproduced the in vivo parent-of-origin-dependent expression of 41 imprinted genes (IGs), including Mest and Cdkn1c known to control corticogenesis. Parent-of-origin-dependent DNA methylation was also conserved at 14 of 18 imprinted differentially methylated regions. The least concordant imprinted locus was Gpr1-Zdbf2, where the aberrant bi-allelic expression of Zdbf2 and Adam23 was concomitant with a gain of methylation on the maternal allele in vitro. Combined, our data argue for a broad conservation of the epigenetic mechanisms at imprinted loci in cortical cells derived from ESCs. We propose that in vitro corticogenesis helps to define the still poorly understood mechanisms that regulate imprinting in the brain and the roles of IGs in cortical development.

Monoallelic Gene Expression in Mammals

Monoallelic expression not due to cis-regulatory sequence polymorphism poses an intriguing problem in epigenetics because it requires the unequal treatment of two segments of DNA that are present in the same nucleus and that can indeed have absolutely identical sequences. Here, I focus on a few recent developments in the field of monoallelic expression that are of particular interest and raise interesting questions for future work. One development is regarding analyses of imprinted genes, in which recent work suggests the possibility that intriguing networks of imprinted genes exist and are important for genetic and physiological studies. Another issue that has been raised in recent years by a number of publications is the question of how skewed allelic expression should be for it to be designated as monoallelic expression and, further, what methods are appropriate or inappropriate for analyzing genomic data to examine allele-specific expression. Perhaps the most exciting recent development in mammalian monoallelic expression is a clever and carefully executed analysis of genetic diversity of autosomal genes subject to random monoallelic expression (RMAE), which provides compelling evidence for distinct evolutionary forces acting on random monoallelically expressed genes.

A Hypothesis Testing Based Method for Normalization and Differential Expression Analysis of RNA-Seq Data

Next-generation sequencing technologies have made RNA sequencing (RNA-seq) a popular choice for measuring gene expression level. To reduce the noise of gene expression measures and compare them between several conditions or samples, normalization is an essential step to adjust for varying sample sequencing depths and other unwanted technical effects. In this paper, we develop a novel global scaling normalization method by employing the available knowledge of housekeeping genes. We formulate the problem from the hypothesis testing perspective and find an optimal scaling factor that minimizes the deviation between the empirical and the nominal type I error. Applying our approach to various simulation studies and real examples, we demonstrate that it is more accurate and robust than the state-of-the-art alternatives in detecting differentially expression genes.

The abundance of cis-acting loci leading to differential allele expression in F1 mice and their relationship to loci harboring genes affecting complex traits

Background

Genome-wide surveys have detected cis-acting quantitative trait loci altering levels of RNA transcripts (RNA-eQTLs) by associating SNV alleles to transcript levels. However, the sensitivity and specificity of detection of cis- expression quantitative trait loci (eQTLs) by genetic approaches, reliant as it is on measurements of transcript levels in recombinant inbred strains or offspring from arranged crosses, is unknown, as is their relationship to QTL’s for complex phenotypes.

Results

We used transcriptome-wide differential allele expression (DAE) to detect cis-eQTLs in forebrain and kidney from reciprocal crosses between three mouse inbred strains, 129S1/SvlmJ, DBA/2J, and CAST/EiJ and C57BL/6 J. Two of these crosses were previously characterized for cis-eQTLs and QTLs for various complex phenotypes by genetic analysis of recombinant inbred (RI) strains. 5.4 %, 1.9 % and 1.5 % of genes assayed in forebrain of B6/129SF1, B6/DBAF1, and B6/CASTF1 mice, respectively, showed differential allelic expression, indicative of cis-acting alleles at these genes. Moreover, the majority of DAE QTLs were observed to be tissue-specific with only a small fraction showing cis-effects in both tissues. Comparing DAE QTLs in F1 mice to cis-eQTLs previously mapped in RI strains we observed that many of the cis-eQTLs were not confirmed by DAE. Additionally several novel DAE-QTLs not identified as cis-eQTLs were identified suggesting that there are differences in sensitivity and specificity for QTL detection between the two methodologies. Strain specific DAE QTLs in B6/DBAF1 mice were located in excess at candidate genes for alcohol use disorders, seizures, and angiogenesis previously implicated by genetic linkage in C57BL/6J × DBA/2JF2 mice or BXD RI strains.

Conclusions

Via a survey for differential allele expression in F1 mice, a substantial proportion of genes were found to have alleles altering expression in cis-acting fashion. Comparing forebrain and kidney, many or most of these alleles were tissue-specific in action. The identification of strain specific DAE QTLs, can assist in assessment of candidate genes located within the large intervals associated with trait QTLs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2922-9) contains supplementary material, which is available to authorized users.

Global assessment of imprinted gene expression in the bovine conceptus by next generation sequencing

Genomic imprinting is an epigenetic mechanism that leads to parental-allele-specific gene expression. Approximately 150 imprinted genes have been identified in humans and mice but less than 30 have been described as imprinted in cattle. For the purpose of de novo identification of imprinted genes in bovine, we determined global monoallelic gene expression in brain, skeletal muscle, liver, kidney and placenta of day ∼105 Bos taurus indicus × Bos taurus taurus F₁ conceptuses using RNA sequencing. To accomplish this, we developed a bioinformatics pipeline to identify parent-specific single nucleotide polymorphism alleles after filtering adenosine to inosine (A-to-I) RNA editing sites. We identified 53 genes subject to monoallelic expression. Twenty three are genes known to be imprinted in the cow and an additional 7 have previously been characterized as imprinted in human and/or mouse that have not been reported as imprinted in cattle. Of the remaining 23 genes, we found that 10 are uncharacterized or unannotated transcripts located in known imprinted clusters, whereas the other 13 genes are distributed throughout the bovine genome and are not close to any known imprinted clusters. To exclude potential cis-eQTL effects on allele expression, we corroborated the parental specificity of monoallelic expression in day 86 Bos taurus taurus × Bos taurus taurus conceptuses and identified 8 novel bovine imprinted genes. Further, we identified 671 candidate A-to-I RNA editing sites and describe random X-inactivation in day 15 bovine extraembryonic membranes. Our results expand the imprinted gene list in bovine and demonstrate that monoallelic gene expression can be the result of cis-eQTL effects.

The evolving landscape of imprinted genes in humans and mice: Conflict among alleles, genes, tissues, and kin

Three recent genome-wide studies in mice and humans have produced the most definitive map to date of genomic imprinting (gene expression that depends on parental origin) by incorporating multiple tissue types and developmental stages. Here, we explore the results of these studies in light of the kinship theory of genomic imprinting, which predicts that imprinting evolves due to differential genetic relatedness between maternal and paternal relatives. The studies produce a list of imprinted genes with around 120-180 in mice and ~100 in humans. The studies agree on broad patterns across mice and humans including the complex patterns of imprinted expression at loci like Igf2 and Grb10. We discuss how the kinship theory provides a powerful framework for hypotheses that can explain these patterns. Finally, since imprinting is rare in the genome despite predictions from the kinship theory that it might be common, we discuss evolutionary factors that could favor biallelic expression.

Testing the kinship theory of intragenomic conflict in honey bees (Apis mellifera)

Sexual reproduction brings genes from two parents (matrigenes and patrigenes) together into one individual. These genes, despite being unrelated, should show nearly perfect cooperation because each gains equally through the production of offspring. However, an individual's matrigenes and patrigenes can have different probabilities of being present in other relatives, so kin selection could act on them differently. Such intragenomic conflict could be implemented by partial or complete silencing (imprinting) of an allele by one of the parents. Evidence supporting this theory is seen in offspring-mother interactions, with patrigenes favoring acquisition of more of the mother's resources if some of the costs fall on half-siblings who do not share the patrigene. The kinship theory of intragenomic conflict is little tested in other contexts, but it predicts that matrigene-patrigene conflict may be rife in social insects. We tested the hypothesis that honey bee worker reproduction is promoted more by patrigenes than matrigenes by comparing across nine reciprocal crosses of two distinct genetic stocks. As predicted, hybrid workers show reproductive trait characteristics of their paternal stock, (indicating enhanced activity of the patrigenes on these traits), greater patrigenic than matrigenic expression, and significantly increased patrigenic-biased expression in reproductive workers. These results support both the general prediction that matrigene-patrigene conflict occurs in social insects and the specific prediction that honey bee worker reproduction is driven more by patrigenes. The success of these predictions suggests that intragenomic conflict may occur in many contexts where matrigenes and patrigenes have different relatednesses to affected kin.

Intersection of genetics and epigenetics in monozygotic twin genomes

As a final function of various epigenetic mechanisms, chromatin regulation is a transcription control process that especially demonstrates active interaction with genetic elements. Thus, chromatin structure has become a principal focus in recent genomics researches that strive to characterize regulatory functions of DNA variants related to diseases or other traits. Although researchers have been focusing on DNA methylation when studying monozygotic (MZ) twins, a great model in epigenetics research, interactions between genetics and epigenetics in chromatin level are expected to be an imperative research trend in the future. In this review, we discuss how the genome, epigenome, and transcriptome of MZ twins can be studied in an integrative manner from this perspective.

Noncanonical Genomic Imprinting Effects in Offspring

Here, we describe an RNA-sequencing (RNA-seq)-based approach that accurately detects even modest maternal or paternal allele expression biases at the tissue level, which we call noncanonical genomic imprinting effects. We profile imprinting in the arcuate nucleus (ARN) and dorsal raphe nucleus of the female mouse brain as well as skeletal muscle (mesodermal) and liver (endodermal). Our study uncovers hundreds of noncanonical autosomal and X-linked imprinting effects. Noncanonical imprinting is highly tissue-specific and enriched in the ARN, but rare in the liver. These effects are reproducible across different genetic backgrounds and associated with allele-specific chromatin. Using in situ hybridization for nascent RNAs, we discover that autosomal noncanonical imprinted genes with a tissue-level allele bias exhibit allele-specific expression effects in subpopulations of neurons in the brain in vivo. We define noncanonical imprinted genes that regulate monoamine signaling and determine that these effects influence the impact of inherited mutations on offspring behavior.

Allelome.PRO, a pipeline to define allele-specific genomic features from high-throughput sequencing data

Detecting allelic biases from high-throughput sequencing data requires an approach that maximises sensitivity while minimizing false positives. Here, we present Allelome.PRO, an automated user-friendly bioinformatics pipeline, which uses high-throughput sequencing data from reciprocal crosses of two genetically distinct mouse strains to detect allele-specific expression and chromatin modifications. Allelome.PRO extends approaches used in previous studies that exclusively analyzed imprinted expression to give a complete picture of the ‘allelome’ by automatically categorising the allelic expression of all genes in a given cell type into imprinted, strain-biased, biallelic or non-informative. Allelome.PRO offers increased sensitivity to analyze lowly expressed transcripts, together with a robust false discovery rate empirically calculated from variation in the sequencing data. We used RNA-seq data from mouse embryonic fibroblasts from F1 reciprocal crosses to determine a biologically relevant allelic ratio cutoff, and define for the first time an entire allelome. Furthermore, we show that Allelome.PRO detects differential enrichment of H3K4me3 over promoters from ChIP-seq data validating the RNA-seq results. This approach can be easily extended to analyze histone marks of active enhancers, or transcription factor binding sites and therefore provides a powerful tool to identify candidate cis regulatory elements genome wide.

Quantitative and functional interrogation of parent-of-origin allelic expression biases in the brain

The maternal and paternal genomes play different roles in mammalian brains as a result of genomic imprinting, an epigenetic regulation leading to differential expression of the parental alleles of some genes. Here we investigate genomic imprinting in the cerebellum using a newly developed Bayesian statistical model that provides unprecedented transcript-level resolution. We uncover 160 imprinted transcripts, including 41 novel and independently validated imprinted genes. Strikingly, many genes exhibit parentally biased--rather than monoallelic--expression, with different magnitudes according to age, organ, and brain region. Developmental changes in parental bias and overall gene expression are strongly correlated, suggesting combined roles in regulating gene dosage. Finally, brain-specific deletion of the paternal, but not maternal, allele of the paternally-biased Bcl-x, (Bcl2l1) results in loss of specific neuron types, supporting the functional significance of parental biases. These findings reveal the remarkable complexity of genomic imprinting, with important implications for understanding the normal and diseased brain.

Genomic imprinting: A missing piece of the Multiple Sclerosis puzzle?

Evidence for parent-of-origin effects in complex diseases such as Multiple Sclerosis (MS) strongly suggests a role for epigenetic mechanisms in their pathogenesis. In this review, we describe the importance of accounting for parent-of-origin when identifying new risk variants for complex diseases and discuss how genomic imprinting, one of the best-characterized epigenetic mechanisms causing parent-of-origin effects, may impact etiology of complex diseases. While the role of imprinted genes in growth and development is well established, the contribution and molecular mechanisms underlying the impact of genomic imprinting in immune functions and inflammatory diseases are still largely unknown. Here we discuss emerging roles of imprinted genes in the regulation of inflammatory responses with a particular focus on the Dlk1 cluster that has been implicated in etiology of experimental MS-like disease and Type 1 Diabetes. Moreover, we speculate on the potential wider impact of imprinting via the action of imprinted microRNAs, which are abundantly present in the Dlk1 locus and predicted to fine-tune important immune functions. Finally, we reflect on how unrelated imprinted genes or imprinted genes together with non-imprinted genes can interact in so-called imprinted gene networks (IGN) and suggest that IGNs could partly explain observed parent-of-origin effects in complex diseases. Unveiling the mechanisms of parent-of-origin effects is therefore likely to teach us not only about the etiology of complex diseases but also about the unknown roles of this fascinating phenomenon underlying uneven genetic contribution from our parents. This article is part of a Directed Issue entitled: Epigenetics dynamics in development and disease.

Genetic conflict reflected in tissue-specific maps of genomic imprinting in human and mouse

Genomic imprinting is an epigenetic process that restricts gene expression to either the maternally or paternally inherited allele. Many theories have been proposed to explain its evolutionary origin, but understanding has been limited by a paucity of data mapping the breadth and dynamics of imprinting within any organism. We generated an atlas of imprinting spanning 33 mouse and 45 human developmental stages and tissues. Nearly all imprinted genes were imprinted in early development and either retained their parent-of-origin expression in adults or lost it completely. Consistent with an evolutionary signature of parental conflict, imprinted genes were enriched for coexpressed pairs of maternally and paternally expressed genes, showed accelerated expression divergence between human and mouse, and were more highly expressed than their non-imprinted orthologs in other species. Our approach demonstrates a general framework for the discovery of imprinting in any species and sheds light on the causes and consequences of genomic imprinting in mammals.

Analyses of allele-specific gene expression in highly divergent mouse crosses identifies pervasive allelic imbalance

Complex human traits are influenced by variation in regulatory DNA through mechanisms that are not fully understood. Since regulatory elements are conserved between humans and mice, a thorough annotation of cis regulatory variants in mice could aid in this process. Here we provide a detailed portrait of mouse gene expression across multiple tissues in a three-way diallel. Greater than 80% of mouse genes have cis regulatory variation. These effects influence complex traits and usually extend to the human ortholog. Further, we estimate that at least one in every thousand SNPs creates a cis regulatory effect. We also observe two types of parent-of-origin effects, including classical imprinting and a novel, global allelic imbalance in favor of the paternal allele. We conclude that, as with humans, pervasive regulatory variation influences complex genetic traits in mice and provide a new resource toward understanding the genetic control of transcription in mammals.

The genetic architecture of the genome-wide transcriptional response to ER stress in the mouse

Endoplasmic reticulum (ER) stress occurs when misfolded proteins accumulate in the ER. The cellular response to ER stress involves complex transcriptional and translational changes, important to the survival of the cell. ER stress is a primary cause and a modifier of many human diseases. A first step to understanding how the ER stress response impacts human disease is to determine how the transcriptional response to ER stress varies among individuals. The genetic diversity of the eight mouse Collaborative Cross (CC) founder strains allowed us to determine how genetic variation impacts the ER stress transcriptional response. We used tunicamycin, a drug commonly used to induce ER stress, to elicit an ER stress response in mouse embryonic fibroblasts (MEFs) derived from the CC founder strains and measured their transcriptional responses. We identified hundreds of genes that differed in response to ER stress across these genetically diverse strains. Strikingly, inflammatory response genes differed most between strains; major canonical ER stress response genes showed relatively invariant responses across strains. To uncover the genetic architecture underlying these strain differences in ER stress response, we measured the transcriptional response to ER stress in MEFs derived from a subset of F1 crosses between the CC founder strains. We found a unique layer of regulatory variation that is only detectable under ER stress conditions. Over 80% of the regulatory variation under ER stress derives from cis-regulatory differences. This is the first study to characterize the genetic variation in ER stress transcriptional response in the laboratory mouse. Our findings indicate that the ER stress transcriptional response is highly variable among strains and arises from genetic variation in individual downstream response genes, rather than major signaling transcription factors. These results have important implications for understanding how genetic variation impacts the ER stress response, an important component of many human diseases.

What does genetics tell us about imprinting and the placenta connection?

Genomic imprinting is an epigenetic gene silencing phenomenon that is specific to eutherians in the vertebrate lineage. The acquisition of both placentation and genomic imprinting has spurred interest in the possible evolutionary link for many years. In this review we examine the genetic evidence and find that while many imprinted domains are anchored by genes required for proper placenta development in a parent of origin fashion, an equal number of imprinted genes have no apparent function that depends on imprinting. Examination of recent data from studies of molecular and genetic mechanisms points to a maternal control of the selection and maintenance of imprint marks, reinforcing the importance of the oocyte in the healthy development of the placenta and fetus.

Characterization of X chromosome inactivation using integrated analysis of whole-exome and mRNA sequencing

In females, X chromosome inactivation (XCI) is an epigenetic, gene dosage compensatory mechanism by inactivation of one copy of X in cells. Random XCI of one of the parental chromosomes results in an approximately equal proportion of cells expressing alleles from either the maternally or paternally inherited active X, and is defined by the XCI ratio. Skewed XCI ratio is suggestive of non-random inactivation, which can play an important role in X-linked genetic conditions. Current methods rely on indirect, semi-quantitative DNA methylation-based assay to estimate XCI ratio. Here we report a direct approach to estimate XCI ratio by integrated, family-trio based whole-exome and mRNA sequencing using phase-by-transmission of alleles coupled with allele-specific expression analysis. We applied this method to in silico data and to a clinical patient with mild cognitive impairment but no clear diagnosis or understanding molecular mechanism underlying the phenotype. Simulation showed that phased and unphased heterozygous allele expression can be used to estimate XCI ratio. Segregation analysis of the patient's exome uncovered a de novo, interstitial, 1.7 Mb deletion on Xp22.31 that originated on the paternally inherited X and previously been associated with heterogeneous, neurological phenotype. Phased, allelic expression data suggested an 83∶20 moderately skewed XCI that favored the expression of the maternally inherited, cytogenetically normal X and suggested that the deleterious affect of the de novo event on the paternal copy may be offset by skewed XCI that favors expression of the wild-type X. This study shows the utility of integrated sequencing approach in XCI ratio estimation.

SNP-guided identification of monoallelic DNA-methylation events from enrichment-based sequencing data

Monoallelic gene expression is typically initiated early in the development of an organism. Dysregulation of monoallelic gene expression has already been linked to several non-Mendelian inherited genetic disorders. In humans, DNA-methylation is deemed to be an important regulator of monoallelic gene expression, but only few examples are known. One important reason is that current, cost-affordable truly genome-wide methods to assess DNA-methylation are based on sequencing post-enrichment. Here, we present a new methodology based on classical population genetic theory, i.e. the Hardy–Weinberg theorem, that combines methylomic data from MethylCap-seq with associated SNP profiles to identify monoallelically methylated loci. Applied on 334 MethylCap-seq samples of very diverse origin, this resulted in the identification of 80 genomic regions featured by monoallelic DNA-methylation. Of these 80 loci, 49 are located in genic regions of which 25 have already been linked to imprinting. Further analysis revealed statistically significant enrichment of these loci in promoter regions, further establishing the relevance and usefulness of the method. Additional validation was done using both 14 whole-genome bisulfite sequencing data sets and 16 mRNA-seq data sets. Importantly, the developed approach can be easily applied to other enrichment-based sequencing technologies, like the ChIP-seq-based identification of monoallelic histone modifications.

Known unknowns for allele-specific expression and genomic imprinting effects

Recent studies have provided evidence for non-canonical imprinting effects that are associated with allele-specific expression biases at the tissue level in mice. These imprinting effects have features that are distinct from canonical imprinting effects that involve allele silencing. Here, I discuss some of the evidence for non-canonical imprinting effects in the context of random X-inactivation and epigenetic allele-specific expression effects on the autosomes. I propose several mechanisms that may underlie non-canonical imprinting effects and outline future directions and approaches to study these effects at the cellular level in vivo. The growing evidence for complex allele-specific expression effects that are cell- and developmental stage-specific has opened a new frontier for study. Currently, the function of these effects and the underlying regulatory mechanisms are largely unknown.

Epigenetics in adaptive evolution and development: the interplay between evolving species and epigenetic mechanisms: extract from Trygve Tollefsbol (ed.) (2011) Handbook of epigenetics--the new molecular and medical genetics. Chapter 26. Amsterdam, USA: Elsevier, pp. 423-446

By comparing epigenetics of current species with fossil records across evolutionary transitions, we can gauge the moment of emergence of some novel mechanisms in evolution, and recognize that epigenetic mechanisms have a bearing on mutation. Understanding the complexity and changeability of these mechanisms, as well as the changes they can effect, is both fascinating and of vital practical benefit. Our most serious pandemics of so-called 'non-communicable' diseases - mental and cardiovascular disorders, obesity and diabetes, rooted in the 'metabolic syndrome' - are evidently related to effects on our evolutionary mechanisms of agricultural and food industrialization, modern lifestyle and diet. Pollution affects us directly as well as indirectly by its destruction of ecologically essential biosystems. Evidently such powerful conditions of existence have epigenetic effects on both our health and our continuing evolution. Such effects are most profound during reproductive and developmental processes, when levels of hormones, as affected by stress particularly, may be due to modern cultures in childbearing such as excessive intervention, separation, maternal distress and disruption of bonding. Mechanisms of genomic imprinting seem likely to throw light on problems in assisted reproductive technology, among other transgenerational effects.

Nonparametric method for detecting imprinting effect using all members of general pedigrees with missing data

Imprinting effects can lead to parent-of-origin patterns in complex human diseases. For a diallelic marker locus, Pedigree Parental-Asymmetry Test (PPAT) and its extension MCPPAT using pedigrees allowing for missing genotypes are simple and powerful for detecting imprinting effects. However, these approaches only take affected offspring into consideration, thus not making full use of the data available. In this paper, we propose Monte Carlo Pedigree Parental-Asymmetry Test using both affected and unaffected (MCPPATu) offsprings, which allows for missing genotypes through Monte Carlo sampling. Simulation studies demonstrate that MCPPATu controls the empirical type I error rate well under the null hypotheses of no parent-of-origin effects. It is also demonstrated that the use of additional information from unaffected offspring and partially observed genotypes in the analysis can greatly improve the statistical power. Indeed, for common diseases, MCPPATu is much more powerful than MCPPAT when all genotypes are observed and the power improvement is even greater when there is missing data. For rarer diseases, there are still substantial power gains with the inclusion of unaffected offspring, although the gains are less impressive compared with those for more common diseases.

An iteration normalization and test method for differential expression analysis of RNA-seq data

Background

Next generation sequencing technologies are powerful new tools for investigating a wide range of biological and medical questions. Statistical and computational methods are key to analyzing massive and complex sequencing data. In order to derive gene expression measures and compare these measures across samples or libraries, we first need to normalize read counts to adjust for varying sample sequencing depths and other potentially technical effects.

Results

In this paper, we develop a normalization method based on iterating median of M-values (IMM) for detecting the differentially expressed (DE) genes. Compared to a previous approach TMM, the IMM method improves the accuracy of DE detection. Simulation studies show that the IMM method outperforms other methods for the sample normalization. We also look into the real data and find that the genes detected by IMM but not by TMM are much more accurate than the genes detected by TMM but not by IMM. What’s more, we discovered that gene UNC5C is highly associated with kidney cancer and so on.

Using next-generation RNA sequencing to identify imprinted genes

Genomic imprinting is manifested as differential allelic expression (DAE) depending on the parent-of-origin. The most direct way to identify imprinted genes is to directly score the DAE in a context where one can identify which parent transmitted each allele. Because many genes display DAE, simply scoring DAE in an individual is not sufficient to identify imprinted genes. In this paper, we outline many technical aspects of a scheme for identification of imprinted genes that makes use of RNA sequencing (RNA-seq) from tissues isolated from F1 offspring derived from the pair of reciprocal crosses. Ideally, the parental lines are from two inbred strains that are not closely related to each other. Aspects of tissue purity, RNA extraction, library preparation and bioinformatic inference of imprinting are all covered. These methods have already been applied in a number of organisms, and one of the most striking results is the evolutionary fluidity with which novel imprinted genes are gained and lost within genomes. The general methodology is also applicable to a wide range of other biological problems that require quantification of allele-specific expression using RNA-seq, such as cis-regulation of gene expression, X chromosome inactivation and random monoallelic expression.

MetaImprint: an information repository of mammalian imprinted genes

Genomic imprinting is a complex genetic and epigenetic phenomenon that plays important roles in mammalian development and diseases. Mammalian imprinted genes have been identified widely by experimental strategies or predicted using computational methods. Systematic information for these genes would be necessary for the identification of novel imprinted genes and the analysis of their regulatory mechanisms and functions. Here, a well-designed information repository, MetaImprint (http://bioinfo.hrbmu.edu.cn/MetaImprint), is presented, which focuses on the collection of information concerning mammalian imprinted genes. The current version of MetaImprint incorporates 539 imprinted genes, including 255 experimentally confirmed genes, and their detailed research courses from eight mammalian species. MetaImprint also hosts genome-wide genetic and epigenetic information of imprinted genes, including imprinting control regions, single nucleotide polymorphisms, non-coding RNAs, DNA methylation and histone modifications. Information related to human diseases and functional annotation was also integrated into MetaImprint. To facilitate data extraction, MetaImprint supports multiple search options, such as by gene ID and disease name. Moreover, a configurable Imprinted Gene Browser was developed to visualize the information on imprinted genes in a genomic context. In addition, an Epigenetic Changes Analysis Tool is provided for online analysis of DNA methylation and histone modification differences of imprinted genes among multiple tissues and cell types. MetaImprint provides a comprehensive information repository of imprinted genes, allowing researchers to investigate systematically the genetic and epigenetic regulatory mechanisms of imprinted genes and their functions in development and diseases.