Epigenetics of Genes Preferentially Expressed in Dissimilar Cell Populations: Myoblasts and Cerebellum.

While studying myoblast methylomes and transcriptomes, we found that CDH15 had a remarkable preference for expression in both myoblasts and cerebellum. To understand how widespread such a relationship was and its epigenetic and biological correlates, we systematically looked for genes with similar transcription profiles and analyzed their DNA methylation and chromatin state and accessibility profiles in many different cell populations. Twenty genes were expressed preferentially in myoblasts and cerebellum (Myob/Cbl genes). Some shared DNA hypo- or hypermethylated regions in myoblasts and cerebellum. Particularly striking was ZNF556, whose promoter is hypomethylated in expressing cells but highly methylated in the many cell populations that do not express the gene. In reporter gene assays, we demonstrated that its promoter's activity is methylation sensitive. The atypical epigenetics of ZNF556 may have originated from its promoter's hypomethylation and selective activation in sperm progenitors and oocytes. Five of the Myob/Cbl genes (KCNJ12, ST8SIA5, ZIC1, VAX2, and EN2) have much higher RNA levels in cerebellum than in myoblasts and displayed myoblast-specific hypermethylation upstream and/or downstream of their promoters that may downmodulate expression. Differential DNA methylation was associated with alternative promoter usage for Myob/Cbl genes MCF2L, DOK7, CNPY1, and ANK1. Myob/Cbl genes PAX3, LBX1, ZNF556, ZIC1, EN2, and VAX2 encode sequence-specific transcription factors, which likely help drive the myoblast and cerebellum specificity of other Myob/Cbl genes. This study extends our understanding of epigenetic/transcription associations related to differentiation and may help elucidate relationships between epigenetic signatures and muscular dystrophies or cerebellar-linked neuropathologies.

Promoter-Adjacent DNA Hypermethylation Can Downmodulate Gene Expression: TBX15 in the Muscle Lineage.

TBX15, which encodes a differentiation-related transcription factor, displays promoter-adjacent DNA hypermethylation in myoblasts and skeletal muscle (psoas) that is absent from non-expressing cells in other lineages. By whole-genome bisulfite sequencing (WGBS) and enzymatic methyl-seq (EM-seq), these hypermethylated regions were found to border both sides of a constitutively unmethylated promoter. To understand the functionality of this DNA hypermethylation, we cloned the differentially methylated sequences (DMRs) in CpG-free reporter vectors and tested them for promoter or enhancer activity upon transient transfection. These cloned regions exhibited strong promoter activity and, when placed upstream of a weak promoter, strong enhancer activity specifically in myoblast host cells. In vitro CpG methylation targeted to the DMR sequences in the plasmids resulted in 86−100% loss of promoter or enhancer activity, depending on the insert sequence. These results as well as chromatin epigenetic and transcription profiles for this gene in various cell types support the hypothesis that DNA hypermethylation immediately upstream and downstream of the unmethylated promoter region suppresses enhancer/extended promoter activity, thereby downmodulating, but not silencing, expression in myoblasts and certain kinds of skeletal muscle. This promoter-border hypermethylation was not found in cell types with a silent TBX15 gene, and these cells, instead, exhibit repressive chromatin in and around the promoter. TBX18, TBX2, TBX3 and TBX1 display TBX15-like hypermethylated DMRs at their promoter borders and preferential expression in myoblasts. Therefore, promoter-adjacent DNA hypermethylation for downmodulating transcription to prevent overexpression may be used more frequently for transcription regulation than currently appreciated.

Similarities and Differences in the Acute-Phase Response to SARS-CoV-2 in Rhesus Macaques and African Green Monkeys.

