When outgroups fail; phylogenomics of rooting the emerging pathogen, Coxiella burnetii

Community ecology and functional potential of bacteria, archaea, eukarya and viruses in Guerrero Negro microbial mat

In this study, the microbial ecology, potential environmental adaptive mechanisms, and the potential evolutionary interlinking of genes between bacterial, archaeal and viral lineages in Guerrero Negro (GN) microbial mat were investigated using metagenomic sequencing across a vertical transect at millimeter scale. The community composition based on unique genes comprised bacteria (98.01%), archaea (1.81%), eukarya (0.07%) and viruses (0.11%). A gene-focused analysis of bacteria archaea, eukarya and viruses showed a vertical partition of the community. The greatest coverages of genes of bacteria and eukarya were detected in first layers, while the highest coverages of genes of archaea and viruses were found in deeper layers. Many genes potentially related to adaptation to the local environment were detected, such as UV radiation, multidrug resistance, oxidative stress, heavy metals, salinity and desiccation. Those genes were found in bacterial, archaeal and viral lineages with 6477, 44, and 1 genes, respectively. The evolutionary histories of those genes were studied using phylogenetic analysis, showing an interlinking between domains in GN mat.

Genetic characterization and description of Leishmania (Leishmania) ellisi sp. nov.: a new human-infecting species from the USA

In a 2018 report, an unusual case of cutaneous leishmaniasis was described in a 72-year-old female patient residing in Arizona, United States of America (USA). Preliminary analysis of the 18S rDNA and glyceraldehyde-3-phosphate dehydrogenase genes supported the conclusion that the Leishmania strain (strain 218-L139) isolated from this case was a novel species, though a complete taxonomic description was not provided. Identification of Leishmania at the species level is critical for clinical management and epidemiologic investigations so it is important that novel human-infecting species are characterized taxonomically and assigned a unique scientific name compliant with the ICZN code. Therefore, we sought to provide a complete taxonomic description of Leishmania strain 218-L139. Phylogenetic analysis of several nuclear loci and partial maxicircle genome sequences supported its position within the subgenus Leishmania and further clarified the distinctness of this new species. Morphological characterization of cultured promastigotes and amastigotes from the original case material is also provided. Thus, we conclude that Leishmania (Leishmania) ellisi is a new cause of autochthonous cutaneous leishmaniasis in the USA.

A critical analysis of the current state of virus taxonomy

Taxonomical classification has preceded evolutionary understanding. For that reason, taxonomy has become a battleground fueled by knowledge gaps, technical limitations, and a priorism. Here we assess the current state of the challenging field, focusing on fallacies that are common in viral classification. We emphasize that viruses are crucial contributors to the genomic and functional makeup of holobionts, organismal communities that behave as units of biological organization. Consequently, viruses cannot be considered taxonomic units because they challenge crucial concepts of organismality and individuality. Instead, they should be considered processes that integrate virions and their hosts into life cycles. Viruses harbor phylogenetic signatures of genetic transfer that compromise monophyly and the validity of deep taxonomic ranks. A focus on building phylogenetic networks using alignment-free methodologies and molecular structure can help mitigate the impasse, at least in part. Finally, structural phylogenomic analysis challenges the polyphyletic scenario of multiple viral origins adopted by virus taxonomy, defeating a polyphyletic origin and supporting instead an ancient cellular origin of viruses. We therefore, prompt abandoning deep ranks and urgently reevaluating the validity of taxonomic units and principles of virus classification.

DISCO+QR: rooting species trees in the presence of GDL and ILS

Motivation

Genes evolve under processes such as gene duplication and loss (GDL), so that gene family trees are multi-copy, as well as incomplete lineage sorting (ILS); both processes produce gene trees that differ from the species tree. The estimation of species trees from sets of gene family trees is challenging, and the estimation of rooted species trees presents additional analytical challenges. Two of the methods developed for this problem are STRIDE, which roots species trees by considering GDL events, and Quintet Rooting (QR), which roots species trees by considering ILS.

Results

We present DISCO+QR, a new approach to rooting species trees that first uses DISCO to address GDL and then uses QR to perform rooting in the presence of ILS. DISCO+QR operates by taking the input gene family trees and decomposing them into single-copy trees using DISCO and then roots the given species tree using the information in the single-copy gene trees using QR. We show that the relative accuracy of STRIDE and DISCO+QR depend on the properties of the dataset (number of species, genes, rate of gene duplication, degree of ILS and gene tree estimation error), and that each provides advantages over the other under some conditions.

Availability and implementation

DISCO and QR are available in github.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Characterization and Identification of Probiotic Features in Lacticaseibacillus Paracasei Using a Comparative Genomic Analysis Approach

Lacticaseibacillus paracasei species are widely used for their health-promoting properties in food and agricultural applications. These bacteria have been isolated from various habitats such as the oral cavity, cereals, vegetables, meats, and dairy products conferring them the ability to consume different carbohydrates. Two subspecies are recognized, Lacticaseibacillus paracasei subsp. paracasei and Lacticaseibacillus paracasei subsp. tolerans according to their acid production from carbohydrates. Some strains are currently used as probiotics. In this study, we performed a comparative genomic analysis of 181 genomes of the Lacticaseibacillus paracasei species to reveal genomic differences at the subspecies level and to reveal adaptive and probiotic features, and special emphasis is given to inulin consumption. No clear distinction at the subspecies level for L. paracasei was shown using a phylogenetic tree with orthologous genes from the core-genome set. In general, a good correlation was observed between genomic distance and isolation origin, suggesting that L. paracasei strains are adapted to their natural habitat, giving rise to genetic differences at the genomic level. A low frequency of undesirable characteristics such as plasmids, prophages, antibiotic resistance genes, absence of virulence factors, and frequent bacteriocin production supports these species being good candidates for use as probiotics. Lastly, we found that the inulin gene cluster in L. paracasei strains seems to differ slightly in the presence or absence of some genes but maintains a core defined by at least three fructose-PTS proteins, one hypothetical protein, and extracellular β-fructosidase. Finally, we conclude that further work has to be done for L. paracasei subspecies classification. Improving outgroup selection criteria is a key factor for their correct subspecies assignation.

