DoMosaics: software for domain arrangement visualization and domain-centric analysis of proteins

Effect of temperature on circadian clock functioning of trees in the context of global warming

Plant survival in a warmer world requires the timely adjustment of biological processes to cyclical changes in the new environment. Circadian oscillators have been proposed to contribute to thermal adaptation and plasticity. However, the influence of temperature on circadian clock performance and its impact on plant behaviour in natural ecosystems are not well-understood. We combined bioinformatics, molecular biology and ecophysiology to investigate the effects of increasing temperatures on the functioning of the circadian clock in two closely related tree species from Patagonian forests that constitute examples of adaptation to different thermal environments based on their altitudinal profiles. Nothofagus pumilio, the species from colder environments, showed a major rearrangement of its transcriptome and reduced ability to maintain rhythmicity at high temperatures compared with Nothofagus obliqua, which inhabits warmer zones. In altitude-swap experiments, N. pumilio, but not N. obliqua, showed limited oscillator function in warmer zones of the forest, and reduced survival and growth. Our findings show that interspecific differences in the influence of temperature on circadian clock performance are associated with preferred thermal niches, and to thermal plasticity of seedlings in natural environments, highlighting the potential role of a resonating oscillator in ecological adaptation to a warming environment.

Genome assembly of the temporary socially parasitic spiny ant Polyrhachis lamellidens and its host Camponotus japonicus

Polyrhachis lamellidens is a temporary socially parasitic ant. The newly mated P. lamellidens queen takes over a colony of several Camponotus ant species and uses the labour of the host workers in the early stages of social parasitism. To facilitate genomic resources for these species, we assembled and annotated the chromosomal genome of P. lamellidens using the 10× Genomics linked-read and Hi-C sequencing, and the draft genome of its host, Camponotus japonicus, using long-read sequencing with the Revio system. The P. lamellidens chromosomal genome assembly is 214.1 Mb, 95.5% BUSCO completeness, and contains 13,703 protein-coding genes. The C. japonicus draft genome assembly is 314.2 Mb, 99.0% BUSCO completeness, and contains 11,207 protein-coding genes. Genome-wide phylogeny and synteny analysis confirmed the phylogenetic position of P. lamellidens and C. japonicus, and a high level of synteny with the genome of both ant species. In addition, P. lamellidens possesses nearly identical chemosensory proteins to its host, C. japonicus, and these genes tended to exhibit higher expression levels in the newly mated queen. The genome assemblies of P. lamellidens and its host C. japonicus provide a valuable resource for the molecular biological and bioinformatic basis for studying the strategy of social parasitism in ants.

Cellulolytic enzymes in Microbulbifer sp. Strain GL-2, a marine fish intestinal bacterium, with emphasis on endo-1,4-β-glucanases Cel5A and Cel8

Cellulose is an abundant biomass on the planet. Various cellulases from environmental microbes have been explored for industrial use of cellulose. Marine fish intestine is of interest as one source of new enzymes. Here, we report the discovery of genes encoding two β-glucosidases (Bgl3A and Bgl3B) and four endo-1,4-β-glucanases (Cel5A, Cel8, Cel5B, and Cel9) as part of the genome sequence of a cellulolytic marine bacterium, Microbulbifer sp. Strain GL-2. Five of these six enzymes (excepting Cel5B) are presumed to localize to the periplasm or outer membrane. Transcriptional analysis demonstrated that all six genes were highly expressed in stationary phase. The transcription was induced by cello-oligosaccharides rather than by glucose, suggesting that the cellulases are produced primarily for nutrient acquisition following initial growth, facilitating the secondary growth phase. We cloned the genes encoding two of the endo-1,4-β-glucanases, Cel5A and Cel8, and purified the corresponding recombinant enzymes following expression in Escherichia coli. The activity of Cel5A was observed across a wide range of temperatures (10-40 ˚C) and pHs (6-8). This pattern differed from those of Cel8 and the commercial cellulase Enthiron, both of which exhibit decreased activities below 30 ˚C and at alkaline pHs. These characteristics suggest that Cel5A might find use in industrial applications. Overall, our results reinforce the hypothesis that marine bacteria remain a possible source of novel cellulolytic activities.

CANDy: Automated analysis of domain architectures in carbohydrate-active enzymes

Carbohydrate-active enzymes (CAZymes) can be found in all domains of life and play a crucial role in metabolic and physiological processes. CAZymes often possess a modular structure, comprising not only catalytic domains but also associated domains such as carbohydrate-binding modules (CBMs) and linker domains. By exploring the modular diversity of CAZy families, catalysts with novel properties can be discovered and further insight in their biological functions and evolutionary relationships can be obtained. Here we present the carbohydrate-active enzyme domain analysis tool (CANDy), an assembly of several novel scripts, tools and databases that allows users to analyze the domain architecture of all protein sequences in a given CAZy family. CANDy's usability is shown on glycoside hydrolase family 48, a small yet underexplored family containing multi-domain enzymes. Our analysis reveals the existence of 35 distinct domain assemblies, including eight known architectures, with the remaining assemblies awaiting characterization. Moreover, we substantiate the occurrence of horizontal gene transfer from prokaryotes to insect orthologs and provide evidence for the subsequent removal of auxiliary domains, likely through a gene fission event. CANDy is available at https://github.com/PyEED/CANDy.

Thermotolerance Mechanism of Fungal GH6 Cellobiohydrolase. Part II. Structural Analysis of Thermotolerant Mutant from the Basidiomycete Phanerochaete chrysosporium

Glycoside hydrolase family 6 cellobiohydrolase (GH6 CBH) is a group of cellulases capable of hydrolyzing crystalline cellulose. However, the synergistic reaction of GH6 CBH with other cellulases is hindered by its relatively low thermotolerance. We previously obtained a thermotolerant double mutant, C240S/C393S, of GH6 CBH from the basidiomycete Phanerochaete chrysosporium (PcCel6A) by replacing the two free cysteine (Cys) residues, C240 and C393, with serine (Yamaguchi et al., J Appl Glycosci. 2020; 67;79-86). In the accompanying paper (Part I; Yamaguchi et al., J Appl Glycosci. 2024; 71: 55-62), we measured the temperature dependence of the activity and folding of C240S/C393S and its single mutants, C240S and C393S, and found that replacement of C393 was the major contributor to the increased thermotolerance of C240S/C393S. Here, in order to investigate the mechanism involved, we crystallized the wild-type and the mutant enzymes and compared their X-ray crystal structures. The overall structures of the wild-type and the three mutant enzymes were similar. However, C240S/C393S had the lowest relative B-factor at both the N-terminal loop (residues 172-177) and the C-terminal loop (residues 390-425). This result suggests that reduced structural fluctuation of the substrate-enclosing loops, possibly due to stronger hydrogen bonding involving C393, could account for the increased thermotolerance of C240S/C393S.