Understanding SARS-CoV-2 immune pathology is critical for the development of effective vaccines and treatments. Here, we employed unbiased serial whole-blood transcriptome profiling by weighted gene network correlation analysis (WGCNA) at pre-specified timepoints of infection to understand SARS-CoV-2-related immune alterations in a cohort of rhesus macaques (RMs) and African green monkeys (AGMs) presenting with varying degrees of pulmonary pathology. We found that the bulk of transcriptional changes occurred at day 3 post-infection and normalized to pre-infection levels by 3 weeks. There was evidence of coordination of transcriptional networks in blood (defined by WGCNA) and the nasopharyngeal SARS-CoV-2 burden as well as the absolute monocyte count. Pathway analysis of gene modules revealed prominent regulation of type I and type II interferon stimulated genes (ISGs) in both RMs and AGMs, with the latter species exhibiting a greater breadth of ISG upregulation. Notably, pathways relating to neutrophil degranulation were enriched in blood of SARS-CoV-2 infected AGMs, but not RMs. Our results elude to hallmark similarities as well as differences in the RM and AGM acute response to SARS-CoV-2 infection, and may help guide the selection of particular NHP species in modeling aspects of COVID-19 disease outcome.

Epigenetics and expression of key genes associated with cardiac fibrosis: NLRP3, MMP2, MMP9, CCN2/CTGF and AGT.

Aims: Excessive inflammatory signaling and pathological remodeling of the extracellular matrix drive cardiac fibrosis and require changes in gene expression. Materials and methods: Using bioinformatics, both tissue-specific expression profiles and epigenomic profiles of some genes critical for cardiac fibrosis were examined, namely, NLRP3, MMP2, MMP9, CCN2/CTGF, AGT (encodes angiotensin II precursors) and hsa-mir-223 (post-transcriptionally regulates NLRP3). Results: In monocytes, neutrophils, fibroblasts, venous cells, liver and brain, enhancers or super-enhancers were found that correlate with high expression of these genes. One enhancer extended into a silent gene neighbor. These enhancers harbored tissue-specific foci of DNA hypomethylation, open chromatin and transcription factor binding. Conclusions: This study identified previously undescribed enhancers containing hypomethylated transcription factor binding subregions that are predicted to regulate expression of these cardiac fibrosis-inducing genes.

Epigenetics of Muscle- and Brain-Specific Expression of KLHL Family Genes.

KLHL and the related KBTBD genes encode components of the Cullin-E3 ubiquitin ligase complex and typically target tissue-specific proteins for degradation, thereby affecting differentiation, homeostasis, metabolism, cell signaling, and the oxidative stress response. Despite their importance in cell function and disease (especially, KLHL40, KLHL41, KBTBD13, KEAP1, and ENC1), previous studies of epigenetic factors that affect transcription were predominantly limited to promoter DNA methylation. Using diverse tissue and cell culture whole-genome profiles, we examined 17 KLHL or KBTBD genes preferentially expressed in skeletal muscle or brain to identify tissue-specific enhancer and promoter chromatin, open chromatin (DNaseI hypersensitivity), and DNA hypomethylation. Sixteen of the 17 genes displayed muscle- or brain-specific enhancer chromatin in their gene bodies, and most exhibited specific intergenic enhancer chromatin as well. Seven genes were embedded in super-enhancers (particularly strong, tissue-specific clusters of enhancers). The enhancer chromatin regions typically displayed foci of DNA hypomethylation at peaks of open chromatin. In addition, we found evidence for an intragenic enhancer in one gene upregulating expression of its neighboring gene, specifically for KLHL40/HHATL and KLHL38/FBXO32 gene pairs. Many KLHL/KBTBD genes had tissue-specific promoter chromatin at their 5' ends, but surprisingly, two (KBTBD11 and KLHL31) had constitutively unmethylated promoter chromatin in their 3' exons that overlaps a retrotransposed KLHL gene. Our findings demonstrate the importance of expanding epigenetic analyses beyond the 5' ends of genes in studies of normal and abnormal gene regulation.

Data showing atherosclerosis-associated differentially methylated regions are often at enhancers.