Correlating Genotyping Data of Coxiella burnetii with Genomic Groups

Coxiella burnetii is a zoonotic pathogen that resides in wild and domesticated animals across the globe and causes a febrile illness, Q fever, in humans. Several distinct genetic lineages or genomic groups have been shown to exist, with evidence for different virulence potential of these lineages. Multispacer Sequence Typing (MST) and Multiple-Locus Variable number tandem repeat Analysis (MLVA) are being used to genotype strains. However, it is unclear how these typing schemes correlate with each other or with the classification into different genomic groups. Here, we created extensive databases for published MLVA and MST genotypes of C. burnetii and analysed the associated metadata, revealing associations between animal host and human disease type. We established a new classification scheme that assigns both MST and MLVA genotypes to a genomic group and which revealed additional sub-lineages in two genomic groups. Finally, we report a novel, rapid genomotyping method for assigning an isolate into a genomic group based on the Cox51 spacer sequence. We conclude that by pooling and streamlining existing datasets, associations between genotype and clinical outcome or host source were identified, which in combination with our novel genomotyping method, should enable an estimation of the disease potential of new C. burnetii isolates.

Pathogen to commensal? Longitudinal within-host population dynamics, evolution, and adaptation during a chronic >16-year Burkholderia pseudomallei infection

Although acute melioidosis is the most common outcome of Burkholderia pseudomallei infection, we have documented a case, P314, where disease severity lessened with time, and the pathogen evolved towards a commensal relationship with the host. In the current study, we used whole-genome sequencing to monitor this long-term symbiotic relationship to better understand B. pseudomallei persistence in P314's sputum despite intensive initial therapeutic regimens. We collected and sequenced 118 B. pseudomallei isolates from P314's airways over a >16-year period, and also sampled the patient's home environment, recovering six closely related B. pseudomallei isolates from the household water system. Using comparative genomics, we identified 126 SNPs in the core genome of the 124 isolates or 162 SNPs/indels when the accessory genome was included. The core SNPs were used to construct a phylogenetic tree, which demonstrated a close relationship between environmental and clinical isolates and detailed within-host evolutionary patterns. The phylogeny had little homoplasy, consistent with a strictly clonal mode of genetic inheritance. Repeated sampling revealed evidence of genetic diversification, but frequent extinctions left only one successful lineage through the first four years and two lineages after that. Overall, the evolution of this population is nonadaptive and best explained by genetic drift. However, some genetic and phenotypic changes are consistent with in situ adaptation. Using a mouse model, P314 isolates caused greatly reduced morbidity and mortality compared to the environmental isolates. Additionally, potentially adaptive phenotypes emerged and included differences in the O-antigen, capsular polysaccharide, motility, and colony morphology. The >13-year co-existence of two long-lived lineages presents interesting hypotheses that can be tested in future studies to provide additional insights into selective pressures, niche differentiation, and microbial adaptation. This unusual melioidosis case presents a rare example of the evolutionary progression towards commensalism by a highly virulent pathogen within a single human host.

Selective whole genome amplification and sequencing of Coxiella burnetii directly from environmental samples

Whole genome sequencing (WGS) is a widely available, inexpensive means of providing a wealth of information about an organism's diversity and evolution. However, WGS for many pathogenic bacteria remain limited because they are difficult, slow and/or dangerous to culture. To avoid culturing, metagenomic sequencing can be performed directly on samples, but the sequencing effort required to characterize low frequency organisms can be expensive. Recently developed methods for selective whole genome amplification (SWGA) can enrich target DNA to provide efficient sequencing. We amplified Coxiella burnetii (a bacterial select agent and human/livestock pathogen) from 3 three environmental samples that were overwhelmed with host DNA. The 68- to 147-fold enrichment of the bacterial sequences provided enough genome coverage for SNP analyses and phylogenetic placement. SWGA is a valuable tool for the study of difficult-to-culture organisms and has the potential to facilitate high-throughput population characterizations as well as targeted epidemiological or forensic investigations.

From cell culture to cynomolgus macaque: infection models show lineage-specific virulence potential of Coxiella burnetii

Coxiella burnetii is an obligate intracellular pathogen that causes the zoonotic disease Q fever in humans, which can occur in either an acute or a chronic form with serious complications. The bacterium has a wide host range, including unicellular organisms, invertebrates, birds and mammals, with livestock representing the most significant reservoir for human infections. Cell culture models have been used to decipher the intracellular lifestyle of C. burnetii, and several infection models, including invertebrates, rodents and non-human primates, are being used to investigate host-pathogen interactions and to identify bacterial virulence factors and vaccine candidates. However, none of the models replicate all aspects of human disease. Furthermore, it is becoming evident that C. burnetii isolates belonging to different lineages exhibit differences in their virulence in these models. Here, we compare the advantages and disadvantages of commonly used infection models and summarize currently available data for lineage-specific virulence.

Extensive genome analysis of Coxiella burnetii reveals limited evolution within genomic groups

BACKGROUND

Coxiella burnetii is a zoonotic pathogen that resides in wild and domesticated animals across the globe and causes a febrile illness, Q fever, in humans. An improved understanding of the genetic diversity of C. burnetii is essential for the development of diagnostics, vaccines and therapeutics, but genotyping data is lacking from many parts of the world. Sporadic outbreaks of Q fever have occurred in the United Kingdom, but the local genetic make-up of C. burnetii has not been studied in detail.

RESULTS

