Epigenetics is one of the most rapidly expanding fields in biomedical research, and the popularity of the high-throughput next-generation sequencing (NGS) highlights the accelerating speed of epigenomics discovery over the past decade. Epigenetics studies the heritable phenotypes resulting from chromatin changes but without alteration on DNA sequence. Epigenetic factors and their interactive network regulate almost all of the fundamental biological procedures, and incorrect epigenetic information may lead to complex diseases. A comprehensive understanding of epigenetic mechanisms, their interactions, and alterations in health and diseases genome widely has become a priority in biological research. Bioinformatics is expected to make a remarkable contribution for this purpose, especially in processing and interpreting the large-scale NGS datasets. In this review, we introduce the epigenetics pioneering achievements in health status and complex diseases; next, we give a systematic review of the epigenomics data generation, summarize public resources and integrative analysis approaches, and finally outline the challenges and future directions in computational epigenomics.

Using accumulating SNP (Single-Nucleotide Polymorphism) data, we performed a genome-wide search for polypeptide hormone ligands showing changes in the mature regions to elucidate genotype/phenotype diversity among various human populations. Neuropeptide S (NPS), a brain peptide hormone highly conserved in vertebrates, has diverse physiological effects on anxiety, fear, hyperactivity, food intake, and sleeping time through its cognate receptor-NPSR. Here, we report a SNP rs4751440 (L(6)-NPS) causing non-synonymous substitution on the 6(th) position (V to L) of the NPS mature peptide region. L(6)-NPS has a higher allele frequency in Europeans than other populations and probably originated from European ancestors ~25,000 yrs ago based on haplotype analysis and Approximate Bayesian Computation. Functional analyses indicate that L(6)-NPS exhibits a significant lower bioactivity than the wild type NPS, with ~20-fold higher EC50 values in the stimulation of NPSR. Additional evolutionary and mutagenesis studies further demonstrate the importance of the valine residue in the 6(th) position for NPS functions. Given the known physiological roles of NPS receptor in inflammatory bowel diseases, asthma pathogenesis, macrophage immune responses, and brain functions, our study provides the basis to elucidate NPS evolution and signaling diversity among human populations.

Stable introns and intronic fragments make up the largest population of RNA in the oocyte nucleus of the frog Xenopus tropicalis. These stable intronic sequence RNAs (sisRNAs) persist through the onset of zygotic transcription when synchronous cell division has ended, and the developing embryo consists of approximately 8000 cells. Despite their abundance, the sequence properties and biological function of sisRNAs are just beginning to be understood.

The occurrence of gene duplication/amplification (GDA) provide potential material for adaptive evolution with environmental stress. Several molecular models have been proposed to explain GDA, recombination via short stretches of sequence similarity plays a crucial role. By screening genomes for such events, we propose a "SRS (short repeated sequence) *N + unit + SRS*N" amplified unit under USCE (unequal sister-chromatid exchange) for tandem amplification mediated by SRS with different repeat numbers in eukaryotes. The amplified units identified from 2131 well-organized amplification events that generate multi gene/element copy amplified with subsequent adaptive evolution in the respective species. Genomic data we analyzed showed dynamic changes among related species or subspecies or plants from different ecotypes/strains. This study clarifies the characteristics of variable copy number SRS on both sides of amplified unit under USCE mechanism, to explain well-organized gene tandem amplification under environmental stress mediated by SRS in all eukaryotes.

The interplay between epigenetic modifications and chromatin structure are integral to our understanding of genome function. Methylation of cytosine (5mC) at CG dinucleotides, traditionally associated with transcriptional repression, is the most highly studied chemical modification of DNA, occurring at over 70% of all CG dinucleotides in the genome. Hypomethylated regions (HMRs) often occur in CG islands (CGIs), however, they also occur outside of CGIs and function as cell-type specific enhancers. During the process of differentiation, reorganization of chromatin and nucleosome arrangement at regulatory regions is thought to occur in order for the establishment of cell-type specific transcriptional programs. However, the specifics regarding the organization of nucleosomes at HMRs and the potential mechanisms regulating nucleosome occupancy in these regions are unknown. Here, we have investigated nucleosome organization around hypomethylated regions (HMRs) identified in two mouse primary cells.

