Comparative genomics and phylogenomics of hemotrophic mycoplasmas

An Update of Bovine Hemoplasmas Based on Phylogenetic and Genomics Analysis

Mycoplasma wenyonii and ‘Candidatus Mycoplasma haemobos’ are bacteria that have been described as significant hemoplasmas that infect cattle worldwide. Currently, three bovine hemoplasma genomes are known. This work aimed to describe the main genomic characteristics and the evolutionary relationships between hemoplasmas, and provide a list of epitopes predicted by immunoinformatics as diagnostic candidates for bovine hemoplasmosis. Thus far, there is no vaccine to prevent this disease that economically impacts cattle production worldwide. Additionally, there is a lack of vaccines against bovine hemoplasmosis. In this work, we performed a genomic characterization of hemoplasmas, including two Mexican strains reported in bovines in the last few years. The generated information is a new scenario about the phylogeny of hemoplasmas. Also, we show genomic features among hemoplasmas that strengthen their characteristic genome plasticity of intracellular lifestyles. Finally, the elucidation of antigenic proteins in Mexican strains represents an opportunity to develop molecular detection methods and diagnoses.

Occurrence of 'Candidatus Mycoplasma haemosuis' in fattening pigs, sows and piglets in Germany using a novel gap-based quantitative real-time PCR assay

Background

The appearance of the novel porcine haemotrophic mycoplasma (HM) species ‘Candidatus Mycoplasma haemosuis’ was reported in apparently healthy but also in clinically sick animals in China, Korea and in a case report from Germany. Outside of Asia, however, nothing further is known about the frequency of ‘Ca. M. haemosuis’ in pigs to date. To investigate the distribution of this novel HM species in Germany, fattening pigs, sows and pre-suckling piglets were examined using a herein developed quantitative real-time PCR assay (qPCR). Because the piglets were sampled before the first colostrum uptake, additional information on a possible vertical transmission from dams to their offspring was obtained.

Results

Our novel qPCR assay successfully detected ‘Ca. M. haemosuis’ in all blood samples from the ‘Ca. M. haemosuis’-infected pigs. No cross-reactivity was detected when DNA from non-target Mycoplasma spp. and other bacterial species representing 10⁵ bacteria/reaction were used as a template. The lower limit of detection of the qPCR was thus 10 gap gene copies per reaction and 2.5 x 10³ genome equivalents (GE) per mL blood.

‘Candidatus M. haemosuis’ was detected by this qPCR in blood samples from a total out of 6.25% sows (13/208), 4.50% pre-suckling piglets (28/622) and 17.50% fattening pigs (35/200). On farm level, 3 out of 21 piglet producing farms (14.28%) and 9 out of 20 fattening farms (45.00%) were positive for ‘Ca. M. haemosuis’. Co-infections with M. suis were evident in all age groups.

Conclusion

‘Candidatus M. haemosuis’ infection is present in German pig farms and the detection of the novel porcine HM species in piglets immediately after birth before colostrum intake indicates vertical transmission. The novel qPCR assay specific for ‘Ca. M. haemosuis’ described herein will be a prerequisite for future studies on the prevalence, epidemiology as well as the clinical and economic impact of ‘Ca. M. haemosuis’ infections.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12917-022-03147-1.

Evaluating the Genetic Capacity of Mycoplasmas for Coenzyme A Biosynthesis in a Search for New Anti-mycoplasma Targets

Mycoplasmas are responsible for a wide range of disease states in both humans and animals, in which their parasitic lifestyle has allowed them to reduce their genome sizes and curtail their biosynthetic capabilities. The subsequent dependence on their host offers a unique opportunity to explore pathways for obtaining and producing cofactors - such as coenzyme A (CoA) - as possible targets for the development of new anti-mycoplasma agents. CoA plays an essential role in energy and fatty acid metabolism and is required for membrane synthesis. However, our current lack of knowledge of the relevance and importance of the CoA biosynthesis pathway in mycoplasmas, and whether it could be bypassed within their pathogenic context, prevents further exploration of the potential of this pathway. In the universal, canonical CoA biosynthesis pathway, five enzymes are responsible for the production of CoA. Given the inconsistent presence of the genes that code for these enzymes across Mycoplasma genomes, this study set out to establish the genetic capacity of mycoplasmas to synthesize their own CoA de novo. Existing functional annotations and sequence, family, motif, and domain analysis of protein products were used to determine the existence of relevant genes in Mycoplasma genomes. We found that most Mycoplasma species do have the genetic capacity to synthesize CoA, but there was a differentiated prevalence of these genes across species. Phylogenetic analysis indicated that the phylogenetic position of a species could not be used to predict its enzyme-encoding gene combinations. Despite this, the final enzyme in the biosynthesis pathway - dephospho-coenzyme A kinase (DPCK) - was found to be the most common among the studied species, suggesting that it has the most potential as a target in the search for new broad-spectrum anti-mycoplasma agents.

