ALBERT

Description: by Yufeng Huang, Chuchu Wang, Yufeng Yao, Xiaoyu Zuo, Shanshan Chen, Chengqi Xu, Hongfu Zhang, Qiulun Lu, Le Chang, Fan Wang, Pengxia Wang, Rongfeng Zhang, Zhenkun Hu, Qixue Song, Xiaowei Yang, Cong Li, Sisi Li, Yuanyuan Zhao, Qin Yang, Dan Yin, Xiaojing Wang, Wenxia Si, Xiuchun Li, Xin Xiong, Dan Wang, Yuan Huang, Chunyan Luo, Jia Li, Jingjing Wang, Jing Chen, Longfei Wang, Li Wang, Meng Han, Jian Ye, Feifei Chen, Jingqiu Liu, Ying Liu, Gang Wu, Bo Yang, Xiang Cheng, Yuhua Liao, Yanxia Wu, Tie Ke, Qiuyun Chen, Xin Tu, Robert Elston, Shaoqi Rao, Yanzong Yang, Yunlong Xia, Qing K. Wang Atrial fibrillation (AF) is the most common cardiac arrhythmia at the clinic. Recent GWAS identified several variants associated with AF, but they account for

Description: by Diane I. Schroeder, Kartika Jayashankar, Kory C. Douglas, Twanda L. Thirkill, Daniel York, Pete J. Dickinson, Lawrence E. Williams, Paul B. Samollow, Pablo J. Ross, Danika L. Bannasch, Gordon C. Douglas, Janine M. LaSalle Over the last 20-80 million years the mammalian placenta has taken on a variety of morphologies through both divergent and convergent evolution. Recently we have shown that the human placenta genome has a unique epigenetic pattern of large partially methylated domains (PMDs) and highly methylated domains (HMDs) with gene body DNA methylation positively correlating with level of gene expression. In order to determine the evolutionary conservation of DNA methylation patterns and transcriptional regulatory programs in the placenta, we performed a genome-wide methylome (MethylC-seq) analysis of human, rhesus macaque, squirrel monkey, mouse, dog, horse, and cow placentas as well as opossum extraembryonic membrane. We found that, similar to human placenta, mammalian placentas and opossum extraembryonic membrane have globally lower levels of methylation compared to somatic tissues. Higher relative gene body methylation was the conserved feature across all mammalian placentas, despite differences in PMD/HMDs and absolute methylation levels. Specifically, higher methylation over the bodies of genes involved in mitosis, vesicle-mediated transport, protein phosphorylation, and chromatin modification was observed compared with the rest of the genome. As in human placenta, higher methylation is associated with higher gene expression and is predictive of genic location across species. Analysis of DNA methylation in oocytes and preimplantation embryos shows a conserved pattern of gene body methylation similar to the placenta. Intriguingly, mouse and cow oocytes and mouse early embryos have PMD/HMDs but their placentas do not, suggesting that PMD/HMDs are a feature of early preimplantation methylation patterns that become lost during placental development in some species and following implantation of the embryo.

Description: by Bin Huang, Yongkui Wang, Wenhao Wang, Juan Chen, Pinglin Lai, Zhongyu Liu, Bo Yan, Song Xu, Zhongmin Zhang, Chun Zeng, Limin Rong, Bin Liu, Daozhang Cai, Dadi Jin, Xiaochun Bai The mechanistic target of rapamycin (mTOR) integrates both intracellular and extracellular signals to regulate cell growth and metabolism. However, the role of mTOR signaling in osteoblast differentiation and bone formation is undefined, and the underlying mechanisms have not been elucidated. Here, we report that activation of mTOR complex 1 (mTORC1) is required for preosteoblast proliferation; however, inactivation of mTORC1 is essential for their differentiation and maturation. Inhibition of mTORC1 prevented preosteoblast proliferation, but enhanced their differentiation in vitro and in mice. Activation of mTORC1 by deletion of tuberous sclerosis 1 ( Tsc1 ) in preosteoblasts produced immature woven bone in mice due to excess proliferation but impaired differentiation and maturation of the cells. The mTORC1-specific inhibitor, rapamycin, restored these in vitro and in vivo phenotypic changes. Mechanistically, mTORC1 prevented osteoblast maturation through activation of the STAT3/p63/Jagged/Notch pathway and downregulation of Runx2. Preosteoblasts with hyperactive mTORC1 reacquired the capacity to fully differentiate and maturate when subjected to inhibition of the Notch pathway. Together, these findings identified the role of mTORC1 in osteoblast formation and established that mTORC1 prevents preosteoblast differentiation and maturation through activation of the Notch pathway.

Description: by David M. Truong, F. Curtis Hewitt, Joseph H. Hanson, Xiaoxia Cui, Alan M. Lambowitz Mobile bacterial group II introns are evolutionary ancestors of spliceosomal introns and retroelements in eukaryotes. They consist of an autocatalytic intron RNA (a “ribozyme”) and an intron-encoded reverse transcriptase, which function together to promote intron integration into new DNA sites by a mechanism termed “retrohoming”. Although mobile group II introns splice and retrohome efficiently in bacteria, all examined thus far function inefficiently in eukaryotes, where their ribozyme activity is limited by low Mg 2+ concentrations, and intron-containing transcripts are subject to nonsense-mediated decay (NMD) and translational repression. Here, by using RNA polymerase II to express a humanized group II intron reverse transcriptase and T7 RNA polymerase to express intron transcripts resistant to NMD, we find that simply supplementing culture medium with Mg 2+ induces the Lactococcus lactis Ll.LtrB intron to retrohome into plasmid and chromosomal sites, the latter at frequencies up to ~0.1%, in viable HEK-293 cells. Surprisingly, under these conditions, the Ll.LtrB intron reverse transcriptase is required for retrohoming but not for RNA splicing as in bacteria. By using a genetic assay for in vivo selections combined with deep sequencing, we identified intron RNA mutations that enhance retrohoming in human cells, but

Description: by María F. Organista, Mercedes Martín, Jesus M. de Celis, Rosa Barrio, Ana López-Varea, Nuria Esteban, Mar Casado, Jose F. de Celis The Drosophila genes spalt major (salm) and spalt-related (salr) encode Zn-finger transcription factors regulated by the Decapentaplegic (Dpp) signalling pathway in the wing imaginal disc. The function of these genes is required for cell survival and proliferation in the central region of the wing disc, and also for vein patterning in the lateral regions. The identification of direct Salm and Salr target genes, and the analysis of their functions, are critical steps towards understanding the genetic control of growth and patterning of the Drosophila wing imaginal disc by the Dpp pathway. To identify candidate Salm/Salr target genes, we have compared the expression profile of salm/salr knockdown wing discs with control discs in microarray experiments. We studied by in situ hybridization the expression pattern of the genes whose mRNA levels varied significantly, and uncovered a complex transcription landscape regulated by the Spalt proteins in the wing disc. Interestingly, candidate Salm/Salr targets include genes which expression is turned off and genes which expression is positively regulated by Salm/Salr. Furthermore, loss-of-function phenotypic analysis of these genes indicates, for a fraction of them, a requirement for wing growth and patterning. The identification and analysis of candidate Salm/Salr target genes opens a new avenue to reconstruct the genetic structure of the wing, linking the activity of the Dpp pathway to the development of this epithelial tissue.

Description: by Alexandra A. Erwin, Mauricio A. Galdos, Michelle L. Wickersheim, Chris C. Harrison, Kendra D. Marr, Jack M. Colicchio, Justin P. Blumenstiel Sexual reproduction allows transposable elements (TEs) to proliferate, leading to rapid divergence between populations and species. A significant outcome of divergence in the TE landscape is evident in hybrid dysgenic syndromes, a strong form of genomic incompatibility that can arise when (TE) family abundance differs between two parents. When TEs inherited from the father are absent in the mother's genome, TEs can become activated in the progeny, causing germline damage and sterility. Studies in Drosophila indicate that dysgenesis can occur when TEs inherited paternally are not matched with a pool of corresponding TE silencing PIWI-interacting RNAs (piRNAs) provisioned by the female germline. Using the D . virilis syndrome of hybrid dysgenesis as a model, we characterize the effects that divergence in TE profile between parents has on offspring. Overall, we show that divergence in the TE landscape is associated with persisting differences in germline TE expression when comparing genetically identical females of reciprocal crosses and these differences are transmitted to the next generation. Moreover, chronic and persisting TE expression coincides with increased levels of genic piRNAs associated with reduced gene expression. Combined with these effects, we further demonstrate that gene expression is idiosyncratically influenced by differences in the genic piRNA profile of the parents that arise though polymorphic TE insertions. Overall, these results support a model in which early germline events in dysgenesis establish a chronic, stable state of both TE and gene expression in the germline that is maintained through adulthood and transmitted to the next generation. This work demonstrates that divergence in the TE profile is associated with diverse piRNA-mediated transgenerational effects on gene expression within populations.

Description: by Caitlin Sedwick

Description: by Tomer Stern, Rona Aviram, Chagai Rot, Tal Galili, Amnon Sharir, Noga Kalish Achrai, Yosi Keller, Ron Shahar, Elazar Zelzer One of the major challenges that developing organs face is scaling, that is, the adjustment of physical proportions during the massive increase in size. Although organ scaling is fundamental for development and function, little is known about the mechanisms that regulate it. Bone superstructures are projections that typically serve for tendon and ligament insertion or articulation and, therefore, their position along the bone is crucial for musculoskeletal functionality. As bones are rigid structures that elongate only from their ends, it is unclear how superstructure positions are regulated during growth to end up in the right locations. Here, we document the process of longitudinal scaling in developing mouse long bones and uncover the mechanism that regulates it. To that end, we performed a computational analysis of hundreds of three-dimensional micro-CT images, using a newly developed method for recovering the morphogenetic sequence of developing bones. Strikingly, analysis revealed that the relative position of all superstructures along the bone is highly preserved during more than a 5-fold increase in length, indicating isometric scaling. It has been suggested that during development, bone superstructures are continuously reconstructed and relocated along the shaft, a process known as drift. Surprisingly, our results showed that most superstructures did not drift at all. Instead, we identified a novel mechanism for bone scaling, whereby each bone exhibits a specific and unique balance between proximal and distal growth rates, which accurately maintains the relative position of its superstructures. Moreover, we show mathematically that this mechanism minimizes the cumulative drift of all superstructures, thereby optimizing the scaling process. Our study reveals a general mechanism for the scaling of developing bones. More broadly, these findings suggest an evolutionary mechanism that facilitates variability in bone morphology by controlling the activity of individual epiphyseal plates.

Description: by Hayato Yanagida, Ariel Gispan, Noam Kadouri, Shelly Rozen, Michal Sharon, Naama Barkai, Dan S. Tawfik Errors in protein synthesis, so-called phenotypic mutations, are orders-of-magnitude more frequent than genetic mutations. Here, we provide direct evidence that alternative protein forms and phenotypic variability derived from translational errors paved the path to genetic, evolutionary adaptations via gene duplication. We explored the evolutionary origins of Saccharomyces cerevisiae IDP3 - an NADP-dependent isocitrate dehydrogenase mediating fatty acids ß -oxidation in the peroxisome. Following the yeast whole genome duplication, IDP3 diverged from a cytosolic ancestral gene by acquisition of a C-terminal peroxisomal targeting signal. We discovered that the pre-duplicated cytosolic IDP s are partially localized to the peroxisome owing to +1 translational frameshifts that bypass the stop codon and unveil cryptic peroxisomal targeting signals within the 3’-UTR. Exploring putative cryptic signals in all 3’-UTRs of yeast genomes, we found that other enzymes related to NADPH production such as pyruvate carboxylase 1 ( PYC1 ) might be prone to peroxisomal localization via cryptic signals. Using laboratory evolution we found that these translational frameshifts are rapidly imprinted via genetic single base deletions occurring within the very same gene location. Further, as exemplified here, the sequences that promote translational frameshifts are also more prone to genetic deletions. Thus, genotypes conferring higher phenotypic variability not only meet immediate challenges by unveiling cryptic 3’-UTR sequences, but also boost the potential for future genetic adaptations.

Description: by Norishige Yamada, Yuko Hasegawa, Minghui Yue, Tomofumi Hamada, Shinichi Nakagawa, Yuya Ogawa To equalize X-linked gene dosage between the sexes in mammalian females, Xist RNA inactivates one of the two X-chromosomes. Here, we report the crucial function of Xist exon 7 in X-inactivation. Xist exon 7 is the second-largest exon with a well-conserved repeat E in eutherian mammals, but its role is often overlooked in X-inactivation. Although female ES cells with a targeted truncation of the Xist exon 7 showed no significant differences in their Xist expression levels and RNA stability from control cells expressing wild-type Xist , compromised localization of Xist RNA and incomplete silencing of X-linked genes on the inactive X-chromosome (Xi) were observed in the exon 7-truncated mutant cells. Furthermore, the interaction between the mutant Xist RNA and hnRNP U required for localization of Xist RNA to the Xi was impaired in the Xist exon 7 truncation mutant cells. Our results suggest that exon 7 of Xist RNA plays an important role for stable Xist RNA localization and silencing of the X-linked genes on the Xi, possibly acting through an interaction with hnRNP U.

Description: by Eliseo Ferrante, Ali Emre Turgut, Edgar Duéñez-Guzmán, Marco Dorigo, Tom Wenseleers Division of labor is ubiquitous in biological systems, as evidenced by various forms of complex task specialization observed in both animal societies and multicellular organisms. Although clearly adaptive, the way in which division of labor first evolved remains enigmatic, as it requires the simultaneous co-occurrence of several complex traits to achieve the required degree of coordination. Recently, evolutionary swarm robotics has emerged as an excellent test bed to study the evolution of coordinated group-level behavior. Here we use this framework for the first time to study the evolutionary origin of behavioral task specialization among groups of identical robots. The scenario we study involves an advanced form of division of labor, common in insect societies and known as “task partitioning”, whereby two sets of tasks have to be carried out in sequence by different individuals. Our results show that task partitioning is favored whenever the environment has features that, when exploited, reduce switching costs and increase the net efficiency of the group, and that an optimal mix of task specialists is achieved most readily when the behavioral repertoires aimed at carrying out the different subtasks are available as pre-adapted building blocks. Nevertheless, we also show for the first time that self-organized task specialization could be evolved entirely from scratch, starting only from basic, low-level behavioral primitives, using a nature-inspired evolutionary method known as Grammatical Evolution. Remarkably, division of labor was achieved merely by selecting on overall group performance, and without providing any prior information on how the global object retrieval task was best divided into smaller subtasks. We discuss the potential of our method for engineering adaptively behaving robot swarms and interpret our results in relation to the likely path that nature took to evolve complex sociality and task specialization.