Atherosclerosis involves phenotypic modulation and transdifferentiation of vascular smooth muscle cells (SMCs). Data are given in tabular or figure format that illustrate genome-wide DNA methylation alterations in atherosclerotic vs. control aorta (athero DMRs). Data based upon publicly available chromatin state profiles are also shown for normal aorta, monocyte, and skeletal muscle tissue-specific DMRs and for aorta-specific chromatin features (enhancer chromatin, promoter chromatin, repressed chromatin, actively transcribed chromatin). Athero hypomethylated and hypermethylated DMRs as well as epigenetic and transcription profiles are described for the following genes: ACTA2, MYH10, MYH11 (SMC-associated genes); SMAD3 (a signaling gene for SMCs and other cell types); CD79B and SH3BP2 (leukocyte-associated genes); and TBX20 and genes in the HOXA, HOXB, HOXC, and HOXD clusters (T-box and homeobox developmental genes). The data reveal strong correlations between athero hypermethylated DMRs and regions of enhancer chromatin in aorta, which are discussed in the linked research article "Atherosclerosis-associated differentially methylated regions can reflect the disease phenotype and are often at enhancers" (M. Lacey et al., 2019).

Atherosclerosis-associated differentially methylated regions can reflect the disease phenotype and are often at enhancers.

BACKGROUND AND AIMS: Atherosclerosis is a widespread and complicated disease involving phenotypic modulation and transdifferentiation of vascular smooth muscle cells (SMCs), the predominant cells in aorta, as well as changes in endothelial cells and infiltrating monocytes. Alterations in DNA methylation are likely to play central roles in these phenotypic changes, just as they do in normal differentiation and cancer. METHODS: We examined genome-wide DNA methylation changes in atherosclerotic aorta using more stringent criteria for differentially methylated regions (DMRs) than in previous studies and compared these DMRs to tissue-specific epigenetic features. RESULTS: We found that disease-linked hypermethylated DMRs account for 85% of the total atherosclerosis-associated DMRs and often overlap aorta-associated enhancer chromatin. These hypermethylated DMRs were associated with functionally different sets of genes compared to atherosclerosis-linked hypomethylated DMRs. The extent and nature of the DMRs could not be explained as direct effects of monocyte/macrophage infiltration. Among the known atherosclerosis- and contractile SMC-related genes that exhibited hypermethylated DMRs at aorta enhancer chromatin were ACTA2 (aorta α2 smooth muscle actin), ELN (elastin), MYOCD (myocardin), C9orf3 (miR-23b and miR-27b host gene), and MYH11 (smooth muscle myosin). Our analyses also suggest a role in atherosclerosis for developmental transcription factor genes having little or no previous association with atherosclerosis, such as NR2F2 (COUP-TFII) and TBX18. CONCLUSIONS: We provide evidence for atherosclerosis-linked DNA methylation changes in aorta SMCs that might help minimize or reverse the standard contractile character of many of these cells by down-modulating aorta SMC-related enhancers, thereby facilitating pro-atherosclerotic phenotypic changes in many SMCs.

Developmentally linked human DNA hypermethylation is associated with down-modulation, repression, and upregulation of transcription.

DNA methylation can affect tissue-specific gene transcription in ways that are difficult to discern from studies focused on genome-wide analyses of differentially methylated regions (DMRs). To elucidate the variety of associations between differentiation-related DNA hypermethylation and transcription, we used available epigenomic and transcriptomic profiles from 38 human cell/tissue types to focus on such relationships in 94 genes linked to hypermethylated DMRs in myoblasts (Mb). For 19 of the genes, promoter-region hypermethylation in Mb (and often a few heterologous cell types) was associated with gene repression but, importantly, DNA hypermethylation was absent in many other repressed samples. In another 24 genes, DNA hypermethylation overlapped cryptic enhancers or super-enhancers and correlated with down-modulated, but not silenced, gene expression. However, such methylation was absent, surprisingly, in both non-expressing samples and highly expressing samples. This suggests that some genes need DMR hypermethylation to help repress cryptic enhancer chromatin only when they are actively transcribed. For another 11 genes, we found an association between intergenic hypermethylated DMRs and positive expression of the gene in Mb. DNA hypermethylation/transcription correlations similar to those of Mb were evident sometimes in diverse tissues, such as aorta and brain. Our findings have implications for the possible involvement of methylated DNA in Duchenne's muscular dystrophy, congenital heart malformations, and cancer. This epigenomic analysis suggests that DNA methylation is not simply the inevitable consequence of changes in gene expression but, instead, is often an active agent for fine-tuning transcription in association with development.

