Phylogenomic approaches reveal a robust time-scale phylogeny of the Terminal Fusarium Clade.

The Terminal Fusarium Clade (TFC) is a group in the Nectriaceae family with agricultural and clinical relevance. In recent years, various phylogenies have been presented in the literature, showing disagreement in the topologies, but only a few studies have conducted analyses on the divergence time scale of the group. Therefore, the evolutionary history of this group is still being determined. This study aimed to understand the evolutionary history of the TFC from a phylogenomic perspective. To achieve this objective, we performed a phylogenomic analysis using the available genomes in GenBank and ran eight different pipelines. We presented a new robust topology of the TFC that differs at some nodes from previous studies. These new relationships allowed us to formulate new hypotheses about the evolutionary history of the TFC. We also inferred new divergence time estimates, which differ from those of previous studies due to topology discordances and taxon sampling. The results suggested an important diversification process in the Neogene period, likely associated with the diversification and predominance of terrestrial ecosystems by angiosperms. In conclusion, we presented a robust time-scale phylogeny that allowed us to formulate new hypotheses regarding the evolutionary history of the TFC.

Improving the Comprehension of Pathogenicity and Phylogeny in 'Candidatus Phytoplasma meliae' through Genome Characterization.

'Candidatus Phytoplasma meliae' is a pathogen associated with chinaberry yellowing disease, which has become a major phytosanitary problem for chinaberry forestry production in Argentina. Despite its economic impact, no genome information of this phytoplasma has been published, which has hindered its characterization at the genomic level. In this study, we used a metagenomics approach to analyze the draft genome of the 'Ca. P. meliae' strain ChTYXIII. The draft assembly consisted of twenty-one contigs with a total length of 751.949 bp, and annotation revealed 669 CDSs, 34 tRNAs, and 1 set of rRNA operons. The metabolic pathways analysis showed that ChTYXIII contains the complete core genes for glycolysis and a functional Sec system for protein translocation. Our phylogenomic analysis based on 133 single-copy genes and genome-to-genome metrics supports the classification as unique 'Ca. P. species' within the MPV clade. We also identified 31 putative effectors, including a homolog to SAP11 and others that have only been described in this pathogen. Our ortholog analysis revealed 37 PMU core genes in the genome of 'Ca. P. meliae' ChTYXIII, leading to the identification of 2 intact PMUs. Our work provides important genomic information for 'Ca. P. meliae' and others phytoplasmas for the 16SrXIII (MPV) group.

Prediction of DNA-Binding Transcription Factors in Bacteria and Archaea Genomes.

DNA-binding transcription factors (TFs) play a central role in the gene expression of all organisms, from viruses to humans, including bacteria and archaea. The role of these proteins is the fate of gene expression in the context of environmental challenges. Because thousands of genomes have been sequenced to date, predictions of the encoded proteins are validated through the use of bioinformatics tools to obtain the necessary experimental, posterior knowledge. In this chapter, we describe three approaches to identify TFs in protein sequences. The first approach integrates the results of sequence comparisons and PFAM assignments, using as reference a manually curated collection of TFs. The second approach considers the prediction of DNA-binding structures, such as the classical helix-turn-helix (HTH); and the third approach considers a deep learning model. We suggest that all approaches must be considered together to increase the possibility of identifying new TFs in bacterial and archaeal genomes.

Prediction of protein architectures involved in the signaling-pathway initiating sporulation in Firmicutes.

OBJECTIVES: Like many other proteins, those belonging to the signal transduction cascade initiating sporulation (Spo0 pathway) have conserved protein domains (Capra and Laub in Annu Rev Microbiol 66:325-47, 2012). Improvements in bioinformatics applications to discover proteins involved in the initiation of the sporulating cascade in newly sequenced genomes is an important task that requires rigorous comparative genomic methods and manual curation to identify endospore-forming bacteria. This note aims to present a collection of predicted proteins involved in the Spo0 pathway found in the proteomes of fully sequenced and manually curated endospore-forming Firmicutes species. This collection may serve as a guide to conduct future experiments in endospore formers in genomic and metagenomic projects. DATA DESCRIPTION: Similar to the report of Davidson et al. (PLoS Genet 14:1-33, 2018), we used Pfam profiles (El-Gebali et al. in Nucleic Acids Res 47:D427-32, 2019) defining each protein and the genomic context surrounding the query gene to predict probable orthologs of the Spo0 pathway in Firmicutes. We present in this note a collection of 325 Firmicutes species organized by phylogenetic class and classified as spore formers, non-spore formers or unknown spore phenotype based on published literature, for which we predicted probable orthologs defining the signal transduction pathway initiating sporulation.

Gene Tags Assessment by Comparative Genomics (GTACG): A User-Friendly Framework for Bacterial Comparative Genomics.