Here, we report whole genome data for nine C. burnetii sequences obtained in the UK. All four genomes of C. burnetii from cattle, as well as one sheep sample, belonged to Multi-spacer sequence type (MST) 20, whereas the goat samples were MST33 (three genomes) and MST32 (one genome), two genotypes that have not been described to be present in the UK to date. We established the phylogenetic relationship between the UK genomes and 67 publically available genomes based on single nucleotide polymorphisms (SNPs) in the core genome, which confirmed tight clustering of strains within genomic groups, but also indicated that sub-groups exist within those groups. Variation is mainly achieved through SNPs, many of which are non-synonymous, thereby confirming that evolution of C. burnetii is based on modification of existing genes. Finally, we discovered genomic-group specific genome content, which supports a model of clonal expansion of previously established genotypes, with large scale dissemination of some of these genotypes across continents being observed.

CONCLUSIONS

The genetic make-up of C. burnetii in the UK is similar to the one in neighboring European countries. As a species, C. burnetii has been considered a clonal pathogen with low genetic diversity at the nucleotide level. Here, we present evidence for significant variation at the protein level between isolates of different genomic groups, which mainly affects secreted and membrane-associated proteins. Our results thereby increase our understanding of the global genetic diversity of C. burnetii and provide new insights into the evolution of this emerging zoonotic pathogen.

Rooting Phylogenies and the Tree of Life While Minimizing Ad Hoc and Auxiliary Assumptions

Phylogenetic methods unearth evolutionary history when supported by three starting points of reason: (1) the continuity axiom begs the existence of a "model" of evolutionary change, (2) the singularity axiom defines the historical ground plan (phylogeny) in which biological entities (taxa) evolve, and (3) the memory axiom demands identification of biological attributes (characters) with historical information. Axiom consequences are interlinked, making the retrodiction enterprise an endeavor of reciprocal fulfillment. In particular, establishing direction of evolutionary change (character polarization) roots phylogenies and enables testing the existence of historical memory (homology). Unfortunately, rooting phylogenies, especially the "tree of life," generally follow narratives instead of integrating empirical and theoretical knowledge of retrodictive exploration. This stems mostly from a focus on molecular sequence analysis and uncertainties about rooting methods. Here, we review available rooting criteria, highlighting the need to minimize both ad hoc and auxiliary assumptions, especially argumentative ad hocness. We show that while the outgroup comparison method has been widely adopted, the generality criterion of nesting and additive phylogenetic change embodied in Weston rule offers the most powerful rooting approach. We also propose a change of focus, from phylogenies that describe the evolution of biological systems to those that describe the evolution of parts of those systems. This weakens violation of character independence, helps formalize the generality criterion of rooting, and provides new ways to study the problem of evolution.

The Importance of Revisiting Legionellales Diversity

Bacteria of the order Legionellales, such as Legionella pneumophila, the agent of Legionnaires' disease, and Coxiella burnetii, the agent of Q fever, are widely recognized as human pathogens. While our view of the Legionellales is often limited to clinical isolates, ecological surveys are continually uncovering new members of the Legionellales that do not fall into the recognized pathogenic species. Here we emphasize that most of these Legionellales are nonpathogenic forms that have evolved symbiotic lifestyles with nonvertebrate hosts. The diversity of nonpathogenic forms remains, however, largely underexplored. We conjecture that its characterization, once contrasted with the data on pathogenic species, will reveal novel highlights on the mechanisms underlying lifestyle transitions of intracellular bacteria, including the emergence of pathogenesis and mutualism, transmission routes, and host specificity.

Comparative pan-genomic analyses of Orientia tsutsugamushi reveal an exceptional model of bacterial evolution driving genomic diversity

Orientia tsutsugamushi, formerly Rickettsia tsutsugamushi, is an obligate intracellular pathogen that causes scrub typhus, an underdiagnosed acute febrile disease with high morbidity. Scrub typhus is transmitted by the larval stage (chigger) of Leptotrombidium mites and is irregularly distributed across endemic regions of Asia, Australia and islands of the western Pacific Ocean. Previous work to understand population genetics in O. tsutsugamushi has been based on sub-genomic sampling methods and whole-genome characterization of two genomes. In this study, we compared 40 genomes from geographically dispersed areas and confirmed patterns of extensive homologous recombination likely driven by transposons, conjugative elements and repetitive sequences. High rates of lateral gene transfer (LGT) among O. tsutsugamushi genomes appear to have effectively eliminated a detectable clonal frame, but not our ability to infer evolutionary relationships and phylogeographical clustering. Pan-genomic comparisons using 31 082 high-quality bacterial genomes from 253 species suggests that genomic duplication in O. tsutsugamushi is almost unparalleled. Unlike other highly recombinant species where the uptake of exogenous DNA largely drives genomic diversity, the pan-genome of O. tsutsugamushi is driven by duplication and divergence. Extensive gene innovation by duplication is most commonly attributed to plants and animals and, in contrast with LGT, is thought to be only a minor evolutionary mechanism for bacteria. The near unprecedented evolutionary characteristics of O. tsutsugamushi, coupled with extensive intra-specific LGT, expand our present understanding of rapid bacterial evolutionary adaptive mechanisms.

Archaea-First and the Co-Evolutionary Diversification of Domains of Life

The origins and evolution of the Archaea, Bacteria, and Eukarya remain controversial. Phylogenomic-wide studies of molecular features that are evolutionarily conserved, such as protein structural domains, suggest Archaea is the first domain of life to diversify from a stem line of descent. This line embodies the last universal common ancestor of cellular life. Here, we propose that ancestors of Euryarchaeota co-evolved with those of Bacteria prior to the diversification of Eukarya. This co-evolutionary scenario is supported by comparative genomic and phylogenomic analyses of the distributions of fold families of domains in the proteomes of free-living organisms, which show horizontal gene recruitments and informational process homologies. It also benefits from the molecular study of cell physiologies responsible for membrane phospholipids, methanogenesis, methane oxidation, cell division, gas vesicles, and the cell wall. Our theory however challenges popular cell fusion and two-domain of life scenarios derived from sequence analysis, demanding phylogenetic reconciliation. Also see the video abstract here: https://youtu.be/9yVWn_Q9faY.

