Putative Regulatory Regions Research Articles

Whole-genome DNA methylomes of tree shrew brains reveal conserved and divergent roles of DNA methylation on sex chromosome regulation

BackgroundThe tree shrew (Tupaia belangeri) is a promising emerging model organism in biomedical studies, notably due to their evolutionary proximity to primates. To enhance our understanding of how DNA methylation is implicated in regulation of gene expression and the X chromosome inactivation (XCI) in tree shrew brains, here we present their first genome-wide, single-base-resolution methylomes integrated with transcriptomes from prefrontal cortices.ResultsGenome-wide relationships between DNA methylation and gene expression are consistent with those in other mammals. Interestingly, we observed a clear and significant global reduction (hypomethylation) of DNA methylation across the entire female X chromosome compared to male X. Female hypomethylation does not directly contribute to the gene silencing of the inactivated X chromosome nor does it significantly drive sex-specific gene expression in tree shrews. However, we identified a putative regulatory region in the 5′ end of the X-inactive-specific transcript (Xist) gene, whose pattern of differential DNA methylation strongly relate to its sex-differential expression in tree shrews. Furthermore, differential methylation of this region is conserved across different species. We also provide evidence suggesting that the observed difference between human and tree shrew X-linked promoter methylation is associated with the difference in genomic CpG contents.ConclusionsOur study offers novel information on genomic DNA methylation of tree shrews as well as insights into the evolution of sex chromosome regulation in mammals. Specifically, we show conserved role of DNA methylation in regulation of Xist expression and propose genomic CpG contents as a factor in driving sex-differential DNA methylation of X-linked promoters.

Common variants increase risk for congenital diaphragmatic hernia within the context of de novo variants

Congenital diaphragmatic hernia (CDH) is a severe congenital anomaly often accompanied by other structural anomalies and/or neurobehavioral manifestations. Rare de novo protein-coding variants and copy-number variations contribute to CDH in the population. However, most individuals with CDH remain genetically undiagnosed. Here, we perform integrated de novo and common-variant analyses using 1,469 CDH individuals, including 1,064 child-parent trios and 6,133 ancestry-matched, unaffected controls for the genome-wide association study. We identify candidate CDH variants in 15 genes, including eight novel genes, through deleterious de novo variants. We further identify two genomic loci contributing to CDH risk through common variants with similar effect sizes among Europeans and Latinx. Both loci are in putative transcriptional regulatory regions of developmental patterning genes. Estimated heritability in common variants is ∼19%. Strikingly, there is no significant difference in estimated polygenic risk scores between isolated and complex CDH or between individuals harboring deleterious de novo variants and individuals without these variants. The data support a polygenic model as part of the CDH genetic architecture.

Variant-to-function mapping of late-onset Alzheimer's disease GWAS signals in human microglial cell models implicates RTFDC1 at the CASS4 locus.

Late-onset Alzheimer's disease (LOAD) research has principally focused on neurons over the years due to their known role in the production of amyloid beta plaques and neurofibrillary tangles. In contrast, recent genomic studies of LOAD have implicated microglia as culprits of the prolonged inflammation exacerbating the neurodegeneration observed in patient brains. Indeed, recent LOAD genome-wide association studies (GWAS) have reported multiple loci near genes related to microglial function, including TREM2, ABI3, and CR1. However, GWAS alone cannot pinpoint underlying causal variants or effector genes at such loci, as most signals reside in non-coding regions of the genome and could presumably confer their influence frequently via long-range regulatory interactions. We elected to carry out a combination of ATAC-seq and high-resolution promoter-focused Capture-C in two human microglial cell models (iPSC-derived microglia and HMC3) in order to physically map interactions between LOAD GWAS-implicated candidate causal variants and their corresponding putative effector genes. Notably, we observed consistent evidence that rs6024870 at the GWAS CASS4 locus contacted the promoter of nearby gene, RTFDC1. We subsequently observed a directionallly consistent decrease in RTFDC1 expression with the the protective minor A allele of rs6024870 via both luciferase assays in HMC3 cells and expression studies in primary human microglia. Through CRISPR-Cas9-mediated deletion of the putative regulatory region harboring rs6024870 in HMC3 cells, we observed increased pro-inflammatory cytokine secretion and decreased DNA double strand break repair related, at least in part, to RTFDC1 expression levels. Our variant-to-function approach therefore reveals that the rs6024870-harboring regulatory element at the LOAD 'CASS4' GWAS locus influences both microglial inflammatory capacity and DNA damage resolution, along with cumulative evidence implicating RTFDC1 as a novel candidate effector gene.