Previously, we identified 8-bps long DNA sequences (8-mers) that localize in human proximal promoters and grouped them into known transcription factor binding sites (TFBS). We now examine split 8-mers consisting of two 4-mers separated by 1-bp to 30-bps (X(4)-N(1-30)-X(4)) to identify pairs of TFBS that localize in proximal promoters at a precise distance. These include two overlapping TFBS: the ETS⇔ETS motif ((C/G)CCGGAAGCGGAA) and the ETS⇔CRE motif ((C/G)CGGAAGTGACGTCAC). The nucleotides in bold are part of both TFBS. Molecular modeling shows that the ETS⇔CRE motif can be bound simultaneously by both the ETS and the B-ZIP domains without protein-protein clashes. The electrophoretic mobility shift assay (EMSA) shows that the ETS protein GABPα and the B-ZIP protein CREB preferentially bind to the ETS⇔CRE motif only when the two TFBS overlap precisely. In contrast, the ETS domain of ETV5 and CREB interfere with each other for binding the ETS⇔CRE. The 11-mer (CGGAAGTGACG), the conserved part of the ETS⇔CRE motif, occurs 226 times in the human genome and 83% are in known regulatory regions. In vivo GABPα and CREB ChIP-seq peaks identified the ETS⇔CRE as the most enriched motif occurring in promoters of genes involved in mRNA processing, cellular catabolic processes, and stress response, suggesting that a specific class of genes is regulated by this composite motif.

The failure of novel therapies effective in preclinical animal models largely reflects the fact that current models do not really mimic the pathological/therapeutic features of glioblastoma (GBM), in which the most effective temozolomide chemoradiotherapy (RT/TMZ) regimen can only slightly extend survival. How to improve RT/TMZ efficacy remains a major challenge in clinic.

In mammals, the cytosine in CG dinucleotides is typically methylated producing 5-methylcytosine (5mC), a chemically less stable form of cytosine that can spontaneously deaminate to thymidine resulting in a T•G mismatched base pair. Unlike other eukaryotes that efficiently repair this mismatched base pair back to C•G, in mammals, 5mCG deamination is mutagenic, sometimes producing TG dinucleotides, explaining the depletion of CG dinucleotides in mammalian genomes. It was suggested that new TG dinucleotides generate genetic diversity that may be critical for evolutionary change. We tested this conjecture by examining the DNA sequence properties of regulatory sequences identified by DNase I hypersensitive sites (DHSs) in human and mouse genomes. We hypothesized that the new TG dinucleotides generate transcription factor binding sites (TFBS) that become tissue-specific DHSs (TS-DHSs). We find that 8-mers containing the CG dinucleotide are enriched in DHSs in both species. However, 8-mers containing a TG and no CG dinucleotide are preferentially enriched in TS-DHSs when compared with 8-mers with neither a TG nor a CG dinucleotide. The most enriched 8-mer with a TG and no CG dinucleotide in tissue-specific regulatory regions in both genomes is the AP-1 motif ( TG: A(C)/GT CA: N), and we find evidence that TG dinucleotides in the AP-1 motif arose from CG dinucleotides. Additional TS-DHS-enriched TFBS containing the TG/CA dinucleotide are the E-Box motif (G CA: GC TG: C), the NF-1 motif (GG CATG: CC), and the GR (glucocorticoid receptor) motif (G-A CATG: T-C). Our results support the suggestion that cytosine methylation is mutagenic in tetrapods producing TG dinucleotides that create TFBS that drive evolution.