Phylogenetic inference of Coxiella burnetii by 16S rRNA gene sequencing

Coxiella burnetii is a human pathogen that causes the serious zoonotic disease Q fever. It is ubiquitous in the environment and due to its wide host range, long-range dispersal potential and classification as a bioterrorism agent, this microorganism is considered an HHS Select Agent. In the event of an outbreak or intentional release, laboratory strain typing methods can contribute to epidemiological investigations, law enforcement investigation and the public health response by providing critical information about the relatedness between C. burnetii isolates collected from different sources. Laboratory cultivation of C. burnetii is both time-consuming and challenging. Availability of strain collections is often limited and while several strain typing methods have been described over the years, a true gold-standard method is still elusive. Building upon epidemiological knowledge from limited, historical strain collections and typing data is essential to more accurately infer C. burnetii phylogeny. Harmonization of auspicious high-resolution laboratory typing techniques is critical to support epidemiological and law enforcement investigation. The single nucleotide polymorphism (SNP) -based genotyping approach offers simplicity, rapidity and robustness. Herein, we demonstrate SNPs identified within 16S rRNA gene sequences can differentiate C. burnetii strains. Using this method, 55 isolates were assigned to six groups based on six polymorphisms. These 16S rRNA SNP-based genotyping results were largely congruent with those obtained by analyzing restriction-endonuclease (RE)-digested DNA separated by SDS-PAGE and by the high-resolution approach based on SNPs within multispacer sequence typing (MST) loci. The SNPs identified within the 16S rRNA gene can be used as targets for the development of additional SNP-based genotyping assays for C. burnetii.

Draft Genome Sequences of the Avirulent Coxiella burnetii Dugway 7D77-80 and Dugway 7E65-68 Strains Isolated from Rodents in Dugway, Utah

Here, we present the draft genome sequences of the Coxiella burnetii Dugway 7D77-80 and Dugway 7E65-68 strains, which were isolated from rodents in Dugway, UT, in the 1950s. The strains reside in a distinct genomic group of C. burnetii and are considered avirulent despite having the largest genomes of the Coxiella genus.

Minimum variance rooting of phylogenetic trees and implications for species tree reconstruction

Phylogenetic trees inferred using commonly-used models of sequence evolution are unrooted, but the root position matters both for interpretation and downstream applications. This issue has been long recognized; however, whether the potential for discordance between the species tree and gene trees impacts methods of rooting a phylogenetic tree has not been extensively studied. In this paper, we introduce a new method of rooting a tree based on its branch length distribution; our method, which minimizes the variance of root to tip distances, is inspired by the traditional midpoint rerooting and is justified when deviations from the strict molecular clock are random. Like midpoint rerooting, the method can be implemented in a linear time algorithm. In extensive simulations that consider discordance between gene trees and the species tree, we show that the new method is more accurate than midpoint rerooting, but its relative accuracy compared to using outgroups to root gene trees depends on the size of the dataset and levels of deviations from the strict clock. We show high levels of error for all methods of rooting estimated gene trees due to factors that include effects of gene tree discordance, deviations from the clock, and gene tree estimation error. Our simulations, however, did not reveal significant differences between two equivalent methods for species tree estimation that use rooted and unrooted input, namely, STAR and NJst. Nevertheless, our results point to limitations of existing scalable rooting methods.

Sequencing, characterization and phylogenomics of the complete mitochondrial genome of Dactylogyrus lamellatus (Monogenea: Dactylogyridae)

Despite the worldwide distribution and pathogenicity of monogenean parasites belonging to the largest helminth genus, Dactylogyrus, there are no complete Dactylogyrinae (subfamily) mitogenomes published to date. In order to fill this knowledge gap, we have sequenced and characterized the complete mitogenome of Dactylogyrus lamellatus, a common parasite on the gills of grass carp (Ctenopharyngodon idella). The circular mitogenome is 15,187 bp in size, containing the standard 22 tRNA genes, 2 rRNA genes, 12 protein-encoding genes and a long non-coding region (NCR). There are two highly repetitive regions in the NCR. We have used concatenated nucleotide sequences of all 36 genes to perform the phylogenetic analysis using Bayesian inference and maximum likelihood approaches. As expected, the two dactylogyrids, D. lamellatus (Dactylogyrinae) and Tetrancistrum nebulosi (Ancyrocephalinae), were closely related to each other. These two formed a sister group with Capsalidae, and this cluster finally formed a further sister group with Gyrodactylidae. Phylogenetic affinity between Dactylogyrinae and Ancyrocephalinae was further confirmed by the similarity in their gene arrangement. The sequencing of the first Dactylogyrinae, along with a more suitable selection of outgroups, has enabled us to infer a much better phylogenetic resolution than recent mitogenomic studies. However, as many lineages of the class Monogenea remain underrepresented or not represented at all, a much larger number of mitogenome sequences will have to be available in order to infer the evolutionary relationships among the monogeneans fully, and with certainty.

Genome-wide characterization of cellulases from the hemi-biotrophic plant pathogen, Bipolaris sorokiniana, reveals the presence of a highly stable GH7 endoglucanase

BACKGROUND

Bipolaris sorokiniana is a filamentous fungus that causes spot blotch disease in cereals like wheat and has severe economic consequences. However, information on the identities and role of the cell wall-degrading enzymes (CWDE) in B. sorokiniana is very limited. Several fungi produce CWDE like glycosyl hydrolases (GHs) that help in host cell invasion. To understand the role of these CWDE in B. sorokiniana, the first step is to identify and annotate all possible genes of the GH families like GH3, GH6, GH7, GH45 and AA9 and then characterize them biochemically.

RESULTS

