Tissue-specific Transcription Research Articles

Tissue-Specific Analysis of Alternative Splicing Events and Differential Isoform Expression in Large Yellow Croaker (Larimichthys crocea) After Cryptocaryon irritans Infection.

The large yellow croaker (Larimichthys crocea) is one of the most important mariculture fish in China. Recently, cryptocaryonosis caused by Cryptocryon irritans infection has brought huge economic losses and threatens the healthy and sustainable development of the L. crocea industry. However, the molecular mechanism and regulation process for L. crocea resistance to C. irritans infection has not been fully researched. Alternative splicing (AS) is an important post-transcriptional regulatory mechanism that allows cells to produce transcriptional and proteomic diversity. The results of AS are tissue dependent, and the expression of tissue-specific transcription subtype genes is determined by AS and transcriptional regulation. However, studies on the tissue specificity of AS events in L. crocea following infection with C. irritans have not been performed. In this study, the L. crocea were artificially infected with C. irritans; their skin and gill were collected at 0h, 24h, 48h, 72h, and 96h post infection. After sequencing and differential expression analysis, a set of 452, 692, 934, 711, 534, and 297 differential alternative splicing(DAS) eventswere identified in 0h, 12h, 24h, 48h, 72h, and 96h post infection respectively. Furthermore, 4160differentially expressed isoforms (DEIs) and 4209 DEI genes were identified from all time point groups. GO enrichment and pathway analysis indicated that many genes of DAS and DEIs were rich in immune-related GO terms and KEGG pathways, such as the Toll and Imd signaling pathway, NOD-like receptor signaling pathway, TNF signaling pathway, and TNF signaling pathway. Among hub DEI genes, alternative splicing-related genes (cwc25, prpf8, and sf3a3), skin function-related gene (fa2h), and oxygen deprivation-related gene (hyo1) were found in DEI genes. This study provided insight into the temporal change of DAS and DEIs between skin and gill of L. crocea against C. irritans infection and revealed that these differences mightplay immune-related roles in the infection process.

Improving reusability along the data life cycle: a regulatory circuits case study

BackgroundIn life sciences, there has been a long-standing effort of standardization and integration of reference datasets and databases. Despite these efforts, many studies data are provided using specific and non-standard formats. This hampers the capacity to reuse the studies data in other pipelines, the capacity to reuse the pipelines results in other studies, and the capacity to enrich the data with additional information. The Regulatory Circuits project is one of the largest efforts for integrating human cell genomics data to predict tissue-specific transcription factor-genes interaction networks. In spite of its success, it exhibits the usual shortcomings limiting its update, its reuse (as a whole or partially), and its extension with new data samples. To address these limitations, the resource has previously been integrated in an RDF triplestore so that TF-gene interaction networks could be generated with two SPARQL queries. However, this triplestore did not store the computed networks and did not integrate metadata about tissues and samples, therefore limiting the reuse of this dataset. In particular, it does not enable to reuse only a portion of Regulatory Circuits if a study focuses on a subset of the tissues, nor to combine the samples described in the datasets with samples from other studies. Overall, these limitations advocate for the design of a complete, flexible and reusable representation of the Regulatory Circuits dataset based on Semantic Web technologies.ResultsWe provide a modular RDF representation of the Regulatory Circuits, called Linked Extended Regulatory Circuits (LERC). It consists in (i) descriptions of biological and experimental context mapped to the references databases, (ii) annotations about TF-gene interactions at the sample level for 808 samples, (iii) annotations about TF-gene interactions at the tissue level for 394 tissues, (iv) metadata connecting the knowledge graphs cited above. LERC is based on a modular organisation into 1,205 RDF named graphs for representing the biological data, the sample-specific and the tissue-specific networks, and the corresponding metadata. In total it contains 3,910,794,050 triples and is available as a SPARQL endpoint.ConclusionThe flexible and modular architecture of LERC supports biologically-relevant SPARQL queries. It allows an easy and fast querying of the resources related to the initial Regulatory Circuits datasets and facilitates its reuse in other studies.Associated websitehttps://regulatorycircuits-lod.genouest.org

