Acinetobacter junii: an emerging One Health pathogen

Acinetobacter junii is an opportunistic human and animal pathogen severely understudied. Here, we conducted the largest genomic epidemiological study on this pathogen to date. Our data show that this bacterium has spread globally. Also, we found that some human and non-human isolates are not well differentiated from one another, implying transmission between clinical and non-clinical, non-human settings. Remarkably, human but also some non-human isolates have clinically important antibiotic resistance genes, and some of these genes are located in plasmids. Given these results, we put forward that A. junii should be considered an emerging One Health problem. In this regard, future molecular epidemiological studies about this species will go beyond human isolates and will consider animal-, plant-, and water-associated environments.

IMPORTANCE

Acinetobacter baumannii is the most well-known species from the genus Acinetobacter. However, other much less studied Acinetobacter species could be important opportunistic pathogens of animals, plants and humans. Here, we conducted the largest genomic epidemiological study of A. junii, which has been described as a source not only of human but also of animal infections. Our analyses show that this bacterium has spread globally and that, in some instances, human and non-human isolates are not well differentiated. Remarkably, some non-human isolates have important antibiotic resistance genes against important antibiotics used in human medicine. Based on our results, we propose that this pathogen must be considered an issue not only for humans but also for veterinary medicine.

Phylogenomic and phenotypic analyses highlight the diversity of antibiotic resistance and virulence in both human and non-human Acinetobacter baumannii

Acinetobacter baumannii is a Gram-negative, opportunistic pathogen that causes infections in the immunocompromised. With a high incidence of muti-drug resistance, carbapenem-resistant A. baumannii is designated as a priority 1 pathogen by the WHO. The current literature has expertly characterized clinical isolates of A. baumannii. As the challenge of these infections has recently been classified as a One Health issue, we set out to explore the diversity of isolates from human and non-clinical sources, such as agricultural surface water, urban streams, various effluents from wastewater treatment plants, and food (tank milk); and, importantly, these isolates came from a wide geographic distribution. Phylogenomic analysis considering almost 200 isolates showed that our diverse set is well-differentiated from the main international clones of A. baumannii. We discovered novel sequence types in both hospital and non-clinical settings and five strains that overexpress the resistance-nodulation-division efflux pump adeIJK without changes in susceptibility reflected by this overexpression. Furthermore, we detected a bla ADC-79 in a non-human isolate despite its sensitivity to all antibiotics. There was no significant differentiation between the virulence profiles of clinical and non-clinical isolates in the Galleria mellonella insect model of virulence, suggesting that virulence is neither dependent on geographic origin nor isolation source. The detection of antibiotic resistance and virulence genes in non-human strains suggests that these isolates may act as a genetic reservoir for clinical strains. This endorses the notion that in order to combat multi-drug-resistant infection caused by A. baumannii, a One Health approach is required, and a deeper understanding of non-clinical strains must be achieved.IMPORTANCEThe global crisis of antibiotic resistance is a silent one. More and more bacteria are becoming resistant to all antibiotics available for treatment, leaving no options remaining. This includes Acinetobacter baumannii. This Gram-negative, opportunistic pathogen shows a high frequency of multi-drug resistance, and many strains are resistant to the last-resort drugs carbapenem and colistin. Research has focused on strains of clinical origin, but there is a knowledge gap regarding virulence traits, particularly how A. baumannii became the notorious pathogen of today. Antibiotic resistance and virulence genes have been detected in strains from animals and environmental locations such as grass and soil. As such, A. baumannii is a One Health concern, which includes the health of humans, animals, and the environment. Thus, in order to truly combat the antibiotic resistance crisis, we need to understand the antibiotic resistance and virulence gene reservoirs of this pathogen under the One Health continuum.

Unsupervised learning and natural language processing highlight research trends in a superbug

Antibiotic-resistance Acinetobacter baumannii is a very important nosocomial pathogen worldwide. Thousands of studies have been conducted about this pathogen. However, there has not been any attempt to use all this information to highlight the research trends concerning this pathogen. Here we use unsupervised learning and natural language processing (NLP), two areas of Artificial Intelligence, to analyse the most extensive database of articles created (5,500+ articles, from 851 different journals, published over 3 decades). K-means clustering found 113 theme clusters and these were defined with representative terms automatically obtained with topic modelling, summarising different research areas. The biggest clusters, all with over 100 articles, are biased toward multidrug resistance, carbapenem resistance, clinical treatment, and nosocomial infections. However, we also found that some research areas, such as ecology and non-human infections, have received very little attention. This approach allowed us to study research themes over time unveiling those of recent interest, such as the use of cefiderocol (a recently approved antibiotic) against A. baumannii. In a broader context, our results show that unsupervised learning, NLP and topic modelling can be used to describe and analyse the research themes for important infectious diseases. This strategy should be very useful to analyse other ESKAPE pathogens or any other pathogens relevant to Public Health.

Assembly collapsing versus heterozygosity oversizing: detection of homokaryotic and heterokaryotic Laccaria trichodermophora strains by hybrid genome assembly

Genome assembly and annotation using short-paired reads is challenging for eukaryotic organisms due to their large size, variable ploidy and large number of repetitive elements. However, the use of single-molecule long reads improves assembly quality (completeness and contiguity), but haplotype duplications still pose assembly challenges. To address the effect of read length on genome assembly quality, gene prediction and annotation, we compared genome assemblers and sequencing technologies with four strains of the ectomycorrhizal fungus Laccaria trichodermophora. By analysing the predicted repertoire of carbohydrate enzymes, we investigated the effects of assembly quality on functional inferences. Libraries were generated using three different sequencing platforms (Illumina Next-Seq, Mi-Seq and PacBio Sequel), and genomes were assembled using single and hybrid assemblies/libraries. Long reads or hybrid assemby resolved the collapsing of repeated regions, but the nuclear heterozygous versions remained unresolved. In dikaryotic fungi, each cell includes two nuclei and each nucleus has differences not only in allelic gene version but also in gene composition and synteny. These heterokaryotic cells produce fragmentation and size overestimation of the genome assembly of each nucleus. Hybrid assembly revealed a wider functional diversity of genomes. Here, several predicted oxidizing activities on glycosyl residues of oligosaccharides and several chitooligosaccharide acetylase activities would have passed unnoticed in short-read assemblies. Also, the size and fragmentation of the genome assembly, in combination with heterozygosity analysis, allowed us to distinguish homokaryotic and heterokaryotic strains isolated from L. trichodermophora fruit bodies.

Pandrug-resistant Acinetobacter baumannii from different clones and regions in Mexico have a similar plasmid carrying the blaOXA-72 gene

Background

Multidrug-resistant Acinetobacter baumannii is a common hospital-acquired pathogen. The increase in antibiotic resistance is commonly due to the acquisition of mobile genetic elements carrying antibiotic resistance genes. To comprehend this, we analyzed the resistome and virulome of Mexican A. baumannii multidrug-resistant isolates.

