Competitive Chemical Reaction Kinetic Model of Nucleosome Assembly Using the Histone Variant H2A.Z and H2A In Vitro.

Nucleosomes not only serve as the basic building blocks for eukaryotic chromatin but also regulate many biological processes, such as DNA replication, repair, and recombination. To modulate gene expression in vivo, the histone variant H2A.Z can be dynamically incorporated into the nucleosome. However, the assembly dynamics of H2A.Z-containing nucleosomes remain elusive. Here, we demonstrate that our previous chemical kinetic model for nucleosome assembly can be extended to H2A.Z-containing nucleosome assembly processes. The efficiency of H2A.Z-containing nucleosome assembly, like that of canonical nucleosome assembly, was also positively correlated with the total histone octamer concentration, reaction rate constant, and reaction time. We expanded the kinetic model to represent the competitive dynamics of H2A and H2A.Z in nucleosome assembly, thus providing a novel method through which to assess the competitive ability of histones to assemble nucleosomes. Based on this model, we confirmed that histone H2A has a higher competitive ability to assemble nucleosomes in vitro than histone H2A.Z. Our competitive kinetic model and experimental results also confirmed that in vitro H2A.Z-containing nucleosome assembly is governed by chemical kinetic principles.

The Metal-binding Protein Atlas (MbPA): An Integrated Database for Curating Metalloproteins in All Aspects.

Metal-binding proteins are essential for the vital activities and engage in their roles by acting in concert with metal cations. MbPA (The Metal-binding Protein Atlas) is the most comprehensive resource up to now dedicated to curating metal-binding proteins. Currently, it contains 106,373 entries and 440,187 sites related to 54 metals and 8169 species. Users can view all metal-binding proteins and species-specific proteins in MbPA. There are also metal-proteomics data that quantitatively describes protein expression in different tissues and organs. By analyzing the data of the amino acid residues at the metal-binding site, it is found that about 80% of the metal ions tend to bind to cysteine, aspartic acid, glutamic acid, and histidine. Moreover, we use Diversity Measure to confirm that the diversity of metal-binding is specific in different area of periodic table, and further elucidate the binding modes of 19 transition metals on 20 amino acids. In addition, MbPA also embraces 6855 potential pathogenic mutations related to metalloprotein. The resource is freely available at http://bioinfor.imu.edu.cn/mbpa.

A cost-effective machine learning-based method for preeclampsia risk assessment and driver genes discovery.

