SECOM: A Novel Hash Seed and Community Detection Based-Approach for Genome-Scale Protein Domain Identification

With rapid advances in the development of DNA sequencing technologies, a plethora of high-throughput genome and proteome data from a diverse spectrum of organisms have been generated. The functional annotation and evolutionary history of proteins are usually inferred from domains predicted from the genome sequences. Traditional database-based domain prediction methods cannot identify novel domains, however, and alignment-based methods, which look for recurring segments in the proteome, are computationally demanding. Here, we propose a novel genome-wide domain prediction method, SECOM. Instead of conducting all-against-all sequence alignment, SECOM first indexes all the proteins in the genome by using a hash seed function. Local similarity can thus be detected and encoded into a graph structure, in which each node represents a protein sequence and each edge weight represents the shared hash seeds between the two nodes. SECOM then formulates the domain prediction problem as an overlapping community-finding problem in this graph. A backward graph percolation algorithm that efficiently identifies the domains is proposed. We tested SECOM on five recently sequenced genomes of aquatic animals. Our tests demonstrated that SECOM was able to identify most of the known domains identified by InterProScan. When compared with the alignment-based method, SECOM showed higher sensitivity in detecting putative novel domains, while it was also three orders of magnitude faster. For example, SECOM was able to predict a novel sponge-specific domain in nucleoside-triphosphatase (NTPases). Furthermore, SECOM discovered two novel domains, likely of bacterial origin, that are taxonomically restricted to sea anemone and hydra. SECOM is an open-source program and available at http://sfb.kaust.edu.sa/Pages/Software.aspx

A genetic screen for modifiers of Drosophila caspase Dcp-1 reveals caspase involvement in autophagy and novel caspase-related genes

BACKGROUND: Caspases are cysteine proteases with essential functions in the apoptotic pathway; their proteolytic activity toward various substrates is associated with the morphological changes of cells. Recent reports have described non-apoptotic functions of caspases, including autophagy. In this report, we searched for novel modifiers of the phenotype of Dcp-1 gain-of-function (GF) animals by screening promoter element- inserted Drosophila melanogaster lines (EP lines). RESULTS: We screened ~15,000 EP lines and identified 72 Dcp-1-interacting genes that were classified into 10 groups based on their functions and pathways: 4 apoptosis signaling genes, 10 autophagy genes, 5 insulin/IGF and TOR signaling pathway genes, 6 MAP kinase and JNK signaling pathway genes, 4 ecdysone signaling genes, 6 ubiquitination genes, 11 various developmental signaling genes, 12 transcription factors, 3 translation factors, and 11 other unclassified genes including 5 functionally undefined genes. Among them, insulin/IGF and TOR signaling pathway, MAP kinase and JNK signaling pathway, and ecdysone signaling are known to be involved in autophagy. Together with the identification of autophagy genes, the results of our screen suggest that autophagy counteracts Dcp-1-induced apoptosis. Consistent with this idea, we show that expression of eGFP-Atg5 rescued the eye phenotype caused by Dcp-1 GF. Paradoxically, we found that over-expression of full-length Dcp-1 induced autophagy, as Atg8b-GFP, an indicator of autophagy, was increased in the eye imaginal discs and in the S2 cell line. Taken together, these data suggest that autophagy suppresses Dcp-1-mediated apoptotic cell death, whereas Dcp-1 positively regulates autophagy, possibly through feedback regulation. CONCLUSIONS: We identified a number of Dcp-1 modifiers that genetically interact with Dcp-1-induced cell death. Our results showing that Dcp-1 and autophagy-related genes influence each other will aid future investigations of the complicated relationships between apoptosis and autophagy

bZIPDB : A database of regulatory information for human bZIP transcription factors