Methods

Six clinical strains of A. baumannii from three Mexican hospitals were sequenced using the Illumina platform, the genomes were assembled with SPAdes and annotated with Prokka. Plasmid SPAdes and MobRecon were used to identify the potential plasmid sequences. Sequence Type (ST) assignation under the MLST Oxford scheme was performed using the PubMLST database. Homologous gene search for known virulent factors was performed using the virulence factor database VFDB and an in silico prediction of the resistome was conducted via the ResFinder databases.

Results

The six strains studied belong to different STs and clonal complexes (CC): two strains were ST208 and one was ST369; these two STs belong to the same lineage CC92, which is part of the international clone (IC) 2. Another two strains were ST758 and one was ST1054, both STs belonging to the same lineage CC636, which is within IC5. The resistome analysis of the six strains identified between 7 to 14 antibiotic resistance genes to different families of drugs, including beta-lactams, aminoglycosides, fluoroquinolones and carbapenems. We detected between 1 to 4 plasmids per strain with sizes from 1,800 bp to 111,044 bp. Two strains from hospitals in Mexico City and Guadalajara had a plasmid each of 10,012 bp pAba78r and pAba79f, respectively, which contained the bla OXA-72 gene. The structure of this plasmid showed the same 13 genes in both strains, but 4 of them were inverted in one of the strains. Finally, the six strains contain 49 identical virulence genes related to immune response evasion, quorum-sensing, and secretion systems, among others.

Conclusion

Resistance to carbapenems due to pAba78r and pAba79f plasmids in Aba pandrug-resistant strains from different geographic areas of Mexico and different clones was detected. Our results provide further evidence that plasmids are highly relevant for the horizontal transfer of antibiotic resistance genes between different clones of A. baumannii.

Genomic epidemiology of Acinetobacter baumannii goes global