Epigenomics analysis of miRNA cis-regulatory elements in pig muscle and fat tissues.

Although large-scale and accurate identification of cis-regulatory elements on pig protein-coding and long non-coding genes has been reported, similar study on pig miRNAs is still lacking. Here, we systematically characterized the cis-regulatory elements of pig miRNAs in muscle and fat by adopting miRNAomes, ChIP-seq, ATAC-seq, RNA-seq and Hi-C data. In total, the cis-regulatory elements of 257 (85.95%) expressed miRNAs including 226 known and 31 novel miRNAs were identified. Especially, the miRNAs associated with super-enhancers, active promoters, and "A" compartment were significantly higher than those associated by typical enhancers, prompters without H3K27ac, and "B" compartment, respectively. The tissue specific transcription factors were the primary determination of core miRNA expression pattern in muscle and fat. Moreover, the miRNA promoters are more evolutionarily conserved than miRNA enhancers, like other type genes. Our study adds additional important information to existing pig epigenetic data and provides essential resource for future in-depth investigation of pig epigenetics.

Targeted DamID in C. elegans reveals a direct role for LIN-22 and NHR-25 in antagonizing the epidermal stem cell fate.

Transcription factors are key players in gene networks controlling cell fate specification during development. In multicellular organisms, they display complex patterns of expression and binding to their targets, hence, tissue specificity is required in the characterization of transcription factor–target interactions. We introduce here targeted DamID (TaDa) as a method for tissue-specific transcription factor target identification in intact Caenorhabditis elegans animals. We use TaDa to recover targets in the epidermis for two factors, the HES1 homolog LIN-22, and the NR5A1/2 nuclear hormone receptor NHR-25. We demonstrate a direct link between LIN-22 and the Wnt signaling pathway through repression of the Frizzled receptor lin-17. We report a direct role for NHR-25 in promoting cell differentiation via repressing the expression of stem cell–promoting GATA factors. Our results expand our understanding of the epidermal gene network and highlight the potential of TaDa to dissect the architecture of tissue-specific gene regulatory networks.

Machine Learning Identifies Six Genetic Variants and Alterations in the Heart Atrial Appendage as Key Contributors to PD Risk Predictivity.

Parkinson’s disease (PD) is a complex neurodegenerative disease with a range of causes and clinical presentations. Over 76 genetic loci (comprising 90 SNPs) have been associated with PD by the most recent GWAS meta-analysis. Most of these PD-associated variants are located in non-coding regions of the genome and it is difficult to understand what they are doing and how they contribute to the aetiology of PD. We hypothesised that PD-associated genetic variants modulate disease risk through tissue-specific expression quantitative trait loci (eQTL) effects. We developed and validated a machine learning approach that integrated tissue-specific eQTL data on known PD-associated genetic variants with PD case and control genotypes from the Wellcome Trust Case Control Consortium. In so doing, our analysis ranked the tissue-specific transcription effects for PD-associated genetic variants and estimated their relative contributions to PD risk. We identified roles for SNPs that are connected with INPP5P, CNTN1, GBA and SNCA in PD. Ranking the variants and tissue-specific eQTL effects contributing most to the machine learning model suggested a key role in the risk of developing PD for two variants (rs7617877 and rs6808178) and eQTL associated transcriptional changes of EAF1-AS1 within the heart atrial appendage. Similarly, effects associated with eQTLs located within the Brain Cerebellum were also recognized to confer major PD risk. These findings were replicated in two additional, independent cohorts (the UK Biobank, and NeuroX) and thus warrant further mechanistic investigations to determine if these transcriptional changes could act as early contributors to PD risk and disease development.

Epigenetics of Mitochondria-Associated Genes in Striated Muscle.