A DELAY OF GERMINATION 1 (DOG1)-like protein regulates spore germination in the moss Physcomitrium patens

DELAY OF GERMINATION 1 is a key regulator of dormancy in flowering plants before seed germination. Bryophytes develop haploid spores with an analogous function to seeds. Here, we investigate whether DOG1 function during germination is conserved between bryophytes and flowering plants and analyse the underlying mechanism of DOG1 action in the moss Physcomitrium patens. Phylogenetic and in silico expression analyses were performed to identify and characterise DOG1 domain-containing genes in P. patens. Germination assays were performed to characterise a Ppdog1-like1 mutant, and replacement with AtDOG1 was carried out. Yeast two-hybrid assays were used to test the interaction of the PpDOG1-like protein with DELLA proteins from P. patens and A. thaliana. P. patens possesses nine DOG1 domain-containing genes. The DOG1-like protein PpDOG1-L1 (Pp3c3_9650) interacts with PpDELLAa and PpDELLAb and the A. thaliana DELLA protein AtRGA in yeast. Protein truncations revealed the DOG1 domain as necessary and sufficient for interaction with PpDELLA proteins. Spores of Ppdog1-l1 mutant germinate faster than wild type, but replacement with AtDOG1 reverses this effect. Our data demonstrate a role for the PpDOG1-LIKE1 protein in moss spore germination, possibly alongside PpDELLAs. This suggests a conserved DOG1 domain function in germination, albeit with differential adaptation of regulatory networks in seed and spore germination.

Systematic Analysis of Diverse Polynucleotide Kinase Clp1 Family Proteins in Eukaryotes: Three Unique Clp1 Proteins of Trypanosoma brucei

The Clp1 family proteins, consisting of the Clp1 and Nol9/Grc3 groups, have polynucleotide kinase (PNK) activity at the 5' end of RNA strands and are important enzymes in the processing of some precursor RNAs. However, it remains unclear how this enzyme family diversified in the eukaryotes. We performed a large-scale molecular evolutionary analysis of the full-length genomes of 358 eukaryotic species to classify the diverse Clp1 family proteins. The average number of Clp1 family proteins in eukaryotes was 2.3 ± 1.0, and most representative species had both Clp1 and Nol9/Grc3 proteins, suggesting that the Clp1 and Nol9/Grc3 groups were already formed in the eukaryotic ancestor by gene duplication. We also detected an average of 4.1 ± 0.4 Clp1 family proteins in members of the protist phylum Euglenozoa. For example, in Trypanosoma brucei, there are three genes of the Clp1 group and one gene of the Nol9/Grc3 group. In the Clp1 group proteins encoded by these three genes, the C-terminal domains have been replaced by unique characteristics domains, so we designated these proteins Tb-Clp1-t1, Tb-Clp1-t2, and Tb-Clp1-t3. Experimental validation showed that only Tb-Clp1-t2 has PNK activity against RNA strands. As in this example, N-terminal and C-terminal domain replacement also contributed to the diversification of the Clp1 family proteins in other eukaryotic species. Our analysis also revealed that the Clp1 family proteins in humans and plants diversified through isoforms created by alternative splicing.

FAS: assessing the similarity between proteins using multi-layered feature architectures

MOTIVATION

Protein sequence comparison is a fundamental element in the bioinformatics toolkit. When sequences are annotated with features such as functional domains, transmembrane domains, low complexity regions or secondary structure elements, the resulting feature architectures allow better informed comparisons. However, many existing schemes for scoring architecture similarities cannot cope with features arising from multiple annotation sources. Those that do fall short in the resolution of overlapping and redundant feature annotations.

RESULTS

Here, we introduce FAS, a scoring method that integrates features from multiple annotation sources in a directed acyclic architecture graph. Redundancies are resolved as part of the architecture comparison by finding the paths through the graphs that maximize the pair-wise architecture similarity. In a large-scale evaluation on more than 10 000 human-yeast ortholog pairs, architecture similarities assessed with FAS are consistently more plausible than those obtained using e-values to resolve overlaps or leaving overlaps unresolved. Three case studies demonstrate the utility of FAS on architecture comparison tasks: benchmarking of orthology assignment software, identification of functionally diverged orthologs, and diagnosing protein architecture changes stemming from faulty gene predictions. With the help of FAS, feature architecture comparisons can now be routinely integrated into these and many other applications.

AVAILABILITY AND IMPLEMENTATION

FAS is available as python package: https://pypi.org/project/greedyFAS/.

Reference Genome Sequences of the Oriental Armyworm, Mythimna separata (Lepidoptera: Noctuidae)

Lepidopteran insects are an important group of animals, including those used as biochemical and physiological model species in the insect and silk industries as well as others that are major agricultural pests. Therefore, the genome sequences of several lepidopteran insects have been reported. The oriental armyworm, Mythimna separata, is an agricultural pest commonly used to study insect immune reactions and interactions with parasitoid wasps as hosts. To improve our understanding of these research topics, reference genome sequences were constructed in the present study. Using long-read and short-read sequence data, de novo assembly and polishing were performed and haplotigs were purged. Subsequently, gene predictions and functional annotations were performed. To search for orthologs of the Toll and Immune Deficiency (IMD) pathways and for C-type lectins, annotation data analysis, BLASTp, and Hummer scans were performed. The M. separata genome is 682 Mbp; its contig N50 was 2.7 Mbp, with 21,970 genes and 24,452 coding sites predicted. All orthologs of the core components of the Toll and IMD pathways and 105 C-type lectins were identified. These results suggest that the genome data were of sufficient quality for use as reference genome data and could contribute to promoting M. separata and lepidopteran research at the molecular and genome levels.

Molecular Evidence of Chemical Disguise by the Socially Parasitic Spiny Ant Polyrhachis lamellidens (Hymenoptera: Formicidae) When Invading a Host Colony