Tissue-specific epigenetics in gene neighborhoods: myogenic transcription factor genes.

Myogenic regulatory factor (MRF) genes, MYOD1, MYOG, MYF6 and MYF5, are critical for the skeletal muscle lineage. Here, we used various epigenome profiles from human myoblasts (Mb), myotubes (Mt), muscle and diverse non-muscle samples to elucidate the involvement of multigene neighborhoods in the regulation of MRF genes. We found more far-distal enhancer chromatin associated with MRF genes in Mb and Mt than previously reported from studies in mice. For the MYF5/MYF6 gene-pair, regions of Mb-associated enhancer chromatin were located throughout the adjacent 236-kb PTPRQ gene even though Mb expressed negligible amounts of PTPRQ mRNA. Some enhancer chromatin regions inside PTPRQ in Mb were also seen in PTPRQ mRNA-expressing non-myogenic cells. This suggests dual-purpose PTPRQ enhancers that upregulate expression of PTPRQ in non-myogenic cells and MYF5/MYF6 in myogenic cells. In contrast, the myogenic enhancer chromatin regions distal to MYOD1 were intergenic and up to 19 kb long. Two of them contain small, known MYOD1 enhancers, and one displayed an unusually high level of 5-hydroxymethylcytosine in a quantitative DNA hydroxymethylation assay. Unexpectedly, three regions of MYOD1-distal enhancer chromatin in Mb and Mt overlapped enhancer chromatin in umbilical vein endothelial cells, which might upregulate a distant gene (PIK3C2A). Lastly, genes surrounding MYOG were preferentially transcribed in Mt, like MYOG itself, and exhibited nearby myogenic enhancer chromatin. These neighboring chromatin regions may be enhancers acting in concert to regulate myogenic expression of multiple adjacent genes. Our findings reveal the very different and complex organization of gene neighborhoods containing closely related transcription factor genes.

Myogenic differential methylation: diverse associations with chromatin structure.

Employing a new algorithm for identifying differentially methylated regions (DMRs) from reduced representation bisulfite sequencing profiles, we identified 1972 hypermethylated and 3250 hypomethylated myogenic DMRs in a comparison of myoblasts (Mb) and myotubes (Mt) with 16 types of nonmuscle cell cultures. DMRs co-localized with a variety of chromatin structures, as deduced from ENCODE whole-genome profiles. Myogenic hypomethylation was highly associated with both weak and strong enhancer-type chromatin, while hypermethylation was infrequently associated with enhancer-type chromatin. Both myogenic hypermethylation and hypomethylation often overlapped weak transcription-type chromatin and Polycomb-repressed-type chromatin. For representative genes, we illustrate relationships between DNA methylation, the local chromatin state, DNaseI hypersensitivity, and gene expression. For example, MARVELD2 exhibited myogenic hypermethylation in transcription-type chromatin that overlapped a silenced promoter in Mb and Mt while TEAD4 had myogenic hypomethylation in intronic subregions displaying enhancer-type or transcription-type chromatin in these cells. For LSP1, alternative promoter usage and active promoter-type chromatin were linked to highly specific myogenic or lymphogenic hypomethylated DMRs. Lastly, despite its myogenesis-associated expression, TBX15 had multiple hypermethylated myogenic DMRs framing its promoter region. This could help explain why TBX15 was previously reported to be underexpressed and, unexpectedly, its promoter undermethylated in placentas exhibiting vascular intrauterine growth restriction.

RNA CoMPASS: a dual approach for pathogen and host transcriptome analysis of RNA-seq datasets.

