Evolutionary dynamics of microRNA target sites across vertebrate evolution

The rapidly evolving X-linked MIR-506 family fine-tunes spermatogenesis to enhance sperm competition

Despite rapid evolution across eutherian mammals, the X-linked MIR-506 family miRNAs are located in a region flanked by two highly conserved protein-coding genes (SLITRK2 and FMR1) on the X chromosome. Intriguingly, these miRNAs are predominantly expressed in the testis, suggesting a potential role in spermatogenesis and male fertility. Here, we report that the X-linked MIR-506 family miRNAs were derived from the MER91C DNA transposons. Selective inactivation of individual miRNAs or clusters caused no discernible defects, but simultaneous ablation of five clusters containing 19 members of the MIR-506 family led to reduced male fertility in mice. Despite normal sperm counts, motility, and morphology, the KO sperm were less competitive than wild-type sperm when subjected to a polyandrous mating scheme. Transcriptomic and bioinformatic analyses revealed that these X-linked MIR-506 family miRNAs, in addition to targeting a set of conserved genes, have more targets that are critical for spermatogenesis and embryonic development during evolution. Our data suggest that the MIR-506 family miRNAs function to enhance sperm competitiveness and reproductive fitness of the male by finetuning gene expression during spermatogenesis.

Convergent and divergent evolution of microRNA-mediated regulation in metazoans

The evolution of microRNAs (miRNAs) has been studied extensively to understand their roles in gene regulation and evolutionary processes. This review focuses on how miRNA-mediated regulation has evolved in bilaterian animals, highlighting both convergent and divergent evolution. Since animals and plants display significant differences in miRNA biogenesis and target recognition, the 'independent origin' hypothesis proposes that miRNA pathways in these groups independently evolved from the RNA interference (RNAi) pathway, leading to modern miRNA repertoires through convergent evolution. However, recent evidence raises the alternative possibility that the miRNA pathway might have already existed in the last common ancestor of eukaryotes, and that the differences in miRNA pathway and miRNA repertoires among animal and plant lineages arise from lineage-specific innovations and losses of miRNA pathways, miRNA acquisition, and loss of miRNAs after eukaryotic divergence. The repertoire of miRNAs has considerably expanded during bilaterian evolution, primarily through de novo creation and duplication processes, generating new miRNAs. Although ancient functionally established miRNAs are rarely lost, many newly emerged miRNAs are transient and lineage specific, following a birth-death evolutionary pattern aligning with the 'out-of-the-testis' and 'transcriptional control' hypotheses. Our focus then shifts to the convergent molecular evolution of miRNAs. We summarize how miRNA clustering and seed mimicry contribute to this phenomenon, and we review how miRNAs from different sources converge to degrade maternal messenger RNAs (mRNAs) during animal development. Additionally, we describe how miRNAs evolve across species due to changes in sequence, seed shifting, arm switching, and spatiotemporal expression patterns, which can result in variations in target sites among orthologous miRNAs across distant strains or species. We also provide a summary of the current understanding regarding how the target sites of orthologous miRNAs can vary across strains or distantly related species. Although many paralogous miRNAs retain their seed or mature sequences after duplication, alterations can occur in the seed or mature sequences or expression patterns of paralogous miRNAs, leading to functional diversification. We discuss our current understanding of the functional divergence between duplicated miRNAs, and illustrate how the functional diversification of duplicated miRNAs impacts target site evolution. By investigating these topics, we aim to enhance our current understanding of the functions and evolutionary dynamics of miRNAs. Additionally, we shed light on the existing challenges in miRNA evolutionary studies, particularly the complexity of deciphering the role of miRNA-mediated regulatory network evolution in shaping gene expression divergence and phenotypic differences among species.

Rapid adaptation of cellular metabolic rate to the MicroRNA complements of mammals and its relevance to the evolution of endothermy

The metabolic efficiency of mammalian cells depends on the attenuation of intrinsic translation noise by microRNAs. We devised a metric of cellular metabolic rate (cMR), rMR/Mexp optimally fit to the number of microRNA families (mirFam), that is robust to variation in mass and sensitive to body temperature (Tb), consistent with the heat dissipation limit theory of Speakman and Król (2010). Using mirFam as predictor, an Ornstein-Uhlenbeck process of stabilizing selection, with an adaptive shift at the divergence of Boreoeutheria, accounted for 95% of the variation in cMR across mammals. Branchwise rates of evolution of cMR, mirFam and Tb concurrently increased 6- to 7-fold at the divergence of Boreoeutheria, independent of mass. Cellular MR variation across placental mammals was also predicted by the sum of model conserved microRNA-target interactions, revealing an unexpected degree of integration of the microRNA-target apparatus into the energy economy of the mammalian cell.

