Improved detection and phylogenetic analysis of plant proteins containing LysM domains.

Plants perceive N-acetyl-d-glucosamine-containing oligosaccharides that play a role in the interaction with bacteria and fungi, through cell-surface receptors containing a tight bundle of three LysM domains in their extracellular region. However, the identification of LysM domains of receptor-like kinases (RLK)/receptor-like proteins (RLP) using sequence based methods has led to some ambiguity, as some proteins have been annotated with only one or two LysM domains. This missing annotation was likely produced by the failure of the LysM hidden Markov model (HMM) from the Pfam database to correctly identify some LysM domains in proteins of plant origin. In this work, we provide improved HMMs for LysM domain detection in plants, that were built from the structural alignment of manually curated LysM domain structures from the Protein Data Bank and AlphaFold Protein Structure Database. Furthermore, we evaluated different sets of ligand-specific HMMs that were able to correctly classify a limited set of fully characterised RLK/Ps by their ligand specificity. In contrast, the phylogenetic analysis of the extracellular region of RLK/Ps, or of their individual LysM domains, was unable to discriminate these proteins by their ligand specificity. The HMMs reported here will allow a more sensitive detection of plant proteins containing LysM domains and help improve their characterisation.

Identification of Ensifer meliloti genes required for survival during peat-based bioinoculant maturation by STM-seq.

Rhizobial inoculants are sold either as rhizobia within a liquid matrix; or as rhizobia adhered to granules composed of peat prill or finely ground peat moss. During the production of peat-based inoculants, a series of physiological changes occur that result in an increased capability of the rhizobia to survive on the seeds. The number of viable rhizobia on preinoculated seeds at the point of sale, however, is often a limiting factor, as is the inefficiency of the inoculant bacteria to compete with the local rhizobia for the host colonization. In the present work, we used STM-seq for the genome-wide screening of Ensifer meliloti mutants affected in the survival during the maturation of peat-based inoculant formulations. Through this approach, we were able to identify a set of mutants whose behavior suggests that persistence in peat inoculants involves a complex phenotype that is connected to diverse cellular activities, mainly related to satisfying the requirements of bacterial nutrition (e.g., carbon sources, ions) and to coping with specific stresses (e.g., oxidative, mutational). These results also provide a base knowledge that could be used to more completely understand the survival mechanisms used by rhizobia during the maturation of peat-based inoculants, as well as for the design and implementation of practical strategies to improve inoculant formulations.

ubiF is involved in acid stress tolerance and symbiotic competitiveness in Rhizobium favelukesii LPU83.

The acidity of soils significantly reduces the productivity of legumes mainly because of the detrimental effects of hydrogen ions on the legume plants, leading to the establishment of an inefficient symbiosis and poor biological nitrogen fixation. We recently reported the analysis of the fully sequenced genome of Rhizobium favelukesii LPU83, an alfalfa-nodulating rhizobium with a remarkable ability to grow, nodulate and compete in acidic conditions. To gain more insight into the genetic mechanisms leading to acid tolerance in R. favelukesii LPU83, we constructed a transposon mutant library and screened for mutants displaying a more acid-sensitive phenotype than the parental strain. We identified mutant Tn833 carrying a single-transposon insertion within LPU83_2531, an uncharacterized short ORF located immediately upstream from ubiF homolog. This gene encodes a protein with an enzymatic activity involved in the biosynthesis of ubiquinone. As the transposon was inserted near the 3' end of LPU83_2531 and these genes are cotranscribed as a part of the same operon, we hypothesized that the phenotype in Tn833 is most likely due to a polar effect on ubiF transcription.We found that a mutant in ubiF was impaired to grow at low pH and other abiotic stresses including 5 mM ascorbate and 0.500 mM Zn2+. Although the ubiF mutant retained the ability to nodulate alfalfa and Phaseolus vulgaris, it was unable to compete with the R. favelukesii LPU83 wild-type strain for nodulation in Medicago sativa and P. vulgaris, suggesting that ubiF is important for competitiveness. Here, we report for the first time an ubiF homolog being essential for nodulation competitiveness and tolerance to specific stresses in rhizobia.

Analysis of SARS-CoV-2 synonymous codon usage evolution throughout the COVID-19 pandemic.