Acinetobacter baumannii is a major public health concern, for which many genomic epidemiology studies have been conducted in the last decade. However, the vast majority of these are local studies focusing on hospitals from one or a few countries. Proper global genomic epidemiology studies are needed if we are to understand the worldwide dissemination of A. baumannii clones. In this regard, a recent study published in mBio is a good step forward. Müller et al. (mBio e2260-23, 2023, https://doi.org/10.1128/mbio.02260-23) sequenced the genomes of 313 carbapenem-resistant A. baumannii isolates from over 100 hospitals in almost 50 countries from Africa, Asia, Europe, and The Americas. With this data set the authors provide an updated view of the global distribution of the major international clones and their carbapenemase genes. Future global genomic epidemiology studies can be enhanced by considering not only human but also non-human isolates, and by considering isolates despite their antibiotic resistance profile.

Acinetobacter baumannii IC2 and IC5 Isolates with Co-Existing blaOXA-143-like and blaOXA-72 and Exhibiting Strong Biofilm Formation in a Mexican Hospital

Acinetobacter baumannii is an opportunistic pathogen responsible for healthcare-associated infections (HAIs) and outbreaks. Antimicrobial resistance mechanisms and virulence factors allow it to survive and spread in the hospital environment. However, the molecular mechanisms of these traits and their association with international clones are frequently unknown in low- and middle-income countries. Here, we analyze the phenotype and genotype of seventy-six HAIs and outbreak-causing A. baumannii isolates from a Mexican hospital over ten years, with special attention to the carbapenem resistome and biofilm formation. The isolates belonged to the global international clone (IC) 2 and the Latin America endemic IC5 and were predominantly extensively drug-resistant (XDR). Oxacillinases were identified as a common source of carbapenem resistance. We noted the presence of the blaOXA-143-like family (not previously described in Mexico), the blaOXA-72 and the blaOXA-398 found in both ICs. A low prevalence of efflux pump overexpression activity associated with carbapenem resistance was observed. Finally, strong biofilm formation was found, and significant biofilm-related genes were identified, including bfmRS, csuA/BABCDE, pgaABCD and ompA. This study provides a comprehensive profile of the carbapenem resistome of A. baumannii isolates belonging to the same pulse type, along with their significant biofilm formation capacity. Furthermore, it contributes to a better understanding of their role in the recurrence of infection and the endemicity of these isolates in a Mexican hospital.

Emergence of ADC-5 Cephalosporinase in environmental Acinetobacter baumannii from a German tank milk with a novel Sequence Type

Bacteria resistant to antibiotics arguably pose the greatest threat to human health in the twenty-first century. One such bacterium that typifies antibiotic resistance is Acinetobacter baumannii . Frequently, hospital strains of A. baumannii display multidrug resistant (MDR) or extensively drug resistant (XDR) phenotypes, often requiring the use of last resort antibiotics for treatment. In addition to hospital settings, A. baumannii has been isolated from many highly divergent sources including wastewater treatment plant effluent, soil, and agricultural run-off with global distribution. However, such isolates remain poorly characterized. In this study, we characterized a strain of A. baumannii, AB341-IK15, isolated from bulk tank milk in Germany that demonstrated resistance to ceftazidime and intermediate resistance to ceftriaxone and piperacillin/tazobactam. Further genetic characterization identified an ADC-5 cephalosporinase, first incidence in an environmental isolate; and an OXA-408 oxacillinase that may contribute to this phenotype. Interestingly, AB341-IK15 is of a novel sequence type. This research underscores the importance of studying isolates of A. baumannii of non-clinical origin to understand the antibiotic resistance and virulence potential of environmental isolates of A. baumannii as well to understand the diversity of this species.

The Acinetobacter baumannii website (Ab-web): a multidisciplinary knowledge hub, communication platform, and workspace

Acinetobacter baumannii is a Gram-negative bacterium increasingly implicated in hospital-acquired infections and outbreaks. Effective prevention and control of such infections are commonly challenged by the frequent emergence of multidrug-resistant strains. Here we introduce Ab-web (https://www.acinetobacterbaumannii.no), the first online platform for sharing expertise on A. baumannii. Ab-web is a species-centric knowledge hub, initially with 10 articles organized into two main sections, 'Overview' and 'Topics', and three themes, 'epidemiology', 'antibiotic resistance', and 'virulence'. The 'workspace' section provides a spot for colleagues to collaborate, build, and manage joint projects. Ab-web is a community-driven initiative amenable to constructive feedback and new ideas.

Acinetobacter baumannii Sampled from Cattle and Pigs Represent Novel Clones

Acinetobacter baumannii is a very important human pathogen. Nonetheless, we know very little about nonhuman isolates of A. baumannii. Here, we determine the genomic identity of 15 Scottish cattle and pig isolates, as well as their antibiotic and virulence genetic determinants, and compare them with 148 genomes from the main human clinical international clones. Our results demonstrate that cattle and pig isolates represent novel clones well separated from the major international clones. Furthermore, these new clones showed fewer antibiotic resistance genes and may have fewer virulence genes than human clinical isolates. IMPORTANCE Over the last decades, huge amounts of information have been obtained for clinical isolates of A. baumannii and the clones they belong to. In contrast, very little is known about the genomic identity and the genomic basis for virulence and resistance of animal isolates. To fulfil this gap, we conducted a genomic epidemiology study of 15 Scottish cattle and pig isolates in the context of almost 150 genomes belonging to the main international clones of A. baumannii. Our findings show that these animal isolates represent novel clones clearly different from the major international clones. Furthermore, these new clones are distinct in nature considering both antibiotic resistance and virulence when compared with their human clinical counterparts.

Reassessing the evolution of the 1-deoxy-D-xylulose 5-phosphate synthase family suggests a possible novel function for the DXS class 3 proteins

The methylerythritol 4-phosphate (MEP) pathway is of paramount importance for generating plastidial isoprenoids. The first enzyme of the MEP pathway, 1-deoxy-D-xylulose-5-phosphate synthase (DXS), catalyzes a flux-controlling step. In plants the DXS gene family is composed of three distinct classes with non-redundant functions. Although the DXS1 and DXS2 subfamilies have been well characterized, the DXS3 subfamily has been considerably understudied. Here, we carried out in silico and functional analyses to better understand the DXS3 class. Our phylogenetic analysis showed high variation in copy number among the different DXS classes, with the apparent absence of DXS1 class in some species. We found that DXS3 subfamily emerged later than DXS1 and DXS2 and it is under less intense purifying selection. Furthermore, in the DXS3 subfamily critical amino acids positions in the thiamine pyrophosphate binding pocket are not conserved. We demonstrated that the DXS3 proteins from Arabidopsis, Maize, and Rice lack functional DXS activity. Moreover, the Arabidopsis DXS3 protein displayed distinctive sub-organellar chloroplast localization not observed in any DXS1 or DXS2 proteins. Co-expression analysis of the DXS3 from Arabidopsis showed that, unlike DXS1 and DXS2 proteins, it co-expresses with genes related to post-embryonic development and reproduction and not with primary metabolism and isoprenoid synthesis.

A novel vieuvirus from multidrug-resistant Acinetobacter baumannii

Bacteriophages are considered the most abundant biological entities on earth, and they are able to modulate the populations of their bacterial hosts. Although the potential of bacteriophages has been accepted as an alternative strategy to combat multidrug-resistant pathogenic bacteria, there still exists a considerable knowledge gap regarding their genetic diversity, which hinders their use as antimicrobial agents. In this study, we undertook a genomic and phylogenetic characterization of the phage Ab11510-phi, which was isolated from a multidrug-resistant Acinetobacter baumannii strain (Ab11510). We found that Ab11510-phi has a narrow host range and belongs to a small group of transposable phages of the genus Vieuvirus that have only been reported to infect Acinetobacter bacteria. Finally, we showed that Ab11510-phi (as well as other vieuvirus phages) has a high level of mosaicism. On a broader level, we demonstrate that comparative genomics and phylogenetic analysis are necessary tools for the proper characterization of phage diversity.

Spirochetes isolated from arthropods constitute a novel genus Entomospira genus novum within the order Spirochaetales

Spirochetal bacteria were successfully isolated from mosquitoes (Culex pipiens, Aedes cinereus) in the Czech Republic between 1999 and 2002. Preliminary 16S rRNA phylogenetic sequence analysis showed that these strains differed significantly from other spirochetal genera within the family Spirochaetaceae and suggested a novel bacterial genus in this family. To obtain more comprehensive genomic information of these isolates, we used Illumina MiSeq and Oxford Nanopore technologies to sequence four genomes of these spirochetes (BR151, BR149, BR193, BR208). The overall size of the genomes varied between 1.68 and 1.78 Mb; the GC content ranged from 38.5 to 45.8%. Draft genomes were compared to 36 publicly available genomes encompassing eight genera from the class Spirochaetes. A phylogeny generated from orthologous genes across all taxa and the percentage of conserved proteins (POCP) confirmed the genus status of these novel spirochetes. The genus Entomospira gen. nov. is proposed with BR151 selected as type species of the genus. For this isolate and the closest related isolate, BR149, we propose the species name Entomospira culicis sp. nov. The two other isolates BR208 and BR193 are named Entomospira nematocera sp. nov. (BR208) and Entomospira entomophilus sp. nov. (BR193). Finally, we discuss their interesting phylogenetic positioning.

Origin and Evolutionary Dynamics of the miR2119 and ADH1 Regulatory Module in Legumes

MicroRNAs are important regulators of gene expression in eukaryotes. Previously, we reported that in Phaseolus vulgaris, the precursor for miR2119 is located in the same gene as miR398a, conceiving a dicistronic MIR gene. Both miRNA precursors are transcribed and processed from a single transcript resulting in two mature microRNAs that regulate the mRNAs encoding ALCOHOL DEHYDROGENASE 1 (ADH1) and COPPER-ZINC SUPEROXIDE DISMUTASE 1 (CSD1). Genes for miR398 are distributed throughout the spermatophytes; however, miR2119 is only found in Leguminosae species, indicating its recent emergence. Here, we used public databases to explore the presence of the miR2119 sequence in several plant species. We found that miR2119 is present only in specific clades within the Papilionoideae subfamily, including important crops used for human consumption and forage. Within this subfamily, MIR2119 and MIR398a are found together as a single gene in the genomes of the Millettioids and Hologalegina. In contrast, in the Dalbergioids MIR2119 is located in a different locus from MIR398a, suggesting this as the ancestral genomic organization. To our knowledge, this is a unique example where two separate MIRNA genes have merged to generate a single polycistronic gene. Phylogenetic analysis of ADH1 gene sequences in the Papilionoideae subfamily revealed duplication events resulting in up to four ADH1 genes in certain species. Notably, the presence of MIR2119 correlates with the conservation of target sites in particular ADH1 genes in each clade. Our results suggest that post-transcriptional regulation of ADH1 genes by miR2119 has contributed to shaping the expansion and divergence of this gene family in the Papilionoideae. Future experimental work on ADH1 regulation by miR2119 in more legume species will help to further understand the evolutionary history of the ADH1 gene family and the relevance of miRNA regulation in this process.

Chromosome Architecture and Gene Content of the Emergent Pathogen Acinetobacter haemolyticus

Acinetobacter haemolyticus is a Gammaproteobacterium that has been involved in serious diseases frequently linked to the nosocomial environment. Most of the strains causing such infections are sensitive to a wide variety of antibiotics, but recent reports indicate that this pathogen is acquiring very efficiently carbapenem-resistance determinants like the blaNDM-1 gene, all over the world. With this work we contribute with a collection set of 31 newly sequenced nosocomial A. haemolyticus isolates. Genome analysis of these sequences and others collected from RefSeq indicates that their chromosomes are organized in 12 syntenic blocks that contain most of the core genome genes. These blocks are separated by hypervariable regions that are rich in unique gene families, but also have signals of horizontal gene transfer. Genes involved in virulence or encoding different secretion systems are located inside syntenic regions and have recombination signals. The relative order of the synthetic blocks along the A. haemolyticus chromosome can change, indicating that they have been subject to several kinds of inversions. Genomes of this microorganism show large differences in gene content even if they are in the same clade. Here we also show that A. haemolyticus has an open pan-genome.

Structure and Evolution of Acinetobacter baumannii Plasmids

Acinetobacter baumannii is an emergent bacterial pathogen that provokes many types of infections in hospitals around the world. The genome of this organism consists of a chromosome and plasmids. These plasmids vary over a wide size range and many of them have been linked to the acquisition of antibiotic-resistance genes. Our bioinformatic analyses indicate that A. baumannii plasmids belong to a small number of plasmid lineages. The general structure of these lineages seems to be very stable and consists not only of genes involved in plasmid maintenance functions but of gene sets encoding poorly characterized proteins, not obviously linked to survival in the hospital setting, and opening the possibility that they improve the parasitic properties of plasmids. An analysis of genes involved in replication, suggests that members of the same plasmid lineage are part of the same plasmid incompatibility group. The same analysis showed the necessity of classifying the Rep proteins in ten new groups, under the scheme proposed by Bertini et al. (2010). Also, we show that some plasmid lineages have the potential capacity to replicate in many bacterial genera including those embracing human pathogen species, while others seem to replicate only within the limits of the Acinetobacter genus. Moreover, some plasmid lineages are widely distributed along the A. baumannii phylogenetic tree. Despite this, a number of them lack genes involved in conjugation or mobilization functions. Interestingly, only 34.6% of the plasmids analyzed here possess antibiotic resistance genes and most of them belong to fourteen plasmid lineages of the twenty one described here. Gene flux between plasmid lineages appears primarily limited to transposable elements, which sometimes carry antibiotic resistance genes. In most plasmid lineages transposable elements and antibiotic resistance genes are secondary acquisitions. Finally, broad host-range plasmids appear to have played a crucial role.

Phylogenomics Reveals Clear Cases of Misclassification and Genus-Wide Phylogenetic Markers for Acinetobacter

The Gram-negative Acinetobacter genus has several species of clear medical relevance. Many fully sequenced genomes belonging to the genus have been published in recent years; however, there has not been a recent attempt to infer the evolutionary history of Acinetobacter with that vast amount of information. Here, through a phylogenomic approach, we established the most up-to-date view of the evolutionary relationships within this genus and highlighted several cases of poor classification, especially for the very closely related species within the Acinetobacter calcoaceticus-Acinetobacter baumannii complex (Acb complex). Furthermore, we determined appropriate phylogenetic markers for this genus and showed that concatenation of the top 13 gives a very decent reflection of the evolutionary relationships for the genus Acinetobacter. The intersection between our top markers and previously defined universal markers is very small. In general, our study shows that, although there seems to be hardly any universal markers, bespoke phylogenomic approaches can be used to infer the phylogeny of different bacterial genera. We expect that ad hoc phylogenomic approaches will be the standard in the years to come and will provide enough information to resolve intricate evolutionary relationships like those observed in the Acb complex.

High Phenotypic and Genotypic Diversity of Enterococcus faecium from Clinical and Commensal Isolates in Third Level Hospital

Background: The use of antimicrobials and myeloablative chemotherapy regimens has promoted multiresistant microorganisms to emerge as nosocomial pathogens, such as vancomycin-resistant Enterococcus faecium (VREfm). We described a polyclonal outbreak of bloodstream infection caused by Efm in a hemato-oncological ward in Mexico. Our aim was to describe the clonal complex (CC) of the Efm strains isolated in the outbreak in comparison with commensal and environmental isolates. Methodology: Sixty Efm clinical, environmental, and commensal strains were included. We constructed a cladogram and a phylogenetic tree using Vitek and Multilocus sequence typing data, respectively. Results: We reported 20 new sequence types (ST), among which 17/43 clinical isolates belonged to CC17. The predominant ST in the clinical strains were ST757, ST1304, ST412, and ST770. Neither environmental nor commensal isolates belonged to CC17. The phylogeny of our collection shows that the majority of the clinical isolates were different from the environmental and commensal isolates, and only a small group of clinical isolates was closely related with environmental and commensal isolates. The cladogram revealed a similar segregation to that of the phylogeny. Conclusions: We found a high diversity among clinical, environmental, and commensal strains in a group of samples in a single hospital. Highest diversity was found between commensal and environmental isolates.

Inaccurate Multilocus Sequence Typing of Acinetobacter baumannii

Multilocus sequence typing has been useful for genotyping pathogens in surveillance and epidemiologic studies. However, it cannot reflect the true relationships of isolates for species with very dynamic genomes. Using a robust genome phylogeny, we demonstrated the limitations of this method for typing Acinetobacter baumannii.

Phylogenomics picks out the par excellence markers for species phylogeny in the genus Staphylococcus

Although genome sequencing has become a very promising approach to conduct microbial taxonomy, few labs have the resources to afford this especially when dealing with data sets of hundreds to thousands of isolates. The goal of this study was to identify the most adequate loci for inferring the phylogeny of the species within the genus Staphylococcus; with the idea that those who cannot afford whole genome sequencing can use these loci to carry out species assignation confidently. We retrieved 177 orthologous groups (OGs) by using a genome-based phylogeny and an average nucleotide identity analysis. The top 26 OGs showed topologies similar to the species tree and the concatenation of them yielded a topology almost identical to that of the species tree. Furthermore, a phylogeny of just the top seven OGs could be used for species assignment. We sequenced four staphylococcus isolates to test the 26 OGs and found that these OGs were far superior to commonly used markers for this genus. On the whole, our procedure allowed identification of the most adequate markers for inferring the phylogeny within the genus Staphylococcus. We anticipate that this approach will be employed for the identification of the most suitable markers for other bacterial genera and can be very helpful to sort out poorly classified genera.

Phylogenomic Classification and the Evolution of Clonal Complex 5 Methicillin-Resistant Staphylococcus aureus in the Western Hemisphere

Clonal complex 5 methicillin-resistant Staphylococcus aureus (CC5-MRSA) includes multiple prevalent clones that cause hospital-associated infections in the Western Hemisphere. Here, we present a phylogenomic study of these MRSA to reveal their phylogeny, spatial and temporal population structure, and the evolution of selected traits. We studied 598 genome sequences, including 409 newly generated sequences, from 11 countries in Central, North, and South America, and references from Asia and Europe. An early-branching CC5-Basal clade is well-dispersed geographically, is methicillin-susceptible and MRSA predominantly of ST5-IV such as the USA800 clone, and includes separate subclades for avian and porcine strains. In the early 1970s and early 1960s, respectively, two clades appeared that subsequently underwent major expansions in the Western Hemisphere: a CC5-I clade in South America and a CC5-II clade largely in Central and North America. The CC5-I clade includes the ST5-I Chilean/Cordobes clone, and the ST228-I South German clone as an early offshoot, but is distinct from other ST5-I clones from Europe that nest within CC5-Basal. The CC5-II clade includes divergent strains of the ST5-II USA100 clone, various other clones, and most known vancomycin-resistant strains of S. aureus, but is distinct from ST5-II strain N315 from Japan that nests within CC5-Basal. The recombination rate of CC5 was much lower than has been reported for other S. aureus genetic backgrounds, which indicates that recurrence of vancomycin resistance in CC5 is not likely due to an enhanced promiscuity. An increased number of antibiotic resistances and decreased number of toxins with distance from the CC5 tree root were observed. Of note, the expansions of the CC5-I and CC5-II clades in the Western Hemisphere were preceded by convergent gains of resistance to fluoroquinolone, macrolide, and lincosamide antibiotics, and convergent losses of the staphylococcal enterotoxin p (sep) gene from the immune evasion gene cluster of phage ϕSa3. Unique losses of surface proteins were also noted for these two clades. In summary, our study has determined the relationships of different clades and clones of CC5 and has revealed genomic changes for increased antibiotic resistance and decreased virulence associated with the expansions of these MRSA in the Western Hemisphere.

Invasive Methicillin-Resistant Staphylococcus aureus USA500 Strains from the U.S. Emerging Infections Program Constitute Three Geographically Distinct Lineages

USA500 isolates are clonal complex 8 (CC8) Staphylococcus aureus strains closely related to the prominent community- and hospital-associated USA300 group. Despite being relatively understudied, USA500 strains cause a significant burden of disease and are the third most common methicillin-resistant S. aureus (MRSA) strains identified in the U.S. Emerging Infections Program (EIP) invasive S. aureus surveillance. To better understand the genetic relationships of the strains, we sequenced the genomes of 539 USA500 MRSA isolates from sterile site infections collected through the EIP between 2005 and 2013 in the United States. USA500 isolates fell into three major clades principally separated by their distribution across different U.S. regions. Clade C1 strains, found principally in the Northeast, were associated with multiple IS256 insertion elements in their genomes and higher levels of antibiotic resistance. C2 was associated with Southern states, and E1 was associated with Western states. C1 and C2 strains all shared a frameshift in the gene encoding AdsA surface-attached surface protein. We propose that the term "USA500" should be used for CC8 strains sharing a recent common ancestor with the C1, C2, and E1 strains but not in the USA300 group.IMPORTANCE In this work, we have removed some of the confusion surrounding the use of the name "USA500," placed USA500 strains in the context of the CC8 group, and developed a strategy for assignment to subclades based on genome sequence. Our new phylogeny of USA300/USA500 will be a reference point for understanding the genetic adaptations that have allowed multiple highly virulent clonal strains to emerge from within CC8 over the past 50 years.

Rapid Gene Turnover as a Significant Source of Genetic Variation in a Recently Seeded Population of a Healthcare-Associated Pathogen

Genome sequencing has been useful to gain an understanding of bacterial evolution. It has been used for studying the phylogeography and/or the impact of mutation and recombination on bacterial populations. However, it has rarely been used to study gene turnover at microevolutionary scales. Here, we sequenced Mexican strains of the human pathogen Acinetobacter baumannii sampled from the same locale over a 3 year period to obtain insights into the microevolutionary dynamics of gene content variability. We found that the Mexican A. baumannii population was recently founded and has been emerging due to a rapid clonal expansion. Furthermore, we noticed that on average the Mexican strains differed from each other by over 300 genes and, notably, this gene content variation has accrued more frequently and faster than the accumulation of mutations. Moreover, due to its rapid pace, gene content variation reflects the phylogeny only at very short periods of time. Additionally, we found that the external branches of the phylogeny had almost 100 more genes than the internal branches. All in all, these results show that rapid gene turnover has been of paramount importance in producing genetic variation within this population and demonstrate the utility of genome sequencing to study alternative forms of genetic variation.

Recurrent evolution of host and vector association in bacteria of the Borrelia burgdorferi sensu lato species complex

Background

The Borrelia burgdorferi sensu lato (s.l.) species complex consists of tick-transmitted bacteria and currently comprises approximately 20 named and proposed genospecies some of which are known to cause Lyme Borreliosis. Species have been defined via genetic distances and ecological niches they occupy. Understanding the evolutionary relationship of species of the complex is fundamental to explaining patterns of speciation. This in turn forms a crucial basis to frame testable hypotheses concerning the underlying processes including host and vector adaptations.

Results

Illumina Technology was used to obtain genome-wide sequence data for 93 strains of 14 named genospecies of the B. burgdorferi species complex and genomic data already published for 18 additional strain (including one new species) was added. Phylogenetic reconstruction based on 114 orthologous single copy genes shows that the genospecies represent clearly distinguishable taxa with recent and still ongoing speciation events apparent in Europe and Asia. The position of Borrelia species in the phylogeny is consistent with host associations constituting a major driver for speciation. Interestingly, the data also demonstrate that vector associations are an additional driver for diversification in this tick-borne species complex. This is particularly obvious in B. bavariensis, a rodent adapted species that has diverged from the bird-associated B. garinii most likely in Asia. It now consists of two populations one of which most probably invaded Europe following adaptation to a new vector (Ixodes ricinus) and currently expands its distribution range.

Conclusions

The results imply that genotypes/species with novel properties regarding host or vector associations have evolved recurrently during the history of the species complex and may emerge at any time. We suggest that the finding of vector associations as a driver for diversification may be a general pattern for tick-borne pathogens. The core genome analysis presented here provides an important source for investigations of the underlying mechanisms of speciation in tick-borne pathogens.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3016-4) contains supplementary material, which is available to authorized users.

Evolution of a Sigma Factor: An All-In-One of Gene Duplication, Horizontal Gene Transfer, Purifying Selection, and Promoter Differentiation

Sigma factors are an essential part of bacterial gene regulation and have been extensively studied as far as their molecular mechanisms and protein structure are concerned. However, their molecular evolution, especially for the alternative sigma factors, is poorly understood. Here, we analyze the evolutionary forces that have shaped the rpoH sigma factors within the alphaproteobacteria. We found that an ancient duplication gave rise to two major groups of rpoH sigma factors and that after this event horizontal gene transfer (HGT) occurred in rpoH₁ group. We also noted that purifying selection has differentially affected distinct parts of the gene; singularly, the gene segment that encodes the region 4.2, which interacts with the −35 motif of the RpoH-dependent genes, has been under relaxed purifying selection. Furthermore, these two major groups are clearly differentiated from one another regarding their promoter selectivity, as rpoH₁ is under the transcriptional control of σ⁷⁰ and σ³², whereas rpoH₂ is under the transcriptional control of σ²⁴. Our results suggest a scenario in which HGT, gene loss, variable purifying selection and clear promoter specialization occurred after the ancestral duplication event. More generally, our study offers insights into the molecular evolution of alternative sigma factors and highlights the importance of analyzing not only the coding regions but also the promoter regions.

Trans-Atlantic exchanges have shaped the population structure of the Lyme disease agent Borrelia burgdorferi sensu stricto

The origin and population structure of Borrelia burgdorferi sensu stricto (s.s.), the agent of Lyme disease, remain obscure. This tick-transmitted bacterial species occurs in both North America and Europe. We sequenced 17 European isolates (representing the most frequently found sequence types in Europe) and compared these with 17 North American strains. We show that trans-Atlantic exchanges have occurred in the evolutionary history of this species and that a European origin of B. burgdorferi s.s. is marginally more likely than a USA origin. The data further suggest that some European human patients may have acquired their infection in North America. We found three distinct genetically differentiated groups: i) the outgroup species Borrelia bissettii, ii) two divergent strains from Europe, and iii) a group composed of strains from both the USA and Europe. Phylogenetic analysis indicated that different genotypes were likely to have been introduced several times into the same area. Our results demonstrate that irrespective of whether B. burgdorferi s.s. originated in Europe or the USA, later trans-Atlantic exchange(s) have occurred and have shaped the population structure of this genospecies. This study clearly shows the utility of next generation sequencing to obtain a better understanding of the phylogeography of this bacterial species.

