Epigenetic footprints: Investigating placental DNA methylation in the context of prenatal exposure to phenols and phthalates

BACKGROUND

Endocrine disrupting compounds (EDCs) such as phthalates and phenols can affect placental functioning and fetal health, potentially via epigenetic modifications. We investigated the associations between pregnancy exposure to synthetic phenols and phthalates estimated from repeated urine sampling and genome wide placental DNA methylation.

METHODS

The study is based on 387 women with placental DNA methylation assessed with Infinium MethylationEPIC arrays and with 7 phenols, 13 phthalates, and two non-phthalate plasticizer metabolites measured in pools of urine samples collected twice during pregnancy. We conducted an exploratory analysis on individual CpGs (EWAS) and differentially methylated regions (DMRs) as well as a candidate analysis focusing on 20 previously identified CpGs. Sex-stratified analyses were also performed.

RESULTS

In the exploratory analysis, when both sexes were studied together no association was observed in the EWAS. In the sex-stratified analysis, 114 individual CpGs (68 in males, 46 in females) were differentially methylated, encompassing 74 genes (36 for males and 38 for females). We additionally identified 28 DMRs in the entire cohort, 40 for females and 42 for males. Associations were mostly positive (for DMRs: 93% positive associations in the entire cohort, 60% in the sex-stratified analysis), with the exception of several associations for bisphenols and DINCH metabolites that were negative. Biomarkers associated with most DMRs were parabens, DEHP, and DiNP metabolite concentrations. Some DMRs encompassed imprinted genes including APC (associated with parabens and DiNP metabolites), GNAS (bisphenols), ZIM2;PEG3;MIMT1 (parabens, monoethyl phthalate), and SGCE;PEG10 (parabens, DINCH metabolites). Terms related to adiposity, lipid and glucose metabolism, and cardiovascular function were among the enriched phenotypes associated with differentially methylated CpGs. The candidate analysis identified one CpG mapping to imprinted LGALS8 gene, negatively associated with ethylparaben.

CONCLUSIONS

By combining improved exposure assessment and extensive placental epigenome coverage, we identified several novel genes associated with the exposure, possibly in a sex-specific manner.

Machine learning on alignment features for parent-of-origin classification of simulated hybrid RNA-seq

BACKGROUND

Parent-of-origin allele-specific gene expression (ASE) can be detected in interspecies hybrids by virtue of RNA sequence variants between the parental haplotypes. ASE is detectable by differential expression analysis (DEA) applied to the counts of RNA-seq read pairs aligned to parental references, but aligners do not always choose the correct parental reference.

RESULTS

We used public data for species that are known to hybridize. We measured our ability to assign RNA-seq read pairs to their proper transcriptome or genome references. We tested software packages that assign each read pair to a reference position and found that they often favored the incorrect species reference. To address this problem, we introduce a post process that extracts alignment features and trains a random forest classifier to choose the better alignment. On each simulated hybrid dataset tested, our machine-learning post-processor achieved higher accuracy than the aligner by itself at choosing the correct parent-of-origin per RNA-seq read pair.

CONCLUSIONS

For the parent-of-origin classification of RNA-seq, machine learning can improve the accuracy of alignment-based methods. This approach could be useful for enhancing ASE detection in interspecies hybrids, though RNA-seq from real hybrids may present challenges not captured by our simulations. We believe this is the first application of machine learning to this problem domain.

Prenatal exposure to triclosan assessed in multiple urine samples and placental DNA methylation

A previous study reported positive associations of maternal urinary concentrations of triclosan, a synthetic phenol with widespread exposure in the general population, with placental DNA methylation of male fetuses. Given the high number of comparisons performed in -omic research, further studies were needed to validate and extend on these findings. Using a cohort of male and female fetuses with repeated maternal urine samples to assess exposure, we studied the associations between triclosan and placental DNA methylation. We assessed triclosan concentrations in two pools of 21 urine samples collected among 395 women from the SEPAGES cohort. We used Infinium Methylation EPIC arrays to measure DNA methylation in placental biopsies collected at delivery. We performed a candidate study restricted to a set of candidate CpGs (n = 500) identified in a previous work as well as an exploratory epigenome-wide association study to investigate the associations between triclosan and differentially methylated probes and regions. Analyses were conducted on the whole population and stratified by child's sex. Mediation analysis was performed to test whether heterogeneity of placental tissue may mediate the observed associations. In the candidate approach, we confirmed 18 triclosan-associated genes when both sexes were considered. After stratification for child's sex, triclosan was associated with 72 genes in females and three in males. Most of the associations were positive and several CpGs mapped to imprinted genes: FBRSL1, KCNQ1, RHOBTB3, and SMOC1. A mediation effect by placental tissue heterogeneity was identified for most of the observed associations. In the exploratory analysis, we identified a few isolated associations in the sex-stratified analysis. In line with a previous study on male placentas, our approach revealed several positive associations between triclosan exposure and placental DNA methylation. Several identified loci mapped to imprinted genes.

Evolution of parent-of-origin effects on placental gene expression in house mice

