Sea-ATI unravels novel vocabularies of plant active cistrome.

The cistrome consists of all cis-acting regulatory elements recognized by transcription factors (TFs). However, only a portion of the cistrome is active for TF binding in a specific tissue. Resolving the active cistrome in plants remains challenging. In this study, we report the assay sequential extraction assisted-active TF identification (sea-ATI), a low-input method that profiles the DNA sequences recognized by TFs in a target tissue. We applied sea-ATI to seven plant tissues to survey their active cistrome and generated 41 motif models, including 15 new models that represent previously unidentified cis-regulatory vocabularies. ATAC-seq and RNA-seq analyses confirmed the functionality of the cis-elements from the new models, in that they are actively bound in vivo, located near the transcription start site, and influence chromatin accessibility and transcription. Furthermore, comparing dimeric WRKY CREs between sea-ATI and DAP-seq libraries revealed that thermodynamics and genetic drifts cooperatively shaped their evolution. Notably, sea-ATI can identify not only positive but also negative regulatory cis-elements, thereby providing unique insights into the functional non-coding genome of plants.

Milletdb: a multi-omics database to accelerate the research of functional genomics and molecular breeding of millets.

Millets are a class of nutrient-rich coarse cereals with high resistance to abiotic stress; thus, they guarantee food security for people living in areas with extreme climatic conditions and provide stress-related genetic resources for other crops. However, no platform is available to provide a comprehensive and systematic multi-omics analysis for millets, which seriously hinders the mining of stress-related genes and the molecular breeding of millets. Here, a free, web-accessible, user-friendly millets multi-omics database platform (Milletdb, http://milletdb.novogene.com) has been developed. The Milletdb contains six millets and their one related species genomes, graph-based pan-genomics of pearl millet, and stress-related multi-omics data, which enable Milletdb to be the most complete millets multi-omics database available. We stored GWAS (genome-wide association study) results of 20 yield-related trait data obtained under three environmental conditions [field (no stress), early drought and late drought] for 2 years in the database, allowing users to identify stress-related genes that support yield improvement. Milletdb can simplify the functional genomics analysis of millets by providing users with 20 different tools (e.g., 'Gene mapping', 'Co-expression', 'KEGG/GO Enrichment' analysis, etc.). On the Milletdb platform, a gene PMA1G03779.1 was identified through 'GWAS', which has the potential to modulate yield and respond to different environmental stresses. Using the tools provided by Milletdb, we found that the stress-related PLATZs TFs (transcription factors) family expands in 87.5% of millet accessions and contributes to vegetative growth and abiotic stress responses. Milletdb can effectively serve researchers in the mining of key genes, genome editing and molecular breeding of millets.

Ectopic Expression of FvVND4c Promotes Secondary Cell Wall Thickening and Flavonoid Accumulation in Fragaria vesca.

Secondary cell wall (SCW) thickening has a significant effect on the growth and development of plants, as well as in the resistance to various biotic and abiotic stresses. Lignin accounts for the strength of SCW. It is synthesized through the phenylpropanoid pathway that also leads to flavonoid synthesis. The coupling strategies for lignin and flavonoid syntheses are diverse in plants. How their syntheses are balanced by transcriptional regulation in fleshy fruits is still unclear. The diploid strawberry (Fragaria vesca) is a model for fleshy fruits research due to its small genome and wide scope of genetic transformation. SCW thickening is regulated by a multilevel transcriptional regulatory network wherein vascular-related NAC domains (VNDs) act as key regulators. In this study, we systematically characterized VNDs in Fragaria vesca and explored their functions. The overexpression of FvVND4c in diploid strawberry fruits resulted in SCW thickening and fruit color changes accompanied with the accumulation of lignin and flavonoids. Genes related to these phenotypes were also induced upon FvVND4c overexpression. Among the induced genes, we found FvMYB46 to be a direct downstream regulator of FvVND4c. The overexpression of FvMYB46 resulted in similar phenotypes as FvVND4c, except for the color change. Transcriptomic analyses suggest that both FvVND4c and FvMYB46 act on phenylpropanoid and flavonoid biosynthesis pathways, and induce lignin synthesis for SCW. These results suggest that FvVND4c and FvMYB46 cooperatively regulate SCW thickening and flavonoid accumulation in Fragaria vesca.

Sequence determinants of human gene regulatory elements.

DNA can determine where and when genes are expressed, but the full set of sequence determinants that control gene expression is unknown. Here, we measured the transcriptional activity of DNA sequences that represent an ~100 times larger sequence space than the human genome using massively parallel reporter assays (MPRAs). Machine learning models revealed that transcription factors (TFs) generally act in an additive manner with weak grammar and that most enhancers increase expression from a promoter by a mechanism that does not appear to involve specific TF-TF interactions. The enhancers themselves can be classified into three types: classical, closed chromatin and chromatin dependent. We also show that few TFs are strongly active in a cell, with most activities being similar between cell types. Individual TFs can have multiple gene regulatory activities, including chromatin opening and enhancing, promoting and determining transcription start site (TSS) activity, consistent with the view that the TF binding motif is the key atomic unit of gene expression.

Binding specificities of human RNA-binding proteins toward structured and linear RNA sequences.

RNA-binding proteins (RBPs) regulate RNA metabolism at multiple levels by affecting splicing of nascent transcripts, RNA folding, base modification, transport, localization, translation, and stability. Despite their central role in RNA function, the RNA-binding specificities of most RBPs remain unknown or incompletely defined. To address this, we have assembled a genome-scale collection of RBPs and their RNA-binding domains (RBDs) and assessed their specificities using high-throughput RNA-SELEX (HTR-SELEX). Approximately 70% of RBPs for which we obtained a motif bound to short linear sequences, whereas ∼30% preferred structured motifs folding into stem-loops. We also found that many RBPs can bind to multiple distinctly different motifs. Analysis of the matches of the motifs in human genomic sequences suggested novel roles for many RBPs. We found that three cytoplasmic proteins-ZC3H12A, ZC3H12B, and ZC3H12C-bound to motifs resembling the splice donor sequence, suggesting that these proteins are involved in degradation of cytoplasmic viral and/or unspliced transcripts. Structural analysis revealed that the RNA motif was not bound by the conventional C3H1 RNA-binding domain of ZC3H12B. Instead, the RNA motif was bound by the ZC3H12B's PilT N terminus (PIN) RNase domain, revealing a potential mechanism by which unconventional RBDs containing active sites or molecule-binding pockets could interact with short, structured RNA molecules. Our collection containing 145 high-resolution binding specificity models for 86 RBPs is the largest systematic resource for the analysis of human RBPs and will greatly facilitate future analysis of the various biological roles of this important class of proteins.

Nucleosome-bound SOX2 and SOX11 structures elucidate pioneer factor function.

'Pioneer' transcription factors are required for stem-cell pluripotency, cell differentiation and cell reprogramming1,2. Pioneer factors can bind nucleosomal DNA to enable gene expression from regions of the genome with closed chromatin. SOX2 is a prominent pioneer factor that is essential for pluripotency and self-renewal of embryonic stem cells3. Here we report cryo-electron microscopy structures of the DNA-binding domains of SOX2 and its close homologue SOX11 bound to nucleosomes. The structures show that SOX factors can bind and locally distort DNA at superhelical location 2. The factors also facilitate detachment of terminal nucleosomal DNA from the histone octamer, which increases DNA accessibility. SOX-factor binding to the nucleosome can also lead to a repositioning of the N-terminal tail of histone H4 that includes residue lysine 16. We speculate that this repositioning is incompatible with higher-order nucleosome stacking, which involves contacts of the H4 tail with a neighbouring nucleosome. Our results indicate that pioneer transcription factors can use binding energy to initiate chromatin opening, and thereby facilitate nucleosome remodelling and subsequent transcription.

The interaction landscape between transcription factors and the nucleosome.

Nucleosomes cover most of the genome and are thought to be displaced by transcription factors in regions that direct gene expression. However, the modes of interaction between transcription factors and nucleosomal DNA remain largely unknown. Here we systematically explore interactions between the nucleosome and 220 transcription factors representing diverse structural families. Consistent with earlier observations, we find that the majority of the studied transcription factors have less access to nucleosomal DNA than to free DNA. The motifs recovered from transcription factors bound to nucleosomal and free DNA are generally similar. However, steric hindrance and scaffolding by the nucleosome result in specific positioning and orientation of the motifs. Many transcription factors preferentially bind close to the end of nucleosomal DNA, or to periodic positions on the solvent-exposed side of the DNA. In addition, several transcription factors usually bind to nucleosomal DNA in a particular orientation. Some transcription factors specifically interact with DNA located at the dyad position at which only one DNA gyre is wound, whereas other transcription factors prefer sites spanning two DNA gyres and bind specifically to each of them. Our work reveals notable differences in the binding of transcription factors to free and nucleosomal DNA, and uncovers a diverse interaction landscape between transcription factors and the nucleosome.

A protein activity assay to measure global transcription factor activity reveals determinants of chromatin accessibility.

No existing method to characterize transcription factor (TF) binding to DNA allows genome-wide measurement of all TF-binding activity in cells. Here we present a massively parallel protein activity assay, active TF identification (ATI), that measures the DNA-binding activity of all TFs in cell or tissue extracts. ATI is based on electrophoretic separation of protein-bound DNA sequences from a highly complex DNA library and subsequent mass-spectrometric identification of the DNA-bound proteins. We applied ATI to four mouse tissues and mouse embryonic stem cells and found that, in a given tissue or cell type, a small set of TFs, which bound to only ∼10 distinct motifs, displayed strong DNA-binding activity. Some of these TFs were found in all cell types, whereas others were specific TFs known to determine cell fate in the analyzed tissue or cell type. We also show that a small number of TFs determined the accessible chromatin landscape of a cell, suggesting that gene regulatory logic may be simpler than previously appreciated.

Two distinct DNA sequences recognized by transcription factors represent enthalpy and entropy optima.

Most transcription factors (TFs) can bind to a population of sequences closely related to a single optimal site. However, some TFs can bind to two distinct sequences that represent two local optima in the Gibbs free energy of binding (ΔG). To determine the molecular mechanism behind this effect, we solved the structures of human HOXB13 and CDX2 bound to their two optimal DNA sequences, CAATAAA and TCGTAAA. Thermodynamic analyses by isothermal titration calorimetry revealed that both sites were bound with similar ΔG. However, the interaction with the CAA sequence was driven by change in enthalpy (ΔH), whereas the TCG site was bound with similar affinity due to smaller loss of entropy (ΔS). This thermodynamic mechanism that leads to at least two local optima likely affects many macromolecular interactions, as ΔG depends on two partially independent variables ΔH and ΔS according to the central equation of thermodynamics, ΔG = ΔH - TΔS.

Effects of magnesium ions and water molecules on the structure of amorphous calcium carbonate: a molecular dynamics study.

Molecular dynamics simulations were conducted to elucidate the effects of Mg(2+) and H2O additives on the structure of amorphous calcium carbonate (ACC). New potential parameters for Mg(2+) ions were developed. The distribution function of the angle formed by three nearest-neighbor atoms was introduced to analyze the short-range local structure of ACC. The simulation indicated that ACC had a weakly ordered local structure resembling the local structure of a CaCO3 crystal. The local structure of pure ACC resembled that of vaterite. The formation of the vaterite-like local structure was hindered by Mg(2+) ions, whereas H2O molecules did not significantly influence the structure of ACC when the fraction of H2O molecules was low. However, when the fraction of H2O was high, the formation of a monohydrocalcite-like local structure was promoted. The effects of the additives on the structure of ACC were verified using the size of the additives and the interaction between the additives and CaCO3. The simulated structure of ACC was compared with the structure of CaCO3 crystals nucleated through the formation of ACC particles in real systems.

Tuning the stability of CaCO3 crystals with magnesium ions for the formation of aragonite thin films on organic polymer templates.

Thin-film growth of aragonite CaCO3 on annealed poly(vinyl alcohol) (PVA) matrices is induced by adding Mg(2+) into a supersaturated solution of CaCO3. Both the growth rate and surface morphology of the aragonite thin films depend upon the concentration of Mg(2+) in the mineralization solution. In the absence of PVA matrices, no thin films are formed, despite the presence of Mg(2+). Molecular dynamics simulation of the CaCO3 precursor suggests that the transition of amorphous calcium carbonate to crystals is suppressed in the presence of Mg(2+). The role for ionic additives in the crystallization of CaCO3 on organic templates obtained in this study may provide useful information for the development of functional hybrid materials.

Subcellular localization of N-deoxyribosyltransferase in Lactobacillus fermentum: cell surface association of an intracellular nucleotide metabolic enzyme.

N-deoxyribosyltransferases are essential enzymes in the nucleotide salvage pathway of lactobacilli. They catalyze the exchange between the purine or pyrimidine bases of 2'-deoxyribonucleosides and free pyrimidine or purine bases. In general, N-deoxyribosyltransferases are referred to as cytoplasmic enzymes, although there is no experimental evidence for this subcellular localization. In this work, the subcellular localization of N-deoxyribosyltransferase II (NTD) from Lactobacillus fermentum was examined by subcellular fractionation, transmission electron microscopy, and fluorescence microscopy. Our results indicate that L. fermentum NTD are distributed not only in the cytoplasm but also on the cell wall surface, and further studies showed that surface-attached NTD can be released into the culture broth and conventional buffers.

Molecular approaches to understand biomineralization of shell nacreous layer.

The nacreous layer of molluskan shells, which consists of highly oriented aragonitic crystals and an organic matrix (including chitin and proteins), is a product of biomineralization. This paper briefly introduces the recent research advances on nacre biomineralization of shells from bivalves and gastropods, which mainly focus on analysis of the micro- and nano-structure and components of shell nacreous layers, and investigations of the characteristics and functions of matrix proteins from nacre. Matrix proteins not only participate in construction of the organic nacre framework, but also control the nucleation and growth of aragonitic crystals, as well as determine the polymorph specificity of calcium carbonate in nacre. Moreover, the inorganic aragonite phase also plays an active role in organizing nacre microstructure. Based on these studies, several models to illustrate the formation mechanism related to lamellar nacre in bivalves, and columnar nacre in gastropods are introduced.

Cloning, characterization and immunolocalization of two subunits of calcineurin from pearl oyster (Pinctada fucata).

Calcineurin (CN), consisting of catalytic subunit (CN A) and regulatory subunit (CN B), is a multifunctional protein involved in many important physiological processes. Here, we cloned two subunits of CN (Pf-CN A and Pf-CN B) from pearl oyster Pinctada fucata and reported, for the first time, its expression patterns in the developmental stages, its enzymatic activity and immunolocalization in various tissues of adult pearl oyster. The Pf-CN A was extensively localized in all the tested tissues including mantle, gonad, digestive gland, gills, adductor muscle, and foot with strong signals detected in gonad, gills, foot, and mantle. Importantly, Pf-CN A was mainly found in the inner epithelial cells of the basal periostracal groove and lateral surface of the inner mantle fold, in which organic macromolecules used for periostracum formation and shell construction are secreted, respectively. In gill, the strong signals were distributed in the epithelial cells of the branchial filaments and the base of gill filaments. All the results suggested that Pf-CN may participate in the development of the pearl oyster and function in many ways in various physiological activities, especially in the shell formation. Our observations could provide some important clues to further understanding of the functions of CN in the oyster.

Cloning and characterization of Prisilkin-39, a novel matrix protein serving a dual role in the prismatic layer formation from the oyster Pinctada fucata.

Molluscs form their shells out of CaCO(3) and a matrix of biomacromolecules. Understanding the role of matrices may shed some light on the mechanism of biomineralization. Here, a 1401-bp full-length cDNA sequence encoding a novel matrix protein was cloned from the mantle of the bivalve oyster, Pinctada fucata. The deduced protein (Prisilkin-39), which has a molecular mass of 39.3 kDa and an isoelectric point of 8.83, was fully characterized, and its role in biomineralization was demonstrated using both in vivo and in vitro crystal growth assays. Prisilkin-39 is a highly repetitive protein with an unusual composition of Gly, Tyr, and Ser residues. Expression of Prisilkin-39 was localized to columnar epithelial cells of the mantle edge, corresponding to the calcitic prismatic layer formation. Immunostaining in situ and immunodetection in vitro revealed the presence of a characteristic pattern of Prisilkin-39 in the organic sheet and in sheaths around the prisms. Prisilkin-39 binds tightly with chitin, an insoluble polysaccharide that forms the highly structured framework of the shell. Antibody injection in vivo resulted in dramatic morphological deformities in the inner shell surface structure, where large amounts of CaCO(3) were deposited in an uncontrolled manner. Moreover, Prisilkin-39 strictly prohibited the precipitation of aragonite in vitro. Taken together, Prisilkin-39 is the first protein shown to have dual function, involved both in the chitinous framework building and in crystal growth regulation during the prismatic layer mineralization. These observations may extend our view on the rare group of basic matrices and their functions during elaboration of the molluscan shell.