While most ant species establish a colony independently, some socially parasitic ants build the foundation of their colony by invading other ant (host) colonies and utilizing their labor forces. Many socially parasitic ants disguise their cuticular hydrocarbon (CHC) profile, which is also known as signature mixture for nestmate discrimination, when invading the host colony. Since the strategy of chemical disguise is widespread in socially parasitic ants, elucidating the mechanism of chemical disguise will promote knowledge on the evolutionary history of social parasitism. However, detailed knowledge is still lacking, as the relevant information has only originated from circumstantial evidence, which was obtained from ecological observations. In this study, we investigated the mechanism of chemical disguise in a new queen of a temporary socially parasitic spiny ant (Polyrhachis lamellidens) by measuring its CHC profile, performing a tracing assay with labeled substances, and analyzing gene expression levels. First, after rubbing behavior was observed against the host workers, the CHC profile in P. lamellidens shifted to pronounced peaks that closely resembling that of the host workers. We also observed a reduction in aggressive behaviors by the host ant against P. lamellidens after rubbing behavior was performed. In addition, P. lamellidens acquired artificially-applied labeling substances from host workers through their rubbing behaviors, while gene expression profiling showed the expression of CHC synthesis-related genes did not change during this behavior. These results suggest that P. lamellidens directly obtains host CHCs through rubbing behavior, and these host CHCs enables P. lamellidens to remain disguised during colony invasion.

Broad thermal spectrum metagenomic laccase with action for dye decolorization and fentin hydroxide treatment

Laccases are multicopper oxidases that act on various phenolic and non-phenolic compounds, enabling numerous applications including xenobiotic bioremediation, biofuel production, drug development, and cosmetic production, and they can be used as additives in the textile and food industries. This wide range of uses makes these enzymes extremely attractive for novel biotechnology applications. Here, we undertook the kinetic characterization of LacMeta, a predicted as homotrimeric (~ 107,93 kDa) small laccase, and demonstrated that this enzyme performs best at an acidic pH (pH 3–5) towards ABTS as substrate and has a broad thermal spectrum (10–60 °C), which can promote high plastic action potential through dynamic environmental temperature fluctuations. This enzyme showed following kinetic parameters: k_cat = 6.377 s⁻¹ ± 0.303, K_m = 4.219 mM, and V_max = 24.43 µM/min (against ABTS as substrate). LacMeta almost completely degraded malachite green (50 mg/mL) in only 2 h. Moreover, the enzyme was able to degrade seven dyes from four distinct classes and it respectively achieved 85% and 83% decolorization of methylene blue and trypan blue with ABTS as the mediator. In addition, LacMeta showed potential for the degradation of two thirds of an agricultural fungicide: fentin hydroxide, thus demonstrating its biotechnological aptitude for bioremediation. The results of this study suggest that LacMeta has potential in textile wastewater treatment and that it could help in the bioremediation of other human/environmental toxins such as pesticides and antibiotic compounds belonging to the same chemical classes as the degraded dyes.

LacMeta is a new two-domain laccase with activity over a wide temperature range

LacMeta maintained 50% activity after 5 months of storage at 4 °C.

Laccase was able to degrade in 2 h the Malachite Green dye, and had the potential to degrade fentin hydroxide

Exploring the multifunctionality of SR proteins

Members of the arginine–serine-rich protein family (SR proteins) are multifunctional RNA-binding proteins that have emerged as key determinants for mRNP formation, identity and fate. They bind to pre-mRNAs early during transcription in the nucleus and accompany bound transcripts until they are translated or degraded in the cytoplasm. SR proteins are mostly known for their essential roles in constitutive splicing and as regulators of alternative splicing. However, many additional activities of individual SR proteins, beyond splicing, have been reported in recent years. We will summarize the different functions of SR proteins and discuss how multifunctionality can be achieved. We will also highlight the difficulties of studying highly versatile SR proteins and propose approaches to disentangle their activities, which is transferrable to other multifunctional RBPs.

Identification of Components of the Hippo Pathway in Hydra and Potential Role of YAP in Cell Division and Differentiation

The Hippo signaling pathway has been shown to be involved in regulating cellular identity, cell/tissue size maintenance and mechanotransduction. The Hippo pathway consists of a kinase cascade which determines the nucleo-cytoplasmic localization of YAP in the cell. YAP is the effector protein in the Hippo pathway, which acts as a transcriptional cofactor for TEAD. Phosphorylation of YAP upon activation of the Hippo pathway prevents it from entering the nucleus and abrogates its function in the transcription of the target genes. In Cnidaria, the information on the regulatory roles of the Hippo pathway is virtually lacking. Here, we report the existence of a complete set of Hippo pathway core components in Hydra for the first time. By studying their phylogeny and domain organization, we report evolutionary conservation of the components of the Hippo pathway. Protein modelling suggested the conservation of YAP-TEAD interaction in Hydra. Further, we characterized the expression pattern of the homologs of yap, hippo, mob and sav in Hydra using whole-mount RNA in situ hybridization and report their possible role in stem cell maintenance. Immunofluorescence assay revealed that Hvul_YAP expressing cells occur in clusters in the body column and are excluded in the terminally differentiated regions. Actively proliferating cells marked by Ki67 exhibit YAP colocalization in their nuclei. Strikingly, a subset of these colocalized cells is actively recruited to the newly developing bud. Disruption of the YAP-TEAD interaction increased the budding rate indicating a critical role of YAP in regulating cell proliferation in Hydra. Collectively, we posit that the Hippo pathway is an essential signaling system in Hydra; its components are ubiquitously expressed in the Hydra body column and play a crucial role in Hydra tissue homeostasis.

Genome Wide Analysis of U-Box E3 Ubiquitin Ligases in Wheat (Triticum aestivum L.)

U-box E3 ligase genes play specific roles in protein degradation by post-translational modification in plant signaling pathways, developmental stages, and stress responses; however, little is known about U-box E3 genes in wheat. We identified 213 U-box E3 genes in wheat based on U-box and other functional domains in their genome sequences. The U-box E3 genes were distributed among 21 chromosomes and most showed high sequence homology with homoeologous U-box E3 genes. Synteny analysis of wheat U-box E3 genes was conducted with other plant species such as Brachypodium distachyon, barley, rice, Triricum uratu, and Aegilops tauschii. A total of 209 RNA-seq samples representing 22 tissue types, from grain, root, leaf, and spike samples across multiple time points, were analyzed for clustering of U-box E3 gene expression during developmental stages, and the genes responded differently in various tissues and developmental stages. In addition, expression analysis of U-box E3 genes under abiotic stress, including drought, heat, and both heat and drought, and cold conditions, was conducted to provide information on U-box E3 gene expression under specific stress conditions. This analysis of U-box E3 genes could provide valuable information to elucidate biological functions for a better understanding of U-box E3 genes in wheat.