Description: by Patrícia Santos-Oliveira, António Correia, Tiago Rodrigues, Teresa M Ribeiro-Rodrigues, Paulo Matafome, Juan Carlos Rodríguez-Manzaneque, Raquel Seiça, Henrique Girão, Rui D. M. Travasso Sprouting angiogenesis, where new blood vessels grow from pre-existing ones, is a complex process where biochemical and mechanical signals regulate endothelial cell proliferation and movement. Therefore, a mathematical description of sprouting angiogenesis has to take into consideration biological signals as well as relevant physical processes, in particular the mechanical interplay between adjacent endothelial cells and the extracellular microenvironment. In this work, we introduce the first phase-field continuous model of sprouting angiogenesis capable of predicting sprout morphology as a function of the elastic properties of the tissues and the traction forces exerted by the cells. The model is very compact, only consisting of three coupled partial differential equations, and has the clear advantage of a reduced number of parameters. This model allows us to describe sprout growth as a function of the cell-cell adhesion forces and the traction force exerted by the sprout tip cell. In the absence of proliferation, we observe that the sprout either achieves a maximum length or, when the traction and adhesion are very large, it breaks. Endothelial cell proliferation alters significantly sprout morphology, and we explore how different types of endothelial cell proliferation regulation are able to determine the shape of the growing sprout. The largest region in parameter space with well formed long and straight sprouts is obtained always when the proliferation is triggered by endothelial cell strain and its rate grows with angiogenic factor concentration. We conclude that in this scenario the tip cell has the role of creating a tension in the cells that follow its lead. On those first stalk cells, this tension produces strain and/or empty spaces, inevitably triggering cell proliferation. The new cells occupy the space behind the tip, the tension decreases, and the process restarts. Our results highlight the ability of mathematical models to suggest relevant hypotheses with respect to the role of forces in sprouting, hence underlining the necessary collaboration between modelling and molecular biology techniques to improve the current state-of-the-art.

Description: by Helen Budworth, Faye R. Harris, Paul Williams, Do Yup Lee, Amy Holt, Jens Pahnke, Bartosz Szczesny, Karina Acevedo-Torres, Sylvette Ayala-Peña, Cynthia T. McMurray Huntington’s Disease (HD) is caused by inheritance of a single disease-length allele harboring an expanded CAG repeat, which continues to expand in somatic tissues with age. The inherited disease allele expresses a toxic protein, and whether further somatic expansion adds to toxicity is unknown. We have created an HD mouse model that resolves the effects of the inherited and somatic expansions. We show here that suppressing somatic expansion substantially delays the onset of disease in littermates that inherit the same disease-length allele. Furthermore, a pharmacological inhibitor, XJB-5-131, inhibits the lengthening of the repeat tracks, and correlates with rescue of motor decline in these animals. The results provide evidence that pharmacological approaches to offset disease progression are possible.

Description: by Sayed-Rzgar Hosseini, Aditya Barve, Andreas Wagner All biological evolution takes place in a space of possible genotypes and their phenotypes. The structure of this space defines the evolutionary potential and limitations of an evolving system. Metabolism is one of the most ancient and fundamental evolving systems, sustaining life by extracting energy from extracellular nutrients. Here we study metabolism’s potential for innovation by analyzing an exhaustive genotype-phenotype map for a space of 10 15 metabolisms that encodes all possible subsets of 51 reactions in central carbon metabolism. Using flux balance analysis, we predict the viability of these metabolisms on 10 different carbon sources which give rise to 1024 potential metabolic phenotypes. Although viable metabolisms with any one phenotype comprise a tiny fraction of genotype space, their absolute numbers exceed 10 9 for some phenotypes. Metabolisms with any one phenotype typically form a single network of genotypes that extends far or all the way through metabolic genotype space, where any two genotypes can be reached from each other through a series of single reaction changes. The minimal distance of genotype networks associated with different phenotypes is small, such that one can reach metabolisms with novel phenotypes – viable on new carbon sources – through one or few genotypic changes. Exceptions to these principles exist for those metabolisms whose complexity (number of reactions) is close to the minimum needed for viability. Increasing metabolic complexity enhances the potential for both evolutionary conservation and evolutionary innovation.

Description: by Roland G. Roberts

Description: by Seth R. Bordenstein, Kevin R. Theis Groundbreaking research on the universality and diversity of microorganisms is now challenging the life sciences to upgrade fundamental theories that once seemed untouchable. To fully appreciate the change that the field is now undergoing, one has to place the epochs and foundational principles of Darwin, Mendel, and the modern synthesis in light of the current advances that are enabling a new vision for the central importance of microbiology. Animals and plants are no longer heralded as autonomous entities but rather as biomolecular networks composed of the host plus its associated microbes, i.e., "holobionts." As such, their collective genomes forge a "hologenome," and models of animal and plant biology that do not account for these intergenomic associations are incomplete. Here, we integrate these concepts into historical and contemporary visions of biology and summarize a predictive and refutable framework for their evaluation. Specifically, we present ten principles that clarify and append what these concepts are and are not, explain how they both support and extend existing theory in the life sciences, and discuss their potential ramifications for the multifaceted approaches of zoology and botany. We anticipate that the conceptual and evidence-based foundation provided in this essay will serve as a roadmap for hypothesis-driven, experimentally validated research on holobionts and their hologenomes, thereby catalyzing the continued fusion of biology's subdisciplines. At a time when symbiotic microbes are recognized as fundamental to all aspects of animal and plant biology, the holobiont and hologenome concepts afford a holistic view of biological complexity that is consistent with the generally reductionist approaches of biology.

Description: by Pengxing Cao, Ada W. C. Yan, Jane M. Heffernan, Stephen Petrie, Robert G. Moss, Louise A. Carolan, Teagan A. Guarnaccia, Anne Kelso, Ian G. Barr, Jodie McVernon, Karen L. Laurie, James M. McCaw Influenza is an infectious disease that primarily attacks the respiratory system. Innate immunity provides both a very early defense to influenza virus invasion and an effective control of viral growth. Previous modelling studies of virus–innate immune response interactions have focused on infection with a single virus and, while improving our understanding of viral and immune dynamics, have been unable to effectively evaluate the relative feasibility of different hypothesised mechanisms of antiviral immunity. In recent experiments, we have applied consecutive exposures to different virus strains in a ferret model, and demonstrated that viruses differed in their ability to induce a state of temporary immunity or viral interference capable of modifying the infection kinetics of the subsequent exposure. These results imply that virus-induced early immune responses may be responsible for the observed viral hierarchy. Here we introduce and analyse a family of within-host models of re-infection viral kinetics which allow for different viruses to stimulate the innate immune response to different degrees. The proposed models differ in their hypothesised mechanisms of action of the non-specific innate immune response. We compare these alternative models in terms of their abilities to reproduce the re-exposure data. Our results show that 1) a model with viral control mediated solely by a virus-resistant state, as commonly considered in the literature, is not able to reproduce the observed viral hierarchy; 2) the synchronised and desynchronised behaviour of consecutive virus infections is highly dependent upon the interval between primary virus and challenge virus exposures and is consistent with virus-dependent stimulation of the innate immune response. Our study provides the first mechanistic explanation for the recently observed influenza viral hierarchies and demonstrates the importance of understanding the host response to multi-strain viral infections. Re-exposure experiments provide a new paradigm in which to study the immune response to influenza and its role in viral control.

Description: by Paul M. Harrison, Laurent Badel, Mark J. Wall, Magnus J. E. Richardson Models of neocortical networks are increasingly including the diversity of excitatory and inhibitory neuronal classes. Significant variability in cellular properties are also seen within a nominal neuronal class and this heterogeneity can be expected to influence the population response and information processing in networks. Recent studies have examined the population and network effects of variability in a particular neuronal parameter with some plausibly chosen distribution. However, the empirical variability and covariance seen across multiple parameters are rarely included, partly due to the lack of data on parameter correlations in forms convenient for model construction. To addess this we quantify the heterogeneity within and between the neocortical pyramidal-cell classes in layers 2/3, 4, and the slender-tufted and thick-tufted pyramidal cells of layer 5 using a combination of intracellular recordings, single-neuron modelling and statistical analyses. From the response to both square-pulse and naturalistic fluctuating stimuli, we examined the class-dependent variance and covariance of electrophysiological parameters and identify the role of the h current in generating parameter correlations. A byproduct of the dynamic I-V method we employed is the straightforward extraction of reduced neuron models from experiment. Empirically these models took the refractory exponential integrate-and-fire form and provide an accurate fit to the perisomatic voltage responses of the diverse pyramidal-cell populations when the class-dependent statistics of the model parameters were respected. By quantifying the parameter statistics we obtained an algorithm which generates populations of model neurons, for each of the four pyramidal-cell classes, that adhere to experimentally observed marginal distributions and parameter correlations. As well as providing this tool, which we hope will be of use for exploring the effects of heterogeneity in neocortical networks, we also provide the code for the dynamic I-V method and make the full electrophysiological data set available.

Description: by Naomi Levy-Strumpf, Meghan Krizus, Hong Zheng, Louise Brown, Joseph G. Culotti Wnt and Netrin signaling regulate diverse essential functions. Using a genetic approach combined with temporal gene expression analysis, we found a regulatory link between the Wnt receptor MOM-5/Frizzled and the UNC-6/Netrin receptor UNC-5. These two receptors play key roles in guiding cell and axon migrations, including the migration of the C . elegans Distal Tip Cells (DTCs). DTCs migrate post-embryonically in three sequential phases: in the first phase along the Antero-Posterior (A/P) axis, in the second, along the Dorso-Ventral (D/V) axis, and in the third, along the A/P axis. Loss of MOM-5/Frizzled function causes third phase A/P polarity reversals of the migrating DTCs. We show that an over-expression of UNC-5 causes similar DTC A/P polarity reversals and that unc-5 deficits markedly suppress the A/P polarity reversals caused by mutations in mom-5/frizzled . This implicates MOM-5/Frizzled as a negative regulator of unc-5 . We provide further evidence that small GTPases mediate MOM-5’s regulation of unc-5 such that one outcome of impaired function of small GTPases like CED-10/Rac and MIG-2/RhoG is an increase in unc-5 function. The work presented here demonstrates the existence of cross talk between components of the Netrin and Wnt signaling pathways and provides further insights into the way guidance signaling mechanisms are integrated to orchestrate directed cell migration.

Description: by Hugo G. Hilton, Paul J. Norman, Neda Nemat-Gorgani, Ana Goyos, Jill A. Hollenbach, Brenna M. Henn, Christopher R. Gignoux, Lisbeth A. Guethlein, Peter Parham Modulating natural killer cell functions in human immunity and reproduction are diverse interactions between the killer cell immunoglobulin-like receptors (KIR) of Natural Killer (NK) cells and HLA class I ligands on the surface of tissue cells. Dominant interactions are between KIR2DL1 and the C2 epitope of HLA-C and between KIR2DL2/3 and the C1 epitope of HLA-C. KhoeSan hunter-gatherers of Southern Africa represent the earliest population divergence known and are the most genetically diverse indigenous people, qualities reflected in their KIR and HLA genes. Of the ten KhoeSan KIR2DL1 alleles, KIR2DL1*022 and KIR2DL1*026 likely originated in the KhoeSan, and later were transmitted at low frequency to the neighboring Zulus through gene flow. These alleles arose by point mutation from other KhoeSan KIR2DL1 alleles that are more widespread globally. Mutation of KIR2DL1*001 gave rise to KIR2DL1*022 , causing loss of C2 recognition and gain of C1 recognition. This makes KIR2DL1*022 a more avid and specific C1 receptor than any KIR2DL2/3 allotype. Mutation of KIR2DL1*012 gave rise to KIR2DL1*026 , causing premature termination of translation at the end of the transmembrane domain. This makes KIR2DL1*026 a membrane-associated receptor that lacks both a cytoplasmic tail and signaling function. At higher frequencies than their parental allotypes, the combined effect of the KhoeSan-specific KIR2DL1*022 and KIR2DL1*026 is to reduce the frequency of strong inhibitory C2 receptors and increase the frequency of strong inhibitory C1 receptors. Because interaction of KIR2DL1 with C2 is associated with risk of pregnancy disorder, these functional changes are potentially advantageous. Whereas all other KhoeSan KIR2DL1 alleles are present on a wide diversity of centromeric KIR haplotypes, KIR2DL1*026 is present on a single KIR haplotype and KIR2DL1*022 is present on two very similar haplotypes. The high linkage disequilibrium across their haplotypes is consistent with a recent emergence for these KIR2DL1 alleles that have distinctive functions.

Description: by Michael S. Behnke, Asis Khan, Elvin J. Lauron, John R. Jimah, Qiuling Wang, Niraj H. Tolia, L. David Sibley Toxoplasma gondii has evolved a number of strategies to evade immune responses in its many hosts. Previous genetic mapping of crosses between clonal type 1, 2, and 3 strains of T . gondii , which are prevalent in Europe and North America, identified two rhoptry proteins, ROP5 and ROP18, that function together to block innate immune mechanisms activated by interferon gamma (IFNg) in murine hosts. However, the contribution of these and other virulence factors in more genetically divergent South American strains is unknown. Here we utilized a cross between the intermediately virulent North American type 2 ME49 strain and the highly virulent South American type 10 VAND strain to map the genetic basis for differences in virulence in the mouse. Quantitative trait locus (QTL) analysis of this new cross identified one peak that spanned the ROP5 locus on chromosome XII. CRISPR-Cas9 mediated deletion of all copies of ROP5 in the VAND strain rendered it avirulent and complementation confirmed that ROP5 is the major virulence factor accounting for differences between type 2 and type 10 strains. To extend these observations to other virulent South American strains representing distinct genetic populations, we knocked out ROP5 in type 8 TgCtBr5 and type 4 TgCtBr18 strains, resulting in complete loss of virulence in both backgrounds. Consistent with this, polymorphisms that show strong signatures of positive selection in ROP5 were shown to correspond to regions known to interface with host immunity factors. Because ROP5 and ROP18 function together to resist innate immune mechanisms, and a significant interaction between them was identified in a two-locus scan, we also assessed the role of ROP18 in the virulence of South American strains. Deletion of ROP18 in South American type 4, 8, and 10 strains resulted in complete attenuation in contrast to a partial loss of virulence seen for ROP18 knockouts in previously described type 1 parasites. These data show that ROP5 and ROP18 are conserved virulence factors in genetically diverse strains from North and South America, suggesting they evolved to resist innate immune defenses in ancestral T . gondii strains, and they have subsequently diversified under positive selection.

Description: by Ariel Afek, Hila Cohen, Shiran Barber-Zucker, Raluca Gordân, David B. Lukatsky Recent genome-wide experiments in different eukaryotic genomes provide an unprecedented view of transcription factor (TF) binding locations and of nucleosome occupancy. These experiments revealed that a large fraction of TF binding events occur in regions where only a small number of specific TF binding sites (TFBSs) have been detected. Furthermore, in vitro protein-DNA binding measurements performed for hundreds of TFs indicate that TFs are bound with wide range of affinities to different DNA sequences that lack known consensus motifs. These observations have thus challenged the classical picture of specific protein-DNA binding and strongly suggest the existence of additional recognition mechanisms that affect protein-DNA binding preferences. We have previously demonstrated that repetitive DNA sequence elements characterized by certain symmetries statistically affect protein-DNA binding preferences. We call this binding mechanism nonconsensus protein-DNA binding in order to emphasize the point that specific consensus TFBSs do not contribute to this effect. In this paper, using the simple statistical mechanics model developed previously, we calculate the nonconsensus protein-DNA binding free energy for the entire C . elegans and D . melanogaster genomes. Using the available chromatin immunoprecipitation followed by sequencing (ChIP-seq) results on TF-DNA binding preferences for ~100 TFs, we show that DNA sequences characterized by low predicted free energy of nonconsensus binding have statistically higher experimental TF occupancy and lower nucleosome occupancy than sequences characterized by high free energy of nonconsensus binding. This is in agreement with our previous analysis performed for the yeast genome. We suggest therefore that nonconsensus protein-DNA binding assists the formation of nucleosome-free regions, as TFs outcompete nucleosomes at genomic locations with enhanced nonconsensus binding. In addition, here we perform a new, large-scale analysis using in vitro TF-DNA preferences obtained from the universal protein binding microarrays (PBM) for ~90 eukaryotic TFs belonging to 22 different DNA-binding domain types. As a result of this new analysis, we conclude that nonconsensus protein-DNA binding is a widespread phenomenon that significantly affects protein-DNA binding preferences and need not require the presence of consensus (specific) TFBSs in order to achieve genome-wide TF-DNA binding specificity.

