UBCG: Up-to-date bacterial core gene set and pipeline for phylogenomic tree reconstruction

Genome-based phylogeny plays a central role in the future taxonomy and phylogenetics of Bacteria and Archaea by replacing 16S rRNA gene phylogeny. The concatenated core gene alignments are frequently used for such a purpose. The bacterial core genes are defined as single-copy, homologous genes that are present in most of the known bacterial species. There have been several studies describing such a gene set, but the number of species considered was rather small. Here we present the up-to-date bacterial core gene set, named UBCG, and software suites to accommodate necessary steps to generate and evaluate phylogenetic trees. The method was successfully used to infer phylogenomic relationship of Escherichia and related taxa and can be used for the set of genomes at any taxonomic ranks of Bacteria. The UBCG pipeline and file viewer are freely available at https://www.ezbiocloud.net/tools/ubcg and https://www.ezbiocloud.net/tools/ubcg_viewer, respectively.     

Oestrogen regulates the expression of cathepsin E-A-like gene through ER$$\upbeta $$ in liver of chicken ( Gallus gallus)

The cathepsin E-A-like, also known as ‘similar to nothepsin’, is a new member of the aspartic protease family, which may take part in processing of egg yolk macromolecules, due to it was identified in the chicken egg-yolk. Previously, studies have suggested that the expression of cathepsin E-A-like increased gradually during sexual maturation of pullets, but the exact regulation mechanism is poorly understood. In this study, to gain insight into the function and regulation mechanism of the gene in egg-laying hen, we cloned the cathepsin E-A-like gene and evaluated its evolutionary origin by using both phylogenetic and syntenic methods. The mode of the gene expression regulation was analysed through stimulating juvenile hens with \(17\upbeta \)-estradiol and chicken embryo hepatocytes with \(17\upbeta \)-estradiol combined with oestrogen receptor antagonists including MPP, ICI 182,780 and tamoxifen. Our results showed that cathepsin E-A-like was an orthologoues gene with nothepsin, which is present in birds but not in mammals. The expression of cathepsin E-A-like significantly increased in a dose-dependent manner after the juvenile hens were treated with \(17\upbeta \)-estradiol (\(P~<~0.05\)). Compared with the \(17\upbeta \)-estradiol treatment group, the expression of cathepsin E-A-like was not significantly changed when the hepatocytes were treated with \(17\upbeta \)-estradiol combined with MPP (\(P~<~0.05\)). In contrast, compared with the \(17\upbeta \)-estradiol combined with MPP treatment group, the expression of cathepsin E-A-like was significantly downregulated when the hepatocytes were treated with \(17\upbeta \)-estradiol combined with tamoxifen or ICI 182,780 (\(P~<~0.05\)). These results demonstrated that cathepsin E-A-like shared the same evolutionary origin with nothepsin. The expression of cathepsin E-A-like was regulated by oestrogen, and the regulative effect was predominantly mediated through ER-\(\upbeta \) in liver of chicken.     

<Emphasis Type="Italic">CellFateScout</Emphasis> – a bioinformatics tool for elucidating small molecule signaling pathways that drive cells in a specific direction

Background

Small molecule effects can be represented by active signaling pathways within functional networks. Identifying these can help to design new strategies to utilize known small molecules, e.g. to trigger specific cellular transformations or to reposition known drugs.

Results

We developed CellFateScout that uses the method of Latent Variables to turn differential high-throughput expression data and a functional network into a list of active signaling pathways. Applying it to Connectivity Map data, i.e., differential expression data describing small molecule effects, we then generated a Human Small Molecule Mechanisms Database. Finally, using a list of active signaling pathways as query, a similarity search can identify small molecules from the database that may trigger these pathways. We validated our approach systematically, using expression data of small molecule perturbations, yielding better predictions than popular bioinformatics tools.

Conclusions

CellFateScout can be used to select small molecules for their desired effects. The CellFateScout Cytoscape plugin, a tutorial and the Human Small Molecule Mechanisms Database are available at https://sourceforge.net/projects/cellfatescout/ under LGPLv2 license.

     

Widespread presence of "bacterial-like" PPP phosphatases in eukaryotes

Background

In eukaryotes, PPP (p rotein p hosphatase P) family is one of the two known protein phosphatase families specific for Ser and Thr. The role of PPP phosphatases in multiple signaling pathways in eukaryotic cell has been extensively studied. Unlike eukaryotic PPP phosphatases, bacterial members of the family have broad substrate specificity or may even be Tyr-specific. Moreover, one group of bacterial PPPs are diadenosine tetraphosphatases, indicating that bacterial PPP phosphatases may not necessarily function as protein phosphatases.