Chromatin plays a critical role in regulating transcription factors (TFs) binding to their canonical transcription factor binding sites (TFBS). Recent studies in vertebrates show that many TFs preferentially bind to genomic regions that are well bound by nucleosomes in vitro. Co-occurring secondary motifs sometimes correlated with functional TFBS.

Ischemic stroke, which accounts for 87% of all strokes, constitutes the leading cause of morbidity and mortality in China. Although the genetics and epigenetics of stroke have been extensively investigated, few studies have examined their relationships at different stages of stroke. This study assessed the characteristics of transcriptome changes at different stages of ischemic stroke using a mouse model of transient middle cerebral artery occlusion (tMCAO) and bioinformatics analyses. Cerebral cortex tissues from tMCAO mice at days 1, 3, 7, 14, and 28 were removed for RNA-Seq and small RNA-Seq library construction, sequencing, and bioinformatics analysis. We identified differentially expressed (DE) genes and miRNAs and revealed an association of the up-regulated or down-regulated DEmiRNAs with the correspondingly altered DEgene targets at each time point. In addition, different biological pathways were activated at different time points; thus, three groups of miRNAs were verified that may represent potential clinical biomarkers corresponding to days 1, 3, and 7 after ischemic stroke. Notably, this represents the first functional association of some of these miRNAs with stroke, e.g., miR-2137, miR-874-5p, and miR-5099. Together, our findings lay the foundation for the transition from a single-point, single-drug stroke treatment approach to multiple-time-point multi-drug combination therapies.

Enhancer of zeste homologue 2 (EZH2), the core member of polycomb repressive complex 2 (PRC2), has multiple splicing modes and performs various physiological functions. However, function and mechanism of alternative splicing at Ezh2 exon 3 in reproduction are unknown.

Significantly decreased H3K4 methylation in oocytes from aged mice indicates the important roles of H3K4 methylation in female reproduction. However, how H3K4 methylation regulates oocyte development remains largely unexplored. In this study, it is demonstrated that oocyte-specific expression of dominant negative mutant H3.3-K4M led to a decrease of the level of H3K4 methylation in mouse oocytes, resulting in reduced transcriptional activity and increased DNA methylation in oocytes, disturbed oocyte developmental potency, and fertility of female mice. The impaired expression of genes regulating mitochondrial functions in H3.3-K4M oocytes, accompanied by mitochondrial abnormalities, is further noticed. Moreover, early embryos from H3.3-K4M oocytes show developmental arrest and reduced zygotic genome activation. Collectively, these results show that H3K4 methylation in oocytes is critical to orchestrating gene expression profile, driving the oocyte developmental program, and ensuring oocyte quality. This study also improves understanding of how histone modifications regulate organelle dynamics in oocytes.

SARS-CoV-2 caused the COVID-19 pandemic. COVID-19 may elevate the risk of cognitive impairment and even cause dementia in infected individuals; it may accelerate cognitive decline in elderly patients with dementia, possibly in Alzheimer's disease (AD) patients. However, the mechanisms underlying the interplay between AD and COVID-19 are still unclear. To investigate the underlying mechanisms and associations between AD progression and SARS-CoV-2 infection, we conducted a series of bioinformatics research into SARS-CoV-2-infected cells, COVID-19 patients, AD patients, and SARS-CoV-2-infected AD patients. We identified the common differentially expressed genes (DEGs) in COVID-19 patients, AD patients, and SARS-CoV-2-infected cells, and these DEGs are enriched in certain pathways, such as immune responses and cytokine storms. We constructed the gene interaction network with the signaling transduction module in the center and identified IRF7, STAT1, STAT2, and OAS1 as the hub genes. We also checked the correlations between several key transcription factors and the SARS-CoV-2 and COVID-19 pathway-related genes. We observed that ACE2 expression is positively correlated with IRF7 expression in AD and coronavirus infections, and interestingly, IRF7 is significantly upregulated in response to different RNA virus infections. Further snRNA-seq analysis indicates that NRGN neurons or endothelial cells may be responsible for the increase in ACE2 and IRF7 expression after SARS-CoV-2 infection. The positive correlation between ACE2 and IRF7 expressions is confirmed in the hippocampal formation (HF) of SARS-CoV-2-infected AD patients. Our findings could contribute to the investigation of the molecular mechanisms underlying the interplay between AD and COVID-19 and to the development of effective therapeutic strategies for AD patients with COVID-19.