Striated muscle has especially large energy demands. We identified 97 genes preferentially expressed in skeletal muscle and heart, but not in aorta, and found significant enrichment for mitochondrial associations among them. We compared the epigenomic and transcriptomic profiles of the 27 genes associated with striated muscle and mitochondria. Many showed strong correlations between their tissue-specific transcription levels, and their tissue-specific promoter, enhancer, or open chromatin as well as their DNA hypomethylation. Their striated muscle-specific enhancer chromatin was inside, upstream, or downstream of the gene, throughout much of the gene as a super-enhancer (CKMT2, SLC25A4, and ACO2), or even overlapping a neighboring gene (COX6A2, COX7A1, and COQ10A). Surprisingly, the 3′ end of the 1.38 Mb PRKN (PARK2) gene (involved in mitophagy and linked to juvenile Parkinson’s disease) displayed skeletal muscle/myoblast-specific enhancer chromatin, a myoblast-specific antisense RNA, as well as brain-specific enhancer chromatin. We also found novel tissue-specific RNAs in brain and embryonic stem cells within PPARGC1A (PGC-1α), which encodes a master transcriptional coregulator for mitochondrial formation and metabolism. The tissue specificity of this gene’s four alternative promoters, including a muscle-associated promoter, correlated with nearby enhancer chromatin and open chromatin. Our in-depth epigenetic examination of these genes revealed previously undescribed tissue-specific enhancer chromatin, intragenic promoters, regions of DNA hypomethylation, and intragenic noncoding RNAs that give new insights into transcription control for this medically important set of genes.

Infection with a hypervirulent strain of Helicobacter pylori primes gastric cells toward intestinal transdifferentiation

BackgroundIntestinal metaplasia, gastric-to-intestinal transdifferentiation, occurs as a result of the misexpression of certain regulatory factors, leading to genetic reprogramming. Here, we have evaluated the H. pylori-induced expression patterns of these candidate genes. MethodsThe expression levels of 1) tissue-specific transcription factors (RUNX3, KLF5, SOX2, SALL4, CDX1 and CDX2), 2) stemness factors (TNFRSF19, LGR5, VIL1) and 3) tissue-specific mucins (MUC5AC, MUC2) were evaluated by quantitative real-time PCR in gastric primary cells (GPCs), in parallel with two gastric cancer (MKN45 and AGS) cell lines, up to 96h following H. pylori infection. ResultsFollowing H. pylori infection of GPCs, RUNX3 declined at 24h post infection (−6.2 ± 0.3) and remained downregulated for up to 96h. Subsequently, overexpression of self-renewal and pluripotency transcription factors, KLF5 (3.6 ± 0.2), SOX2 (7.6 ± 0.5) and SALL4 (4.3 ± 0.2) occurred. The expression of TNFRSF19 and LGR5, demonstrated opposing trends, with an early rise of the former (4.5 ± 0.3) at 8h, and a simultaneous fall of the latter (-1.8 ± 0.5). This trend was reversed at 96h, with the decline in TNFRSF19 (−5.5 ± 0.2), and escalation of LGR5 (2.6 ± 0.2) and VIL1 (1.8 ± 0.3). Ultimately, CDX1 and CDX2 were upregulated by 1.9 and 4.7-fold, respectively. The above scenario was, variably observed in MKN45 and AGS cells. ConclusionOur data suggests an interdependent gene regulatory network, induced by H. pylori infection. This interaction begins with the downregulation of RUNX3, upregulation of self-renewal and pluripotency transcription factors, KLF5, SOX2 and SALL4, leading to the downregulation of TNFRSF19, upregulation of LGR5 and aberrant expression of intestine-specific transcription factors, potentially facilitating the process of gastric-to-intestinal transdifferentiation.

Dynamic Expression and Chromatin Incorporation of ACT and CREM Transcription Factors in Testis Tissues of Infertile Men.