Results

We describe the presence in eukaryotes of three groups of expressed genes encoding "non-conventional" phosphatases of the PPP family. These enzymes are more closely related to bacterial PPP phosphatases than to the known eukaryotic members of the family. One group, found exclusively in land plants, is most closely related to PPP phosphatases from some α-Proteobacteria, including Rhizobiales, Rhodobacterales and Rhodospirillaceae. This group is therefore termed Rhi zobiales / Rh odobacterales / Rh odospirillaceae-l ike ph osphatases, or Rhilphs. Phosphatases of the other group are found in Viridiplantae, Rhodophyta, Trypanosomatidae, Plasmodium and some fungi. They are structurally related to phosphatases from psychrophilic bacteria Shewanella and Colwellia, and are termed She wanella-l ike ph osphatases, or Shelphs. Phosphatases of the third group are distantly related to ApaH, bacterial diadenosine tetraphosphatases, and are termed A paH-l ike ph osphatases, or Alphs. Patchy distribution of Alphs in animals, plants, fungi, diatoms and kinetoplasts suggests that these phosphatases were present in the common ancestor of eukaryotes but were independently lost in many lineages. Rhilphs, Shelphs and Alphs form PPP clades, as divergent from "conventional" eukaryotic PPP phosphatases as they are from each other and from major bacterial clades. In addition, comparison of primary structures revealed a previously unrecognised (I/L/V)D(S/T)G motif, conserved in all bacterial and "bacterial-like" eukaryotic PPPs, but not in "conventional" eukaryotic and archaeal PPPs.

Conclusions

Our findings demonstrate that many eukaryotes possess diverse "bacterial-like" PPP phosphatases, the enzymatic characteristics, physiological roles and precise evolutionary history of which have yet to be determined.

     

Interspecific responses of wild African carnivores to odour of 3-mercapto-3-methylbutanol,a component of wildcat and leopard urine

The scent of 3-mercapto-3-methylbutanol (3-M-3-MB), a volatile component of leopard (Panthera pardus) and domestic cat (Felis silvestris catus) urine, released at about 10 ng/s from slow-release dispensers, elicited scent-marking from African civet (Civettictis civetta), small-spotted genet (Genetta genetta) and slender mongoose (Galerella sanguinea), as well as African wildcat (F. s. cafra). A female leopard was apparently repelled by the scent. The scent-marking and scent-rubbing by species other than African wildcats and leopards were unexpected and have important implications for the design of studies to investigate chemical communication between wild mammals and the use of camera traps to estimate animal numbers. Videos showing the behaviours referred to in this article are available at; http://www.momo-p.com/showdetail-e.php?movieid=momo161223fs01a; http://www.momo-p.com/showdetail-e.php?movieid=momo161223gs01a; http://www.momo-p.com/showdetail-e.php?movieid=momo161223gg01a.     

Studies in the taxonomy of the Polypores II

The following genera are redefined:Albatrellus S. F. Gray,Heterobasidion Bref.,Haploporus Bond. et Sing. ex Sing.,Fomitopsis P. Karst. andRigidoporus Murrill two new subgenera are described:Polyporus subgen.Dendropolyporus Pouz. (type:Polyporus umbellatus) andRigidoporus subgen.Neooxyporus Pouz. (type:Polyporus latemarginatus); the genusOxyporus (Bourd. etGalz.)Donk is classified as a subgenus of the genusRigidoporus,Murrill and the generaBjerkandera P. Karst. andLeptoporus quél. are classified as subgenera of the genusTyromyces P. Karst. The new subfamilyAlbatrelloideae Pouz. (genera:Albatrellus andGrifola) is described and 14 new specific combinations are made. The new genusIrpicodon Pouz. (type:Irpex pendulus) is proposed.     

ClpXP affects the cell metabolism of <Emphasis Type="Italic">Salmonella typhimurium</Emphasis> partially in an RpoS-dependent manner

Introduction

ClpXP protease is an important proteolytic system in Salmonella enterica serovar typhimurium (S. typhimurium). Inactivation of ClpXP by deletion of clpP resulted in overproduction of RpoS and a growth defect phenotype. Only one report has indicated that deleting rpoS can restore the growth of a S. typhimurium clpP mutant to the wild-type level. Whether overproduction of RpoS is responsible for the growth deficiency resulting from clpP disruption and how ClpXP affects the cell metabolism of S. typhimurium remain to be elucidated.

