SOX2 O-GlcNAcylation alters its protein-protein interactions and genomic occupancy to modulate gene expression in pluripotent cells.

The transcription factor SOX2 is central in establishing and maintaining pluripotency. The processes that modulate SOX2 activity to promote pluripotency are not well understood. Here, we show SOX2 is O-GlcNAc modified in its transactivation domain during reprogramming and in mouse embryonic stem cells (mESCs). Upon induction of differentiation SOX2 O-GlcNAcylation at serine 248 is decreased. Replacing wild type with an O-GlcNAc-deficient SOX2 (S248A) increases reprogramming efficiency. ESCs with O-GlcNAc-deficient SOX2 exhibit alterations in gene expression. This change correlates with altered protein-protein interactions and genomic occupancy of the O-GlcNAc-deficient SOX2 compared to wild type. In addition, SOX2 O-GlcNAcylation impairs the SOX2-PARP1 interaction, which has been shown to regulate ESC self-renewal. These findings show that SOX2 activity is modulated by O-GlcNAc, and provide a novel regulatory mechanism for this crucial pluripotency transcription factor.

Bat Accelerated Regions Identify a Bat Forelimb Specific Enhancer in the HoxD Locus.

The molecular events leading to the development of the bat wing remain largely unknown, and are thought to be caused, in part, by changes in gene expression during limb development. These expression changes could be instigated by variations in gene regulatory enhancers. Here, we used a comparative genomics approach to identify regions that evolved rapidly in the bat ancestor, but are highly conserved in other vertebrates. We discovered 166 bat accelerated regions (BARs) that overlap H3K27ac and p300 ChIP-seq peaks in developing mouse limbs. Using a mouse enhancer assay, we show that five Myotis lucifugus BARs drive gene expression in the developing mouse limb, with the majority showing differential enhancer activity compared to the mouse orthologous BAR sequences. These include BAR116, which is located telomeric to the HoxD cluster and had robust forelimb expression for the M. lucifugus sequence and no activity for the mouse sequence at embryonic day 12.5. Developing limb expression analysis of Hoxd10-Hoxd13 in Miniopterus natalensis bats showed a high-forelimb weak-hindlimb expression for Hoxd10-Hoxd11, similar to the expression trend observed for M. lucifugus BAR116 in mice, suggesting that it could be involved in the regulation of the bat HoxD complex. Combined, our results highlight novel regulatory regions that could be instrumental for the morphological differences leading to the development of the bat wing.

Transcriptomic and epigenomic characterization of the developing bat wing.

Bats are the only mammals capable of powered flight, but little is known about the genetic determinants that shape their wings. Here we generated a genome for Miniopterus natalensis and performed RNA-seq and ChIP-seq (H3K27ac and H3K27me3) analyses on its developing forelimb and hindlimb autopods at sequential embryonic stages to decipher the molecular events that underlie bat wing development. Over 7,000 genes and several long noncoding RNAs, including Tbx5-as1 and Hottip, were differentially expressed between forelimb and hindlimb, and across different stages. ChIP-seq analysis identified thousands of regions that are differentially modified in forelimb and hindlimb. Comparative genomics found 2,796 bat-accelerated regions within H3K27ac peaks, several of which cluster near limb-associated genes. Pathway analyses highlighted multiple ribosomal proteins and known limb patterning signaling pathways as differentially regulated and implicated increased forelimb mesenchymal condensation in differential growth. In combination, our work outlines multiple genetic components that likely contribute to bat wing formation, providing insights into this morphological innovation.

motifDiverge: a model for assessing the statistical significance of gene regulatory motif divergence between two DNA sequences.

Next-generation sequencing technology enables the identification of thousands of gene regulatory sequences in many cell types and organisms. We consider the problem of testing if two such sequences differ in their number of binding site motifs for a given transcription factor (TF) protein. Binding site motifs impart regulatory function by providing TFs the opportunity to bind to genomic elements and thereby affect the expression of nearby genes. Evolutionary changes to such functional DNA are hypothesized to be major contributors to phenotypic diversity within and between species; but despite the importance of TF motifs for gene expression, no method exists to test for motif loss or gain. Assuming that motif counts are Binomially distributed, and allowing for dependencies between motif instances in evolutionarily related sequences, we derive the probability mass function of the difference in motif counts between two nucleotide sequences. We provide a method to numerically estimate this distribution from genomic data and show through simulations that our estimator is accurate. Finally, we introduce the R package motifDiverge that implements our methodology and illustrate its application to gene regulatory enhancers identified by a mouse developmental time course experiment. While this study was motivated by analysis of regulatory motifs, our results can be applied to any problem involving two correlated Bernoulli trials.

Human-specific transcriptional networks in the brain.