Species-specific diversity of novel bacterial lineages and differential abundance of predicted pathways for toxic compound degradation in scorpion gut microbiota

Scorpions are considered 'living fossils' that have conserved ancestral anatomical features and have adapted to numerous habitats. However, their gut microbiota diversity has not been studied. Here, we characterized the gut microbiota of two scorpion species, Vaejovis smithi and Centruroides limpidus. Our results indicate that scorpion gut microbiota is species-specific and that food deprivation reduces bacterial diversity. 16S rRNA gene phylogenetic analysis revealed novel bacterial lineages showing a low level of sequence identity to any known bacteria. Furthermore, these novel bacterial lineages were each restricted to a different scorpion species. Additionally, our results of the predicted metagenomic profiles revealed a core set of pathways that were highly abundant in both species, and mostly related to amino acid, carbohydrate, vitamin and cofactor metabolism. Notably, the food-deprived V. smithi shotgun metagenome matched almost completely the metabolic features of the prediction. Finally, comparisons among predicted metagenomic profiles showed that toxic compound degradation pathways were more abundant in recently captured C. limpidus scorpions. This study gives a first insight into the scorpion gut microbiota and provides a reference for future studies on the gut microbiota from other arachnid species.

RNA-Seq analysis of the multipartite genome of Rhizobium etli CE3 shows different replicon contributions under heat and saline shock