We confirmed and annotated the homologs of GH3, GH6, GH7, GH45 and AA9 enzymes in the B. sorokiniana genome using the sequence and domain features of these families. Quantitative real-time PCR analyses of these homologs revealed that the transcripts of the BsGH7-3 (3rd homolog of the GH 7 family in B. sorokiniana) were most abundant. BsGH7-3, the gene of BsGH7-3, was thus cloned into pPICZαC Pichia pastoris vector and expressed in X33 P. pastoris host to be characterized. BsGH7-3 enzyme showed a temperature optimum of 60 °C and a pH_opt of 8.1. BsGH7-3 was identified to be an endoglucanase based on its broad substrate specificity and structural comparisons with other such endoglucanases. BsGH7-3 has a very long half-life and retains 100% activity even in the presence of 4 M NaCl, 4 M KCl and 20% (v/v) ionic liquids. The enzyme activity is stimulated up to fivefold in the presence of Mn⁺² and Fe⁺² without any deleterious effects on enzyme thermostability.

CONCLUSIONS

Here we reanalysed the B. sorokiniana genome and selected one GH7 enzyme for further characterization. The present work demonstrates that BsGH7-3 is an endoglucanase with a long half-life and no loss in activity in the presence of denaturants like salt and ionic liquids, and lays the foundation towards exploring the Bipolaris genome for other cell wall-degrading enzymes.

Phylogenetic and Evolutionary Analysis of Plant ARGONAUTES

Comparative sequence analysis is widely used for the reconstruction of phylogeny and for understanding the evolutionary history of gene families. Here, we describe the methodologies to reconstruct the phylogenetic and evolutionary history of a gene family across genomes with a focus on the ARGONAUTE (AGO) family of proteins in plants. The method described here may easily be adapted for studying molecular evolution of a wide variety of gene families. We enlist methods as well as parameters for the collection of molecular data (nucleic acids and peptides), preparation of datasets, and selection of evolutionary models and various methods for the phylogenetic and evolutionary analysis, such as maximum likelihood and Bayesian inference.

Coxiella burnetii in central Italy: novel genotypes are circulating in cattle and goats