Proteome-Scale Detection of Differential Conservation Patterns at Protein and Subprotein Levels with BLUR

In the multiomics era, comparative genomics studies based on gene repertoire comparison are increasingly used to investigate evolutionary histories of species, to study genotype–phenotype relations, species adaptation to various environments, or to predict gene function using phylogenetic profiling. However, comparisons of orthologs have highlighted the prevalence of sequence plasticity among species, showing the benefits of combining protein and subprotein levels of analysis to allow for a more comprehensive study of genotype/phenotype correlations. In this article, we introduce a new approach called BLUR (BLAST Unexpected Ranking), capable of detecting genotype divergence or specialization between two related clades at different levels: gain/loss of proteins but also of subprotein regions. These regions can correspond to known domains, uncharacterized regions, or even small motifs. Our method was created to allow two types of research strategies: 1) the comparison of two groups of species with no previous knowledge, with the aim of predicting phenotype differences or specializations between close species or 2) the study of specific phenotypes by comparing species that present the phenotype of interest with species that do not. We designed a website to facilitate the use of BLUR with a possibility of in-depth analysis of the results with various tools, such as functional enrichments, protein–protein interaction networks, and multiple sequence alignments. We applied our method to the study of two different biological pathways and to the comparison of several groups of close species, all with very promising results. BLUR is freely available at http://lbgi.fr/blur/.

Transcriptomic profiling of the medicinal plant Clitoria ternatea: identification of potential genes in cyclotide biosynthesis

Clitoria ternatea a perennial climber of the Fabaceae family, is well known for its agricultural and medical applications. It is also currently the only known member of the Fabaceae family that produces abundant amounts of the ultra-stable macrocyclic peptides, cyclotides, across all tissues. Cyclotides are a class of gene-encoded, disulphide-rich, macrocyclic peptides (26–37 residues) acting as defensive metabolites in several plant species. Previous transcriptomic studies have demonstrated the genetic origin of cyclotides from the Fabaceae plant family to be embedded in the albumin-1 genes, unlike its counterparts in other plant families. However, the complete mechanism of its biosynthesis and the repertoire of enzymes involved in cyclotide folding and processing remains to be understood. In this study, using RNA-Seq data and de novo transcriptome assembly of Clitoria ternatea, we have identified 71 precursor genes of cyclotides. Out of 71 unique cyclotide precursor genes obtained, 51 sequences display unique cyclotide domains, of which 26 are novel cyclotide sequences, arising from four individual tissues. MALDI-TOF mass spectrometry analysis of fractions from different tissue extracts, coupled with precursor protein sequences obtained from transcriptomic data, established the cyclotide diversity in this plant species. Special focus in this study has also been on identifying possible enzymes responsible for proper folding and processing of cyclotides in the cell. Transcriptomic mining for oxidative folding enzymes such as protein-disulphide isomerases (PDI), ER oxidoreductin-1 (ERO1) and peptidylprolyl cis-trans isomerases (PPIases)/cyclophilins, and their levels of expression are also reported. In particular, it was observed that the CtPDI genes formed plant-specific clusters among PDI genes as compared to those from other plant species. Collectively, this work provides insights into the biogenesis of the medicinally important cyclotides and establishes the expression of certain key enzymes participating in peptide biosynthesis. Also, several novel cyclotide sequences are reported and precursor sequences are analysed in detail. In the absence of a published reference genome, a comprehensive transcriptomics approach was adopted to provide an overview of diverse properties and constituents of C. ternatea.

Large-Scale Molecular Evolutionary Analysis Uncovers a Variety of Polynucleotide Kinase Clp1 Family Proteins in the Three Domains of Life

Clp1, a polyribonucleotide 5′-hydroxyl kinase in eukaryotes, is involved in pretRNA splicing and mRNA 3′-end formation. Enzymes similar in amino acid sequence to Clp1, Nol9, and Grc3, are present in some eukaryotes and are involved in prerRNA processing. However, our knowledge of how these Clp1 family proteins evolved and diversified is limited. We conducted a large-scale molecular evolutionary analysis of the Clp1 family proteins in all living organisms for which protein sequences are available in public databases. The phylogenetic distribution and frequencies of the Clp1 family proteins were investigated in complete genomes of Bacteria, Archaea and Eukarya. In total, 3,557 Clp1 family proteins were detected in the three domains of life, Bacteria, Archaea, and Eukarya. Many were from Archaea and Eukarya, but a few were found in restricted, phylogenetically diverse bacterial species. The domain structures of the Clp1 family proteins also differed among the three domains of life. Although the proteins were, on average, 555 amino acids long (range, 196–2,728), 122 large proteins with >1,000 amino acids were detected in eukaryotes. These novel proteins contain the conserved Clp1 polynucleotide kinase domain and various other functional domains. Of these proteins, >80% were from Fungi or Protostomia. The polyribonucleotide kinase activity of Thermus scotoductus Clp1 (Ts-Clp1) was characterized experimentally. Ts-Clp1 preferentially phosphorylates single-stranded RNA oligonucleotides (Km value for ATP, 2.5 µM), or single-stranded DNA at higher enzyme concentrations. We propose a comprehensive assessment of the diversification of the Clp1 family proteins and the molecular evolution of their functional domains.

Molecular signature of an ancient organizer regulated by Wnt/β-catenin signalling during primary body axis patterning in Hydra

Wnt/β-catenin signalling has been shown to play a critical role during head organizer formation in Hydra. Here, we characterized the Wnt signalling regulatory network involved in formation of the head organizer. We found that Wnt signalling regulates genes that are important in tissue morphogenesis. We identified that majority of transcription factors (TFs) regulated by Wnt/β-catenin signalling belong to the homeodomain and forkhead families. Silencing of Margin, one of the Wnt regulated homeodomain TFs, results in loss of the ectopic tentacle phenotype typically seen upon activation of Wnt signalling. Furthermore, we show that the Margin promoter is directly bound and regulated by β-catenin. Ectopic expression of Margin in zebrafish embryos results in body axis abnormalities suggesting that Margin plays a role in axis patterning. Our findings suggest that homeobox TFs came under the regulatory umbrella of Wnt/β-catenin signalling presumably resulting in the evolution of primary body axis in animal phyla.

Co-expression of synaptic genes in the sponge Amphimedon queenslandica uncovers ancient neural submodules