The Rapidly Evolving X-linked miR-506 Family Finetunes Spermatogenesis to Enhance Sperm Competition

Despite rapid evolution across eutherian mammals, the X-linked miR-506 family miRNAs are located in a region flanked by two highly conserved protein-coding genes (Slitrk2 and Fmr1) on the X chromosome. Intriguingly, these miRNAs are predominantly expressed in the testis, suggesting a potential role in spermatogenesis and male fertility. Here, we report that the X-linked miR-506 family miRNAs were derived from the MER91C DNA transposons. Selective inactivation of individual miRNAs or clusters caused no discernable defects, but simultaneous ablation of five clusters containing nineteen members of the miR-506 family led to reduced male fertility in mice. Despite normal sperm counts, motility and morphology, the KO sperm were less competitive than wild-type sperm when subjected to a polyandrous mating scheme. Transcriptomic and bioinformatic analyses revealed that these X-linked miR-506 family miRNAs, in addition to targeting a set of conserved genes, have more targets that are critical for spermatogenesis and embryonic development during evolution. Our data suggest that the miR-506 family miRNAs function to enhance sperm competitiveness and reproductive fitness of the male by finetuning gene expression during spermatogenesis.

The miR-181 family: Wide-ranging pathophysiological effects on cell fate and function

MicroRNAs (miRNAs) are epigenetic regulators that can target and inhibit translation of multiple mRNAs within a given cell type. As such, a number of different pathways and networks may be modulated as a result. In fact, miRNAs are known to regulate many cellular processes including differentiation, proliferation, inflammation, and metabolism. This review focuses on the miR-181 family and provides information from the published literature on the role of miR-181 homologs in regulating a range of activities in different cell types and tissues. Of note, we have not included details on miR-181 expression and function in the context of cancer since this is a broad topic area requiring independent review. Instead, we have focused on describing the function and mechanism of miR-181 family members on differentiation toward a number of cell lineages in various non-neoplastic conditions (e.g., immune/hematopoietic cells, osteoblasts, osteoclasts, chondrocytes, adipocytes). We have also provided information on how modulation of miR-181 homologs can have positive effects on disease states such as cardiac abnormalities, pulmonary arterial hypertension, thrombosis, osteoarthritis, and vascular inflammation. In this context, we have used some examples of FDA-approved drugs that modulate miR-181 expression. We conclude by discussing some common mechanisms by which miR-181 homologs appear to regulate a number of different cellular processes and how targeting specific miR-181 family members may lead to attractive therapeutic approaches to treat a number of human disease or repair conditions, including those associated with the aging process.

Repetitive DNA is Functional and Encodes Parts of the Non-Coding RNA Repertoire

This is a commentary on the article by Eviatar Nevo and Kexin Li entitled "Sympatric Speciation in Mole Rats and Wild Barley and Their Genome Repeatome Evolution: A Commentary", published recently in Advanced Genetics.

Evolution of miRNA binding sites and regulatory networks in cichlids

The divergence of regulatory regions and gene regulatory network (GRN) rewiring is a key driver of cichlid phenotypic diversity. However, the contribution of miRNA binding site turnover has yet to be linked to GRN evolution across cichlids. Here, we extend our previous studies by analysing the selective constraints driving evolution of miRNA and transcription factor (TF) binding sites of target genes, to infer instances of cichlid GRN rewiring associated with regulatory binding site turnover. Comparative analyses identified increased species-specific networks that are functionally associated to traits of cichlid phenotypic diversity. The evolutionary rewiring is associated with differential models of miRNA and TF binding site turnover, driven by a high proportion of fast-evolving polymorphic sites in adaptive trait genes compared to subsets of random genes. Positive selection acting upon discrete mutations in these regulatory regions is likely to be an important mechanism in rewiring GRNs in rapidly radiating cichlids. Regulatory variants of functionally associated miRNA and TF binding sites of visual opsin genes differentially segregate according to phylogeny and ecology of Lake Malawi species, identifying both rewired e.g. clade-specific and conserved network motifs of adaptive trait associated GRNs. Our approach revealed several novel candidate regulators, regulatory regions and three-node motifs across cichlid genomes with previously reported associations to known adaptive evolutionary traits.

Association of miRNA targetome variants in LAMC1 and GNB3 genes with colorectal cancer and obesity