Hemotropic mycoplasmas (hemoplasmas) in wild boars, hunting dogs, and hunters from two Brazilian regions

Haemotropic mycoplasmas (haemoplasmas) are small pleomorphic bacteria infecting erythrocytes of several mammalian species, including human beings. No study to date has focused on the risk of bacteria exposure in hunting activities, particularly in natural environments of highly tick-infested areas. Accordingly, the present study aimed to assess haemoplasma occurrence in the complex encompassing wild boars, hunting dogs and hunters of Brazil. A total of 38/65 (58.5%) wild boars and 94/159 (59.1%) dogs were positive by qPCR for at least one haemoplasma. All 25 hunters were negative. Dogs with high hunting frequency were 2.4 more likely to be infected. Sequencing revealed a probable novel haemoplasma species in wild boars. Although exposure to haemoplasma species was present, the study herein found no evidence of cross-species transmission.

Worldwide occurrence of haemoplasmas in wildlife: Insights into the patterns of infection, transmission, pathology and zoonotic potential

Haemotropic mycoplasmas (haemoplasmas) have increasingly attracted the attention of wildlife disease researchers due to a combination of wide host range, high prevalence and genetic diversity. A systematic review identified 75 articles that investigated haemoplasma infection in wildlife by molecular methods (chiefly targeting partial 16S rRNA gene sequences), which included 131 host genera across six orders. Studies were less common in the Eastern Hemisphere (especially Africa and Asia) and more frequent in the Artiodactyla and Carnivora. Meta-analysis showed that infection prevalence did not vary by geographic region nor host order, but wild hosts showed significantly higher prevalence than captive hosts. Using a taxonomically flexible machine learning algorithm, we also found vampire bats and cervids to have greater prevalence, whereas mink, a subclade of vesper bats, and true foxes all had lower prevalence compared to the remaining sampled mammal phylogeny. Haemoplasma genotype and nucleotide diversity varied little among wild mammals but were marginally lower in primates and bats. Coinfection with more than one haemoplasma species or genotype was always confirmed when assessed. Risk factors of infection identified were sociality, age, males and high trophic levels, and both prevalence and diversity were often higher in undisturbed environments. Haemoplasmas likely use different and concurrent transmission routes and typically display enzootic dynamics when wild populations are studied longitudinally. Haemoplasma pathology is poorly known in wildlife but appears subclinical. Candidatus Mycoplasma haematohominis, which causes disease in humans, probably has it natural host in bats. Haemoplasmas can serve as a model system in ecological and evolutionary studies, and future research on these pathogens in wildlife must focus on increasing the geographic range and taxa of studies and elucidating pathology, transmission and zoonotic potential. To facilitate such work, we recommend using universal PCR primers or NGS protocols to detect novel haemoplasmas and other genetic markers to differentiate among species and infer cross-species transmission.

Persistence of Functional Protein Domains in Mycoplasma Species and their Role in Host Specificity and Synthetic Minimal Life

Mycoplasmas are the smallest self-replicating organisms and obligate parasites of a specific vertebrate host. An in-depth analysis of the functional capabilities of mycoplasma species is fundamental to understand how some of simplest forms of life on Earth succeeded in subverting complex hosts with highly sophisticated immune systems. In this study we present a genome-scale comparison, focused on identification of functional protein domains, of 80 publically available mycoplasma genomes which were consistently re-annotated using a standardized annotation pipeline embedded in a semantic framework to keep track of the data provenance. We examined the pan- and core-domainome and studied predicted functional capability in relation to host specificity and phylogenetic distance. We show that the pan- and core-domainome of mycoplasma species is closed. A comparison with the proteome of the “minimal” synthetic bacterium JCVI-Syn3.0 allowed us to classify domains and proteins essential for minimal life. Many of those essential protein domains, essential Domains of Unknown Function (DUFs) and essential hypothetical proteins are not persistent across mycoplasma genomes suggesting that mycoplasma species support alternative domain configurations that bypass their essentiality. Based on the protein domain composition, we could separate mycoplasma species infecting blood and tissue. For selected genomes of tissue infecting mycoplasmas, we could also predict whether the host is ruminant, pig or human. Functionally closely related mycoplasma species, which have a highly similar protein domain repertoire, but different hosts could not be separated. This study provides a concise overview of the functional capabilities of mycoplasma species, which can be used as a basis to further understand host-pathogen interaction or to design synthetic minimal life.

