Phylogenetically informed logic relationships improve detection of biological network organization

Background: A "phylogenetic profile" refers to the presence or absence of a gene across a set of organisms, and it has been proven valuable for understanding gene functional relationships and network organization. Despite this success, few studies have attempted to search beyond just pairwise relationships among genes. Here we search for logic relationships involving three genes, and explore its potential application in gene network analyses. Results: Taking advantage of a phylogenetic matrix constructed from the large orthologs database Roundup, we invented a method to create balanced profiles for individual triplets of genes that guarantee equal weight on the different phylogenetic scenarios of coevolution between genes. When we applied this idea to LAPP, the method to search for logic triplets of genes, the balanced profiles resulted in significant performance improvement and the discovery of hundreds of thousands more putative triplets than unadjusted profiles. We found that logic triplets detected biological network organization and identified key proteins and their functions, ranging from neighbouring proteins in local pathways, to well separated proteins in the whole pathway, and to the interactions among different pathways at the system level. Finally, our case study suggested that the directionality in a logic relationship and the profile of a triplet could disclose the connectivity between the triplet and surrounding networks. Conclusion: Balanced profiles are superior to the raw profiles employed by traditional methods of phylogenetic profiling in searching for high order gene sets. Gene triplets can provide valuable information in detection of biological network organization and identification of key genes at different levels of cellular interaction

Detecting biological network organization and functional gene orthologs

SUMMARY: We developed a package TripletSearch to compute relationships within triplets of genes based on Roundup, an orthologous gene database containing >1500 genomes. These relationships, derived from the coevolution of genes, provide valuable information in the detection of biological network organization from the local to the system level, in the inference of protein functions and in the identification of functional orthologs. To run the computation, users need to provide the GI IDs of the genes of interest

Roundup 2.0: Enabling Comparative Genomics for over 1800 Genomes

Summary: Roundup is an online database of gene orthologs for over 1800 genomes, including 226 Eukaryota, 1447 Bacteria, 113 Archaea, and 21 Viruses. Orthologs are inferred using the Reciprocal Smallest Distance algorithm. Users may query Roundup for single-linkage clusters of orthologous genes based on any group of genomes. Annotated query results may be viewed in a variety of ways including as clusters of orthologs and as phylogenetic profiles. Genomic results may be downloaded in formats suitable for functional as well as phylogenetic analysis, including the recent OrthoXML standard. In addition, gene IDs can be retrieved using FASTA sequence search. All orthology results and source code are freely available

EAVN Astrometry toward the Extreme Outer Galaxy: Kinematic distance with the proper motion of G034.84-00.95

We aim to reveal the structure and kinematics of the Outer-Scutum-Centaurus (OSC) arm located on the far side of the Milky Way through very long baseline interferometry (VLBI) astrometry using KaVA, which is composed of KVN (Korean VLBI Network) and VERA (VLBI Exploration of Radio Astrometry). We report the proper motion of a 22 GHz H$_{2}$O maser source, which is associated with the star-forming region G034.84$-$00.95, to be ($\mu_{\alpha} \rm{cos}\delta$, $\mu_{\delta}$) = ($-$1.61$\pm$0.18, $-$4.29$\pm$0.16) mas yr$^{-1}$ in equatorial coordinates (J2000). We estimate the 2D kinematic distance to the source to be 18.6$\pm$1.0 kpc, which is derived from the variance-weighted average of kinematic distances with LSR velocity and the Galactic-longitude component of the measured proper motion. Our result places the source in the OSC arm and implies that G034.84$-$00.95 is moving away from the Galactic plane with a vertical velocity of $-$38$\pm$16 km s$^{-1}$. Since the H I supershell GS033+06$-$49 is located at a kinematic distance roughly equal to that of G034.84$-$00.95, it is expected that gas circulation occurs between the outer Galactic disk around G034.84$-$00.95 with a Galactocentric distance of 12.8$^{+1.0}_{-0.9}$ kpc and halo. We evaluate possible origins of the fast vertical motion of G034.84$-$00.95, which are (1) supernova explosions and (2) cloud collisions with the Galactic disk. However, neither of the possibilities are matched with the results of VLBI astrometry as well as spatial distributions of H II regions and H I gas.Comment: Accepted for publication in PASJ. 14 figures; 8 table

Trichomonas Transmembrane Cyclases Result from Massive Gene Duplication and Concomitant Development of Pseudogenes

Trichomonas vaginalis is the only medically important protist (single-cell eukaryote) that is sexually transmitted. The ∼160-Mb Trichomonas genome contains more predicted protein-encoding genes (∼60,000) than the human genome. To begin to understand why there are so many copies of some genes, we chose here to study a large family of genes encoding unique transmembrane cyclases. Our most important results include the following. More than 100 transmembrane cyclase genes do not result from chromosomal duplications, because for the most part only the coding regions of the genes, rather than flanking sequences, are duplicated. Almost half of the transmembrane cyclase genes are pseudogenes, and these pseudogenes are polymorphic among laboratory strains of Trichomonas. Messenger RNAs for numerous transmembrane cyclases are expressed simultaneously, and representative cyclase domains have adenylyl cyclase activity. In summary, the large family of Trichomonas genes encoding transmembrane adenylyl cyclases results from massive gene duplication and concomitant development of pseudogenes