Background

Colorectal cancer (CRC) is one of the most common obesity‐associated cancers. Inflammation is also considered the most important factor between obesity and CRC. This study aimed to examine miRNAs binding sites variants on inflammatory genes identified using bioinformatics and systematic approach on clinical samples that were collected from CRC patients and controls.

Methods

The candidate variants related to CRC inflammatory genes were obtained from genome‐wide association studies and their population‐specific haplotypes. The variants were analyzed according to their genomic position on the miRNA targetome. Targetome variants in inflammation‐related genes were selected for genetic association study by TaqMan genotyping assay.

Results

The GG genotype of rs7473 decreased the risk of obesity (p < 0.05). Heterozygous genotype (GA) of rs1547715 decreased the risk of CRC (p < 0.05). In the rs7473/rs1547715 genotype and haplotype, the frequencies of AA/GA and GG/AA lessened in CRC and obesity, respectively (p < 0.05).

Conclusions

The variants of rs7473 and rs1547715 were associated with obesity and CRC, respectively. The above‐mentioned associations could be made based on the interactions of these variants with miRNAs.

Novel biomarkers of microRNAs in gastric cancer; an overview from diagnosis to treatment

The fourth frequent disease in the world and the second cause of cancer-related death is gastric cancer (GC). In this way, over 80% of diagnoses are made in the middle to advanced degrees of the disease, underscoring the requirement for innovative biomarkers that can be identified quickly. Meaningly, biomarkers that can complement endoscopic diagnosis and be used to detect patients with a high risk of GC are desperately needed. These biomarkers will allow for the accurate prediction of therapy response and prognosis in GC patients, as well as the development of an optimal treatment strategy for each individual. Conspicoiusly, microRNAs (miRNAs) and small noncoding RNA regulates the expression of target mRNA and thereby modifies critical biological mechanisms. According to the data, abnormally miRNAs expression in GC is linked to tumor growth, carcinogenesis, aggression and distant metastasis. Importantly, miRNA expression patterns and next-generation sequencing (NGS) can also be applied to analyze kinds of tissues and cancers. Given the high death rates and poor prognosis of GC, and the absence of a clinical diagnostic factor that is adequately sensitive to GC, research into novel sensitive and specific markers for GC diagnosis is critical. In this review,we evaluate the latest research findings that suggest the feasibility and clinical utility of miRNAs in GC.

Binding sites of miRNA on the overexpressed genes of oral cancer using 7mer-seed match

The microRNAs having a length of ~ 19-22 nucleotides are the small, non-coding RNAs. The evolution of microRNAs in many disorders may hold the key to tackle complex challenges. Oral cancer belongs to the group of head and neck cancer. It occurs in the mouth region that appears as an ulcer. In this study, we collected information on the overexpressed genes of oral cancer. The coding sequences of the genes were derived from NCBI and the entire set of human microRNAs present in miRBASE 21 was retrieved. The human microRNAs that can target the overexpressed genes of oral cancer were determined with the aid of our in-house software. The interaction between microRNAs and the overexpressed genes was evaluated with 7mer-m8 model of microRNA targeting. The genes DKK1 and APLN paired with only one miRNA i.e., miR-447 and miR-6087, respectively. But the genes INHBA and MMP1 were found to be targeted by 2 miRNAs, while the genes FN1, FAP, TGFPI, COL4A1, COL4A2, and LOXL2 were found to be targeted by 16, 5, 9, 18, 29, and 11 miRNAs. Subsequently, several measures such as free energy, translation efficiency, and cosine similarity metric were used to estimate the binding process. It was found that the target region's stability was higher than the upstream and downstream zones. The overexpressed genes' GC contents were calculated, revealing that the codons in target miRNA region were overall GC rich as well as GC3 rich. Lastly, gene ontology was performed to better understand each gene's involvement in biological processes, molecular function, and cellular component. Our study showed the role of microRNAs in gene repression, which could possibly aid in the prognosis and diagnosis of oral cancer.

3' untranslated regions of Marburg and Ebola virus mRNAs possess negative regulators of translation that are modulated by ADAR1 editing