Description: by Alexander Suh, Linnéa Smeds, Hans Ellegren The diversification of neoavian birds is one of the most rapid adaptive radiations of extant organisms. Recent whole-genome sequence analyses have much improved the resolution of the neoavian radiation and suggest concurrence with the Cretaceous-Paleogene (K-Pg) boundary, yet the causes of the remaining genome-level irresolvabilities appear unclear. Here we show that genome-level analyses of 2,118 retrotransposon presence/absence markers converge at a largely consistent Neoaves phylogeny and detect a highly differential temporal prevalence of incomplete lineage sorting (ILS), i.e., the persistence of ancestral genetic variation as polymorphisms during speciation events. We found that ILS-derived incongruences are spread over the genome and involve 35% and 34% of the analyzed loci on the autosomes and the Z chromosome, respectively. Surprisingly, Neoaves diversification comprises three adaptive radiations, an initial near-K-Pg super-radiation with highly discordant phylogenetic signals from near-simultaneous speciation events, followed by two post-K-Pg radiations of core landbirds and core waterbirds with much less pronounced ILS. We provide evidence that, given the extreme level of up to 100% ILS per branch in super-radiations, particularly rapid speciation events may neither resemble a fully bifurcating tree nor are they resolvable as such. As a consequence, their complex demographic history is more accurately represented as local networks within a species tree.

Description: by Yu Hasegawa, Deanne Taylor, Dmitry A. Ovchinnikov, Ernst J. Wolvetang, Laurence de Torrenté, Jessica C. Mar An analysis of gene expression variability can provide an insightful window into how regulatory control is distributed across the transcriptome. In a single cell analysis, the inter-cellular variability of gene expression measures the consistency of transcript copy numbers observed between cells in the same population. Application of these ideas to the study of early human embryonic development may reveal important insights into the transcriptional programs controlling this process, based on which components are most tightly regulated. Using a published single cell RNA-seq data set of human embryos collected at four-cell, eight-cell, morula and blastocyst stages, we identified genes with the most stable, invariant expression across all four developmental stages. Stably-expressed genes were found to be enriched for those sharing indispensable features, including essentiality, haploinsufficiency, and ubiquitous expression. The stable genes were less likely to be associated with loss-of-function variant genes or human recessive disease genes affected by a DNA copy number variant deletion, suggesting that stable genes have a functional impact on the regulation of some of the basic cellular processes. Genes with low expression variability at early stages of development are involved in regulation of DNA methylation, responses to hypoxia and telomerase activity, whereas by the blastocyst stage, low-variability genes are enriched for metabolic processes as well as telomerase signaling. Based on changes in expression variability, we identified a putative set of gene expression markers of morulae and blastocyst stages. Experimental validation of a blastocyst-expressed variability marker demonstrated that HDDC2 plays a role in the maintenance of pluripotency in human ES and iPS cells. Collectively our analyses identified new regulators involved in human embryonic development that would have otherwise been missed using methods that focus on assessment of the average expression levels; in doing so, we highlight the value of studying expression variability for single cell RNA-seq data.

Description: by The PLOS Computational Biology Staff

Description: by James Tamerius, Cécile Viboud, Jeffrey Shaman, Gerardo Chowell While a relationship between environmental forcing and influenza transmission has been established in inter-pandemic seasons, the drivers of pandemic influenza remain debated. In particular, school effects may predominate in pandemic seasons marked by an atypical concentration of cases among children. For the 2009 A/H1N1 pandemic, Mexico is a particularly interesting case study due to its broad geographic extent encompassing temperate and tropical regions, well-documented regional variation in the occurrence of pandemic outbreaks, and coincidence of several school breaks during the pandemic period. Here we fit a series of transmission models to daily laboratory-confirmed influenza data in 32 Mexican states using MCMC approaches, considering a meta-population framework or the absence of spatial coupling between states. We use these models to explore the effect of environmental, school–related and travel factors on the generation of spatially-heterogeneous pandemic waves. We find that the spatial structure of the pandemic is best understood by the interplay between regional differences in specific humidity (explaining the occurrence of pandemic activity towards the end of the school term in late May-June 2009 in more humid southeastern states), school vacations (preventing influenza transmission during July-August in all states), and regional differences in residual susceptibility (resulting in large outbreaks in early fall 2009 in central and northern Mexico that had yet to experience fully-developed outbreaks). Our results are in line with the concept that very high levels of specific humidity, as present during summer in southeastern Mexico, favor influenza transmission, and that school cycles are a strong determinant of pandemic wave timing.

Description: by Alireza Alemi, Carlo Baldassi, Nicolas Brunel, Riccardo Zecchina Understanding the theoretical foundations of how memories are encoded and retrieved in neural populations is a central challenge in neuroscience. A popular theoretical scenario for modeling memory function is the attractor neural network scenario, whose prototype is the Hopfield model. The model simplicity and the locality of the synaptic update rules come at the cost of a poor storage capacity, compared with the capacity achieved with perceptron learning algorithms. Here, by transforming the perceptron learning rule, we present an online learning rule for a recurrent neural network that achieves near-maximal storage capacity without an explicit supervisory error signal, relying only upon locally accessible information. The fully-connected network consists of excitatory binary neurons with plastic recurrent connections and non-plastic inhibitory feedback stabilizing the network dynamics; the memory patterns to be memorized are presented online as strong afferent currents, producing a bimodal distribution for the neuron synaptic inputs. Synapses corresponding to active inputs are modified as a function of the value of the local fields with respect to three thresholds. Above the highest threshold, and below the lowest threshold, no plasticity occurs. In between these two thresholds, potentiation/depression occurs when the local field is above/below an intermediate threshold. We simulated and analyzed a network of binary neurons implementing this rule and measured its storage capacity for different sizes of the basins of attraction. The storage capacity obtained through numerical simulations is shown to be close to the value predicted by analytical calculations. We also measured the dependence of capacity on the strength of external inputs. Finally, we quantified the statistics of the resulting synaptic connectivity matrix, and found that both the fraction of zero weight synapses and the degree of symmetry of the weight matrix increase with the number of stored patterns.

Description: by Heather A. Flores, Jaclyn E. Bubnell, Charles F. Aquadro, Daniel A. Barbash Many reproductive proteins from diverse taxa evolve rapidly and adaptively. These proteins are typically involved in late stages of reproduction such as sperm development and fertilization, and are more often functional in males than females. Surprisingly, many germline stem cell (GSC) regulatory genes, which are essential for the earliest stages of reproduction, also evolve adaptively in Drosophila. One example is the bag of marbles ( bam ) gene, which is required for GSC differentiation and germline cyst development in females and for regulating mitotic divisions and entry to spermatocyte differentiation in males. Here we show that the extensive divergence of bam between Drosophila melanogaster and D . simulans affects bam function in females but has no apparent effect in males. We further find that infection with Wolbachia pipientis , an endosymbiotic bacterium that can affect host reproduction through various mechanisms, partially suppresses female sterility caused by bam mutations in D . melanogaster and interacts differentially with bam orthologs from D . melanogaster and D . simulans . We propose that the adaptive evolution of bam has been driven at least in part by the long-term interactions between Drosophila species and Wolbachia . More generally, we suggest that microbial infections of the germline may explain the unexpected pattern of evolution of several GSC regulatory genes.

Description: by Gila Arad, Revital Levy, Iris Nasie, Dalia Hillman, Ziv Rotfogel, Uri Barash, Emmanuelle Supper, Tomer Shpilka, Adi Minis, Raymond Kaempfer

Description: by Diana M. Mitchell, Craig B. Stevens, Ruth A. Frey, Samuel S. Hunter, Ryuichi Ashino, Shoji Kawamura, Deborah L. Stenkamp The signaling molecule retinoic acid (RA) regulates rod and cone photoreceptor fate, differentiation, and survival. Here we elucidate the role of RA in differential regulation of the tandemly-duplicated long wavelength-sensitive ( LWS ) cone opsin genes. Zebrafish embryos were treated with RA from 48 hours post-fertilization (hpf) to 75 hpf, and RNA was isolated from eyes for microarray analysis. ~170 genes showed significantly altered expression, including several transcription factors and components of cellular signaling pathways. Of interest, the LWS1 opsin gene was strongly upregulated by RA. LWS1 is the upstream member of the tandemly duplicated LWS opsin array and is normally not expressed embryonically. Embryos treated with RA 48 hpf to 100 hpf or beyond showed significant reductions in LWS2 -expressing cones in favor of LWS1 -expressing cones. The LWS reporter line, LWS-PAC(H) provided evidence that individual LWS cones switched from LWS2 to LWS1 expression in response to RA. The RA signaling reporter line, RARE : YFP indicated that increased RA signaling in cones was associated with this opsin switch, and experimental reduction of RA signaling in larvae at the normal time of onset of LWS1 expression significantly inhibited LWS1 expression. A role for endogenous RA signaling in regulating differential expression of the LWS genes in postmitotic cones was further supported by the presence of an RA signaling domain in ventral retina of juvenile zebrafish that coincided with a ventral zone of LWS1 expression. This is the first evidence that an extracellular signal may regulate differential expression of opsin genes in a tandemly duplicated array.

Description: by Su Deng, Ingo Bothe, Mary K. Baylies The formation of multinucleated muscle cells through cell-cell fusion is a conserved process from fruit flies to humans. Numerous studies have shown the importance of Arp2/3, its regulators, and branched actin for the formation of an actin structure, the F-actin focus, at the fusion site. This F-actin focus forms the core of an invasive podosome-like structure that is required for myoblast fusion. In this study, we find that the formin Diaphanous (Dia), which nucleates and facilitates the elongation of actin filaments, is essential for Drosophila myoblast fusion. Following cell recognition and adhesion, Dia is enriched at the myoblast fusion site, concomitant with, and having the same dynamics as, the F-actin focus. Through analysis of Dia loss-of-function conditions using mutant alleles but particularly a dominant negative Dia transgene, we demonstrate that reduction in Dia activity in myoblasts leads to a fusion block. Significantly, no actin focus is detected, and neither branched actin regulators, SCAR or WASp, accumulate at the fusion site when Dia levels are reduced. Expression of constitutively active Dia also causes a fusion block that is associated with an increase in highly dynamic filopodia, altered actin turnover rates and F-actin distribution, and mislocalization of SCAR and WASp at the fusion site. Together our data indicate that Dia plays two roles during invasive podosome formation at the fusion site: it dictates the level of linear F-actin polymerization, and it is required for appropriate branched actin polymerization via localization of SCAR and WASp. These studies provide new insight to the mechanisms of cell-cell fusion, the relationship between different regulators of actin polymerization, and invasive podosome formation that occurs in normal development and in disease.

Description: by Jun Ho Ji, Young Lyun Oh, Mineui Hong, Jae Won Yun, Hyun-Woo Lee, DeokGeun Kim, Yongick Ji, Duk-Hwan Kim, Woong-Yang Park, Hyun-Tae Shin, Kyoung-Mee Kim, Myung-Ju Ahn, Keunchil Park, Jong-Mu Sun The genetic landscape of medullary thyroid cancer (MTC) is not yet fully understood, although some oncogenic mutations have been identified. To explore genetic profiles of MTCs, formalin-fixed, paraffin-embedded tumor tissues from MTC patients were assayed on the Ion AmpliSeq Cancer Panel v2. Eighty-four sporadic MTC samples and 36 paired normal thyroid tissues were successfully sequenced. We discovered 101 hotspot mutations in 18 genes in the 84 MTC tissue samples. The most common mutation was in the ret proto-oncogene, which occurred in 47 cases followed by mutations in genes encoding Harvey rat sarcoma viral oncogene homolog (N = 14), serine/threonine kinase 11 (N = 11), v-kit Hardy-Zuckerman 4 feline sarcoma viral oncogene homolog (N = 6), mutL homolog 1 (N = 4), Kiesten rat sarcoma viral oncogene homolog (N = 3) and MET proto-oncogene (N = 3). We also evaluated anaplastic lymphoma kinase ( ALK ) rearrangement by immunohistochemistry and break-apart fluorescence in situ hybridization (FISH). Two of 98 screened cases were positive for ALK FISH. To identify the genomic breakpoint and 5’ fusion partner of ALK , customized targeted cancer panel sequencing was performed using DNA from tumor samples of the two patients. Glutamine:fructose-6-phosphate transaminase 1 ( GFPT1)-ALK and echinoderm microtubule-associated protein-like 4 ( EML4 ) -ALK fusions were identified. Additional PCR analysis, followed by Sanger sequencing, confirmed the GFPT1 - ALK fusion, indicating that the fusion is a result of intra-chromosomal translocation or deletion. Notably, a metastatic MTC case harboring the EML4-ALK fusion showed a dramatic response to an ALK inhibitor, crizotinib. In conclusion, we found several genetic mutations in MTC and are the first to identify ALK fusions in MTC. Our results suggest that the EML4-ALK fusion in MTC may be a potential driver mutation and a valid target of ALK inhibitors. Furthermore, the GFPT1-ALK fusion may be a potential candidate for molecular target therapy.

Description: by Sander Land, Steven A. Niederer Biophysical models of cardiac tension development provide a succinct representation of our understanding of force generation in the heart. The link between protein kinetics and interactions that gives rise to high cooperativity is not yet fully explained from experiments or previous biophysical models. We propose a biophysical ODE-based representation of cross-bridge (XB), tropomyosin and troponin within a contractile regulatory unit (RU) to investigate the mechanisms behind cooperative activation, as well as the role of cooperativity in dynamic tension generation across different species. The model includes cooperative interactions between regulatory units (RU-RU), between crossbridges (XB-XB), as well more complex interactions between crossbridges and regulatory units (XB-RU interactions). For the steady-state force-calcium relationship, our framework predicts that: (1) XB-RU effects are key in shifting the half-activation value of the force-calcium relationship towards lower [Ca 2+ ], but have only small effects on cooperativity. (2) XB-XB effects approximately double the duty ratio of myosin, but do not significantly affect cooperativity. (3) RU-RU effects derived from the long-range action of tropomyosin are a major factor in cooperative activation, with each additional unblocked RU increasing the rate of additional RU’s unblocking. (4) Myosin affinity for short (1–4 RU) unblocked stretches of actin of is very low, and the resulting suppression of force at low [Ca 2+ ] is a major contributor in the biphasic force-calcium relationship. We also reproduce isometric tension development across mouse, rat and human at physiological temperature and pacing rate, and conclude that species differences require only changes in myosin affinity and troponin I/troponin C affinity. Furthermore, we show that the calcium dependence of the rate of tension redevelopment k tr is explained by transient blocking of RU’s by a temporary decrease in XB-RU effects.

Description: by Jonas Paulsen, Odin Gramstad, Philippe Collas The three-dimensional (3D) structure of the genome is important for orchestration of gene expression and cell differentiation. While mapping genomes in 3D has for a long time been elusive, recent adaptations of high-throughput sequencing to chromosome conformation capture (3C) techniques, allows for genome-wide structural characterization for the first time. However, reconstruction of "consensus" 3D genomes from 3C-based data is a challenging problem, since the data are aggregated over millions of cells. Recent single-cell adaptations to the 3C-technique, however, allow for non-aggregated structural assessment of genome structure, but data suffer from sparse and noisy interaction sampling. We present a manifold based optimization (MBO) approach for the reconstruction of 3D genome structure from chromosomal contact data. We show that MBO is able to reconstruct 3D structures based on the chromosomal contacts, imposing fewer structural violations than comparable methods. Additionally, MBO is suitable for efficient high-throughput reconstruction of large systems, such as entire genomes, allowing for comparative studies of genomic structure across cell-lines and different species.