Ureaplasma diversum Genome Provides New Insights about the Interaction of the Surface Molecules of This Bacterium with the Host

Whole genome sequencing and analyses of Ureaplasma diversum ATCC 49782 was undertaken as a step towards understanding U. diversum biology and pathogenicity. The complete genome showed 973,501 bp in a single circular chromosome, with 28.2% of G+C content. A total of 782 coding DNA sequences (CDSs), and 6 rRNA and 32 tRNA genes were predicted and annotated. The metabolic pathways are identical to other human ureaplasmas, including the production of ATP via hydrolysis of the urea. Genes related to pathogenicity, such as urease, phospholipase, hemolysin, and a Mycoplasma Ig binding protein (MIB)-Mycoplasma Ig protease (MIP) system were identified. More interestingly, a large number of genes (n = 40) encoding surface molecules were annotated in the genome (lipoproteins, multiple-banded antigen like protein, membrane nuclease lipoprotein and variable surface antigens lipoprotein). In addition, a gene encoding glycosyltransferase was also found. This enzyme has been associated with the production of capsule in mycoplasmas and ureaplasma. We then sought to detect the presence of a capsule in this organism. A polysaccharide capsule from 11 to 17 nm of U. diversum was observed trough electron microscopy and using specific dyes. This structure contained arabinose, xylose, mannose, galactose and glucose. In order to understand the inflammatory response against these surface molecules, we evaluated the response of murine macrophages J774 against viable and non-viable U. diversum. As with viable bacteria, non-viable bacteria were capable of promoting a significant inflammatory response by activation of Toll like receptor 2 (TLR2), indicating that surface molecules are important for the activation of inflammatory response. Furthermore, a cascade of genes related to the inflammasome pathway of macrophages was also up-regulated during infection with viable organisms when compared to non-infected cells. In conclusion, U. diversum has a typical ureaplasma genome and metabolism, and its surface molecules, including the identified capsular material, represent major components of the organism immunopathogenesis.

Functional Basis of Microorganism Classification

Correctly identifying nearest “neighbors” of a given microorganism is important in industrial and clinical applications where close relationships imply similar treatment. Microbial classification based on similarity of physiological and genetic organism traits (polyphasic similarity) is experimentally difficult and, arguably, subjective. Evolutionary relatedness, inferred from phylogenetic markers, facilitates classification but does not guarantee functional identity between members of the same taxon or lack of similarity between different taxa. Using over thirteen hundred sequenced bacterial genomes, we built a novel function-based microorganism classification scheme, functional-repertoire similarity-based organism network (FuSiON; flattened to fusion). Our scheme is phenetic, based on a network of quantitatively defined organism relationships across the known prokaryotic space. It correlates significantly with the current taxonomy, but the observed discrepancies reveal both (1) the inconsistency of functional diversity levels among different taxa and (2) an (unsurprising) bias towards prioritizing, for classification purposes, relatively minor traits of particular interest to humans. Our dynamic network-based organism classification is independent of the arbitrary pairwise organism similarity cut-offs traditionally applied to establish taxonomic identity. Instead, it reveals natural, functionally defined organism groupings and is thus robust in handling organism diversity. Additionally, fusion can use organism meta-data to highlight the specific environmental factors that drive microbial diversification. Our approach provides a complementary view to cladistic assignments and holds important clues for further exploration of microbial lifestyles. Fusion is a more practical fit for biomedical, industrial, and ecological applications, as many of these rely on understanding the functional capabilities of the microbes in their environment and are less concerned with phylogenetic descent.

Taxonomic classification of microorganisms according to similarity is important for industrial and clinical applications where close relationships imply similar uses and/or treatments. Current microbial taxonomy is phylogeny-guided, i.e., the organisms are grouped based on their evolutionary relationships, defined by vertical inheritance of genetic information from mother to daughter cells. Microbes, however, are capable of horizontal gene transfer (HGT). Thus, the current taxonomic assignments cannot guarantee genome-encoded molecular functional similarity; i.e. two microbes of the same taxonomic group inhabiting different environments may be very different—just as your cousin may be more different from you than your unrelated best friend. Our work establishes a computational framework for comparison of microorganisms based on their molecular functionality. In our functional-repertoire similarity-based organism network (FuSiON; flattened to fusion) representation, organisms can be consistently assigned to groups based on a quantitative measure of their functional similarities. Our approach highlights the specific environmental factor(s) that explain the functional differences between groups of microorganism. Fusion is a more practical choice for biomedical, industrial, and ecological applications, as many of these rely on understanding the functional capabilities of the microbes in their environment.