High-throughput RNA sequencing (RNA-seq) has become an instrumental assay for the analysis of multiple aspects of an organism's transcriptome. Further, the analysis of a biological specimen's associated microbiome can also be performed using RNA-seq data and this application is gaining interest in the scientific community. There are many existing bioinformatics tools designed for analysis and visualization of transcriptome data. Despite the availability of an array of next generation sequencing (NGS) analysis tools, the analysis of RNA-seq data sets poses a challenge for many biomedical researchers who are not familiar with command-line tools. Here we present RNA CoMPASS, a comprehensive RNA-seq analysis pipeline for the simultaneous analysis of transcriptomes and metatranscriptomes from diverse biological specimens. RNA CoMPASS leverages existing tools and parallel computing technology to facilitate the analysis of even very large datasets. RNA CoMPASS has a web-based graphical user interface with intrinsic queuing to control a distributed computational pipeline. RNA CoMPASS was evaluated by analyzing RNA-seq data sets from 45 B-cell samples. Twenty-two of these samples were derived from lymphoblastoid cell lines (LCLs) generated by the infection of naïve B-cells with the Epstein Barr virus (EBV), while another 23 samples were derived from Burkitt's lymphomas (BL), some of which arose in part through infection with EBV. Appropriately, RNA CoMPASS identified EBV in all LCLs and in a fraction of the BLs. Cluster analysis of the human transcriptome component of the RNA CoMPASS output clearly separated the BLs (which have a germinal center-like phenotype) from the LCLs (which have a blast-like phenotype) with evidence of activated MYC signaling and lower interferon and NF-kB signaling in the BLs. Together, this analysis illustrates the utility of RNA CoMPASS in the simultaneous analysis of transcriptome and metatranscriptome data. RNA CoMPASS is freely available at http://rnacompass.sourceforge.net/.

Modeling, simulation and analysis of methylation profiles from reduced representation bisulfite sequencing experiments.

The ENCODE project has funded the generation of a diverse collection of methylation profiles using reduced representation bisulfite sequencing (RRBS) technology, enabling the analysis of epigenetic variation on a genomic scale at single-site resolution. A standard application of RRBS experiments is in the location of differentially methylated regions (DMRs) between two sets of samples. Despite numerous publications reporting DMRs identified from RRBS datasets, there have been no formal analyses of the effects of experimental and biological factors on the performance of existing or newly developed analytical methods. These factors include variable read coverage, differing group sample sizes across genomic regions, uneven spacing between CpG dinucleotide sites, and correlation in methylation levels among sites in close proximity. To better understand the interplay among technical and biological variables in the analysis of RRBS methylation profiles, we have developed an algorithm for the generation of experimentally realistic RRBS datasets. Applying insights derived from our simulation studies, we present a novel procedure that can identify DMRs spanning as few as three CpG sites with both high sensitivity and specificity. Using RRBS data from muscle vs. non-muscle cell cultures as an example, we demonstrate that our method reveals many more DMRs that are likely to be of biological significance than previous methods.

DNA methylation and differentiation: HOX genes in muscle cells.

BACKGROUND: Tight regulation of homeobox genes is essential for vertebrate development. In a study of genome-wide differential methylation, we recently found that homeobox genes, including those in the HOX gene clusters, were highly overrepresented among the genes with hypermethylation in the skeletal muscle lineage. Methylation was analyzed by reduced representation bisulfite sequencing (RRBS) of postnatal myoblasts, myotubes and adult skeletal muscle tissue and 30 types of non-muscle-cell cultures or tissues. RESULTS: In this study, we found that myogenic hypermethylation was present in specific subregions of all four HOX gene clusters and was associated with various chromatin epigenetic features. Although the 3' half of the HOXD cluster was silenced and enriched in polycomb repression-associated H3 lysine 27 trimethylation in most examined cell types, including myoblasts and myotubes, myogenic samples were unusual in also displaying much DNA methylation in this region. In contrast, both HOXA and HOXC clusters displayed myogenic hypermethylation bordering a central region containing many genes preferentially expressed in myogenic progenitor cells and consisting largely of chromatin with modifications typical of promoters and enhancers in these cells. A particularly interesting example of myogenic hypermethylation was HOTAIR, a HOXC noncoding RNA gene, which can silence HOXD genes in trans via recruitment of polycomb proteins. In myogenic progenitor cells, the preferential expression of HOTAIR was associated with hypermethylation immediately downstream of the gene. Other HOX gene regions also displayed myogenic DNA hypermethylation despite being moderately expressed in myogenic cells. Analysis of representative myogenic hypermethylated sites for 5-hydroxymethylcytosine revealed little or none of this base, except for an intragenic site in HOXB5 which was specifically enriched in this base in skeletal muscle tissue, whereas myoblasts had predominantly 5-methylcytosine at the same CpG site. CONCLUSIONS: Our results suggest that myogenic hypermethylation of HOX genes helps fine-tune HOX sense and antisense gene expression through effects on 5' promoters, intragenic and intergenic enhancers and internal promoters. Myogenic hypermethylation might also affect the relative abundance of different RNA isoforms, facilitate transcription termination, help stop the spread of activation-associated chromatin domains and stabilize repressive chromatin structures.