Objectives

The aim of this study is to investigate the effect of ClpXP on cell metabolism of S. typhimurium and explore the possible co-effect of RpoS associated with ClpXP in cell metabolism.

Method

We constructed a clpP rpoS double deletion mutant TT-19 (ΔclpP ΔrpoS TT-1) using a two-step phage transduction technique. We then compared the metabolite fingerprints of Salmonella rpoS deletion mutant TT-14 (ΔrpoS TT-1), clpP deletion mutant TT-16 (ΔclpP TT-1), and clpP rpoS double deletion mutant TT-19 (ΔclpP ΔrpoS TT-1) with those of the wild-type strain TT-1 by using gas chromatography coupled with mass spectrometry (GC–MS).

Results

Deletion of rpoS recovered only a part of the growth of Salmonella clpP mutant. Further metabolome analysis indicated that clpP disruption changed the levels of 16 extra- and 19 intracellular substances, while the extracellular concentrations of 4 compounds (serine, l-5-oxoproline, l-glutamic acid, and l-tryptophan) and intracellular concentrations of 10 compounds (l-isoleucine, glycine, serine, l-methionine, l-phenylalanine, malic acid, citric acid, urea, putrescine, and 6-hydroxypurine) returned to their wild-type levels when rpoS was also deleted.

Conclusion

ClpXP affects the cell metabolism of S. typhimurium partially in an RpoS-dependent manner.

     

Rational design and analysis of an <Emphasis Type="Italic">Escherichia coli</Emphasis> strain for high-efficiency tryptophan production

l-tryptophan (l-trp) is a precursor of various bioactive components and has great pharmaceutical interest. However, due to the requirement of several precursors and complex regulation of the pathways involved, the development of an efficient l-trp production strain is challenging. In this study, Escherichia coli (E. coli) strain KW001 was designed to overexpress the l-trp operator sequences (trpEDCBA) and 3-deoxy-D-arabinoheptulosonate-7-phosphate synthase (aroG ^fbr ). To further improve the production of l-trp, pyruvate kinase (pykF) and the phosphotransferase system HPr (ptsH) were deleted after inactivation of repression (trpR) and attenuation (attenuator) to produce strain KW006. To overcome the relatively slow growth and to increase the transport rate of glucose, strain KW018 was generated by combinatorial regulation of glucokinase (galP) and galactose permease (glk) expression. To reduce the production of acetic acid, strain KW023 was created by repressive regulation of phosphate acetyltransferase (pta) expression. In conclusion, strain KW023 efficiently produced 39.7 g/L of l-trp with a conversion rate of 16.7% and a productivity of 1.6 g/L/h in a 5 L fed-batch fermentation system.     

Importing statistical measures into Artemis enhances gene identification in the <Emphasis Type="Italic">Leishmania</Emphasis> genome project

Background

Seattle Biomedical Research Institute (SBRI) as part of the Leishmania Genome Network (LGN) is sequencing chromosomes of the trypanosomatid protozoan species Leishmania major. At SBRI, chromosomal sequence is annotated using a combination of trained and untrained non-consensus gene-prediction algorithms with ARTEMIS, an annotation platform with rich and user-friendly interfaces.

Results

Here we describe a methodology used to import results from three different protein-coding gene-prediction algorithms (GLIMMER, TESTCODE and GENESCAN) into the ARTEMIS sequence viewer and annotation tool. Comparison of these methods, along with the CODON USAGE algorithm built into ARTEMIS, shows the importance of combining methods to more accurately annotate the L. major genomic sequence.

Conclusion

An improvised and powerful tool for gene prediction has been developed by importing data from widely-used algorithms into an existing annotation platform. This approach is especially fruitful in the Leishmania genome project where there is large proportion of novel genes requiring manual annotation.

     

Gene Ontology annotation of the rice blast fungus, <Emphasis Type="Italic">Magnaporthe oryzae</Emphasis>

Background

Magnaporthe oryzae, the causal agent of blast disease of rice, is the most destructive disease of rice worldwide. The genome of this fungal pathogen has been sequenced and an automated annotation has recently been updated to Version 6 http://www.broad.mit.edu/annotation/genome/magnaporthe_grisea/MultiDownloads.html. However, a comprehensive manual curation remains to be performed. Gene Ontology (GO) annotation is a valuable means of assigning functional information using standardized vocabulary. We report an overview of the GO annotation for Version 5 of M. oryzae genome assembly.