The mammalian placenta is a hotspot for the evolution of genomic imprinting, a form of gene regulation that involves the parent-specific epigenetic silencing of one allele. Imprinted genes are central to placental development and are thought to contribute to the evolution of reproductive barriers between species. However, it is unclear how rapidly imprinting evolves or how functional specialization among placental tissues influences the evolution of imprinted expression. We compared parent-of-origin expression bias across functionally distinct placental layers sampled from reciprocal crosses within three closely related lineages of mice ( Mus ). Using genome-wide gene expression and DNA methylation data from fetal and maternal tissues, we developed an analytical strategy to minimize pervasive bias introduced by maternal contamination of placenta samples. We corroborated imprinted expression at 42 known imprinted genes and identified five candidate imprinted genes showing parent-of-origin specific expression and DNA methylation. Paternally-biased expression was enriched in the labyrinth zone, a layer specialized in nutrient transfer, and maternally-biased genes were enriched in the junctional zone, which specializes in modulation of maternal physiology. Differentially methylated regions were predominantly determined through epigenetic modification of the maternal genome and were associated with both maternally- and paternally-biased gene expression. Lastly, comparisons between lineages revealed a small set of co-regulated genes showing rapid divergence in expression levels and imprinted status in the M. m. domesticus lineage. Together, our results reveal important links between core functional elements of placental biology and the evolution of imprinted gene expression among closely related rodent species.

dbEmbryo multi-omics database for analyses of synergistic regulation in early mammalian embryo development

During early mammalian embryo development, different epigenetic marks undergo reprogramming and play crucial roles in the mediation of gene expression. Currently, several databases provide multi-omics information on early embryos. However, how interconnected epigenetic markers function together to coordinate the expression of the genetic code in a spatiotemporal manner remains difficult to analyze, markedly limiting scientific and clinical research. Here, we present dbEmbryo, an integrated and interactive multi-omics database for human and mouse early embryos. dbEmbryo integrates data on gene expression, DNA methylation, histone modifications, chromatin accessibility, and higher-order chromatin structure profiles for human and mouse early embryos. It incorporates customized analysis tools, such as "multi-omics visualization," "Gene&Peak annotation," "ZGA gene cluster," "cis-regulation," "synergistic regulation," "promoter signal enrichment," and "3D genome." Users can retrieve gene expression and epigenetic profile patterns to analyze synergistic changes across different early embryo developmental stages. We showed the uniqueness of dbEmbryo among extant databases containing data on early embryo development and provided an overview. Using dbEmbryo, we obtained a phase-separated model of transcriptional control during early embryo development. dbEmbryo offers web-based analytical tools and a comprehensive resource for biologists and clinicians to decipher molecular regulatory mechanisms of human and mouse early embryo development.

Allele-specific H3K9me3 and DNA methylation co-marked CpG-rich regions serve as potential imprinting control regions in pre-implantation embryo

Parental DNA methylation and histone modifications undergo distinct global reprogramming in mammalian pre-implantation embryos, but the landscape of epigenetic crosstalk and its effects on embryogenesis are largely unknown. Here we comprehensively analyse the association between DNA methylation and H3K9me3 reprogramming in mouse pre-implantation embryos and reveal that CpG-rich genomic loci with high H3K9me3 signal and DNA methylation level (CHM) are hotspots of DNA methylation maintenance during pre-implantation embryogenesis. We further profile the allele-specific epigenetic map with unprecedented resolution in gynogenetic and androgenetic embryos, respectively, and identify 1,279 allele-specific CHMs, including 19 known imprinting control regions (ICRs). Our study suggests that 22 ICR-like regions (ICRLRs) may regulate allele-specific transcription similarly to known ICRs, and five of them are confirmed to be important for mouse embryo development. Taken together, our study reveals the widespread existence of allele-specific CHMs and largely extends the scope of allele-specific regulation in mammalian pre-implantation embryos.

Dynamic nucleosome organization after fertilization reveals regulatory factors for mouse zygotic genome activation

Chromatin remodeling is essential for epigenome reprogramming after fertilization. However, the underlying mechanisms of chromatin remodeling remain to be explored. Here, we investigated the dynamic changes in nucleosome occupancy and positioning in pronucleus-stage zygotes using ultra low-input MNase-seq. We observed distinct features of inheritance and reconstruction of nucleosome positioning in both paternal and maternal genomes. Genome-wide de novo nucleosome occupancy in the paternal genome was observed as early as 1 h after the injection of sperm into ooplasm. The nucleosome positioning pattern was continually rebuilt to form nucleosome-depleted regions (NDRs) at promoters and transcription factor (TF) binding sites with differential dynamics in paternal and maternal genomes. NDRs formed more quickly on the promoters of genes involved in zygotic genome activation (ZGA), and this formation is closely linked to histone acetylation, but not transcription elongation or DNA replication. Importantly, we found that NDR establishment on the binding motifs of specific TFs might be associated with their potential pioneer functions in ZGA. Further investigations suggested that the predicted factors MLX and RFX1 played important roles in regulating minor and major ZGA, respectively. Our data not only elucidate the nucleosome positioning dynamics in both male and female pronuclei following fertilization, but also provide an efficient method for identifying key transcription regulators during development.