SARS-CoV-2, the seventh coronavirus known to infect humans, can cause severe life-threatening respiratory pathologies. To better understand SARS-CoV-2 evolution, genome-wide analyses have been made, including the general characterization of its codons usage profile. Here we present a bioinformatic analysis of the evolution of SARS-CoV-2 codon usage over time using complete genomes collected since December 2019. Our results show that SARS-CoV-2 codon usage pattern is antagonistic to, and it is getting farther away from that of the human host. Further, a selection of deoptimized codons over time, which was accompanied by a decrease in both the codon adaptation index and the effective number of codons, was observed. All together, these findings suggest that SARS-CoV-2 could be evolving, at least from the perspective of the synonymous codon usage, to become less pathogenic.

Easy identification of insertion sequence mobilization events in related bacterial strains with ISCompare.

Bacterial genomes are composed of core and accessory genomes. The first is composed of housekeeping and essential genes, while the second is highly enriched in mobile genetic elements, including transposable elements (TEs). Insertion sequences (ISs), the smallest TEs, have an important role in genome evolution, and contribute to bacterial genome plasticity and adaptability. ISs can spread in a genome, presenting different locations in nearly related strains, and producing phenotypic variations. Few tools are available which can identify differentially located ISs (DLISs) on assembled genomes. Here, we introduce ISCompare, a new program to profile IS mobilization events in related bacterial strains using complete or draft genome assemblies. ISCompare was validated using artificial genomes with simulated random IS insertions and real sequences, achieving the same or better results than other available tools, with the advantage that ISCompare can analyze multiple ISs at the same time and outputs a list of candidate DLISs. ISCompare provides an easy and straightforward approach to look for differentially located ISs on bacterial genomes.

Comparative Analysis of Three Bradyrhizobium diazoefficiens Genomes Show Specific Mutations Acquired during Selection for a Higher Motility Phenotype and Adaption to Laboratory Conditions.

Microbial genomes are being extensively studied using next-generation sequencing technologies in order to understand the changes that occur under different selection regimes. In this work, the number and type of mutations that have occurred in three Bradyrhizobium diazoefficiens USDA 110T strains under laboratory conditions and during selection for a more motile phenotypic variant were analyzed. Most of the mutations found in both processes consisted of single nucleotide polymorphisms, single nucleotide deletions or insertions. In the case of adaptation to laboratory conditions, half of the changes occurred within intergenic regions, and around 80% were insertions. When the more motile phenotypic variant was evaluated, eight single nucleotide polymorphisms and an 11-bp deletion were found, although none of them was directly related to known motility or chemotaxis genes. Two mutants were constructed to evaluate the 11-bp deletion affecting the alpha subunit of 2-oxoacid:acceptor oxidoreductase (AAV28_RS30705-blr6743). The results showed that this single deletion was not responsible for the enhanced motility phenotype. IMPORTANCE The genetic and genomic changes that occur under laboratory conditions in Bradyrhizobium diazoefficiens genomes remain poorly studied. Only a few genome sequences of this important nitrogen-fixing species are available, and there are no genome-wide comparative analyses of related strains. In the present work, we sequenced and compared the genomes of strains derived from a parent strain, B. diazoefficiens USDA 110, that has undergone processes of repeated culture in the laboratory environment, or phenotypic selection toward antibiotic resistance and enhanced motility. Our results represent the first analysis in B. diazoefficiens that provides insights into the specific mutations that are acquired during these processes.

The two-component system ActJK is involved in acid stress tolerance and symbiosis in Sinorhizobium meliloti.

The nitrogen-fixing α-proteobacterium Sinorhizobium meliloti genome codifies at least 50 response regulator (RR) proteins mediating different and, in many cases, unknown processes. RR-mutant library screening allowed us to identify genes potentially implicated in survival to acid conditions. actJ mutation resulted in a strain with reduced growth rate under mildly acidic conditions as well as a lower capacity to tolerate a sudden shift to lethal acidic conditions compared with the parental strain. Mutation of the downstream gene actK, which encodes for a histidine kinase, showed a similar phenotype in acidic environments suggesting a functional two-component system. Interestingly, even though nodulation kinetics, quantity, and macroscopic morphology of Medicago sativa nodules were not affected in actJ and actK mutants, ActK was required to express the wild-type nitrogen fixation phenotype and ActJK was necessary for full bacteroid development and nodule occupancy. The actJK regulatory system presented here provides insights into an evolutionary process in rhizobium adaptation to acidic environments and suggests that actJK-controlled functions are crucial for optimal symbiosis development.

