De novo deciphering three-dimensional chromatin interaction and topological domains by wavelet transformation of epigenetic profiles

Defining chromatin interaction frequencies and topological domains is a great challenge for the annotations of genome structures. Although the chromosome conformation capture (3C) and its derivative methods have been developed for exploring the global interactome, they are limited by high experimental complexity and costs. Here we describe a novel computational method, called CITD, for de novo prediction of the chromatin interaction map by integrating histone modification data. We used the public epigenomic data from human fibroblast IMR90 cell and embryonic stem cell (H1) to develop and test CITD, which can not only successfully reconstruct the chromatin interaction frequencies discovered by the Hi-C technology, but also provide additional novel details of chromosomal organizations. We predicted the chromatin interaction frequencies, topological domains and their states (e.g. active or repressive) for 98 additional cell types from Roadmap Epigenomics and ENCODE projects. A total of 131 protein-coding genes located near 78 preserved boundaries among 100 cell types are found to be significantly enriched in functional categories of the nucleosome organization and chromatin assembly. CITD and its predicted results can be used for complementing the topological domains derived from limited Hi-C data and facilitating the understanding of spatial principles underlying the chromosomal organization.

Zinc treatment increases the titre of 'Candidatus Liberibacter asiaticus' in huanglongbing-affected citrus plants while affecting the bacterial microbiomes

AIMS

Huanglongbing (HLB)-affected citrus often display zinc deficiency symptoms. In this study, supplemental zinc was applied to citrus to determine its effect on Candidatus Liberibacter asiaticus (Las) titre, HLB symptoms, and leaf microbiome.

METHODS AND RESULTS

HLB-affected citrus were treated with various amounts of zinc. The treatments promoted Las growth and affected microbiomes in citrus leaves. Phylochip(™) -based results indicated that 5475 of over 50 000 known Operational Taxonomic Units (OTUs) in 52 phyla were detected in the midribs of HLB-affected citrus, of which Proteobacteria was the most abundant, followed by Firmicutes and Actinobacteria. In comparison, the microbiomes of zinc-treated diseased plants had overall more OTUs with higher amounts of Proteobacteria, but decreased percentages of Firmicutes and Actinobacteria. In addition, more OTUs of siderophore-producing bacteria were present. Only zinc-sensitive Staphylococcaceae had higher OTU's in the diseased plants without zinc treatments.

CONCLUSIONS

Although HLB-affected citrus appear zinc deficient, zinc amendments increased the pathogen levels and shifted the microbiome.

SIGNIFICANCE AND IMPACT OF THE STUDY

HLB is currently the most devastating disease of citrus worldwide. Zinc is often applied to HLB-affected citrus due to zinc deficiency symptoms. This study provided new insights into the potential effects of zinc on HLB and the microbial ecology of citrus.

Neurobehavioural abnormalities induced by repeated exposure of neonatal rats to sevoflurane can be aggravated by social isolation and enrichment deprivation initiated after exposure to the anaesthetic

BACKGROUND

We tested the hypothesis that developmental effects of repeated neonatal exposure to sevoflurane in rats are exacerbated by stressful experiences received later in life.

METHODS

Sprague-Dawley male rats received sequential exposures to 3% sevoflurane for two h on postnatal days (P) six, seven, and eight. After weaning at P21, rats were housed either in pairs in an enriched environment (EE) or singly in an enrichment-deprived environment (an adverse environment, AE). The hippocampal concentrations of brain-derived neurotrophic factor (BDNF), and synaptic markers were assessed at P8 and P53. The dentate gyrus neural progenitor proliferation was evaluated at P11 and P53 after administration of bromodeoyuridine (BrdU) at P8 to P10 and at P22 to P27, respectively. Neurobehavioural evaluations were performed at P49 to P53.

RESULTS

Repeated sevoflurane exposure acutely reduced concentrations of BDNF, synaptic markers and neural progenitor proliferation. The sevoflurane group housed in the AE conditions (sevoflurane+AE) had decreased concentrations of BDNF and synaptic markers, and survival of new granule cells and impaired cognitive function compared with the control+AE, control+EE, and sevoflurane+EE groups. The neurobehavioural parameters in the sevoflurane+EE and control+EE groups were similar.

CONCLUSIONS

Neurocognitive abnormalities induced by repeated neonatal exposure to sevoflurane can be aggravated by stressful conditions such as social isolation and enrichment deprivation.

Validation of SCT Methylation as a Hallmark Biomarker for Lung Cancers

INTRODUCTION

The human secretin gene (SCT) encodes secretin, a hormone with limited tissue distribution. Analysis of the 450k methylation array data in The Cancer Genome Atlas (TCGA) indicated that the SCT promoter region is differentially hypermethylated in lung cancer. Our purpose was to validate SCT methylation as a potential biomarker for lung cancer.

METHODS

We analyzed data from TCGA and developed and applied SCT-specific bisulfite DNA sequencing and quantitative methylation-specific polymerase chain reaction assays.

RESULTS

The analyses of TCGA 450K data for 801 samples showed that SCT hypermethylation has an area under the curve (AUC) value greater than 0.98 that can be used to distinguish lung adenocarcinomas or squamous cell carcinomas from nonmalignant lung tissue. Bisulfite sequencing of lung cancer cell lines and normal blood cells allowed us to confirm that SCT methylation is highly discriminative. By applying a quantitative methylation-specific polymerase chain reaction assay, we found that SCT hypermethylation is frequently detected in all major subtypes of malignant non-small cell lung cancer (AUC = 0.92, n = 108) and small cell lung cancer (AUC = 0.93, n = 40) but is less frequent in lung carcinoids (AUC = 0.54, n = 20). SCT hypermethylation appeared in samples of lung carcinoma in situ during multistage pathogenesis and increased in invasive samples. Further analyses of TCGA 450k data showed that SCT hypermethylation is highly discriminative in most other types of malignant tumors but less frequent in low-grade malignant tumors. The only normal tissue with a high level of methylation was the placenta.

CONCLUSIONS

Our findings demonstrated that SCT methylation is a highly discriminative biomarker for lung and other malignant tumors, is less frequent in low-grade malignant tumors (including lung carcinoids), and appears at the carcinoma in situ stage.

3CPET: finding co-factor complexes from ChIA-PET data using a hierarchical Dirichlet process

Various efforts have been made to elucidate the cooperating proteins involved in maintaining chromatin interactions; however, many are still unknown. Here, we present 3CPET, a tool based on a non-parametric Bayesian approach, to infer the set of the most probable protein complexes involved in maintaining chromatin interactions and the regions that they may control, making it a valuable downstream analysis tool in chromatin conformation studies. 3CPET does so by combining data from ChIA-PET, transcription factor binding sites, and protein interactions. 3CPET results show biologically significant and accurate predictions when validated against experimental and simulation data.

Multi-modulus algorithm based on global artificial fish swarm intelligent optimization of DNA encoding sequences

Aimed to address the defects of the large mean square error (MSE), and the slow convergence speed in equalizing the multi-modulus signals of the constant modulus algorithm (CMA), a multi-modulus algorithm (MMA) based on global artificial fish swarm (GAFS) intelligent optimization of DNA encoding sequences (GAFS-DNA-MMA) was proposed. To improve the convergence rate and reduce the MSE, this proposed algorithm adopted an encoding method based on DNA nucleotide chains to provide a possible solution to the problem. Furthermore, the GAFS algorithm, with its fast convergence and global search ability, was used to find the best sequence. The real and imaginary parts of the initial optimal weight vector of MMA were obtained through DNA coding of the best sequence. The simulation results show that the proposed algorithm has a faster convergence speed and smaller MSE in comparison with the CMA, the MMA, and the AFS-DNA-MMA.

