Codon bias, tRNA pools and horizontal gene transfer

Intraspecific and interspecific variations in the synonymous codon usage in mitochondrial genomes of 8 pleurotus strains

In this study, we investigated the codon bias of twelve mitochondrial core protein coding genes (PCGs) in eight Pleurotus strains, two of which are from the same species. The results revealed that the codons of all Pleurotus strains had a preference for ending in A/T. Furthermore, the correlation between codon base compositions and codon adaptation index (CAI), codon bias index (CBI) and frequency of optimal codons (FOP) indices was also detected, implying the influence of base composition on codon bias. The two P. ostreatus species were found to have differences in various base bias indicators. The average effective number of codons (ENC) of mitochondrial core PCGs of Pleurotus was found to be less than 35, indicating strong codon preference of mitochondrial core PCGs of Pleurotus. The neutrality plot analysis and PR2-Bias plot analysis further suggested that natural selection plays an important role in Pleurotus codon bias. Additionally, six to ten optimal codons (ΔRSCU > 0.08 and RSCU > 1) were identified in eight Pleurotus strains, with UGU and ACU being the most widely used optimal codons in Pleurotus. Finally, based on the combined mitochondrial sequence and RSCU value, the genetic relationship between different Pleurotus strains was deduced, showing large variations between them. This research has improved our understanding of synonymous codon usage characteristics and evolution of this important fungal group.

Transcriptional landscape of the pMP7017 megaplasmid and its impact on the Bifidobacterium breve UCC2003 transcriptome

The 190 kb megaplasmid pMP7017 of Bifidobacterium breve JCM7017 represents the first conjugative and largest plasmid characterised within this genus to date. In the current study, we adopted an integrated approach combining transcriptomics, whole genome comparative analysis and metagenomic data mining to understand the biology of pMP7017 and related megaplasmids, and to assess the impact of plasmid-carriage on the host strain. The data generated revealed variations within basic features of promoter elements which correlate with a high level of transcription on the plasmid and highlight the transcriptional activity of genes encoding both offensive and defensive adaptations, including a Type IIL restriction-modification system, an anti-restriction system and four Type II toxin-antitoxin systems. Furthermore, a highly transcribed tmRNA, which likely provides translational support to the host strain, was identified, making pMP7017 the first plasmid of the Bifidobacterium genus and the smallest plasmid known to express a tmRNA. Analyses of synteny and variability among pMP7017 and related plasmids indicate substantial diversity in gene organisation and accessory gene cargo highlighting diverse (co-)evolution and potential host-specific rearrangements and adaptations. Systematic analysis of the codon usage profile of transcriptionally active pMP7017-encoded genes suggests that pMP7017 originated from (sub)species of Bifidobacterium longum. Furthermore, mining of metagenomic data suggests the presence of pMP7017-homologues in ~10% of microbiome samples, mostly infants and/or mothers from various geographical locations. Comparative transcriptome analysis of the B. breve UCC2003 chromosome in the presence or absence of pMP7017 revealed differential expression of genes representing 8% of the total gene pool. Genes involved in genetic information processing were exclusively upregulated, while altered expression of genes involved in biofilm production and polysaccharide biosynthesis was also observed.

Synonymous and Nonsynonymous Substitutions in Dictyostelium discoideum Ammonium Transporter amtA Are Necessary for Functional Complementation in Saccharomyces cerevisiae