Differences in gastric carcinoma microenvironment stratify according to EBV infection intensity: implications for possible immune adjuvant therapy.

Epstein-Barr virus (EBV) is associated with roughly 10% of gastric carcinomas worldwide (EBVaGC). Although previous investigations provide a strong link between EBV and gastric carcinomas, these studies were performed using selected EBV gene probes. Using a cohort of gastric carcinoma RNA-seq data sets from The Cancer Genome Atlas (TCGA), we performed a quantitative and global assessment of EBV gene expression in gastric carcinomas and assessed EBV associated cellular pathway alterations. EBV transcripts were detected in 17% of samples but these samples varied significantly in EBV coverage depth. In four samples with the highest EBV coverage (hiEBVaGC - high EBV associated gastric carcinoma), transcripts from the BamHI A region comprised the majority of EBV reads. Expression of LMP2, and to a lesser extent, LMP1 were also observed as was evidence of abortive lytic replication. Analysis of cellular gene expression indicated significant immune cell infiltration and a predominant IFNG response in samples expressing high levels of EBV transcripts relative to samples expressing low or no EBV transcripts. Despite the apparent immune cell infiltration, high levels of the cytotoxic T-cell (CTL) and natural killer (NK) cell inhibitor, IDO1, was observed in the hiEBVaGCs samples suggesting an active tolerance inducing pathway in this subgroup. These results were confirmed in a separate cohort of 21 Vietnamese gastric carcinoma samples using qRT-PCR and on tissue samples using in situ hybridization and immunohistochemistry. Lastly, a panel of tumor suppressors and candidate oncogenes were expressed at lower levels in hiEBVaGC versus EBV-low and EBV-negative gastric cancers suggesting the direct regulation of tumor pathways by EBV.

Early de novo DNA methylation and prolonged demethylation in the muscle lineage.

Myogenic cell cultures derived from muscle biopsies are excellent models for human cell differentiation. We report the first comprehensive analysis of myogenesis-specific DNA hyper- and hypo-methylation throughout the genome for human muscle progenitor cells (both myoblasts and myotubes) and skeletal muscle tissue vs. 30 non-muscle samples using reduced representation bisulfite sequencing. We also focused on four genes with extensive hyper- or hypo-methylation in the muscle lineage (PAX3, TBX1, MYH7B/MIR499 and OBSCN) to compare DNA methylation, DNaseI hypersensitivity, histone modification, and CTCF binding profiles. We found that myogenic hypermethylation was strongly associated with homeobox or T-box genes and muscle hypomethylation with contractile fiber genes. Nonetheless, there was no simple relationship between differential gene expression and myogenic differential methylation, rather only for subsets of these genes, such as contractile fiber genes. Skeletal muscle retained ~30% of the hypomethylated sites but only ~3% of hypermethylated sites seen in myogenic progenitor cells. By enzymatic assays, skeletal muscle was 2-fold enriched globally in genomic 5-hydroxymethylcytosine (5-hmC) vs. myoblasts or myotubes and was the only sample type enriched in 5-hmC at tested myogenic hypermethylated sites in PAX3/CCDC140 andTBX1. TET1 and TET2 RNAs, which are involved in generation of 5-hmC and DNA demethylation, were strongly upregulated in myoblasts and myotubes. Our findings implicate de novo methylation predominantly before the myoblast stage and demethylation before and after the myotube stage in control of transcription and co-transcriptional RNA processing. They also suggest that, in muscle, TET1 or TET2 are involved in active demethylation and in formation of stable 5-hmC residues.