The imprinted Igf2-Igf2r axis is critical for matching placental microvasculature expansion to fetal growth

In all eutherian mammals, growth of the fetus is dependent upon a functional placenta, but whether and how the latter adapts to putative fetal signals is currently unknown. Here, we demonstrate, through fetal, endothelial, hematopoietic, and trophoblast-specific genetic manipulations in the mouse, that endothelial and fetus-derived IGF2 is required for the continuous expansion of the feto-placental microvasculature in late pregnancy. The angiocrine effects of IGF2 on placental microvasculature expansion are mediated, in part, through IGF2R and angiopoietin-Tie2/TEK signaling. Additionally, IGF2 exerts IGF2R-ERK1/2-dependent pro-proliferative and angiogenic effects on primary feto-placental endothelial cells ex vivo. Endothelial and fetus-derived IGF2 also plays an important role in trophoblast morphogenesis, acting through Gcm1 and Synb. Thus, our study reveals a direct role for the imprinted Igf2-Igf2r axis on matching placental development to fetal growth and establishes the principle that hormone-like signals from the fetus play important roles in controlling placental microvasculature and trophoblast morphogenesis.

Pregnancy exposure to synthetic phenols and placental DNA methylation - An epigenome-wide association study in male infants from the EDEN cohort

In utero exposure to environmental chemicals, such as synthetic phenols, may alter DNA methylation in different tissues, including placenta - a critical organ for fetal development. We studied associations between prenatal urinary biomarker concentrations of synthetic phenols and placental DNA methylation. Our study involved 202 mother-son pairs from the French EDEN cohort. Nine phenols were measured in spot urine samples collected between 22 and 29 gestational weeks. We performed DNA methylation analysis of the fetal side of placental tissues using the IlluminaHM450 BeadChips. We evaluated methylation changes of individual CpGs in an adjusted epigenome-wide association study (EWAS) and identified differentially methylated regions (DMRs). We performed mediation analysis to test whether placental tissue heterogeneity mediated the association between urinary phenol concentrations and DNA methylation. We identified 46 significant DMRs (≥5 CpGs) associated with triclosan (37 DMRs), 2,4-dichlorophenol (3), benzophenone-3 (3), methyl- (2) and propylparaben (1). All but 2 DMRs were positively associated with phenol concentrations. Out of the 46 identified DMRs, 7 (6 for triclosan) encompassed imprinted genes (APC, FOXG1, GNAS, GNASAS, MIR886, PEG10, SGCE), which represented a significant enrichment. Other identified DMRs encompassed genes encoding proteins responsible for cell signaling, transmembrane transport, cell adhesion, inflammatory, apoptotic and immunological response, genes encoding transcription factors, histones, tumor suppressors, genes involved in tumorigenesis and several cancer risk biomarkers. Mediation analysis suggested that placental cell heterogeneity may partly explain these associations. This is the first study describing the genome-wide modifications of placental DNA methylation associated with pregnancy exposure to synthetic phenols or their precursors. Our results suggest that cell heterogeneity might mediate the effects of triclosan exposure on placental DNA methylation. Additionally, the enrichment of imprinted genes within the DMRs suggests mechanisms by which certain exposures, mainly to triclosan, could affect fetal development.

Bioinformatic tools for DNA methylation and histone modification: A survey

Epigenetic inheritance occurs due to different mechanisms such as chromatin and histone modifications, DNA methylation and processes mediated by non-coding RNAs. It leads to changes in gene expressions and the emergence of new traits in different organisms in many diseases such as cancer. Recent advances in experimental methods led to the identification of epigenetic target sites in various organisms. Computational approaches have enabled us to analyze mass data produced by these methods. Next-generation sequencing (NGS) methods have been broadly used to identify these target sites and their patterns. By using these patterns, the emergence of diseases could be prognosticated. In this study, target site prediction tools for two major epigenetic mechanisms comprising histone modification and DNA methylation are reviewed. Publicly accessible databases are reviewed as well. Some suggestions regarding the state-of-the-art methods and databases have been made, including examining patterns of epigenetic changes that are important in epigenotypes detection.

Immediate and durable effects of maternal tobacco consumption alter placental DNA methylation in enhancer and imprinted gene-containing regions

BACKGROUND

Although exposure to cigarette smoking during pregnancy has been associated with alterations of DNA methylation in the cord blood or placental cells, whether such exposure before pregnancy could induce epigenetic alterations in the placenta of former smokers has never been investigated.

METHODS

Our approach combined the analysis of placenta epigenomic (ENCODE) data with newly generated DNA methylation data obtained from 568 pregnant women, the largest cohort to date, either actively smoking during their pregnancy or formerly exposed to tobacco smoking.

RESULTS