ObjectiveActivator of CREM in the testis (ACT) is a tissue specific transcription factor which activates cAMP responsive element modulator (CREM), a key transcription factor in differentiation of round spermatids into mature spermatozoa. They bind to CRE region in the promoters of transition protein genes (TNP1, TNP2) and protamine genes (PRM1 and PRM2), which are essential for sperm chromatin compaction, and regulates their transcription. This study was conducted to consider the expression of ACT and CREM and their regulatory roles on the expression of PRM1, PRM2, TNP1 and TNP2 genes in testis tissues of infertile men.Materials and MethodsIn this case-control study, testicular biopsies were collected from 40 infertile men and classified into three groups: obstructive azoospermia (OA, n=10, positive control), round spermatid maturation arrest (SMA, n=20), Sertoli cell-only syndrome (SCOS, n=10, negative control group). Using quantitative real-time polymerase chain reaction (PCR), the expression profile of ACT, CREM, TNP1, TNP2, PRM1 and PRM2 genes were assessed in testicular samples and incorporation of ACT and CREM proteins on the promoters of PRM1, PRM2, TNP1 and TNP2 genes were also evaluated by ChIP-real time PCR. ResultsOur results demonstrated significant decrease in the expression levels of ACT, CREM and in their incorporations on their target genes in SMA group in comparison to control groups (P<0.05).ConclusionThese data confirm that there is low expression and incorporation of ACT and CREM and of their target genes in infertilities which are associated with post-meiotic arrest.

Satb2 acts as a gatekeeper for major developmental transitions during early vertebrate embryogenesis

Zygotic genome activation (ZGA) initiates regionalized transcription underlying distinct cellular identities. ZGA is dependent upon dynamic chromatin architecture sculpted by conserved DNA-binding proteins. However, the direct mechanistic link between the onset of ZGA and the tissue-specific transcription remains unclear. Here, we have addressed the involvement of chromatin organizer Satb2 in orchestrating both processes during zebrafish embryogenesis. Integrative analysis of transcriptome, genome-wide occupancy and chromatin accessibility reveals contrasting molecular activities of maternally deposited and zygotically synthesized Satb2. Maternal Satb2 prevents premature transcription of zygotic genes by influencing the interplay between the pluripotency factors. By contrast, zygotic Satb2 activates transcription of the same group of genes during neural crest development and organogenesis. Thus, our comparative analysis of maternal versus zygotic function of Satb2 underscores how these antithetical activities are temporally coordinated and functionally implemented highlighting the evolutionary implications of the biphasic and bimodal regulation of landmark developmental transitions by a single determinant.

Thyroid Transcription Factor-1: Structure, Expression, Function and Its Relationship with Disease.

Thyroid transcription factor-1 (TTF-1/NKx2.1) is a member of the NKx2 tissue-specific transcription factor family, which is expressed in thyroid follicle, parathyroid gland, alveolar epithelium, and diencephalon which originated from ectoderm, and participates in the differentiation, development, and functional maintenance of the above organs. Recent studies have shown that the abnormal expression of TTF-1 is closely related to the occurrence of a variety of human diseases and can be used as a potential new target for the diagnosis and treatment of related diseases. In this article, in order to strengthen the systematic understanding of TTF-1 and promote the progress of related research, we reviewed the structure, expression regulation, biological functions of TTF-1, and its role in the occurrence and development of human-related clinical diseases. Meanwhile, we prospect the future research direction of TTF-1, which might ultimately contribute to the understanding of the pathogenesis of related clinical diseases and the development of new prevention and treatment strategies.

Genome-wide identification of the DNA-binding one zinc finger (Dof) transcription factor gene family and their potential functioning in nitrogen remobilization in tea plant (Camellia sinensis L.)

