Search

Summary Androgenetic haploid embryonic stem cells (AG-haESCs) hold great promise for exploring gene functions and generating gene-edited semi-cloned (SC) mice. However, the high incidence of self-diploidization and low efficiency of SC mouse production are major obstacles preventing widespread use of these cells. Moreover, although SC mice generation could be greatly improved by knocking out the differentially methylated regions of two imprinted genes, 50% of the SC mice did not survive into adulthood. Here, we found that the genome-wide DNA methylation level in AG-haESCs is extremely low. Subsequently, downregulation of both de novo methyltransferase Dnmt3b and other methylation-related genes was determined to be responsible for DNA hypomethylation. We further demonstrated that ectopic expression of Dnmt3b in AG-haESCs could effectively improve DNA methylation level, and the high incidence of self-diploidization could be markedly rescued. More importantly, the developmental potential of SC embryos was improved, and most SC mice could survive into adulthood., Graphical Abstract, Highlights • Dnmt3b restores methylation of hypomethylated AG-haESCs • Dnmt3b suppresses self-diploidization of AG-haESCs • Dnmt3b promotes developmental potential of SC embryos • Dnmt3b improves SC mice survival to adulthood, Ectopic expression of Dnmt3b could rescue DNA methylation level in repetitive sequences of hypomethylated AG-haESCs, suppress high incidence of self-diploidization, and promote developmental potential of SC embryos, and most SC mice could survive into adulthood.

// Bo Wang 1, * , Wei Du 2, * , Yutao Jia 1 , Xiaobai Zhang 1 , Guorui Ma 1 1 Department of Paediatrics, The First Affiliated Hospital of Henan University of Science and Technology, Luoyang, Henan, China 2 Department of Medical Laboratory, Luoyang Central Hospital, Luoyang, Henan, China * These authors have contributed equally to this work Correspondence to: Guorui Ma, email: aquacity@126.com Keywords: type 1 diabetes, cytotoxic T-lymphocyte associated protein 4, polymorphism, risk, meta-analysis Received: October 05, 2016 Accepted: December 13, 2016 Published: January 02, 2017 ABSTRACT Type 1 diabetes (T1D) is a heritable disease associated with multiple genetic variants. This systematic review and meta-analysis assessed the correlation between cytotoxic T-lymphocyte-associated protein 4(CTLA-4) +49A/G polymorphisms and the risk of T1D in children. The random effects model was used to estimate the related odds ratios (ORs) and 95% confidence intervals (CIs). Trial sequential analysis (TSA) was used to determine whether the currently available evidence was sufficient and conclusive. Our results indicated that CTLA-4 gene polymorphisms significantly increased the risk of childhood T1D in an allelic model (G vs. A: OR=1.33, 95%CI=1.19-1.48; I 2 =44.0% and P =0.001for heterogeneity) and a codominant model (GG vs. AA: OR=1.75, 95%CI=1.37-2.24; I 2 =57.5% and P =0.001for heterogeneity; GA vs. AA: OR=1.26, 95%CI=1.09-1.46; I 2 =40.4% and P =0.036for heterogeneity). Subgroup analysis results indicated that the ORs were higher in the Asian population (OR allelic model =1.60, OR GG vs. AA =2.46 and OR GA vs. AA =1.58) than the Caucasian population (OR allelic model ==1.24, OR GG vs. AA =1.55 and OR GA vs. AA =1.19). The TSA results indicated that the evidence of the effect was sufficient. In conclusion, CTLA4 +49A/G polymorphisms increased the risk of T1D in children, and CTLA4 +49A/G can be considered to be a genetic marker for T1D in children.

SummaryThe core pluripotency factor Oct4 plays key roles in somatic cell reprogramming through transcriptional control. Here, we profile Oct4 occupancy, epigenetic changes, and gene expression in reprogramming. We find that Oct4 binds in a hierarchical manner to target sites with primed epigenetic modifications. Oct4 binding is temporally continuous and seldom switches between bound and unbound. Oct4 occupancy in most of promoters is maintained throughout the entire reprogramming process. In contrast, somatic cell-specific enhancers are silenced in the early and intermediate stages, whereas stem cell-specific enhancers are activated in the late stage in parallel with cell fate transition. Both epigenetic remodeling and Oct4 binding contribute to the hyperdynamic enhancer signature transitions. The hierarchical Oct4 bindings are associated with distinct functional themes at different stages. Collectively, our results provide a comprehensive molecular roadmap of Oct4 binding in concert with epigenetic rearrangements and rich resources for future reprogramming studies.