Description: by Hiroo Kenzaki, Shoji Takada Nucleosomes, basic units of chromatin, are known to show spontaneous DNA unwrapping dynamics that are crucial for transcriptional activation, but its structural details are yet to be elucidated. Here, employing a coarse-grained molecular model that captures residue-level structural details up to histone tails, we simulated equilibrium fluctuations and forced unwrapping of single nucleosomes at various conditions. The equilibrium simulations showed spontaneous unwrapping from outer DNA and subsequent rewrapping dynamics, which are in good agreement with experiments. We found several distinct partially unwrapped states of nucleosomes, as well as reversible transitions among these states. At a low salt concentration, histone tails tend to sit in the concave cleft between the histone octamer and DNA, tightening the nucleosome. At a higher salt concentration, the tails tend to bound to the outer side of DNA or be expanded outwards, which led to higher degree of unwrapping. Of the four types of histone tails, H3 and H2B tail dynamics are markedly correlated with partial unwrapping of DNA, and, moreover, their contributions were distinct. Acetylation in histone tails was simply mimicked by changing their charges, which enhanced the unwrapping, especially markedly for H3 and H2B tails.

Description: by Lyne Khair, Richard E. Baker, Erin K. Linehan, Carol E. Schrader, Janet Stavnezer Activation-induced cytidine deaminase (AID) is required for initiation of Ig class switch recombination (CSR) and somatic hypermutation (SHM) of antibody genes during immune responses. AID has also been shown to induce chromosomal translocations, mutations, and DNA double-strand breaks (DSBs) involving non-Ig genes in activated B cells. To determine what makes a DNA site a target for AID-induced DSBs, we identify off-target DSBs induced by AID by performing chromatin immunoprecipitation (ChIP) for Nbs1, a protein that binds DSBs, followed by deep sequencing (ChIP-Seq). We detect and characterize hundreds of off-target AID-dependent DSBs. Two types of tandem repeats are highly enriched within the Nbs1-binding sites: long CA repeats, which can form Z-DNA, and tandem pentamers containing the AID target hotspot WGCW. These tandem repeats are not nearly as enriched at AID-independent DSBs, which we also identified. Msh2, a component of the mismatch repair pathway and important for genome stability, increases off-target DSBs, similar to its effect on Ig switch region DSBs, which are required intermediates during CSR. Most of the off-target DSBs are two-ended, consistent with generation during G1 phase, similar to DSBs in Ig switch regions. However, a minority are one-ended, presumably due to conversion of single-strand breaks to DSBs during replication. One-ended DSBs are repaired by processes involving homologous recombination, including break-induced replication repair, which can lead to genome instability. Off-target DSBs, especially those present during S phase, can lead to chromosomal translocations, deletions and gene amplifications, resulting in the high frequency of B cell lymphomas derived from cells that express or have expressed AID.

Description: by Nicola De Maio, Chieh-Hsi Wu, Kathleen M O’Reilly, Daniel Wilson Phylogeographic methods aim to infer migration trends and the history of sampled lineages from genetic data. Applications of phylogeography are broad, and in the context of pathogens include the reconstruction of transmission histories and the origin and emergence of outbreaks. Phylogeographic inference based on bottom-up population genetics models is computationally expensive, and as a result faster alternatives based on the evolution of discrete traits have become popular. In this paper, we show that inference of migration rates and root locations based on discrete trait models is extremely unreliable and sensitive to biased sampling. To address this problem, we introduce BASTA (BAyesian STructured coalescent Approximation), a new approach implemented in BEAST2 that combines the accuracy of methods based on the structured coalescent with the computational efficiency required to handle more than just few populations. We illustrate the potentially severe implications of poor model choice for phylogeographic analyses by investigating the zoonotic transmission of Ebola virus. Whereas the structured coalescent analysis correctly infers that successive human Ebola outbreaks have been seeded by a large unsampled non-human reservoir population, the discrete trait analysis implausibly concludes that undetected human-to-human transmission has allowed the virus to persist over the past four decades. As genomics takes on an increasingly prominent role informing the control and prevention of infectious diseases, it will be vital that phylogeographic inference provides robust insights into transmission history.

Description: by Nuno M. Oliveira, Esteban Martinez-Garcia, Joao Xavier, William M. Durham, Roberto Kolter, Wook Kim, Kevin R. Foster

Description: by Sebastian Bitzer, Jelle Bruineberg, Stefan J. Kiebel Even for simple perceptual decisions, the mechanisms that the brain employs are still under debate. Although current consensus states that the brain accumulates evidence extracted from noisy sensory information, open questions remain about how this simple model relates to other perceptual phenomena such as flexibility in decisions, decision-dependent modulation of sensory gain, or confidence about a decision. We propose a novel approach of how perceptual decisions are made by combining two influential formalisms into a new model. Specifically, we embed an attractor model of decision making into a probabilistic framework that models decision making as Bayesian inference. We show that the new model can explain decision making behaviour by fitting it to experimental data. In addition, the new model combines for the first time three important features: First, the model can update decisions in response to switches in the underlying stimulus. Second, the probabilistic formulation accounts for top-down effects that may explain recent experimental findings of decision-related gain modulation of sensory neurons. Finally, the model computes an explicit measure of confidence which we relate to recent experimental evidence for confidence computations in perceptual decision tasks.

Description: by Richard Robinson

Description: by Isidoro Feliciello, Ivana Akrap, Đurđica Ugarković Non-coding repetitive DNAs have been proposed to perform a gene regulatory role, however for tandemly repeated satellite DNA no such role was defined until now. Here we provide the first evidence for a role of satellite DNA in the modulation of gene expression under specific environmental conditions. The major satellite DNA TCAST1 in the beetle Tribolium castaneum is preferentially located within pericentromeric heterochromatin but is also dispersed as single repeats or short arrays in the vicinity of protein-coding genes within euchromatin. Our results show enhanced suppression of activity of TCAST1-associated genes and slower recovery of their activity after long-term heat stress relative to the same genes without associated TCAST1 satellite DNA elements. The level of gene suppression is not influenced by the distance of TCAST1 elements from the associated genes up to 40 kb from the genes’ transcription start sites, but it does depend on the copy number of TCAST1 repeats within an element, being stronger for the higher number of copies. The enhanced gene suppression correlates with the enrichment of the repressive histone marks H3K9me2/3 at dispersed TCAST1 elements and their flanking regions as well as with increased expression of TCAST1 satellite DNA. The results reveal transient, RNAi based heterochromatin formation at dispersed TCAST1 repeats and their proximal regions as a mechanism responsible for enhanced silencing of TCAST1-associated genes. Differences in the pattern of distribution of TCAST1 elements contribute to gene expression diversity among T . castaneum strains after long-term heat stress and might have an impact on adaptation to different environmental conditions.

Description: by Patrick Fischer, Martina K. La Rosa, Adriana Schulz, Anette Preiss, Anja C. Nagel In multicellular organisms, growth and proliferation is adjusted to nutritional conditions by a complex signaling network. The Insulin receptor/target of rapamycin (InR/TOR) signaling cascade plays a pivotal role in nutrient dependent growth regulation in Drosophila and mammals alike. Here we identify Cyclin G (CycG) as a regulator of growth and metabolism in Drosophila . CycG mutants have a reduced body size and weight and show signs of starvation accompanied by a disturbed fat metabolism. InR/TOR signaling activity is impaired in cycG mutants, combined with a reduced phosphorylation status of the kinase Akt1 and the downstream factors S6-kinase and eukaryotic translation initiation factor 4E binding protein (4E-BP). Moreover, the expression and accumulation of Drosophila insulin like peptides (dILPs) is disturbed in cycG mutant brains. Using a reporter assay, we show that the activity of one of the first effectors of InR signaling, Phosphoinositide 3-kinase (PI3K92E), is unaffected in cycG mutants. However, the metabolic defects and weight loss in cycG mutants were rescued by overexpression of Akt1 specifically in the fat body and by mutants in widerborst ( wdb ), the B'-subunit of the phosphatase PP2A, known to downregulate Akt1 by dephosphorylation. Together, our data suggest that CycG acts at the level of Akt1 to regulate growth and metabolism via PP2A in Drosophila .

Description: by Kfir Baruch Umansky, Yael Gruenbaum-Cohen, Michael Tsoory, Ester Feldmesser, Dalia Goldenberg, Ori Brenner, Yoram Groner Following myonecrosis, muscle satellite cells proliferate, differentiate and fuse, creating new myofibers. The Runx1 transcription factor is not expressed in naïve developing muscle or in adult muscle tissue. However, it is highly expressed in muscles exposed to myopathic damage yet, the role of Runx1 in muscle regeneration is completely unknown. Our study of Runx1 function in the muscle’s response to myonecrosis reveals that this transcription factor is activated and cooperates with the MyoD and AP-1/c-Jun transcription factors to drive the transcription program of muscle regeneration. Mice lacking dystrophin and muscle Runx1 ( mdx - /Runx1 f/f ), exhibit impaired muscle regeneration leading to age-dependent muscle waste, gradual decrease in motor capabilities and a shortened lifespan. Runx1-deficient primary myoblasts are arrested at cell cycle G 1 and consequently differentiate. Such premature differentiation disrupts the myoblasts’ normal proliferation/differentiation balance, reduces the number and size of regenerating myofibers and impairs muscle regeneration. Our combined Runx1-dependent gene expression, ChIP-seq, ATAC-seq and histone H3K4me1/H3K27ac modification analyses revealed a subset of Runx1-regulated genes that are co-occupied by MyoD and c-Jun in mdx - /Runx1 f/f muscle. The data provide unique insights into the transcriptional program driving muscle regeneration and implicate Runx1 as an important participant in the pathology of muscle wasting diseases.

Description: by Xiao-Jun Xie, Fu-Ning Hsu, Xinsheng Gao, Wu Xu, Jian-Quan Ni, Yue Xing, Liying Huang, Hao-Ching Hsiao, Haiyan Zheng, Chenguang Wang, Yani Zheng, Alus M. Xiaoli, Fajun Yang, Sarah E. Bondos, Jun-Yuan Ji The steroid hormone ecdysone and its receptor (EcR) play critical roles in orchestrating developmental transitions in arthropods. However, the mechanism by which EcR integrates nutritional and developmental cues to correctly activate transcription remains poorly understood. Here, we show that EcR-dependent transcription, and thus, developmental timing in Drosophila , is regulated by CDK8 and its regulatory partner Cyclin C (CycC), and the level of CDK8 is affected by nutrient availability. We observed that cdk8 and cycC mutants resemble EcR mutants and EcR-target genes are systematically down-regulated in both mutants. Indeed, the ability of the EcR-Ultraspiracle (USP) heterodimer to bind to polytene chromosomes and the promoters of EcR target genes is also diminished. Mass spectrometry analysis of proteins that co-immunoprecipitate with EcR and USP identified multiple Mediator subunits, including CDK8 and CycC. Consistently, CDK8-CycC interacts with EcR-USP in vivo; in particular, CDK8 and Med14 can directly interact with the AF1 domain of EcR. These results suggest that CDK8-CycC may serve as transcriptional cofactors for EcR-dependent transcription. During the larval–pupal transition, the levels of CDK8 protein positively correlate with EcR and USP levels, but inversely correlate with the activity of sterol regulatory element binding protein (SREBP), the master regulator of intracellular lipid homeostasis. Likewise, starvation of early third instar larvae precociously increases the levels of CDK8, EcR and USP, yet down-regulates SREBP activity. Conversely, refeeding the starved larvae strongly reduces CDK8 levels but increases SREBP activity. Importantly, these changes correlate with the timing for the larval–pupal transition. Taken together, these results suggest that CDK8-CycC links nutrient intake to developmental transitions (EcR activity) and fat metabolism (SREBP activity) during the larval–pupal transition.

Description: by Blair B. Madison, Arjun N. Jeganathan, Rei Mizuno, Monte M. Winslow, Antoni Castells, Miriam Cuatrecasas, Anil K. Rustgi Let-7 miRNAs comprise one of the largest and most highly expressed family of miRNAs among vertebrates, and is critical for promoting differentiation, regulating metabolism, inhibiting cellular proliferation, and repressing carcinogenesis in a variety of tissues. The large size of the Let-7 family of miRNAs has complicated the development of mutant animal models. Here we describe the comprehensive repression of all Let-7 miRNAs in the intestinal epithelium via low-level tissue-specific expression of the Lin28b RNA-binding protein and a conditional knockout of the MirLet7c-2/Mirlet7b locus. This ablation of Let-7 triggers the development of intestinal adenocarcinomas concomitant with reduced survival. Analysis of both mouse and human intestinal cancer specimens reveals that stem cell markers were significantly associated with loss of Let-7 miRNA expression, and that a number of Let-7 targets were elevated, including Hmga1 and Hmga2. Functional studies in 3-D enteroids revealed that Hmga2 is necessary and sufficient to mediate many characteristics of Let-7 depletion, namely accelerating cell cycle progression and enhancing a stem cell phenotype. In addition, inactivation of a single Hmga2 allele in the mouse intestine epithelium significantly represses tumorigenesis driven by Lin28b. In aggregate, we conclude that Let-7 depletion drives a stem cell phenotype and the development of intestinal cancer, primarily via Hmga2.

Description: by Andrew M. Phillips, Rebecca A. Calvo, Daniel B. Kearns Flagella are assembled sequentially from the inside-out with morphogenetic checkpoints that enforce the temporal order of subunit addition. Here we show that flagellar basal bodies fail to proceed to hook assembly at high frequency in the absence of the monotopic protein SwrB of Bacillus subtilis . Genetic suppressor analysis indicates that SwrB activates the flagellar type III secretion export apparatus by the membrane protein FliP. Furthermore, mutants defective in the flagellar C-ring phenocopy the absence of SwrB for reduced hook frequency and C-ring defects may be bypassed either by SwrB overexpression or by a gain-of-function allele in the polymerization domain of FliG. We conclude that SwrB enhances the probability that the flagellar basal body adopts a conformation proficient for secretion to ensure that rod and hook subunits are not secreted in the absence of a suitable platform on which to polymerize.

Description: by Malachi Griffith, Jason R. Walker, Nicholas C. Spies, Benjamin J. Ainscough, Obi L. Griffith Massively parallel RNA sequencing (RNA-seq) has rapidly become the assay of choice for interrogating RNA transcript abundance and diversity. This article provides a detailed introduction to fundamental RNA-seq molecular biology and informatics concepts. We make available open-access RNA-seq tutorials that cover cloud computing, tool installation, relevant file formats, reference genomes, transcriptome annotations, quality-control strategies, expression, differential expression, and alternative splicing analysis methods. These tutorials and additional training resources are accompanied by complete analysis pipelines and test datasets made available without encumbrance at www.rnaseq.wiki.

Description: by Arjun Bharioke, Dmitri B. Chklovskii Neurons must faithfully encode signals that can vary over many orders of magnitude despite having only limited dynamic ranges. For a correlated signal, this dynamic range constraint can be relieved by subtracting away components of the signal that can be predicted from the past, a strategy known as predictive coding, that relies on learning the input statistics. However, the statistics of input natural signals can also vary over very short time scales e.g., following saccades across a visual scene. To maintain a reduced transmission cost to signals with rapidly varying statistics, neuronal circuits implementing predictive coding must also rapidly adapt their properties. Experimentally, in different sensory modalities, sensory neurons have shown such adaptations within 100 ms of an input change. Here, we show first that linear neurons connected in a feedback inhibitory circuit can implement predictive coding. We then show that adding a rectification nonlinearity to such a feedback inhibitory circuit allows it to automatically adapt and approximate the performance of an optimal linear predictive coding network, over a wide range of inputs, while keeping its underlying temporal and synaptic properties unchanged. We demonstrate that the resulting changes to the linearized temporal filters of this nonlinear network match the fast adaptations observed experimentally in different sensory modalities, in different vertebrate species. Therefore, the nonlinear feedback inhibitory network can provide automatic adaptation to fast varying signals, maintaining the dynamic range necessary for accurate neuronal transmission of natural inputs.