Genotyping of bacteria is critical for diagnosis, treatment, and epidemiological surveillance. Coxiella burnetii, the etiological agent of Q fever, has been recognized to have a potential for bioterrorism purposes. Because few serosurveys have been conducted in Italy, there is still limited information about the distribution of this pathogen in natural conditions. In this paper, we describe the genotyping of C. burnetii strains by multispacer sequence typing (MST) detected in cattle and goat farms in the Abruzzi region of Italy. Biological samples (milk, aborted fetus) positive for C. burnetii DNA were sequenced in the spacer regions and compared with those already publicly available ( http://ifr48.timone.univ-mrs.fr/MST_Coxiella/mst/group_detail ). The MST profile of C. burnetii detected in milk samples demonstrated the presence of a new allele, whereas the C. burnetii spacer sequences from fetus and milk goat samples displayed a new allelic combination. The results suggest the circulation of novel genotypes of C. burnetii in Italy.

Massive dispersal of Coxiella burnetii among cattle across the United States

Q-fever is an underreported disease caused by the bacterium Coxiella burnetii, which is highly infectious and has the ability to disperse great distances. It is a completely clonal pathogen with low genetic diversity and requires whole-genome analysis to identify discriminating features among closely related isolates. C. burnetii, and in particular one genotype (ST20), is commonly found in cow's milk across the entire dairy industry of the USA. This single genotype dominance is suggestive of host-specific adaptation, rapid dispersal and persistence within cattle. We used a comparative genomic approach to identify SNPs for high-resolution and high-throughput genotyping assays to better describe the dispersal of ST20 across the USA. We genotyped 507 ST20 cow milk samples and discovered three subgenotypes, all of which were present across the entire country and over the complete time period studied. Only one of these sub-genotypes was observed in a single dairy herd. The temporal and geographic distribution of these sub-genotypes is consistent with a model of large-scale, rapid, frequent and continuous dissemination on a continental scale. The distribution of subgenotypes is not consistent with wind-based dispersal alone, and it is likely that animal husbandry and transportation practices, including pooling of milk from multiple herds, have also shaped the patterns. On the scale of an entire country, there appear to be few barriers to rapid, frequent and large-scale dissemination of the ST20 subgenotypes.

Right on Q: genetics begin to unravel Coxiella burnetii host cell interactions

Invasion of macrophages and replication within an acidic and degradative phagolysosome-like vacuole are essential for disease pathogenesis by Coxiella burnetii, the bacterial agent of human Q fever. Previous experimental constraints imposed by the obligate intracellular nature of Coxiella limited knowledge of pathogen strategies that promote infection. Fortunately, new genetic tools facilitated by axenic culture now allow allelic exchange and transposon mutagenesis approaches for virulence gene discovery. Phenotypic screens have illuminated the critical importance of Coxiella's type 4B secretion system in host cell subversion and discovered genes encoding translocated effector proteins that manipulate critical infection events. Here, we highlight the cellular microbiology and genetics of Coxiella and how recent technical advances now make Coxiella a model organism to study macrophage parasitism.

Genotyping and Axenic Growth of Coxiella burnetii Isolates Found in the United States Environment

Coxiella burnetii is a gram-negative bacterium that is the etiologic agent of the zoonotic disease Q fever. Common reservoirs of C. burnetii include sheep, goats, and cattle. These animals shed C. burnetii into the environment, and humans are infected by inhalation of aerosols. A survey of 1622 environmental samples taken across the United States in 2006-2008 found that 23.8% of the samples contained C. burnetii DNA. To identify the strains circulating in the U.S. environment, DNA from these environmental samples was genotyped using an SNP-based approach to derive sequence types (ST) that are also compatible with multispacer sequence typing methods. Three different sequence types were observed in 31 samples taken from 19 locations. ST8 was associated with goats and ST20 with dairy cattle. ST16/26 was detected in locations with exposure to various animals and also in locations with no direct animal contact. Viable isolates were obtained for all three sequence types, but only the ST20 and ST16/26 isolates grew in acidified citrate cysteine medium (ACCM)-2 axenic media. Examination of a variety of isolates with different sequence types showed that ST8 and closely related isolates did not grow in ACCM-2. These results suggest that a limited number of C. burnetii sequence types are circulating in the U.S. environment and these strains have close associations with specific reservoir species. Growth in ACCM-2 may not be suitable for isolation of many C. burnetii strains.

A New Orthology Assessment Method for Phylogenomic Data: Unrooted Phylogenetic Orthology

Current sequencing technologies are making available unprecedented amounts of genetic data for a large variety of species including nonmodel organisms. Although many phylogenomic surveys spend considerable time finding orthologs from the wealth of sequence data, these results do not transcend the original study and after being processed for specific phylogenetic purposes these orthologs do not become stable orthology hypotheses. We describe a procedure to detect and document the phylogenetic distribution of orthologs allowing researchers to use this information to guide selection of loci best suited to test specific evolutionary questions. At the core of this pipeline is a new phylogenetic orthology method that is neither affected by the position of the root nor requires explicit assignment of outgroups. We discuss the properties of this new orthology assessment method and exemplify its utility for phylogenomics using a small insects dataset. In addition, we exemplify the pipeline to identify and document stable orthologs for the group of orb-weaving spiders (Araneoidea) using RNAseq data. The scripts used in this study, along with sample files and additional documentation, are available at https://github.com/ballesterus/UPhO.

Rooting Trees, Methods for

Identifying the root of a phylogenetic tree is important because incorrectly rooted phylogenetic trees may mislead evolutionary and taxonomic inferences. Many techniques for inferring the root have been proposed, but each has shortcomings that may make it inappropriate for any particular dataset. Here we outline the various ways to root phylogenetic trees, which include: outgroup, midpoint rooting, molecular clock rooting, and Bayesian molecular clock rooting. In addition, we discuss the pros and cons and also list software availability for each of the rooting methods.

[Not Available]

Q fever is a zoonosis of worldwide distribution with the exception of New Zealand. It is caused by an intracellular bacterium, Coxiella burnetii. The disease often goes underdiagnosed because the main manifestation of its acute form is a general self-limiting flu-like syndrome. The Dutch epidemics renewed attention to this disease, which was less considered before. This review summarizes the description of C. burnetii (taxonomy, intracellular cycle, and genome) and Q fever disease (description, diagnosis, epidemiology, and pathogenesis). Finally, vaccination in humans and animals is also considered.

Coxiella-like endosymbiont in argasid ticks (Ornithodoros muesebecki) from a Socotra Cormorant colony in Umm Al Quwain, United Arab Emirates

Coxiella burnetii is a pathogen causing Q fever in domestic animals and humans. Seabirds have been implicated as possible reservoirs of this bacterium in the Arabian Gulf and in the Western Indian Ocean. Recently, Coxiella species closely related to C. burnetii was detected from ticks collected from oil rigs used as roosting areas by Socotra Cormorants (Phalacrocorax nigrogularis) in the western Arabian Gulf. We collected ticks from the largest breeding colony of Socotra Cormorants in the United Arab Emirates on the eastern extreme of the species' breeding range to determine the prevalence of C. burnetii and evaluate its role as a wild reservoir. All ticks were identified as Ornithodoros muesebecki and genomic DNA was extracted from larval and nymph/adult tick pools. Multiplex PCR tests were performed targeting three C. burnetii specific genes. C. burnetii was not detected although a Coxiella-like endosymbiont was identified that was closely related to Coxiella symbionts from Ornithodoros capensis ticks. Because domestic and wild ungulates are the primary source of C. burnetii, we suggest that the presence of free-ranging, native and non-native ungulates in some off-shore islands in the Arabian Gulf could disseminate C. burnetii to seabirds. More comprehensive studies on seabird colonies are needed to better understand the diversity and prevalence of Coxiella symbionts and to establish if C. burnetii is endemic on some of these islands.

Pan-genomic analysis to redefine species and subspecies based on quantum discontinuous variation: the Klebsiella paradigm

Background

Various methods are currently used to define species and are based on the phylogenetic marker 16S ribosomal RNA gene sequence, DNA-DNA hybridization and DNA GC content. However, these are restricted genetic tools and showed significant limitations.

Results

In this work, we describe an alternative method to build taxonomy by analyzing the pan-genome composition of different species of the Klebsiella genus. Klebsiella species are Gram-negative bacilli belonging to the large Enterobacteriaceae family. Interestingly, when comparing the core/pan-genome ratio; we found a clear discontinuous variation that can define a new species.

Conclusions

Using this pan-genomic approach, we showed that Klebsiella pneumoniae subsp. ozaenae and Klebsiella pneumoniae subsp. rhinoscleromatis are species of the Klebsiella genus, rather than subspecies of Klebsiella pneumoniae. This pan-genomic analysis, helped to develop a new tool for defining species introducing a quantic perspective for taxonomy.

Reviewers

This article was reviewed by William Martin, Pierre Pontarotti and Pere Puigbo (nominated by Dr Yuri Wolf).

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-015-0085-2) contains supplementary material, which is available to authorized users.

Using Whole Genome Analysis to Examine Recombination across Diverse Sequence Types of Staphylococcus aureus