Transcriptional and epigenetic characterization of a new in vitro platform to model the formation of human pharyngeal endoderm

BackgroundThe Pharyngeal Endoderm (PE) is an extremely relevant developmental tissue, serving as the progenitor for the esophagus, parathyroids, thyroids, lungs, and thymus. While several studies have highlighted the importance of PE cells, a detailed transcriptional and epigenetic characterization of this important developmental stage is still missing, especially in humans, due to technical and ethical constraints pertaining to its early formation.ResultsHere we fill this knowledge gap by developing an in vitro protocol for the derivation of PE-like cells from human Embryonic Stem Cells (hESCs) and by providing an integrated multi-omics characterization. Our PE-like cells robustly express PE markers and are transcriptionally homogenous and similar to in vivo mouse PE cells. In addition, we define their epigenetic landscape and dynamic changes in response to Retinoic Acid by combining ATAC-Seq and ChIP-Seq of histone modifications. The integration of multiple high-throughput datasets leads to the identification of new putative regulatory regions and to the inference of a Retinoic Acid-centered transcription factor network orchestrating the development of PE-like cells.ConclusionsBy combining hESCs differentiation with computational genomics, our work reveals the epigenetic dynamics that occur during human PE differentiation, providing a solid resource and foundation for research focused on the development of PE derivatives and the modeling of their developmental defects in genetic syndromes.

Whole-genome DNA methylomes of Tree shrew brains reveal conserved and divergent roles of DNA methylation on sex chromosome regulation.

The tree shrew (Tupaia belangeri) is a promising emerging model organism in biomedical studies, notably due to their evolutionary proximity to primates. To enhance our understanding of how DNA methylation is implicated in regulation of gene expression and the X chromosome inactivation (XCI) in tree shrew brains, here we present their first genome-wide, single-base-resolution methylomes integrated with transcriptomes from prefrontal cortices. We discovered both divergent and conserved features of tree shrew DNA methylation compared to that of other mammals. DNA methylation levels of promoter and gene body regions are negatively correlated with gene expression, consistent with patterns in other mammalian brains studied. Comparing DNA methylation patterns of the female and male X chromosomes, we observed a clear and significant global reduction (hypomethylation) of DNA methylation across the entire X chromosome in females. Our data suggests that the female X hypomethylation does not directly contribute to the gene silencing of the inactivated X chromosome nor does it significantly drive sex-specific gene expression of tree shrews. However, we identified a putative regulatory region in the 5' end of the X inactive specific transcript (Xist) gene, a key gene for XCI, whose pattern of differential DNA methylation strongly relate to its differential expression between male and female tree shrews. We show that differential methylation of this region is conserved across different species. Moreover, we provide evidence suggesting that the observed difference between human and tree shrew X-linked promoter methylation is associated with the difference in genomic CpG contents. Our study offers novel information on genomic DNA methylation of tree shrews, as well as insights into the evolution of X chromosome regulation in mammals.

Predicting gene expression state and prioritizing putative enhancers using 5hmC signal