The synapse is a complex cellular module crucial to the functioning of neurons. It evolved largely through the exaptation of pre-existing smaller submodules, each of which are comprised of ancient sets of proteins that are conserved in modern animals and other eukaryotes. Although these ancient submodules themselves have non-neural roles, it has been hypothesized that they may mediate environmental sensing behaviors in aneural animals, such as sponges. Here we identify orthologues in the sponge Amphimedon queenslandica of genes encoding synaptic submodules in neural animals, and analyse their cell-type specific and developmental expression to determine their potential to be co-regulated. We find that genes comprising certain synaptic submodules, including those involved in vesicle trafficking, calcium-regulation and scaffolding of postsynaptic receptor clusters, are co-expressed in adult choanocytes and during metamorphosis. Although these submodules may contribute to sensory roles in this cell type and this life cycle stage, total synaptic gene co-expression profiles do not support the existence of a functional synapse in A. queenslandica. The lack of evidence for the co-regulation of genes necessary for pre- and post-synaptic functioning in A. queenslandica suggests that sponges, and perhaps the last common ancestor of sponges and other extant animals, had the ability to promulgate sensory inputs without complete synapse-like functionalities. The differential co-expression of multiple synaptic submodule genes in sponge choanocytes, which have sensory and feeding roles, however, is consistent with the metazoan ancestor minimally being able to undergo exo- and endocytosis in a controlled and localized manner.

The RNA-binding ubiquitin ligase MKRN1 functions in ribosome-associated quality control of poly(A) translation

BACKGROUND

Cells have evolved quality control mechanisms to ensure protein homeostasis by detecting and degrading aberrant mRNAs and proteins. A common source of aberrant mRNAs is premature polyadenylation, which can result in non-functional protein products. Translating ribosomes that encounter poly(A) sequences are terminally stalled, followed by ribosome recycling and decay of the truncated nascent polypeptide via ribosome-associated quality control.

RESULTS

Here, we demonstrate that the conserved RNA-binding E3 ubiquitin ligase Makorin Ring Finger Protein 1 (MKRN1) promotes ribosome stalling at poly(A) sequences during ribosome-associated quality control. We show that MKRN1 directly binds to the cytoplasmic poly(A)-binding protein (PABPC1) and associates with polysomes. MKRN1 is positioned upstream of poly(A) tails in mRNAs in a PABPC1-dependent manner. Ubiquitin remnant profiling and in vitro ubiquitylation assays uncover PABPC1 and ribosomal protein RPS10 as direct ubiquitylation substrates of MKRN1.

CONCLUSIONS

We propose that MKRN1 mediates the recognition of poly(A) tails to prevent the production of erroneous proteins from prematurely polyadenylated transcripts, thereby maintaining proteome integrity.

PhyloProfile: dynamic visualization and exploration of multi-layered phylogenetic profiles

Summary

Phylogenetic profiles form the basis for tracing proteins and their functions across species and through time. Novel genome sequences nowadays often represent species from the remotest corner of the tree of life. Thus, phylogenetic profiling becomes increasingly important for functionally annotating this data and to integrate it into a comprehensive view on organismal evolution. To strengthen the link between the sharing of a gene across species and of the corresponding function, it is meanwhile common to complement phylogenetic profiles with additional information, such as domain architecture similarities between orthologs, or pairwise similarities of other protein features. However, there are few visualization tools that facilitate an intuitive integration of these various information layers. Here, we present PhyloProfile, an R-based tool to visualize, explore and analyze multi-layered phylogenetic profiles.

Availability and implementation

PhyloProfile is available as open source code under the MIT license at https://github.com/BIONF/phyloprofile. An online version for testing PhyloProfile and for small to medium-scale analyses is available at http://applbio.biologie.uni-frankfurt.de/phyloprofile.

The Evolutionary Traceability of a Protein

Orthologs document the evolution of genes and metabolic capacities encoded in extant and ancient genomes. However, the similarity between orthologs decays with time, and ultimately it becomes insufficient to infer common ancestry. This leaves ancient gene set reconstructions incomplete and distorted to an unknown extent. Here we introduce the “evolutionary traceability” as a measure that quantifies, for each protein, the evolutionary distance beyond which the sensitivity of the ortholog search becomes limiting. Using yeast, we show that genes that were thought to date back to the last universal common ancestor are of high traceability. Their functions mostly involve catalysis, ion transport, and ribonucleoprotein complex assembly. In turn, the fraction of yeast genes whose traceability is not sufficient to infer their presence in last universal common ancestor is enriched for regulatory functions. Computing the traceabilities of genes that have been experimentally characterized as being essential for a self-replicating cell reveals that many of the genes that lack orthologs outside bacteria have low traceability. This leaves open whether their orthologs in the eukaryotic and archaeal domains have been overlooked. Looking at the example of REC8, a protein essential for chromosome cohesion, we demonstrate how a traceability-informed adjustment of the search sensitivity identifies hitherto missed orthologs in the fast-evolving microsporidia. Taken together, the evolutionary traceability helps to differentiate between true absence and nondetection of orthologs, and thus improves our understanding about the evolutionary conservation of functional protein networks. “protTrace,” a software tool for computing evolutionary traceability, is freely available at https://github.com/BIONF/protTrace.git; last accessed February 10, 2019.

Evolution of Protein Domain Architectures

This chapter reviews current research on how protein domain architectures evolve. We begin by summarizing work on the phylogenetic distribution of proteins, as this will directly impact which domain architectures can be formed in different species. Studies relating domain family size to occurrence have shown that they generally follow power law distributions, both within genomes and larger evolutionary groups. These findings were subsequently extended to multi-domain architectures. Genome evolution models that have been suggested to explain the shape of these distributions are reviewed, as well as evidence for selective pressure to expand certain domain families more than others. Each domain has an intrinsic combinatorial propensity, and the effects of this have been studied using measures of domain versatility or promiscuity. Next, we study the principles of protein domain architecture evolution and how these have been inferred from distributions of extant domain arrangements. Following this, we review inferences of ancestral domain architecture and the conclusions concerning domain architecture evolution mechanisms that can be drawn from these. Finally, we examine whether all known cases of a given domain architecture can be assumed to have a single common origin (monophyly) or have evolved convergently (polyphyly). We end by a discussion of some available tools for computational analysis or exploitation of protein domain architectures and their evolution.

Genome-wide analysis of the soybean CRK-family and transcriptional regulation by biotic stress signals triggering plant immunity