Background

Regulation of transcription is essential for any organism and Rhizobium etli (a multi-replicon, nitrogen-fixing symbiotic bacterium) is no exception. This bacterium is commonly found in the rhizosphere (free-living) or inside of root-nodules of the common bean (Phaseolus vulgaris) in a symbiotic relationship. Abiotic stresses, such as high soil temperatures and salinity, compromise the genetic stability of R. etli and therefore its symbiotic interaction with P. vulgaris. However, it is still unclear which genes are up- or down-regulated to cope with these stress conditions. The aim of this study was to identify the genes and non-coding RNAs (ncRNAs) that are differentially expressed under heat and saline shock, as well as the promoter regions of the up-regulated loci.

Results

Analysing the heat and saline shock responses of R. etli CE3 through RNA-Seq, we identified 756 and 392 differentially expressed genes, respectively, and 106 were up-regulated under both conditions. Notably, the set of genes over-expressed under either condition was preferentially encoded on plasmids, although this observation was more significant for the heat shock response. In contrast, during either saline shock or heat shock, the down-regulated genes were principally chromosomally encoded. Our functional analysis shows that genes encoding chaperone proteins were up-regulated during the heat shock response, whereas genes involved in the metabolism of compatible solutes were up-regulated following saline shock. Furthermore, we identified thirteen and nine ncRNAs that were differentially expressed under heat and saline shock, respectively, as well as eleven ncRNAs that had not been previously identified. Finally, using an in silico analysis, we studied the promoter motifs in all of the non-coding regions associated with the genes and ncRNAs up-regulated under both conditions.