Repeated extragenic palindromes (REPs) in the enterobacterial genomes are usually composed of individual palindromic units separated by linker sequences. A total of 355 annotated REPs are distributed along the Escherichia coli genome. RNA sequence (RNAseq) analysis showed that almost 80% of the REPs in E. coli are transcribed. The DNA sequence of REP325 showed that it is a cluster of six repeats, each with two palindromic units capable of forming cruciform structures in supercoiled DNA. Here, we report that components of the REP325 element and at least one of its RNA products play a role in bacterial nucleoid DNA condensation. These RNA not only are present in the purified nucleoid but bind to the bacterial nucleoid-associated HU protein as revealed by RNA IP followed by microarray analysis (RIP-Chip) assays. Deletion of REP325 resulted in a dramatic increase of the nucleoid size as observed using transmission electron microscopy (TEM), and expression of one of the REP325 RNAs, nucleoid-associated noncoding RNA 4 (naRNA4), from a plasmid restored the wild-type condensed structure. Independently, chromosome conformation capture (3C) analysis demonstrated physical connections among various REP elements around the chromosome. These connections are dependent in some way upon the presence of HU and the REP325 element; deletion of HU genes and/or the REP325 element removed the connections. Finally, naRNA4 together with HU condensed DNA in vitro by connecting REP325 or other DNA sequences that contain cruciform structures in a pairwise manner as observed by atomic force microscopy (AFM). On the basis of our results, we propose molecular models to explain connections of remote cruciform structures mediated by HU and naRNA4.

An altered pattern of epigenetic modifications, such as DNA methylation and histone modification, is critical to many common human diseases, including cancer. Recently, mitochondrial DNA (mtDNA) was reported to be associated with tumorigenesis through epigenetic regulation of methylation patterns. One of the promising approaches to study DNA methylation and CpG islands (CGIs) is sequencing and analysis of clones derived from the physical library generated by methyl-CpG-binding domain proteins and restriction enzyme MseI. In this study, we observed that the most redundant sequences of 349 clones in a human CGI library were all generated from the human mitochondrial genome. Further analysis indicated that there was a 5,845-bp DNA transfer from mtDNA to chromosome 1, and all the clones should be the products of a 510-bp MseI fragment, which contained a putative CGI of 270 bp. The 510-bp fragment was annotated as part of cytochrome c oxidase subunit II (COXII), and phylogenetic analysis of homologous sequences containing COXII showed three DNA transfer events from mtDNA to nuclear genome, one of which underwent secondary transfer events between different chromosomes. These results may further our understanding of how the mtDNA regulates DNA methylation in the nucleus.

To achieve proper spatiotemporal control of gene expression, transcription factors cooperatively assemble onto specific DNA sequences. The ETS domain protein monomer of GABPα and the B-ZIP domain protein dimer of CREB1 cooperatively bind DNA only when the ETS ((C)/GCGGAA GT: ) and CRE ( GT: GACGTCAC) motifs overlap precisely, producing the ETS↔CRE motif ((C)/GCGGAA GT: GACGTCAC). We designed a Protein Binding Microarray (PBM) with 60-bp DNAs containing four identical sectors, each with 177,440 features that explore the cooperative interactions between GABPα and CREB1 upon binding the ETS↔CRE motif. The DNA sequences include all 15-mers of the form (C)/GCGGA--CG-, the ETS↔CRE motif, and all single nucleotide polymorphisms (SNPs), and occurrences in the human and mouse genomes. CREB1 enhanced GABPα binding to the canonical ETS↔CRE motif CCGGAAGT two-fold, and up to 23-fold for several SNPs at the beginning and end of the ETS motif, which is suggestive of two separate and distinct allosteric mechanisms of cooperative binding. We show that the ETS-CRE array data can be used to identify regions likely cooperatively bound by GABPα and CREB1 in vivo, and demonstrate their ability to identify human genetic variants that might inhibit cooperative binding.