Reconfigurable hybrid interface for molecular marker diagnostics and in-situ reporting

Combinations of molecular signals such as transcription factors and microRNAs in cells are a reliable indicator of multi-gene disorders. A system capable of detecting these conditions in-situ may be used as a tool for diagnosis and monitoring of disease. Here, we engineer genetic circuits that sense endogenous levels of the androgen receptor (AR), the glucocorticoid receptor (GR), and the microRNA hsa-miR-21 (miR-21) in cervical cancer cells (HeLa). Furthermore, using the mediator molecule human chorionic gonadotropin (hCG), we interface the intracellular information to enzyme-linked immunosorbent assay (ELISA) test strips. We demonstrate that this hybrid genetic circuit and test-strip interface can accommodate combinatorial, low-cost, and in-situ reporting, a versatile profiling tool.

Fast dimension reduction and integrative clustering of multi-omics data using low-rank approximation: application to cancer molecular classification

Background

One major goal of large-scale cancer omics study is to identify molecular subtypes for more accurate cancer diagnoses and treatments. To deal with high-dimensional cancer multi-omics data, a promising strategy is to find an effective low-dimensional subspace of the original data and then cluster cancer samples in the reduced subspace. However, due to data-type diversity and big data volume, few methods can integrative and efficiently find the principal low-dimensional manifold of the high-dimensional cancer multi-omics data.

Results

In this study, we proposed a novel low-rank approximation based integrative probabilistic model to fast find the shared principal subspace across multiple data types: the convexity of the low-rank regularized likelihood function of the probabilistic model ensures efficient and stable model fitting. Candidate molecular subtypes can be identified by unsupervised clustering hundreds of cancer samples in the reduced low-dimensional subspace. On testing datasets, our method LRAcluster (low-rank approximation based multi-omics data clustering) runs much faster with better clustering performances than the existing method. Then, we applied LRAcluster on large-scale cancer multi-omics data from TCGA. The pan-cancer analysis results show that the cancers of different tissue origins are generally grouped as independent clusters, except squamous-like carcinomas. While the single cancer type analysis suggests that the omics data have different subtyping abilities for different cancer types.

Conclusions

LRAcluster is a very useful method for fast dimension reduction and unsupervised clustering of large-scale multi-omics data. LRAcluster is implemented in R and freely available via http://bioinfo.au.tsinghua.edu.cn/software/lracluster/.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2223-8) contains supplementary material, which is available to authorized users.

NONCODE 2016: an informative and valuable data source of long non-coding RNAs