Ammonium transporters are present in all three domains of life. They have undergone extensive horizontal gene transfer (HGT), gene duplication, and functional diversification and therefore offer an excellent paradigm to study protein evolution. We attempted to complement a mep1Δmep2Δmep3Δ strain of Saccharomyces cerevisiae (triple-deletion strain), which otherwise cannot grow on ammonium as a sole nitrogen source at concentrations of <3 mM, with amtA of Dictyostelium discoideum, an orthologue of S. cerevisiae MEP2. We observed that amtA did not complement the triple-deletion strain of S. cerevisiae for growth on low-ammonium medium. We isolated two mutant derivatives of amtA (amtA M1 and amtA M2) from a PCR-generated mutant plasmid library that complemented the triple-deletion strain of S. cerevisiae. amtA M1 bears three nonsynonymous and two synonymous substitutions, which are necessary for its functionality. amtA M2 bears two nonsynonymous substitutions and one synonymous substitution, all of which are necessary for functionality. Interestingly, AmtA M1 transports ammonium but does not confer methylamine toxicity, while AmtA M2 transports ammonium and confers methylamine toxicity, demonstrating functional diversification. Preliminary biochemical analyses indicated that the mutants differ in their conformations as well as their mechanisms of ammonium transport. These intriguing results clearly point out that protein evolution cannot be fathomed by studying nonsynonymous and synonymous substitutions in isolation. The above-described observations have significant implications for various facets of biological processes and are discussed in detail. IMPORTANCE Functional diversification following gene duplication is one of the major driving forces of protein evolution. While the role of nonsynonymous substitutions in the functional diversification of proteins is well recognized, knowledge of the role of synonymous substitutions in protein evolution is in its infancy. Using functional complementation, we isolated two functional alleles of the D. discoideum ammonium transporter gene (amtA), which otherwise does not function in S. cerevisiae as an ammonium transporters. One of them is an ammonium transporter, while the other is an ammonium transporter that also confers methylammonium (ammonium analogue) toxicity, suggesting functional diversification. Surprisingly, both alleles require a combination of synonymous and nonsynonymous substitutions for their functionality. These results bring out a hitherto-unknown pathway of protein evolution and pave the way for not only understanding protein evolution but also interpreting single nucleotide polymorphisms (SNPs).

Bipartite Genomes in Enterobacterales: Independent Origins of Chromids, Elevated Openness and Donors of Horizontally Transferred Genes

Multipartite bacteria have one chromosome and one or more chromid. Chromids are believed to have properties that enhance genomic flexibility, making them a favored integration site for new genes. However, the mechanism by which chromosomes and chromids jointly contribute to this flexibility is not clear. To shed light on this, we analyzed the openness of chromosomes and chromids of the two bacteria, Vibrio and Pseudoalteromonas, both which belong to the Enterobacterales order of Gammaproteobacteria, and compared the genomic openness with that of monopartite genomes in the same order. We applied pangenome analysis, codon usage analysis and the HGTector software to detect horizontally transferred genes. Our findings suggest that the chromids of Vibrio and Pseudoalteromonas originated from two separate plasmid acquisition events. Bipartite genomes were found to be more open compared to monopartite. We found that the shell and cloud pangene categories drive the openness of bipartite genomes in Vibrio and Pseudoalteromonas. Based on this and our two recent studies, we propose a hypothesis that explains how chromids and the chromosome terminus region contribute to the genomic plasticity of bipartite genomes.

Machine learning classifiers predict key genomic and evolutionary traits across the kingdoms of life

In this study, we investigate how an organism's codon usage bias can serve as a predictor and classifier of various genomic and evolutionary traits across the domains of life. We perform secondary analysis of existing genetic datasets to build several AI/machine learning models. When trained on codon usage patterns of nearly 13,000 organisms, our models accurately predict the organelle of origin and taxonomic identity of nucleotide samples. We extend our analysis to identify the most influential codons for phylogenetic prediction with a custom feature ranking ensemble. Our results suggest that the genetic code can be utilized to train accurate classifiers of taxonomic and phylogenetic features. We then apply this classification framework to open reading frame (ORF) detection. Our statistical model assesses all possible ORFs in a nucleotide sample and rejects or deems them plausible based on the codon usage distribution. Our dataset and analyses are made publicly available on GitHub and the UCI ML Repository to facilitate open-source reproducibility and community engagement.

Tracking of Diversity and Evolution in the Brown Rot Fungi Monilinia fructicola, Monilinia fructigena, and Monilinia laxa

