Search

[Background]: Systems biologists study interaction data to understand the behaviour of whole cell systems, and their environment, at a molecular level. In order to effectively achieve this goal, it is critical that researchers have high quality interaction datasets available to them, in a standard data format, and also a suite of tools with which to analyse such data and form experimentally testable hypotheses from them. The PSI-MI XML standard interchange format was initially published in 2004, and expanded in 2007 to enable the download and interchange of molecular interaction data. PSI-XML2.5 was designed to describe experimental data and to date has fulfilled this basic requirement. However, new use cases have arisen that the format cannot properly accommodate. These include data abstracted from more than one publication such as allosteric/cooperative interactions and protein complexes, dynamic interactions and the need to link kinetic and affinity data to specific mutational changes. [Results]: The Molecular Interaction workgroup of the HUPO-PSI has extended the existing, well-used XML interchange format for molecular interaction data to meet new use cases and enable the capture of new data types, following extensive community consultation. PSI-MI XML3.0 expands the capabilities of the format beyond simple experimental data, with a concomitant update of the tool suite which serves this format. The format has been implemented by key data producers such as the International Molecular Exchange (IMEx) Consortium of protein interaction databases and the Complex Portal. [Conclusions]: PSI-MI XML3.0 has been developed by the data producers, data users, tool developers and database providers who constitute the PSI-MI workgroup. This group now actively supports PSI-MI XML2.5 as the main interchange format for experimental data, PSI-MI XML3.0 which additionally handles more complex data types, and the simpler, tab-delimited MITAB2.5, 2.6 and 2.7 for rapid parsing and download., MD, MK, AS, JS, JH and YY were funded by BBSRC MIDAS grant (BB/L024179/1), this grant provided the funds for the design of PSI-MI XML3.0 and its implementation by the IntAct database. KVR was funded by European Commission (FP7-HEALTH-2009-242129 SyBoSS), LL by ELIXIR-IIB, the Italian Node of the European ELIXIR infrastructure, IJ was funded by Ontario Research Fund (GL2–01-030, #34876) and Canada Research Chair Program (#225404), DJL by EMBL Australia and FP7-HEALTH2011-278568, SRB and NTM by Fondation pour la Recherche Médicale (grant n° DBI20141231336) and by the French Institute of Bioinformatics (2015 call), NHC and RCL by British Heart Foundation (RG/13/5/30112), GC by the European Research Council (Grant Agreement 32274), CC was funded by the Wellcome Trust (grant numbers 103139, 063412, 203149) and LS by National Institutes of Health.

Summary InterMine is a data integration warehouse and analysis software system developed for large and complex biological data sets. Designed for integrative analysis, it can be accessed through a user-friendly web interface. For bioinformaticians, extensive web services as well as programming interfaces for most common scripting languages support access to all features. The web interface includes a useful identifier look-up system, and both simple and sophisticated search options. Interactive results tables enable exploration, and data can be filtered, summarized, and browsed. A set of graphical analysis tools provide a rich environment for data exploration including statistical enrichment of sets of genes or other entities. InterMine databases have been developed for the major model organisms, budding yeast, nematode worm, fruit fly, zebrafish, mouse, and rat together with a newly developed human database. Here, we describe how this has facilitated interoperation and development of cross-organism analysis tools and reports. InterMine as a data exploration and analysis tool is also described. All the InterMine-based systems described in this article are resources freely available to the scientific community. genesis 53:547–560, 2015. © 2015 Wiley Periodicals, Inc.

Bdelloid rotifers are microscopic invertebrates that have existed for millions of years apparently without sex or meiosis. They inhabit a variety of temporary and permanent freshwater habitats globally, and many species are remarkably tolerant of desiccation. Bdelloids offer an opportunity to better understand the evolution of sex and recombination, but previous work has emphasized desiccation as the cause of several unusual genomic features in this group. Here, we evaluate the relative effects of asexuality and desiccation tolerance on genome evolution by comparing whole genome sequences for three bdelloid species:Adineta ricciae(desiccation tolerant),Rotaria macruraandRotaria magnacalcarata(both desiccation intolerant) to the only published bdelloid genome to date, that ofAdineta vaga(also desiccation tolerant). We find that tetraploidy is conserved among all four bdelloid species, but homologous divergence in obligately aquaticRotariagenomes is low, well within the range observed between alleles in obligately sexual, diploid animals. In addition, we find that homologous regions inA. ricciaeare largely collinear and do not form palindromic repeats as observed in the publishedA. vagaassembly. These findings are contrary to current understanding of the role of desiccation in shaping the bdelloid genome, and indicate that various features interpreted as genomic evidence for long-term ameiotic evolution are not general to all bdelloid species, even within the same genus. Finally, we substantiate previous findings of high levels of horizontally transferred non-metazoan genes encoded in both desiccating and non-desiccating bdelloid species, and show that this is a unique feature of bdelloids among related animal phyla. Comparisons within bdelloids and to other desiccation-tolerant animals, however, again call into question the purported role of desiccation in horizontal transfer.

Poly(A) signals located at the 3' end of eukaryotic genes drive cleavage and polyadenylation at the same end of pre-mRNA. Although these sequences are expected only at the 3' end of genes, we found that strong poly(A) signals are also predicted within the 5' untranslated regions (UTRs) of many Drosophila melanogaster mRNAs. Most of these 5' poly(A) signals have little influence on the processing of the endogenous transcripts, but they are very active when placed at the 3' end of reporter genes. In investigating these unexpected observations, we discovered that both these novel poly(A) signals and standard poly(A) signals become functionally silent when they are positioned close to transcription start sites in either Drosophila or human cells. This indicates that the stage when the poly(A) signal emerges from the polymerase II (Pol II) transcription complex determines whether a putative poly(A) signal is recognized as functional. The data suggest that this mechanism, which probably prevents cryptic poly(A) signals from causing premature transcription termination, depends on low Ser2 phosphorylation of the C-terminal domain of Pol II and inefficient recruitment of processing factors.

The anamniote lateral line system, comprising mechanosensory neuromasts and electrosensory ampullary organs, is a useful model for investigating the developmental and evolutionary diversification of different organs and cell types. Zebrafish neuromast development is increasingly well understood, but neither zebrafish nor Xenopus is electroreceptive and our molecular understanding of ampullary organ development is rudimentary. We have used RNA-seq to generate a lateral line-enriched gene-set from late-larval paddlefish (Polyodon spathula). Validation of a subset reveals expression in developing ampullary organs of transcription factor genes critical for hair cell development, and genes essential for glutamate release at hair cell ribbon synapses, suggesting close developmental, physiological and evolutionary links between non-teleost electroreceptors and hair cells. We identify an ampullary organ-specific proneural transcription factor, and candidates for the voltage-sensing L-type Ca channel and rectifying K channel predicted from skate (cartilaginous fish) ampullary organ electrophysiology. Overall, our results illuminate ampullary organ development, physiology and evolution.

Exosomes derived from all nephron segments are present in human urine, where their functionality is incompletely understood. Most studies have focused on biomarker discovery rather than exosome function. Through sequencing we identified the miRNA repertoire of urinary exosomes from healthy volunteers; 276 mature miRNAs and 345 pre-miRNAs were identified (43%/7% of reads). Among the most abundant were members of the miR-10, miR-30 and let-7 families. Targets for the identified miRNAs were predicted using five different databases; genes encoding membrane transporters and their regulators were enriched, highlighting the possibility that these miRNAs could modulate key renal tubular functions in a paracrine manner. As proof of concept, cultured renal epithelial cells were exposed to urinary exosomes and cellular exosomal uptake was confirmed; thereafter, reduced levels of the potassium channel ROMK and kinases SGK1 and WNK1 were observed in a human collecting duct cell line, while SPAK was unaltered. In proximal tubular cells, mRNA levels of the amino acid transporter gene SLC38A2 were diminished and reflected in a significant decrement of its encoded protein SNAT2. Protein levels of the kinase SGK1 did not change. Thus we demonstrated a novel potential function for miRNA in urinary exosomes.

ThaleMine (https://apps.araport.org/thalemine/) is a comprehensive data warehouse that integrates a wide array of genomic information of the model plant Arabidopsis thaliana. The data collection currently includes the latest structural and functional annotation from the Araport11 update, the Col-0 genome sequence, RNA-seq and array expression, co-expression, protein interactions, homologs, pathways, publications, alleles, germplasm and phenotypes. The data are collected from a wide variety of public resources. Users can browse gene-specific data through Gene Report pages, identify and create gene lists based on experiments or indexed keywords, and run GO enrichment analysis to investigate the biological significance of selected gene sets. Developed by the Arabidopsis Information Portal project (Araport, https://www.araport.org/), ThaleMine uses the InterMine software framework, which builds well-structured data, and provides powerful data query and analysis functionality. The warehoused data can be accessed by users via graphical interfaces, as well as programmatically via web-services. Here we describe recent developments in ThaleMine including new features and extensions, and discuss future improvements. InterMine has been broadly adopted by the model organism research community including nematode, rat, mouse, zebrafish, budding yeast, the modENCODE project, as well as being used for human data. ThaleMine is the first InterMine developed for a plant model. As additional new plant InterMines are developed by the legume and other plant research communities, the potential of cross-organism integrative data analysis will be further enabled.

Despite half a century of research, the biology of dinoflagellates remains enigmatic: they defy many functional and genetic traits attributed to typical eukaryotic cells. Genomic approaches to study dinoflagellates are often stymied due to their large, multi-gigabase genomes. Members of the genus Symbiodinium are photosynthetic endosymbionts of stony corals that provide the foundation of coral reef ecosystems. Their smaller genome sizes provide an opportunity to interrogate evolution and functionality of dinoflagellate genomes and endosymbiosis. We sequenced the genome of the ancestral Symbiodinium microadriaticum and compared it to the genomes of the more derived Symbiodinium minutum and Symbiodinium kawagutii and eukaryote model systems as well as transcriptomes from other dinoflagellates. Comparative analyses of genome and transcriptome protein sets show that all dinoflagellates, not only Symbiodinium, possess significantly more transmembrane transporters involved in the exchange of amino acids, lipids, and glycerol than other eukaryotes. Importantly, we find that only Symbiodinium harbor an extensive transporter repertoire associated with the provisioning of carbon and nitrogen. Analyses of these transporters show species-specific expansions, which provides a genomic basis to explain differential compatibilities to an array of hosts and environments, and highlights the putative importance of gene duplications as an evolutionary mechanism in dinoflagellates and Symbiodinium.

Systems biology approaches to the integrative study of cells, organs and organisms offer the best means of understanding in a holistic manner the diversity of molecular assays that can be now be implemented in a high throughput manner. Such assays can sample the genome, epigenome, proteome, metabolome and microbiome contemporaneously, allowing us for the first time to perform a complete analysis of physiological activity. The central problem remains empowering the scientific community to actually implement such an integration, across seemingly diverse data types and measurements. One promising solution is to apply semantic techniques on a self-consistent and implicitly correct ontological representation of these data types. In this paper we describe how we have applied one such solution, based around the InterMine data warehouse platform which uses as its basis the Sequence Ontology, to facilitate a systems biology analysis of virulence in the apicomplexan pathogen $Toxoplasma~gondii$, a common parasite that infects up to half the worlds population, with acute pathogenic risks for immuno-compromised individuals or pregnant mothers. Our solution, which we named `toxoMine', has provided both a platform for our collaborators to perform such integrative analyses and also opportunities for such cyberinfrastructure to be further developed, particularly to take advantage of possible semantic similarities of value to knowledge discovery in the Omics enterprise. We discuss these opportunities in the context of further enhancing the capabilities of this powerful integrative platform., To appear in the Proceedings of the 2016 IEEE Tenth International Conference on Semantic Computing (ICSC 2016)

Poly(A) signals located at the 3' end of eukaryotic genes drive cleavage and polyadenylation at the same end of pre-mRNA. Although these sequences are expected only at the 3' end of genes, we found that strong poly(A) signals are also predicted within the 5' untranslated regions (UTRs) of many Drosophila melanogaster mRNAs. Most of these 5' poly(A) signals have little influence on the processing of the endogenous transcripts, but they are very active when placed at the 3' end of reporter genes. In investigating these unexpected observations, we discovered that both these novel poly(A) signals and standard poly(A) signals become functionally silent when they are positioned close to transcription start sites in either Drosophila or human cells. This indicates that the stage when the poly(A) signal emerges from the polymerase II (Pol II) transcription complex determines whether a putative poly(A) signal is recognized as functional. The data suggest that this mechanism, which probably prevents cryptic poly(A) signals from causing premature transcription termination, depends on low Ser2 phosphorylation of the C-terminal domain of Pol II and inefficient recruitment of processing factors.

From Genome to Regulatory Networks For biologists, having a genome in hand is only the beginning—much more investigation is still needed to characterize how the genome is used to help to produce a functional organism (see the Perspective by Blaxter ). In this vein, Gerstein et al. (p. 1775 ) summarize for the Caenorhabditis elegans genome, and The modENCODE Consortium (p. 1787 ) summarize for the Drosophila melanogaster genome, full transcriptome analyses over developmental stages, genome-wide identification of transcription factor binding sites, and high-resolution maps of chromatin organization. Both studies identified regions of the nematode and fly genomes that show highly occupied targets (or HOT) regions where DNA was bound by more than 15 of the transcription factors analyzed and the expression of related genes were characterized. Overall, the studies provide insights into the organization, structure, and function of the two genomes and provide basic information needed to guide and correlate both focused and genome-wide studies.

The primary objective of the Human Genome Project was to produce high-quality sequences not just for the human genome but also for those of the chief model organisms: Escherichia coli, yeast (Saccharomyces cerevisiae), worm (Caenorhabditis elegans), fly (Drosophila melanogaster) and mouse (Mus musculus). Free access to the resultant data has prompted much biological research, including development of a map of common human genetic variants (the International HapMap Project)1, expression profiling of healthy and diseased cells2 and in-depth studies of many individual genes. These genome sequences have enabled researchers to carry out genetic and functional genomic studies not previously possible, revealing new biological insights with broad relevance across the animal kingdom 3, 4.

BACKGROUND: A fundamental concept in biology is that heritable material, DNA, is passed from parent to offspring, a process called vertical gene transfer. An alternative mechanism of gene acquisition is through horizontal gene transfer (HGT), which involves movement of genetic material between different species. HGT is well-known in single-celled organisms such as bacteria, but its existence in higher organisms, including animals, is less well established, and is controversial in humans. RESULTS: We have taken advantage of the recent availability of a sufficient number of high-quality genomes and associated transcriptomes to carry out a detailed examination of HGT in 26 animal species (10 primates, 12 flies and four nematodes) and a simplified analysis in a further 14 vertebrates. Genome-wide comparative and phylogenetic analyses show that HGT in animals typically gives rise to tens or hundreds of active 'foreign' genes, largely concerned with metabolism. Our analyses suggest that while fruit flies and nematodes have continued to acquire foreign genes throughout their evolution, humans and other primates have gained relatively few since their common ancestor. We also resolve the controversy surrounding previous evidence of HGT in humans and provide at least 33 new examples of horizontally acquired genes. CONCLUSIONS: We argue that HGT has occurred, and continues to occur, on a previously unsuspected scale in metazoans and is likely to have contributed to biochemical diversification during animal evolution.

The EBI Complex Portal is a manually curated, unifying resource of macromolecular complexes from model organisms. Each entry has a unique and stable identifier and links participating molecule to their unique reference database (UniProt, ChEBI and RNAcentral). Each complex is annotated with information about their stoichiometry, topology and structural assembly, function, complex-centric Gene Ontology terms and evidence codes. Complexes are extensively cross-referenced to ChEMBL, EMDB, Experimental Factor Ontology, Intenz, MatrixDB, the PDB and Reactome. Bespoke visualisation tools for the general topology and stoichiometry, crystal structures, molecular reactions and gene expression data are provided. All data is open-source and available in PSI-MI xml2.5 and xml3.0 standard formats, MI-JSON and tab-delineated ComplexTAB format.

Summary: The DAGViewer BioJS component is a reusable JavaScript component made available as part of the BioJS project and intended to be used to display graphs of structured data, with a particular emphasis on Directed Acyclic Graphs (DAGs). It enables users to embed representations of graphs of data, such as ontologies or phylogenetic trees, in hyper-text documents (HTML). This component is generic, since it is capable (given the appropriate configuration) of displaying any kind of data that is organised as a graph. The features of this component which are useful for examining and filtering large and complex graphs are described.Availability: http://github.com/alexkalderimis/dag-viewer-biojs; http://github.com/biojs/biojs; http://dx.doi.org/10.5281/zenodo.8303.

Position-dependent gene expression is a critical aspect of the development and behaviour of multicellular organisms. It requires a complex series of interactions to occur between different cell types in addition to intracellular signalling cascades. We used Escherichia coli to study the properties of an artificial signalling system at the interface between two expanding cell populations. We genetically engineered one population to produce a diffusible acyl-homoserine lactone (AHL) signal, and another population to respond to it. Our experiments demonstrate how such a signal can be used to reproducibly generate simple visible patterns with high accuracy in swimming agar. The producing and responding cassettes of two such signalling systems can be linked to produce a symmetric interface for bidirectional communication that can be used to visualise molecular logic. Intracellular feedback between these two cassettes would then create a framework for self-organised patterning of higher complexity. Adapting the experiments of Basu et al. (Basu et al., 2005) using cell motility, rather than a differential response to AHL concentrations as a way to define zones of response, we noted how the interaction of sender and receiver cell populations on a swimming plate could lead to complex pattern formation. Equipping highly motile strains such as E. coli MC1000 with AHL-mediated auto-inducing systems based on Vibrio fischeri luxI/luxR and Pseudomonas aeruginosa lasI/lasR cassettes would allow the amplification of a response to an AHL signal and its propagation. We designed and synthesised codon-optimised auto-inducing luxI/R and lasI/R cassettes as optimal gene expression is crucial for the generation of robust patterns. We still have to complete and test the entire genetic circuitry, although by modelling the system we were able to demonstrate its feasibility.

Coral reefs are major contributors to marine biodiversity. However, they are in rapid decline due to global environmental changes such as rising sea surface temperatures, ocean acidification, and pollution. Genomic and transcriptomic analyses have broadened our understanding of coral biology, but a study of the microRNA (miRNA) repertoire of corals is missing. miRNAs constitute a class of small non-coding RNAs of ∼22 nt in size that play crucial roles in development, metabolism, and stress response in plants and animals alike. In this study, we examined the coral Stylophora pistillata for the presence of miRNAs and the corresponding core protein machinery required for their processing and function. Based on small RNA sequencing, we present evidence for 31 bona fide microRNAs, 5 of which (miR-100, miR-2022, miR-2023, miR-2030, and miR-2036) are conserved in other metazoans. Homologues of Argonaute, Piwi, Dicer, Drosha, Pasha, and HEN1 were identified in the transcriptome of S. pistillata based on strong sequence conservation with known RNAi proteins, with additional support derived from phylogenetic trees. Examination of putative miRNA gene targets indicates potential roles in development, metabolism, immunity, and biomineralisation for several of the microRNAs. Here, we present first evidence of a functional RNAi machinery and five conserved miRNAs in S. pistillata, implying that miRNAs play a role in organismal biology of scleractinian corals. Analysis of predicted miRNA target genes in S. pistillata suggests potential roles of miRNAs in symbiosis and coral calcification. Given the importance of miRNAs in regulating gene expression in other metazoans, further expression analyses of small non-coding RNAs in transcriptional studies of corals should be informative about miRNA-affected processes and pathways. published

We have isolated and sequenced a cosmid clone from the compact genome of the Japanese pufferfish (Fugu rubripes) containing portions of three genes that have the same order as in human. The gene order is microtubule-associated protein (MAP-2), myosin light chain (MYL-1), and carbamoyl phosphate synthetase (CPS III). The intron-exon organization of Fugu CPS III is identical with that of rat CPS I, although the equivalent genomic fragments of rat and Fugu CPS span 87.9 and 21 kb, respectively. This is the first report of a piscine CPS III genomic structure and predicts a close evolutionary link between CPS III and CPS I. The 8-kb intergenic region between MYL-1 and CPS gave no clear areas of transcription factor-binding sites by pairwise comparison with shark or rat CPS promoter regions. However, there was a match with the rat myosin light chain 2 (MLC-2) gene promoter and a MyoD transcription factor-binding site 874 bp upstream of the MYL-1 gene.

We have sequenced a contiguous 284, 495-bp segment of DNA extending from the terminal (TTAGGG)n repeats of the short arm of chromosome 16, providing a full description of the transition from telomeric through subtelomeric DNA to sequences that are unique to the chromosome. To complement and extend analysis of the primary sequence, we have characterized mRNA transcripts, patterns of DNA methylation and DNase I sensitivity. Together with previous data these studies describe in detail the structural and functional organization of a human telomeric region.

The breast cancer susceptibility gene BRCA2 encodes a protein of 3418 amino acids which does not exhibit substantial sequence similarity to any other protein in the public databases. A dot matrix comparison of BRCA2 with itself revealed an eight times repeated motif in the segment of the protein encoded by exon 11. As a preliminary test of the hypothesis that these motifs are functionally significant, we have sequenced exon 11 of BRCA2 in six mammals. An alignment of the predicted protein sequences shows that, overall, the motifs have been conserved while much of the intervening sequences has diverged. These data support the notion that the BRC motifs are important in BRCA2 function. There is, however, considerable interspecies variation within certain motif units, raising the possibility of redundancy and that not all of the repeats are required for the normal function of BRCA2.

CORRIGENDUM: InterMOD: integrated data and tools for the unification of model organism research

Model organisms are widely used for understanding basic biology and have significantly contributed to the study of human disease. In recent years, genomic analysis has provided extensive evidence of widespread conservation of gene sequence and function amongst eukaryotes, allowing insights from model organisms to help decipher gene function in a wider range of species. The InterMOD consortium is developing an infrastructure based around the InterMine data warehouse system to integrate genomic and functional data from a number of key model organisms, leading the way to improved cross-species research. So far including budding yeast, nematode worm, fruit fly, zebrafish, rat and mouse, the project has set up data warehouses, synchronized data models and created analysis tools and links between data from different species. The project unites a number of major model organism databases, improving both the consistency and accessibility of comparative research, to the benefit of the wider scientific community.

The Huntington's disease (HD) gene encodes a novel protein with as yet no known function. In order to identify the functionally important domains of this protein, we have cloned and sequenced the homologue of the HD gene in the pufferfish, Fugu rubripes. The Fugu HD gene spans only 23 kb of genomic DNA, compared to the 170 kb human gene, and yet all 67 exons are conserved. The first coding exon, the site of the disease-causing triplet repeat, is highly conserved. However, the glutamine repeat in Fugu consists of just four residues. We also show that gene order may be conserved over longer stretches of the two genomes. Our work describes a detailed example of sequence comparison between human and Fugu, and illustrates the power of the pufferfish genome as a model system in the analysis of human genes.

Bdelloid rotifers are microinvertebrates with unique characteristics: they have survived tens of millions of years without sexual reproduction; they withstand extreme desiccation by undergoing anhydrobiosis; and they tolerate very high levels of ionizing radiation. Recent evidence suggests that subtelomeric regions of the bdelloid genome contain sequences originating from other organisms by horizontal gene transfer (HGT), of which some are known to be transcribed. However, the extent to which foreign gene expression plays a role in bdelloid physiology is unknown. We address this in the first large scale analysis of the transcriptome of the bdelloid Adineta ricciae: cDNA libraries from hydrated and desiccated bdelloids were subjected to massively parallel sequencing and assembled transcripts compared against the UniProtKB database by blastx to identify their putative products. Of ∼29,000 matched transcripts, ∼10% were inferred from blastx matches to be horizontally acquired, mainly from eubacteria but also from fungi, protists, and algae. After allowing for possible sources of error, the rate of HGT is at least 8%–9%, a level significantly higher than other invertebrates. We verified their foreign nature by phylogenetic analysis and by demonstrating linkage of foreign genes with metazoan genes in the bdelloid genome. Approximately 80% of horizontally acquired genes expressed in bdelloids code for enzymes, and these represent 39% of enzymes in identified pathways. Many enzymes encoded by foreign genes enhance biochemistry in bdelloids compared to other metazoans, for example, by potentiating toxin degradation or generation of antioxidants and key metabolites. They also supplement, and occasionally potentially replace, existing metazoan functions. Bdelloid rotifers therefore express horizontally acquired genes on a scale unprecedented in animals, and foreign genes make a profound contribution to their metabolism. This represents a potential mechanism for ancient asexuals to adapt rapidly to changing environments and thereby persist over long evolutionary time periods in the absence of sex., Author Summary Bdelloid rotifers are tiny invertebrates with unusual characteristics: they withstand stresses, such as desiccation and high levels of ionising radiation, that kill other animals, and they have survived over millions of years without sexual reproduction, which contradicts theories on the evolutionary advantages of sex. In this study, we investigate another bizarre feature of bdelloids, namely their ability to acquire genes from other organisms in a process known as horizontal gene transfer (HGT). We show that HGT happens on an unprecedented scale in bdelloids: approximately 10% of active genes are “foreign,” mostly originating from bacteria and other simple organisms like fungi and algae, but now functioning as bdelloid genes. About 80% of foreign genes code for enzymes, and these make a major contribution to bdelloid biochemistry: 39% of enzyme activities have a foreign contribution, and in 23% of cases the activity in question is uniquely specified by a foreign gene. This indicates biochemistry, such as toxin degradation and antioxidant generation, that is unknown in other animals and that is expected to improve the “robustness” of the bdelloid. It also represents a possible mechanism for survival without sex, by diversification of functional capacity and even replacement of defective genes by foreign counterparts.

Background BioHackathon 2010 was the third in a series of meetings hosted by the Database Center for Life Sciences (DBCLS) in Tokyo, Japan. The overall goal of the BioHackathon series is to improve the quality and accessibility of life science research data on the Web by bringing together representatives from public databases, analytical tool providers, and cyber-infrastructure researchers to jointly tackle important challenges in the area of in silico biological research. Results The theme of BioHackathon 2010 was the 'Semantic Web', and all attendees gathered with the shared goal of producing Semantic Web data from their respective resources, and/or consuming or interacting those data using their tools and interfaces. We discussed on topics including guidelines for designing semantic data and interoperability of resources. We consequently developed tools and clients for analysis and visualization. Conclusion We provide a meeting report from BioHackathon 2010, in which we describe the discussions, decisions, and breakthroughs made as we moved towards compliance with Semantic Web technologies - from source provider, through middleware, to the end-consumer.

Background Bdelloid rotifers are microscopic animals that have apparently survived without sex for millions of years and are able to survive desiccation at all life stages through a process called anhydrobiosis. Both of these characteristics are believed to have played a role in shaping several unusual features of bdelloid genomes discovered in recent years. Studies into the impact of asexuality and anhydrobiosis on bdelloid genomes have focused on understanding gene copy number. Here we investigate copy number and sequence divergence in alpha tubulin. Alpha tubulin is conserved and normally present in low copy numbers in animals, but multiplication of alpha tubulin copies has occurred in animals adapted to extreme environments, such as cold-adapted Antarctic fish. Using cloning and sequencing we compared alpha tubulin copy variation in four species of bdelloid rotifers and four species of monogonont rotifers, which are facultatively sexual and cannot survive desiccation as adults. Results were verified using transcriptome data from one bdelloid species, Adineta ricciae. Results In common with the typical pattern for animals, monogonont rotifers contain either one or two copies of alpha tubulin, but bdelloid species contain between 11 and 13 different copies, distributed across five classes. Approximately half of the copies form a highly conserved group that vary by only 1.1% amino acid pairwise divergence with each other and with the monogonont copies. The other copies have divergent amino acid sequences that evolved significantly faster between classes than within them, relative to synonymous changes, and vary in predicted biochemical properties. Copies of each class were expressed under the laboratory conditions used to construct the transcriptome. Conclusions Our findings are consistent with recent evidence that bdelloids are degenerate tetraploids and that functional divergence of ancestral copies of genes has occurred, but show how further duplication events in the ancestor of bdelloids led to proliferation in both conserved and functionally divergent copies of this gene.

Background: A central issue in genome analysis is the identification and characterization of coding regions. Estimating the coding complexity of vertebrate genomes by measuring the kinetic complexity of mRNA populations and by sequence analysis of cDNAs is limited by the fact that any given source of mRNA represents a very biased sample of all genes. Exon trapping is a method that enables the identification of genes irrespective of their transcriptional status. Results Exons were trapped from the entire mouse genome, and the resulting fragments cloned. About 7% of a random sample of exons taken from this library have significant structural homology or sequence similarity to previously sequenced genes. Using cDNAs derived from several stages of mouse development, evidence for expression of about 62% of this sample of exons was found. These data suggest that the great majority of ‘exons' in the library are derived from genes. We estimate that the fraction of the genome contained in trapped exons is 2.4%; this corresponds to a sequence complexity of about 72 megabases. Conclusion The library of exons trapped from the entire mouse genome probably represents one of the least biased and most comprehensive libraries of mouse coding regions, and should therefore prove very useful for finding genes during genome mapping and sequencing.

In an effort to comprehensively characterize the functional elements within the genomes of the important model organisms Drosophila melanogaster and Caenorhabditis elegans, the NHGRI model organism Encyclopaedia of DNA Elements (modENCODE) consortium has generated an enormous library of genomic data along with detailed, structured information on all aspects of the experiments. The modMine database (http://intermine.modencode.org) described here has been built by the modENCODE Data Coordination Center to allow the broader research community to (i) search for and download data sets of interest among the thousands generated by modENCODE; (ii) access the data in an integrated form together with non-modENCODE data sets; and (iii) facilitate fine-grained analysis of the above data. The sophisticated search features are possible because of the collection of extensive experimental metadata by the consortium. Interfaces are provided to allow both biologists and bioinformaticians to exploit these rich modENCODE data sets now available via modMine.

Background The Drosophila melanogaster genome contains 29 serpin genes, 12 as single transcripts and 17 within 6 gene clusters. Many of these serpins have a conserved "hinge" motif characteristic of active proteinase inhibitors. However, a substantial proportion (42%) lacks this motif and represents non-inhibitory serpin-fold proteins of unknown function. Currently, it is not known whether orthologous, inhibitory serpin genes retain the same target proteinase specificity within the Drosophilid lineage, nor whether they give rise to non-inhibitory serpin-fold proteins or other, more diverged, proteins. Results We collated 188 orthologues to the D. melanogaster serpins from the other 11 Drosophilid genomes and used synteny to find further family members, raising the total to 226, or 71% of the number of orthologues expected assuming complete conservation across all 12 Drosophilid species. In general the sequence constraints on the serpin-fold itself are loose. The critical Reactive Centre Loop (RCL) sequence, including the target proteinase cleavage site, is strongly conserved in inhibitory serpins, although there are 3 exceptional sets of orthologues in which the evolutionary constraints are looser. Conversely, the RCL of non-inhibitory serpin orthologues is less conserved, with 3 exceptions that presumably bind to conserved partner molecules. We derive a consensus hinge motif, for Drosophilid inhibitory serpins, which differs somewhat from that of the vertebrate consensus. Three gene clusters appear to have originated in the melanogaster subgroup, Spn28D, Spn77B and Spn88E, each containing one inhibitory serpin orthologue that is present in all Drosophilids. In addition, the Spn100A transcript appears to represent a novel serpin-derived fold. Conclusion In general, inhibitory serpins rarely change their range of proteinase targets, except by a duplication/divergence mechanism. Non-inhibitory serpins appear to derive from inhibitory serpins, but not the reverse. The conservation of different family members varied widely across the 12 sequenced Drosophilid genomes. An approach considering synteny as well as homology was important to find the largest set of orthologues.

A gene located within the intron of a larger gene is an uncommon arrangement in any species. Few of these nested gene arrangements have been explored from an evolutionary perspective. Here we report a phylogenetic analysis of kayak (kay) and fos intron gene (fig), a divergently transcribed gene located in a kay intron, utilizing 12 Drosophila species. The evolutionary relationship between these genes is of interest because kay is the homolog of the proto-oncogene c-fos whose function is modulated by serine/threonine phosphorylation and fig is a predicted PP2C phosphatase specific for serine/threonine residues. We found that, despite an extraordinary level of diversification in the intron–exon structure of kay (11 inversions and six independent exon losses), the nested arrangement of kay and fig is conserved in all species. A genomewide analysis of protein-coding nested gene pairs revealed that ∼20% of nested pairs in D. melanogaster are also nested in D. pseudoobscura and D. virilis. A phylogenetic examination of fig revealed that there are three subfamilies of PP2C phosphatases in all 12 species of Drosophila. Overall, our phylogenetic and genomewide analyses suggest that the nested arrangement of kay and fig may be due to a functional relationship between them.

This novel web-based database provides unique accessibility and querying of integrated genomic and proteomic data for Drosophila and Anopheles., FlyMine is a data warehouse that addresses one of the important challenges of modern biology: how to integrate and make use of the diversity and volume of current biological data. Its main focus is genomic and proteomics data for Drosophila and other insects. It provides web access to integrated data at a number of different levels, from simple browsing to construction of complex queries, which can be executed on either single items or lists.

We have sequenced and compared DNA from the ends of three human chromosomes: 4p, 16p and 22q. In all cases the pro-terminal regions are subdivided by degenerate (TTAGGG)n repeats into distal and proximal sub-domains with entirely different patterns of homology to other chromosome ends. The distal regions contain numerous, short (

Abstract 3525 Leukaemic transformation is driven by aberrant transcriptional programs often in combination with abnormal proliferative signalling. These programs operate in normal hematopoiesis where they are involved in hematopoietic stem cell (HSC) proliferation and maintenance. ERG is a component of normal and leukemic stem cell signatures and high ERG expression has been proposed as a risk factor for poor prognosis in acute myeloid leukemia (AML). However, mechanisms that underlie ERG expression in AML and how its expression relates to leukemic stemness are unknown. We report that ERG expression in AML is associated with activity of the ERG+85 stem cell enhancer (SCE) and a heptad of transcription factors that combinatorially regulate genes in normal HSCs. Gene expression signatures derived from ERG+85 stem cell enhancer (Fig A) and heptad activity (Fig B) predict clinical outcome in a cytogenetically normal cohort of AML (CN-AML) patients when ERG expression alone fails. The heptad signature is an independent risk factor for poor overall and event-free survival (Fig C). There were no long-term survivors amongst patients with a heptad signature, FLT3 mutations and wild-type NPM1 (Fig D) pointing to a hitherto unappreciated link between aberrant signaling and transcriptional mediators of hematopoietic stem cell identity. In two independent cohorts, the heptad signature was as closely associated with wild-type NPM1 AML as the HOX signature was with mutant NPM1 AML (Fig E–F) suggestive of a collective role for these transcription factors in mediating the leukemic signature in the former. Taken together, these results show that key transcriptional regulators cooperate in establishing stem cell signatures in leukemic cells and that the underlying spectrum of somatic mutations contributes to the development of these signatures and modulate their influence on clinical outcome. Disclosures: No relevant conflicts of interest to declare.

MicroRNA (miRNA) dysregulation has been associated with numerous cancers including breast cancer. The dysregulation of miRNAs in cancer has been shown to perturb various pathways, with oncogenic effects. Here we investigate the relationship between dysregulated miRNAs and pathways involved in breast cancer by integrating miRNA and mRNA expression data. From a list of dysregulated miRNAs, we started by selecting the subset that appear to be regulating genes differentially expressed in breast cancer vs. normal tissue. Individually and as a group, this subset was found to target several canonical oncogenic pathways including the p53 signalling pathway, MAPK signalling pathway, TGFβ signalling pathway, focal adhesion and cell cycle progression. These results suggest that the dysregulation of miRNAs in breast cancer not only results in widespread changes to gene expression, but also the dysregulation of key oncogenic pathways.

The HMR E silencer represses transcription of silent mating-type genes in the budding yeast Saccharomyces cerevisiae and contains three redundant regulatory elements A, E and B (ref. 1). The A element contains the 11 base pair consensus sequence that is essential for the firing of DNA replication origins. A multisubunit protein called the origin recognition complex (ORC) binds specifically to this consensus sequence within yeast origins in vitro and in vivo. We isolated mutants in A element-mediated silencing and report here that one of the genes we identified, RRR1, encodes ORC2, the 72K subunit of ORC. RRR1/ORC2 is an essential gene, but the rrr1-316 allele, which is viable, is defective in the replication of nuclear DNA and the maintenance of the 2-microns episomal DNA. This is, to our knowledge, the first genetic evidence that ORC is involved in DNA replication and silencing.

Repression of the yeast silent mating type loci requires cis -acting sequences located over 1 kb from the regulated promoters. One of these sites, a "silencer," exhibits enhancer-like distance- and orientation-independence. The silencer demonstrates both autonomous replication sequence (ARS) activity and a centromere-like segregation function, suggesting roles for DNA replication and segregation in transcriptional repression. Here we identify three sequences (A, E, and B) involved both in repression and in either ARS or segregation activity. The sequences are functionally redundant: no one is essential for complete transcriptional control, but mutations in any two inactivate the silencer. Surprisingly, elements E and B can each activate transcription from heterologous promoters, and E shows striking homology to several yeast upstream activation sequences. Therefore, sequences individually involved in replication, segregation, and transcriptional activation can, at the silencer, efficiently repress transcription.

SummaryStem cells have the remarkable ability to give rise to both self-renewing and differentiating daughter cells. Drosophila neural stem cells segregate cell-fate determinants from the self-renewing cell to the differentiating daughter at each division. Here, we show that one such determinant, the homeodomain transcription factor Prospero, regulates the choice between stem cell self-renewal and differentiation. We have identified the in vivo targets of Prospero throughout the entire genome. We show that Prospero represses genes required for self-renewal, such as stem cell fate genes and cell-cycle genes. Surprisingly, Prospero is also required to activate genes for terminal differentiation. We further show that in the absence of Prospero, differentiating daughters revert to a stem cell-like fate: they express markers of self-renewal, exhibit increased proliferation, and fail to differentiate. These results define a blueprint for the transition from stem cell self-renewal to terminal differentiation.

Homeotic genes confer identity to the different segments of Drosophila. These genes are expressed in many cell types over long periods of time. To determine when the homeotic genes are required for specific developmental events we have expressed the Ultrabithorax, abdominal-A and Abdominal-Bm proteins at different times during development using the GAL4 targeting technique. We find that early transient homeotic gene expression has no lasting effects on the differentiation of the larval epidermis, but it switches the fate of other cell types irreversibly (e.g. the spiracle primordia). We describe one cell type in the peripheral nervous system that makes sequential, independent responses to homeotic gene expression. We also provide evidence that supports the hypothesis of in vivo competition between the bithorax complex proteins for the regulation of their down-stream targets.