This strategy resulted in several major findings. First, among the 203 differentially methylated regions (DMRs) identified by the epigenome-wide association study, 152 showed "reversible" alterations of DNA methylation, only present in the placenta of current smokers, whereas 26 were also found altered in former smokers, whose placenta had not been exposed directly to cigarette smoking. Although the absolute methylation changes were smaller than those observed in other contexts, such as in some congenital diseases, the observed alterations were consistent within each DMR. This observation was further supported by a demethylation of LINE-1 sequences in the placentas of both current (beta-coefficient (β) (95% confidence interval (CI)), - 0.004 (- 0.008; 0.001)) and former smokers (β (95% CI), - 0.006 (- 0.011; - 0.001)) compared to nonsmokers. Second, the 203 DMRs were enriched in epigenetic marks corresponding to enhancer regions, including monomethylation of lysine 4 and acetylation of lysine 27 of histone H3 (respectively H3K4me1 and H3K27ac). Third, smoking-associated DMRs were also found near and/or overlapping 10 imprinted genes containing regions (corresponding to 16 genes), notably including the NNAT, SGCE/PEG10, and H19/MIR675 loci.

CONCLUSIONS

Our results pointing towards genomic regions containing the imprinted genes as well as enhancers as preferential targets suggest mechanisms by which tobacco could directly impact the fetus and future child. The persistence of significant DNA methylation changes in the placenta of former smokers supports the hypothesis of an "epigenetic memory" of exposure to cigarette smoking before pregnancy. This observation not only is conceptually revolutionary, but these results also bring crucial information in terms of public health concerning potential long-term detrimental effects of smoking in women.

Genomic Imprinting As a Window into Human Language Evolution

Humans spend large portions of their time and energy talking to one another, yet it remains unclear whether this activity is primarily selfish or altruistic. Here, it is shown how parent-of-origin specific gene expression-or "genomic imprinting"-may provide an answer to this question. First, it is shown why, regarding language, only altruistic or selfish scenarios are expected. Second, it is pointed out that an individual's maternal-origin and paternal-origin genes may have different evolutionary interests regarding investment into language, and that this intragenomic conflict may drive genomic imprinting which-as the direction of imprint depends upon whether investment into language is relatively selfish or altruistic-may be used to discriminate between these two possibilities. Third, predictions concerning the impact of various mutations and epimutations at imprinted loci on language pathologies are derived. In doing so, a framework is developed that highlights avenues for using intragenomic conflicts to investigate the evolutionary drivers of language.

Identification and comprehensive characterization of lncRNAs with copy number variations and their driving transcriptional perturbed subpathways reveal functional significance for cancer