BackgroundLike its parent base 5-methylcytosine (5mC), 5-hydroxymethylcytosine (5hmC) is a direct epigenetic modification of cytosines in the context of CpG dinucleotides. 5hmC is the most abundant oxidized form of 5mC, generated through the action of TET dioxygenases at gene bodies of actively-transcribed genes and at active or lineage-specific enhancers. Although such enrichments are reported for 5hmC, to date, predictive models of gene expression state or putative regulatory regions for genes using 5hmC have not been developed.ResultsHere, by using only 5hmC enrichment in genic regions and their vicinity, we develop neural network models that predict gene expression state across 49 cell types. We show that our deep neural network models distinguish high vs low expression state utilizing only 5hmC levels and these predictive models generalize to unseen cell types. Further, in order to leverage 5hmC signal in distal enhancers for expression prediction, we employ an Activity-by-Contact model and also develop a graph convolutional neural network model with both utilizing Hi-C data and 5hmC enrichment to prioritize enhancer-promoter links. These approaches identify known and novel putative enhancers for key genes in multiple immune cell subsets.ConclusionsOur work highlights the importance of 5hmC in gene regulation through proximal and distal mechanisms and provides a framework to link it to genome function. With the recent advances in 6-letter DNA sequencing by short and long-read techniques, profiling of 5mC and 5hmC may be done routinely in the near future, hence, providing a broad range of applications for the methods developed here.

Protein disulfide isomerase-mediated transcriptional upregulation of Nox1 contributes to vascular dysfunction in hypertension.

Nox1 signaling is a causal key element in arterial hypertension. Recently, we identified protein disulfide isomerase A1 (PDI) as a novel regulatory protein that regulates Nox1 signaling in VSMCs. Spontaneously hypertensive rats (SHR) have increased levels of PDI in mesenteric resistance arteries compared with Wistar controls; however, its consequences remain unclear. Herein, we investigated the role of PDI in mediating Nox1 transcriptional upregulation and its effects on vascular dysfunction in hypertension. We demonstrate that PDI contributes to the development of hypertension via enhanced transcriptional upregulation of Nox1 in vascular smooth muscle cells (VSMCs). We show for the first time that PDI sulfenylation by hydrogen peroxide contributes to EGFR activation in hypertension via increased shedding of epidermal growth factor-like ligands. PDI also increases intracellular calcium levels, and contractile responses induced by ANG II. PDI silencing or pharmacological inhibition in VSMCs significantly decreases EGFR activation and Nox1 transcription. Overexpression of PDI in VSMCs enhances ANG II-induced EGFR activation and ATF1 translocation to the nucleus. Mechanistically, PDI increases ATF1-induced Nox1 transcription and enhances the contractile responses to ANG II. Herein we show that ATF1 binding to Nox1 transcription putative regulatory regions is augmented by PDI. Altogether, we provide evidence that HB-EGF in SHR resistance vessels promotes the nuclear translocation of ATF1, under the control of PDI, and thereby induces Nox1 gene expression and increases vascular reactivity. Thus, PDI acts as a thiol redox-dependent enhancer of vascular dysfunction in hypertension and could represent a novel therapeutic target for the treatment of this disease.

179 Identification of Regulatory Mutations Affecting Complex Traits