Drosophila neural development undergoes extensive chromatin remodeling and precise epigenetic regulation. However, the roles of chromatin remodeling in establishment and maintenance of cell identity during cell fate transition remain enigmatic. Here, we compared the changes in gene expression, as well as the dynamics of nucleosome positioning and key histone modifications between the four major neural cell types during Drosophila neural development. We find that the neural progenitors can be separated from the terminally differentiated cells based on their gene expression profiles, whereas nucleosome distribution in the flanking regions of transcription start sites fails to identify the relationships between the progenitors and the differentiated cells. H3K27me3 signal in promoters and enhancers can not only distinguish the progenitors from the differentiated cells but also identify the differentiation path of the neural stem cells (NSCs) to the intermediate progenitor cells to the glial cells. In contrast, H3K9ac signal fails to identify the differentiation path, although it activates distinct sets of genes with neuron-specific and glia-related functions during the differentiation of the NSCs into neurons and glia, respectively. Together, our study provides novel insights into the crucial roles of chromatin remodeling in determining cell type during Drosophila neural development. Keywords: Nucleosome, Histone modification, Neural stem cell, Neuron, Glia

It has long been known that trypanosomes regulate mitochondrial biogenesis during the life cycle of the parasite; however, the mitochondrial protein inventory (MitoCarta) and its regulation remain unknown. We present a novel computational method for genome-wide prediction of mitochondrial proteins using a support vector machine-based classifier with approximately 90% prediction accuracy. Using this method, we predicted the mitochondrial localization of 468 proteins with high confidence and have experimentally verified the localization of a subset of these proteins. We then applied a recently developed parallel sequencing technology to determine the expression profiles and the splicing patterns of a total of 1065 predicted MitoCarta transcripts during the development of the parasite, and showed that 435 of the transcripts significantly changed their expressions while 630 remain unchanged in any of the three life stages analyzed. Furthermore, we identified 298 alternatively splicing events, a small subset of which could lead to dual localization of the corresponding proteins.

Previous phylogenetic analyses have led to incongruent evolutionary relationships between tree shrews and other suborders of Euarchontoglires. What caused the incongruence remains elusive. In this study, we identified 6845 orthologous genes between seventeen placental mammals. Tree shrews and Primates were monophyletic in the phylogenetic trees derived from the first or/and second codon positions whereas tree shrews and Glires formed a monophyly in the trees derived from the third or all codon positions. The same topology was obtained in the phylogeny inference using the slowly and fast evolving genes, respectively. This incongruence was likely attributed to the fast substitution rate in tree shrews and Glires. Notably, sequence GC content only was not informative to resolve the controversial phylogenetic relationships between tree shrews, Glires, and Primates. Finally, estimation in the confidence of the tree selection strongly supported the phylogenetic affiliation of tree shrews to Primates as a monophyly.

Neurons are a key component of the nervous system and differentiate from multipotent neural stem cells (NSCs). Chromatin remodeling has a critical role in the differentiation process. However, its in vivo epigenetic regulatory role remains unknown. We show here that nucleosome depletion regions (NDRs) form in both proximal promoters and distal enhancers during NSCs differentiating into neurons in the early Drosophila embryonic development. NDR formation in the regulatory regions involves nucleosome shift and eviction. Nucleosome occupancy in promoter NDRs is inversely proportional to the gene activity. Genes with promoter NDR formation during differentiation are enriched for functions related to neuron development and maturation. Active histone-modification signals (H3K4me3 and H3K9ac) in promoters are gained in neurons in two modes: de novo establishment to high levels or increase from the existing levels in NSCs. The gene sets corresponding to the two modes have different neuron-related functions. Dynamic changes of H3K27ac and H3K9ac signals in enhancers and promoters synergistically repress genes associated with neural stem or progenitor cell-related pluripotency and upregulate genes associated with neuron projection morphogenesis, neuron differentiation, and so on. Our results offer new insights into chromatin remodeling during in vivo neuron development and lay a foundation for its epigenetic regulatory mechanism study of other lineage specification.