NONCODE (http://www.bioinfo.org/noncode/) is an interactive database that aims to present the most complete collection and annotation of non-coding RNAs, especially long non-coding RNAs (lncRNAs). The recently reduced cost of RNA sequencing has produced an explosion of newly identified data. Revolutionary third-generation sequencing methods have also contributed to more accurate annotations. Accumulative experimental data also provides more comprehensive knowledge of lncRNA functions. In this update, NONCODE has added six new species, bringing the total to 16 species altogether. The lncRNAs in NONCODE have increased from 210 831 to 527,336. For human and mouse, the lncRNA numbers are 167,150 and 130,558, respectively. NONCODE 2016 has also introduced three important new features: (i) conservation annotation; (ii) the relationships between lncRNAs and diseases; and (iii) an interface to choose high-quality datasets through predicted scores, literature support and long-read sequencing method support. NONCODE is also accessible through http://www.noncode.org/.

Histone Deacetylases Positively Regulate Transcription through the Elongation Machinery

Transcription elongation regulates the expression of many genes, including oncogenes. Histone deacetylase (HDAC) inhibitors (HDACIs) block elongation, suggesting that HDACs are involved in gene activation. To understand this, we analyzed nascent transcription and elongation factor binding genome-wide after perturbation of elongation with small molecule inhibitors. We found that HDACI-mediated repression requires heat shock protein 90 (HSP90) activity. HDACIs promote the association of RNA polymerase II (RNAP2) and negative elongation factor (NELF), a complex stabilized by HSP90, at the same genomic sites. Additionally, HDACIs redistribute bromodomain-containing protein 4 (BRD4), a key elongation factor involved in enhancer activity. BRD4 binds to newly acetylated sites, and its occupancy at promoters and enhancers is reduced. Furthermore, HDACIs reduce enhancer activity, as measured by enhancer RNA production. Therefore, HDACs are required for limiting acetylation in gene bodies and intergenic regions. This facilitates the binding of elongation factors to properly acetylated promoters and enhancers for efficient elongation.

Functional and association studies of the cholesteryl ester transfer protein (CETP) gene in a Wannan Black pig model

Some polymorphisms of the human CETP gene are causally and significantly associated with serum lipids levels; however, the information regarding this gene in pigs is sparse. To evaluate the effects of CETP on blood lipid traits and fat deposition in pig, porcine CETP tissue expression patterns were observed by quantitative real-time polymerase chain reaction (qPCR) first. High expression was detected in liver, spleen, gluteus medius (GM) muscle and backfat. A de novo polymorphism (AF333037:g.795C>T) in the intron 1 region of porcine CETP was identified. This polymorphism was further genotyped by direct sequencing of the PCR products of 390 Wannan Black pigs, a Chinese native breed population. Association analyses at 45 and 300 days of age revealed highly significant associations between CETP genotypes and serum lipid traits. Furthermore, this polymorphism was proved to be associated with differences in liver CETP mRNA levels: pigs at 300 days of age with the TT genotype had higher levels than did those with other genotypes (P = 0.021). Additionally, analysis at 300 days of age showed that GM CETP mRNA expression correlated positively with serum lipids levels as well as with carcass backfat thickness and intramuscular fat content in GM. These results indicate that CETP is involved in serum, adipose and muscle lipid metabolism in pigs. The mechanisms underlying such relationships and their functional implications are worthy of further research.

Genome-wide analysis of the response to nitric oxide in uropathogenic Escherichia coli CFT073

Uropathogenic Escherchia coli (UPEC) is the causative agent of urinary tract infections. Nitric oxide (NO) is a toxic water-soluble gas that is encountered by UPEC in the urinary tract. Therefore, UPEC probably requires mechanisms to detoxify NO in the host environment. Thus far, flavohaemoglobin (Hmp), an NO denitrosylase, is the only demonstrated NO detoxification system in UPEC. Here we show that, in E. coli strain CFT073, the NADH-dependent NO reductase flavorubredoxin (FlRd) also plays a major role in NO scavenging. We generated a mutant that lacks all known and candidate NO detoxification pathways (Hmp, FlRd and the respiratory nitrite reductase, NrfA). When grown and assayed anaerobically, this mutant expresses an NO-inducible NO scavenging activity, pointing to the existence of a novel detoxification mechanism. Expression of this activity is inducible by both NO and nitrate, and the enzyme is membrane-associated. Genome-wide transcriptional profiling of UPEC grown under anaerobic conditions in the presence of nitrate (as a source of NO) highlighted various aspects of the response of the pathogen to nitrate and NO. Several virulence-associated genes are upregulated, suggesting that host-derived NO is a potential regulator of UPEC virulence. Chromatin immunoprecipitation and sequencing was used to evaluate the NsrR regulon in CFT073. We identified 49 NsrR binding sites in promoter regions in the CFT073 genome, 29 of which were not previously identified in E. coli K-12. NsrR may regulate some CFT073 genes that do not have homologues in E. coli K-12.

SCT Promoter Methylation is a Highly Discriminative Biomarker for Lung and Many Other Cancers

Aberrant DNA methylation has long been implicated in cancers. In this work we present a highly discriminative DNA methylation biomarker for non-small cell lung cancers and fourteen other cancers. Based on 69 NSCLC cell lines and 257 cancer-free lung tissues we identified a CpG island in SCT gene promoter which was verified by qMSP experiment in 15 NSCLC cell lines and 3 immortalized human respiratory epithelium cells. In addition, we found that SCT promoter was methylated in 23 cancer cell lines involving >10 cancer types profiled by ENCODE. We found that SCT promoter is hyper-methylated in primary tumors from TCGA lung cancer cohort. Additionally, we found that SCT promoter is methylated at high frequencies in fifteen malignancies and is not methylated in~1000 non-cancerous tissues across >30 organ types. Our study indicates that SCT promoter methylation is a highly discriminative biomarker for lung and many other cancers.

CRISPR Inversion of CTCF Sites Alters Genome Topology and Enhancer/Promoter Function

CTCF and the associated cohesin complex play a central role in insulator function and higher-order chromatin organization of mammalian genomes. Recent studies identified a correlation between the orientation of CTCF-binding sites (CBSs) and chromatin loops. To test the functional significance of this observation, we combined CRISPR/Cas9-based genomic-DNA-fragment editing with chromosome-conformation-capture experiments to show that the location and relative orientations of CBSs determine the specificity of long-range chromatin looping in mammalian genomes, using protocadherin (Pcdh) and β-globin as model genes. Inversion of CBS elements within the Pcdh enhancer reconfigures the topology of chromatin loops between the distal enhancer and target promoters and alters gene-expression patterns. Thus, although enhancers can function in an orientation-independent manner in reporter assays, in the native chromosome context, the orientation of at least some enhancers carrying CBSs can determine both the architecture of topological chromatin domains and enhancer/promoter specificity. These findings reveal how 3D chromosome architecture can be encoded by linear genome sequences.

MICC: an R package for identifying chromatin interactions from ChIA-PET data

ChIA-PET is rapidly emerging as an important experimental approach to detect chromatin long-range interactions at high resolution. Here, we present Model based Interaction Calling from ChIA-PET data (MICC), an easy-to-use R package to detect chromatin interactions from ChIA-PET sequencing data. By applying a Bayesian mixture model to systematically remove random ligation and random collision noise, MICC could identify chromatin interactions with a significantly higher sensitivity than existing methods at the same false discovery rate.

AVAILABILITY AND IMPLEMENTATION

http://bioinfo.au.tsinghua.edu.cn/member/xwwang/MICCusage

CONTACT

michael.zhang@utdallas.edu or xwwang@tsinghua.edu.cn.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Differential connectivity of splicing activators and repressors to the human spliceosome

Background

During spliceosome assembly, protein-protein interactions (PPI) are sequentially formed and disrupted to accommodate the spatial requirements of pre-mRNA substrate recognition and catalysis. Splicing activators and repressors, such as SR proteins and hnRNPs, modulate spliceosome assembly and regulate alternative splicing. However, it remains unclear how they differentially interact with the core spliceosome to perform their functions.

Results

Here, we investigate the protein connectivity of SR and hnRNP proteins to the core spliceosome using probabilistic network reconstruction based on the integration of interactome and gene expression data. We validate our model by immunoprecipitation and mass spectrometry of the prototypical splicing factors SRSF1 and hnRNPA1. Network analysis reveals that a factor’s properties as an activator or repressor can be predicted from its overall connectivity to the rest of the spliceosome. In addition, we discover and experimentally validate PPIs between the oncoprotein SRSF1 and members of the anti-tumor drug target SF3 complex. Our findings suggest that activators promote the formation of PPIs between spliceosomal sub-complexes, whereas repressors mostly operate through protein-RNA interactions.

Conclusions

This study demonstrates that combining in-silico modeling with biochemistry can significantly advance the understanding of structure and function relationships in the human spliceosome.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0682-5) contains supplementary material, which is available to authorized users.

Molecular cloning, characterization and expression analysis of a Doublesex gene from Daphnia carinata (Crustacea: Cladocera) during different reproductive stages

To better understand the reproductive transformation mechanism of Daphnia carinata, a Doublesex (Dsx) gene was cloned based on rapid amplification of cDNA ends (RACE), and was designated DapcaDsx2. Next, we compared similarities and assumed homology based on deduced amino acid sequences. It showed 97.52, 87.94, and 85.11% identity to orthologous genes in D. magna, D. pulex, and D. galeata respectively. Phylogenetic analysis revealed that DapcaDsx2 clustered in the same class, and was evolutionarily more distant to sequences from other species. qRT-PCR showed that DapcaDsx2 was most abundantly expressed during sexual reproduction (P < 0.05). Using digoxigenin-labeled RNA probes, we studied DapcaDsx2 expression in parthenogenetic and sexual females by whole-mount in situ hybridization. The results revealed that DapcaDsx2 was mainly expressed in the second antennae and several sites of the ventral carapace, whereas higher expression levels were found in sexual than in parthenogenetic females. This suggests that the DapcaDsx2 gene is involved in switching modes of reproduction and in sexual differentiation.

Small C-terminal Domain Phosphatase 3 Dephosphorylates the Linker Sites of Receptor-regulated Smads (R-Smads) to Ensure Transforming Growth Factor β (TGFβ)-mediated Germ Layer Induction in Xenopus Embryos

Germ layer induction is one of the earliest events shortly after fertilization that initiates body formation of vertebrate embryos. In Xenopus, the maternally deposited transcriptional factor VegT promotes the expression of zygotic Nodal/Activin ligands that further form a morphogen gradient along the vegetal-animal axis and trigger the induction of the three germ layers. Here we found that SCP3 (small C-terminal domain phosphatase 3) is maternally expressed and vegetally enriched in Xenopus embryos and is essential for the timely induction of germ layers. SCP3 is required for the full activation of Nodal/Activin and bone morphogenetic protein signals and functions via dephosphorylation in the linker regions of receptor-regulated Smads. Consistently, the linker regions of receptor-regulated Smads are heavily phosphorylated in fertilized eggs, and this phosphorylation is gradually removed when embryos approach the midblastula transition. Knockdown of maternal SCP3 attenuates these dephosphorylation events and the activation of Nodal/Activin and bone morphogenetic protein signals after midblastula transition. This study thus suggested that the maternal SCP3 serves as a vegetally enriched, intrinsic factor to ensure a prepared status of Smads for their activation by the upcoming ligands during germ layer induction of Xenopus embryos.

Design and bioinformatics analysis of genome-wide CLIP experiments