Abstract Understanding genetic bases controlling complex traits of economic importance is fundamental to improve breeding programs. Recent efforts from FAANG consortium resulted in significant progress in the identification of regulatory regions in farm animals. However, most studies in cattle have focused on Bos taurus, and limited information is available for Bos indicus. Considering that Bos indicus is the predominant subspecies in tropical regions, the objective of this work was to identify regulatory variants associated with meat quality by combining information from whole genome sequence (WGS), gene expression (RNA-seq), and open chromatin regions (ATAC-seq). For genome wide association study (GWAS), an imputed panel of SNPs was obtained from WGS of 26 sires, a Bovine HD 770K panel of 374 progenies, and SNPs called from RNA-seq data of Longissimus thoracis muscle (LT) of 192 progenies. The SNPs were identified using GATK4, and the dataset included almost 4.5 million variants that were pruned to 553,581 tag-SNPs (PLINK, r2 = 0.8). ATAC-seq data was obtained from duplicated samples of LT collected from two animals. Phenotype data included back fat thickness (BFT), intramuscular fat (IMF), ribeye area (REA) and shear force (SF). Association of genotypes (553,581 SNPs) and gene expression (12,991 genes), performed in Matrix eQTL resulted in 36,916 cis-eQTL and 14,408 trans-eQTL (FDR < 0.05). For GWAS analysis of the eQTLs with phenotypic traits, we used a linear model in PLINK program (considering population stratification, contemporary group and hot carcass weight). GWAS resulted in 24 eQTLs associated with IMF, three with BFT, 12 with REA, and eight with SF (FDR < 0.10). ATAC-seq analysis (using nfcore/atacseq pipeline implemented in Nextflow) resulted in 33,734 consensus peaks with a mean size of 2.4 kb. The overlap analyses among GWAS, eQTL, and ATAC-seq positions resulted in two, four and one putative regulatory regions that control gene expression and are associated with BFT, IMF and SF, respectively. Functional studies are still necessary to establish if identified mutations are the causal mutation. In conclusion, integration of GWAS, RNAseq, and ATAC-Seq data can result in identification of regulatory mutations associated with phenotypes of economic importance.

Allele specific binding of histone modifications and a transcription factor does not predict allele specific expression in correlated ChIP-seq peak-exon pairs

Allele specific expression (ASE) is widespread in many species including cows. Therefore, regulatory regions which control gene expression should show cis-regulatory variation which mirrors this differential expression within the animal. ChIP-seq peaks for histone modifications and transcription factors measure activity at functional regions and the height of some peaks have been shown to correlate across tissues with the expression of particular genes, suggesting these peaks are putative regulatory regions. In this study we identified ASE in the bovine genome in multiple tissues and investigated whether ChIP-seq peaks for four histone modifications and the transcription factor CTCF show allele specific binding (ASB) differences in the same tissues. We then investigate whether peak height and gene expression, which correlates across tissues, also correlates within the animal by investigating whether the direction of ASB in putative regulatory regions, mirrors that of the ASE in the genes they are putatively regulating. We found that ASE and ASB were widespread in the bovine genome but vary in extent between tissues. However, even when the height of a peak was positively correlated across tissues with expression of an exon, ASE of the exon and ASB of the peak were in the same direction only half the time. A likely explanation for this finding is that the correlations between peak height and exon expression do not indicate that the height of the peak causes the extent of exon expression, at least in some cases.

Evolutionary dynamics of genome size and content during the adaptive radiation of Heliconiini butterflies

Heliconius butterflies, a speciose genus of Müllerian mimics, represent a classic example of an adaptive radiation that includes a range of derived dietary, life history, physiological and neural traits. However, key lineages within the genus, and across the broader Heliconiini tribe, lack genomic resources, limiting our understanding of how adaptive and neutral processes shaped genome evolution during their radiation. Here, we generate highly contiguous genome assemblies for nine Heliconiini, 29 additional reference-assembled genomes, and improve 10 existing assemblies. Altogether, we provide a dataset of annotated genomes for a total of 63 species, including 58 species within the Heliconiini tribe. We use this extensive dataset to generate a robust and dated heliconiine phylogeny, describe major patterns of introgression, explore the evolution of genome architecture, and the genomic basis of key innovations in this enigmatic group, including an assessment of the evolution of putative regulatory regions at the Heliconius stem. Our work illustrates how the increased resolution provided by such dense genomic sampling improves our power to generate and test gene-phenotype hypotheses, and precisely characterize how genomes evolve.

Conservation of the restricted expression of Brassicaceae Bsister-like genes in seeds requires a transposable element in Arabidopsis thaliana.