Monilinia species are among the most devastating fungi worldwide as they cause brown rot and blossom blight on fruit trees. To understand the molecular bases of their pathogenic lifestyles, we compared the newly assembled genomes of single strains of Monilinia fructicola, M. fructigena and M. laxa, with those of Botrytis cinerea and Sclerotinia sclerotiorum, as the closest species within Sclerotiniaceae. Phylogenomic analysis of orthologous proteins and syntenic investigation suggest that M. laxa is closer to M. fructigena than M. fructicola, and is closest to the other investigated Sclerotiniaceae species. This indicates that M. laxa was the earliest result of the speciation process. Distinct evolutionary profiles were observed for transposable elements (TEs). M. fructicola and M. laxa showed older bursts of TE insertions, which were affected (mainly in M. fructicola) by repeat-induced point (RIP) mutation gene silencing mechanisms. These suggested frequent occurrence of the sexual process in M. fructicola. More recent TE expansion linked with low RIP action was observed in M. fructigena, with very little in S. sclerotiorum and B. cinerea. The detection of active non-syntenic TEs is indicative of horizontal gene transfer and has resulted in alterations in specific gene functions. Analysis of candidate effectors, biosynthetic gene clusters for secondary metabolites and carbohydrate-active enzymes, indicated that Monilinia genus has multiple virulence mechanisms to infect host plants, including toxins, cell-death elicitor, putative virulence factors and cell-wall-degrading enzymes. Some species-specific pathogenic factors might explain differences in terms of host plant and organ preferences between M. fructigena and the other two Monilinia species.

Translational demand is not a major source of plasmid-associated fitness costs

Plasmids are key drivers of bacterial evolution because they are crucial agents for the horizontal transfer of adaptive traits, such as antibiotic resistance. Most plasmids entail a metabolic burden that reduces the fitness of their host if there is no selection for plasmid-encoded genes. It has been hypothesized that the translational demand imposed by plasmid-encoded genes is a major mechanism driving the fitness cost of plasmids. Plasmid-encoded genes typically present a different codon usage from host chromosomal genes. As a consequence, the translation of plasmid-encoded genes might sequestrate ribosomes on plasmid transcripts, overwhelming the translation machinery of the cell. However, the pervasiveness and origins of the translation-derived costs of plasmids are yet to be assessed. Here, we systematically altered translation efficiency in the host cell to disentangle the fitness effects produced by six natural antibiotic resistance plasmids. We show that limiting translation efficiency either by reducing the number of available ribosomes or their processivity does not increase plasmid costs. Overall, our results suggest that ribosomal paucity is not a major contributor to plasmid fitness costs. This article is part of the theme issue 'The secret lives of microbial mobile genetic elements'.

Gene Validation and Remodelling Using Proteogenomics of Phytophthora cinnamomi, the Causal Agent of Dieback

Phytophthora cinnamomi is a pathogenic oomycete that causes plant dieback disease across a range of natural ecosystems and in many agriculturally important crops on a global scale. An annotated draft genome sequence is publicly available (JGI Mycocosm) and suggests 26,131 gene models. In this study, soluble mycelial, extracellular (secretome), and zoospore proteins of P. cinnamomi were exploited to refine the genome by correcting gene annotations and discovering novel genes. By implementing the diverse set of sub-proteomes into a generated proteogenomics pipeline, we were able to improve the P. cinnamomi genome annotation. Liquid chromatography mass spectrometry was used to obtain high confidence peptides with spectral matching to both the annotated genome and a generated 6-frame translation. Two thousand seven hundred sixty-four annotations from the draft genome were confirmed by spectral matching. Using a proteogenomic pipeline, mass spectra were used to edit the P. cinnamomi genome and allowed identification of 23 new gene models and 60 edited gene features using high confidence peptides obtained by mass spectrometry, suggesting a rate of incorrect annotations of 3% of the detectable proteome. The novel features were further validated by total peptide support, alongside functional analysis including the use of Gene Ontology and functional domain identification. We demonstrated the use of spectral data in combination with our proteogenomics pipeline can be used to improve the genome annotation of important plant diseases and identify missed genes. This study presents the first use of spectral data to edit and manually annotate an oomycete pathogen.