Cysteine-rich receptor-like kinases (CRKs) are transmembrane proteins that exhibit ectodomains containing the domain of unknown function 26 (DUF26). The CRKs form a large subfamily of receptor-like kinases in plants, and their possible functions remain to be elucidated. Several lines of evidence suggest that CRKs play important roles in plant defense responses to environmental stress, including plant immunity. We performed a genome-wide analysis of CRK encoding genes in soybean (Glycine max). We found 91 GmCRKs distributed in 16 chromosomes, and identified several tandem and segmental duplications, which influenced the expansion of this gene family. According to our phylogenetic analysis, GmCRKs are grouped in four clades. Furthermore, 12% of the members exhibited GmCRKs with a duplicated bi-modular organization of the ectodomains, containing four DUF26 domains. Expression analysis of GmCRKs was performed by exploring publicly available databases, and by RT-qPCR analysis of selected genes in soybean leaves responding to biotic stress signals. GmCRKs exhibited diverse expression patterns in leaves, stems, roots, and other tissues. Some of them were highly expressed in only one type of tissue, suggesting predominant roles in specific tissues. Furthermore, several GmCRKs were induced with PAMPs, DAMPs and the pathogens Phakopsora pachyrhizi and Phytophthora sojae. Expression profiles of several GmCRKs encoding highly similar proteins exhibited antagonist modes of regulation. The results suggest a fine-tuning control of GmCRKs transcriptional regulation in response to external stimuli, including PAMPs and DAMPs. This study offers a comprehensive view of the GmCRKs family in soybean, and provides a foundation for evolutionary and functional analysis of this family of plant proteins involved in the perception of pathogens and activation of plant immunity.

Evolution of asexual and sexual reproduction in the aspergilli

Aspergillus nidulans has long-been used as a model organism to gain insights into the genetic basis of asexual and sexual developmental processes both in other members of the genus Aspergillus, and filamentous fungi in general. Paradigms have been established concerning the regulatory mechanisms of conidial development. However, recent studies have shown considerable genome divergence in the fungal kingdom, questioning the general applicability of findings from Aspergillus, and certain longstanding evolutionary theories have been questioned. The phylogenetic distribution of key regulatory elements of asexual reproduction in A. nidulans was investigated in a broad taxonomic range of fungi. This revealed that some proteins were well conserved in the Pezizomycotina (e.g. AbaA, FlbA, FluG, NsdD, MedA, and some velvet proteins), suggesting similar developmental roles. However, other elements (e.g. BrlA) had a more restricted distribution solely in the Eurotiomycetes, and it appears that the genetic control of sporulation seems to be more complex in the aspergilli than in some other taxonomic groups of the Pezizomycotina. The evolution of the velvet protein family is discussed based on the history of expansion and contraction events in the early divergent fungi. Heterologous expression of the A. nidulans abaA gene in Monascus ruber failed to induce development of complete conidiophores as seen in the aspergilli, but did result in increased conidial production. The absence of many components of the asexual developmental pathway from members of the Saccharomycotina supports the hypothesis that differences in the complexity of their spore formation is due in part to the increased diversity of the sporulation machinery evident in the Pezizomycotina. Investigations were also made into the evolution of sex and sexuality in the aspergilli. MAT loci were identified from the heterothallic Aspergillus (Emericella) heterothallicus and Aspergillus (Neosartorya) fennelliae and the homothallic Aspergillus pseudoglaucus (=Eurotium repens). A consistent architecture of the MAT locus was seen in these and other heterothallic aspergilli whereas much variation was seen in the arrangement of MAT loci in homothallic aspergilli. This suggested that it is most likely that the common ancestor of the aspergilli exhibited a heterothallic breeding system. Finally, the supposed prevalence of asexuality in the aspergilli was examined. Investigations were made using A. clavatus as a representative 'asexual' species. It was possible to induce a sexual cycle in A. clavatus given the correct MAT1-1 and MAT1-2 partners and environmental conditions, with recombination confirmed utilising molecular markers. This indicated that sexual reproduction might be possible in many supposedly asexual aspergilli and beyond, providing general insights into the nature of asexuality in fungi.

Identification and comparison of key RNA interference machinery from western corn rootworm, fall armyworm, and southern green stink bug

RNA interference (RNAi)-based technology shows great potential for use in agriculture, particularly for management of costly insect pests. In the decade since the insecticidal effects of environmentally-introduced RNA were first reported, this treatment has been applied to several types of insect pests. Through the course of those efforts, it has become apparent that different insects exhibit a range of sensitivity to environmentally-introduced RNAs. The variation in responses across insect is not well-understood, with differences in the underlying RNAi mechanisms being one explanation. This study evaluates eight proteins among three agricultural pests whose responses to environmental RNAi are known to differ: western corn rootworm (Diabrotica virgifera virgifera), fall armyworm (Spodoptera frugiperda), and southern green stink bug (Nezara viridula). These proteins have been identified in various organisms as centrally involved in facilitating the microRNA- and small interfering-RNA-mediated interference responses. Various bioinformatics tools, as well as gene expression profiling, were used to identify and evaluate putative homologues for characteristics that may contribute to the differing responses of these insects, such as the absence of critical functional domains within expressed sequences, the absence of entire gene sequences, or unusually low or undetectable expression of critical genes. Though many similarities were observed, the number of isoforms and expression levels of double-stranded RNA-binding and argonaute proteins varied across insect. Differences among key RNAi machinery genes of these three pests may impact the function of their RNAi pathways, and therefore, their respective responses to exogenous RNAs.

The structure of the S-layer of Clostridium difficile

The nosocomially acquired pathogen Clostridium difficile is the primary causative agent of antibiotic associated diarrhoea and causes tens of thousands of deaths globally each year. C. difficile presents a paracrystalline protein array on the surface of the cell known as an S-layer. S-layers have been demonstrated to possess a wide range of important functions, which, combined with their inherent accessibility, makes them a promising drug target. The unusually complex S-layer of C. difficile is primarily comprised of the high- and low- molecular weight S-layer proteins, HMW SLP and LMW SLP, formed from the cleavage of the S-layer precursor protein, SlpA, but may also contain up to 28 SlpA paralogues. A model of how the S-layer functions as a whole is required if it is to be exploited in fighting the bacterium. Here, we provide a summary of what is known about the S-layer of C. difficile and each of the paralogues and, considering some of the domains present, suggest potential roles for them.

Studying AMPK in an Evolutionary Context