Changes in transcription factor binding sites (TFBSs) can alter the spatio-temporal expression pattern and transcript abundance of genes. Loss and gain of TFBSs were shown to cause shifts in expression patterns in numerous cases. However, we know little about the evolution of extended regulatory sequences incorporating many TFBSs. We compare, across the crucifers (Brassicaceae, cabbage family), the sequences between the translated regions of Arabidopsis Bsister- (ABS-) like MADS-box genes (including paralogous GOA-like genes) and the next gene upstream, as an example of family wide evolution of putative regulatory regions (PURRs). ABS-like genes are essential for integument development of ovules and endothelium formation in seeds of Arabidopsis thaliana. A combination of motif based GO-enrichment and reporter gene analysis using A. thaliana as common trans-regulatory environment allows analysis of selected Brassicaceae Bsister gene PURRs. Comparison of TFBS of transcriptionally active ABS-like genes with those of transcriptionally largely inactive GOA-like genes shows that the number of in silico predicted TF binding sites (TFBSs) is similar between paralogs, emphasizing the importance of experimental verification for in silico characterization of TFBS activity and analysis of their evolution. Further, our data show highly conserved expression of Brassicaceae ABS-like genes almost exclusively in the chalazal region of ovules. The Arabidopsis-specific insertion of a transposable element (TE) into the ABS PURRs is required for stabilizing this spatially restricted expression, while other Brassicaceae achieve chalaza-specific expression without TE insertion. We hypothesize that the chalaza-specific expression of ABS is regulated by CREs provided by the TE.

Animal Models for Studying the Genetic Basis of Metabolic Regulation

Modern genetics has developed methods to modify the expression of genes in animals to study the factors responsible for the tissue-specific expression and hormonal and dietary regulation of metabolic processes. As these methods are applied to genes that code for critical proteins in metabolic pathways, a new insight into the control of metabolism is emerging. There are three general approaches currently in use. First, is the introduction of genes into the germ line to create transgenic animal models in which the gene of interest is over-expressed. This model is of particular value for promoter analysis because it is possible to introduce specific mutations into a putative regulatory region of a transgene and study its transcriptional control in the intact animal. Second, the developmental function of a gene product and its effect on various metabolic processes in a mouse can be directly determined by deleting a gene of interest by homologous recombination. Gene “knockout” mice are currently available with deletions in the genes for a variety of transcription factors and other biologically active proteins, permitting a critical analysis of the proteins responsible for the metabolic patterning of the animal. Third, the metabolic role of a gene of interest in a specific tissue can be studied by ablating its mRNA by the introduction of a transgene that codes for antisense mRNA targeted against the gene transcript. Because it is possible to use a transgene with a tissue-specific promoter, this procedure allows the isolation of the metabolic effect to a selected tissue in the transgenic animal. Taken together, these procedures provide a unique set of metabolic models for an in-depth study of metabolic regulation. This review will present examples of selected animal models currently available and will outline the challenge these animals present for investigators in the nutritional sciences.

Patterns of Heterochromatin Transitions Linked to Changes in the Expression of Plasmodium falciparum Clonally Variant Genes.

The survival of malaria parasites in the changing human blood environment largely depends on their ability to alter gene expression by epigenetic mechanisms. The active state of Plasmodium falciparum clonally variant genes (CVGs) is associated with euchromatin characterized by the histone mark H3K9ac, whereas the silenced state is characterized by H3K9me3-based heterochromatin. Expression switches are linked to euchromatin-heterochromatin transitions, but these transitions have not been characterized for the majority of CVGs. To define the heterochromatin distribution patterns associated with the alternative transcriptional states of CVGs, we compared H3K9me3 occupancy at a genome-wide level among several parasite subclones of the same genetic background that differed in the transcriptional state of many CVGs. We found that de novo heterochromatin formation or the complete disruption of a heterochromatin domain is a relatively rare event, and for the majority of CVGs, expression switches can be explained by the expansion or retraction of heterochromatin domains. We identified different modalities of heterochromatin changes linked to transcriptional differences, but despite this complexity, heterochromatin distribution patterns generally enable the prediction of the transcriptional state of specific CVGs. We also found that in some subclones, several var genes were simultaneously in an active state. Furthermore, the heterochromatin levels in the putative regulatory region of the gdv1 antisense noncoding RNA, a regulator of sexual commitment, varied between parasite lines with different sexual conversion rates. IMPORTANCE The malaria parasite P. falciparum is responsible for more than half a million deaths every year. P. falciparum clonally variant genes (CVGs) mediate fundamental host-parasite interactions and play a key role in parasite adaptation to fluctuations in the conditions of the human host. The expression of CVGs is regulated at the epigenetic level by changes in the distribution of a type of chromatin called heterochromatin. Here, we describe at a genome-wide level the changes in the heterochromatin distribution associated with the different transcriptional states of CVGs. Our results also reveal a likely role for heterochromatin at a particular locus in determining the parasite investment in transmission to mosquitoes. Additionally, this data set will enable the prediction of the transcriptional state of CVGs from epigenomic data, which is important for the study of parasite adaptation to the conditions of the host in natural malaria infections.