Numerous studies have shown that copy number variation (CNV) in lncRNA regions play critical roles in the initiation and progression of cancer. However, our knowledge about their functionalities is still limited. Here, we firstly provided a computational method to identify lncRNAs with copy number variation (lncRNAs-CNV) and their driving transcriptional perturbed subpathways by integrating multidimensional omics data of cancer. The high reliability and accuracy of our method have been demonstrated. Then, the method was applied to 14 cancer types, and a comprehensive characterization and analysis was performed. LncRNAs-CNV had high specificity in cancers, and those with high CNV level may perturb broad biological functions. Some core subpathways and cancer hallmarks widely perturbed by lncRNAs-CNV were revealed. Moreover, subpathways highlighted the functional diversity of lncRNAs-CNV in various cancers. Survival analysis indicated that functional lncRNAs-CNV could be candidate prognostic biomarkers for clinical applications, such as ST7-AS1, CDKN2B-AS1 and EGFR-AS1. In addition, cascade responses and a functional crosstalk model among lncRNAs-CNV, impacted genes, driving subpathways and cancer hallmarks were proposed for understanding the driving mechanism of lncRNAs-CNV. Finally, we developed a user-friendly web interface-LncCASE (http://bio-bigdata.hrbmu.edu.cn/LncCASE/) for exploring lncRNAs-CNV and their driving subpathways in various cancer types. Our study identified and systematically characterized lncRNAs-CNV and their driving subpathways and presented valuable resources for investigating the functionalities of non-coding variations and the mechanisms of tumorigenesis.

Pervasive Inter-Individual Variation in Allele-Specific Expression in Monozygotic Twins

Despite being developed from one zygote, heterokaryotypic monozygotic (MZ) co-twins exhibit discordant karyotypes. Epigenomic studies in biological samples from heterokaryotypic MZ co-twins are of the most significant value for assessing the effects on gene- and allele-specific expression of an extranumerary chromosomal copy or structural chromosomal disparities in otherwise nearly identical germline genetic contributions. Here, we use RNA-Seq data from existing repositories to establish within-pair correlations for the breadth and magnitude of allele-specific expression (ASE) in heterokaryotypic MZ co-twins discordant for trisomy 21 and maternal 21q inheritance, as well as homokaryotypic co-twins. We show that there is a genome-wide disparity at ASE sites between the heterokaryotypic MZ co-twins. Although most of the disparity corresponds to changes in the magnitude of biallelic imbalance, ASE sites switching from either strictly monoallelic to biallelic imbalance or the reverse occur in few genes that are known or predicted to be imprinted, subject to X-chromosome inactivation or A-to-I(G) RNA edited. We also uncovered comparable ASE differences between homokaryotypic MZ twins. The extent of ASE discordance in MZ twins (2.7%) was about 10-fold lower than the expected between pairs of unrelated, non-twin males or females. The results indicate that the observed within-pair dissimilarities in breadth and magnitude of ASE sites in the heterokaryotypic MZ co-twins could not solely be attributable to the aneuploidy and the missing allelic heritability at 21q.

Development of a method for identifying and functionally analyzing allele-specific DNA methylation based on BS-seq data

Aim: To comprehensively identify allele-specific DNA methylation (ASM) at the genome-wide level. Methods: Here, we propose a new method, called GeneASM, to identify ASM using high-throughput bisulfite sequencing data in the absence of haplotype information. Results: A total of 2194 allele-specific DNA methylated genes were identified in the GM12878 lymphocyte lineage using GeneASM. These genes are mainly enriched in cell cytoplasm function, subcellular component movement or cellular linkages. GM12878 methylated DNA immunoprecipitation sequencing, and methylation sensitive restriction enzyme sequencing data were used to evaluate ASM. The relationship between ASM and disease was further analyzed using the The Cancer Genome Atlas (TCGA) data of lung adenocarcinoma (LUAD), and whole genome bisulfite sequencing data. Conclusion: GeneASM, which recognizes ASM by high-throughput bisulfite sequencing and heterozygous single-nucleotide polymorphisms, provides new perspective for studying genomic imprinting.

Characterization of bovine (Bos taurus) imprinted genes from genomic to amino acid attributes by data mining approaches

Genomic imprinting results in monoallelic expression of genes in mammals and flowering plants. Understanding the function of imprinted genes improves our knowledge of the regulatory processes in the genome. In this study, we have employed classification and clustering algorithms with attribute weighting to specify the unique attributes of both imprinted (monoallelic) and biallelic expressed genes. We have obtained characteristics of 22 known monoallelically expressed (imprinted) and 8 biallelic expressed genes that have been experimentally validated alongside 208 randomly selected genes in bovine (Bos taurus). Attribute weighting methods and various supervised and unsupervised algorithms in machine learning were applied. Unique characteristics were discovered and used to distinguish mono and biallelic expressed genes from each other in bovine. To obtain the accuracy of classification, 10-fold cross-validation with concerning each combination of attribute weighting (feature selection) and machine learning algorithms, was used. Our approach was able to accurately predict mono and biallelic genes using the genomics and proteomics attributes.

EmExplorer: a database for exploring time activation of gene expression in mammalian embryos

Understanding early development offers a striking opportunity to investigate genetic disease, stem cell and assisted reproductive technology. Recent advances in high-throughput sequencing technology have led to the rising influx of omics data, which have rapidly boosted our understanding of mammalian developmental mechanisms. Here, we review the database EmExplorer (a database for exploring time activation of gene expression in mammalian embryos), which systematically organizes the genes from development-related pathways, and which we have already established and continue to update it. The current version of EmExplorer incorporates over 26 000 genes obtained from 306 functional pathways in five species. The function annotations of development-related genes were also integrated into EmExplorer. To facilitate data extraction, the database also contains the following information. (i) The dynamic expression values for each development stage are matched to the corresponding genes. (ii) A two-layer search tool which supports multi-option searching, such as by official symbol, pathway name and function annotation. The returned entries can directly link to the analysis results for the corresponding gene or pathway in the analysis module. (iii) The analysis module provides different gene comparisons at the multi-species level and functional pathway level, which shows the species specificity and stage specificity at the gene or pathway level. (iv) The analysis based on the hypergeometric distribution test reveals the enrichment of gene functions at a particular stage of one organism's pathway. (v) The browser is designed for users with ambiguous searching goals and greatly helps new users to get a general idea of the contents of the database. (vi) The experimentally validated pathways are manually curated and shown on the home page. EmExplorer will be helpful for elucidating early developmental mechanisms and exploring time activation genes. EmExplorer is freely available at http://bioinfor.imu.edu.cn/emexplorer.

Assessment of maternal and parent of origin effects in genetic variation of economic traits in Iranian native fowl

1. The objective of the study was to investigate the influence of maternal and parent of origin effects (POE) on genetic variation of Iranian native fowl on economic traits. 2. Studied traits were body weights at birth (BW0), at eight (BW8) and 12 weeks of age (BW12), age (ASM) and weight at sexual maturity (WSM), egg number (EN) and average egg weight (AEW). 3. Several models, including additive, maternal additive genetics, permanent environmental effects and POE were compared using Wombat software. Bayesian Information Criterion (BIC) was used to identify the best model for each trait. The chance of reranking of birds between models was investigated using Spearman correlation and Wilcoxon rank test. 4. Based on the best model, direct heritability estimates for BW0, BW8, BW12, ASM, WSM, EN and AEW traits were 0.05, 0.21, 0.23, 0.30, 0.39, 0.22 and 0.38, respectively. Proportion of variance due to paternal POE for BW8 was 4% and proportion of variance due to maternal POE for BW12 was 5%. 5. Estimated maternal heritability for BW0 was 0.30 and for BW8 and BW12 were 0.00 and 0.01, respectively, which shows that maternal heritability was reduced by age. 6. Based on the results, considering POE for BW8 and BW12 and maternal genetic effects for BW0 improved the accuracy of estimations and avoid reranking of birds for these traits.

Tissue-Specific Monoallelic Expression of Bovine AXL is Associated with DNA Methylation of Promoter DMR

The AXL protein is a receptor tyrosine kinase and is often implicated in proliferation, migration and therapy resistance in various cancers. The AXL gene in humans is maternally expressed and paternally imprinted with differentially methylated regions (DMR) surrounding the promoter region. However, the imprinting status and epigenetic regulation of AXL gene in cattle remain unclear. Therefore, we explored the molecular structure along with the patterns of allelic expression and DNA methylation of the bovine AXL gene. First, the complete cDNA sequence of bovine AXL was gathered by Sanger method, from transcripts obtained from RT-PCR, 5' and 3' -RACE. In silico BLAST alignments showed that the longest mRNA sequence of bovine AXL consists of 19 exons and encodes a protein of 887 amino acids. We further analyzed the allelic expression of bovine AXL by employing single-nucleotide polymorphism (SNP)-based sequencing method. A SNP site (GenBank Accession no: rs210020651) found in exon 7 allowed us to distinguish the two parental alleles. Monoallelic expression of AXL was observed in four adult bovine tissues (heart, liver, spleen and fat), while biallelic expression was found in the other adult tissues such as the lung, kidney, muscle, brain and placenta. To determine whether the DNA methylation played a role in the tissue-specific imprinting of bovine AXL, we performed bisulfite sequencing of two regions: region 1 was a CpG island (CGI) in AXL promoter, mapping to 643 bp upstream of the transcription start site of AXL 5'-v1 transcripts, while region two was homologous to the region of human AXL DMR, with 10 CpG sites overlapping the first translation start site (TSS1) of bovine AXL. In region 2, DNA from both monoallelic and biallelic expressed tissues were mostly found to be completely unmethylated. However, tissue-specific differential methylation patterns were found in monoallelic expressed tissues such as the heart and liver while hypomethylation was noted in the promoter CpG island in biallelic expressed tissues such as the lung. These observations demonstrated that the tissue-specific monoallelic expression of bovine AXL is dependent on the DNA methylation of its promoter region.

Effects of maternal nutrition on the expression of genomic imprinted genes in ovine fetuses

Genomic imprinting is an epigenetic phenomenon of differential allelic expression based on parental origin. To date, 263 imprinted genes have been identified among all investigated mammalian species. However, only 21 have been described in sheep, of which 11 are annotated in the current ovine genome. Here, we aim to i) use DNA/RNA high throughput sequencing to identify new monoallelically expressed and imprinted genes in day 135 ovine fetuses and ii) determine whether maternal diet (100%, 60%, or 140% of National Research Council Total Digestible Nutrients) influences expression of imprinted genes. We also reported strategies to solve technical challenges in the data analysis pipeline. We identified 80 monoallelically expressed, 13 new putative imprinted genes, and five known imprinted genes in sheep using the 263 genes stated above as a guide. Sanger sequencing confirmed allelic expression of seven genes, CASD1, COPG2, DIRAS3, INPP5F, PLAGL1, PPP1R9A, and SLC22A18. Among the 13 putative imprinted genes, five were localized in the known sheep imprinting domains of MEST on chromosome 4, DLK1/GTL2 on chromosome 18 and KCNQ1 on chromosome 21, and three were in a novel sheep imprinted cluster on chromosome 4, known in other species as PEG10/SGCE. The expression of DIRAS3, IGF2, PHLDA2, and SLC22A18 was altered by maternal diet, albeit without allelic expression reversal. Together, our results expanded the list of sheep imprinted genes to 34 and demonstrated that while the expression levels of four imprinted genes were changed by maternal diet, the allelic expression patterns were un-changed for all imprinted genes studied.