Methods

A similarity-based (i.e., computational) GO annotation with manual review was conducted, which was then integrated with a literature-based GO annotation with computational assistance. For similarity-based GO annotation a stringent reciprocal best hits method was used to identify similarity between predicted proteins of M. oryzae and GO proteins from multiple organisms with published associations to GO terms. Significant alignment pairs were manually reviewed. Functional assignments were further cross-validated with manually reviewed data, conserved domains, or data determined by wet lab experiments. Additionally, biological appropriateness of the functional assignments was manually checked.

Results

In total, 6,286 proteins received GO term assignment via the homology-based annotation, including 2,870 hypothetical proteins. Literature-based experimental evidence, such as microarray, MPSS, T-DNA insertion mutation, or gene knockout mutation, resulted in 2,810 proteins being annotated with GO terms. Of these, 1,673 proteins were annotated with new terms developed for Plant-Associated Microbe Gene Ontology (PAMGO). In addition, 67 experiment-determined secreted proteins were annotated with PAMGO terms. Integration of the two data sets resulted in 7,412 proteins (57%) being annotated with 1,957 distinct and specific GO terms. Unannotated proteins were assigned to the 3 root terms. The Version 5 GO annotation is publically queryable via the GO site http://amigo.geneontology.org/cgi-bin/amigo/go.cgi. Additionally, the genome of M. oryzae is constantly being refined and updated as new information is incorporated. For the latest GO annotation of Version 6 genome, please visit our website http://scotland.fgl.ncsu.edu/smeng/GoAnnotationMagnaporthegrisea.html. The preliminary GO annotation of Version 6 genome is placed at a local MySql database that is publically queryable via a user-friendly interface Adhoc Query System.

Conclusion

Our analysis provides comprehensive and robust GO annotations of the M. oryzae genome assemblies that will be solid foundations for further functional interrogation of M. oryzae.

     

Time-consistent reconciliation maps and forbidden time travel

Background

In the absence of horizontal gene transfer it is possible to reconstruct the history of gene families from empirically determined orthology relations, which are equivalent to event-labeled gene trees. Knowledge of the event labels considerably simplifies the problem of reconciling a gene tree T with a species trees S, relative to the reconciliation problem without prior knowledge of the event types. It is well-known that optimal reconciliations in the unlabeled case may violate time-consistency and thus are not biologically feasible. Here we investigate the mathematical structure of the event labeled reconciliation problem with horizontal transfer.

Results

We investigate the issue of time-consistency for the event-labeled version of the reconciliation problem, provide a convenient axiomatic framework, and derive a complete characterization of time-consistent reconciliations. This characterization depends on certain weak conditions on the event-labeled gene trees that reflect conditions under which evolutionary events are observable at least in principle. We give an \(\mathcal {O}(|V(T)|\log (|V(S)|))\)-time algorithm to decide whether a time-consistent reconciliation map exists. It does not require the construction of explicit timing maps, but relies entirely on the comparably easy task of checking whether a small auxiliary graph is acyclic. The algorithms are implemented in C++ using the boost graph library and are freely available at https://github.com/Nojgaard/tc-recon.

Significance

The combinatorial characterization of time consistency and thus biologically feasible reconciliation is an important step towards the inference of gene family histories with horizontal transfer from orthology data, i.e., without presupposed gene and species trees. The fast algorithm to decide time consistency is useful in a broader context because it constitutes an attractive component for all tools that address tree reconciliation problems.

     

Complexity and algorithms for copy-number evolution problems

Background

Cancer is an evolutionary process characterized by the accumulation of somatic mutations in a population of cells that form a tumor. One frequent type of mutations is copy number aberrations, which alter the number of copies of genomic regions. The number of copies of each position along a chromosome constitutes the chromosome’s copy-number profile. Understanding how such profiles evolve in cancer can assist in both diagnosis and prognosis.

Results

We model the evolution of a tumor by segmental deletions and amplifications, and gauge distance from profile \(\mathbf {a}\) to \(\mathbf {b}\) by the minimum number of events needed to transform \(\mathbf {a}\) into \(\mathbf {b}\). Given two profiles, our first problem aims to find a parental profile that minimizes the sum of distances to its children. Given k profiles, the second, more general problem, seeks a phylogenetic tree, whose k leaves are labeled by the k given profiles and whose internal vertices are labeled by ancestral profiles such that the sum of edge distances is minimum.