Whole-genome sequencing of the Akata and Mutu Epstein-Barr virus strains.

Using a simple viral genome enrichment approach, we report the de novo assembly of the Akata and Mutu Epstein-Barr virus (EBV) genomes from a single lane of next-generation sequencing (NGS) reads. The Akata and Mutu viral genomes are type I EBV strains of approximately 171 kb in length. Evidence for genome heterogeneity was found for the Akata but not for the Mutu strain. A comparative analysis of Akata with another four completely sequenced EBV strains, B95-8/Raji, AG876, Mutu, and GD1, demonstrated that the Akata strain is most closely related to the GD1 strain and exhibits the greatest divergence from the type II strain, AG876. A global comparison of latent and lytic gene sequences showed that the four latency genes, EBNA2, EBNA3A, EBNA3B, and EBNA3C, are uniquely defining of type I and type II strain differences. Within type I strains, LMP1, the latency gene, is among the most divergent of all EBV genes, with three insertion or deletion loci in its CTAR2 and CTAR3 signaling domains. Analysis of the BHLF1 and LF3 genes showed that the reading frames identified in the B95-8/Raji genome are not conserved in Akata (or Mutu, for BHLF1), suggesting a primarily non-protein-coding function in EBV's life cycle. The Akata and Mutu viral-genome sequences should be a useful resource for homology-based functional prediction and for molecular studies, such as PCR, RNA-seq, recombineering, and transcriptome studies. As an illustration, we identified novel RNA-editing events in ebv-miR-BART6 antisense transcripts using the Akata and Mutu reference genomes.

Gene expression during normal and FSHD myogenesis.

BACKGROUND: Facioscapulohumeral muscular dystrophy (FSHD) is a dominant disease linked to contraction of an array of tandem 3.3-kb repeats (D4Z4) at 4q35. Within each repeat unit is a gene, DUX4, that can encode a protein containing two homeodomains. A DUX4 transcript derived from the last repeat unit in a contracted array is associated with pathogenesis but it is unclear how. METHODS: Using exon-based microarrays, the expression profiles of myogenic precursor cells were determined. Both undifferentiated myoblasts and myoblasts differentiated to myotubes derived from FSHD patients and controls were studied after immunocytochemical verification of the quality of the cultures. To further our understanding of FSHD and normal myogenesis, the expression profiles obtained were compared to those of 19 non-muscle cell types analyzed by identical methods. RESULTS: Many of the ~17,000 examined genes were differentially expressed (>2-fold, p<0.01) in control myoblasts or myotubes vs. non-muscle cells (2185 and 3006, respectively) or in FSHD vs. control myoblasts or myotubes (295 and 797, respectively). Surprisingly, despite the morphologically normal differentiation of FSHD myoblasts to myotubes, most of the disease-related dysregulation was seen as dampening of normal myogenesis-specific expression changes, including in genes for muscle structure, mitochondrial function, stress responses, and signal transduction. Other classes of genes, including those encoding extracellular matrix or pro-inflammatory proteins, were upregulated in FSHD myogenic cells independent of an inverse myogenesis association. Importantly, the disease-linked DUX4 RNA isoform was detected by RT-PCR in FSHD myoblast and myotube preparations only at extremely low levels. Unique insights into myogenesis-specific gene expression were also obtained. For example, all four Argonaute genes involved in RNA-silencing were significantly upregulated during normal (but not FSHD) myogenesis relative to non-muscle cell types. CONCLUSIONS: DUX4's pathogenic effect in FSHD may occur transiently at or before the stage of myoblast formation to establish a cascade of gene dysregulation. This contrasts with the current emphasis on toxic effects of experimentally upregulated DUX4 expression at the myoblast or myotube stages. Our model could explain why DUX4's inappropriate expression was barely detectable in myoblasts and myotubes but nonetheless linked to FSHD.

Duration learning for analysis of nanopore ionic current blockades.

