SUMMARY N6-methyl-adenosine (m6A) is the most abundant modification on messenger RNAs and is linked to human diseases, but its functions in mammalian development are poorly understood. Here we reveal the evolutionary conservation and function of m6A by mapping the m6A methylome in mouse and human embryonic stem cells. Thousands of messenger and long noncoding RNAs show conserved m6A modification, including transcripts encoding core pluripotency transcription factors. m6A is enriched over 3′ untranslated regions at defined sequence motifs, and marks unstable transcripts, including transcripts turned over upon differentiation. Genetic inactivation or depletion of mouse and human Mettl3, one of the m6A methylases, led to m6A erasure on select target genes, prolonged Nanog expression upon differentiation, and impaired ESC’s exit from self-renewal towards differentiation into several lineages in vitro and in vivo. Thus, m6A is a mark of transcriptome flexibility required for stem cells to differentiate to specific lineages.

Development of the AMOEBA (Atomic Multipole Optimized Energetics for Biomolecular Simulation) force field for proteins is presented. The current version (AMOEBA-2013) utilizes permanent electrostatic multipole moments through the quadrupole at each atom, and explicitly treats polarization effects in various chemical and physical environments. The atomic multipole electrostatic parameters for each amino acid residue type are derived from high-level gas phase quantum mechanical calculations via a consistent and extensible protocol. Molecular polarizability is modeled via a Thole-style damped interactive induction model based upon distributed atomic polarizabilities. Inter- and intramolecular polarization is treated in a consistent fashion via the Thole model. The intramolecular polarization model ensures transferability of electrostatic parameters among different conformations, as demonstrated by the agreement between QM and AMOEBA electrostatic potentials, and dipole moments of dipeptides. The backbone and side chain torsional parameters were determined by comparing to gas-phase QM (RI-TRIM MP2/CBS) conformational energies of dipeptides and to statistical distributions from the Protein Data Bank. Molecular dynamics simulations are reported for short peptides in explicit water to examine their conformational properties in solution. Overall the calculated conformational free energies and J-coupling constants are consistent with PDB statistics and experimental NMR results, respectively. In addition, the experimental crystal structures of a number of proteins are well maintained during molecular dynamics (MD) simulation. While further calculations are necessary to fully validate the force field, initial results suggest the AMOEBA polarizable multipole force field is able to describe the structure and energetics of peptides and proteins, in both gas-phase and solution environments.

In parallel to the genetic code for protein synthesis, a second layer of information is embedded in all RNA transcripts in the form of RNA structure. RNA structure influences practically every step in the gene expression program1. Yet the nature of most RNA structures or effects of sequence variation on structure are not known. Here we report the initial landscape and variation of RNA secondary structures (RSS) in a human family Trio, providing a comprehensive RSS map of human coding and noncoding RNAs. We identify unique RSS signatures that demarcate open reading frames, splicing junctions, and define authentic microRNA binding sites. Comparison of native deproteinized RNA isolated from cells versus refolded purified RNA suggests that the majority of the RSS information is encoded within RNA sequence. Over 1900 transcribed single nucleotide variants (~15% of all transcribed SNVs) alter local RNA structure. We discover simple sequence and spacing rules that determine the ability of point mutations to impact RSS. Selective depletion of RiboSNitches versus structurally synonymous variants at precise locations suggests selection for specific RNA shapes at thousands of sites, including 3’UTRs, binding sites of miRNAs and RNA binding proteins genome-wide. These results highlight the potentially broad contribution of RNA structure and its variation to gene regulation.

Pseudogenes are thought to be inactive gene sequences, but recent evidence of extensive pseudogene transcription raised the question of potential function. Here we discover and characterize the sets of mouse lncRNAs induced by inflammatory signaling via TNFα. TNFα regulates hundreds of lncRNAs, including 54 pseudogene lncRNAs, several of which show exquisitely selective expression in response to specific cytokines and microbial components in a NF-κB-dependent manner. Lethe, a pseudogene lncRNA, is selectively induced by proinflammatory cytokines via NF-κB or glucocorticoid receptor agonist, and functions in negative feedback signaling to NF-κB. Lethe interacts with NF-κB subunit RelA to inhibit RelA DNA binding and target gene activation. Lethe level decreases with organismal age, a physiological state associated with increased NF-κB activity. These findings suggest that expression of pseudogenes lncRNAs are actively regulated and constitute functional regulators of inflammatory signaling.DOI: http://dx.doi.org/10.7554/eLife.00762.001