Conclusions

Our data suggest that the replicon contribution is different for different stress responses and that the heat shock response is more complex than the saline shock response. In general, this work exemplifies how strategies that not only consider differentially regulated genes but also regulatory elements of the stress response provide a more comprehensive view of bacterial gene regulation.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-770) contains supplementary material, which is available to authorized users.

Genetic relationships of Vibrio parahaemolyticus isolates from clinical, human carrier, and environmental sources in Thailand, determined by multilocus sequence analysis

Vibrio parahaemolyticus is a seafood-borne pathogenic bacterium that is a major cause of gastroenteritis worldwide. We investigated the genetic and evolutionary relationships of 101 V. parahaemolyticus isolates originating from clinical, human carrier, and various environmental and seafood production sources in Thailand using multilocus sequence analysis. The isolates were recovered from clinical samples (n = 15), healthy human carriers (n = 18), various types of fresh seafood (n = 18), frozen shrimp (n = 16), fresh-farmed shrimp tissue (n = 18), and shrimp farm water (n = 16). Phylogenetic analysis revealed a high degree of genetic diversity within the V. parahaemolyticus population, although isolates recovered from clinical samples and from farmed shrimp and water samples represented distinct clusters. The tight clustering of the clinical isolates suggests that disease-causing isolates are not a random sample of the environmental reservoir, although the source of infection remains unclear. Extensive serotypic diversity occurred among isolates representing the same sequence types and recovered from the same source at the same time. These findings suggest that the O- and K-antigen-encoding loci are subject to exceptionally high rates of recombination. There was also strong evidence of interspecies horizontal gene transfer and intragenic recombination involving the recA locus in a large proportion of isolates. As the majority of the intragenic recombinational exchanges involving recA occurred among clinical and carrier isolates, it is possible that the human intestinal tract serves as a potential reservoir of donor and recipient strains that is promoting horizontal DNA transfer, driving evolutionary change, and leading to the emergence of new, potentially pathogenic strains.