The past decades have witnessed a surge of discoveries revealing RNA regulation as a central player in cellular processes. RNAs are regulated by RNA-binding proteins (RBPs) at all post-transcriptional stages, including splicing, transportation, stabilization and translation. Defects in the functions of these RBPs underlie a broad spectrum of human pathologies. Systematic identification of RBP functional targets is among the key biomedical research questions and provides a new direction for drug discovery. The advent of cross-linking immunoprecipitation coupled with high-throughput sequencing (genome-wide CLIP) technology has recently enabled the investigation of genome-wide RBP–RNA binding at single base-pair resolution. This technology has evolved through the development of three distinct versions: HITS-CLIP, PAR-CLIP and iCLIP. Meanwhile, numerous bioinformatics pipelines for handling the genome-wide CLIP data have also been developed. In this review, we discuss the genome-wide CLIP technology and focus on bioinformatics analysis. Specifically, we compare the strengths and weaknesses, as well as the scopes, of various bioinformatics tools. To assist readers in choosing optimal procedures for their analysis, we also review experimental design and procedures that affect bioinformatics analyses.

ChIP-Array 2: integrating multiple omics data to construct gene regulatory networks

Transcription factors (TFs) play an important role in gene regulation. The interconnections among TFs, chromatin interactions, epigenetic marks and cis-regulatory elements form a complex gene transcription apparatus. Our previous work, ChIP-Array, combined TF binding and transcriptome data to construct gene regulatory networks (GRNs). Here we present an enhanced version, ChIP-Array 2, to integrate additional types of omics data including long-range chromatin interaction, open chromatin region and histone modification data to dissect more comprehensive GRNs involving diverse regulatory components. Moreover, we substantially extended our motif database for human, mouse, rat, fruit fly, worm, yeast and Arabidopsis, and curated large amount of omics data for users to select as input or backend support. With ChIP-Array 2, we compiled a library containing regulatory networks of 18 TFs/chromatin modifiers in mouse embryonic stem cell (mESC). The web server and the mESC library are publicly free and accessible athttp://jjwanglab.org/chip-array.

Integrative analysis of 111 reference human epigenomes

The reference human genome sequence set the stage for studies of genetic variation and its association with human disease, but epigenomic studies lack a similar reference. To address this need, the NIH Roadmap Epigenomics Consortium generated the largest collection so far of human epigenomes for primary cells and tissues. Here we describe the integrative analysis of 111 reference human epigenomes generated as part of the programme, profiled for histone modification patterns, DNA accessibility, DNA methylation and RNA expression. We establish global maps of regulatory elements, define regulatory modules of coordinated activity, and their likely activators and repressors. We show that disease- and trait-associated genetic variants are enriched in tissue-specific epigenomic marks, revealing biologically relevant cell types for diverse human traits, and providing a resource for interpreting the molecular basis of human disease. Our results demonstrate the central role of epigenomic information for understanding gene regulation, cellular differentiation and human disease.

Silencing ECHS1 attenuates the proliferation and induces the autophagy of hepatocellular carcinoma via impairing cell metabolism and activating AMPK

Hepatocellular carcinoma (HCC) is among the most common cancers in the world with a low survival rate. Our previous study showed Short chain enoyl-CoA hydratase (ECHS1) could bind to HBsAg (HBs) and that ECHS1's localization in mitochondria induced HepG2 cell apoptosis. However, the role of the ECHS1 in energy metabolism and autophagy during hepatocellular carcinoma development remains undefined. We aimed to determine what ECHS1 does to energy metabolism and its effects on HCC progression. We performed CCK-8, EdU assays in hepatocellular carcinoma cell lines (HepG2 and HuH7) with stable ECHS1 knock-down. ATP and NADP+/NADPH levels were measured using an colorimetric assay. Our data demonstrated that ECHS1 silencing inhibited cell proliferation and induced autophagy. ECHS1 knockdown did not increase fatty acid synthesis, but decreased cellular ATP. This resulted in AMP-activated protein kinase (AMPK) activation and induced HCC cell autophagy. Our results showed that silencing ECHS1 to attenuate proliferation and induce autophagy may make it a novel cancer therapy target.

Epigenetic modification of MiR-429 promotes liver tumour-initiating cell properties by targeting Rb binding protein 4

OBJECTIVE

Liver tumour-initiating cells (T-ICs) are critical for hepatocarcinogenesis. However, the underlying mechanism regulating the function of liver T-ICs remains unclear.

METHODS

Tissue microarrays containing 242 hepatocellular carcinoma (HCC) samples were used for prognostic analysis. Magnetically activated cell sorting was used to isolate epithelial cell adhesion molecule (EPCAM)-positive cells. The gene expressions affected by miR-429 were determined by arrays. Co-immunoprecipitation was used to study interactions among retinoblastoma protein (RB1), Rb binding protein 4 (RBBP4) and E2F transcription factor 1 (E2F1). The DNA methylation status in CpG islands was detected by quantitative methylation analysis. miRNAs in microvesicles were isolated by a syringe filter system.

RESULTS

The significant prognosis factor miR-429 was upregulated in HCC tissues and also in primary liver T-ICs isolated from clinical samples. The enrichment of miR-429 in EPCAM+ T-ICs contributed to hepatocyte self-renewal, malignant proliferation, chemoresistance and tumorigenicity. A novel functional axis involving miR-429, RBBP4, E2F1 and POU class 5 homeobox 1 (POU5F1 or OCT4) governing the regulation of liver EPCAM+ T-ICs was established in vitro and in vivo. The molecular mechanism regulating miR-429 expression, involving four abnormal hypomethylated sites upstream of the miR-200b/miR-200a/miR-429 cluster, was first defined in both EPCAM+ liver T-ICs and very early-stage HCC tissues. miR-429 secreted by high-expressing cells has the potential to become a proactive signalling molecule to mediate intercellular communication.

CONCLUSIONS

Epigenetic modification of miR-429 can manipulate liver T-ICs by targeting the RBBP4/E2F1/OCT4 axis. This miRNA might be targeted to inactivate T-ICs, thus providing a novel strategy for HCC prevention and treatment.

A common set of distinct features that characterize noncoding RNAs across multiple species

To find signature features shared by various ncRNA sub-types and characterize novel ncRNAs, we have developed a method, RNAfeature, to investigate >600 sets of genomic and epigenomic data with various evolutionary and biophysical scores. RNAfeature utilizes a fine-tuned intra-species wrapper algorithm that is followed by a novel feature selection strategy across species. It considers long distance effect of certain features (e.g. histone modification at the promoter region). We finally narrow down on 10 informative features (including sequences, structures, expression profiles and epigenetic signals). These features are complementary to each other and as a whole can accurately distinguish canonical ncRNAs from CDSs and UTRs (accuracies: >92% in human, mouse, worm and fly). Moreover, the feature pattern is conserved across multiple species. For instance, the supervised 10-feature model derived from animal species can predict ncRNAs in Arabidopsis (accuracy: 82%). Subsequently, we integrate the 10 features to define a set of noncoding potential scores, which can identify, evaluate and characterize novel noncoding RNAs. The score covers all transcribed regions (including unconserved ncRNAs), without requiring assembly of the full-length transcripts. Importantly, the noncoding potential allows us to identify and characterize potential functional domains with feature patterns similar to canonical ncRNAs (e.g. tRNA, snRNA, miRNA, etc) on ∼70% of human long ncRNAs (lncRNAs).

Genome-wide map of regulatory interactions in the human genome