Staphylococcus aureus is an important clinical pathogen worldwide and understanding this organism's phylogeny and, in particular, the role of recombination, is important both to understand the overall spread of virulent lineages and to characterize outbreaks. To further elucidate the phylogeny of S. aureus, 35 diverse strains were sequenced using whole genome sequencing. In addition, 29 publicly available whole genome sequences were included to create a single nucleotide polymorphism (SNP)-based phylogenetic tree encompassing 11 distinct lineages. All strains of a particular sequence type fell into the same clade with clear groupings of the major clonal complexes of CC8, CC5, CC30, CC45 and CC1. Using a novel analysis method, we plotted the homoplasy density and SNP density across the whole genome and found evidence of recombination throughout the entire chromosome, but when we examined individual clonal lineages we found very little recombination. However, when we analyzed three branches of multiple lineages, we saw intermediate and differing levels of recombination between them. These data demonstrate that in S. aureus, recombination occurs across major lineages that subsequently expand in a clonal manner. Estimated mutation rates for the CC8 and CC5 lineages were different from each other. While the CC8 lineage rate was similar to previous studies, the CC5 lineage was 100-fold greater. Fifty known virulence genes were screened in all genomes in silico to determine their distribution across major clades. Thirty-three genes were present variably across clades, most of which were not constrained by ancestry, indicating horizontal gene transfer or gene loss.

The Recent Evolution of a Maternally-Inherited Endosymbiont of Ticks Led to the Emergence of the Q Fever Pathogen, Coxiella burnetii

Q fever is a highly infectious disease with a worldwide distribution. Its causative agent, the intracellular bacterium Coxiella burnetii, infects a variety of vertebrate species, including humans. Its evolutionary origin remains almost entirely unknown and uncertainty persists regarding the identity and lifestyle of its ancestors. A few tick species were recently found to harbor maternally-inherited Coxiella-like organisms engaged in symbiotic interactions, but their relationships to the Q fever pathogen remain unclear. Here, we extensively sampled ticks, identifying new and atypical Coxiella strains from 40 of 58 examined species, and used this data to infer the evolutionary processes leading to the emergence of C. burnetii. Phylogenetic analyses of multi-locus typing and whole-genome sequencing data revealed that Coxiella-like organisms represent an ancient and monophyletic group allied to ticks. Remarkably, all known C. burnetii strains originate within this group and are the descendants of a Coxiella-like progenitor hosted by ticks. Using both colony-reared and field-collected gravid females, we further establish the presence of highly efficient maternal transmission of these Coxiella-like organisms in four examined tick species, a pattern coherent with an endosymbiotic lifestyle. Our laboratory culture assays also showed that these Coxiella-like organisms were not amenable to culture in the vertebrate cell environment, suggesting different metabolic requirements compared to C. burnetii. Altogether, this corpus of data demonstrates that C. burnetii recently evolved from an inherited symbiont of ticks which succeeded in infecting vertebrate cells, likely by the acquisition of novel virulence factors.

How virulent infectious diseases emerge from non-pathogenic organisms is a challenging question. Here, we address this evolutionary issue in the case of Q fever. Its causative agent, the intracellular bacterium Coxiella burnetii, is extremely infectious to humans and a variety of animals. However, uncertainty persists regarding its evolutionary origin, including the identity and lifestyle of its ancestors. In this article, we show that C. burnetii arose from a rare evolutionary transformation of a maternally-inherited endosymbiont of ticks into a specialized and virulent pathogen of vertebrates. While arthropod symbionts are typically transmitted maternally and thought not to be infectious to vertebrates, we establish here that one Coxiella symbiont has evolved the necessary adaptations to exploit the vertebrate cell, leading to the emergence of Q fever.

Genomic epidemiology of the Haitian cholera outbreak: a single introduction followed by rapid, extensive, and continued spread characterized the onset of the epidemic

For centuries, cholera has been one of the most feared diseases. The causative agent Vibrio cholerae is a waterborne Gram-negative enteric pathogen eliciting a severe watery diarrheal disease. In October 2010, the seventh pandemic reached Haiti, a country that had not experienced cholera for more than a century. By using whole-genome sequence typing and mapping strategies of 116 serotype O1 strains from global sources, including 44 Haitian genomes, we present a detailed reconstructed evolutionary history of the seventh pandemic with a focus on the Haitian outbreak. We catalogued subtle genomic alterations at the nucleotide level in the genome core and architectural rearrangements from whole-genome map comparisons. Isolates closely related to the Haitian isolates caused several recent outbreaks in southern Asia. This study provides evidence for a single-source introduction of cholera from Nepal into Haiti followed by rapid, extensive, and continued clonal expansion. The phylogeographic patterns in both southern Asia and Haiti argue for the rapid dissemination of V. cholerae across the landscape necessitating real-time surveillance efforts to complement the whole-genome epidemiological analysis. As eradication efforts move forward, phylogeographic knowledge will be important for identifying persistent sources and monitoring success at regional levels. The results of molecular and epidemiological analyses of this outbreak suggest that an indigenous Haitian source of V. cholerae is unlikely and that an indigenous source has not contributed to the genomic evolution of this clade.

In this genomic epidemiology study, we have applied high-resolution whole-genome-based sequence typing methodologies on a comprehensive set of genome sequences that have become available in the aftermath of the Haitian cholera epidemic. These sequence resources enabled us to reassess the degree of genomic heterogeneity within the Vibrio cholerae O1 serotype and to refine boundaries and evolutionary relationships. The established phylogenomic framework showed how outbreak isolates fit into the global phylogeographic patterns compared to a comprehensive globally and temporally diverse strain collection and provides strong molecular evidence that points to a nonindigenous source of the 2010 Haitian cholera outbreak and refines epidemiological standards used in outbreak investigations for outbreak inclusion/exclusion following the concept of genomic epidemiology. The generated phylogenomic data have major public health relevance in translating sequence-based information to assist in future diagnostic, epidemiological, surveillance, and forensic studies of cholera.

The genome of Coxiella burnetii Z3055, a clone linked to the Netherlands Q fever outbreaks, provides evidence for the role of drift in the emergence of epidemic clones