The presence of the blood-brain barrier (BBB) uniquely distinguishes the brain from other organs, and various brain pathologies, including cancer, can disrupt or breach the BBB. The specific implications of BBB alterations in glioma have not been sufficiently clarified.

Recent progress in high-throughput technologies has greatly contributed to the development of DNA methylation profiling. Although there are several reports that describe methylome detection of whole genome bisulfite sequencing, the high cost and heavy demand on bioinformatics analysis prevents its extensive application. Thus, current strategies for the study of mammalian DNA methylomes is still based primarily on genome-wide methylated DNA enrichment combined with DNA microarray detection or sequencing. Methylated DNA enrichment is a key step in a microarray based genome-wide methylation profiling study, and even for future high-throughput sequencing based methylome analysis.

To evaluate the effect of CG methylation on DNA binding of sequence-specific B-ZIP transcription factors (TFs) in a high-throughput manner, we enzymatically methylated the cytosine in the CG dinucleotide on protein binding microarrays. Two Agilent DNA array designs were used. One contained 40,000 features using de Bruijn sequences where each 8-mer occurs 32 times in various positions in the DNA sequence. The second contained 180,000 features with each CG containing 8-mer occurring three times. The first design was better for identification of binding motifs, while the second was better for quantification. Using this novel technology, we show that CG methylation enhanced binding for CEBPA and CEBPB and inhibited binding for CREB, ATF4, JUN, JUND, CEBPD, and CEBPG. The CEBPB|ATF4 heterodimer bound a novel motif CGAT|GCAA 10-fold better when methylated. The electrophoretic mobility shift assay (EMSA) confirmed these results. CEBPB ChIP-seq data using primary female mouse dermal fibroblasts with 50× methylome coverage for each strand indicate that the methylated sequences well-bound on the arrays are also bound in vivo. CEBPB bound 39% of the methylated canonical 10-mers ATTGC|GCAAT in the mouse genome. After ATF4 protein induction by thapsigargin which results in ER stress, CEBPB binds methylated CGAT|GCAA in vivo, recapitulating what was observed on the arrays. This methodology can be used to identify new methylated DNA sequences preferentially bound by TFs, which may be functional in vivo.

The Silkworm Knowledgebase (SilkDB) is a web-based repository for the curation, integration and study of silkworm genetic and genomic data. With the recent accomplishment of a approximately 6X draft genome sequence of the domestic silkworm (Bombyx mori), SilkDB provides an integrated representation of the large-scale, genome-wide sequence assembly, cDNAs, clusters of expressed sequence tags (ESTs), transposable elements (TEs), mutants, single nucleotide polymorphisms (SNPs) and functional annotations of genes with assignments to InterPro domains and Gene Ontology (GO) terms. SilkDB also hosts a set of ESTs from Bombyx mandarina, a wild progenitor of B.mori, and a collection of genes from other Lepidoptera. Comparative analysis results between the domestic and wild silkworm, between B.mori and other Lepidoptera, and between B.mori and the two sequenced insects, fruitfly and mosquito, are displayed by using B.mori genome sequence as a reference framework. Designed as a basic platform, SilkDB strives to provide a comprehensive knowledgebase about the silkworm and present the silkworm genome and related information in systematic and graphical ways for the convenience of in-depth comparative studies. SilkDB is publicly accessible at http://silkworm.genomics.org.cn.