Genetic variation in histone modifications and gene expression identifies regulatory variants in the mammary gland of cattle.

Causal variants for complex traits, such as eQTL are often found in non-coding regions of the genome, where they are hypothesised to influence phenotypes by regulating gene expression. Many regulatory regions are marked by histone modifications, which can be assayed by chromatin immunoprecipitation followed by sequencing (ChIP-seq). Sequence reads from ChIP-seq form peaks at putative regulatory regions, which may reflect the amount of regulatory activity at this region. Therefore, eQTL which are also associated with differences in histone modifications are excellent candidate causal variants. We assayed the histone modifications H3K4Me3, H3K4Me1 and H3K27ac and mRNA in the mammary gland of up to 400 animals. We identified QTL for peak height (histone QTL), exon expression (eeQTL), allele specific expression (aseQTL) and allele specific binding (asbQTL). By intersecting these results, we identify variants which may influence gene expression by altering regulatory regions of the genome, and may be causal variants for other traits. Lastly, we find that these variants are found in putative transcription factor binding sites, identifying a mechanism for the effect of many eQTL. We find that allele specific and traditional QTL analysis often identify the same genetic variants and provide evidence that many eQTL are regulatory variants which alter activity at regulatory regions of the bovine genome. Our work provides methodological and biological updates on how regulatory mechanisms interplay at multi-omics levels.

Identification of NHLRC1 as a Novel AKT Activator from a Lung Cancer Epigenome-Wide Association Study (EWAS).

Changes in DNA methylation identified by epigenome-wide association studies (EWAS) have been recently linked to increased lung cancer risk. However, the cellular effects of these differentially methylated positions (DMPs) are often unclear. Therefore, we investigated top differentially methylated positions identified from an EWAS study. This included a putative regulatory region of NHLRC1. Hypomethylation of this gene was recently linked with decreased survival rates in lung cancer patients. HumanMethylation450 BeadChip array (450K) analysis was performed on 66 lung cancer case-control pairs from the European Prospective Investigation into Cancer and Nutrition Heidelberg lung cancer EWAS (EPIC HD) cohort. DMPs identified in these pre-diagnostic blood samples were then investigated for differential DNA methylation in lung tumor versus adjacent normal lung tissue from The Cancer Genome Atlas (TCGA) and replicated in two independent lung tumor versus adjacent normal tissue replication sets with MassARRAY. The EPIC HD top hypermethylated DMP cg06646708 was found to be a hypomethylated region in multiple data sets of lung tumor versus adjacent normal tissue. Hypomethylation within this region caused increased mRNA transcription of the closest gene NHLRC1 in lung tumors. In functional assays, we demonstrate attenuated proliferation, viability, migration, and invasion upon NHLRC1 knock-down in lung cancer cells. Furthermore, diminished AKT phosphorylation at serine 473 causing expression of pro-apoptotic AKT-repressed genes was detected in these knock-down experiments. In conclusion, this study demonstrates the powerful potential for discovery of novel functional mechanisms in oncogenesis based on EWAS DNA methylation data. NHLRC1 holds promise as a new prognostic biomarker for lung cancer survival and prognosis, as well as a target for novel treatment strategies in lung cancer patients.