<p>Abstract</p> <p>Background</p> <p>Basic region-leucine zipper (bZIP) proteins are a class of transcription factors (TFs) that play diverse roles in eukaryotes. Malfunctions in these proteins lead to cancer and various other diseases. For detailed characterization of these TFs, further public resources are required.</p> <p>Description</p> <p>We constructed a database, designated bZIPDB, containing information on 49 human bZIP TFs, by means of automated literature collection and manual curation. bZIPDB aims to provide public data required for deciphering the gene regulatory network of the human bZIP family, e.g., evaluation or reference information for the identification of regulatory modules. The resources provided by bZIPDB include (1) protein interaction data including direct binding, phosphorylation and functional associations between bZIP TFs and other cellular proteins, along with other types of interactions, (2) bZIP TF-target gene relationships, (3) the cellular network of bZIP TFs in particular cell lines, and (4) gene information and ontology. In the current version of the database, 721 protein interactions and 560 TF-target gene relationships are recorded. bZIPDB is annually updated for the newly discovered information.</p> <p>Conclusion</p> <p>bZIPDB is a repository of detailed regulatory information for human bZIP TFs that is collected and processed from the literature, designed to facilitate analysis of this protein family. bZIPDB is available for public use at <url>http://biosoft.kaist.ac.kr/bzipdb</url>.</p

Molecular Response to Extreme Summer Temperatures Differs Between Two Genetically Differentiated Populations of a Coral Reef Fish

Extreme thermal events are increasing in frequency and duration as the climate continues to warm, with potential detrimental effects on marine organisms. However, the effects of heatwaves may differ among geographically separated populations depending on their capacity for thermal plasticity. Here, we compared the response to simulated summer heatwave temperatures (+1.5 and +3.0°C above average) in two populations of a coral reef damselfish with different capacities for thermal plasticity. We found that the more thermally tolerant population had greater plasticity of gene expression and had significantly more downregulated genes, which may provide more energy to repair damage associated with thermal stress and to maintain basic functions at these extreme temperatures. In contrast, the thermally sensitive population exhibited higher basal levels of heat shock proteins and had three times fewer changes in gene expression overall. The limited changes in gene regulation suggest that individuals have reduced genome plasticity to tolerate thermal fluctuations and consequently may not have enough energy to repair damage and resume cellular homeostasis at extreme temperatures. Thus, we have identified the molecular signatures of how two genetically distinct fish populations cope with an extreme thermal event, and why they differ in their capacity for thermal plasticity

The chromosome-scale genome assembly of the yellowtail clownfish Amphiprion clarkii provides insights into the melanic pigmentation of anemonefish

Anemonefish are an emerging group of model organisms for studying genetic, ecological, evolutionary, and developmental traits of coral reef fish. The yellowtail clownfish Amphiprion clarkii possesses species-specific characteristics such as inter-species co-habitation, high intra-species color variation, no anemone specificity, and a broad geographic distribution, that can increase our understanding of anemonefish evolutionary history, behavioral strategies, fish-anemone symbiosis, and color pattern evolution. Despite its position as an emerging model species, the genome of A. clarkii is yet to be published. Using PacBio long-read sequencing and Hi-C chromatin capture technology, we generated a high-quality chromosome-scale genome assembly initially comprised of 1, 840 contigs with an N50 of 1, 203, 211 bp. These contigs were successfully anchored into 24 chromosomes of 843, 582, 782 bp and annotated with 25, 050 protein-coding genes encompassing 97.0% of conserved actinopterygian genes, making the quality and completeness of this genome the highest among all published anemonefish genomes to date. Transcriptomic analysis identified tissue-specific gene expression patterns, with the brain and optic lobe having the largest number of expressed genes. Further analyses revealed higher copy numbers of erbb3b (a gene involved in melanocyte development) in A. clarkii compared with other anemonefish, thus suggesting a possible link between erbb3b and the natural melanism polymorphism observed in A. clarkii. The publication of this high-quality genome, along with A. clarkii's many unique traits, position this species as an ideal model organism for addressing scientific questions across a range of disciplines

