Algorithms in nature: the convergence of systems biology and computational thinking

Biologists rely on computational methods to analyze and integrate large data sets, while several computational methods were inspired by the high-level design principles of biological systems. This Perspectives discusses the recent convergence of these two ways of thinking

Small RNA populations revealed by blocking rRNA fragments in Drosophila melanogaster reproductive tissues

RNA interference (RNAi) is a complex and highly conserved regulatory mechanism mediated via small RNAs (sRNAs). Recent technical advances in high throughput sequencing have enabled an increasingly detailed analysis of sRNA abundances and profiles in specific body parts and tissues. This enables investigations of the localized roles of microRNAs (miRNAs) and small interfering RNAs (siRNAs). However, variation in the proportions of non-coding RNAs in the samples being compared can hinder these analyses. Specific tissues may vary significantly in the proportions of fragments of longer non-coding RNAs (such as ribosomal RNA or transfer RNA) present, potentially reflecting tissue-specific differences in biological functions. For example, in Drosophila, some tissues contain a highly abundant 30nt rRNA fragment (the 2S rRNA) as well as abundant 5’ and 3’ terminal rRNA fragments. These can pose difficulties for the construction of sRNA libraries as they can swamp the sequencing space and obscure sRNA abundances. Here we addressed this problem and present a modified “rRNA blocking” protocol for the construction of high-definition (HD) adapter sRNA libraries, in D. melanogaster reproductive tissues. The results showed that 2S rRNAs targeted by blocking oligos were reduced from >80% to < 0.01% total reads. In addition, the use of multiple rRNA blocking oligos to bind the most abundant rRNA fragments allowed us to reveal the underlying sRNA populations at increased resolution. Side-by-side comparisons of sequencing libraries of blocked and non-blocked samples revealed that rRNA blocking did not change the miRNA populations present, but instead enhanced their abundances. We suggest that this rRNA blocking procedure offers the potential to improve the in-depth analysis of differentially expressed sRNAs within and across different tissues

The genomes of two key bumblebee species with primitive eusocial organization

Background: The shift from solitary to social behavior is one of the major evolutionary transitions. Primitively eusocial bumblebees are uniquely placed to illuminate the evolution of highly eusocial insect societies. Bumblebees are also invaluable natural and agricultural pollinators, and there is widespread concern over recent population declines in some species. High-quality genomic data will inform key aspects of bumblebee biology, including susceptibility to implicated population viability threats. Results: We report the high quality draft genome sequences of Bombus terrestris and Bombus impatiens, two ecologically dominant bumblebees and widely utilized study species. Comparing these new genomes to those of the highly eusocial honeybee Apis mellifera and other Hymenoptera, we identify deeply conserved similarities, as well as novelties key to the biology of these organisms. Some honeybee genome features thought to underpin advanced eusociality are also present in bumblebees, indicating an earlier evolution in the bee lineage. Xenobiotic detoxification and immune genes are similarly depauperate in bumblebees and honeybees, and multiple categories of genes linked to social organization, including development and behavior, show high conservation. Key differences identified include a bias in bumblebee chemoreception towards gustation from olfaction, and striking differences in microRNAs, potentially responsible for gene regulation underlying social and other traits. Conclusions: These two bumblebee genomes provide a foundation for post-genomic research on these key pollinators and insect societies. Overall, gene repertoires suggest that the route to advanced eusociality in bees was mediated by many small changes in many genes and processes, and not by notable expansion or depauperation

Differential and converging molecular mechanisms of antidepressants' action in the hippocampal dentate gyrus

Major depression is a highly prevalent, multidimensional disorder. Although several classes of antidepressants (ADs) are currently available, treatment efficacy is limited, and relapse rates are high; thus, there is a need to find better therapeutic strategies. Neuroplastic changes in brain regions such as the hippocampal dentate gyrus (DG) accompany depression and its amelioration with ADs. In this study, the unpredictable chronic mild stress (uCMS) rat model of depression was used to determine the molecular mediators of chronic stress and the targets of four ADs with different pharmacological profiles (fluoxetine, imipramine, tianeptine, and agomelatine) in the hippocampal DG. All ADs, except agomelatine, reversed the depression-like behavior and neuroplastic changes produced by uCMS. Chronic stress induced significant molecular changes that were generally reversed by fluoxetine, imipramine, and tianeptine. Fluoxetine primarily acted on neurons to reduce the expression of pro-inflammatory response genes and increased a set of genes involved in cell metabolism. Similarities were found between the molecular actions and targets of imipramine and tianeptine that activated pathways related to cellular protection. Agomelatine presented a unique profile, with pronounced effects on genes related to Rho-GTPase-related pathways in oligodendrocytes and neurons. These differential molecular signatures of ADs studied contribute to our understanding of the processes implicated in the onset and treatment of depression-like symptoms.Patricia Patricio, Antonio Mateus-Pinheiro, Monica Morais, and Nuno Dinis Alves received fellowships from the Portuguese Foundation for Science and Technology (FCT). Michal Korostynski and Marcin Piechota were funded by the POIG De-Me-Ter 3.1 and NCN 2011/03/D/NZ3/01686 grants. This study was co-funded by the Life and Health Sciences Research Institute (ICVS) and ON. 2-O NOVO NORTE-North Portugal Regional Operational Programme 2007/2013, of the National Strategic Reference Framework (NSRF) 2007/ 2013, through the European Regional Development Fund (ERDF) and by the SwitchBox Consortium (Contract FP7-Health-F2-2010-259772 from the European Union). The authors declare no conflict of interest