The identification of oppA gene homologues as part of the oligopeptide transport system in mycoplasmas

The lack of an annotated oppA gene as part of many oligopeptide permease (opp) operons has questioned the necessity of the oligopeptide-binding domain (OppA) as a part of the Opp transport system in mycoplasmas. This study investigated the occurrence of an oppA gene as part of the oppBCDF operon in 42 mycoplasma genomes. Except for hemoplasma, all mycoplasmas were found to possess one or more copies of the oppBCDF operon and with the help of similarity searches their oppA genes could be identified. Phylogenetic analysis of the combined OppABCDF amino acid sequences allowed them to be grouped into three types. Each type has a unique set of conserved motifs, which are likely to reflect substrate preference and adaption strategies. Our approach allowed the identification of oppA gene homologues for all mycoplasma opp operons and thereby provides a method for re-evaluating the current annotation of oppA genes in mycoplasma genomes.

Ends of the line for tmRNA-SmpB

Genes for the RNA tmRNA and protein SmpB, partners in the trans-translation process that rescues stalled ribosomes, have previously been found in all bacteria and some organelles. During a major update of The tmRNA Website (relocated to http://bioinformatics.sandia.gov/tmrna), including addition of an SmpB sequence database, we found some bacteria that lack functionally significant regions of SmpB. Three groups with reduced genomes have lost the central loop of SmpB, which is thought to improve alanylation and EF-Tu activation: Carsonella, Hodgkinia, and the hemoplasmas (hemotropic Mycoplasma). Carsonella has also lost the SmpB C-terminal tail, thought to stimulate the decoding center of the ribosome. We validate recent identification of tmRNA homologs in oomycete mitochondria by finding partner genes from oomycete nuclei that target SmpB to the mitochondrion. We have moreover identified through exhaustive search a small number of complete, but often highly derived, bacterial genomes that appear to lack a functional copy of either the tmRNA or SmpB gene (but not both). One Carsonella isolate exhibits complete degradation of the tmRNA gene sequence yet its smpB shows no evidence for relaxed selective constraint, relative to other genes in the genome. After loss of the SmpB central loop in the hemoplasmas, one subclade apparently lost tmRNA. Carsonella also exhibits gene overlap such that tmRNA maturation should produce a non-stop smpB mRNA. At least some of the tmRNA/SmpB-deficient strains appear to further lack the ArfA and ArfB backup systems for ribosome rescue. The most frequent neighbors of smpB are the tmRNA gene, a ratA/rnfH unit, and the gene for RNaseR, a known physical and functional partner of tmRNA-SmpB.

Microscopy and genomic analysis of Mycoplasma parvum strain Indiana

Mycoplasma parvum [Eperythrozoon parvum] is the second hemotrophic mycoplasma (hemoplasma) described in pigs. Unlike M. suis, its closest phylogenetic relative, M. parvum, is considered a non-pathogenic bacterium in this host species. Natural infection of a domestic, 6-month-old splenectomized pig with M. parvum strain Indiana is described herein. Light and scanning electron microscopy of the bacteria were performed in addition to whole genome sequencing, analysis, and comparison to the genome of M. suis strain Illinois. Neither clinical signs nor anemia were observed during the infection. Microscopy analyses revealed coccoid to rod- shaped organisms varying from 0.2 to 0.5 μm; they were observed individually or in short chains by both light and electron microscopy, however less than 30% of the red blood cells were infected at peak bacteremia. The single circular chromosome of M. parvum was only 564 395 bp, smaller than M. genitalium, previously considered the tiniest member of the Mollicutes. Its general genomic features were similar to others in this class and species circumscription was verified by phylogenomic analysis. A gene-by-gene comparison between M. suis and M. parvum revealed all protein coding sequences (CDS) with assigned functions were shared, including metabolic functions, transporters and putative virulence factors. However, the number of CDS in paralogous gene families was remarkably different with about half as many paralogs in M. parvum. The differences in paralogous genes may be implicated in the different pathogenic potential of these two species, however variable gene expression may also play a role. Both are areas of ongoing investigation.