The filovirus family includes deadly pathogens such as Ebola virus (EBOV) and Marburg virus (MARV). A substantial portion of filovirus genomes encode 5′ and 3′ untranslated regions (UTRs) of viral mRNAs. Select viral genomic RNA sequences corresponding to 3′ UTRs are prone to editing by adenosine deaminase acting on RNA 1 (ADAR1). A reporter mRNA approach, in which different 5′ or 3′ UTRs were inserted into luciferase-encoding mRNAs, demonstrates that MARV 3′ UTRs yield different levels of reporter gene expression, suggesting modulation of translation. The modulation occurs in cells unable to produce microRNAs (miRNAs) and can be recapitulated in a MARV minigenome assay. Deletion mutants identified negative regulatory regions at the ends of the MARV nucleoprotein (NP) and large protein (L) 3′ UTRs. Apparent ADAR1 editing mutants were previously identified within the MARV NP 3′ UTR. Introduction of these changes into the MARV nucleoprotein (NP) 3′ UTR or deletion of the region targeted for editing enhances translation, as indicated by reporter assays and polysome analysis. In addition, the parental NP 3′ UTR, but not the edited or deletion mutant NP 3′ UTRs, induces a type I interferon (IFN) response upon transfection into cells. Because some EBOV isolates from the West Africa outbreak exhibited ADAR1 editing of the viral protein of 40 kDa (VP40) 3′ UTR, VP40 3′ UTRs with parental and edited sequences were similarly assayed. The EBOV VP40 3′ UTR edits also enhanced translation, but neither the wild-type nor the edited 3′ UTRs induced IFN. These findings implicate filoviral mRNA 3′ UTRs as negative regulators of translation that can be inactivated by innate immune responses that induce ADAR1.

IMPORTANCE UTRs comprise a large percentage of filovirus genomes and are apparent targets of editing by ADAR1, an enzyme with pro- and antiviral activities. However, the functional significance of the UTRs and ADAR1 editing has been uncertain. This study demonstrates that MARV and EBOV 3′ UTRs can modulate translation, in some cases negatively. ADAR1 editing or deletion of select regions within the translation suppressing 3′ UTRs relieves the negative effects of the UTRs. These data indicate that filovirus 3′ UTRs contain translation regulatory elements that are modulated by activation of ADAR1, suggesting a complex interplay between filovirus gene expression and innate immunity.

miRNA: local guardians of presynaptic function in plasticity and disease

Environmental fitness is an essential component of animal survival. Fitness is achieved through responsive physiological plasticity of tissues across the entire body, and particularly in the nervous system. At the molecular level, neural plasticity is mediated via gene-environmental interactions whereby developmental cues and experience dependent input adapt neuronal function to ever changing demands. To this end, neuronal gene regulation must be coupled to changes in neural activity. Seminal discoveries of the 20th century demonstrated neural activity modifies gene expression through calcium-dependent gene transcription. Building on this model, recent work over the last two decades shows that mRNA products of transcriptional programming continue to be regulated in the neuron through the activity-dependent post-transcriptional action of microRNAs (miRNAs). miRNAs are special post-transcriptional regulators that can tune gene expression within the spatial and temporal requirements of synaptic compartments. This mode of gene regulation has proven to be essential for synaptic function and plasticity as miRNA loss of function is highly associated with neural disease. In this review we will discuss current perspective on the link between presynaptic plasticity and miRNA biogenesis in the neuron.

Direct molecular evidence for an ancient, conserved developmental toolkit controlling post-transcriptional gene regulation in land plants

In plants, miRNA production is orchestrated by a suite of proteins that control transcription of the pri-miRNA gene, post-transcriptional processing and nuclear export of the mature miRNA. Post-transcriptional processing of miRNAs is controlled by a pair of physically interacting proteins, hyponastic leaves 1 (HYL1) and Dicer-like 1 (DCL1). However, the evolutionary history and structural basis of the HYL1–DCL1 interaction is unknown. Here we use ancestral sequence reconstruction and functional characterization of ancestral HYL1 in vitro and in Arabidopsis thaliana to better understand the origin and evolution of the HYL1–DCL1 interaction and its impact on miRNA production and plant development. We found the ancestral plant HYL1 evolved high affinity for both double-stranded RNA (dsRNA) and its DCL1 partner before the divergence of mosses from seed plants (∼500 Ma), and these high-affinity interactions remained largely conserved throughout plant evolutionary history. Structural modeling and molecular binding experiments suggest that the second of two dsRNA-binding motifs (DSRMs) in HYL1 may interact tightly with the first of two C-terminal DCL1 DSRMs to mediate the HYL1–DCL1 physical interaction necessary for efficient miRNA production. Transgenic expression of the nearly 200 Ma-old ancestral flowering-plant HYL1 in A. thaliana was sufficient to rescue many key aspects of plant development disrupted by HYL1⁻ knockout and restored near-native miRNA production, suggesting that the functional partnership of HYL1–DCL1 originated very early in and was strongly conserved throughout the evolutionary history of terrestrial plants. Overall, our results are consistent with a model in which miRNA-based gene regulation evolved as part of a conserved plant “developmental toolkit.”