Description: by Ata Ghavidel, Kunal Baxi, Vladimir Ignatchenko, Martin Prusinkiewicz, Terra G. Arnason, Thomas Kislinger, Carlos E. Carvalho, Troy A. A. Harkness Proliferating eukaryotic cells undergo a finite number of cell divisions before irreversibly exiting mitosis. Yet pathways that normally limit the number of cell divisions remain poorly characterized. Here we describe a screen of a collection of 3762 single gene mutants in the yeast Saccharomyces cerevisiae , accounting for 2/3 of annotated yeast ORFs, to search for mutants that undergo an atypically high number of cell divisions. Many of the potential longevity genes map to cellular processes not previously implicated in mitotic senescence, suggesting that regulatory mechanisms governing mitotic exit may be broader than currently anticipated. We focused on an ER-Golgi gene cluster isolated in this screen to determine how these ubiquitous organelles integrate into mitotic longevity. We report that a chronic increase in ER protein load signals an expansion in the assembly of autophagosomes in an Ire1-independent manner, accelerates trafficking of high molecular weight protein aggregates from the cytoplasm to the vacuoles, and leads to a profound enhancement of daughter cell production. We demonstrate that this catabolic network is evolutionarily conserved, as it also extends reproductive lifespan in the nematode Caenorhabditis elegans . Our data provide evidence that catabolism of protein aggregates, a natural byproduct of high protein synthesis and turn over in dividing cells, is among the drivers of mitotic longevity in eukaryotes.

Description: by Roland G. Roberts

Description: by Kenneth H. Wolfe Whole-genome duplications (WGDs) are rare evolutionary events with profound consequences. They double an organism’s genetic content, immediately creating a reproductive barrier between it and its ancestors and providing raw material for the divergence of gene functions between paralogs. Almost all eukaryotic genome sequences bear evidence of ancient WGDs, but the causes of these events and the timing of intermediate steps have been difficult to discern. One of the best-characterized WGDs occurred in the lineage leading to the baker’s yeast Saccharomyces cerevisiae . Marcet-Houben and Gabaldón now show that, rather than simply doubling the DNA of a single ancestor, the yeast WGD likely involved mating between two different ancestral species followed by a doubling of the genome to restore fertility.

Description: by Murat Alp, Vipan K. Parihar, Charles L. Limoli, Francis A. Cucinotta In this work, a stochastic computational model of microscopic energy deposition events is used to study for the first time damage to irradiated neuronal cells of the mouse hippocampus. An extensive library of radiation tracks for different particle types is created to score energy deposition in small voxels and volume segments describing a neuron’s morphology that later are sampled for given particle fluence or dose. Methods included the construction of in silico mouse hippocampal granule cells from neuromorpho.org with spine and filopodia segments stochastically distributed along the dendritic branches. The model is tested with high-energy 56 Fe, 12 C, and 1 H particles and electrons. Results indicate that the tree-like structure of the neuronal morphology and the microscopic dose deposition of distinct particles may lead to different outcomes when cellular injury is assessed, leading to differences in structural damage for the same absorbed dose. The significance of the microscopic dose in neuron components is to introduce specific local and global modes of cellular injury that likely contribute to spine, filopodia, and dendrite pruning, impacting cognition and possibly the collapse of the neuron. Results show that the heterogeneity of heavy particle tracks at low doses, compared to the more uniform dose distribution of electrons, juxtaposed with neuron morphology make it necessary to model the spatial dose painting for specific neuronal components. Going forward, this work can directly support the development of biophysical models of the modifications of spine and dendritic morphology observed after low dose charged particle irradiation by providing accurate descriptions of the underlying physical insults to complex neuron structures at the nano-meter scale.

Description: by Brinda Vallat, Carlos Madrid-Aliste, Andras Fiser Predicting the three-dimensional structure of proteins from their amino acid sequences remains a challenging problem in molecular biology. While the current structural coverage of proteins is almost exclusively provided by template-based techniques, the modeling of the rest of the protein sequences increasingly require template-free methods. However, template-free modeling methods are much less reliable and are usually applicable for smaller proteins, leaving much space for improvement. We present here a novel computational method that uses a library of supersecondary structure fragments, known as Smotifs, to model protein structures. The library of Smotifs has saturated over time, providing a theoretical foundation for efficient modeling. The method relies on weak sequence signals from remotely related protein structures to create a library of Smotif fragments specific to the target protein sequence. This Smotif library is exploited in a fragment assembly protocol to sample decoys, which are assessed by a composite scoring function. Since the Smotif fragments are larger in size compared to the ones used in other fragment-based methods, the proposed modeling algorithm, SmotifTF, can employ an exhaustive sampling during decoy assembly. SmotifTF successfully predicts the overall fold of the target proteins in about 50% of the test cases and performs competitively when compared to other state of the art prediction methods, especially when sequence signal to remote homologs is diminishing. Smotif-based modeling is complementary to current prediction methods and provides a promising direction in addressing the structure prediction problem, especially when targeting larger proteins for modeling.

Description: by Elisabetta Menna, Andrea Disanza, Cinzia Cagnoli, Ursula Schenk, Giuliana Gelsomino, Emanuela Frittoli, Maud Hertzog, Nina Offenhauser, Corinna Sawallisch, Hans-Jürgen Kreienkamp, Frank B. Gertler, Pier Paolo Di Fiore, Giorgio Scita, Michela Matteoli

Description: by Janelle Weaver

Description: by Jeffrey E. Markowitz, William A. Liberti, Grigori Guitchounts, Tarciso Velho, Carlos Lois, Timothy J. Gardner Time-locked sequences of neural activity can be found throughout the vertebrate forebrain in various species and behavioral contexts. From “time cells” in the hippocampus of rodents to cortical activity controlling movement, temporal sequence generation is integral to many forms of learned behavior. However, the mechanisms underlying sequence generation are not well known. Here, we describe a spatial and temporal organization of the songbird premotor cortical microcircuit that supports sparse sequences of neural activity. Multi-channel electrophysiology and calcium imaging reveal that neural activity in premotor cortex is correlated with a length scale of 100 µm. Within this length scale, basal-ganglia–projecting excitatory neurons, on average, fire at a specific phase of a local 30 Hz network rhythm. These results show that premotor cortical activity is inhomogeneous in time and space, and that a mesoscopic dynamical pattern underlies the generation of the neural sequences controlling song.

Description: by Saumya Gupta, Aparna Radhakrishnan, Pandu Raharja-Liu, Gen Lin, Lars M. Steinmetz, Julien Gagneur, Himanshu Sinha Even with identification of multiple causal genetic variants for common human diseases, understanding the molecular processes mediating the causal variants’ effect on the disease remains a challenge. This understanding is crucial for the development of therapeutic strategies to prevent and treat disease. While static profiling of gene expression is primarily used to get insights into the biological bases of diseases, it makes differentiating the causative from the correlative effects difficult, as the dynamics of the underlying biological processes are not monitored. Using yeast as a model, we studied genome-wide gene expression dynamics in the presence of a causal variant as the sole genetic determinant, and performed allele-specific functional validation to delineate the causal effects of the genetic variant on the phenotype. Here, we characterized the precise genetic effects of a functional MKT1 allelic variant in sporulation efficiency variation. A mathematical model describing meiotic landmark events and conditional activation of MKT1 expression during sporulation specified an early meiotic role of this variant. By analyzing the early meiotic genome-wide transcriptional response, we demonstrate an MKT1 -dependent role of novel modulators, namely, RTG1/3 , regulators of mitochondrial retrograde signaling, and DAL82 , regulator of nitrogen starvation, in additively effecting sporulation efficiency. In the presence of functional MKT1 allele, better respiration during early sporulation was observed, which was dependent on the mitochondrial retrograde regulator, RTG3 . Furthermore, our approach showed that MKT1 contributes to sporulation independent of Puf3, an RNA-binding protein that steady-state transcription profiling studies have suggested to mediate MKT1 -pleiotropic effects during mitotic growth. These results uncover interesting regulatory links between meiosis and mitochondrial retrograde signaling. In this study, we highlight the advantage of analyzing allele-specific transcriptional dynamics of mediating genes. Applications in higher eukaryotes can be valuable for inferring causal molecular pathways underlying complex dynamic processes, such as development, physiology and disease progression.

Description: by Elpida Fragouli, Katharina Spath, Samer Alfarawati, Fiona Kaper, Andrew Craig, Claude-Edouard Michel, Felix Kokocinski, Jacques Cohen, Santiago Munne, Dagan Wells Mitochondria play a vital role in embryo development. They are the principal site of energy production and have various other critical cellular functions. Despite the importance of this organelle, little is known about the extent of variation in mitochondrial DNA (mtDNA) between individual human embryos prior to implantation. This study investigated the biological and clinical relevance of the quantity of mtDNA in 379 embryos. These were examined via a combination of microarray comparative genomic hybridisation (aCGH), quantitative PCR and next generation sequencing (NGS), providing information on chromosomal status, amount of mtDNA, and presence of mutations in the mitochondrial genome. The quantity of mtDNA was significantly higher in embryos from older women (P=0.003). Additionally, mtDNA levels were elevated in aneuploid embryos, independent of age (P=0.025). Assessment of clinical outcomes after transfer of euploid embryos to the uterus revealed that blastocysts that successfully implanted tended to contain lower mtDNA quantities than those failing to implant (P=0.007). Importantly, an mtDNA quantity threshold was established, above which implantation was never observed. Subsequently, the predictive value of this threshold was confirmed in an independent blinded prospective study, indicating that abnormal mtDNA levels are present in 30% of non-implanting euploid embryos, but are not seen in embryos forming a viable pregnancy. NGS did not reveal any increase in mutation in blastocysts with elevated mtDNA levels. The results of this study suggest that increased mtDNA may be related to elevated metabolism and are associated with reduced viability, a possibility consistent with the ‘quiet embryo’ hypothesis. Importantly, the findings suggest a potential role for mitochondria in female reproductive aging and the genesis of aneuploidy. Of clinical significance, we propose that mtDNA content represents a novel biomarker with potential value for in vitro fertilisation (IVF) treatment, revealing chromosomally normal blastocysts incapable of producing a viable pregnancy.

Description: by Simon van Vliet, Martin Ackermann Multicellular eukaryotes can perform functions that exceed the possibilities of an individual cell. These functions emerge through interactions between differentiated cells that are precisely arranged in space. Bacteria also form multicellular collectives that consist of differentiated but genetically identical cells. How does the functionality of these collectives depend on the spatial arrangement of the differentiated bacteria? In a previous issue of PLOS Biology , van Gestel and colleagues reported an elegant example of how the spatial arrangement of differentiated cells gives rise to collective behavior in Bacillus subtilus colonies, further demonstrating the similarity of bacterial collectives to higher multicellular organisms.

Description: by Ingvild Flåtten, Solveig Fossum-Raunehaug, Riikka Taipale, Silje Martinsen, Kirsten Skarstad The bacterial replication cycle is driven by the DnaA protein which cycles between the active ATP-bound form and the inactive ADP-bound form. It has been suggested that DnaA also is the main controller of initiation frequency. Initiation is thought to occur when enough ATP-DnaA has accumulated. In this work we have performed cell cycle analysis of cells that contain a surplus of ATP-DnaA and asked whether initiation then occurs earlier. It does not. Cells with more than a 50% increase in the concentration of ATP-DnaA showed no changes in timing of replication. We suggest that although ATP-DnaA is the main actor in initiation of replication, its accumulation does not control the time of initiation. ATP-DnaA is the motor that drives the initiation process, but other factors will be required for the exact timing of initiation in response to the cell’s environment. We also investigated the in vivo roles of datA dependent DnaA inactivation (DDAH) and the DnaA-binding protein DiaA. Loss of DDAH affected the cell cycle machinery only during slow growth and made it sensitive to the concentration of DiaA protein. The result indicates that compromised cell cycle machines perform in a less robust manner.

Description: by Po-Wei Chen, Luis L. Fonseca, Yusuf A. Hannun, Eberhard O. Voit The article demonstrates that computational modeling has the capacity to convert metabolic snapshots, taken sequentially over time, into a description of cellular, dynamic strategies. The specific application is a detailed analysis of a set of actions with which Saccharomyces cerevisiae responds to heat stress. Using time dependent metabolic concentration data, we use a combination of mathematical modeling, reverse engineering, and optimization to infer dynamic changes in enzyme activities within the sphingolipid pathway. The details of the sphingolipid responses to heat stress are important, because they guide some of the longer-term alterations in gene expression, with which the cells adapt to the increased temperature. The analysis indicates that all enzyme activities in the system are affected and that the shapes of the time trends in activities depend on the fatty-acyl CoA chain lengths of the different ceramide species in the system.

Description: by Ran Lin, Danielle Rittenhouse, Katelyn Sweeney, Prasanth Potluri, Douglas C. Wallace The heavy consumption of ethanol can lead to alcohol use disorders (AUDs) which impact patients, their families, and societies. Yet the genetic and physiological factors that predispose humans to AUDs remain unclear. One hypothesis is that alterations in mitochondrial function modulate neuronal sensitivity to ethanol exposure. Using Drosophila genetics we report that inactivation of the mitochondrial outer membrane translocator protein 18kDa (TSPO), also known as the peripheral benzodiazepine receptor, affects ethanol sedation and tolerance in male flies. Knockdown of dTSPO in adult male neurons results in increased sensitivity to ethanol sedation, and this effect requires the dTSPO depletion-mediated increase in reactive oxygen species (ROS) production and inhibition of caspase activity in fly heads. Systemic loss of dTSPO in male flies blocks the development of tolerance to repeated ethanol exposures, an effect that is not seen when dTSPO is only inactivated in neurons. Female flies are naturally more sensitive to ethanol than males, and female fly heads have strikingly lower levels of dTSPO mRNA than males. Hence, mitochondrial TSPO function plays an important role in ethanol sensitivity and tolerance. Since a large array of benzodiazepine analogues have been developed that interact with the peripheral benzodiazepine receptor, the mitochondrial TSPO might provide an important new target for treating AUDs.

Description: by Christine E. Cucinotta, Alexandria N. Young, Kristin M. Klucevsek, Karen M. Arndt Eukaryotes regulate gene expression and other nuclear processes through the posttranslational modification of histones. In S . cerevisiae , the mono-ubiquitylation of histone H2B on lysine 123 (H2B K123ub) affects nucleosome stability, broadly influences gene expression and other DNA-templated processes, and is a prerequisite for additional conserved histone modifications that are associated with active transcription, namely the methylation of lysine residues in H3. While the enzymes that promote these chromatin marks are known, regions of the nucleosome required for the recruitment of these enzymes are undefined. To identify histone residues required for H2B K123ub, we exploited a functional interaction between the ubiquitin-protein ligase, Rkr1/Ltn1, and H2B K123ub in S . cerevisiae . Specifically, we performed a synthetic lethal screen with cells lacking RKR1 and a comprehensive library of H2A and H2B residue substitutions, and identified H2A residues that are required for H2B K123ub. Many of these residues map to the nucleosome acidic patch. The substitutions in the acidic patch confer varying histone modification defects downstream of H2B K123ub, indicating that this region contributes differentially to multiple histone modifications. Interestingly, substitutions in the acidic patch result in decreased recruitment of H2B K123ub machinery to active genes and defects in transcription elongation and termination. Together, our findings reveal a role for the nucleosome acidic patch in recruitment of histone modification machinery and maintenance of transcriptional integrity.