Genomics research has produced an exponential amount of data. However, the genetic knowledge pertaining to certain phenotypic characteristics is lacking. Also, a considerable part of these genomes have coding sequences (CDSs) with unknown functions, posing additional challenges to researchers. Phylogenetically close microorganisms share much of their CDSs, and certain phenotypes unique to a set of microorganisms may be the result of the genes found exclusively in those microorganisms. This study presents the GTACG framework, an easy-to-use tool for identifying in the subgroups of bacterial genomes whose microorganisms have common phenotypic characteristics, to find data that differentiates them from other associated genomes in a simple and fast way. The GTACG analysis is based on the formation of homologous CDS clusters from local alignments. The front-end is easy to use, and the installation packages have been developed to enable users lacking knowledge of programming languages or bioinformatics analyze high-throughput data using the tool. The validation of the GTACG framework has been carried out based on a case report involving a set of 161 genomes from the Xanthomonadaceae family, in which 19 families of orthologous proteins were found in 90% of the plant-associated genomes, allowing the identification of the proteins potentially associated with adaptation and virulence in plant tissue. The results show the potential use of GTACG in the search for new targets for molecular studies, and GTACG can be used as a research tool by biologists who lack advanced knowledge in the use of computational tools for bacterial comparative genomics.

Study on the inhibitory effect of furafylline and troleandomycin in the 7-methoxyresorufin-O-demethylase and nifedipine oxidase activities in hepatic microsomes from four poultry species using high-performance liquid chromatography coupled with fluorescence and ultraviolet detection.

The present study reports the in vitro studies with furafylline and troleandomycin (TAO) as specific inhibitors of activities 7-methoxyresorufin-O-demethylase (MROD) and nifedipine oxidase, catalyzed by cytochrome P450 1 A2 (CYP1 A2) and 3A4 human enzymes, respectively, in hepatic microsomes of quail, duck, turkey and chicken. The results suggest that in chicken and quail the MROD activity is carried out by orthologs CYP1 A4 and 1 A5, meanwhile in duck and turkey by a CYP1 A5 ortholog. The nifedipine oxidase activity is carried out by orthologs of the CYP3A family in the four bird species. The use of furafylline and TAO significantly decreased these activities (P < 0.05) and suggested that the biotransformation of resorufin methyl ether (RME) may be related to more than one avian ortholog.

Gene Phylogenies and Orthologous Groups.

This chapter covers the theory and practice of ortholog gene set computation. In the theoretical part we give detailed and formal descriptions of the relevant concepts. We also cover the topic of graph-based clustering as a tool to compute ortholog gene sets. In the second part we provide an overview of practical considerations intended for researchers who need to determine orthologous genes from a collection of annotated genomes, briefly describing some of the most popular programs and resources currently available for this task.

The Transcriptional Regulators of the CRP Family Regulate Different Essential Bacterial Functions and Can Be Inherited Vertically and Horizontally.

One of the best-studied transcriptional regulatory proteins in bacteria is the Escherichia coli catabolite repressor protein (CRP) that when complexed with 3'-5'-cyclic AMP (cAMP) changes its conformation and interacts with specific DNA-sequences. CRP DNA-binding can result in positive or negative regulation of gene expression depending on the position of its interaction with respect to RNA polymerase binding site. The aim of this work is to review the biological role and phylogenetic relations that some members of the CRP family of transcriptional regulators (also known as cAMP receptor protein family) have in different bacterial species. This work is not intended to give an exhaustive revision of bacterial CRP-orthologs, but to provide examples of the role that these proteins play in the expression of genes that are fundamental for the life style of some bacterial species. We highlight the conservation of their structural characteristics and of their binding to conserved-DNA sequences, in contrast to their very diverse repertoire of gene activation. CRP activates a wide variety of fundamental genes for the biological characteristic of each bacterial species, which in several instances form part of their core-genome (defined as the gene sequences present in all members of a bacterial species). We present evidence that support the fact that some of the transcriptional regulators that belong to the CRP family in different bacterial species, and some of the genes that are regulated by them, can be inherited by horizontal gene transfer. These data are discussed in the framework of bacterial evolution models.

Multiple isoforms for the catalytic subunit of PKA in the basal fungal lineage Mucor circinelloides.

Protein kinase A (PKA) activity is involved in dimorphism of the basal fungal lineage Mucor. From the recently sequenced genome of Mucor circinelloides we could predict ten catalytic subunits of PKA. From sequence alignment and structural prediction we conclude that the catalytic core of the isoforms is conserved, and the difference between them resides in their amino termini. This high number of isoforms is maintained in the subdivision Mucoromycotina. Each paralogue, when compared to the ones form other fungi is more homologous to one of its orthologs than to its paralogs. All of these fungal isoforms cannot be included in the class I or II in which fungal protein kinases have been classified. mRNA levels for each isoform were measured during aerobic and anaerobic growth. The expression of each isoform is differential and associated to a particular growth stage. We reanalyzed the sequence of PKAC (GI 20218944), the only cloned sequence available until now for a catalytic subunit of M. circinelloides. PKAC cannot be classified as a PKA because of its difference in the conserved C-tail; it shares with PKB a conserved C2 domain in the N-terminus. No catalytic activity could be measured for this protein nor predicted bioinformatically. It can thus be classified as a pseudokinase. Its importance can not be underestimated since it is expressed at the mRNA level in different stages of growth, and its deletion is lethal.