Evolutionary responses to codon usage of horizontally transferred genes in Pseudomonas aeruginosa: gene retention, amelioration and compensatory evolution

Prokaryote genome evolution is characterized by the frequent gain of genes through horizontal gene transfer (HGT). For a gene, being horizontally transferred can represent a strong change in its genomic and physiological context. If the codon usage of a transferred gene deviates from that of the receiving organism, the fitness benefits it provides can be reduced due to a mismatch with the expression machinery. Consequently, transferred genes with a deviating codon usage can be selected against or elicit evolutionary responses that enhance their integration, such as gene amelioration and compensatory evolution. Within bacterial species, the extent and relative importance of these different mechanisms has never been considered altogether. In this study, a phylogeny-based method was used to investigate the occurrence of these different evolutionary responses in Pseudomonas aeruginosa. Selection on codon usage of genes acquired through HGT was observed over evolutionary time, with the overall codon usage converging towards that of the core genome. Gene amelioration, through the accumulation of synonymous mutations after HGT, did not seem to systematically affect transferred genes. This pattern therefore seemed to be mainly driven by selective retention of transferred genes with an initial codon usage similar to that of the core genes. Additionally, variation in the copy number of tRNA genes was often associated with the acquisition of genes for which the observed variation could enhance their expression. This provides evidence that compensatory evolution might be an important mechanism for the integration of horizontally transferred genes.

Comparative phylo-pangenomics reveals generalist lifestyles in representative Acinetobacter species and proposes candidate gene markers for species identification

Acinetobacter species have the potential to invade and colonize immunocompromised patients, therefore being well-known as opportunistic pathogens. Among these bacteria, the species of the Acinetobacter calcoaceticus-Acinetobacter baumannii "complex" (Acb members) emerge as the main often isolated bacteria in clinical specimens. The unequivocal taxonomy is crucial to correctly identify these species and associated with comparative genomic analyses aids to understand their life-styles as well. In this study, all publicly available Acinetobacter species at the date of this study preparation were analyzed. The results revealed that the Acb members are in fact a complex when phenotypic methods are confronted, while for comparative and phylogenomics analyses this term is misleading, since they composed a monophyletic group instead. Nine best gene markers (response regulator, recJ, recG, phosphomannomutase, pepSY, monovalent cation/H + antiporter subunit D, mnmE, glnE, and bamA) were selected for identification of Acinetobacter species. Moreover, representative strains of each species were split according their isolation sources in the categories: environmental, human, insect and non-human vertebrate. Neither niche-specific genome signature nor niche-associated functional and pathogenic potential were associated with their isolation source, meaning it is not the main force acting on Acinetobacter adaptation in a given niche and corroborating that their ubiquitous distribution is a reflex of their generalist life-styles.

Phylogenomic Analyses of Non-Dikarya Fungi Supports Horizontal Gene Transfer Driving Diversification of Secondary Metabolism in the Amphibian Gastrointestinal Symbiont, Basidiobolus

Research into secondary metabolism (SM) production by fungi has resulted in the discovery of diverse, biologically active compounds with significant medicinal applications. The fungi rich in SM production are taxonomically concentrated in the subkingdom Dikarya, which comprises the phyla Ascomycota and Basidiomycota. Here, we explore the potential for SM production in Mucoromycota and Zoopagomycota, two phyla of nonflagellated fungi that are not members of Dikarya, by predicting and identifying core genes and gene clusters involved in SM. The majority of non-Dikarya have few genes and gene clusters involved in SM production except for the amphibian gut symbionts in the genus Basidiobolus. Basidiobolus genomes exhibit an enrichment of SM genes involved in siderophore, surfactin-like, and terpene cyclase production, all these with evidence of constitutive gene expression. Gene expression and chemical assays also confirm that Basidiobolus has significant siderophore activity. The expansion of SMs in Basidiobolus are partially due to horizontal gene transfer from bacteria, likely as a consequence of its ecology as an amphibian gut endosymbiont.