Understanding human-specific patterns of brain gene expression and regulation can provide key insights into human brain evolution and speciation. Here, we use next-generation sequencing, and Illumina and Affymetrix microarray platforms, to compare the transcriptome of human, chimpanzee, and macaque telencephalon. Our analysis reveals a predominance of genes differentially expressed within human frontal lobe and a striking increase in transcriptional complexity specific to the human lineage in the frontal lobe. In contrast, caudate nucleus gene expression is highly conserved. We also identify gene coexpression signatures related to either neuronal processes or neuropsychiatric diseases, including a human-specific module with CLOCK as its hub gene and another module enriched for neuronal morphological processes and genes coexpressed with FOXP2, a gene important for language evolution. These data demonstrate that transcriptional networks have undergone evolutionary remodeling even within a given brain region, providing a window through which to view the foundation of uniquely human cognitive capacities.

RBFOX1 regulates both splicing and transcriptional networks in human neuronal development.

RNA splicing plays a critical role in the programming of neuronal differentiation and, consequently, normal human neurodevelopment, and its disruption may underlie neurodevelopmental and neuropsychiatric disorders. The RNA-binding protein, fox-1 homolog (RBFOX1; also termed A2BP1 or FOX1), is a neuron-specific splicing factor predicted to regulate neuronal splicing networks clinically implicated in neurodevelopmental disease, including autism spectrum disorder (ASD), but only a few targets have been experimentally identified. We used RNA sequencing to identify the RBFOX1 splicing network at a genome-wide level in primary human neural stem cells during differentiation. We observe that RBFOX1 regulates a wide range of alternative splicing events implicated in neuronal development and maturation, including transcription factors, other splicing factors and synaptic proteins. Downstream alterations in gene expression define an additional transcriptional network regulated by RBFOX1 involved in neurodevelopmental pathways remarkably parallel to those affected by splicing. Several of these differentially expressed genes are further implicated in ASD and related neurodevelopmental diseases. Weighted gene co-expression network analysis demonstrates a high degree of connectivity among these disease-related genes, highlighting RBFOX1 as a key factor coordinating the regulation of both neurodevelopmentally important alternative splicing events and clinically relevant neuronal transcriptional programs in the development of human neurons.

A geometric constraint, the head-to-tail exclusion rule, may be the basis for the isolated-pentagon rule in fullerenes with more than 60 vertices.

Carbon atoms self-assemble into the famous soccer-ball shaped Buckminsterfullerene (C(60)), the smallest fullerene cage that obeys the isolated-pentagon rule (IPR). Carbon atoms self-assemble into larger (n > 60 vertices) empty cages as well-but only the few that obey the IPR-and at least 1 small fullerene (n 60, the head-to-tail exclusion rule permits only (and all) fullerene cages and nanotubes that obey the IPR. We therefore suggest that self-assembly that obeys the IPR may be explained by the head-to-tail exclusion rule, a geometric constraint.

The physical basis for the head-to-tail rule that excludes most fullerene cages from self-assembly.

In the companion article, we proposed that fullerene cages with head-to-tail dihedral angle discrepancies do not self-assemble. Here we show why. If an edge abuts a pentagon at one end and a hexagon at the other, the dihedral angle about the edge increases, producing a dihedral angle discrepancy (DAD) vector. The DADs about all five/six edges of a central pentagonal/hexagonal face are determined by the identities-pentagon or hexagon-of its five/six surrounding faces. Each "Ring"-central face plus specific surrounding faces-may have zero, two, or four edges with DAD. In most Rings, the nonplanarity induced by DADs is shared among surrounding faces. However, in a Ring that has DADs arranged head of one to tail of another, the nonplanarity cannot be shared, so some surrounding faces would be especially nonplanar. Because the head-to-tail exclusion rule is an implicit geometric constraint, the rule may operate either by imposing a kinetic barrier that prevents assembly of certain Rings or by imposing an energy cost that makes those Rings unlikely to last in an equilibrium circumstance. Since Rings with head-to-tail DADs would be unlikely to self-assemble or last, fullerene cages with those Rings would be unlikely to self-assemble.

Gastrointestinal leukocytoclastic vasculitis: an adverse effect of sirolimus.

An 18-yr-old Hispanic female with end-stage renal disease secondary to chronic glomerulonephritis of unknown etiology underwent cadaveric renal transplantation. She was placed on a steroid-free protocol with tacrolimus and mycophenolate mofetil (MMF) for maintenance immunosuppression. Approximately 8 months post-transplantation, MMF was replaced by sirolimus (SRL) because of persistent leukopenia. Four months after the initiation of SRL, the patient began to experience chronic, constant periumbilical abdominal pain in the absence of vomiting, diarrhea or melena. Esophagogastroduodenoscopy and CT scans revealed significant diffuse mucosal thickening of the antrum, duodenum, and jejunum; leukocytoclastic vasculitis was identified on antral biopsy. A repeat biopsy after reduction of sirolimus dose by 50% over 6 months showed mild chronic inflammation of stomach and duodenum with some improvement in abdominal pain. Discontinuation of SRL and replacement with low dose MMF resulted in complete resolution of pain and normalization of gastrointestinal anatomy by imaging studies within 2 months. In light of this report, drug-induced leukocytoclastic vasculitis caused by SRL should be considered in the differential diagnosis of chronic abdominal pain in a patient with organ transplantation.