Description: by Fong Kuan Wong, Ji-Feng Fei, Felipe Mora-Bermúdez, Elena Taverna, Christiane Haffner, Jun Fu, Konstantinos Anastassiadis, A. Francis Stewart, Wieland B. Huttner The evolutionary expansion of the neocortex in mammals has been linked to enlargement of the subventricular zone (SVZ) and increased proliferative capacity of basal progenitors (BPs), notably basal radial glia (bRG). The transcription factor Pax6 is known to be highly expressed in primate, but not mouse, BPs. Here, we demonstrate that sustaining Pax6 expression selectively in BP-genic apical radial glia (aRG) and their BP progeny of embryonic mouse neocortex suffices to induce primate-like progenitor behaviour. Specifically, we conditionally expressed Pax6 by in utero electroporation using a novel, Tis21 –CreER T2 mouse line. This expression altered aRG cleavage plane orientation to promote bRG generation, increased cell-cycle re-entry of BPs, and ultimately increased upper-layer neuron production. Upper-layer neuron production was also increased in double-transgenic mouse embryos with sustained Pax6 expression in the neurogenic lineage. Strikingly, increased BPs existed not only in the SVZ but also in the intermediate zone of the neocortex of these double-transgenic mouse embryos. In mutant mouse embryos lacking functional Pax6, the proportion of bRG among BPs was reduced. Our data identify specific Pax6 effects in BPs and imply that sustaining this Pax6 function in BPs could be a key aspect of SVZ enlargement and, consequently, the evolutionary expansion of the neocortex.

Description: by Jason M. Glanz, Courtney R. Kraus, Matthew F. Daley The recent United States measles epidemic has sparked another contentious national discussion about childhood vaccination. A growing number of parents are expressing concerns about the safety of vaccines, often fueled by misinformation from the internet, books, and other nonmedical sources. Many of these concerned parents are choosing to refuse or delay childhood vaccines, placing their children and surrounding communities at risk for serious diseases that are nearly 100% preventable with vaccination. Between 10% and 15% of parents are asking physicians to space out the timing of vaccines, which often poses an ethical dilemma for physicians. This trend reflects a tension between personal liberty and public health, as parents fight to control the decisions that affect the health of their children and public health officials strive to maintain high immunization rates to prevent outbreaks of vaccine-preventable diseases. Interventions to address this emerging public health issue are needed. We describe a framework by which web-based interventions can be used to help parents make evidence-based decisions about childhood vaccinations.

Description: by Marina Marcet-Houben, Toni Gabaldón Whole-genome duplications have shaped the genomes of several vertebrate, plant, and fungal lineages. Earlier studies have focused on establishing when these events occurred and on elucidating their functional and evolutionary consequences, but we still lack sufficient understanding of how genome duplications first originated. We used phylogenomics to study the ancient genome duplication occurred in the yeast Saccharomyces cerevisiae lineage and found compelling evidence for the existence of a contemporaneous interspecies hybridization. We propose that the genome doubling was a direct consequence of this hybridization and that it served to provide stability to the recently formed allopolyploid. This scenario provides a mechanism for the origin of this ancient duplication and the lineage that originated from it and brings a new perspective to the interpretation of the origin and consequences of whole-genome duplications.

Description: by Pengyi Yang, Xiaofeng Zheng, Vivek Jayaswal, Guang Hu, Jean Yee Hwa Yang, Raja Jothi Cell signaling underlies transcription/epigenetic control of a vast majority of cell-fate decisions. A key goal in cell signaling studies is to identify the set of kinases that underlie key signaling events. In a typical phosphoproteomics study, phosphorylation sites (substrates) of active kinases are quantified proteome-wide. By analyzing the activities of phosphorylation sites over a time-course, the temporal dynamics of signaling cascades can be elucidated. Since many substrates of a given kinase have similar temporal kinetics, clustering phosphorylation sites into distinctive clusters can facilitate identification of their respective kinases. Here we present a knowledge-based CLUster Evaluation (CLUE) approach for identifying the most informative partitioning of a given temporal phosphoproteomics data. Our approach utilizes prior knowledge, annotated kinase-substrate relationships mined from literature and curated databases, to first generate biologically meaningful partitioning of the phosphorylation sites and then determine key kinases associated with each cluster. We demonstrate the utility of the proposed approach on two time-series phosphoproteomics datasets and identify key kinases associated with human embryonic stem cell differentiation and insulin signaling pathway. The proposed approach will be a valuable resource in the identification and characterizing of signaling networks from phosphoproteomics data.

Description: by Deborah A. Striegel, Manami Hara, Vipul Periwal Pancreatic islets of Langerhans consist of endocrine cells, primarily α, β and δ cells, which secrete glucagon, insulin, and somatostatin, respectively, to regulate plasma glucose. β cells form irregular locally connected clusters within islets that act in concert to secrete insulin upon glucose stimulation. Due to the central functional significance of this local connectivity in the placement of β cells in an islet, it is important to characterize it quantitatively. However, quantification of the seemingly stochastic cytoarchitecture of β cells in an islet requires mathematical methods that can capture topological connectivity in the entire β-cell population in an islet. Graph theory provides such a framework. Using large-scale imaging data for thousands of islets containing hundreds of thousands of cells in human organ donor pancreata, we show that quantitative graph characteristics differ between control and type 2 diabetic islets. Further insight into the processes that shape and maintain this architecture is obtained by formulating a stochastic theory of β-cell rearrangement in whole islets, just as the normal equilibrium distribution of the Ornstein-Uhlenbeck process can be viewed as the result of the interplay between a random walk and a linear restoring force. Requiring that rearrangements maintain the observed quantitative topological graph characteristics strongly constrained possible processes. Our results suggest that β-cell rearrangement is dependent on its connectivity in order to maintain an optimal cluster size in both normal and T2D islets.

Description: by Marcel Ander, Sivaraman Subramaniam, Karim Fahmy, A. Francis Stewart, Erik Schäffer Repair of DNA breaks by single-strand annealing (SSA) is a major mechanism for the maintenance of genomic integrity. SSA is promoted by proteins (single-strand-annealing proteins [SSAPs]), such as eukaryotic RAD52 and λ phage Redβ. These proteins use a short single-stranded region to find sequence identity and initiate homologous recombination. However, it is unclear how SSAPs detect homology and catalyze annealing. Using single-molecule experiments, we provide evidence that homology is recognized by Redβ monomers that weakly hold single DNA strands together. Once annealing begins, dimerization of Redβ clamps the double-stranded region and nucleates nucleoprotein filament growth. In this manner, DNA clamping ensures and secures a successful detection for DNA sequence homology. The clamp is characterized by a structural change of Redβ and a remarkable stability against force up to 200 pN. Our findings not only present a detailed explanation for SSAP action but also identify the DNA clamp as a very stable, noncovalent, DNA–protein interaction.

Description: by Kenta Funayama, Genki Minamisawa, Nobuyoshi Matsumoto, Hiroshi Ban, Allen W. Chan, Norio Matsuki, Timothy H. Murphy, Yuji Ikegaya Animals are constantly exposed to the time-varying visual world. Because visual perception is modulated by immediately prior visual experience, visual cortical neurons may register recent visual history into a specific form of offline activity and link it to later visual input. To examine how preceding visual inputs interact with upcoming information at the single neuron level, we designed a simple stimulation protocol in which a brief, orientated flashing stimulus was subsequently coupled to visual stimuli with identical or different features. Using in vivo whole-cell patch-clamp recording and functional two-photon calcium imaging from the primary visual cortex (V1) of awake mice, we discovered that a flash of sinusoidal grating per se induces an early, transient activation as well as a long-delayed reactivation in V1 neurons. This late response, which started hundreds of milliseconds after the flash and persisted for approximately 2 s, was also observed in human V1 electroencephalogram. When another drifting grating stimulus arrived during the late response, the V1 neurons exhibited a sublinear, but apparently increased response, especially to the same grating orientation. In behavioral tests of mice and humans, the flashing stimulation enhanced the detection power of the identically orientated visual stimulation only when the second stimulation was presented during the time window of the late response. Therefore, V1 late responses likely provide a neural basis for admixing temporally separated stimuli and extracting identical features in time-varying visual environments.

Description: by Ghanim Ullah, Yina Wei, Markus A Dahlem, Martin Wechselberger, Steven J Schiff Cell volume changes are ubiquitous in normal and pathological activity of the brain. Nevertheless, we know little of how cell volume affects neuronal dynamics. We here performed the first detailed study of the effects of cell volume on neuronal dynamics. By incorporating cell swelling together with dynamic ion concentrations and oxygen supply into Hodgkin-Huxley type spiking dynamics, we demonstrate the spontaneous transition between epileptic seizure and spreading depression states as the cell swells and contracts in response to changes in osmotic pressure. Our use of volume as an order parameter further revealed a dynamical definition for the experimentally described physiological ceiling that separates seizure from spreading depression, as well as predicted a second ceiling that demarcates spreading depression from anoxic depolarization. Our model highlights the neuroprotective role of glial K buffering against seizures and spreading depression, and provides novel insights into anoxic depolarization and the relevant cell swelling during ischemia. We argue that the dynamics of seizures, spreading depression, and anoxic depolarization lie along a continuum of the repertoire of the neuron membrane that can be understood only when the dynamic ion concentrations, oxygen homeostasis,and cell swelling in response to osmotic pressure are taken into consideration. Our results demonstrate the feasibility of a unified framework for a wide range of neuronal behaviors that may be of substantial importance in the understanding of and potentially developing universal intervention strategies for these pathological states.

Description: by John R. Houser, Craig Barnhart, Daniel R. Boutz, Sean M. Carroll, Aurko Dasgupta, Joshua K. Michener, Brittany D. Needham, Ophelia Papoulas, Viswanadham Sridhara, Dariya K. Sydykova, Christopher J. Marx, M. Stephen Trent, Jeffrey E. Barrick, Edward M. Marcotte, Claus O. Wilke How do bacteria regulate their cellular physiology in response to starvation? Here, we present a detailed characterization of Escherichia coli growth and starvation over a time-course lasting two weeks. We have measured multiple cellular components, including RNA and proteins at deep genomic coverage, as well as lipid modifications and flux through central metabolism. Our study focuses on the physiological response of E . coli in stationary phase as a result of being starved for glucose, not on the genetic adaptation of E . coli to utilize alternative nutrients. In our analysis, we have taken advantage of the temporal correlations within and among RNA and protein abundances to identify systematic trends in gene regulation. Specifically, we have developed a general computational strategy for classifying expression-profile time courses into distinct categories in an unbiased manner. We have also developed, from dynamic models of gene expression, a framework to characterize protein degradation patterns based on the observed temporal relationships between mRNA and protein abundances. By comparing and contrasting our transcriptomic and proteomic data, we have identified several broad physiological trends in the E . coli starvation response. Strikingly, mRNAs are widely down-regulated in response to glucose starvation, presumably as a strategy for reducing new protein synthesis. By contrast, protein abundances display more varied responses. The abundances of many proteins involved in energy-intensive processes mirror the corresponding mRNA profiles while proteins involved in nutrient metabolism remain abundant even though their corresponding mRNAs are down-regulated.

Description: by Shaun S. Sanders, Dale D. O. Martin, Stefanie L. Butland, Mathieu Lavallée-Adam, Diego Calzolari, Chris Kay, John R. Yates, Michael R. Hayden Palmitoylation involves the reversible posttranslational addition of palmitate to cysteines and promotes membrane binding and subcellular localization. Recent advancements in the detection and identification of palmitoylated proteins have led to multiple palmitoylation proteomics studies but these datasets are contained within large supplemental tables, making downstream analysis and data mining time-consuming and difficult. Consequently, we curated the data from 15 palmitoylation proteomics studies into one compendium containing 1,838 genes encoding palmitoylated proteins; representing approximately 10% of the genome. Enrichment analysis revealed highly significant enrichments for Gene Ontology biological processes, pathway maps, and process networks related to the nervous system. Strikingly, 41% of synaptic genes encode a palmitoylated protein in the compendium. The top disease associations included cancers and diseases and disorders of the nervous system, with Schizophrenia, HD, and pancreatic ductal carcinoma among the top five, suggesting that aberrant palmitoylation may play a pivotal role in the balance of cell death and survival. This compendium provides a much-needed resource for cell biologists and the palmitoylation field, providing new perspectives for cancer and neurodegeneration.

Description: by Liat Rockah-Shmuel, Ágnes Tóth-Petróczy, Dan S. Tawfik Systematic mappings of the effects of protein mutations are becoming increasingly popular. Unexpectedly, these experiments often find that proteins are tolerant to most amino acid substitutions, including substitutions in positions that are highly conserved in nature. To obtain a more realistic distribution of the effects of protein mutations, we applied a laboratory drift comprising 17 rounds of random mutagenesis and selection of M.HaeIII, a DNA methyltransferase. During this drift, multiple mutations gradually accumulated. Deep sequencing of the drifted gene ensembles allowed determination of the relative effects of all possible single nucleotide mutations. Despite being averaged across many different genetic backgrounds, about 67% of all nonsynonymous, missense mutations were evidently deleterious, and an additional 16% were likely to be deleterious. In the early generations, the frequency of most deleterious mutations remained high. However, by the 17th generation, their frequency was consistently reduced, and those remaining were accepted alongside compensatory mutations. The tolerance to mutations measured in this laboratory drift correlated with sequence exchanges seen in M.HaeIII’s natural orthologs. The biophysical constraints dictating purging in nature and in this laboratory drift also seemed to overlap. Our experiment therefore provides an improved method for measuring the effects of protein mutations that more closely replicates the natural evolutionary forces, and thereby a more realistic view of the mutational space of proteins.

Description: by Julien Y. Dutheil, Kasper Munch, Kiwoong Nam, Thomas Mailund, Mikkel H. Schierup The human and chimpanzee X chromosomes are less divergent than expected based on autosomal divergence. We study incomplete lineage sorting patterns between humans, chimpanzees and gorillas to show that this low divergence can be entirely explained by megabase-sized regions comprising one-third of the X chromosome, where polymorphism in the human-chimpanzee ancestral species was severely reduced. We show that background selection can explain at most 10% of this reduction of diversity in the ancestor. Instead, we show that several strong selective sweeps in the ancestral species can explain it. We also report evidence of population specific sweeps in extant humans that overlap the regions of low diversity in the ancestral species. These regions further correspond to chromosomal sections shown to be devoid of Neanderthal introgression into modern humans. This suggests that the same X-linked regions that undergo selective sweeps are among the first to form reproductive barriers between diverging species. We hypothesize that meiotic drive is the underlying mechanism causing these two observations.

Description: by Greicy H. Goto, Sevil Zencir, Yukinori Hirano, Hiroo Ogi, Andreas Ivessa, Katsunori Sugimoto Telomeres, the ends of linear eukaryotic chromosomes, have a specialized chromatin structure that provides a stable chromosomal terminus. In budding yeast Rap1 protein binds to telomeric TG repeat and negatively regulates telomere length. Here we show that binding of multiple Rap1 proteins stimulates DNA double-stranded break (DSB) induction at both telomeric and non-telomeric regions. Consistent with the role of DSB induction, Rap1 stimulates nearby recombination events in a dosage-dependent manner. Rap1 recruits Rif1 and Rif2 to telomeres, but neither Rif1 nor Rif2 is required for DSB induction. Rap1-mediated DSB induction involves replication fork progression but inactivation of checkpoint kinase Mec1 does not affect DSB induction. Rap1 tethering shortens artificially elongated telomeres in parallel with telomerase inhibition, and this telomere shortening does not require homologous recombination. These results suggest that Rap1 contributes to telomere homeostasis by promoting chromosome breakage.