Polymorphic Imprinting of SLC38A4 Gene in Bovine Placenta

Imprinted genes are characterized by monoallelic expression that is dependent on parental origin. Comparative analysis of imprinted genes between species is a powerful tool for understanding the biological significance of genomic imprinting. The slc38a4 gene encodes a neutral amino acid transporter and is identified as imprinted in mice. In this study, the imprinting status of SLC38A4 was assessed in bovine adult tissues and placenta using a polymorphism-based approach. Results indicate that SLC38A4 is not imprinted in eight adult bovine tissues including heart, liver, spleen, lung, kidney, muscle, fat, and brain. It was interesting to note that SLC38A4 showed polymorphic status in five heterogeneous placentas, with three exhibiting paternal monoallelic expression and two exhibiting biallelic expression. Monoallelic expression of imprinted genes is generally associated with allele-specific differentially methylation regions (DMRs) of CpG islands (CGIs)-encompassed promoter; therefore, the DNA methylation statuses of three CGIs in the SLC38A4 promoter and exon 1 region were tested in three placentas (two exhibiting paternal monoallelic and one showing biallelic expression of SLC38A4) and their corresponding paternal sperms. Unexpectedly, extreme hypomethylation (< 3%) of the DNA was observed in all the three detected placentas and their corresponding paternal sperms. The absence of DMR in bovine SLC38A4 promoter region implied that DNA methylation of these three CGIs does not directly or indirectly affect the polymorphic imprinting of SLC38A4 in bovine placenta. This suggested other epigenetic features other than DNA methylation are needed in regulating the imprinting of bovine SLC38A4, which is different from that of mouse with respect to a DMR existence at the mouse's slc38a4 promoter region. Although further work is needed, this first characterization of polymorphic imprinting status of SLC38A4 in cattle placenta provides valuable information on investigating the genomic imprinting phenomenon itself.