SUMMARY Long noncoding RNAs (lncRNAs) are thought to be prevalent regulators of gene expression, but the consequences of lncRNA inactivation in vivo are mostly unknown. Here we show that targeted deletion of mouse Hotair lncRNA leads to de-repression of hundreds of genes, resulting in homeotic transformation of the spine and malformation of metacarpal-carpal bones. RNA-seq and conditional inactivation reveal an ongoing requirement of Hotair to repress HoxD genes and several imprinted loci such as Dlk1-Meg3 and Igf2-H19, without affecting imprinting choice. Hotair binds to both Polycomb repressive complex 2 that methylates histone H3 at lysine 27 (H3K27) and Lsd1 complex that demethylates histone H3 at lysine 4 (H3K4) in vivo. Hotair inactivation causes H3K4me3 gain and, to a lesser extent, H3K27me3 loss at target genes. These results reveal the function and mechanisms of Hotair lncRNA to enforce silent chromatin state at Hox and additional genes.

Ionogels offer great potential for diverse electric applications. However, it remains challenging to fabricate high-performance ionogels with both good mechanical strength and high conductivity. Here, a new kind of transparent ionogel with both good mechanical strength and high conductivity is designed via locking a kind of free ionic liquid (IL), i.e., 1-ethyl-3-methylimidazolium dicyanamide ([EMIm][DCA]), into charged poly(2-acrylamido-2-methyl-1-propanesulfonic acid) (PAMPS)-based double networks. On the one hand, the charged PAMPS double network provides good mechanical strength and excellent recovery property. On the other hand, the free [EMIm][DCA] locked in the charged double network through electrostatic interaction offers ionic conductivity as high as ≈1.7-2.4 S m at 25 °C. It is demonstrated that the designed ionogel can be successfully used for a flexible skin sensor even under harsh conditions. Considering the rationally designed chemical structures of ILs and the diversity of charged polymer networks, it is envisioned that this strategy can be extended to a broad range of polymer systems. Moreover, functional components such as conducting polymers, 0D nanoparticles, 1D nanowires, and 2D nanosheets can be introduced into the polymer systems to fabricate diverse novel ionogels with unique functions. It is believed that this design principle will provide a new opportunity to construct next-generation multifunctional ionogels.

An approach combining fluorescence-activated cell sorting and high-throughput DNA sequencing (FACS-seq) was employed to determine the efficiency of start codon recognition for all possible translation initiation sites (TIS) utilizing AUG start codons. Using FACS-seq, we measured translation from a genetic reporter library representing all 65,536 possible TIS sequences spanning the −6 to +5 positions. We found that the motif RYMRMVAUGGC enhanced start codon recognition and translation efficiency. However, dinucleotide interactions, which cannot be conveyed by a single motif, were also important for modeling TIS efficiency. Our dataset combined with modeling allowed us to predict genome-wide translation initiation efficiency for all mRNA transcripts. Additionally, we screened somatic TIS mutations associated with tumorigenesis to identify candidate driver mutations consistent with known tumor expression patterns. Finally, we implemented a quantitative leaky scanning model to predict alternative initiation sites that produce truncated protein isoforms and compared predictions with ribosome footprint profiling data. The comprehensive analysis of the TIS sequence space enables quantitative predictions of translation initiation based on genome sequence.

Animals respond to environmental threats, e.g. looming visual stimuli, with innate defensive behaviors such as escape and freezing. The key neural circuits that participate in the generation of such dimorphic defensive behaviors remain unclear. Here we show that the dimorphic behavioral patterns triggered by looming visual stimuli are mediated by parvalbumin-positive (PV+) projection neurons in mouse superior colliculus (SC). Two distinct groups of SC PV+ neurons form divergent pathways to transmit threat-relevant visual signals to neurons in the parabigeminal nucleus (PBGN) and lateral posterior thalamic nucleus (LPTN). Activations of PV+ SC-PBGN and SC-LPTN pathways mimic the dimorphic defensive behaviors. The PBGN and LPTN neurons are co-activated by looming visual stimuli. Bilateral inactivation of either nucleus results in the defensive behavior dominated by the other nucleus. Together, these data suggest that the SC orchestrates dimorphic defensive behaviors through two separate tectofugal pathways that may have interactions.