Conclusions

For the former problem we give a pseudo-polynomial dynamic programming algorithm that is linear in the profile length, and an integer linear program formulation. For the latter problem we show it is NP-hard and give an integer linear program formulation that scales to practical problem instance sizes. We assess the efficiency and quality of our algorithms on simulated instances.

Availability

https://github.com/raphael-group/CNT-ILP

     

Flechtensystematische Studien IV. Die GattungGyalidea Lett

Im Artikel wird die ursprüngliche Gattungsdiagnose derGyalidea Lett. um die Merkmale der Pilzkomponente ergänzt. Die angenommenen verwandtschaftlichen Beziehungen derGyalidea zu der FamilieAsterothyriaceae werden diskutiert. Alle bisher festgestellten Arten mit Übersichten der revidierten Proben werden kurz beschrieben. Neue Taxa und Kombinationen:Gyalidea dodgei spec. n.,G. epiphylla spec. n.,G. mayaguezensis spec. n.,G. portoricensis spec. n.,G. fritzei (Stein) comb. n.,G. fritzei var.rivularis (Eitn.) comb. n.,G. hyalinescens (Nyl.) comb. n.,G. lecideopsis var.convarians (Nyl.) comb. n.,G. lecideopsis var.stigmatoides (Nyl.) comb. n.,G. mexicana (B. de Lesd.) comb. n. undG. subscutellaris (Vězda) comb. n.     

Mating and aggressive behaviour of <Emphasis Type="Italic">Brachyrhaphis olomina</Emphasis> (Cyprinodontiformes: Poeciliidae)

Despite the increasing interest in the use of intromittent male genitalia and coercive mating behaviour in poeciliids, detailed studies of the mating behaviour of most species in this family are lacking. We describe here the mating and aggressive behaviours of Brachyrhaphis olomina, and correlate them with the condition of the female’s ovum and embryos (immature, mature and pregnant). B. olomina performed a wide range of aggressive (sidle spread, tail beating, coordinate) and mating behaviours (approximation, touch, lateral display, touch-lateral display). Some behaviours (e.g. tail beating) are shared with other poeciliids, but two sexual behaviours (touch and lateral display) and one aggressive (coordinate) behaviour may be unique to B. olomina and were not reported in a previous study. Differences in male behaviour when paired with a female with mature ovum (more mating displays, no agonistic movements) suggest that males detect the female’s reproductive condition from some distance. The distinctive nature of mating behaviour in B. olomina highlights the importance of studying different species to have a better understanding of the evolution of mating and aggressive behaviours in poeciliids. Digital video images related to the article are available at http://www.momo-p.com/showdetail-e.php?movieid=momo170720bo01a, http://www.momo-p.com/showdetail-e.php?movieid=momo170720bo02a and http://www.momo-p.com/showdetail-e.php?movieid=momo170720bo03a.     

The variability of critical species ofLotus L. in Iran and in the neighbouring countries I. the circle ofL. corniculatus L. andL. gebelia Vent.

The region of Iran, Iraq, Afghanistan and the neighbouring countries is important for some groups of the speciesLotus L., especially those of the circle ofL. corniculatus L. andL. gebelia Vent. The first group is represented by the speciesL. corniculatus L. with 4 subspecies (3 of which are important for this region), andL. tenuis Waldst. etKit. which here attains the eastern boundary of the continuous area of distribution, and by the eastern speciesL. krylovii Schischk. etSerg. andL. rechingeri Chrtková-?ertová. The second group is represented by the speciesL. gebelia Vent.,L. michauxianus Ser. in DC. andL. libanoticus Boiss. their areas of distribution covering mostly those regions. Most of the species show considerable variability within the species.     

Differentiation of species in Gonioloboceras

The type specimen ofGonioloboceras goniolobum (Meek), rediscovered by Spath in the British Museum, is the foundation for a more accurate comparative study of this and other species ofGonioloboceras.Gonioloboceras described asG. goniolobum byElias in 1938 is differentiated asGonioloboceras schmidti, new species. Suture sets (new term) for several growth stages inG. goniolobum (Meek),G. welleriSmith,G. schmidtiElias, G.eliasiMiller &Owen, andG. asiaticumLibrovitch are assembled and used for differentiation of the species.The Kazakhstan goniatite faunule containingG. asiaticum is considered of very late Pennsylvanian age.