Description: by Laureline Logiaco, René Quilodran, Emmanuel Procyk, Angelo Arleo The frontal cortex controls behavioral adaptation in environments governed by complex rules. Many studies have established the relevance of firing rate modulation after informative events signaling whether and how to update the behavioral policy. However, whether the spatiotemporal features of these neuronal activities contribute to encoding imminent behavioral updates remains unclear. We investigated this issue in the dorsal anterior cingulate cortex (dACC) of monkeys while they adapted their behavior based on their memory of feedback from past choices. We analyzed spike trains of both single units and pairs of simultaneously recorded neurons using an algorithm that emulates different biologically plausible decoding circuits. This method permits the assessment of the performance of both spike-count and spike-timing sensitive decoders. In response to the feedback, single neurons emitted stereotypical spike trains whose temporal structure identified informative events with higher accuracy than mere spike count. The optimal decoding time scale was in the range of 70–200 ms, which is significantly shorter than the memory time scale required by the behavioral task. Importantly, the temporal spiking patterns of single units were predictive of the monkeys’ behavioral response time. Furthermore, some features of these spiking patterns often varied between jointly recorded neurons. All together, our results suggest that dACC drives behavioral adaptation through complex spatiotemporal spike coding. They also indicate that downstream networks, which decode dACC feedback signals, are unlikely to act as mere neural integrators.

Description: by Shuai Yuan, H. Richard Johnston, Guosheng Zhang, Yun Li, Yi-Juan Hu, Zhaohui S. Qin With rapid decline of the sequencing cost, researchers today rush to embrace whole genome sequencing (WGS), or whole exome sequencing (WES) approach as the next powerful tool for relating genetic variants to human diseases and phenotypes. A fundamental step in analyzing WGS and WES data is mapping short sequencing reads back to the reference genome. This is an important issue because incorrectly mapped reads affect the downstream variant discovery, genotype calling and association analysis. Although many read mapping algorithms have been developed, the majority of them uses the universal reference genome and do not take sequence variants into consideration. Given that genetic variants are ubiquitous, it is highly desirable if they can be factored into the read mapping procedure. In this work, we developed a novel strategy that utilizes genotypes obtained a priori to customize the universal haploid reference genome into a personalized diploid reference genome. The new strategy is implemented in a program named RefEditor. When applying RefEditor to real data, we achieved encouraging improvements in read mapping, variant discovery and genotype calling. Compared to standard approaches, RefEditor can significantly increase genotype calling consistency (from 43% to 61% at 4X coverage; from 82% to 92% at 20X coverage) and reduce Mendelian inconsistency across various sequencing depths. Because many WGS and WES studies are conducted on cohorts that have been genotyped using array-based genotyping platforms previously or concurrently, we believe the proposed strategy will be of high value in practice, which can also be applied to the scenario where multiple NGS experiments are conducted on the same cohort. The RefEditor sources are available at https://github.com/superyuan/refeditor.

Description: by James J. Collins, Xiaowen Hou, Elena V. Romanova, Bramwell G. Lambrus, Claire M. Miller, Amir Saberi, Jonathan V. Sweedler, Phillip A. Newmark

Description: by Alexey A. Gritsenko, Marc Hulsman, Marcel J. T. Reinders, Dick de Ridder Translation of RNA to protein is a core process for any living organism. While for some steps of this process the effect on protein production is understood, a holistic understanding of translation still remains elusive. In silico modelling is a promising approach for elucidating the process of protein synthesis. Although a number of computational models of the process have been proposed, their application is limited by the assumptions they make. Ribosome profiling (RP), a relatively new sequencing-based technique capable of recording snapshots of the locations of actively translating ribosomes, is a promising source of information for deriving unbiased data-driven translation models. However, quantitative analysis of RP data is challenging due to high measurement variance and the inability to discriminate between the number of ribosomes measured on a gene and their speed of translation. We propose a solution in the form of a novel multi-scale interpretation of RP data that allows for deriving models with translation dynamics extracted from the snapshots. We demonstrate the usefulness of this approach by simultaneously determining for the first time per-codon translation elongation and per-gene translation initiation rates of Saccharomyces cerevisiae from RP data for two versions of the Totally Asymmetric Exclusion Process (TASEP) model of translation. We do this in an unbiased fashion, by fitting the models using only RP data with a novel optimization scheme based on Monte Carlo simulation to keep the problem tractable. The fitted models match the data significantly better than existing models and their predictions show better agreement with several independent protein abundance datasets than existing models. Results additionally indicate that the tRNA pool adaptation hypothesis is incomplete, with evidence suggesting that tRNA post-transcriptional modifications and codon context may play a role in determining codon elongation rates.

Description: by Sourav Roy, Tusar T. Saha, Lisa Johnson, Bo Zhao, Jisu Ha, Kevin P. White, Thomas Girke, Zhen Zou, Alexander S. Raikhel In multicellular organisms, development, growth and reproduction require coordinated expression of numerous functional and regulatory genes. Insects, in addition to being the most speciose animal group with enormous biological and economical significance, represent outstanding model organisms for studying regulation of synchronized gene expression due to their rapid development and reproduction. Disease-transmitting female mosquitoes have adapted uniquely for ingestion and utilization of the huge blood meal required for swift reproductive events to complete egg development within a 72-h period. We investigated the network of regulatory factors mediating sequential gene expression in the fat body, a multifunctional organ analogous to the vertebrate liver and adipose tissue, of the female Aedes aegypti mosquito. Transcriptomic and bioinformatics analyses revealed that ~7500 transcripts are differentially expressed in four sequential waves during the 72-h reproductive period. A combination of RNA-interference gene-silencing and in-vitro organ culture identified the major regulators for each of these waves. Amino acids (AAs) regulate the first wave of gene activation between 3 h and 12 h post-blood meal (PBM). During the second wave, between 12 h and 36 h, most genes are highly upregulated by a synergistic action of AAs, 20-hydroxyecdysone (20E) and the Ecdysone-Receptor (EcR). Between 36 h and 48 h, the third wave of gene activation—regulated mainly by HR3—occurs. Juvenile Hormone (JH) and its receptor Methoprene-Tolerant (Met) are major regulators for the final wave between 48 h and 72 h. Each of these key regulators also has repressive effects on one or more gene sets. Our study provides a better understanding of the complexity of the regulatory mechanisms related to temporal coordination of gene expression during reproduction. We have detected the novel function of 20E/EcR responsible for transcriptional repression. This study also reveals the previously unidentified large-scale effects of HR3 and JH/Met on transcriptional regulation during the termination of vitellogenesis and remodeling of the fat body.

Description: by Jintao Li, Yu Zhao, Huangwei Chu, Likai Wang, Yanru Fu, Ping Liu, Narayana Upadhyaya, Chunli Chen, Tongmin Mou, Yuqi Feng, Prakash Kumar, Jian Xu Little is known about how the size of meristem cells is regulated and whether it participates in the control of meristem size in plants. Here, we report our findings on shoebox ( shb ), a mild gibberellin (GA) deficient rice mutant that has a short root meristem size. Quantitative analysis of cortical cell length and number indicates that shb has shorter, rather than fewer, cells in the root meristem until around the fifth day after sowing, from which the number of cortical cells is also reduced. These defects can be either corrected by exogenous application of bioactive GA or induced in wild-type roots by a dose-dependent inhibitory effect of paclobutrazol on GA biosynthesis, suggesting that GA deficiency is the primary cause of shb mutant phenotypes. SHB encodes an AP2/ERF transcription factor that directly activates transcription of the GA biosynthesis gene KS1 . Thus, root meristem size in rice is modulated by SHB-mediated GA biosynthesis that regulates the elongation and proliferation of meristem cells in a developmental stage-specific manner.

Description: by Prashant R. Desai, Lasse van Wijlick, Dagmar Kurtz, Mateusz Juchimiuk, Joachim F. Ernst Candida albicans is a common commensal in the human gut but in predisposed patients it can become an important human fungal pathogen. As a commensal, C . albicans adapts to low-oxygen conditions and represses its hyphal development by the transcription factor Efg1, which under normoxia activates filamentation. The repressive hypoxic but not the normoxic function of Efg1 required its unmodified N-terminus, was prevented by phosphomimetic residues at normoxic phosphorylation sites T179 and T206 and occurred only at temperatures ≤35°C. Genome-wide binding sites for native Efg1 identified 300 hypoxia-specific target genes, which overlapped partially with hypoxic binding sites for Ace2, a known positive regulator of hypoxic filamentation. Transcriptional analyses revealed that EFG1 , ACE2 and their identified targets BCR1 and BRG1 encode an interconnected regulatory hub, in which Efg1/Bcr1 act as negative and Ace2/Brg1 act as positive regulators of gene expression under hypoxia. In this circuit, the hypoxic function of Ace2 was stimulated by elevated CO 2 levels. The hyperfilamentous phenotype of efg1 and bcr1 mutants depended on Ace2/Brg1 regulators and required increased expression of genes encoding Cek1 MAP kinase and its downstream target Cph1. The intricate temperature-dependent regulatory mechanisms under hypoxia suggest that C . albicans restricts hyphal morphogenesis in oxygen-poor body niches, possibly to persist as a commensal in the human host.

Description: by Emmanuel Haillot, Maria Dolores Molina, François Lapraz, Thierry Lepage Specification of the dorsal-ventral axis in the highly regulative sea urchin embryo critically relies on the zygotic expression of nodal , but whether maternal factors provide the initial spatial cue to orient this axis is not known. Although redox gradients have been proposed to entrain the dorsal-ventral axis by acting upstream of nodal , manipulating the activity of redox gradients only has modest consequences, suggesting that other factors are responsible for orienting nodal expression and defining the dorsal-ventral axis. Here we uncover the function of Panda, a maternally provided transforming growth factor beta (TGF-β) ligand that requires the activin receptor-like kinases (Alk) Alk3/6 and Alk1/2 receptors to break the radial symmetry of the embryo and orient the dorsal-ventral axis by restricting nodal expression. We found that the double inhibition of the bone morphogenetic protein (BMP) type I receptors Alk3/6 and Alk1/2 causes a phenotype dramatically more severe than the BMP2/4 loss-of-function phenotype, leading to extreme ventralization of the embryo through massive ectopic expression of nodal , suggesting that an unidentified signal acting through BMP type I receptors cooperates with BMP2/4 to restrict nodal expression. We identified this ligand as the product of maternal Panda mRNA. Double inactivation of panda and bmp2/4 led to extreme ventralization, mimicking the phenotype caused by inactivation of the two BMP receptors. Inhibition of maternal panda mRNA translation disrupted the early spatial restriction of nodal , leading to persistent massive ectopic expression of nodal on the dorsal side despite the presence of Lefty. Phylogenetic analysis indicates that Panda is not a prototypical BMP ligand but a member of a subfamily of TGF-β distantly related to Inhibins, Lefty, and TGF-β that includes Maverick from Drosophila and GDF15 from vertebrates. Indeed, overexpression of Panda does not appear to directly or strongly activate phosphoSmad1/5/8 signaling, suggesting that although this TGF-β may require Alk1/2 and/or Alk3/6 to antagonize nodal expression, it may do so by sequestering a factor essential for Nodal signaling, by activating a non-Smad pathway downstream of the type I receptors, or by activating extremely low levels of pSmad1/5/8. We provide evidence that, although panda mRNA is broadly distributed in the early embryo, local expression of panda mRNA efficiently orients the dorsal-ventral axis and that Panda activity is required locally in the early embryo to specify this axis. Taken together, these findings demonstrate that maternal panda mRNA is both necessary and sufficient to orient the dorsal-ventral axis. These results therefore provide evidence that in the highly regulative sea urchin embryo, the activity of spatially restricted maternal factors regulates patterning along the dorsal-ventral axis.

Description: by Michelle K. Scotland, Justin M. H. Heltzel, James E. Kath, Jung-Suk Choi, Anthony J. Berdis, Joseph J. Loparo, Mark D. Sutton Translesion DNA synthesis (TLS) by specialized DNA polymerases (Pols) is a conserved mechanism for tolerating replication blocking DNA lesions. The actions of TLS Pols are managed in part by ring-shaped sliding clamp proteins. In addition to catalyzing TLS, altered expression of TLS Pols impedes cellular growth. The goal of this study was to define the relationship between the physiological function of Escherichia coli Pol IV in TLS and its ability to impede growth when overproduced. To this end, 13 novel Pol IV mutants were identified that failed to impede growth. Subsequent analysis of these mutants suggest that overproduced levels of Pol IV inhibit E . coli growth by gaining inappropriate access to the replication fork via a Pol III-Pol IV switch that is mechanistically similar to that used under physiological conditions to coordinate Pol IV-catalyzed TLS with Pol III-catalyzed replication. Detailed analysis of one mutant, Pol IV-T120P, and two previously described Pol IV mutants impaired for interaction with either the rim (Pol IV R ) or the cleft (Pol IV C ) of the β sliding clamp revealed novel insights into the mechanism of the Pol III-Pol IV switch. Specifically, Pol IV-T120P retained complete catalytic activity in vitro but, like Pol IV R and Pol IV C , failed to support Pol IV TLS function in vivo . Notably, the T120P mutation abrogated a biochemical interaction of Pol IV with Pol III that was required for Pol III-Pol IV switching. Taken together, these results support a model in which Pol III-Pol IV switching involves interaction of Pol IV with Pol III, as well as the β clamp rim and cleft. Moreover, they provide strong support for the view that Pol III-Pol IV switching represents a vitally important mechanism for regulating TLS in vivo by managing access of Pol IV to the DNA.

Description: by Johannes Raffler, Nele Friedrich, Matthias Arnold, Tim Kacprowski, Rico Rueedi, Elisabeth Altmaier, Sven Bergmann, Kathrin Budde, Christian Gieger, Georg Homuth, Maik Pietzner, Werner Römisch-Margl, Konstantin Strauch, Henry Völzke, Melanie Waldenberger, Henri Wallaschofski, Matthias Nauck, Uwe Völker, Gabi Kastenmüller, Karsten Suhre Genome-wide association studies with metabolic traits (mGWAS) uncovered many genetic variants that influence human metabolism. These genetically influenced metabotypes (GIMs) contribute to our metabolic individuality, our capacity to respond to environmental challenges, and our susceptibility to specific diseases. While metabolic homeostasis in blood is a well investigated topic in large mGWAS with over 150 known loci, metabolic detoxification through urinary excretion has only been addressed by few small mGWAS with only 11 associated loci so far. Here we report the largest mGWAS to date, combining targeted and non-targeted 1 H NMR analysis of urine samples from 3,861 participants of the SHIP-0 cohort and 1,691 subjects of the KORA F4 cohort. We identified and replicated 22 loci with significant associations with urinary traits, 15 of which are new ( HIBCH , CPS1 , AGXT , XYLB , TKT , ETNPPL , SLC6A19 , DMGDH , SLC36A2 , GLDC , SLC6A13 , ACSM3 , SLC5A11 , PNMT , SLC13A3 ). Two-thirds of the urinary loci also have a metabolite association in blood. For all but one of the 6 loci where significant associations target the same metabolite in blood and urine, the genetic effects have the same direction in both fluids. In contrast, for the SLC5A11 locus, we found increased levels of myo -inositol in urine whereas mGWAS in blood reported decreased levels for the same genetic variant. This might indicate less effective re-absorption of myo -inositol in the kidneys of carriers. In summary, our study more than doubles the number of known loci that influence urinary phenotypes. It thus allows novel insights into the relationship between blood homeostasis and its regulation through excretion. The newly discovered loci also include variants previously linked to chronic kidney disease ( CPS1 , SLC6A13 ), pulmonary hypertension ( CPS1 ), and ischemic stroke ( XYLB ). By establishing connections from gene to disease via metabolic traits our results provide novel hypotheses about molecular mechanisms involved in the etiology of diseases.