Chronic kidney disease is becoming a global public health problem, which will usually cause uremia at the end stage of chronic kidney failure. So far, kidney transplant is the most effective and proper therapy for uremia, however, the short supply of matched donor kidney has been a persistent bottleneck for transplantation. HLA matching of HLA-A, -B and -DRB1 loci is very important for the allocation of kidney transplants. In this study, we investigated genotypes of HLA-A, -B and -DRB1 loci based on 1,464 uremia patients and 10,000 unrelated healthy individuals in Henan province of China, and compared the frequency distribution of these HLA alleles and corresponding haplotypes between patient and healthy groups. We detected 23 HLA-A, 49 HLA-B and 17 HLA-DRB1 alleles in total. The predominant alleles of HLA-A, -B and -DRB1 loci in patients are the same as those in healthy group. The seven most frequent alleles account for about 87%, 50%, and 77% at HLA-A, -B and -DRB1 loci, respectively. The haplotypes (combinations of HLA-A, -B, and -DRB1) with significantly different frequency between patients and controls mostly account for less than 1%. Overall, this suggests that HLA matching is not a potential difficulty for kidney transplant of uremia patients. However, three of the top seven frequent HLA-DRB1 alleles have a significantly different distribution in patients and controls, while only one alleles for HLA-B and zero for HLA-A loci. These HLA-DRB1 alleles may be closely associated with uremia. This study sheds new lights on the composition and difference of HLA genotypes in uremia patients and healthy populations in Central China that can serve as a guide to HLA matching for kidney transplants and a resource for HLA typing-related studies.

MicroRNAs (miRNAs) are a class of small non-coding RNAs discovered in recent years, which are found to play important regulatory roles in various organisms. As the number of experimentally validated miRNAs is rapidly increasing, systematic analysis on the characteristics of these known miRNAs is necessary and indispensable, especially for computational prediction of new miRNAs. We extensively analyzed precursor sequences for all experimentally validated mature miRNAs in metazoan species, focusing on the characteristics at the level of primary sequences and secondary structures. An observation over the secondary structures of 2729 miRNA precursors (pre-miRNAs) reveals that these hairpin structures can be approximately classified into two types: one with a hairpin loop, and the other with multiple loops. Interestingly, the two types of pre-miRNAs show significant differences in both sequence and structure characteristics, and our study indicates that separate consideration on each type of pre-miRNAs is more reasonable, especially in computational prediction. Besides, we develop a new criterion called mAMFE which shows robust discriminative power in distinguishing pre-miRNAs against other RNAs, thus can potentially serve as a discriminative feature in prediction of new pre-miRNAs.

ATP-dependent chromatin remodeling complexes regulate nucleosome organizations. In Drosophila, gene Brm encodes the core Brahma complex, the ATPase subunit of SWI/SNF class of chromatin remodelers. Its role in modulating the nucleosome landscape in vivo is unclear. In this study, we knocked down Brm in Drosophila third instar larvae to explore the changes in nucleosome profiles and global gene transcription. The results show that Brm knockdown leads to nucleosome occupancy changes throughout the entire genome with a bias in occupancy decrease. In contrast, the knockdown has limited impacts on nucleosome position shift. The knockdown also alters another important physical property of nucleosome positioning, fuzziness. Nucleosome position shift, gain or loss and fuzziness changes are all enriched in promoter regions. Nucleosome arrays around the 5′ ends of genes are reorganized in five patterns as a result of Brm knockdown. Intriguingly, the concomitant changes in the genes adjacent to the Brahma-dependent remodeling regions have important roles in development and morphogenesis. Further analyses reveal abundance of AT-rich motifs for transcription factors in the remodeling regions.

Polybrominated diphenyl ethers (PBDEs) are widely used as flame-retardant additives in consumer and household products and can escape into the environment over time. PBDEs have become a global environmental organic pollutant due to the properties of persistence, toxicity, and bioaccumulation. The well-studied toxic effects of PBDEs mainly include thyroid hormone disruption and neurotoxicity. There is no consistent conclusions on the carcinogenic potential of PBDEs to date. Here, we explored the toxic effects of BDE-209 on human embryonic kidney cells (HEK293T). The comparison of the gene expression profiles of HEK293T cells with BDE-209 treatment and the negative control found that BDE-209 exposure may alter nucleosome organization through significantly changing the expression of histone gene clusters. The remodeled chromatin structure could further disturb systemic lupus erythematosus as one of the toxic effects of BDE-209. Additionally, gene sets of different cancer modules are positively correlated with BDE-209 exposure. This suggests that BDE-209 has carcinogenic potential for a variety of tumors. Collectively, BDE-209 has a broader toxicity not limited to disruption of thyroid hormone-related biological processes. Notably, the toxic effects of BDE-209 dissolved in dimethyl sulfoxide (DMSO) is not the simply additive effects of BDE-209 and DMSO alone.