The tempo of human childhood: a maternal foot on the accelerator, a paternal foot on the brake

Relative to the life history of other great apes, that of humans is characterized by early weaning and short interbirth intervals (IBIs). We propose that in modern humans, birth until adrenarche, or the rise in adrenal androgens, developmentally corresponds to the period from birth until weaning in great apes and ancestral hominins. According to this hypothesis, humans achieved short IBIs by subdividing ancestral infancy into a nurseling phase, during which offspring fed at the breast, and a weanling phase, during which offspring fed specially prepared foods. Imprinted genes influence the timing of human weaning and adrenarche, with paternally expressed genes promoting delays in childhood maturation and maternally expressed genes promoting accelerated maturation. These observations suggest that the tempo of human development has been shaped by consequences for the fitness of kin, with faster development increasing maternal fitness at a cost to child fitness. The effects of imprinted genes suggest that the duration of the juvenile period (adrenarche until puberty) has also been shaped by evolutionary conflicts within the family.

Maternal 5mCpG Imprints at the PARD6G-AS1 and GCSAML Differentially Methylated Regions Are Decoupled From Parent-of-Origin Expression Effects in Multiple Human Tissues

A hallmark of imprinted genes in mammals is the occurrence of parent-of-origin-dependent asymmetry of DNA cytosine methylation (5^mC) of alleles at CpG islands (CGIs) in their promoter regions. This 5^mCpG asymmetry between the parental alleles creates allele-specific imprinted differentially methylated regions (iDMRs). iDMRs are often coupled to the transcriptional repression of the methylated allele and the activation of the unmethylated allele in a tissue-specific, developmental-stage-specific and/or isoform-specific fashion. iDMRs function as regulatory platforms, built through the recruitment of chemical modifications to histones to achieve differential, parent-of-origin-dependent chromatin segmentation states. Here, we used a comparative computational data mining approach to identify 125 novel constitutive candidate iDMRs that integrate the maximal number of allele-specific methylation region records overlapping CGIs in human methylomes. Twenty-nine candidate iDMRs display gametic 5^mCpG asymmetry, and another 96 are candidate secondary iDMRs. We established the maternal origin of the 5^mCpG imprints of one gametic (PARD6G-AS1) and one secondary (GCSAML) iDMRs. We also found a constitutively hemimethylated, nonimprinted domain at the PWWP2AP1 promoter CGI with oocyte-derived methylation asymmetry. Given that the 5^mCpG level at the iDMRs is not a sufficient criterion to predict active or silent locus states and that iDMRs can regulate genes from a distance of more than 1 Mb, we used RNA-Seq experiments from the Genotype-Tissue Expression project and public archives to assess the transcriptional expression profiles of SNPs across 4.6 Mb spans around the novel maternal iDMRs. We showed that PARD6G-AS1 and GCSAML are expressed biallelically in multiple tissues. We found evidence of tissue-specific monoallelic expression of ZNF124 and OR2L13, located 363 kb upstream and 419 kb downstream, respectively, of the GCSAML iDMR. We hypothesize that the GCSAML iDMR regulates the tissue-specific, monoallelic expression of ZNF124 but not of OR2L13. We annotated the non-coding epigenomic marks in the two maternal iDMRs using data from the Roadmap Epigenomics project and showed that the PARD6G-AS1 and GCSAML iDMRs achieve contrasting activation and repression chromatin segmentations. Lastly, we found that the maternal 5^mCpG imprints are perturbed in several hematopoietic cancers. We conclude that the maternal 5^mCpG imprints at PARD6G-AS1 and GCSAML iDMRs are decoupled from parent-of-origin transcriptional expression effects in multiple tissues.

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.

DNA methylation dynamics: identification and functional annotation

DNA methylation is an epigenetic modification of cytosines that undergoes dynamic changes in a temporal, spatial and cell-type-specific manner. Recent advances in technology have permitted the profiling of high-throughput methylomes in large numbers of biological samples. Various computational tools have been developed to identify and analyze DNA methylation dynamics in a variety of critical biological processes. As DNA methylation is becoming increasingly viewed as a dynamic process, the mechanisms governing DNA methylation dynamics and its roles in the transcriptional regulatory network are of great interest. It has been reported that DNA methylation dynamics plays essential roles in multiple biological processes, including development and cancer. As a functional event, the dynamics of DNA methylation have become increasingly relevant to many researchers. Here, we review state-of-the-art advances at three levels (genome-wide identification, regulatory mechanism investigation and the functional annotation) in the field of DNA methylation dynamics, as well as the future perspective of DNA methylation dynamics.