EBV miRNAs are potent effectors of tumor cell transcriptome remodeling in promoting immune escape

The Epstein Barr virus (EBV) contributes to the tumor phenotype through a limited set of primarily non-coding viral RNAs, including 31 mature miRNAs. Here we investigated the impact of EBV miRNAs on remodeling the tumor cell transcriptome. Strikingly, EBV miRNAs displayed exceptionally abundant expression in primary EBV-associated Burkitt’s Lymphomas (BLs) and Gastric Carcinomas (GCs). To investigate viral miRNA targeting, we used the high-resolution approach, CLASH in GC and BL cell models. Affinity constant calculations of targeting efficacies for CLASH hits showed that viral miRNAs bind their targets more effectively than their host counterparts, as did Kaposi’s sarcoma-associated herpesvirus (KSHV) and murine gammaherpesvirus 68 (MHV68) miRNAs. Using public BL and GC RNA-seq datasets, we found that high EBV miRNA targeting efficacies translates to enhanced reduction of target expression. Pathway analysis of high efficacy EBV miRNA targets showed enrichment for innate and adaptive immune responses. Inhibition of the immune response by EBV miRNAs was functionally validated in vivo through the finding of inverse correlations between EBV miRNAs and immune cell infiltration and T-cell diversity in BL and GC datasets. Together, this study demonstrates that EBV miRNAs are potent effectors of the tumor transcriptome that play a role in suppressing host immune response.

Burkitt’s Lymphoma and gastric cancer are both associated with EBV, a prolific DNA tumor virus that latently resides in nearly all human beings. Despite mostly restricting viral gene expression to noncoding RNAs, EBV has important influences on the fitness of infected tumor cells. Here, we show that the miRNA class of viral noncoding RNAs are a major viral contributor to remodeling the tumor cell regulatory machinery in patient tumor samples. First, an assessment of miRNA expression in clinical tumor samples showed that EBV miRNAs are expressed at unexpectedly high levels relative to cell miRNAs. Using a highly specific miRNA target identification approach, CLASH, we comprehensively identified both viral and cellular miRNA targets and the relative abundance of each miRNA-mRNA interaction. We also show that viral miRNAs bind to and alter the expression of their mRNA targets more effectively than their cellular miRNA counterparts. Pathway analysis of the most effectively targeted mRNAs revealed enrichment of immune signaling pathways and we show a corresponding inverse correlation between EBV miRNA expression and infiltrating immune cells in EBV positive primary tumors. Altogether, this study shows that EBV miRNAs are key regulators of the tumor cell phenotype and the immune cell microenvironment.

Variability in porcine microRNA genes and its association with mRNA expression and lipid phenotypes

BACKGROUND

Mature microRNAs (miRNAs) play an important role in repressing the expression of a wide range of mRNAs. The presence of polymorphic sites in miRNA genes and their corresponding 3'UTR binding sites can disrupt canonical conserved miRNA-mRNA pairings, and thus modify gene expression patterns. However, to date such polymorphic sites in miRNA genes and their association with gene expression phenotypes and complex traits are poorly characterized in pigs.

RESULTS

By analyzing whole-genome sequences from 120 pigs and wild boars from Europe and Asia, we identified 285 single nucleotide polymorphisms (SNPs) that map to miRNA loci, and 109,724 SNPs that are located in predicted 7mer-m8 miRNA binding sites within porcine 3'UTR. In porcine miRNA genes, SNP density is reduced compared with their flanking non-miRNA regions. By sequencing the genomes of five Duroc boars, we identified 12 miRNA SNPs that were subsequently genotyped in their offspring (N = 345, Lipgen population). Association analyses of miRNA SNPs with 38 lipid-related traits and hepatic and muscle microarray expression phenotypes recorded in the Lipgen population were performed. The most relevant detected association was between the genotype of the rs319154814 (G/A) SNP located in the apical loop of the ssc-miR-326 hairpin precursor and PPP1CC mRNA levels in the liver (q-value = 0.058). This result was subsequently confirmed by qPCR (P-value = 0.027). The rs319154814 (G/A) genotype was also associated with several fatty acid composition traits.

CONCLUSIONS