Increasing evidence suggests that interactions between regulatory genomic elements play an important role in regulating gene expression. We generated a genome-wide interaction map of regulatory elements in human cells (ENCODE tier 1 cells, K562, GM12878) using Chromatin Interaction Analysis by Paired-End Tag sequencing (ChIA-PET) experiments targeting six broadly distributed factors. Bound regions covered 80% of DNase I hypersensitive sites including 99.7% of TSS and 98% of enhancers. Correlating this map with ChIP-seq and RNA-seq data sets revealed cohesin, CTCF, and ZNF143 as key components of three-dimensional chromatin structure and revealed how the distal chromatin state affects gene transcription. Comparison of interactions between cell types revealed that enhancer-promoter interactions were highly cell-type-specific. Construction and comparison of distal and proximal regulatory networks revealed stark differences in structure and biological function. Proximal binding events are enriched at genes with housekeeping functions, while distal binding events interact with genes involved in dynamic biological processes including response to stimulus. This study reveals new mechanistic and functional insights into regulatory region organization in the nucleus.

Gene module based regulator inference identifying miR-139 as a tumor suppressor in colorectal cancer

Colorectal cancer is one of the most commonly diagnosed cancer types worldwide. Identification of the key regulators of the altered biological networks is crucial for understanding the complex molecular mechanisms of colorectal cancer. We proposed a gene module based approach to infer key miRNAs regulating the major gene network alterations in cancer tissues. By integrating gene differential expression and co-expression information with a protein-protein interaction network, the differential gene expression modules, which captured the major gene network changes, were identified for colorectal cancer. Then, several key miRNAs, which extensively regulate the gene modules, were inferred by analyzing their target gene enrichment in the modules. Among the inferred candidates, three miRNAs, miR-101, miR-124 and miR-139, are frequently down-regulated in colorectal cancers. The following computational and experimental analyses demonstrate that miR-139 can inhibit cell proliferation and cell cycle G1/S transition. A known oncogene ETS1, a key transcription factor in the gene module, was experimentally verified as a novel target of miR-139. miR-139 was found to be significantly down-regulated in early pathological cancer stages and its expression remained at very low levels in advanced stages. These results indicate that miR-139, inferred by the gene module based approach, should be a key tumor suppressor in early cancer development.

First Report of Myrothecium roridum Causing Leaf Spot of Anubias barteri var. nana in Mainland China

Anubias barteri var. nana (belonging to family Araceae) is one of the most popular ornamental plants for aquaria. It has slow growth and small attractive ornamental leaves. It originated from West Africa, was first described in 1860 by Heinrich Wilhelm Schott (1), and is now commercially planted widely in southern China. A new disease of A. barteri var. nana was observed on leaves. Approximately 30% of the plants in an ornamental nursery in Baiyun district, Gaungzhou, China (113°22'45.15″ E, 23°23'41.15″ N) in July 2013 were found to be infected. The pathogen caused necrotic spots that progressed to form holes on leaves of the infected plants, negatively affecting their ornamental market value. Eventually the diseased plant died. White fungal fruiting bodies and black sporodochia were observed on the surfaces of the symptomatic leaves. A single-spore isolate (myr5) was obtained from the concentric lesions and cultured on potato dextrose agar (PDA). The floccose fungal colonies were white to buff, and black conidiomata were often visible on the surface after 25 days of incubation. The conidia had rounded ends and the average size of conidia was 5.98 ± 0.15 × 2.24 ± 0.08 μm. The rDNA internal transcribed spacer (ITS) region of isolate myr5 was amplified by PCR using ITS1 and ITS4 primer pairs. The amplified product was sequenced and deposited in GenBank (Accession No. KJ572115) and showed 99% identity to Myrothecium roridum isolates BBA 71015 (AJ302001) and BBA 67679 (AJ301995) (unpublished), and 100% identity to the isolate myr2-2 from Dieffenbachia picta 'Camilla' in Taiwan (2). On the basis of morphological and molecular characterization, the fungus causing leaf spot on A. barteri var. nana plants was identified as M. roridum Tode ex Fr (2,3). To confirm the etiology of the disease, Koch's postulates was performed in a greenhouse (28 ± 2°C) using 3-day-old cultures of isolate myr5 and fungal spore suspensions of 1 × 10⁵ conidia ml^-1 containing 0.05% of Tween 20. Sixteen healthy leaves, two each from eight plants, were infiltrated on two different parts close to the midribs with 100 μl of the fungal spore using needleless syringes. Water infiltration was used as the control treatment. Water-soaked brown lesions appeared on leaves 7 days after inoculation, followed by the development of dark concentric rings within the necrotic areas on the surface of the inoculated leaves after 15 days of incubation. These symptoms were similar to those in the naturally infected aquarium plants. No symptoms were observed on any of the water micro-infiltrated plants. The fungus was re-isolated from the inoculated plants, but not from control plants. To our knowledge, this is the first report of Myrothecium leaf spot caused by M. roridum on A. barteri var. nana in mainland China. References: (1) W. Crusio. Primitiae Africanae XII 79 (14):1-48, 1979. (2) C. F. Hong et al. Plant Dis. 97:1253, 2013. (3) M. Tulloch. Mycol. Pap. 130:1-42, 1972.

Exome sequencing identifies a novel frameshift mutation of MYO6 as the cause of autosomal dominant nonsyndromic hearing loss in a Chinese family

Autosomal dominant types of nonsyndromic hearing loss (ADNSHL) are typically postlingual in onset and progressive. High genetic heterogeneity, late onset age, and possible confounding due to nongenetic factors hinder the timely molecular diagnoses for most patients. In this study, exome sequencing was applied to investigate a large Chinese family segregating ADNSHL in which we initially failed to find strong evidence of linkage to any locus by whole-genome linkage analysis. Two affected family members were selected for sequencing. We identified two novel mutations disrupting known ADNSHL genes and shared by the sequenced samples: c.328C>A in COCH (DFNA9) resulting in a p.Q110K substitution and a deletion c. 2814_2815delAA in MYO6 (DFNA22) causing a frameshift alteration p.R939Tfs*2. The pathogenicity of novel coding variants in ADNSHL genes was carefully evaluated by analysis of co-segregation with phenotype in the pedigree and in light of established genotype-phenotype correlations. The frameshift deletion in MYO6 was confirmed as the causative variant for this pedigree, whereas the missense mutation in COCH had no clinical significance. The results allowed us to retrospectively identify the phenocopy in one patient that contributed to the negative finding in the linkage scan. Our clinical data also supported the emerging genotype-phenotype correlation for DFNA22.

Hsa-miR-1246, hsa-miR-320a and hsa-miR-196b-5p inhibitors can reduce the cytotoxicity of Ebola virus glycoprotein in vitro

Ebola virus (EBOV) causes a highly lethal hemorrhagic fever syndrome in humans and has been associated with mortality rates of up to 91% in Zaire, the most lethal strain. Though the viral envelope glycoprotein (GP) mediates widespread inflammation and cellular damage, these changes have mainly focused on alterations at the protein level, the role of microRNAs (miRNAs) in the molecular pathogenesis underlying this lethal disease is not fully understood. Here, we report that the mi-RNAs hsa-miR-1246, hsa-miR-320a and hsa-miR-196b-5p were induced in human umbilical vein endothelial cells (HUVECs) following expression of EBOV GP. Among the proteins encoded by predicted targets of these miRNAs, the adhesion-related molecules tissue factor pathway inhibitor (TFPI), dystroglycan1 (DAG1) and the caspase 8 and FADD-like apoptosis regulator (CFLAR) were significantly downregulated in EBOV GP-expressing HUVECs. Moreover, inhibition of hsa-miR-1246, hsa-miR-320a and hsa-miR-196b-5p, or overexpression of TFPI, DAG1 and CFLAR rescued the cell viability that was induced by EBOV GP. Our results provide a novel molecular basis for EBOV pathogenesis and may contribute to the development of strategies to protect against future EBOV pandemics.