Universal evolutionary selection for high dimensional silent patterns of information hidden in the redundancy of viral genetic code

Motivation

Understanding how viruses co-evolve with their hosts and adapt various genomic level strategies in order to ensure their fitness may have essential implications in unveiling the secrets of viral evolution, and in developing new vaccines and therapeutic approaches. Here, based on a novel genomic analysis of 2625 different viruses and 439 corresponding host organisms, we provide evidence of universal evolutionary selection for high dimensional ‘silent’ patterns of information hidden in the redundancy of viral genetic code.

Results

Our model suggests that long substrings of nucleotides in the coding regions of viruses from all classes, often also repeat in the corresponding viral hosts from all domains of life. Selection for these substrings cannot be explained only by such phenomena as codon usage bias, horizontal gene transfer and the encoded proteins. Genes encoding structural proteins responsible for building the core of the viral particles were found to include more host-repeating substrings, and these substrings tend to appear in the middle parts of the viral coding regions. In addition, in human viruses these substrings tend to be enriched with motives related to transcription factors and RNA binding proteins. The host-repeating substrings are possibly related to the evolutionary pressure on the viruses to effectively interact with host's intracellular factors and to efficiently escape from the host's immune system.

Supplementary information

Supplementary data are available at Bioinformatics online.

Complete Assembly of the Genome of an Acidovorax citrulli Strain Reveals a Naturally Occurring Plasmid in This Species

Acidovorax citrulli is the causal agent of bacterial fruit blotch (BFB), a serious threat to cucurbit crop production worldwide. Based on genetic and phenotypic properties, A. citrulli strains are divided into two major groups: group I strains have been generally isolated from melon and other non-watermelon cucurbits, while group II strains are closely associated with watermelon. In a previous study, we reported the genome of the group I model strain, M6. At that time, the M6 genome was sequenced by MiSeq Illumina technology, with reads assembled into 139 contigs. Here, we report the assembly of the M6 genome following sequencing with PacBio technology. This approach not only allowed full assembly of the M6 genome, but it also revealed the occurrence of a ∼53 kb plasmid. The M6 plasmid, named pACM6, was further confirmed by plasmid extraction, Southern-blot analysis of restricted fragments and obtention of M6-derivative cured strains. pACM6 occurs at low copy numbers (average of ∼4.1 ± 1.3 chromosome equivalents) in A. citrulli M6 and contains 63 open reading frames (ORFs), most of which (55.6%) encoding hypothetical proteins. The plasmid contains several genes encoding type IV secretion components, and typical plasmid-borne genes involved in plasmid maintenance, replication and transfer. The plasmid also carries an operon encoding homologs of a Fic-VbhA toxin-antitoxin (TA) module. Transcriptome data from A. citrulli M6 revealed that, under the tested conditions, the genes encoding the components of this TA system are among the highest expressed genes in pACM6. Whether this TA module plays a role in pACM6 maintenance is still to be determined. Leaf infiltration and seed transmission assays revealed that, under tested conditions, the loss of pACM6 did not affect the virulence of A. citrulli M6. We also show that pACM6 or similar plasmids are present in several group I strains, but absent in all tested group II strains of A. citrulli.

Natural tuning of restriction endonuclease synthesis by cluster of rare arginine codons

Restriction–modification (R-M) systems are highly widespread among bacteria and archaea, and they appear to play a pivotal role in modulating horizontal gene transfer, as well as in protecting the host organism against viruses and other invasive DNA particles. Type II R-M systems specify two independent enzymes: a restriction endonuclease (REase) and protective DNA methyltransferase (MTase). If the cell is to survive, the counteracting activities as toxin and antitoxin, must be finely balanced in vivo. The molecular basis of this regulatory process remains unclear and current searches for regulatory elements in R-M modules are focused mainly at the transcription step. In this report, we show new aspects of REase control that are linked to translation. We used the EcoVIII R-M system as a model. Both, the REase and MTase genes for this R-M system contain an unusually high number of rare arginine codons (AGA and AGG) when compared to the rest of the E. coli K-12 genome. Clusters of these codons near the N-terminus of the REase greatly affect the translational efficiency. Changing these to higher frequency codons for E. coli (CGC) improves the REase synthesis, making the R-M system more potent to defend its host against bacteriophages. However, this improved efficiency in synthesis reduces host fitness due to increased autorestriction. We hypothesize that expression of the endonuclease gene can be modulated depending on the host genetic context and we propose a novel post-transcriptional mode of R–M system regulation that alleviates the potential lethal action of the restriction enzyme.