A genomic portrait of the emergence, evolution, and global spread of a methicillin-resistant Staphylococcus aureus pandemic

The widespread use of antibiotics in association with high-density clinical care has driven the emergence of drug-resistant bacteria that are adapted to thrive in hospitalized patients. Of particular concern are globally disseminated methicillin-resistant Staphylococcus aureus (MRSA) clones that cause outbreaks and epidemics associated with health care. The most rapidly spreading and tenacious health-care-associated clone in Europe currently is EMRSA-15, which was first detected in the UK in the early 1990s and subsequently spread throughout Europe and beyond. Using phylogenomic methods to analyze the genome sequences for 193 S. aureus isolates, we were able to show that the current pandemic population of EMRSA-15 descends from a health-care-associated MRSA epidemic that spread throughout England in the 1980s, which had itself previously emerged from a primarily community-associated methicillin-sensitive population. The emergence of fluoroquinolone resistance in this EMRSA-15 subclone in the English Midlands during the mid-1980s appears to have played a key role in triggering pandemic spread, and occurred shortly after the first clinical trials of this drug. Genome-based coalescence analysis estimated that the population of this subclone over the last 20 yr has grown four times faster than its progenitor. Using comparative genomic analysis we identified the molecular genetic basis of 99.8% of the antimicrobial resistance phenotypes of the isolates, highlighting the potential of pathogen genome sequencing as a diagnostic tool. We document the genetic changes associated with adaptation to the hospital environment and with increasing drug resistance over time, and how MRSA evolution likely has been influenced by country-specific drug use regimens.

Phylogeographic variation in recombination rates within a global clone of methicillin-resistant Staphylococcus aureus

BACKGROUND

Next-generation sequencing (NGS) is a powerful tool for understanding both patterns of descent over time and space (phylogeography) and the molecular processes underpinning genome divergence in pathogenic bacteria. Here, we describe a synthesis between these perspectives by employing a recently developed Bayesian approach, BRATNextGen, for detecting recombination on an expanded NGS dataset of the globally disseminated methicillin-resistant Staphylococcus aureus (MRSA) clone ST239.

RESULTS

The data confirm strong geographical clustering at continental, national and city scales and demonstrate that the rate of recombination varies significantly between phylogeographic sub-groups representing independent introductions from Europe. These differences are most striking when mobile non-core genes are included, but remain apparent even when only considering the stable core genome. The monophyletic ST239 sub-group corresponding to isolates from South America shows heightened recombination, the sub-group predominantly from Asia shows an intermediate level, and a very low level of recombination is noted in a third sub-group representing a large collection from Turkey.

CONCLUSIONS

We show that the rapid global dissemination of a single pathogenic bacterial clone results in local variation in measured recombination rates. Possible explanatory variables include the size and time since emergence of each defined sub-population (as determined by the sampling frame), variation in transmission dynamics due to host movement, and changes in the bacterial genome affecting the propensity for recombination.

The impact of recombination on dN/dS within recently emerged bacterial clones

The development of next-generation sequencing platforms is set to reveal an unprecedented level of detail on short-term molecular evolutionary processes in bacteria. Here we re-analyse genome-wide single nucleotide polymorphism (SNP) datasets for recently emerged clones of methicillin resistant Staphylococcus aureus (MRSA) and Clostridium difficile. We note a highly significant enrichment of synonymous SNPs in those genes which have been affected by recombination, i.e. those genes on mobile elements designated “non-core” (in the case of S. aureus), or those core genes which have been affected by homologous replacements (S. aureus and C. difficile). This observation suggests that the previously documented decrease in dN/dS over time in bacteria applies not only to genomes of differing levels of divergence overall, but also to horizontally acquired genes of differing levels of divergence within a single genome. We also consider the role of increased drift acting on recently emerged, highly specialised clones, and the impact of recombination on selection at linked sites. This work has implications for a wide range of genomic analyses.

As bacteria diversify, many of the nucleotide changes that emerge will render the cell slightly less competitive, and these mutations will tend to be removed by natural selection. However, this purging process does not happen instantaneously, and this delay allows deleterious mutations to survive in the population long enough to be sampled. Genomes at the very initial stages of diversification therefore exhibit a relatively high proportion of slightly deleterious mutations and, as most of these are non-synonymous mutations, this is manifest as a high dN/dS ratio. However, the effective population size will also impact on this ratio, as will selection operating on neighbouring SNPs. Here we examine the distribution of synonymous and non-synonymous SNPs within recently emerged clones of two important nosocomial pathogens, methicillin resistant Staphylococcus aureus (MRSA) and Clostridium difficile. In both species, we note a much higher proportion of synonymous changes in those single nucleotide polymorphisms (SNPs) likely to have emerged through recombination compared to de novo mutations. We argue that this effect is explained by the very recent emergence of the mutational SNPs combined with a reduction in the efficiency of selection due to niche specialisation.

Transcriptional profiling of fetal hypothalamic TRH neurons

Background

During murine hypothalamic development, different neuroendocrine cell phenotypes are generated in overlapping periods; this suggests that cell-type specific developmental programs operate to achieve complete maturation. A balance between programs that include cell proliferation, cell cycle withdrawal as well as epigenetic regulation of gene expression characterizes neurogenesis. Thyrotropin releasing hormone (TRH) is a peptide that regulates energy homeostasis and autonomic responses. To better understand the molecular mechanisms underlying TRH neuron development, we performed a genome wide study of its transcriptome during fetal hypothalamic development.

Results

In primary cultures, TRH cells constitute 2% of the total fetal hypothalamic cell population. To purify these cells, we took advantage of the fact that the segment spanning -774 to +84 bp of the Trh gene regulatory region confers specific expression of the green fluorescent protein (GFP) in the TRH cells. Transfected TRH cells were purified by fluorescence activated cell sorting, various cell preparations pooled, and their transcriptome compared to that of GFP- hypothalamic cells. TRH cells undergoing the terminal phase of differentiation, expressed genes implicated in protein biosynthesis, intracellular signaling and transcriptional control. Among the transcription-associated transcripts, we identified the transcription factors Klf4, Klf10 and Atf3, which were previously uncharacterized within the hypothalamus.

Conclusion

To our knowledge, this is one of the first reports identifying transcripts with a potentially important role during the development of a specific hypothalamic neuronal phenotype. This genome-scale study forms a rational foundation for identifying genes that might participate in the development and function of hypothalamic TRH neurons.

Horizontal gene transfer and diverse functional constrains within a common replication-partitioning system in Alphaproteobacteria: the repABC operon

Background

The repABC plasmid family, which is extensively present within Alphaproteobacteria, and some secondary chromosomes of the Rhizobiales have the particular feature that all the elements involved in replication and partitioning reside within one transcriptional unit, the repABC operon. Given the functional interactions among the elements of the repABC operon, and the fact that they all reside in the same operon, a common evolutionary history would be expected if the entire operon had been horizontally transferred. Here, we tested whether there is a common evolutionary history within the repABC operon. We further examined different incompatibility groups in terms of their differentiation and degree of adaptation to their host.