Murine In vitro Memory T Cell Differentiation

Upon pathogen encounter, naïve CD8⁺ T cells are primed and undergo massive clonal expansion. A fraction of effector CD8⁺ T cells remains during the contraction phase and differentiate into memory T cells critical for mounting robust recall responses in response to secondary infection. Low frequency of memory T cells in vivo is a major obstacle to investigate their functional aspects including migration capacity and genetic regulation. Here, we describe detailed protocol for memory T cell differentiation developed by von Andrian's group to generate large number of CD44^hiCD62L^hi antigen-specific memory T cells in vitro.

First Report of Phoma sp. Causing Twisting and Curling of Crown Leaves of Sugarcane in Mainland of China

Sugarcane is a major sugar and the leading energy crop worldwide and Guangxi is the largest sugarcane production area in China (2011 Sugar Annual China, www.gain.fas.usda.gov). During survey of sugarcane crops in September 2012 and June 2013, ~5 to 10% of sugarcane (cvs. FN-40 and ROC22) planted in Chongzuo and Laibing, Guangxi Zhuang Autonomous Region, P. R. China, had twisted and curling symptoms of crown leaves similar to sugarcane Pokkah boeng disease (caused by Fusarium moniliforme Sheldon). The symptoms started appearing as yellowing on midribs and leaf margins that spread further to the entire leaf, along with twisting and curling of crown leaves. The symptomatic leaf tissues (5 × 5 mm) were surface-sterilized by 0.1% HgCl₂ solution for 30 s, followed by rinsing three times in sterile water, placed on potato dextrose agar (PDA), and then incubated in darkness at 28°C. After 3 days of incubation, the isolated fungal colony appeared as white villous, spherical, radial, and dense colorless mycelium from the top, while it was reddish-brown at the bottom and later became grayish. Chlamydospores were also observed with a diameter of 5 to 10 μm and were dark brown, unicellular, intercalary, and smooth. The binucleate hyphae were colorless and transparent. Pycnidia appeared on the colonies after 20 days, and were dark brown, subglobose, and 150 to 230 μm in diameter, and the conidia were ~3 to 7 × 2.5 to 6 μm, unicellular, colorless, and ovoid to oval. The fungal isolates from the symptomatic leaves were obtained and pathogenicity was evaluated. Conidial suspensions (10⁷ CFU/ml) of the single isolate from FN-40 were micro-injected into 20 sugarcane seedlings of cultivar FN-40. Another 20 seedlings were injected with water without conidia as control. The inoculated plants were grown in a growth chamber at 28°C with a 16-h photoperiod. Twisted and curly symptoms similar to the field appeared on the inoculated leaves at 10 days after inoculation, while the control leaves remained asymptomatic. The fungus was re-isolated and identified. Genomic DNA from the cultured fungal isolate was extracted with a modified Fungal DNA Midi Kit (Omega Bio-Tek, Inc., Norcross, GA), and amplified using fungus-conserved primer sequences (ITS1: 5'-TCCGTAGGTGAACCTGCGG-3' and ITS4: 5'-TCCTCCGCTTATTGATATGC-3'). The consensus rDNA-internal transcribed spacer sequence (GenBank Accession No. KC524502) was 100% identical with 97% coverage to the ITS sequence from Phoma sp. 3. TMS-2011 (HQ631000.1) in GenBank (3). The fungus Phoma sp. was identified on the basis of morphological characteristics (2,4) and the ITS sequence of rDNA (1,3). Disease caused by Phoma sp. has been reported earlier on sugarcane from Pakistan, Hawaii, and Taiwan, causing leaf blight and curling (2,4). However, to the best of our knowledge, this is the first report of Phoma sp. causing twisting and curling of crown leaves of sugarcane in mainland of China. References: (1) M. M. Aveskamp et al. Stud. Mycol. 65:1, 2010. (2) A. Sanguino and H. Tokeski. ISSCT Proc. 17:1555, 1980. (3) P. Shrestha et al. Appl. Environ. Microbiol. 77:5490, 2011. (4) Z. N. Wang. Rep. Taiwan Sugar Res. Inst. 129:1, 1980.

Nucleosome eviction and multiple co-factor binding predict estrogen-receptor-alpha-associated long-range interactions

Many enhancers regulate their target genes via long-distance interactions. High-throughput experiments like ChIA-PET have been developed to map such largely cell-type-specific interactions between cis-regulatory elements genome-widely. In this study, we integrated multiple types of data in order to reveal the general hidden patterns embedded in the ChIA-PET data. We found characteristic distance features related to promoter-promoter, enhancer-enhancer and insulator-insulator interactions. Although a protein may have many binding sites along the genome, our hypothesis is that those sites that share certain open chromatin structure can accommodate relatively larger protein complex consisting of specific regulatory and 'bridging' factors, and may be more likely to form robust long-range deoxyribonucleic acid (DNA) loops. This hypothesis was validated in the estrogen receptor alpha (ERα) ChIA-PET data. An efficient classifier was built to predict ERα-associated long-range interactions solely from the related ChIP-seq data, hence linking distal ERα-dependent enhancers to their target genes. We further applied the classifier to generate additional novel interactions, which were undetected in the original ChIA-PET paper but were validated by other independent experiments. Our work provides a new insight into the long-range chromatin interactions through deeper and integrative ChIA-PET data analysis and demonstrates DNA looping predictability from ordinary ChIP-seq data.

Chd5 orchestrates chromatin remodelling during sperm development

One of the most remarkable chromatin remodelling processes occurs during spermiogenesis, the post-meiotic phase of sperm development during which histones are replaced with sperm-specific protamines to repackage the genome into the highly compact chromatin structure of mature sperm. Here we identify Chromodomain helicase DNA binding protein 5 (Chd5) as a master regulator of the histone-to-protamine chromatin remodelling process. Chd5 deficiency leads to defective sperm chromatin compaction and male infertility in mice, mirroring the observation of low CHD5 expression in testes of infertile men. Chd5 orchestrates a cascade of molecular events required for histone removal and replacement, including histone 4 (H4) hyperacetylation, histone variant expression, nucleosome eviction and DNA damage repair. Chd5 deficiency also perturbs expression of transition proteins (Tnp1/Tnp2) and protamines (Prm1/2). These findings define Chd5 as a multi-faceted mediator of histone-to-protamine replacement and depict the cascade of molecular events underlying this process of extensive chromatin remodelling.

Assembly and validation of versatile transcription activator-like effector libraries