Genome-Wide Analysis Reveals Ancestral Lack of Seventeen Different tRNAs and Clade-Specific Loss of tRNA-CNNs in Archaea

Transfer RNA (tRNA) is a category of RNAs that specifically decode messenger RNAs (mRNAs) into proteins by recognizing a set of 61 codons commonly adopted by different life domains. The composition and abundance of tRNAs play critical roles in shaping codon usage and pairing bias, which subsequently modulate mRNA translation efficiency and accuracy. Over the past few decades, effort has been concentrated on evaluating the specificity and redundancy of different tRNA families. However, the mechanism and processes underlying tRNA evolution have only rarely been investigated. In this study, by surveying tRNA genes in 167 completely sequenced genomes, we systematically investigated the composition and evolution of tRNAs in Archaea from a phylogenetic perspective. Our data revealed that archaeal genomes are compact in both tRNA types and copy number. Generally, no more than 44 different types of tRNA are present in archaeal genomes to decode the 61 canonical codons, and most of them have only one gene copy per genome. Among them, tRNA-Met was significantly overrepresented, with an average of three copies per genome. In contrast, the tRNA-UAU and 16 tRNAs with A-starting anticodons (tRNA-ANNs) were rarely detected in all archaeal genomes. The conspicuous absence of these tRNAs across the archaeal phylogeny suggests they might have not been evolved in the common ancestor of Archaea, rather than have lost independently from different clades. Furthermore, widespread absence of tRNA-CNNs in the Methanococcales and Methanobacteriales genomes indicates convergent loss of these tRNAs in the two clades. This clade-specific tRNA loss may be attributing to the reductive evolution of their genomes. Our data suggest that the current tRNA profiles in Archaea are contributed not only by the ancestral tRNA composition, but also by differential maintenance and loss of redundant tRNAs.

GC3-biased gene domains in mammalian genomes

Motivation: Synonymous codon usage bias has been shown to be correlated with many genomic features among different organisms. However, the biological significance of codon bias with respect to gene function and genome organization remains unclear.

Results: Guanine and cytosine content at the third codon position (GC3) could be used as a good indicator of codon bias. Here, we used relative GC3 bias values to compare the strength of GC3 bias of genes in human and mouse. We reported, for the first time, that GC3-rich and GC3-poor gene products might have distinct sub-cellular spatial distributions. Moreover, we extended the view of genomic gene domains and identified conserved GC3 biased gene domains along chromosomes. Our results indicated that similar GC3 biased genes might be co-translated in specific spatial regions to share local translational machineries, and that GC3 could be involved in the organization of genome architecture.

Availability and implementation: Source code is available upon request from the authors.

Contact:zhaozh@nic.bmi.ac.cn or zany1983@gmail.com

Supplementary information:Supplementary data are available at Bioinformatics online.

Optimization of translation profiles enhances protein expression and solubility