BACKGROUND: The placenta, as a unique exchange organ between mother and fetus, is essential for successful human pregnancy and fetal health. Preeclampsia (PE) caused by placental dysfunction contributes to both maternal and infant morbidity and mortality. Accurate identification of PE patients plays a vital role in the formulation of treatment plans. However, the traditional clinical methods of PE have a high misdiagnosis rate. RESULTS: Here, we first designed a computational biology method that used single-cell transcriptome (scRNA-seq) of healthy pregnancy (38 wk) and early-onset PE (28-32 wk) to identify pathological cell subpopulations and predict PE risk. Based on machine learning methods and feature selection techniques, we observed that the Tuning ReliefF (TURF) score hybrid with XGBoost (TURF_XGB) achieved optimal performance, with 92.61% accuracy and 92.46% recall for classifying nine cell subpopulations of healthy placentas. Biological landscapes of placenta heterogeneity could be mapped by the 110 marker genes screened by TURF_XGB, which revealed the superiority of the TURF feature mining. Moreover, we processed the PE dataset with LASSO to obtain 497 biomarkers. Integration analysis of the above two gene sets revealed that dendritic cells were closely associated with early-onset PE, and C1QB and C1QC might drive preeclampsia by mediating inflammation. In addition, an ensemble model-based risk stratification card was developed to classify preeclampsia patients, and its area under the receiver operating characteristic curve (AUC) could reach 0.99. For broader accessibility, we designed an accessible online web server ( http://bioinfor.imu.edu.cn/placenta ). CONCLUSION: Single-cell transcriptome-based preeclampsia risk assessment using an ensemble machine learning framework is a valuable asset for clinical decision-making. C1QB and C1QC may be involved in the development and progression of early-onset PE by affecting the complement and coagulation cascades pathway that mediate inflammation, which has important implications for better understanding the pathogenesis of PE.

EmAtlas: a comprehensive atlas for exploring spatiotemporal activation in mammalian embryogenesis.

The emerging importance of embryonic development research rapidly increases the volume for a professional resource related to multi-omics data. However, the lack of global embryogenesis repository and systematic analysis tools limits the preceding in stem cell research, human congenital diseases and assisted reproduction. Here, we developed the EmAtlas, which collects the most comprehensive multi-omics data and provides multi-scale tools to explore spatiotemporal activation during mammalian embryogenesis. EmAtlas contains data on multiple types of gene expression, chromatin accessibility, DNA methylation, nucleosome occupancy, histone modifications, and transcription factors, which displays the complete spatiotemporal landscape in mouse and human across several time points, involving gametogenesis, preimplantation, even fetus and neonate, and each tissue involves various cell types. To characterize signatures involved in the tissue, cell, genome, gene and protein levels during mammalian embryogenesis, analysis tools on these five scales were developed. Additionally, we proposed EmRanger to deliver extensive development-related biological background annotations. Users can utilize these tools to analyze, browse, visualize, and download data owing to the user-friendly interface. EmAtlas is freely accessible at http://bioinfor.imu.edu.cn/ematlas.

Temporal Dynamic Analysis of Alternative Splicing During Embryonic Development in Zebrafish.

Alternative splicing is pervasive in mammalian genomes and involved in embryo development, whereas research on crosstalk of alternative splicing and embryo development was largely restricted to mouse and human and the alternative splicing regulation during embryogenesis in zebrafish remained unclear. We constructed the alternative splicing atlas at 18 time-course stages covering maternal-to-zygotic transition, gastrulation, somitogenesis, pharyngula stages, and post-fertilization in zebrafish. The differential alternative splicing events between different developmental stages were detected. The results indicated that abundance alternative splicing and differential alternative splicing events are dynamically changed and remarkably abundant during the maternal-to-zygotic transition process. Based on gene expression profiles, we found splicing factors are expressed with specificity of developmental stage and largely expressed during the maternal-to-zygotic transition process. The better performance of cluster analysis was achieved based on the inclusion level of alternative splicing. The biological function analysis uncovered the important roles of alternative splicing during embryogenesis. The identification of isoform switches of alternative splicing provided a new insight into mining the regulated mechanism of transcript isoforms, which always is hidden by gene expression. In conclusion, we inferred that alternative splicing activation is synchronized with zygotic genome activation and discovered that alternative splicing is coupled with transcription during embryo development in zebrafish. We also unveiled that the temporal expression dynamics of splicing factors during embryo development, especially co-orthologous splicing factors. Furthermore, we proposed that the inclusion level of alternative splicing events can be employed for cluster analysis as a novel parameter. This work will provide a deeper insight into the regulation of alternative splicing during embryogenesis in zebrafish.

RaacFold: a webserver for 3D visualization and analysis of protein structure by using reduced amino acid alphabets.

Protein structure exhibits greater complexity and diversity than DNA structure, and usually affects the interpretation of the function, interactions and biological annotations. Reduced amino acid alphabets (Raaa) exhibit a powerful ability to decrease protein complexity and identify functional conserved regions, which motivated us to create RaacFold. The RaacFold provides 687 reduced amino acid clusters (Raac) based on 58 reduction methods and offers three analysis tools: Protein Analysis, Align Analysis, and Multi Analysis. The Protein Analysis and Align Analysis provide reduced representations of sequence-structure according to physicochemical similarities and computational biology strategies. With the simplified representations, the protein structure can be viewed more concise and clearer to capture biological insight than the unreduced structure. Thus, the design of artificial protein will be more convenient, and redundant interference is avoided. In addition, Multi Analysis allows users to explore biophysical variation and conservation in the evolution of protein structure and function. This supplies important information for the identification and exploration of the nonhomologous functions of paralogs. Simultaneously, RaacFold provides powerful 2D and 3D rendering performance with advanced parameters for sequences, structures, and related annotations. RaacFold is freely available at http://bioinfor.imu.edu.cn/raacfold.

Nucleosome Assembly and Disassembly in vitro Are Governed by Chemical Kinetic Principles.

As the elementary unit of eukaryotic chromatin, nucleosomes in vivo are highly dynamic in many biological processes, such as DNA replication, repair, recombination, or transcription, to allow the necessary factors to gain access to their substrate. The dynamic mechanism of nucleosome assembly and disassembly has not been well described thus far. We proposed a chemical kinetic model of nucleosome assembly and disassembly in vitro. In the model, the efficiency of nucleosome assembly was positively correlated with the total concentration of histone octamer, reaction rate constant and reaction time. All the corollaries of the model were well verified for the Widom 601 sequence and the six artificially synthesized DNA sequences, named CS1-CS6, by using the salt dialysis method in vitro. The reaction rate constant in the model may be used as a new parameter to evaluate the nucleosome reconstitution ability with DNAs. Nucleosome disassembly experiments for the Widom 601 sequence detected by Förster resonance energy transfer (FRET) and fluorescence thermal shift (FTS) assays demonstrated that nucleosome disassembly is the inverse process of assembly and can be described as three distinct stages: opening phase of the (H2A-H2B) dimer/(H3-H4)2 tetramer interface, release phase of the H2A-H2B dimers from (H3-H4)2 tetramer/DNA and removal phase of the (H3-H4)2 tetramer from DNA. Our kinetic model of nucleosome assembly and disassembly allows to confirm that nucleosome assembly and disassembly in vitro are governed by chemical kinetic principles.

A Genome-Wide Analysis of the Gene Expression and Alternative Splicing Events in a Whole-Body Hypoxic Preconditioning Mouse Model.

Exposure to specific doses of hypoxia can trigger endogenous neuroprotective and neuroplastic mechanisms of the central nervous system. These molecular mechanisms, together referred to as hypoxic preconditioning (HPC), remain poorly understood. In the present study, we applied RNA sequencing and bioinformatics analyses to study HPC in a whole-body HPC mouse model. The preconditioned (H4) and control (H0) groups showed 605 differentially expressed genes (DEGs), of which 263 were upregulated and 342 were downregulated. Gene Ontology enrichment analysis indicated that these DEGs were enriched in several biological processes, including metabolic stress and angiogenesis. The Kyoto Encyclopedia of Genes and Genomes enrichment analysis showed that the FOXO and Notch signaling pathways were involved in hypoxic tolerance and protection during HPC. Furthermore, 117 differential alternative splicing events (DASEs) were identified, with exon skipping being the dominant one (48.51%). Repeated exposure to systemic hypoxia promoted skipping of exon 7 in Edrf1 and exon 9 or 13 in Lrrc45. This study expands the understanding of the endogenous protective mechanisms of HPC and the DASEs that occur during HPC.

A deformation energy model reveals sequence-dependent property of nucleosome positioning.

We present a deformation energy model for predicting nucleosome positioning, in which a position-dependent structural parameter set derived from crystal structures of nucleosomes was used to calculate the DNA deformation energy. The model is successful in predicting nucleosome occupancy genome-wide in budding yeast, nucleosome free energy, and rotational positioning of nucleosomes. Our model also indicates that the genomic regions underlying the MNase-sensitive nucleosomes in budding yeast have high deformation energy and, consequently, low nucleosome-forming ability, while the MNase-sensitive non-histone particles are characterized by much lower DNA deformation energy and high nucleosome preference. In addition, we also revealed that remodelers, SNF2 and RSC8, are likely to act in chromatin remodeling by binding to broad nucleosome-depleted regions that are intrinsically favorable for nucleosome positioning. Our data support the important role of position-dependent physical properties of DNA in nucleosome positioning.

Transcriptome profile analysis reveals a silica-induced immune response and fibrosis in a silicosis rat model.

Silicosis is a type of pneumoconiosis caused by the inhalation of silica dust. It is characterized by inflammation and fibrosis of the lung. Although many studies have reported that crystalline silica-inhalation into the lung initiates the immune response, activating effector cells and triggering the inflammatory cascade with subsequent elaboration of the extracellular matrix and fibrosis, the mechanism of silicosis pathogenesis remains unclear. In the present study, we established a silica inhalation-induced silicosis rat model validated by histological and cytokine analyses. RNA-seq and bioinformatic analyses showed that 600 genes were upregulated and 537 genes were downregulated in the silica-treated group. GO enrichment analysis indicates that these differentially expressed genes are enriched in several biological processes including immune response and organism remodeling. KEGG enrichment analysis showed that 53 enriched pathways were mainly associated with human diseases, immune response, signal transduction, and fibrosis process. Since alternative splicing of pre-mRNAs is also essential for the regulation of gene expression, we identified several alternative pre-mRNA splicing events in the fibrotic process. This study will provide a foundation to understand the molecular mechanism of the pulmonary fibrosis caused by silica.

Dynamic Alternative Splicing During Mouse Preimplantation Embryo Development.

The mechanism of alternative pre-mRNA splicing (AS) during preimplantation development is largely unknown. In order to capture the dynamic changes of AS occurring during embryogenesis, we carried out bioinformatics analysis based on scRNA-seq data over the time-course preimplantation development in mouse. We detected numerous previously-unreported differentially expressed genes at specific developmental stages and investigated the nature of AS at both minor and major zygotic genome activation (ZGA). The AS and differential AS atlas over preimplantation development were established. The differentially alternatively spliced genes (DASGs) are likely to be key splicing factors (SFs) during preimplantation development. We also demonstrated that there is a regulatory cascade of AS events in which some key SFs are regulated by differentially AS of their own gene transcripts. Moreover, 212 isoform switches (ISs) during preimplantation development were detected, which may be critical for decoding the mechanism of early embryogenesis. Importantly, we uncovered that zygotic AS activation (ZASA) is in conformity with ZGA and revealed that AS is coupled with transcription during preimplantation development. Our results may provide a deeper insight into the regulation of early embryogenesis.

Deform-nu: A DNA Deformation Energy-Based Predictor for Nucleosome Positioning.

The structure and function of chromatin can be regulated through positioning patterns of nucleosomes. DNA-based processes are regulated via nucleosomes. Therefore, it is significant to determine nucleosome positions in DNA-based processes. A deformation energy model was proposed to predict nucleosome positions in our previous study. A free web server based on the model (http://lin-group.cn/server/deform-nu/) was firstly established to estimate the occupancy and rotational positioning of nucleosomes in the study. Then, the performance of the model was verified by several examples. The results indicated that nucleosome positioning relied on the physical properties of DNA, such as deformation energy.

The affinity of DNA sequences containing R5Y5 motif and TA repeats with 10.5-bp periodicity to histone octamer in vitro.

Nucleosome positioning along the genome is partially determined by the intrinsic DNA sequence preferences on histone. RRRRRYYYYY (R5Y5, R = Purine and Y = Pyrimidine) motif in nucleosome DNA, which was presented based on several theoretical models by Trifonov et al., might be a facilitating sequence pattern for nucleosome assembly. However, there is not a high conformity experimental evidence to support the concept that R5Y5 motif is a key element for the determination of nucleosome positioning. In this work, the ability of the canonical, H2A.Z- and H3.3-containing octamers to assemble nucleosome on DNA templates containing R5Y5 motif and TA repeats within 10.5-bp periodicity was investigated by using salt-dialysis method in vitro. The results showed that the10.5-bp periodical distributions of both R5Y5 motif and TA repeats along DNA templates can significantly promote canonical nucleosome assembly and may be key sequence factors for canonical nucleosome assembly. Compared with TA repeats within 10.5-bp periodicity, R5Y5 motif in DNA templates did not elevate H2A.Z- and H3.3-containing nucleosome formation efficiency in vitro. This result indicates that R5Y5 motif probably isn't a pivotal factor to regulate nucleosome assembly on histone variants. It is speculated that the regulatory mechanism of nucleosome assembly is different between canonical and variant histone. These conclusions can provide a deeper insight on the mechanism of nucleosome positioning. Communicated by Ramaswamy H. Sarma.

The implication of DNA bending energy for nucleosome positioning and sliding.

Nucleosome not only directly affects cellular processes, such as DNA replication, recombination, and transcription, but also severs as a fundamentally important target of epigenetic modifications. Our previous study indicated that the bending property of DNA is important in nucleosome formation, particularly in predicting the dyad positions of nucleosomes on a DNA segment. Here, we investigated the role of bending energy in nucleosome positioning and sliding in depth to decipher sequence-directed mechanism. The results show that bending energy is a good physical index to predict the free energy in the process of nucleosome reconstitution in vitro. Our data also imply that there are at least 20% of the nucleosomes in budding yeast do not adopt canonical positioning, in which underlying sequences wrapped around histones are structurally symmetric. We also revealed distinct patterns of bending energy profile for distinctly organized chromatin structures, such as well-positioned nucleosomes, fuzzy nucleosomes, and linker regions and discussed nucleosome sliding in terms of bending energy. We proposed that the stability of a nucleosome is positively correlated with the strength of the bending anisotropy of DNA segment, and both accessibility and directionality of nucleosome sliding is likely to be modulated by diverse patterns of DNA bending energy profile.

MiasDB: A Database of Molecular Interactions Associated with Alternative Splicing of Human Pre-mRNAs.

Alternative splicing (AS) is pervasive in human multi-exon genes and is a major contributor to expansion of the transcriptome and proteome diversity. The accurate recognition of alternative splice sites is regulated by information contained in networks of protein-protein and protein-RNA interactions. However, the mechanisms leading to splice site selection are not fully understood. Although numerous databases have been built to describe AS, molecular interaction databases associated with AS have only recently emerged. In this study, we present a new database, MiasDB, that provides a description of molecular interactions associated with human AS events. This database covers 938 interactions between human splicing factors, RNA elements, transcription factors, kinases and modified histones for 173 human AS events. Every entry includes the interaction partners, interaction type, experimental methods, AS type, tissue specificity or disease-relevant information, a simple description of the functionally tested interaction in the AS event and references. The database can be queried easily using a web server (http://47.88.84.236/Miasdb). We display some interaction figures for several genes. With this database, users can view the regulation network describing AS events for 12 given genes.

A deformation energy-based model for predicting nucleosome dyads and occupancy.

Nucleosome plays an essential role in various cellular processes, such as DNA replication, recombination, and transcription. Hence, it is important to decode the mechanism of nucleosome positioning and identify nucleosome positions in the genome. In this paper, we present a model for predicting nucleosome positioning based on DNA deformation, in which both bending and shearing of the nucleosomal DNA are considered. The model successfully predicted the dyad positions of nucleosomes assembled in vitro and the in vitro map of nucleosomes in Saccharomyces cerevisiae. Applying the model to Caenorhabditis elegans and Drosophila melanogaster, we achieved satisfactory results. Our data also show that shearing energy of nucleosomal DNA outperforms bending energy in nucleosome occupancy prediction and the ability to predict nucleosome dyad positions is attributed to bending energy that is associated with rotational positioning of nucleosomes.

Influence factors on the correlations between expression levels of neighboring pattern genes.

Some genes tend to cluster and be co-expressed. Multiple factors affect gene co-expression. In this study, we investigated the relationships between multiple factors and the correlations of expression levels of neighboring genes, which were divided into four kinds of pattern genes and one type of non-pattern gene. Our results indicate that the correlation between expression levels of neighboring non-pattern genes is related to multiple factors with the exception of transcriptional orientations of neighboring genes. The correlation between expression levels of neighboring specific genes or neighboring repressed genes is likely to be dependent on the co-functions of neighboring genes. The correlation between expression levels of neighboring housekeeping genes is associated with histone modifications in intergenic regions.

Using weighted features to predict recombination hotspots in Saccharomyces cerevisiae.

Characterization and accurate prediction of recombination hotspots and coldspots have crucial implications for the mechanism of recombination. Several models have predicted recombination hot/cold spots successfully, but there is still much room for improvement. We present a novel classifier in which k-mer frequency, physical and thermodynamic properties of DNA sequences are incorporated in the form of weighted features. Applying the classifier to recombination hot/cold ORFs in Saccharomyces cerevisiae, we achieved an accuracy of 90%, which is ~5% higher than existing methods, such as iRSpot-PseDNC, IDQD and Random Forest. The model also predicted non-ORF recombination hot/cold spots sequences in S. cerevisiae with high accuracy. A broad applicability of the model in the field of classification is expected.

GAA triplet-repeats cause nucleosome depletion in the human genome.

Although there have been many investigations into how trinucleotide repeats affect nucleosome formation and local chromatin structure, the nucleosome positioning of GAA triplet-repeats in the human genome has remained elusive. In this work, the nucleosome occupancy around GAA triplet-repeats across the human genome was computed statistically. The results showed a nucleosome-depleted region in the vicinity of GAA triplet-repeats in activated and resting CD4(+) T cells. Furthermore, the A-tract was frequently adjacent to the upstream region of GAA triplet-repeats and could enhance the depletion surrounding GAA triplet-repeats. In vitro chromatin reconstitution assays with GAA-containing plasmids also demonstrated that the inserted GAA triplet-repeats destabilized the ability of recombinant plasmids to assemble nucleosomes. Our results suggested that GAA triplet-repeats have lower affinity to histones and can change local nucleosome positioning. These findings may be helpful for understanding the mechanism of Friedreich's ataxia, which is associated with GAA triplet-repeats at the chromatin level.

Genome-wide characterization and prediction of Arabidopsis thaliana replication origins.

Identification of replication origins is crucial for the faithful duplication of genomic DNA. The frequencies of single nucleotides and dinucleotides, GC/AT bias and GC/AT profile in the vicinity of Arabidopsis thaliana replication origins were analyzed in the present work. The guanine content or cytosine content is higher in origin of replication (Ori) than in non-Ori. The SS (S=G or C) dinucleotides are favoured in Ori whereas WW (W=A or T) dinucleotides are favoured in non-Ori. GC/AT bias and GC/AT profile in Ori are significantly different from that in non-Ori. Furthermore, by inputting DNA sequence features into support vector machine, we distinguished between the Ori and non-Ori regions in A. thaliana. The total prediction accuracy is about 69.5% as evaluated by the 10-fold cross-validation. This result suggested that apart from DNA sequence, deciphering the selection of replication origin must integrate many other factors including nucleosome positioning, DNA methylation, histone modification, etc. In addition, by comparing predictive performance we found that the predictive accuracy of SVM using sequence features on the context of WS language is significantly better than that of RY language. Furthermore, the same conclusion was also obtained in S. cerevisiae and D. melanogaster.