The DNA-binding one zinc finger (Dof) proteins are a type of transcription factor (TF) specific to plants that are involved in various aspects of their biological processes. Although the Dof TF has been identified in a variety of species, its systemic analysis in tea (Camellia sinensis L.) is lacking and its potential role in this plant's nitrogen remobilization remains unclear. Here, we carried out a genome-wide identification of Dofs in tea and identified 41 CsDofs at the chromosome level. Their genetic structure, conserved motif, and synteny relations were analyzed and visualized, using information extracted from a newly available genome database. Phylogenetic analysis of Dofs from seven different plant species revealed that Dof members could be clustered into six groups. Heatmap analysis showed that CsDofs have similar expression patterns among different species of genus Camellia, and tissue-specific transcription patterns among tea tissues. Expression patterns of CsDofs in mature leaves during the development of new shoots revealed that most of CsDofs peaked on March 24 when the new shoots developed into one bud and two leaves. Furthermore, the expression pattern of CsDof13/16/29/33/34/35/36/39 was significantly correlated with that of genes involved in N metabolism, such as GDH or GS, indicating these CsDofs might play roles in regulating the expression of genes participating in N metabolism of tea. Taken together, this study provides information on the genome-wide identification of the CsDofs’ family in tea, and identifies putative candidates possible involved in N remobilization in mature leaves during the development of new shoots in spring.

Identifying the lungs as a susceptible site for allele-specific regulatory changes associated with type 1 diabetes risk

Type 1 diabetes (T1D) etiology is complex. We developed a machine learning approach that ranked the tissue-specific transcription regulatory effects for T1D SNPs and estimated their relative contributions to conversion to T1D by integrating case and control genotypes (Wellcome Trust Case Control Consortium and UK Biobank) with tissue-specific expression quantitative trait loci (eQTL) data. Here we show an eQTL (rs6679677) associated with changes to AP4B1-AS1 transcript levels in lung tissue makes the largest gene regulatory contribution to the risk of T1D development. Luciferase reporter assays confirmed allele-specific enhancer activity for the rs6679677 tagged locus in lung epithelial cells (i.e. A549 cells; C > A reduces expression, p = 0.005). Our results identify tissue-specific eQTLs for SNPs associated with T1D. The strongest tissue-specific eQTL effects were in the lung and may help explain associations between respiratory infections and risk of islet autoantibody seroconversion in young children.

Identification of the transcription factor MAZ as a regulator of erythropoiesis

AbstractErythropoiesis requires a combination of ubiquitous and tissue-specific transcription factors (TFs). Here, through DNA affinity purification followed by mass spectrometry, we have identified the widely expressed protein MAZ (Myc-associated zinc finger) as a TF that binds to the promoter of the erythroid-specific human α-globin gene. Genome-wide mapping in primary human erythroid cells revealed that MAZ also occupies active promoters as well as GATA1-bound enhancer elements of key erythroid genes. Consistent with an important role during erythropoiesis, knockdown of MAZ reduces α-globin expression in K562 cells and impairs differentiation in primary human erythroid cells. Genetic variants in the MAZ locus are associated with changes in clinically important human erythroid traits. Taken together, these findings reveal the zinc-finger TF MAZ to be a previously unrecognized regulator of the erythroid differentiation program.

Overexpression of OsERF83, a Vascular Tissue-Specific Transcription Factor Gene, Confers Drought Tolerance in Rice.

Abiotic stresses severely affect plant growth and productivity. To cope with abiotic stresses, plants have evolved tolerance mechanisms that are tightly regulated by reprogramming transcription factors (TFs). APETALA2/ethylene-responsive factor (AP2/ERF) transcription factors are known to play an important role in various abiotic stresses. However, our understanding of the molecular mechanisms remains incomplete. In this study, we identified the role of OsERF83, a member of the AP2/ERF transcription factor family, in response to drought stress. OsERF83 is a transcription factor localized to the nucleus and induced in response to various abiotic stresses, such as drought and abscisic acid (ABA). Overexpression of OsERF83 in transgenic plants (OsERF83OX) significantly increased drought tolerance, with higher photochemical efficiency in rice. OsERF83OX was also associated with growth retardation, with reduced grain yields under normal growth conditions. OsERF83 is predominantly expressed in the vascular tissue of all organs. Transcriptome analysis revealed that OsERF83 regulates drought response genes, which are related to the transporter (OsNPF8.10, OsNPF8.17, OsLH1), lignin biosynthesis (OsLAC17, OsLAC10, CAD8D), terpenoid synthesis (OsTPS33, OsTPS14, OsTPS3), cytochrome P450 family (Oscyp71Z4, CYP76M10), and abiotic stress-related genes (OsSAP, OsLEA14, PCC13-62). OsERF83 also up-regulates biotic stress-associated genes, including PATHOGENESIS-RELATED PROTEIN (PR), WALL-ASSOCIATED KINASE (WAK), CELLULOSE SYNTHASE-LIKE PROTEIN E1 (CslE1), and LYSM RECEPTOR-LIKE KINASE (RLK) genes. Our results provide new insight into the multiple roles of OsERF83 in the cross-talk between abiotic and biotic stress signaling pathways.