The AMPK protein kinase forms the heart of a complex network controlling the metabolic activities in a eukaryotic cell. Unraveling the steps by which this pathway evolved from its primordial roots in the last eukaryotic common ancestor to its present status in contemporary species has the potential to shed light on the evolution of eukaryotes. A homolog search for the proteins interacting in this pathway is considerably straightforward. However, interpreting the results, when reconstructing the evolutionary history of the pathway over larger evolutionary distances, bears a number of pitfalls. With this in mind, we present a protocol to trace a metabolic pathway across contemporary species and backward in evolutionary time. Alongside the individual analysis steps, we provide guidelines for data interpretation generalizing beyond the analysis of AMPK.

Origin and Evolution of the Sponge Aggregation Factor Gene Family

Although discriminating self from nonself is a cardinal animal trait, metazoan allorecognition genes do not appear to be homologous. Here, we characterize the Aggregation Factor (AF) gene family, which encodes putative allorecognition factors in the demosponge Amphimedon queenslandica, and trace its evolution across 24 sponge (Porifera) species. The AF locus in Amphimedon is comprised of a cluster of five similar genes that encode Calx-beta and Von Willebrand domains and a newly defined Wreath domain, and are highly polymorphic. Further AF variance appears to be generated through individualistic patterns of RNA editing. The AF gene family varies between poriferans, with protein sequences and domains diagnostic of the AF family being present in Amphimedon and other demosponges, but absent from other sponge classes. Within the demosponges, AFs vary widely with no two species having the same AF repertoire or domain organization. The evolution of AFs suggests that their diversification occurs via high allelism, and the continual and rapid gain, loss and shuffling of domains over evolutionary time. Given the marked differences in metazoan allorecognition genes, we propose the rapid evolution of AFs in sponges provides a model for understanding the extensive diversification of self-nonself recognition systems in the animal kingdom.

Aquerium: A web application for comparative exploration of domain-based protein occurrences on the taxonomically clustered genome tree

Gene duplication and loss are major driving forces in evolution. While many important genomic resources provide information on gene presence, there is a lack of tools giving equal importance to presence and absence information as well as web platforms enabling easy visual comparison of multiple domain-based protein occurrences at once. Here, we present Aquerium, a platform for visualizing genomic presence and absence of biomolecules with a focus on protein domain architectures. The web server offers advanced domain organization querying against the database of pre-computed domains for ∼26,000 organisms and it can be utilized for identification of evolutionary events, such as fusion, disassociation, duplication, and shuffling of protein domains. The tool also allows alternative inputs of custom entries or BLASTP results for visualization. Aquerium will be a useful tool for biologists who perform comparative genomic and evolutionary analyses. The web server is freely accessible at http://aquerium.utk.edu. Proteins 2016; 85:72-77. © 2016 Wiley Periodicals, Inc.

Phylogenetic profiles of all membrane transport proteins of the malaria parasite highlight new drug targets

In order to combat the on-going malaria epidemic, discovery of new drug targets remains vital. Proteins that are essential to survival and specific to malaria parasites are key candidates. To survive within host cells, the parasites need to acquire nutrients and dispose of waste products across multiple membranes. Additionally, like all eukaryotes, they must redistribute ions and organic molecules between their various internal membrane bound compartments. Membrane transport proteins mediate all of these processes and are considered important mediators of drug resistance as well as drug targets in their own right. Recently, using advanced experimental genetic approaches and streamlined life cycle profiling, we generated a large collection of Plasmodium berghei gene deletion mutants and assigned essential gene functions, highlighting potential targets for prophylactic, therapeutic, and transmission-blocking anti-malarial drugs. Here, we present a comprehensive orthology assignment of all Plasmodium falciparum putative membrane transport proteins and provide a detailed overview of the associated essential gene functions obtained through experimental genetics studies in human and murine model parasites. Furthermore, we discuss the phylogeny of selected potential drug targets identified in our functional screen. We extensively discuss the results in the context of the functional assignments obtained using gene targeting available to date.

Emergence of de novo proteins from 'dark genomic matter' by 'grow slow and moult'

Proteins are the workhorses of the cell and, over billions of years, they have evolved an amazing plethora of extremely diverse and versatile structures with equally diverse functions. Evolutionary emergence of new proteins and transitions between existing ones are believed to be rare or even impossible. However, recent advances in comparative genomics have repeatedly called some 10%-30% of all genes without any detectable similarity to existing proteins. Even after careful scrutiny, some of those orphan genes contain protein coding reading frames with detectable transcription and translation. Thus some proteins seem to have emerged from previously non-coding 'dark genomic matter'. These 'de novo' proteins tend to be disordered, fast evolving, weakly expressed but also rapidly assuming novel and physiologically important functions. Here we review mechanisms by which 'de novo' proteins might be created, under which circumstances they may become fixed and why they are elusive. We propose a 'grow slow and moult' model in which first a reading frame is extended, coding for an initially disordered and non-globular appendage which, over time, becomes more structured and may also become associated with other proteins.

TreeDom: a graphical web tool for analysing domain architecture evolution

Summary: We present TreeDom, a web tool for graphically analysing the evolutionary history of domains in multi-domain proteins. Individual domains on the same protein chain may have distinct evolutionary histories, which is important to grasp in order to understand protein function. For instance, it may be important to know whether a domain was duplicated recently or long ago, to know the origin of inserted domains, or to know the pattern of domain loss within a protein family. TreeDom uses the Pfam database as the source of domain annotations, and displays these on a sequence tree. An advantage of TreeDom is that the user can limit the analysis to N sequences that are most similar to a query, or provide a list of sequence IDs to include. Using the Pfam alignment of the selected sequences, a tree is built and displayed together with the domain architecture of each sequence.

Availablility and implementation:  http://TreeDom.sbc.su.se

Contact:  Erik.Sonnhammer@scilifelab.se

The language of the protein universe

Proteins, the main cell machinery which play a major role in nearly every cellular process, have always been a central focus in biology. We live in the post-genomic era, and inferring information from massive data sets is a steadily growing universal challenge. The increasing availability of fully sequenced genomes can be regarded as the 'Rosetta Stone' of the protein universe, allowing the understanding of genomes and their evolution, just as the original Rosetta Stone allowed Champollion to decipher the ancient Egyptian hieroglyphics. In this review, we consider aspects of the protein domain architectures repertoire that are closely related to those of human languages and aim to provide some insights about the language of proteins.

Acinetobacter baumannii Virulence Is Mediated by the Concerted Action of Three Phospholipases D