Tbx6 Regulates Left/Right Patterning in Mouse Embryos through Effects on Nodal Cilia and Perinodal Signaling

Background: The determination of left/right body axis during early embryogenesis sets up a developmental cascade that coordinates the development of the viscera and is essential to the correct placement and alignment of organ systems and vasculature. Defective left-right patterning can lead to congenital cardiac malformations, vascular anomalies and other serious health problems. Here we describe a novel role for the T-box transcription factor gene Tbx6 in left/right body axis determination in the mouse. Results: Embryos lacking Tbx6 show randomized embryo turning and heart looping. Our results point to multiple mechanisms for this effect. First, Dll1, a direct target of Tbx6, is down regulated around the node in Tbx6 mutants and there is a subsequent decrease in nodal signaling, which is required for laterality determination. Secondly, in spite of a lack of expression of Tbx6 in the node, we document a profound effect of the Tbx6 mutation on the morphology and motility of nodal cilia. This results in the loss of asymmetric calcium signaling at the periphery of the node, suggesting that unidirectional nodal flow is disrupted. To carry out these studies, we devised a novel method for direct labeling and live imaging cilia in vivo using a genetically-encoded fluorescent protein fusion that labels tubulin, combined with laser point scanning confocal microscopy for direct visualization of cilia movement. Conclusions: We conclude that the transcription factor gene Tbx6 is essential for correct left/right axis determination in th

Narrowing down the conditions for extinction of Pavlovian feature-positive discriminations in humans

The aim of this study was to delineate the minimal conditions for extinction of Pavlovian modulation in humans. Previous experiments at our lab showed that, after X-- A+/A- acquisition training, X- trials did not extinguish differential X-- A+/A- responding, while X-- A- trials did. Additionally, X-- A- extinction training seemed only to extinguish differential X-- A+/A- responding, while leaving differential responding on a concurrently trained Y [Symbol: see text] B+/B- discrimination intact. It thus seemed that the X-- A+/A- discrimination can only be extinguished by X-- A- extinction trials. (Rescorla, Journal of Experimental Psychology: Animal Behavior Processes 12, 16-24, 1986), on the other hand, found that the minimal conditions for extinction were broader in pigeons: Namely, he found that an acquired X-- A+/A- discrimination could be extinguished by presenting the original feature X in combination with a different target (B) that was minimally trained as an exciter. We thus wanted to examine whether this was also the case in humans. We found that nonreinforced X-- B- presentations did not abolish discriminative X-- A/A responding when target B was a nonreinforced stimulus. Nonreinforced X-- B- trials did extinguish the X-- A+/A- discrimination when target B had previously been trained as a target for modulation (X-- B+/B- or Y [Symbol: see text] B+/B- training) or as a reinforced exciter (B+). Our results thusf parallel and extend those in nonhuman animals (Rescorla, Journal of Experimental Psychology: Animal Behavior Processes 12, 16-24, 1986)

Criteria for the definition of Pituitary Tumor Centers of Excellence (PTCOE): A Pituitary Society Statement

Introduction With the goal of generate uniform criteria among centers dealing with pituitary tumors and to enhance patient care, the Pituitary Society decided to generate criteria for developing Pituitary Tumors Centers of Excellence (PTCOE). Methods To develop that task, a group of ten experts served as a Task Force and through two years of iterative work an initial draft was elaborated. This draft was discussed, modified and finally approved by the Board of Directors of the Pituitary Society. Such document was presented and debated at a specific session of the Congress of the Pituitary Society, Orlando 2017, and suggestions were incorporated. Finally the document was distributed to a large group of global experts that introduced further modifications with final endorsement. Results After five years of iterative work a document with the ideal criteria for a PTCOE is presented. Conclusions Acknowledging that very few centers in the world, if any, likely fulfill the requirements here presented, the document may be a tool to guide improvements of care delivery to patients with pituitary disorders. All these criteria must be accommodated to the regulations and organization of Health of a given country

From Dynamic Expression Patterns to Boundary Formation in the Presomitic Mesoderm

The segmentation of the vertebrate body is laid down during early embryogenesis. The formation of signaling gradients, the periodic expression of genes of the Notch-, Fgf- and Wnt-pathways and their interplay in the unsegmented presomitic mesoderm (PSM) precedes the rhythmic budding of nascent somites at its anterior end, which later develops into epithelialized structures, the somites. Although many in silico models describing partial aspects of somitogenesis already exist, simulations of a complete causal chain from gene expression in the growth zone via the interaction of multiple cells to segmentation are rare. Here, we present an enhanced gene regulatory network (GRN) for mice in a simulation program that models the growing PSM by many virtual cells and integrates WNT3A and FGF8 gradient formation, periodic gene expression and Delta/Notch signaling. Assuming Hes7 as core of the somitogenesis clock and LFNG as modulator, we postulate a negative feedback of HES7 on Dll1 leading to an oscillating Dll1 expression as seen in vivo. Furthermore, we are able to simulate the experimentally observed wave of activated NOTCH (NICD) as a result of the interactions in the GRN. We esteem our model as robust for a wide range of parameter values with the Hes7 mRNA and protein decays exerting a strong influence on the core oscillator. Moreover, our model predicts interference between Hes1 and HES7 oscillators when their intrinsic frequencies differ. In conclusion, we have built a comprehensive model of somitogenesis with HES7 as core oscillator that is able to reproduce many experimentally observed data in mice