Tissue-specific 5-hydroxymethylcytosine landscape of the human genome

5-Hydroxymethylcytosine (5hmC) is an important epigenetic mark that regulates gene expression. Charting the landscape of 5hmC in human tissues is fundamental to understanding its regulatory functions. Here, we systematically profiled the whole-genome 5hmC landscape at single-base resolution for 19 types of human tissues. We found that 5hmC preferentially decorates gene bodies and outperforms gene body 5mC in reflecting gene expression. Approximately one-third of 5hmC peaks are tissue-specific differentially-hydroxymethylated regions (tsDhMRs), which are deposited in regions that potentially regulate the expression of nearby tissue-specific functional genes. In addition, tsDhMRs are enriched with tissue-specific transcription factors and may rewire tissue-specific gene expression networks. Moreover, tsDhMRs are associated with single-nucleotide polymorphisms identified by genome-wide association studies and are linked to tissue-specific phenotypes and diseases. Collectively, our results show the tissue-specific 5hmC landscape of the human genome and demonstrate that 5hmC serves as a fundamental regulatory element affecting tissue-specific gene expression programs and functions.

Phosphorylation of Pit-1 by cyclin-dependent kinase 5 at serine 126 is associated with cell proliferation and poor prognosis in prolactinomas

Abstract Pit-1 (POU1F1) is a POU-homeodomain transcription factor, and it is one of the most important tissue-specific transcription factors in pituitary development. Cyclin-dependent kinase 5 (CDK5) is a protein kinase that can phosphorylate many key transcription factors, but the mechanism under which CDK5 phosphorylates Pit-1 is unclear. To investigate whether CDK5 can regulate cell proliferation and promote hormone secretion through phosphorylation of Ser126-Pit-1 in prolactinomas, we generated an antibody that specifically recognizes phosphorylated serine at position 126 of Pit-1 (Ser126-Pit-1). We used western blotting to detect the level of Pit-1 phosphorylation and observed the proliferation and apoptosis of GH3 cells with different levels of Pit-1 phosphorylation by clone formation experiments, cell viability assays, and flow cytometry. ELISA was used to measure the level of PRL in the supernatant of GH3 cells. Tissue microarrays and immunohistochemistry were used to evaluate the expression of the phosphorylation level of Ser126-Pit-1 (pSer126-Pit-1) in prolactinomas. Our data indicated that Ser126-Pit-1 is specifically phosphorylated by CDK5 and high-level pSer-126-Pit-1 can promote cell proliferation and PRL secretion. In addition, a higher level of pSer-126-Pit-1 correlates with a worse prognosis in patients with prolactinoma. Our results show that CDK5 mediated Ser126-Pit-1 phosphorylation and regulated prolactinoma progression and PRL secretion.

The Role of lncRNAs in the Stem Phenotype of Pancreatic Ductal Adenocarcinoma.

Pancreatic ductal adenocarcinoma is one of the deadliest tumors. This neoplasia is characterized by an important cellular and phenotypic heterogeneity. In particular, it has been shown that at least two subtypes can be found: basal-like, which presents stem-like properties, and classical. Cancer stem cells have been isolated and characterized from these tumors, showing their dependance on general and tissue-specific stem transcription factors and signaling pathways. Nevertheless, little is known about their tissue microenvironment and cell non-autonomous regulators, such as long-non-coding RNAs. (lncRNAs). In this review, we summarize the current knowledge about the positive and negative effects of lncRNAs in the stemness phenotype of pancreatic ductal adenocarcinoma cancer (PDAC).