Use of Giardia, which appears to have a single nucleotide-sugar transporter for UDP-GlcNAc, to identify the UDP-Glc transporter of Entamoeba

Nucleotide-sugar transporters (NSTs) transport activated sugars (e.g. UDP-GlcNAc) from the cytosol to the lumen of the endoplasmic reticulum or Golgi apparatus where they are used to make glycoproteins and glycolipids. UDP-Glc is an important component of the N-glycan-dependent quality control (QC) system for protein folding. Because Entamoeba has this QC system while Giardia does not, we hypothesized that transfected Giardia might be used to identify the UDP-Glc transporter of Entamoeba. Here we show Giardia membranes transport UDP-GlcNAc and have apyrases, which hydrolyze nucleoside-diphosphates to make the antiporter nucleoside-monophosphate. The only NST of Giardia (GlNst), which we could identify, transports UDP-GlcNAc in transfected Saccharomyces and is present in perinuclear and peripheral vesicles and increases in expression during encystation. Entamoeba membranes transport three nucleotide-sugars (UDP-Gal, UDP-GlcNAc, and UDP-Glc), and Entamoeba has three NSTs, one of which has been shown previously to transport UDP-Gal (EhNst1). Here we show recombinant EhNst2 transports UDP-Glc in transfected Giardia, while recombinant EhNst3 transports UDP-GlcNAc in transfected Saccharomyces. In summary, all three NSTs of Entamoeba and the single NST of Giardia have been molecularly characterized, and transfected Giardia provides a new system for testing heterologous UDP-Glc transporters

Giardia, Entamoeba, and Trichomonas Enzymes Activate Metronidazole (Nitroreductases) and Inactivate Metronidazole (Nitroimidazole Reductases) ▿ †

Infections with Giardia lamblia, Entamoeba histolytica, and Trichomonas vaginalis, which cause diarrhea, dysentery, and vaginitis, respectively, are each treated with metronidazole. Here we show that Giardia, Entamoeba, and Trichomonas have oxygen-insensitive nitroreductase (ntr) genes which are homologous to those genes that have nonsense mutations in metronidazole-resistant Helicobacter pylori isolates. Entamoeba and Trichomonas also have nim genes which are homologous to those genes expressed in metronidazole-resistant Bacteroides fragilis isolates. Recombinant Giardia, Entamoeba, and Trichomonas nitroreductases used NADH rather than the NADPH used by Helicobacter, and two recombinant Entamoeba nitroreductases increased the metronidazole sensitivity of transformed Escherichia coli strains. Conversely, the recombinant nitroimidazole reductases (NIMs) of Entamoeba and Trichmonas conferred very strong metronidazole resistance to transformed bacteria. The Ehntr1 gene of the genome project HM-1:IMSS strain of Entamoeba histolytica had a nonsense mutation, and the same nonsense mutation was present in 3 of 22 clinical isolates of Entamoeba. While ntr and nim mRNAs were variably expressed by cultured Entamoeba and Trichomonas isolates, there was no relationship to metronidazole sensitivity. We conclude that microaerophilic protists have bacterium-like enzymes capable of activating metronidazole (nitroreductases) and inactivating metronidazole (NIMs). While Entamoeba and Trichomonas displayed some of the changes (nonsense mutations and gene overexpression) associated with metronidazole resistance in bacteria, these changes did not confer metronidazole resistance to the microaerophilic protists examined here

The diversity of dolichol-linked precursors to Asn-linked glycans likely results from secondary loss of sets of glycosyltransferases

The vast majority of eukaryotes (fungi, plants, animals, slime mold, and euglena) synthesize Asn-linked glycans (Alg) by means of a lipid-linked precursor dolichol-PP-GlcNAc2Man9Glc3. Knowledge of this pathway is important because defects in the glycosyltransferases (Alg1-Alg12 and others not yet identified), which make dolichol-PP-glycans, lead to numerous congenital disorders of glycosylation. Here we used bioinformatic and experimental methods to characterize Alg glycosyltransferases and dolichol-PP-glycans of diverse protists, including many human pathogens, with the following major conclusions. First, it is demonstrated that common ancestry is a useful method of predicting the Alg glycosyltransferase inventory of each eukaryote. Second, in the vast majority of cases, this inventory accurately predicts the dolichol-PP-glycans observed. Third, Alg glycosyltransferases are missing in sets from each organism (e.g., all of the glycosyltransferases that add glucose and mannose are absent from Giardia and Plasmodium). Fourth, dolichol-PP-GlcNAc2Man5 (present in Entamoeba and Trichomonas) and dolichol-PP- and N-linked GlcNAc2 (present in Giardia) have not been identified previously in wild-type organisms. Finally, the present diversity of protist and fungal dolichol-PP-linked glycans appears to result from secondary loss of glycosyltransferases from a common ancestor that contained the complete set of Alg glycosyltransferases