Results

We did not find a single evolutionary history within the repABC operon. Each protein had a particular phylogeny, horizontal gene transfer events of the individual genes within the operon were detected, and different functional constraints were found within and between the Rep proteins. When different repABC operons coexisted in the same genome, they were well differentiated from one another. Finally, we found different levels of adaptation to the host genome within and between repABC operons coexisting in the same species.

Conclusion

Horizontal gene transfer with conservation of the repABC operon structure provides a highly dynamic operon in which each member of this operon has its own evolutionary dynamics. In addition, it seems that different incompatibility groups present in the same species have different degrees of adaptation to their host genomes, in proportion to the amount of time the incompatibility group has coexisted with the host genome.

The Rhizobium etli RpoH1 and RpoH2 sigma factors are involved in different stress responses

The physiological role and transcriptional expression of Rhizobium etli sigma factors rpoH1 and rpoH2 are reported in this work. Both rpoH1 and rpoH2 were able to complement the temperature-sensitive phenotype of an Escherichia coli rpoH mutant. The R. etli rpoH1 mutant was sensitive to heat shock, sodium hypochlorite and hydrogen peroxide, whereas the rpoH2 mutant was sensitive to NaCl and sucrose. The rpoH2 rpoH1 double mutant had increased sensitivity to heat shock and oxidative stress when compared with the rpoH1 single mutant. This suggests that in R. etli, RpoH1 is the main heat-shock sigma factor, but a more complete protective response could be achieved with the participation of RpoH2. Conversely, RpoH2 is involved in osmotic tolerance. In symbiosis with bean plants, the R. etli rpoH1 and rpoH2 rpoH1 mutants still elicited nodule formation, but exhibited reduced nitrogenase activity and bacterial viability in early and late symbiosis compared with nodules produced by rpoH2 mutants and wild-type strains. In addition, nodules formed by R. etli rpoH1 and rpoH2 rpoH1 mutants showed premature senescence. It was also determined that fixNf and fixKf expression was affected in rpoH1 mutants. Both rpoH genes were induced under microaerobic conditions and in the stationary growth phase, but not in response to heat shock. Analysis of the upstream region of rpoH1 revealed a sigma70 and a probable sigmaE promoter, whereas in rpoH2, one probable sigmaE-dependent promoter was detected. In conclusion, the two RpoH proteins operate under different stress conditions, RpoH1 in heat-shock and oxidative responses, and RpoH2 in osmotic tolerance.

Factors affecting the concordance between orthologous gene trees and species tree in bacteria

BACKGROUND

As originally defined, orthologous genes implied a reflection of the history of the species. In recent years, many studies have examined the concordance between orthologous gene trees and species trees in bacteria. These studies have produced contradictory results that may have been influenced by orthologous gene misidentification and artefactual phylogenetic reconstructions. Here, using a method that allows the detection and exclusion of false positives during identification of orthologous genes, we address the question of whether putative orthologous genes within bacteria really reflect the history of the species.

RESULTS

We identified a set of 370 orthologous genes from the bacterial order Rhizobiales. Although manifesting strong vertical signal, almost every orthologous gene had a distinct phylogeny, and the most common topology among the orthologous gene trees did not correspond with the best estimate of the species tree. However, each orthologous gene tree shared an average of 70% of its bipartitions with the best estimate of the species tree. Stochastic error related to gene size affected the concordance between the best estimated of the species tree and the orthologous gene trees, although this effect was weak and distributed unevenly among the functional categories. The nodes showing the greatest discordance were those defined by the shortest internal branches in the best estimated of the species tree. Moreover, a clear bias was evident with respect to the function of the orthologous genes, and the degree of divergence among the orthologous genes appeared to be related to their functional classification.

CONCLUSION

Orthologous genes do not reflect the history of the species when taken as individual markers, but they do when taken as a whole. Stochastic error affected the concordance of orthologous genes with the species tree, albeit weakly. We conclude that two important biological causes of discordance among orthologous genes are incomplete lineage sorting and functional restriction.

A common genomic framework for a diverse assembly of plasmids in the symbiotic nitrogen fixing bacteria

This work centres on the genomic comparisons of two closely-related nitrogen-fixing symbiotic bacteria, Rhizobium leguminosarum biovar viciae 3841 and Rhizobium etli CFN42. These strains maintain a stable genomic core that is also common to other rhizobia species plus a very variable and significant accessory component. The chromosomes are highly syntenic, whereas plasmids are related by fewer syntenic blocks and have mosaic structures. The pairs of plasmids p42f-pRL12, p42e-pRL11 and p42b-pRL9 as well large parts of p42c with pRL10 are shown to be similar, whereas the symbiotic plasmids (p42d and pRL10) are structurally unrelated and seem to follow distinct evolutionary paths. Even though purifying selection is acting on the whole genome, the accessory component is evolving more rapidly. This component is constituted largely for proteins for transport of diverse metabolites and elements of external origin. The present analysis allows us to conclude that a heterogeneous and quickly diversifying group of plasmids co-exists in a common genomic framework.

Rapid evolutionary change of common bean (Phaseolus vulgaris L) plastome, and the genomic diversification of legume chloroplasts

Background

Fabaceae (legumes) is one of the largest families of flowering plants, and some members are important crops. In contrast to what we know about their great diversity or economic importance, our knowledge at the genomic level of chloroplast genomes (cpDNAs or plastomes) for these crops is limited.

Results

We sequenced the complete genome of the common bean (Phaseolus vulgaris cv. Negro Jamapa) chloroplast. The plastome of P. vulgaris is a 150,285 bp circular molecule. It has gene content similar to that of other legume plastomes, but contains two pseudogenes, rpl33 and rps16. A distinct inversion occurred at the junction points of trnH-GUG/rpl14 and rps19/rps8, as in adzuki bean [1]. These two pseudogenes and the inversion were confirmed in 10 varieties representing the two domestication centers of the bean. Genomic comparative analysis indicated that inversions generally occur in legume plastomes and the magnitude and localization of insertions/deletions (indels) also vary. The analysis of repeat sequences demonstrated that patterns and sequences of tandem repeats had an important impact on sequence diversification between legume plastomes and tandem repeats did not belong to dispersed repeats. Interestingly, P. vulgaris plastome had higher evolutionary rates of change on both genomic and gene levels than G. max, which could be the consequence of pressure from both mutation and natural selection.

Conclusion

Legume chloroplast genomes are widely diversified in gene content, gene order, indel structure, abundance and localization of repetitive sequences, intracellular sequence exchange and evolutionary rates. The P. vulgaris plastome is a rapidly evolving genome.