Description: by Santiago Schnell

Description: by Yuzuru Imai, Yoshito Kobayashi, Tsuyoshi Inoshita, Hongrui Meng, Taku Arano, Kengo Uemura, Takeshi Asano, Kenji Yoshimi, Chang-Liang Zhang, Gen Matsumoto, Toshiyuki Ohtsuka, Ryoichiro Kageyama, Hiroshi Kiyonari, Go Shioi, Nobuyuki Nukina, Nobutaka Hattori, Ryosuke Takahashi Leucine-rich repeat kinase 2 (LRRK2) is a key molecule in the pathogenesis of familial and idiopathic Parkinson’s disease (PD). We have identified two novel LRRK2-associated proteins, a HECT-type ubiquitin ligase, HERC2, and an adaptor-like protein with six repeated Neuralized domains, NEURL4. LRRK2 binds to NEURL4 and HERC2 via the LRRK2 Ras of complex proteins (ROC) domain and NEURL4, respectively. HERC2 and NEURL4 link LRRK2 to the cellular vesicle transport pathway and Notch signaling, through which the LRRK2 complex promotes the recycling of the Notch ligand Delta-like 1 (Dll1)/Delta (Dl) through the modulation of endosomal trafficking. This process negatively regulates Notch signaling through cis -inhibition by stabilizing Dll1/Dl, which accelerates neural stem cell differentiation and modulates the function and survival of differentiated dopaminergic neurons. These effects are strengthened by the R1441G ROC domain-mutant of LRRK2. These findings suggest that the alteration of Notch signaling in mature neurons is a component of PD etiology linked to LRRK2 .

Description: by Richard Robinson

Description: by Alexander Ullrich, Mathias A. Böhme, Johannes Schöneberg, Harald Depner, Stephan J. Sigrist, Frank Noé Synaptic vesicle fusion is mediated by SNARE proteins forming in between synaptic vesicle (v-SNARE) and plasma membrane (t-SNARE), one of which is Syntaxin-1A. Although exocytosis mainly occurs at active zones, Syntaxin-1A appears to cover the entire neuronal membrane. By using STED super-resolution light microscopy and image analysis of Drosophila neuro-muscular junctions, we show that Syntaxin-1A clusters are more abundant and have an increased size at active zones. A computational particle-based model of syntaxin cluster formation and dynamics is developed. The model is parametrized to reproduce Syntaxin cluster-size distributions found by STED analysis, and successfully reproduces existing FRAP results. The model shows that the neuronal membrane is adjusted in a way to strike a balance between having most syntaxins stored in large clusters, while still keeping a mobile fraction of syntaxins free or in small clusters that can efficiently search the membrane or be traded between clusters. This balance is subtle and can be shifted toward almost no clustering and almost complete clustering by modifying the syntaxin interaction energy on the order of only 1 k B T. This capability appears to be exploited at active zones. The larger active-zone syntaxin clusters are more stable and provide regions of high docking and fusion capability, whereas the smaller clusters outside may serve as flexible reserve pool or sites of spontaneous ectopic release.

Description: by Malika Ihle, Bart Kempenaers, Wolfgang Forstmeier Research on mate choice has primarily focused on preferences for quality indicators, assuming that all individuals show consensus about who is the most attractive. However, in some species, mating preferences seem largely individual-specific, suggesting that they might target genetic or behavioral compatibility. Few studies have quantified the fitness consequences of allowing versus preventing such idiosyncratic mate choice. Here, we report on an experiment that controls for variation in overall partner quality and show that zebra finch ( Taeniopygia guttata ) pairs that resulted from free mate choice achieved a 37% higher reproductive success than pairs that were forced to mate. Cross-fostering of freshly laid eggs showed that embryo mortality (before hatching) primarily depended on the identity of the genetic parents, whereas offspring mortality during the rearing period depended on foster-parent identity. Therefore, preventing mate choice should lead to an increase in embryo mortality if mate choice targets genetic compatibility (for embryo viability), and to an increase in offspring mortality if mate choice targets behavioral compatibility (for better rearing). We found that pairs from both treatments showed equal rates of embryo mortality, but chosen pairs were better at raising offspring. These results thus support the behavioral, but not the genetic, compatibility hypothesis. Further exploratory analyses reveal several differences in behavior and fitness components between “free-choice” and “forced” pairs.

Description: by Stephan Köhler, Friederike Schmid, Giovanni Settanni Fibrinogen is a serum multi-chain protein which, when activated, aggregates to form fibrin, one of the main components of a blood clot. Fibrinolysis controls blood clot dissolution through the action of the enzyme plasmin, which cleaves fibrin at specific locations. Although the main biochemical factors involved in fibrin formation and lysis have been identified, a clear mechanistic picture of how these processes take place is not available yet. This picture would be instrumental, for example, for the design of improved thrombolytic or anti-haemorrhagic strategies, as well as, materials with improved biocompatibility. Here, we present extensive molecular dynamics simulations of fibrinogen which reveal large bending motions centered at a hinge point in the coiled-coil regions of the molecule. This feature, likely conserved across vertebrates according to our analysis, suggests an explanation for the mechanism of exposure to lysis of the plasmin cleavage sites on fibrinogen coiled-coil region. It also explains the conformational variability of fibrinogen observed during its adsorption on inorganic surfaces and it is supposed to play a major role in the determination of the hydrodynamic properties of fibrinogen. In addition the simulations suggest how the dynamics of the D region of fibrinogen may contribute to the allosteric regulation of the blood coagulation cascade through a dynamic coupling between the a- and b-holes, important for fibrin polymerization, and the integrin binding site P1.

Description: by Sergei Maslov, Kim Sneppen Populations of species in ecosystems are often constrained by availability of resources within their environment. In effect this means that a growth of one population, needs to be balanced by comparable reduction in populations of others. In neutral models of biodiversity all populations are assumed to change incrementally due to stochastic births and deaths of individuals. Here we propose and model another redistribution mechanism driven by abrupt and severe reduction in size of the population of a single species freeing up resources for the remaining ones. This mechanism may be relevant e.g. for communities of bacteria, with strain-specific collapses caused e.g. by invading bacteriophages, or for other ecosystems where infectious diseases play an important role. The emergent dynamics of our system is characterized by cyclic ‘‘diversity waves’’ triggered by collapses of globally dominating populations. The population diversity peaks at the beginning of each wave and exponentially decreases afterwards. Species abundances have bimodal time-aggregated distribution with the lower peak formed by populations of recently collapsed or newly introduced species while the upper peak - species that has not yet collapsed in the current wave. In most waves both upper and lower peaks are composed of several smaller peaks. This self-organized hierarchical peak structure has a long-term memory transmitted across several waves. It gives rise to a scale-free tail of the time-aggregated population distribution with a universal exponent of 1.7. We show that diversity wave dynamics is robust with respect to variations in the rules of our model such as diffusion between multiple environments, species-specific growth and extinction rates, and bet-hedging strategies.

Description: by Eva Mejia-Ramirez, Oliver Limbo, Petra Langerak, Paul Russell Abstract Phosphorylation of histone H2AX by ATM and ATR establishes a chromatin recruitment platform for DNA damage response proteins. Phospho-H2AX (γH2AX) has been most intensively studied in the context of DNA double-strand breaks caused by exogenous clastogens, but recent studies suggest that DNA replication stress also triggers formation of γH2A (ortholog of γH2AX) in Schizosaccharomyces pombe . Here, a focused genetic screen in fission yeast reveals that γH2A is critical when there are defects in Replication Factor C (RFC), which loads proliferating cell nuclear antigen (PCNA) clamp onto duplex DNA. Surprisingly Chk1, Cds1/Chk2 and the Rad9-Hus1-Rad1 checkpoint clamp, which are crucial for surviving many genotoxins, are fully dispensable in RFC-defective cells. Immunoblot analysis confirms that Rad9-Hus1-Rad1 is not required for formation of γH2A by Rad3/ATR in S-phase. Defects in DNA polymerase epsilon, which binds PCNA in the replisome, also create an acute need for γH2A. These requirements for γH2A were traced to its role in docking with Brc1, which is a 6-BRCT-domain protein that is structurally related to budding yeast Rtt107 and mammalian PTIP. Brc1, which localizes at stalled replication forks by binding γH2A, prevents aberrant formation of Replication Protein A (RPA) foci in RFC-impaired cells, suggesting that Brc1-coated chromatin stabilizes replisomes when PCNA or DNA polymerase availability limits DNA synthesis.

Description: by Song Chen, Niels Lory, Johannes Stauber, Ute Hoecker The Arabidopsis COP1/SPA E3 ubiquitin ligase is a key negative regulator that represses light signaling in darkness by targeting transcription factors involved in the light response for degradation. The COP1/SPA complex consists of COP1 and members of the four-member SPA protein family (SPA1-SPA4). Genetic analysis indicated that COP1/SPA2 function is particularly strongly repressed by light when compared to complexes carrying the other three SPAs, thereby promoting a light response after exposure of plants to extremely low light. Here, we show that the SPA2 protein is degraded within 5–15 min after exposure of dark-grown seedlings to a pulse of light. Phytochrome photoreceptors are required for the rapid degradation of SPA2 in red, far-red and also in blue light, whereas cryptochromes are not involved in the rapid, blue light-induced reduction in SPA2 protein levels. These results uncover a photoreceptor-specific mechanism of light-induced inhibition of COP1/SPA2 function. Phytochrome A (phyA) is required for the severe blue light responsiveness of spa triple mutants expressing only SPA2, thus confirming the important role of phyA in downregulating SPA2 function in blue light. In blue light, SPA2 forms a complex with cryptochrome 1 (cry1), but not with cryptochrome 2 (cry2) in vivo , indicating that the lack of a rapid blue light response of the SPA2 protein is only in part caused by a failure to interact with cryptochromes. Since SPA1 interacts with both cry1 and cry2, these results provide first molecular evidence that the light-regulation of different SPA proteins diverged during evolution. SPA2 degradation in the light requires COP1 and the COP1-interacting coiled-coil domain of SPA2, supporting that SPA2 is ubiquitinated by COP1. We propose that light perceived by phytochromes causes a switch in the ubiquitination activity of COP1/SPA2 from ubiquitinating downstream substrates to ubiquitinating SPA2, which subsequently causes a repression of COP1/SPA2 function.

Description: by The PLOS Genetics Editors

Description: by Hannah Edwards, Charlotte M. Deane Several protein structure classification schemes exist that partition the protein universe into structural units called folds. Yet these schemes do not discuss how these units sit relative to each other in a global structure space. In this paper we construct networks that describe such global relationships between folds in the form of structural bridges. We generate these networks using four different structural alignment methods across multiple score thresholds. The networks constructed using the different methods remain a similar distance apart regardless of the probability threshold defining a structural bridge. This suggests that at least some structural bridges are method specific and that any attempt to build a picture of structural space should not be reliant on a single structural superposition method. Despite these differences all representations agree on an organisation of fold space into five principal community structures: all- α , all- β sandwiches, all- β barrels, α / β and α + β . We project estimated fold ages onto the networks and find that not only are the pairings of unconnected folds associated with higher age differences than bridged folds, but this difference increases with the number of networks displaying an edge. We also examine different centrality measures for folds within the networks and how these relate to fold age. While these measures interpret the central core of fold space in varied ways they all identify the disposition of ancestral folds to fall within this core and that of the more recently evolved structures to provide the peripheral landscape. These findings suggest that evolutionary information is encoded along these structural bridges. Finally, we identify four highly central pivotal folds representing dominant topological features which act as key attractors within our landscapes.

Description: by Ayal Lavi, Omri Perez, Uri Ashery Neuronal microcircuits generate oscillatory activity, which has been linked to basic functions such as sleep, learning and sensorimotor gating. Although synaptic release processes are well known for their ability to shape the interaction between neurons in microcircuits, most computational models do not simulate the synaptic transmission process directly and hence cannot explain how changes in synaptic parameters alter neuronal network activity. In this paper, we present a novel neuronal network model that incorporates presynaptic release mechanisms, such as vesicle pool dynamics and calcium-dependent release probability, to model the spontaneous activity of neuronal networks. The model, which is based on modified leaky integrate-and-fire neurons, generates spontaneous network activity patterns, which are similar to experimental data and robust under changes in the model's primary gain parameters such as excitatory postsynaptic potential and connectivity ratio. Furthermore, it reliably recreates experimental findings and provides mechanistic explanations for data obtained from microelectrode array recordings, such as network burst termination and the effects of pharmacological and genetic manipulations. The model demonstrates how elevated asynchronous release, but not spontaneous release, synchronizes neuronal network activity and reveals that asynchronous release enhances utilization of the recycling vesicle pool to induce the network effect. The model further predicts a positive correlation between vesicle priming at the single-neuron level and burst frequency at the network level; this prediction is supported by experimental findings. Thus, the model is utilized to reveal how synaptic release processes at the neuronal level govern activity patterns and synchronization at the network level.

Description: by Trevor L. Cameron, Katrina M. Bell, Irma L. Gresshoff, Lisa Sampurno, Lorna Mullan, Joerg Ermann, Laurie H. Glimcher, Raymond P. Boot-Handford, John F. Bateman Schmid metaphyseal chondrodysplasia (MCDS) involves dwarfism and growth plate cartilage hypertrophic zone expansion resulting from dominant mutations in the hypertrophic zone collagen, Col10a1 . Mouse models phenocopying MCDS through the expression of an exogenous misfolding protein in the endoplasmic reticulum (ER) in hypertrophic chondrocytes have demonstrated the central importance of ER stress in the pathology of MCDS. The resultant unfolded protein response (UPR) in affected chondrocytes involved activation of canonical ER stress sensors, IRE1, ATF6, and PERK with the downstream effect of disrupted chondrocyte differentiation. Here, we investigated the role of the highly conserved IRE1/XBP1 pathway in the pathology of MCDS. Mice with a MCDS collagen X p.N617K knock-in mutation ( ColX N617K ) were crossed with mice in which Xbp1 was inactivated specifically in cartilage ( Xbp1 CartΔEx2 ), generating the compound mutant, C/X . The severity of dwarfism and hypertrophic zone expansion in C/X did not differ significantly from ColX N617K , revealing surprising redundancy for the IRE1/XBP1 UPR pathway in the pathology of MCDS. Transcriptomic analyses of hypertrophic zone cartilage identified differentially expressed gene cohorts in MCDS that are pathologically relevant (XBP1-independent) or pathologically redundant (XBP1-dependent). XBP1-independent gene expression changes included large-scale transcriptional attenuation of genes encoding secreted proteins and disrupted differentiation from proliferative to hypertrophic chondrocytes. Moreover, these changes were consistent with disruption of C/EBP-β, a master regulator of chondrocyte differentiation, by CHOP, a transcription factor downstream of PERK that inhibits C/EBP proteins, and down-regulation of C/EBP-β transcriptional co-factors, GADD45-β and RUNX2. Thus we propose that the pathology of MCDS is underpinned by XBP1 independent UPR-induced dysregulation of C/EBP-β-mediated chondrocyte differentiation. Our data suggest that modulation of C/EBP-β activity in MCDS chondrocytes may offer therapeutic opportunities.

Description: by Richard Robinson