Molecular cloning, characterization, and steady-state levels of retinol-binding protein (rbp) genes in response to dietary soybean oil in brown trout (Salmo trutta)

The cDNA of rbp1a, rbp2a, rbp2b, rbp5, and rbp7b genes was cloned and characterized in brown trout, and their mRNA expressions in response to different levels of dietary soybean oil (SBO) were determined along with their growth parameters, muscle fatty acid (FA) profile, and hepatic retinol level. Dietary SBO did not affect the growth parameters, survival, and hepatosomatic index. The FA composition and retinol content of fish reflected those of experimental diets. The retinol-binding proteins (Rbps; except Rbp5) contained lipocalin domains. The phylogenetic analysis showed that brown trout rbps are orthologs of the rbp/RBPs from other vertebrates; this was also supported by identity/similarity analysis. The rbp genes were exhibited to have conserved synteny with zebrafish and human rbps/RBPs. The tissue-specific transcription of duplicated rbp2a and rbp2b genes varied among all the tested tissues. Moreover, the diet including 100% SBO affected the rbp2a and rbp2b mRNA transcriptions differently, and this was explained by the divergence of the transcriptional regulation site of these genes due to the teleost-specific whole genome duplication. rbp1b and rbp7a were lost in the brown trout genome most likely because of nonfunctionalization. The mRNA expression of all the rbp genes studied was affected by dietary SBO. Therefore, their mRNA expressions could be modulated by dietary FAs and retinol content.

Functional annotations of three domestic animal genomes provide vital resources for comparative and agricultural research

Gene regulatory elements are central drivers of phenotypic variation and thus of critical importance towards understanding the genetics of complex traits. The Functional Annotation of Animal Genomes consortium was formed to collaboratively annotate the functional elements in animal genomes, starting with domesticated animals. Here we present an expansive collection of datasets from eight diverse tissues in three important agricultural species: chicken (Gallus gallus), pig (Sus scrofa), and cattle (Bos taurus). Comparative analysis of these datasets and those from the human and mouse Encyclopedia of DNA Elements projects reveal that a core set of regulatory elements are functionally conserved independent of divergence between species, and that tissue-specific transcription factor occupancy at regulatory elements and their predicted target genes are also conserved. These datasets represent a unique opportunity for the emerging field of comparative epigenomics, as well as the agricultural research community, including species that are globally important food resources.

Genetic Drift Dominates Genome-Wide Regulatory Evolution Following an Ancient Whole-Genome Duplication in Atlantic Salmon.

Whole-genome duplications (WGD) have been considered as springboards that potentiate lineage diversification through increasing functional redundancy. Divergence in gene regulatory elements is a central mechanism for evolutionary diversification, yet the patterns and processes governing regulatory divergence following events that lead to massive functional redundancy, such as WGD, remain largely unknown. We studied the patterns of divergence and strength of natural selection on regulatory elements in the Atlantic salmon (Salmo salar) genome, which has undergone WGD 100–80 Ma. Using ChIPmentation, we first show that H3K27ac, a histone modification typical to enhancers and promoters, is associated with genic regions, tissue-specific transcription factor binding motifs, and with gene transcription levels in immature testes. Divergence in transcription between duplicated genes from WGD (ohnologs) correlated with difference in the number of proximal regulatory elements, but not with promoter elements, suggesting that functional divergence between ohnologs after WGD is mainly driven by enhancers. By comparing H3K27ac regions between duplicated genome blocks, we further show that a longer polyploid state post-WGD has constrained regulatory divergence. Patterns of genetic diversity across natural populations inferred from resequencing indicate that recent evolutionary pressures on H3K27ac regions are dominated by largely neutral evolution. In sum, our results suggest that post-WGD functional redundancy in regulatory elements continues to have an impact on the evolution of the salmon genome, promoting largely neutral evolution of regulatory elements despite their association with transcription levels. These results highlight a case where genome-wide regulatory evolution following an ancient WGD is dominated by genetic drift.