Coxiella burnetii is a pathogen causing Q fever. The aim of our work was to study Z3055, a strain that is genotypically related to the strain causing the Netherlands outbreak. We compared Z3055 to 5 other completed genomes available in GenBank. We calculated the blast score ratio (BSR) to analyze genetic differences among the strains. The ratio core genome/pangenome was 98% likely other bacteria with closed pangenomes. Differences between Z3055 and the reference NMI consisted only of point mutations and insertion/deletion (INDELs). Non-synonymous mutations significantly increased in genes coding for membrane proteins (16/156 vs 103/1757, bilateral Chi(2) test, p<0.05), ankyrin repeat domains containing proteins (2/9 vs 117/1904, bilateral Chi(2) test, p<0.05), transcription factors (7/53 vs 112/1860, bilateral Chi(2) test, p<0.05) and translation proteins (15/144 vs 109/1655, bilateral Chi(2) test, p<0.05). The evolution of this strain may have been driven by mutations in critical genes.

Phylogeography of Bacillus anthracis in the country of Georgia shows evidence of population structuring and is dissimilar to other regional genotypes

Sequence analyses and subtyping of Bacillus anthracis strains from Georgia reveal a single distinct lineage (Aust94) that is ecologically established. Phylogeographic analysis and comparisons to a global collection reveals a clade that is mostly restricted to Georgia. Within this clade, many groups are found around the country, however at least one subclade is only found in the eastern part. This pattern suggests that dispersal into and out of Georgia has been rare and despite historical dispersion within the country, for at least for one lineage, current spread is limited.

Genotyping of Coxiella burnetii from domestic ruminants and human in Hungary: indication of various genotypes

BACKGROUND

Information about the genotypic characteristic of Coxiella burnetii from Hungary is lacking. The aim of this study is to describe the genetic diversity of C. burnetii in Hungary and compare genotypes with those found elsewhere. A total of 12 samples: (cattle, n = 6, sheep, n = 5 and human, n = 1) collected from across Hungary were studied by a 10-loci multispacer sequence typing (MST) and 6-loci multiple-locus variable-number of tandem repeat analysis (MLVA). Phylogenetic relationships among MST genotypes show how these Hungarian samples are related to others collected around the world.

RESULTS

Three MST genotypes were identified: sequence type (ST) 20 has also been identified in ruminants from other European countries and the USA, ST28 was previously identified in Kazakhstan, and the proposed ST37 is novel. All MST genotypes yielded different MLVA genotypes and three different MLVA genotypes were identified within ST20 samples alone. Two novel MLVA types 0-9-5-5-6-2 (AG) and 0-8-4-5-6-2 (AF) (Ms23-Ms24-Ms27-Ms28-Ms33-Ms34) were defined in the ovine materials correlated with ST28 and ST37. Samples from different parts of the phylogenetic tree were associated with different hosts, suggesting host-specific adaptations.

CONCLUSIONS

Even with the limited number of samples analysed, this study revealed high genetic diversity among C. burnetii in Hungary. Understanding the background genetic diversity will be essential in identifying and controlling outbreaks.

High prevalence and two dominant host-specific genotypes of Coxiella burnetii in U.S. milk

Background

Coxiella burnetii causes Q fever in humans and Coxiellosis in animals; symptoms range from general malaise to fever, pneumonia, endocarditis and death. Livestock are a significant source of human infection as they shed C. burnetii cells in birth tissues, milk, urine and feces. Although prevalence of C. burnetii is high, few Q fever cases are reported in the U.S. and we have a limited understanding of their connectedness due to difficulties in genotyping. Here, we develop canonical SNP genotyping assays to evaluate spatial and temporal relationships among C. burnetii environmental samples and compare them across studies. Given the genotypic diversity of historical collections, we hypothesized that the current enzootic of Coxiellosis is caused by multiple circulating genotypes. We collected A) 23 milk samples from a single bovine herd, B) 134 commercial bovine and caprine milk samples from across the U.S., and C) 400 bovine and caprine samples from six milk processing plants over three years.

Results

We detected C. burnetii DNA in 96% of samples with no variance over time. We genotyped 88.5% of positive samples; bovine milk contained only a single genotype (ST20) and caprine milk was dominated by a second type (mostly ST8).

Conclusions

The high prevalence and lack of genotypic diversity is consistent with a model of rapid spread and persistence. The segregation of genotypes between host species is indicative of species-specific adaptations or dissemination barriers and may offer insights into the relative lack of human cases and characterizing genotypes.

Eight new genomes and synthetic controls increase the accessibility of rapid melt-MAMA SNP typing of Coxiella burnetii

The case rate of Q fever in Europe has increased dramatically in recent years, mainly because of an epidemic in the Netherlands in 2009. Consequently, there is a need for more extensive genetic characterization of the disease agent Coxiella burnetii in order to better understand the epidemiology and spread of this disease. Genome reference data are essential for this purpose, but only thirteen genome sequences are currently available. Current methods for typing C. burnetii are criticized for having problems in comparing results across laboratories, require the use of genomic control DNA, and/or rely on markers in highly variable regions. We developed in this work a method for single nucleotide polymorphism (SNP) typing of C. burnetii isolates and tissue samples based on new assays targeting ten phylogenetically stable synonymous canonical SNPs (canSNPs). These canSNPs represent previously known phylogenetic branches and were here identified from sequence comparisons of twenty-one C. burnetii genomes, eight of which were sequenced in this work. Importantly, synthetic control templates were developed, to make the method useful to laboratories lacking genomic control DNA. An analysis of twenty-one C. burnetii genomes confirmed that the species exhibits high sequence identity. Most of its SNPs (7,493/7,559 shared by >1 genome) follow a clonal inheritance pattern and are therefore stable phylogenetic typing markers. The assays were validated using twenty-six genetically diverse C. burnetii isolates and three tissue samples from small ruminants infected during the epidemic in the Netherlands. Each sample was assigned to a clade. Synthetic controls (vector and PCR amplified) gave identical results compared to the corresponding genomic controls and are viable alternatives to genomic DNA. The results from the described method indicate that it could be useful for cheap and rapid disease source tracking at non-specialized laboratories, which requires accurate genotyping, assay accessibility and inter-laboratory comparisons.