Implementation of homology based and non-homology based computational methods for the identification and annotation of orphan enzymes: using Mycobacterium tuberculosis H37Rv as a case study

Functional prediction of proteins from the human gut archaeome

The human gastrointestinal tract contains diverse microbial communities, including archaea. Among them, Methanobrevibacter smithii represents a highly active and clinically relevant methanogenic archaeon, being involved in gastrointestinal disorders, such as inflammatory bowel disease and obesity. Herein, we present an integrated approach using sequence and structure information to improve the annotation of M. smithii proteins using advanced protein structure prediction and annotation tools, such as AlphaFold2, trRosetta, ProFunc, and DeepFri. Of an initial set of 873 481 archaeal proteins, we found 707 754 proteins exclusively present in the human gut. Having analysed archaeal proteins together with 87 282 994 bacterial proteins, we identified unique archaeal proteins and archaeal-bacterial homologs. We then predicted and characterized functional domains and structures of 73 unique and homologous archaeal protein clusters linked the human gut and M. smithii. We refined annotations based on the predicted structures, extending existing sequence similarity-based annotations. We identified gut-specific archaeal proteins that may be involved in defense mechanisms, virulence, adhesion, and the degradation of toxic substances. Interestingly, we identified potential glycosyltransferases that could be associated with N-linked and O-glycosylation. Additionally, we found preliminary evidence for interdomain horizontal gene transfer between Clostridia species and M. smithii, which includes sporulation Stage V proteins AE and AD. Our study broadens the understanding of archaeal biology, particularly M. smithii, and highlights the importance of considering both sequence and structure for the prediction of protein function.

An informatic workflow for the enhanced annotation of excretory/secretory proteins of Haemonchus contortus

Major advances in genomic and associated technologies have demanded reliable bioinformatic tools and workflows for the annotation of genes and their products via comparative analyses using well-curated reference data sets, accessible in public repositories. However, the accurate in silico annotation of molecules (proteins) encoded in organisms (e.g., multicellular parasites) which are evolutionarily distant from those for which these extensive reference data sets are available, including invertebrate model organisms (e.g., Caenorhabditis elegans - free-living nematode, and Drosophila melanogaster - the vinegar fly) and vertebrate species (e.g., Homo sapiens and Mus musculus), remains a major challenge. Here, we constructed an informatic workflow for the enhanced annotation of biologically-important, excretory/secretory (ES) proteins ("secretome") encoded in the genome of a parasitic roundworm, called Haemonchus contortus (commonly known as the barber's pole worm). We critically evaluated the performance of five distinct methods, refined some of them, and then combined the use of all five methods to comprehensively annotate ES proteins, according to gene ontology, biological pathways and/or metabolic (enzymatic) processes. Then, using optimised parameter settings, we applied this workflow to comprehensively annotate 2591 of all 3353 proteins (77.3%) in the secretome of H. contortus. This result is a substantial improvement (10-25%) over previous annotations using individual, "off-the-shelf" algorithms and default settings, indicating the ready applicability of the present, refined workflow to gene/protein sequence data sets from a wide range of organisms in the Tree-of-Life.

Molecular Insight into Mycobacterium tuberculosis Resistance to Nitrofuranyl Amides Gained through Metagenomics-like Analysis of Spontaneous Mutants

We performed synthesis of new nitrofuranyl amides and investigated their anti-TB activity and primary genetic response of mycobacteria through whole-genome sequencing (WGS) of spontaneous resistant mutants. The in vitro activity was assessed on reference strain Mycobacterium tuberculosis H37Rv. The most active compound 11 was used for in vitro selection of spontaneous resistant mutants. The same mutations in six genes were detected in bacterial cultures grown under increased concentrations of 11 (2×, 4×, 8× MIC). The mutant positions were presented as mixed wild type and mutant alleles while increasing the concentration of the compound led to the semi-proportional and significant increase in mutant alleles. The identified genes belong to different categories and pathways. Some of them were previously reported as mediating drug resistance or drug tolerance, and counteracting oxidative and nitrosative stress, in particular: Rv0224c, fbiC, iniA, and Rv1592c. Gene-set interaction analysis revealed a certain weak interaction for gene pairs Rv1592–Rv1639c and Rv1592–Rv0224c. To conclude, this study experimentally demonstrated a multifaceted primary genetic response of M. tuberculosis to the action of nitrofurans. All three 11-treated subcultures independently presented the same six SNPs, which suggests their non-random occurrence and likely causative relationship between compound action and possible resistance mechanism.

Identification of a novel gene required for competitive growth at high temperature in the thermotolerant yeast Kluyveromyces marxianus

It is important to understand the basis of thermotolerance in yeasts to broaden their application in industrial biotechnology. The capacity to run bioprocesses at temperatures above 40 °C is of great interest but this is beyond the growth range of most of the commonly used yeast species. In contrast, some industrial yeasts such as Kluyveromyces marxianus can grow at temperatures of 45 °C or higher. Such species are valuable for direct use in industrial biotechnology and as a vehicle to study the genetic and physiological basis of yeast thermotolerance. In previous work, we reported that evolutionarily young genes disproportionately changed expression when yeast were growing under stressful conditions and postulated that such genes could be important for long-term adaptation to stress. Here, we tested this hypothesis in K. marxianus by identifying and studying species-specific genes that showed increased expression during high-temperature growth. Twelve such genes were identified and 11 were successfully inactivated using CRISPR-mediated mutagenesis. One gene, KLMX_70384, is required for competitive growth at high temperature, supporting the hypothesis that evolutionary young genes could play roles in adaptation to harsh environments. KLMX_70384 is predicted to encode an 83 aa peptide, and RNA sequencing and ribo-sequencing were used to confirm transcription and translation of the gene. The precise function of KLMX_70384 remains unknown but some features are suggestive of RNA-binding activity. The gene is located in what was previously considered an intergenic region of the genome, which lacks homologues in other yeasts or in databases. Overall, the data support the hypothesis that genes that arose de novo in K. marxianus after the speciation event that separated K. marxianus and K. lactis contribute to some of its unique traits.