First transcriptome assembly of a newly discovered vent mussel, Gigantidas vrijenhoeki, at Onnuri Vent Field on the northern Central Indian Ridge

This is the first report of a transcriptome assembly of a newly discovered hydrothermal vent mussel, Gigantidas vrijenhoeki (Bivalvia: Mytilidae), on the Central Indian Ridge. Gigantidas vrijenhoeki was identified from material collected at the newly discovered Onnuri Vent Field (OVF) on the Central Indian Ridge in 2018, and was reported as a new species, distinct from another dominant hydrothermal vent mussel, Bathymodiolus marisindicus, in 2020. We sequenced the transcriptome of G. vrijenhoeki using the Illumina HiSeq X System. De novo assembly and analysis of the coding regions predicted 25,405 genes, 84.76% of which was annotated by public databases. The transcriptome of G. vrijenhoeki will be a valuable resource in studying the ecological and biological characteristics of this new species, which is distinct from other deep-sea mussels. These data should also support the investigation of the relationship between the environmental conditions of hydrothermal vents and the unique distribution of G. vrijenhoeki in the OVF of the Central Indian Ridge

Shedding light on cell compartmentation in the candidate phylum Poribacteria by high resolution visualisation and transcriptional profiling

Assigning functions to uncultivated environmental microorganisms continues to be a challenging endeavour. Here, we present a new microscopy protocol for fluorescence in situ hybridisation-correlative light and electron microscopy (FISH-CLEM) that enabled, to our knowledge for the first time, the identification of single cells within their complex microenvironment at electron microscopy resolution. Members of the candidate phylum Poribacteria, common and uncultivated symbionts of marine sponges, were used towards this goal. Cellular 3D reconstructions revealed bipolar, spherical granules of low electron density, which likely represent carbon reserves. Poribacterial activity profiles were retrieved from prokaryotic enriched sponge metatranscriptomes using simulation-based optimised mapping. We observed high transcriptional activity for proteins related to bacterial microcompartments (BMC) and we resolved their subcellular localisation by combining FISH-CLEM with immunohistochemistry (IHC) on ultra-thin sponge tissue sections. In terms of functional relevance, we propose that the BMC-A region may be involved in 1,2-propanediol degradation. The FISH-IHC-CLEM approach was proven an effective toolkit to combine -omics approaches with functional studies and it should be widely applicable in environmental microbiology

The Draft Genome of an Octocoral, Dendronephthya gigantea

Coral reefs composed of stony corals are threatened by global marine environmental changes. However, soft coral communities of octocorallian species, appear more resilient. The genomes of several cnidarians species have been published, including from stony corals, sea anemones, and hydra. To fill the phylogenetic gap for octocoral species of cnidarians, we sequenced the octocoral, Dendronephthya gigantea, a nonsymbiotic soft coral, commonly known as the carnation coral. The D. gigantea genome size is similar to 276 Mb. A high-quality genome assembly was constructed from PacBio long reads (29.85 Gb with 108x coverage) and Illumina short paired-end reads (35.54 Gb with 128x coverage) resulting in the highest N50 value (1.4 Mb) reported thus far among cnidarian genomes. About 12% of the genome is repetitive elements and contained 28,879 predicted protein-coding genes. This gene set is composed of 94% complete BUSCO ortholog benchmark genes, which is the second highest value among the cnidarians, indicating high quality. Based on molecular phylogenetic analysis, octocoral and hexacoral divergence times were estimated at 544 MYA. There is a clear difference in Hox gene composition between these species: unlike hexacorals, the Antp superclass Evx gene was absent in D. gigantea. Here, we present the first genome assembly of a nonsymbiotic octocoral, D. gigantea to aid in the comparative genomic analysis of cnidarians, including stony and soft corals, both symbiotic and nonsymbiotic. The D. gigantea genome may also provide clues to mechanisms of differential coping between the soft and stony corals in response to scenarios of global warming