BACKGROUND: Ionic current blockade signal processing, for use in nanopore detection, offers a promising new way to analyze single molecule properties, with potential implications for DNA sequencing. The alpha-Hemolysin transmembrane channel interacts with a translocating molecule in a nontrivial way, frequently evidenced by a complex ionic flow blockade pattern. Typically, recorded current blockade signals have several levels of blockade, with various durations, all obeying a fixed statistical profile for a given molecule. Hidden Markov Model (HMM) based duration learning experiments on artificial two-level Gaussian blockade signals helped us to identify proper modeling framework. We then apply our framework to the real multi-level DNA hairpin blockade signal. RESULTS: The identified upper level blockade state is observed with durations that are geometrically distributed (consistent with an a physical decay process for remaining in any given state). We show that mixture of convolution chains of geometrically distributed states is better for presenting multimodal long-tailed duration phenomena. Based on learned HMM profiles we are able to classify 9 base-pair DNA hairpins with accuracy up to 99.5% on signals from same-day experiments. CONCLUSION: We have demonstrated several implementations for de novo estimation of duration distribution probability density function with HMM framework and applied our model topology to the real data. The proposed design could be handy in molecular analysis based on nanopore current blockade signal.

A novel, fast, HMM-with-Duration implementation - for application with a new, pattern recognition informed, nanopore detector.

BACKGROUND: Hidden Markov Models (HMMs) provide an excellent means for structure identification and feature extraction on stochastic sequential data. An HMM-with-Duration (HMMwD) is an HMM that can also exactly model the hidden-label length (recurrence) distributions - while the regular HMM will impose a best-fit geometric distribution in its modeling/representation. RESULTS: A Novel, Fast, HMM-with-Duration (HMMwD) Implementation is presented, and experimental results are shown that demonstrate its performance on two-state synthetic data designed to model Nanopore Detector Data. The HMMwD experimental results are compared to (i) the ideal model and to (ii) the conventional HMM. Its accuracy is clearly an improvement over the standard HMM, and matches that of the ideal solution in many cases where the standard HMM does not. Computationally, the new HMMwD has all the speed advantages of the conventional (simpler) HMM implementation. In preliminary work shown here, HMM feature extraction is then used to establish the first pattern recognition-informed (PRI) sampling control of a Nanopore Detector Device (on a "live" data-stream). CONCLUSION: The improved accuracy of the new HMMwD implementation, at the same order of computational cost as the standard HMM, is an important augmentation for applications in gene structure identification and channel current analysis, especially PRI sampling control, for example, where speed is essential. The PRI experiment was designed to inherit the high accuracy of the well characterized and distinctive blockades of the DNA hairpin molecules used as controls (or blockade "test-probes"). For this test set, the accuracy inherited is 99.9%.

Cheminformatics methods for novel nanopore analysis of HIV DNA termini.

BACKGROUND: Channel current feature extraction methods, using Hidden Markov Models (HMMs) have been designed for tracking individual-molecule conformational changes. This information is derived from observation of changes in ionic channel current blockade "signal" upon that molecule's interaction with (and occlusion of) a single nanometer-scale channel in a "nanopore detector". In effect, a nanopore detector transduces single molecule events into channel current blockades. HMM analysis tools described are used to help systematically explore DNA dinucleotide flexibility, with particular focus on HIV's highly conserved (and highly flexible/reactive) viral DNA termini. One of the most critical stages in HIV's attack is the binding between viral DNA and the retroviral integrase, which is influenced by the dynamic-coupling induced high flexibility of a CA/TG dinucleotide positioned precisely two base-pairs from the blunt terminus of the duplex viral DNA. This suggests the study of a family of such CA/TG dinucleotide molecules via nanopore measurement and cheminformatics analysis. RESULTS: HMMs are used for level identification on the current blockades, HMM/EM with boosted variance emissions are used for level projection pre-processing, and time-domain FSAs are used to parse the level-projected waveform for kinetic information. The observed state kinetics of the DNA hairpins containing the CA/TG dinucleotide provides clear evidence for HIV's selection of a peculiarly flexible/interactive DNA terminus.