Novel tnpR-based transposable promoter traps suitable for RIVET studies in different gram-negative bacteria.

The preparation of plasmid-borne RIVET libraries can be troublesome when high genomic coverages are needed. We present here the construction and functional validation of a new set of miniTn5 promoter traps to generate tnpR-based RIVET libraries. The ability to generate tnpR transcriptional fusions by transposition will significantly facilitate the setup of RIVET studies in those bacteria where Tn5 transposition is operative.

Genetic and functional characterization of a yet-unclassified rhizobial Dtr (DNA-transfer-and-replication) region from a ubiquitous plasmid conjugal system present in Sinorhizobium meliloti, in Sinorhizobium medicae, and in other nonrhizobial Gram-negative bacteria.

Rhizobia are Gram-negative bacteria that live in soils and associate with leguminous plants to establish nitrogen-fixing symbioses. The ability of these bacteria to undergo horizontal gene transfer (HGT) is thought to be one of the main features to explain both the origin of their symbiotic life-style and the plasticity and dynamics of their genomes. In our laboratory we have previously characterized at the species level the non-pSym plasmid mobilome in Sinorhizobium meliloti, the symbiont of Medicago spp., and have found a high incidence of conjugal activity in many plasmids (Pistorio et al., 2008). In this work we characterized the Dtr (DNA-transfer-and-replication) region of one of those plasmids, pSmeLPU88b. This mobilization region was found to represent a previously unclassified Dtr type in rhizobia (hereafter type-IV), highly ubiquitous in S. meliloti and found in other genera of Gram-negative bacteria as well; including Agrobacterium, Ochrobactrum, and Chelativorans. The oriT of the type-IV Dtr described here could be located by function within a DNA fragment of 278 bp, between the divergent genes parA and mobC. The phylogenetic analysis of the cognate relaxase MobZ indicated that this protein groups close to the previously defined MOB(P3) and MOB(P4) type of enzymes, but is located in a separate and novel cluster that we have designated MOB(P0). Noteworthy, MOB(P0) and MOB(P4) relaxases were frequently associated with plasmids present in rhizospheric soil bacteria. A comparison of the nod-gene locations with the phylogenetic topology of the rhizobial relaxases revealed that the symbiotic genes are found on diverse plasmids bearing any of the four Dtr types, thus indicating that pSym plasmids are not specifically associated with any particular mobilization system. Finally, we demonstrated that the type-IV Dtr promoted the mobilization of plasmids from S. meliloti to Sinorhizobium medicae as well as from these rhizobia to other bacteria by means of their own helper functions. The results present an as-yet-unclassified and seemingly ubiquitous conjugal system that provides a mechanistic support for the HGT between sympatric rhizobia of Medicago roots, and between other soil and rhizospheric bacteria.

Conjugal properties of the Sinorhizobium meliloti plasmid mobilome.

The biology and biochemistry of plasmid transfer in soil bacteria is currently under active investigation because of its central role in prokaryote adaptation and evolution. In this work, we examined the conjugal properties of the cryptic plasmids present in a collection of the N(2)-fixing legume-symbiont Sinorhizobium meliloti. The study was performed on 65 S. meliloti isolates recovered from 25 humic soils of Argentina, which were grouped into 22 plasmid-profile types [i.e. plasmid operational taxonomic units (OTUs)]. The cumulative Shannon index calculated for the observed plasmid profiles showed a clear saturation plateau, thus indicating an adequate representation of the S. meliloti plasmid-profile types in the isolates studied. The results show that isolates of nearly 14% of the plasmid OTUs hosted transmissible plasmids and that isolates of 29% of the plasmid OTUs were able to retransfer the previously characterized mobilizable-cryptic plasmid pSmeLPU88b to a third recipient strain. It is noteworthy that isolates belonging to 14% of the plasmid OTUs proved to be refractory to the entrance of the model plasmid pSmeLPU88b, suggesting either the presence of surface exclusion phenomena or the occurrence of restriction incompatibility with the incoming replicon. Incompatibility for replication between resident plasmids and plasmid pSmeLPU88b was observed in c. 20% of the OTUs. The results reported here reveal a widespread compatibility among the conjugal functions of the cryptic plasmids in S. meliloti, and this fact, together with the observed high proportion of existing donor genotypes, points to the extrachromosomal compartment of the species as being an extremely active plasmid mobilome.