The ability to perturb individual genes in genome-wide experiments has been instrumental in unraveling cellular and disease properties. Here we introduce, describe the assembly, and demonstrate the use of comprehensive and versatile transcription activator-like effector (TALE) libraries. As a proof of principle, we built an 11-mer library that covers all possible combinations of the nucleotides that determine the TALE-DNA binding specificity. We demonstrate the versatility of the methodology by constructing a constraint library, customized to bind to a known p53 motif. To verify the functionality in assays, we applied the 11-mer library in yeast-one-hybrid screens to discover TALEs that activate human SCN9A and miR-34b respectively. Additionally, we performed a genome-wide screen using the complete 11-mer library to confirm known genes that confer cycloheximide resistance in yeast. Considering the highly modular nature of TALEs and the versatility and ease of constructing these libraries we envision broad implications for high-throughput genomic assays.

Three-dimensional visualization of rat brain microvasculature following permanent focal ischaemia by synchrotron radiation

OBJECTIVE

Identifying morphological changes that occur in microvessels under both normal and ischaemic conditions is crucial for understanding and treating stroke. However, conventional imaging techniques are not able to detect microvessels on a micron or sub-micron scale without angiography. In the present study, synchrotron radiation (SR)-based X-ray in-line phase contrast imaging (ILPCI) was used to acquire high-resolution and high-contrast images of rat brain tissues in both normal and ischaemic states.

METHODS

ILPCI was performed at the Shanghai Synchrotron Radiation Facility, Shanghai, China, without the use of contrast agents. CT slices were reformatted and then converted into three-dimensional (3D) reconstruction images to analyse subtle details of the cerebral microvascular network.

RESULTS

By using ILPCI, brain vessels up to 11.8 μm in diameter were resolved. The number of cortical and penetrating arteries detected were found to undergo a remarkable decrease within the infarct area. 3 days after permanent ischaemia, vascular masses were also observed in the peripheral region of the infarcts.

CONCLUSION

SR-based ILPCI-CT can serve as a powerful tool to accurately visualize brain microvasculature. The morphological parameters of blood vessels in both CT slices and 3D reconstructions were determined, and this approach has great potential for providing an effective diagnosis and evaluation for rehabilitation therapy for stroke.

ADVANCES IN KNOWLEDGE

In the absence of contrast agent, the 3D morphologies of the brain microvasculature in normal and stroke rats were obtained using SR-based ILPCI. SR imaging is a sensitive and promising method which can be used to explore primary brain function.

Murine in vivo CD8+ T Cell Killing Assay

Antigen-specific killing ability of effector CD8⁺ T cells is critical for protective immunity against infection. Here, we describe in vivo cytotoxic T cell assay to examine effector function of antigen-specific CD8⁺ T cells. Mice infected with Listeria monocytogenes (L. monocytogenes) expressing chicken ovalbumin as a model antigen mount ovalbumin-specific CD8⁺ T cell responses. Effector CD8⁺ T cell function in vivo is determined by mixed transfer of OVA peptide-pulsed target cells with control target cells into the previously immunized mice. Difference in CFSE expression levels clearly marks two distinct populations: Antigen-pulsed target cells-CFSE^low vs. unpulsed target cells-CFSE^hi. The frequencies between antigen-pulsed target cells and control target cells are used as readouts of antigen-specific killing.

MIROR: a method for cell-type specific microRNA occupancy rate prediction

This work provides a novel method to quantitatively predict miRNA–mRNA interactions in a specific cell type.

Characterizing the strand-specific distribution of non-CpG methylation in human pluripotent cells

DNA methylation is an important defense and regulatory mechanism. In mammals, most DNA methylation occurs at CpG sites, and asymmetric non-CpG methylation has only been detected at appreciable levels in a few cell types. We are the first to systematically study the strand-specific distribution of non-CpG methylation. With the divide-and-compare strategy, we show that CHG and CHH methylation are not intrinsically different in human embryonic stem cells (ESCs) and induced pluripotent stem cells (iPSCs). We also find that non-CpG methylation is skewed between the two strands in introns, especially at intron boundaries and in highly expressed genes. Controlling for the proximal sequences of non-CpG sites, we show that the skew of non-CpG methylation in introns is mainly guided by sequence skew. By studying subgroups of transposable elements, we also found that non-CpG methylation is distributed in a strand-specific manner in both short interspersed nuclear elements (SINE) and long interspersed nuclear elements (LINE), but not in long terminal repeats (LTR). Finally, we show that on the antisense strand of Alus, a non-CpG site just downstream of the A-box is highly methylated. Together, the divide-and-compare strategy leads us to identify regions with strand-specific distributions of non-CpG methylation in humans.

BS-Seeker2: a versatile aligning pipeline for bisulfite sequencing data

BACKGROUND

DNA methylation is an important epigenetic modification involved in many biological processes. Bisulfite treatment coupled with high-throughput sequencing provides an effective approach for studying genome-wide DNA methylation at base resolution. Libraries such as whole genome bisulfite sequencing (WGBS) and reduced represented bisulfite sequencing (RRBS) are widely used for generating DNA methylomes, demanding efficient and versatile tools for aligning bisulfite sequencing data.

RESULTS

We have developed BS-Seeker2, an updated version of BS Seeker, as a full pipeline for mapping bisulfite sequencing data and generating DNA methylomes. BS-Seeker2 improves mappability over existing aligners by using local alignment. It can also map reads from RRBS library by building special indexes with improved efficiency and accuracy. Moreover, BS-Seeker2 provides additional function for filtering out reads with incomplete bisulfite conversion, which is useful in minimizing the overestimation of DNA methylation levels. We also defined CGmap and ATCGmap file formats for full representations of DNA methylomes, as part of the outputs of BS-Seeker2 pipeline together with BAM and WIG files.

CONCLUSIONS

Our evaluations on the performance show that BS-Seeker2 works efficiently and accurately for both WGBS data and RRBS data. BS-Seeker2 is freely available at http://pellegrini.mcdb.ucla.edu/BS_Seeker2/ and the Galaxy server.

DNA methylation data analysis and its application to cancer research

With the rapid development of genome-wide high-throughput technologies, including expression arrays, SNP arrays and next-generation sequencing platforms, enormous amounts of molecular data have been generated and deposited in the public domain. The application of computational approaches is required to yield biological insights from this enormous, ever-growing resource. A particularly interesting subset of these resources is related to epigenetic regulation, with DNA methylation being the most abundant data type. In this paper, we will focus on the analysis of DNA methylation data and its application to cancer studies. We first briefly review the molecular techniques that generate such data, much of which has been obtained with the use of the most recent version of Infinium HumanMethylation450 BeadChip(®) technology (Illumina, CA, USA). We describe the coverage of the methylome by this technique. Several examples of data mining are provided. However, it should be understood that reliance on a single aspect of epigenetics has its limitations. In the not too distant future, these defects may be rectified, providing scientists with previously unavailable opportunities to explore in detail the role of epigenetics in cancer and other disease states.

Integrated omics study delineates the dynamics of lipid droplets in Rhodococcus opacus PD630