Identification of candidate enhancers controlling the transcriptome during the formation of interphalangeal joints

The formation of the synovial joint begins with the visible emergence of a stripe of densely packed mesenchymal cells located between distal ends of the developing skeletal anlagen called the interzone. Recently the transcriptome of the early synovial joint was reported. Knowledge about enhancers would complement these data and lead to a better understanding of the control of gene transcription at the onset of joint development. Using ChIP-sequencing we have mapped the H3-signatures H3K27ac and H3K4me1 to locate regulatory elements specific for the interzone and adjacent phalange, respectively. This one-stage atlas of candidate enhancers (CEs) was used to map the association between these respective joint tissue specific CEs and biological processes. Subsequently, integrative analysis of transcriptomic data and CEs identified new putative regulatory elements of genes expressed in interzone (e.g., GDF5, BMP2 and DACT2) and phalange (e.g., MATN1, HAPLN1 and SNAI1). We also linked such CEs to genes known as crucial in synovial joint hypermobility and osteoarthritis, as well as phalange malformations. These analyses show that the CE atlas can serve as resource for identifying, and as starting point for experimentally validating, putative disease-causing genomic regulatory regions in patients with synovial joint dysfunctions and/or phalange disorders, and enhancer-controlled synovial joint and phalange formation.

Transcription Factor Affinity at Dorsal Root Ganglia Enhancers after Nerve Injury

Misregulation of gene expression in persistent pain supports a role for context-dependent regulation of gene expression in primary sensory neurons through coordinated transcription factor binding to cis-regulatory elements. In this study, we used a multiomic approach to identify regulatory regions and transcription factors (TFs) that may be involved in the development of nerve injury-induced persistent pain. Adult Sprague Dawley rats were randomized to Chronic Constriction Injury (CCI) to the sciatic nerve or no treatment (naïve). The ipsilateral L4-L6 dorsal root ganglia (DRG)s were removed on Day 14 and used for ChIP-seq for H3K4me1, ATAC-seq, and RNA-seq. Computational integration of this sequencing data generated a comprehensive map of chromatin accessibility changes in the DRG after CCI at cis-regulatory elements. We then used luciferase assays to determine the functional significance of these sequences and determined their ability to bind to TFs using a custom human TF array and EMSA. In vitro validation of binding was confirmed by ChIP-PCR and physiological relevance informed by shRNA knockdown in the 50B11 (nociceptor) cell line. Of 58,446 putative regulatory regions identified genome-wide, 2145 (3.67%) were differentially accessible after CCI. Gene ontology analysis identified molecular functions and biological processes associated with neuronal activation and synaptic signaling. Luciferase assays confirmed their regulatory potential in 50B11 cells which suggests that altered chromatin structure at these regions contributes to the development of mechanical hypersensitivity after nerve injury. Several TFs showed significant binding to one of these regions. CEBPG binding was then confirmed by EMSA and ChIP-qPCR with a fold enrichment of 2.72 ±0.47 over IgG. RNA-seq then identified pathways altered by CEBPG shRNA knockdown. Our data provides a comprehensive map of chromatin accessibility changes in the DRG after CCI and emphasizes the importance of chromatin structure at cis-regulatory elements in persistent pain. This study was supported by grants from National Institutes of Health (Bethesda, Maryland, USA) F32NR015728 (KES), KL2TR003108 (KES), NS110598 (YG), NS117761 (YG), R01GM118760 (SDT), the Arkansas Children's Research Institute (KES), the Arkansas Breast Cancer Research Program (KES) as well as a seed grant from the Johns Hopkins Blaustein Pain Research Fund (SDT).

OsMADS14 and NF-YB1 cooperate in the direct activation of OsAGPL2 and Waxy during starch synthesis in rice endosperm.