mRNA is translated with a non-uniform speed that actively coordinates co-translational folding of protein domains. Using structure-based homology we identified the structural domains in epoxide hydrolases (EHs) and introduced slow-translating codons to delineate the translation of single domains. These changes in translation speed dramatically improved the solubility of two EHs of metagenomic origin in Escherichia coli. Conversely, the importance of transient attenuation for the folding, and consequently solubility, of EH was evidenced with a member of the EH family from Agrobacterium radiobacter, which partitions in the soluble fraction when expressed in E. coli. Synonymous substitutions of codons shaping the slow-transiting regions to fast-translating codons render this protein insoluble. Furthermore, we show that low protein yield can be enhanced by decreasing the free folding energy of the initial 5’-coding region, which can disrupt mRNA secondary structure and enhance ribosomal loading. This study provides direct experimental evidence that mRNA is not a mere messenger for translation of codons into amino acids but bears an additional layer of information for folding, solubility and expression level of the encoded protein. Furthermore, it provides a general frame on how to modulate and fine-tune gene expression of a target protein.

Antibiotic-free selection in biotherapeutics: now and forever

The continuously improving sophistication of molecular engineering techniques gives access to novel classes of bio-therapeutics and new challenges for their production in full respect of the strengthening regulations. Among these biologic agents are DNA based vaccines or gene therapy products and to a lesser extent genetically engineered live vaccines or delivery vehicles. The use of antibiotic-based selection, frequently associated with genetic manipulation of microorganism is currently undergoing a profound metamorphosis with the implementation and diversification of alternative selection means. This short review will present examples of alternatives to antibiotic selection and their context of application to highlight their ineluctable invasion of the bio-therapeutic world.

Genetic determinants of reutericyclin biosynthesis in Lactobacillus reuteri

Reutericyclin is a unique antimicrobial tetramic acid produced by some strains of Lactobacillus reuteri. This study aimed to identify the genetic determinants of reutericyclin biosynthesis. Comparisons of the genomes of reutericyclin-producing L. reuteri strains with those of non-reutericyclin-producing strains identified a genomic island of 14 open reading frames (ORFs) including genes coding for a nonribosomal peptide synthetase (NRPS), a polyketide synthase (PKS), homologues of PhlA, PhlB, and PhlC, and putative transport and regulatory proteins. The protein encoded by rtcN is composed of a condensation domain, an adenylation domain likely specific for d-leucine, and a thiolation domain. rtcK codes for a PKS that is composed of a ketosynthase domain, an acyl-carrier protein domain, and a thioesterase domain. The products of rtcA, rtcB, and rtcC are homologous to the diacetylphloroglucinol-biosynthetic proteins PhlABC and may acetylate the tetramic acid moiety produced by RtcN and RtcK, forming reutericyclin. Deletion of rtcN or rtcABC in L. reuteri TMW1.656 abrogated reutericyclin production but did not affect resistance to reutericyclin. Genes coding for transport and regulatory proteins could be deleted only in the reutericyclin-negative L. reuteri strain TMW1.656ΔrtcN, and these deletions eliminated reutericyclin resistance. The genomic analyses suggest that the reutericyclin genomic island was horizontally acquired from an unknown source during a unique event. The combination of PhlABC homologues with both an NRPS and a PKS has also been identified in the lactic acid bacteria Streptococcus mutans and Lactobacillus plantarum, suggesting that the genes in these organisms and those in L. reuteri share an evolutionary origin.

Environmental shaping of codon usage and functional adaptation across microbial communities

Microbial communities represent the largest portion of the Earth's biomass. Metagenomics projects use high-throughput sequencing to survey these communities and shed light on genetic capabilities that enable microbes to inhabit every corner of the biosphere. Metagenome studies are generally based on (i) classifying and ranking functions of identified genes; and (ii) estimating the phyletic distribution of constituent microbial species. To understand microbial communities at the systems level, it is necessary to extend these studies beyond the species' boundaries and capture higher levels of metabolic complexity. We evaluated 11 metagenome samples and demonstrated that microbes inhabiting the same ecological niche share common preferences for synonymous codons, regardless of their phylogeny. By exploring concepts of translational optimization through codon usage adaptation, we demonstrated that community-wide bias in codon usage can be used as a prediction tool for lifestyle-specific genes across the entire microbial community, effectively considering microbial communities as meta-genomes. These findings set up a 'functional metagenomics' platform for the identification of genes relevant for adaptations of entire microbial communities to environments. Our results provide valuable arguments in defining the concept of microbial species through the context of their interactions within the community.