Bisphenol A (BPA) is an environmental endocrine disruptor which has been detected in human bodies. Many studies have implied that BPA exposure is harmful to human health. Previous studies mainly focused on BPA effects on estrogen receptor (ER)-positive cells. Genome-wide impacts of BPA on gene expression in ER-negative cells is unclear. In this study, we performed RNA-seq to characterize BPA-induced cellular and molecular impacts on ER-negative HEK293 cells. The microscopic observation showed that low-dose BPA exposure did not affect cell viability and morphology. Gene expression profiling analysis identified a list of differentially expressed genes in response to BPA exposure in HEK293 cells. These genes were involved in variable important biological processes including ion transport, cysteine metabolic process, apoptosis, DNA damage repair, etc. Notably, BPA up-regulated the expression of ERCC5 encoding a DNA endonuclease for nucleotide-excision repair. Further electrochemical experiment showed that BPA induced significant DNA damage in ER-positive MCF-7 cells but not in ER-negative HEK293 cells. Collectively, our study revealed that ER-negative HEK293 cells employed mechanisms in response to BPA exposure different from ER-positive cells.

microDoR Human is a web server tool based on machine learning algorithm to predict Human miRNA-mediated gene silencing: mRNA degradation or translational repression. For a miRNA-mRNA interaction, microDoR Human identifies whether the target mRNA is cleaved or translationally inhibite. microDoR Human is expected to serve as a useful tool for understanding of miRNA-mediated gene regulation. The microDoR Human 1.0 can be available in http://reprod.njmu.edu.cn/microdor /.

Computational analysis of microarray data has provided an effective way to identify disease-related genes. Traditional disease gene selection methods from microarray data such as statistical test always focus on differentially expressed genes in different samples by individual gene prioritization. These traditional methods might miss differentially coexpressed (DCE) gene subsets because they ignore the interaction between genes. In this paper, MIClique algorithm is proposed to identify DEC gene subsets based on mutual information and clique analysis. Mutual information is used to measure the coexpression relationship between each pair of genes in two different kinds of samples. Clique analysis is a commonly used method in biological network, which generally represents biological module of similar function. By applying the MIClique algorithm to real gene expression data, some DEC gene subsets which correlated under one experimental condition but uncorrelated under another condition are detected from the graph of colon dataset and leukemia dataset.

Protein turnover metabolism plays important roles in cell cycle progression, signal transduction, and differentiation. Those proteins with short half-lives are involved in various regulatory processes. To better understand the regulation of cell process, it is important to study the key sequence-derived factors affecting short-lived protein degradation. Until now, most of protein half-lives are still unknown due to the difficulties of traditional experimental methods in measuring protein half-lives in human cells. To investigate the molecular determinants that affect short-lived proteins, a computational method was proposed in this work to recognize short-lived proteins based on sequence-derived features in human cells. In this study, we have systematically analyzed many features that perhaps correlated with short-lived protein degradation. It is found that a large fraction of proteins with signal peptides and transmembrane regions in human cells are of short half-lives. We have constructed an SVM-based classifier to recognize short-lived proteins, due to the fact that short-lived proteins play pivotal roles in the control of various cellular processes. By employing the SVM model on human dataset, we achieved 80.8% average sensitivity and 79.8% average specificity, respectively, on ten testing dataset (TE1-TE10). We also obtained 89.9%, 99% and 83.9% of average accuracy on an independent validation datasets iTE1, iTE2 and iTE3 respectively. The approach proposed in this paper provides a valuable alternative for recognizing the short-lived proteins in human cells, and is more accurate than the traditional N-end rule. Furthermore, the web server SProtP (http://reprod.njmu.edu.cn/sprotp) has been developed and is freely available for users.