A bioinformatic survey of distribution, conservation, and probable functions of LuxR solo regulators in bacteria.

LuxR solo transcriptional regulators contain both an autoinducer binding domain (ABD; N-terminal) and a DNA binding Helix-Turn-Helix domain (HTH; C-terminal), but are not associated with a cognate N-acyl homoserine lactone (AHL) synthase coding gene in the same genome. Although a few LuxR solos have been characterized, their distributions as well as their role in bacterial signal perception and other processes are poorly understood. In this study we have carried out a systematic survey of distribution of all ABD containing LuxR transcriptional regulators (QS domain LuxRs) available in the InterPro database (IPR005143), and identified those lacking a cognate AHL synthase. These LuxR solos were then analyzed regarding their taxonomical distribution, predicted functions of neighboring genes and the presence of complete AHL-QS systems in the genomes that carry them. Our analyses reveal the presence of one or multiple predicted LuxR solos in many proteobacterial genomes carrying QS domain LuxRs, some of them harboring genes for one or more AHL-QS circuits. The presence of LuxR solos in bacteria occupying diverse environments suggests potential ecological functions for these proteins beyond AHL and interkingdom signaling. Based on gene context and the conservation levels of invariant amino acids of ABD, we have classified LuxR solos into functionally meaningful groups or putative orthologs. Surprisingly, putative LuxR solos were also found in a few non-proteobacterial genomes which are not known to carry AHL-QS systems. Multiple predicted LuxR solos in the same genome appeared to have different levels of conservation of invariant amino acid residues of ABD questioning their binding to AHLs. In summary, this study provides a detailed overview of distribution of LuxR solos and their probable roles in bacteria with genome sequence information.

Non conserved residues between Cqm1 and Aam1 mosquito α-glucosidases are critical for the capacity of Cqm1 to bind the Binary toxin from Lysinibacillus sphaericus.

The Binary (Bin) toxin from the entomopathogenic bacterium Lysinibacillus sphaericus acts on larvae of the culicid Culex quinquefasciatus through its binding to Cqm1, a midgut-bound α-glucosidase. Specific binding by the BinB subunit to the Cqm1 receptor is essential for toxicity however the toxin is unable to bind to the Cqm1 ortholog from the refractory species Aedes aegypti (Aam1). Here, to investigate the molecular basis for the interaction between Cqm1 and BinB, recombinant Cqm1 and Aam1 were first expressed as soluble forms in Sf9 cells. The two proteins were found to display the same glycosilation patterns and BinB binding properties as the native α-glucosidases. Chimeric constructs were then generated through the exchange of reciprocal fragments between the corresponding cqm1 and aam1 cDNAs. Subsequent expression and binding experiments defined a Cqm1 segment encompassing residues S129 and A312 as critical for the interaction with BinB. Through site directed mutagenesis experiments, replacing specific sets of residues from Cqm1 with those of Aam1, the 159GG160 doublet was required for this interaction. Molecular modeling mapped these residues to an exposed loop within the Cqm1's structure, compatible with a target site for BinB and providing a possible explanation for its lack of binding to Aam1.

Chromosomal localization of selected SCARs converted from new RAPD, ISSR and R-ISSR markers linked to rye (Secale cereale L. ) tolerance to nutrient deprivation stress identified using bulk segregant analysis

Background: An adaptive mechanism in plant roots is initiated in the event of nitrogen and potassium deficiency, and it facilitates the active uptake of these elements in order to ensure plant growth and survival in stress conditions. Signaling and transduction of signals in response to changing nitrogen and potassium concentrations is a complex process, affected by interactions between various gene expression products, and often subjected to modifications. Results: In order to identify genotypic differences between phenotypes of two populations of recombinant inbred rye lines (153/79-1 x Ot1-3 and Ot0-6 x Ot1-3) in response to nutrition stress caused by nitrogen and potassium deficiency at the seedling stage, bulk segregant analysis was utilized. Identification of genotypic differences between and within pooled DNA samples involved 424 RAPD, 120 ISSR primers and 50 combinations of R-ISSR. Identified markers were sequenced and converted to SCAR, attributing to them unique ESTs annotations, and chromosomal ones to selected localizations. Significant relationships with the examined trait were described for nine and eight RAPD markers, four and five ISSR, one and three R-ISSR markers for population 153/79-1 x Ot1-3 and Ot0-6 x Ot1-3, respectively. Sequences identified for the rye genome were characterized by a uniqueness and a similarity to the sequence of aquaporin PIP1, a gene encoding protein related to the function of the transcription factor in plant response to iron deficiency and the putative ethylene-responsive transcription factor, cytosolic acetyl-CoA carboxylase, HvHKT1 transporter, as well as HCBT proteins. Conclusion: Identified molecular markers differentiating rye genotypes of extreme response of root system on nitrogen and potassium deficiency play a significant role in systemic plant response to stress, including stress caused by nitrogen and potassium deficiency. They may constitute a system facilitating selection, and together with the material they are described in, they may be a starting point for research on mechanisms of sensing and transduction of signal across the plant.