Acinetobacter baumannii causes a broad range of opportunistic infections in humans. Its success as an emerging pathogen is due to a combination of increasing antibiotic resistance, environmental persistence and adaptation to the human host. To date very little is known about the molecular basis of the latter. Here we demonstrate that A. baumannii can use phosphatidylcholine, an integral part of human cell membranes, as sole carbon and energy source. We report on the identification of three phospholipases belonging to the PLD superfamily. PLD1 and PLD2 appear restricted to the bacteria and display the general features of bacterial phospholipases D. They possess two PLDc_2 PFAM domains each encompassing the HxKx₄Dx₆GS/GGxN (HKD) motif necessary for forming the catalytic core. The third candidate, PLD3, is found in bacteria as well as in eukaryotes and harbours only one PLDc_2 PFAM domain and one conserved HKD motif, which however do not overlap. Employing a markerless mutagenesis system for A. baumannii ATCC 19606^T, we generated a full set of PLD knock-out mutants. Galleria mellonella infection studies as well as invasion experiments using A549 human lung epithelial cells revealed that the three PLDs act in a concerted manner as virulence factors and are playing an important role in host cell invasion.

The Rise and Fall of TRP-N, an Ancient Family of Mechanogated Ion Channels, in Metazoa

Mechanoreception, the sensing of mechanical forces, is an ancient means of orientation and communication and tightly linked to the evolution of motile animals. In flies, the transient-receptor-potential N protein (TRP-N) was found to be a cilia-associated mechanoreceptor. TRP-N belongs to a large and diverse family of ion channels. Its unusually long N-terminal repeat of 28 ankyrin domains presumably acts as the gating spring by which mechanical energy induces channel gating. We analyzed the evolutionary origins and possible diversification of TRP-N. Using a custom-made set of highly discriminative sequence profiles we scanned a representative set of metazoan genomes and subsequently corrected several gene models. We find that, contrary to other ion channel families, TRP-N is remarkably conserved in its domain arrangements and copy number (1) in all Bilateria except for amniotes, even in the wake of several whole-genome duplications. TRP-N is absent in Porifera but present in Ctenophora and Placozoa. Exceptional multiplications of TRP-N occurred in Cnidaria, independently along the Hydra and the Nematostella lineage. Molecular signals of subfunctionalization can be attributed to different mechanisms of activation of the gating spring. In Hydra this is further supported by in situ hybridization and immune staining, suggesting that at least three paralogs adapted to nematocyte discharge, which is key for predation and defense. We propose that these new candidate proteins help explain the sensory complexity of Cnidaria which has been previously observed but so far has lacked a molecular underpinning. Also, the ancient appearance of TRP-N supports a common origin of important components of the nervous systems in Ctenophores, Cnidaria, and Bilateria.

Domain similarity based orthology detection

Background

Orthologous protein detection software mostly uses pairwise comparisons of amino-acid sequences to assert whether two proteins are orthologous or not. Accordingly, when the number of sequences for comparison increases, the number of comparisons to compute grows in a quadratic order. A current challenge of bioinformatic research, especially when taking into account the increasing number of sequenced organisms available, is to make this ever-growing number of comparisons computationally feasible in a reasonable amount of time. We propose to speed up the detection of orthologous proteins by using strings of domains to characterize the proteins.

Results

We present two new protein similarity measures, a cosine and a maximal weight matching score based on domain content similarity, and new software, named porthoDom. The qualities of the cosine and the maximal weight matching similarity measures are compared against curated datasets. The measures show that domain content similarities are able to correctly group proteins into their families. Accordingly, the cosine similarity measure is used inside porthoDom, the wrapper developed for proteinortho. porthoDom makes use of domain content similarity measures to group proteins together before searching for orthologs. By using domains instead of amino acid sequences, the reduction of the search space decreases the computational complexity of an all-against-all sequence comparison.

Conclusion

We demonstrate that representing and comparing proteins as strings of discrete domains, i.e. as a concatenation of their unique identifiers, allows a drastic simplification of search space. porthoDom has the advantage of speeding up orthology detection while maintaining a degree of accuracy similar to proteinortho. The implementation of porthoDom is released using python and C++ languages and is available under the GNU GPL licence 3 at http://www.bornberglab.org/pages/porthoda.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0570-8) contains supplementary material, which is available to authorized users.

LacSubPred: predicting subtypes of Laccases, an important lignin metabolism-related enzyme class, using in silico approaches

Background

Laccases (E.C. 1.10.3.2) are multi-copper oxidases that have gained importance in many industries such as biofuels, pulp production, textile dye bleaching, bioremediation, and food production. Their usefulness stems from the ability to act on a diverse range of phenolic compounds such as o-/p-quinols, aminophenols, polyphenols, polyamines, aryl diamines, and aromatic thiols. Despite acting on a wide range of compounds as a family, individual Laccases often exhibit distinctive and varied substrate ranges. This is likely due to Laccases involvement in many metabolic roles across diverse taxa. Classification systems for multi-copper oxidases have been developed using multiple sequence alignments, however, these systems seem to largely follow species taxonomy rather than substrate ranges, enzyme properties, or specific function. It has been suggested that the roles and substrates of various Laccases are related to their optimal pH. This is consistent with the observation that fungal Laccases usually prefer acidic conditions, whereas plant and bacterial Laccases prefer basic conditions. Based on these observations, we hypothesize that a descriptor-based unsupervised learning system could generate homology independent classification system for better describing the functional properties of Laccases.

Results

In this study, we first utilized unsupervised learning approach to develop a novel homology independent Laccase classification system. From the descriptors considered, physicochemical properties showed the best performance. Physicochemical properties divided the Laccases into twelve subtypes. Analysis of the clusters using a t-test revealed that the majority of the physicochemical descriptors had statistically significant differences between the classes. Feature selection identified the most important features as negatively charges residues, the peptide isoelectric point, and acidic or amidic residues. Secondly, to allow for classification of new Laccases, a supervised learning system was developed from the clusters. The models showed high performance with an overall accuracy of 99.03%, error of 0.49%, MCC of 0.9367, precision of 94.20%, sensitivity of 94.20%, and specificity of 99.47% in a 5-fold cross-validation test. In an independent test, our models still provide a high accuracy of 97.98%, error rate of 1.02%, MCC of 0.8678, precision of 87.88%, sensitivity of 87.88% and specificity of 98.90%.

Conclusion

This study provides a useful classification system for better understanding of Laccases from their physicochemical properties perspective. We also developed a publically available web tool for the characterization of Laccase protein sequences (http://lacsubpred.bioinfo.ucr.edu/). Finally, the programs used in the study are made available for researchers interested in applying the system to other enzyme classes (https://github.com/tweirick/SubClPred).