Systematic identification and annotation of human methylation marks based on bisulfite sequencing methylomes reveals distinct roles of cell type-specific hypomethylation in the regulation of cell identity genes

DNA methylation is a key epigenetic mark that is critical for gene regulation in multicellular eukaryotes. Although various human cell types may have the same genome, these cells have different methylomes. The systematic identification and characterization of methylation marks across cell types are crucial to understand the complex regulatory network for cell fate determination. In this study, we proposed an entropy-based framework termed SMART to integrate the whole genome bisulfite sequencing methylomes across 42 human tissues/cells and identified 757 887 genome segments. Nearly 75% of the segments showed uniform methylation across all cell types. From the remaining 25% of the segments, we identified cell type-specific hypo/hypermethylation marks that were specifically hypo/hypermethylated in a minority of cell types using a statistical approach and presented an atlas of the human methylation marks. Further analysis revealed that the cell type-specific hypomethylation marks were enriched through H3K27ac and transcription factor binding sites in cell type-specific manner. In particular, we observed that the cell type-specific hypomethylation marks are associated with the cell type-specific super-enhancers that drive the expression of cell identity genes. This framework provides a complementary, functional annotation of the human genome and helps to elucidate the critical features and functions of cell type-specific hypomethylation.

Genomic imprinting in psoriasis and atopic dermatitis: A review

Genomic imprinting is a genetic process where only one allele of a particular gene is expressed in a parent-of-origin dependent manner. Epigenetic changes in the DNA, such as methylation or acetylation of histones, are primarily thought to be responsible for silencing of the imprinted allele. Recently, global CpG methylation changes have been identified in psoriatic skin in comparison to normal skin, particularly near genes known to be upregulated in psoriasis such as KYNU, OAS2, and SERPINB3. Furthermore, imprinting has been associated with multi-chromosomal human disease, including diabetes and multiple sclerosis. This paper is the first to review the clinical and genetic evidence that exists in the literature for the association between imprinting and general skin disorders, including atopic dermatitis and psoriatic disease. Atopy was found to have evidence of imprinting on chromosomes 6, 11, 14, and 13. The β subunit of the IgE receptor on chromosome 11q12-13 may be imprinted. Psoriatic disease may be related to imprinting effects on chromosome 6 for psoriasis and 16 for psoriatic arthritis.

Genomic imprinting effects on complex traits in domesticated animal species

Monoallelically expressed genes that exert their phenotypic effect in a parent-of-origin specific manner are considered to be subject to genomic imprinting, the most well understood form of epigenetic regulation of gene expression in mammals. The observed differences in allele specific gene expression for imprinted genes are not attributable to differences in DNA sequence information, but to specific chemical modifications of DNA and chromatin proteins. Since the discovery of genomic imprinting some three decades ago, over 100 imprinted mammalian genes have been identified and considerable advances have been made in uncovering the molecular mechanisms regulating imprinted gene expression. While most genomic imprinting studies have focused on mouse models and human biomedical disorders, recent work has highlighted the contributions of imprinted genes to complex trait variation in domestic livestock species. Consequently, greater understanding of genomic imprinting and its effect on agriculturally important traits is predicted to have major implications for the future of animal breeding and husbandry. In this review, we discuss genomic imprinting in mammals with particular emphasis on domestic livestock species and consider how this information can be used in animal breeding research and genetic improvement programs.

Identification and characterization of long intergenic non-coding RNAs related to mouse liver development

Long non-coding RNAs (lncRNAs) have been studied extensively over the last few years. Liver is an important organ that plays a crucial role in glucose metabolism and homeostasis; however, there are few reports of the identification and functional characterization of lncRNAs with important roles in liver development. Therefore, it is necessary to systematically identify lncRNAs that are involved in liver development. In this paper, we assembled the transcriptome using published RNA-seq data across three mouse liver developmental stages and identified 4,882 putative long intergenic non-coding RNAs (lincRNAs) expressed in at least one of the investigated stages. Combining these with Ensembl lincRNAs, we established a reference catalog of 6,602 transcribed lincRNAs in the mouse liver. We then analyzed all the lincRNAs in this reference catalog systematically and revealed that liver lincRNAs carry different genomic signatures from protein-coding genes, while the putative lincRNAs are generally comparable with known Ensembl lincRNAs. In addition, putative lincRNAs are functionally associated with essential biological processes, including RNA splicing, protein localization and fatty acid metabolic process, implying that they may play an important role in regulating liver development. The validation of selected lincRNAs that are specifically expressed in developing liver tissues further suggested the effectiveness of our approach. Our study shows that lincRNAs that are differentially expressed during three liver developmental stages could have important regulatory roles in liver development. The identified putative lincRNAs are a valuable resource for further functional studies.