Rhodococcus opacus strain PD630 (R. opacus PD630), is an oleaginous bacterium, and also is one of few prokaryotic organisms that contain lipid droplets (LDs). LD is an important organelle for lipid storage but also intercellular communication regarding energy metabolism, and yet is a poorly understood cellular organelle. To understand the dynamics of LD using a simple model organism, we conducted a series of comprehensive omics studies of R. opacus PD630 including complete genome, transcriptome and proteome analysis. The genome of R. opacus PD630 encodes 8947 genes that are significantly enriched in the lipid transport, synthesis and metabolic, indicating a super ability of carbon source biosynthesis and catabolism. The comparative transcriptome analysis from three culture conditions revealed the landscape of gene-altered expressions responsible for lipid accumulation. The LD proteomes further identified the proteins that mediate lipid synthesis, storage and other biological functions. Integrating these three omics uncovered 177 proteins that may be involved in lipid metabolism and LD dynamics. A LD structure-like protein LPD06283 was further verified to affect the LD morphology. Our omics studies provide not only a first integrated omics study of prokaryotic LD organelle, but also a systematic platform for facilitating further prokaryotic LD research and biofuel development.

FastDMA: an infinium humanmethylation450 beadchip analyzer

DNA methylation is vital for many essential biological processes and human diseases. Illumina Infinium HumanMethylation450 Beadchip is a recently developed platform studying genome-wide DNA methylation state on more than 480,000 CpG sites and a few CHG sites with high data quality. To analyze the data of this promising platform, we developed FastDMA which can be used to identify significantly differentially methylated probes. Besides single probe analysis, FastDMA can also do region-based analysis for identifying the differentially methylated region (DMRs). A uniformed statistical model, analysis of covariance (ANCOVA), is used to achieve all the analyses in FastDMA. We apply FastDMA on three large-scale DNA methylation datasets from The Cancer Genome Atlas (TCGA) and find many differentially methylated genomic sites in different types of cancer. On the testing datasets, FastDMA shows much higher computational efficiency than current tools. FastDMA can benefit the data analyses of large-scale DNA methylation studies with an integrative pipeline and a high computational efficiency. The software is freely available via http://bioinfo.au.tsinghua.edu.cn/software/fastdma/.

The transcription factor Foxo1 controls central-memory CD8+ T cell responses to infection

Memory T cells protect hosts from pathogen reinfection, but how these cells emerge from a pool of antigen-experienced T cells is unclear. Here, we show that mice lacking the transcription factor Foxo1 in activated CD8+ T cells have defective secondary, but not primary, responses to Listeria monocytogenes infection. Compared to short-lived effector T cells, memory-precursor T cells expressed higher amounts of Foxo1, which promoted their generation and maintenance. Chromatin immunoprecipitation sequencing revealed the transcription factor Tcf7 and the chemokine receptor Ccr7 as Foxo1-bound target genes, which have critical functions in central-memory T cell differentiation and trafficking. These findings demonstrate that Foxo1 is selectively incorporated into the genetic program that regulates memory CD8+ T cell responses to infection.

ProteoMirExpress: inferring microRNA and protein-centered regulatory networks from high-throughput proteomic and mRNA expression data

MicroRNAs (miRNAs) regulate gene expression through translational repression and RNA degradation. Recently developed high-throughput proteomic methods measure gene expression changes at protein level and therefore can reveal the direct effects of miRNAs' translational repression. Here, we present a web server, ProteoMirExpress, that integrates proteomic and mRNA expression data together to infer miRNA-centered regulatory networks. With both types of high-throughput data from the users, ProteoMirExpress is able to discover not only miRNA targets that have decreased mRNA, but also subgroups of targets with suppressed proteins whose mRNAs are not significantly changed or with decreased mRNA whose proteins are not significantly changed, which are usually ignored by most current methods. Furthermore, both direct and indirect targets of miRNAs can be detected. Therefore, ProteoMirExpress provides more comprehensive miRNA-centered regulatory networks. We used several published data to assess the quality of our inferred networks and prove the value of our server. ProteoMirExpress is available online, with free access to academic users.

Population dynamics of cancer cells with cell state conversions

Cancer stem cell (CSC) theory suggests a cell-lineage structure in tumor cells in which CSCs are capable of giving rise to the other non-stem cancer cells (NSCCs) but not vice versa. However, an alternative scenario of bidirectional interconversions between CSCs and NSCCs was proposed very recently. Here we present a general population model of cancer cells by integrating conventional cell divisions with direct conversions between different cell states, namely, not only can CSCs differentiate into NSCCs by asymmetric cell division, NSCCs can also dedifferentiate into CSCs by cell state conversion. Our theoretical model is validated when applying the model to recent experimental data. It is also found that the transient increase in CSCs proportion initiated from the purified NSCCs subpopulation cannot be well predicted by the conventional CSC model where the conversion from NSCCs to CSCs is forbidden, implying that the cell state conversion is required especially for the transient dynamics. The theoretical analysis also gives the condition such that our general model can be equivalently reduced into a simple Markov chain with only cell state transitions keeping the same cell proportion dynamics.

Computational identification of active enhancers in model organisms

As a class of cis-regulatory elements, enhancers were first identified as the genomic regions that are able to markedly increase the transcription of genes nearly 30years ago. Enhancers can regulate gene expression in a cell-type specific and developmental stage specific manner. Although experimental technologies have been developed to identify enhancers genome-wide, the design principle of the regulatory elements and the way they rewire the transcriptional regulatory network tempo-spatially are far from clear. At present, developing predictive methods for enhancers, particularly for the cell-type specific activity of enhancers, is central to computational biology. In this review, we survey the current computational approaches for active enhancer prediction and discuss future directions.

Epigenomic analysis of multilineage differentiation of human embryonic stem cells

Epigenetic mechanisms have been proposed to play crucial roles in mammalian development, but their precise functions are only partially understood. To investigate epigenetic regulation of embryonic development, we differentiated human embryonic stem cells into mesendoderm, neural progenitor cells, trophoblast-like cells, and mesenchymal stem cells and systematically characterized DNA methylation, chromatin modifications, and the transcriptome in each lineage. We found that promoters that are active in early developmental stages tend to be CG rich and mainly engage H3K27me3 upon silencing in nonexpressing lineages. By contrast, promoters for genes expressed preferentially at later stages are often CG poor and primarily employ DNA methylation upon repression. Interestingly, the early developmental regulatory genes are often located in large genomic domains that are generally devoid of DNA methylation in most lineages, which we termed DNA methylation valleys (DMVs). Our results suggest that distinct epigenetic mechanisms regulate early and late stages of ES cell differentiation.

OLego: fast and sensitive mapping of spliced mRNA-Seq reads using small seeds

A crucial step in analyzing mRNA-Seq data is to accurately and efficiently map hundreds of millions of reads to the reference genome and exon junctions. Here we present OLego, an algorithm specifically designed for de novo mapping of spliced mRNA-Seq reads. OLego adopts a multiple-seed-and-extend scheme, and does not rely on a separate external aligner. It achieves high sensitivity of junction detection by strategic searches with small seeds (∼14 nt for mammalian genomes). To improve accuracy and resolve ambiguous mapping at junctions, OLego uses a built-in statistical model to score exon junctions by splice-site strength and intron size. Burrows–Wheeler transform is used in multiple steps of the algorithm to efficiently map seeds, locate junctions and identify small exons. OLego is implemented in C++ with fully multithreaded execution, and allows fast processing of large-scale data. We systematically evaluated the performance of OLego in comparison with published tools using both simulated and real data. OLego demonstrated better sensitivity, higher or comparable accuracy and substantially improved speed. OLego also identified hundreds of novel micro-exons (<30 nt) in the mouse transcriptome, many of which are phylogenetically conserved and can be validated experimentally in vivo. OLego is freely available at http://zhanglab.c2b2.columbia.edu/index.php/OLego.