Analysis of the use of codon pairs in the HE gene of the ISA virus shows a correlation between bias in HPR codon-pair use and mortality rates caused by the virus

BACKGROUND

Segment 6 of the ISA virus codes for hemoagglutinin-esterase (HE). This segment is highly variable, with more than 26 variants identified. The major variation is observed in what is called the high polymorphism region (HPR). The role of the different HPR zones in the viral cycle or evolution remains unknown. However viruses that present the HPR0 are avirulent, while viruses with important deletions in this region have been responsible for outbreaks with high mortality rates. In this work, using bioinformatic tools, we examined the influence of different HPRs on the adaptation of HE genes to the host translational machinery and the relationship to observed virulence.

METHODS

Translational efficiency of HE genes and their HPR were estimated analyzing codon-pair bias (CPB), adaptation to host codon use (codon adaptation index-CAI) and the adaptation to available tRNAs (tAI). These values were correlated with reported mortality for the respective ISA virus and the ΔG of RNA folding. tRNA abundance was inferred from tRNA gene numbers identified in the Salmo salar genome using tRNAScan-SE. Statistical correlation between data was performed using a non-parametric test.

RESULTS

We found that HPR0 contains zones with codon pairs of low frequency and low availability of tRNA with respect to salmon codon-pair usage, suggesting that HPR modifies HE translational efficiency. Although calculating tAI was impossible because one third of tRNAs (~60.000) were tRNA-ala, translational efficiency measured by CPB shows that as HPR size increases, the CPB value of the HE gene decreases (P = 2x10⁻⁷, ρ = -0.675, n = 63) and that these values correlate positively with the mortality rates caused by the virus (ρ = 0.829, P = 2x10⁻⁷, n = 11). The mortality associated with different virus isolates or their corresponding HPR sizes were not related with the ΔG of HPR RNA folding, suggesting that the secondary structure of HPR RNA does not modify virulence.

CONCLUSIONS

Our results suggest that HPR size affects the efficiency of gene translation, which modulates the virulence of the virus by a mechanism similar to that observed in production of live attenuated vaccines through deoptimization of codon-pair usage.

The rhizome of the multidrug-resistant Enterobacter aerogenes genome reveals how new "killer bugs" are created because of a sympatric lifestyle

Here, we sequenced the 5,419,609 bp circular genome of an Enterobacter aerogenes clinical isolate that killed a patient and was resistant to almost all current antibiotics (except gentamicin) commonly used to treat Enterobacterial infections, including colistin. Genomic and phylogenetic analyses explain the discrepancies of this bacterium and show that its core genome originates from another genus, Klebsiella. Atypical characteristics of this bacterium (i.e., motility, presence of ornithine decarboxylase, and lack of urease activity) are attributed to genomic mosaicism, by acquisition of additional genes, such as the complete 60,582 bp flagellar assembly operon acquired "en bloc" from the genus Serratia. The genealogic tree of the 162,202 bp multidrug-resistant conjugative plasmid shows that it is a chimera of transposons and integrative conjugative elements from various bacterial origins, resembling a rhizome. Moreover, we demonstrate biologically that a G53S mutation in the pmrA gene results in colistin resistance. E. aerogenes has a large RNA population comprising 8 rRNA operons and 87 cognate tRNAs that have the ability to translate transferred genes that use different codons, as exemplified by the significantly different codon usage between genes from the core genome and the "mobilome." On the basis of our findings, the evolution of this bacterium to become a "killer bug" with new genomic repertoires was from three criteria that are "opportunity, power, and usage" to indicate a sympatric lifestyle: "opportunity" to meet other bacteria and exchange foreign sequences since this bacteria was similar to sympatric bacteria; "power" to integrate these foreign sequences such as the acquisition of several mobile genetic elements (plasmids, integrative conjugative element, prophages, transposons, flagellar assembly system, etc.) found in his genome; and "usage" to have the ability to translate these sequences including those from rare codons to serve as a translator of foreign languages.