Our findings show a reduced variability of porcine miRNA genes, which is consistent with strong purifying selection, particularly in the seed region that plays a critical role in miRNA binding. Although it is generally assumed that SNPs mapping to the seed region are those with the most pronounced consequences on mRNA expression, we show that a SNP mapping to the apical region of ssc-miR-326 is significantly associated with hepatic mRNA levels of the PPP1CC gene, one of its predicted targets. Although experimental confirmation of such an interaction is reported in humans but not in pigs, this result highlights the need to further investigate the functional effects of miRNA polymorphisms that are located outside the seed region on gene expression in pigs.

EBV miRNAs are potent effectors of tumor cell transcriptome remodeling in promoting immune escape

The Epstein Barr virus (EBV) contributes to the tumor phenotype through a limited set of primarily non-coding viral RNAs, including 31 mature miRNAs. Here we investigated the impact of EBV miRNAs on remodeling the tumor cell transcriptome. Strikingly, EBV miRNAs displayed exceptionally abundant expression in primary EBV-associated Burkitt's Lymphomas (BLs) and Gastric Carcinomas (GCs). To investigate viral miRNA targeting, we used the high-resolution approach, CLASH in GC and BL cell models. Affinity constant calculations of targeting efficacies for CLASH hits showed that viral miRNAs bind their targets more effectively than their host counterparts, as did Kaposi's sarcoma-associated herpesvirus (KSHV) and murine gammaherpesvirus 68 (MHV68) miRNAs. Using public BL and GC RNA-seq datasets, we found that high EBV miRNA targeting efficacies translates to enhanced reduction of target expression. Pathway analysis of high efficacy EBV miRNA targets showed enrichment for innate and adaptive immune responses. Inhibition of the immune response by EBV miRNAs was functionally validated in vivo through the finding of inverse correlations between EBV miRNAs and immune cell infiltration and T-cell diversity in BL and GC datasets. Together, this study demonstrates that EBV miRNAs are potent effectors of the tumor transcriptome that play a role in suppressing host immune response.

Silencing lung cancer genes using miRNAs identified by 7mer-seed matching

Lung cancer (LC) is the main cause of cancer-associated deaths in both men and women globally with a very high mortality rate. The microRNAs (miRNAs) are a class of noncoding RNAs consisting of 18-25 nucleotides. They inhibit translation of protein through binding to complementary target mRNAs. The non-coding miRNAs are recognized as potent biomarkers for detection, development and treatment of malignancy. In this study, we screened a set of 12 genes over expressed in small cell lung cancer, non small cell lung cancer and the genes involved in both categories and their binding sites for human miRNAs as no work was reported yet. Screening of human miRNAs revealed that a few genes showed numerous miRNA binding sites. Free energy values of mRNA sequences revealed that they might acquire compact folded structure causing complexity for miRNAs to interact. GC content in the target site was relatively higher than that of their flanks. It was observed through analysis of cosine similarity metric and compAI parameters that the genes related to lung cancer were encoded with non optimal codons and thus might be translationally less efficient for producing polypeptides. Gene ontology analysis was carried out to understand the diverse functions of these 12 genes.

Comparative epigenetics in animal physiology: An emerging frontier

The unprecedented access to annotated genomes now facilitates the investigation of the molecular basis of epigenetic phenomena in phenotypically diverse animals. In this critical review, we describe the roles of molecular epigenetic mechanisms in regulating mitotically and meiotically stable spatiotemporal gene expression, phenomena that provide the molecular foundation for the intra-, inter-, and trans-generational emergence of physiological phenotypes. By focusing principally on emerging comparative epigenetic roles of DNA-level and transcriptome-level epigenetic mark dynamics in the emergence of phenotypes, we highlight the relationship between evolutionary conservation and innovation of specific epigenetic pathways, and their interplay as a priority for future study. This comparative approach is expected to significantly advance our understanding of epigenetic phenomena, as animals show a diverse array of strategies to epigenetically modify physiological responses. Additionally, we review recent technological advances in the field of molecular epigenetics (single-cell epigenomics and transcriptomics and editing of epigenetic marks) in order to (1) investigate environmental and endogenous factor dependent epigenetic mark dynamics in an integrative manner; (2) functionally test the contribution of specific epigenetic marks for animal phenotypes via genome and transcript-editing tools. Finally, we describe advantages and limitations of emerging animal models, which under the Krogh principle, may be particularly useful in the advancement of comparative epigenomics and its potential translational applications in animal science, ecotoxicology, ecophysiology, climate change science and wild-life conservation, as well as organismal health.