Starch synthesis makes a dramatic contribution to the yield and nutritional value of cereal crops. Although several starch synthesis enzymes and related regulators have been reported, the underlying regulatory mechanisms of starch synthesis remain largely unknown. OsMADS14 is a FRUITFULL (FUL)-like MADS-box gene in rice (Oryza sativa). Here we show that two null mutations of OsMADS14 result in a shrunken and chalky grain phenotype. It is caused by obviously defective compound starch granules and a significantly reduced content of both total starch and amylose in the endosperm. Transcriptomic profiling analyses revealed that the loss-of-function of OsMADS14 leads to significantly downregulated expression of many core starch synthesis genes, including OsAGPL2 and Waxy. Both in vitro and in vivo assays demonstrate that the OsMADS14 protein directly binds to stretches of DNA with a CArG-box consensus in the putative regulatory regions of OsAGPL2 and Waxy. Protein-protein interaction experiments also suggest that OsMADS14 interacts with nuclear factor NF-YB1 to promote the transcription of OsAGPL2 and Waxy. Our study thus demonstrates that OsMADS14 plays an essential role in the synthesis of storage starch and provides novel insights into the underlying molecular mechanism that may be used to improve rice cultivars by molecular breeding.

Testis-specific expression pattern of dmrt1 and its putative regulatory region in the sea urchin (Mesocentrotus nudus).

Sea urchin (Mseocentrotus nudus) is an economically important mariculture species in several Asian countries. The growth rate and immunocompetence differ by sex in this species. However, the mechanisms of sex determination in M. nudus have remain unclear. In the present study, we focus on the dmrt1 gene of M. nudus (Mndmrt1) to investigate its dynamic expression pattern during different developmental stages. Real-time quantitative PCR (RT-qPCR) revealed that Mndmrt1 exhibits testis-specific expression and undetectable during the whole embryogenesis. With the development of ontogenetic, Mndmrt1 transcripts are first detected at 9months post-fertilization (mpf). In addition, both the transcripts and protein of Mndmrt1 gene were specifically expressed in spermatogonia and spermatocytes, indicating that it might be a male germ cells marker in sea urchin. Significantly, the 1441bp promoter sequence of Mndmrt1 gene was obtained by DNA walking, and one positive regulatory region at -1197/ -968 in the promoter, as well as one negative regulatory region at -1441/ -1198 have been identified by promoter activity analysis. Moreover, two regulatory regions contain multiple putative binding sites for transcription factors, including Sp1, Egr1, Sox5, CEBP, GATA and SRY. These findings suggest that Mndmrt1 may be related to testis differentiation and spermatogenesis in sea urchin and will provide an insight into understanding the regulatory mechanism of the dmrt1 gene.

Profiling of open chromatin in developing pig (Sus scrofa) muscle to identify regulatory regions.

There is very little information about how the genome is regulated in domestic pigs (Sus scrofa). This lack of knowledge hinders efforts to define and predict the effects of genetic variants in pig breeding programs. To address this knowledge gap, we need to identify regulatory sequences in the pig genome starting with regions of open chromatin. We used the “Improved Protocol for the Assay for Transposase-Accessible Chromatin (Omni-ATAC-Seq)” to identify putative regulatory regions in flash-frozen semitendinosus muscle from 24 male piglets. We collected samples from the smallest-, average-, and largest-sized male piglets from each litter through five developmental time points. Of the 4661 ATAC-Seq peaks identified that represent regions of open chromatin, >50% were within 1 kb of known transcription start sites. Differential read count analysis revealed 377 ATAC-Seq defined genomic regions where chromatin accessibility differed significantly across developmental time points. We found regions of open chromatin associated with downregulation of genes involved in muscle development that were present in small-sized fetal piglets but absent in large-sized fetal piglets at day 90 of gestation. The dataset that we have generated provides a resource for studies of genome regulation in pigs and contributes valuable functional annotation information to filter genetic variants for use in genomic selection in pig breeding programs.