Comparative genomics of Mycoplasma: analysis of conserved essential genes and diversity of the pan-genome

Whole-Genome Sequencing and Phenotypic Analysis of Streptococcus equi subsp. zooepidemicus Sequence Type 147 Isolated from China

Streptococcus equi subsp. zooepidemicus (S. zooepidemicus) is one of the important zoonotic and opportunistic pathogens. In recent years, there has been growing evidence that supports the potential role of S. zooepidemicus in severe diseases in horses and other animals, including humans. Furthermore, the clinical isolation and drug resistance rates of S. zooepidemicus have been increasing yearly, leading to interest in its in-depth genomic analysis. In order to deepen the understanding of the S. zooepidemicus characteristics and genomic features, we investigated the genomic islands, mobile genetic elements, virulence and resistance genes, and phenotype of S. zooepidemicus strain ZHZ 211 (ST147), isolated from an equine farm in China. We obtained a 2.18 Mb, high-quality chromosome and found eight genomic islands. According to a comparative genomic investigation with other reference strains, ZHZ 211 has more virulence factors, like an iron uptake system, adherence, exoenzymes, and antiphagocytosis. More interestingly, ZHZ 211 has acquired a mobile genetic element (MGE), prophage Ph01, which was found to be in the chromosome of this strain and included two hyaluronidase (hyl) genes, important virulence factors of the strain. Moreover, two transposons and two virulence (virD4) genes were found to be located in the same genome island of ZHZ 211. In vitro phenotypic results showed that ZHZ 211 grows faster and is resistant to clarithromycin, enrofloxacin, and sulfonamides. The higher biofilm-forming capabilities of ZHZ 211 may provide a competitive advantage for survival in its niche. The results expand our understanding of the genomic, pathogenicity, and resistance characterization of Streptococcus zooepidemicus and facilitate further exploration of its molecular pathogenic mechanism.

Recent advances in the characterization of essential genes and development of a database of essential genes

Over the past few decades, there has been a significant interest in the study of essential genes, which are crucial for the survival of an organism under specific environmental conditions and thus have practical applications in the fields of synthetic biology and medicine. An increasing amount of experimental data on essential genes has been obtained with the continuous development of technological methods. Meanwhile, various computational prediction methods, related databases and web servers have emerged accordingly. To facilitate the study of essential genes, we have established a database of essential genes (DEG), which has become popular with continuous updates to facilitate essential gene feature analysis and prediction, drug and vaccine development, as well as artificial genome design and construction. In this article, we summarized the studies of essential genes, overviewed the relevant databases, and discussed their practical applications. Furthermore, we provided an overview of the main applications of DEG and conducted comprehensive analyses based on its latest version. However, it should be noted that the essential gene is a dynamic concept instead of a binary one, which presents both opportunities and challenges for their future development.

Mycoplasmas as Host Pantropic and Specific Pathogens: Clinical Implications, Gene Transfer, Virulence Factors, and Future Perspectives

Mycoplasmas as economically important and pantropic pathogens can cause similar clinical diseases in different hosts by eluding host defense and establishing their niches despite their limited metabolic capacities. Besides, enormous undiscovered virulence has a fundamental role in the pathogenesis of pathogenic mycoplasmas. On the other hand, they are host-specific pathogens with some highly pathogenic members that can colonize a vast number of habitats. Reshuffling mycoplasmas genetic information and evolving rapidly is a way to avoid their host's immune system. However, currently, only a few control measures exist against some mycoplasmosis which are far from satisfaction. This review aimed to provide an updated insight into the state of mycoplasmas as pathogens by summarizing and analyzing the comprehensive progress, current challenge, and future perspectives of mycoplasmas. It covers clinical implications of mycoplasmas in humans and domestic and wild animals, virulence-related factors, the process of gene transfer and its crucial prospects, the current application and future perspectives of nanotechnology for diagnosing and curing mycoplasmosis, Mycoplasma vaccination, and protective immunity. Several questions remain unanswered and are recommended to pay close attention to. The findings would be helpful to develop new strategies for basic and applied research on mycoplasmas and facilitate the control of mycoplasmosis for humans and various species of animals.

Metal utilization in genome-reduced bacteria: Do human mycoplasmas rely on iron?

Mycoplasmas are parasitic bacteria with streamlined genomes and complex nutritional requirements. Although iron is vital for almost all organisms, its utilization by mycoplasmas is controversial. Despite its minimalist nature, mycoplasmas can survive and persist within the host, where iron availability is rigorously restricted through nutritional immunity. In this review, we describe the putative iron-enzymes, transporters, and metalloregulators of four relevant human mycoplasmas. This work brings in light critical differences in the mycoplasma-iron interplay. Mycoplasma penetrans, the species with the largest genome (1.36 Mb), shows a more classic repertoire of iron-related proteins, including different enzymes using iron-sulfur clusters as well as iron storage and transport systems. In contrast, the iron requirement is less apparent in the three species with markedly reduced genomes, Mycoplasma genitalium (0.58 Mb), Mycoplasma hominis (0.67 Mb) and Mycoplasma pneumoniae (0.82 Mb), as they exhibit only a few proteins possibly involved in iron homeostasis. The multiple facets of iron metabolism in mycoplasmas illustrate the remarkable evolutive potential of these minimal organisms when facing nutritional immunity and question the dependence of several human-infecting species for iron. Collectively, our data contribute to better understand the unique biology and infective strategies of these successful pathogens.

A pan-genome method to determine core regions of the Bacillus subtilis and Escherichia coli genomes

Background: Synthetic engineering of bacteria to produce industrial products is a burgeoning field of research and application. In order to optimize genome design, designers need to understand which genes are essential, which are optimal for growth, and locations in the genome that will be tolerated by the organism when inserting engineered cassettes.

Methods: We present a pan-genome based method for the identification of core regions in a genome that are strongly conserved at the species level.

Results: We show that the core regions determined by our method contain all or almost all essential genes. This demonstrates the accuracy of our method as essential genes should be core genes. We show that we outperform previous methods by this measure. We also explain why there are exceptions to this rule for our method.

Conclusions: We assert that synthetic engineers should avoid deleting or inserting into these core regions unless they understand and are manipulating the function of the genes in that region. Similarly, if the designer wishes to streamline the genome, non-core regions and in particular low penetrance genes would be good targets for deletion. Care should be taken to remove entire cassettes with similar penetrance of the genes within cassettes as they may harbor toxin/antitoxin genes which need to be removed in tandem. The bioinformatic approach introduced here saves considerable time and effort relative to knockout studies on single isolates of a given species and captures a broad understanding of the conservation of genes that are core to a species.

Attenuation of Mycoplasma hyopneumoniae Strain ES-2 and Comparative Genomic Analysis of ES-2 and Its Attenuated Form ES-2L

Mycoplasma hyopneumoniae causes swine respiratory disease worldwide. Due to the difficulty of isolating and cultivating M. hyopneumoniae, very few attenuated strains have been successfully isolated, which hampers the development of attenuated vaccines. In order to produce an attenuated M. hyopneumoniae strain, we used the highly virulent M. hyopneumoniae strain ES-2, which was serially passaged in vitro 200 times to produce the attenuated strain ES-2L, and its virulence was evidenced to be low in an animal experiment. In order to elucidate the mechanisms underlying virulence attenuation, we performed whole-genome sequencing of both strains and conducted comparative genomic analyses of strain ES-2 and its attenuated form ES-2L. Strain ES-2L showed three large fragment deletion regions including a total of 18 deleted genes, compared with strain ES-2. Analysis of single-nucleotide polymorphisms (SNPs) and indels indicated that 22 dels were located in 19 predicted coding sequences. In addition to these indels, 348 single-nucleotide variations (SNVs) were identified between strains ES-2L and ES-2. These SNVs mapped to 99 genes where they appeared to induce amino acid substitutions and translation stops. The deleted genes and SNVs may be associated with decreased virulence of strain ES-2L. Our work provides a foundation for further examining virulence factors of M. hyopneumoniae and for the development of attenuated vaccines.

Development of New PCR Assay with SYBR Green I for Detection of Mycoplasma, Acholeplasma, and Ureaplasma sp. in Cell Cultures

Mycoplasma, Acholeplasma, and Ureaplasma sp. are atypical bacteria responsible for in vitro cell culture contaminations that can warp the results. These bacteria also cause human and animal infections and may lead to chronic diseases. In developed polymerase chain reaction (PCR) in this study a quantitative PCR with SYBR Green I fluorochrome was applied to facilitate the Mycoplasma, Acholeplasma, and Ureaplasma sp. DNA detection and identification. Screening Test-1 v.1 (triplex qPCR) allowed for the detection of 11 species. Test-1 v.2 (three single qPCRs) pre-identified three subgroups, allowing for the reduction of using single qPCRs in Test-2 for species identification. The range of both tests was consistent with pharmacopeial requirements for microbial quality control of mammal cells and included detection of M. arginini, M. orale, M. hyorhinis, M. fermentans, M. genitalium, M. hominis, M. pneumoniae, M. salivarium, M. pirum, A. laidlawii, and U. urealyticum. Limit of detection values varied between 125–300 and 50–100 number of copies per milliliter in Test-1 and Test-2, respectively. Test-1 and Test-2 showed fully concordant results, allowed for time-saving detection and/or identification of selected species from Mycoplasma, Acholeplasma, and Ureaplasma in tested cell cultures.

A pan-genome method to determine core regions of the Bacillus subtilis and Escherichia coli genomes

Background: Synthetic engineering of bacteria to produce industrial products is a burgeoning field of research and application. In order to optimize genome design, designers need to understand which genes are essential, which are optimal for growth, and locations in the genome that will be tolerated by the organism when inserting engineered cassettes. Methods: We present a pan-genome based method for the identification of core regions in a genome that are strongly conserved at the species level. Results: We show that the core regions determined by our method contain all or almost all essential genes. This demonstrates the accuracy of our method as essential genes should be core genes. We show that we outperform previous methods by this measure. We also explain why there are exceptions to this rule for our method. Conclusions: We assert that synthetic engineers should avoid deleting or inserting into these core regions unless they understand and are manipulating the function of the genes in that region. Similarly, if the designer wishes to streamline the genome, non-core regions and in particular low penetrance genes would be good targets for deletion. Care should be taken to remove entire cassettes with similar penetrance of the genes within cassettes as they may harbor toxin/antitoxin genes which need to be removed in tandem. The bioinformatic approach introduced here saves considerable time and effort relative to knockout studies on single isolates of a given species and captures a broad understanding of the conservation of genes that are core to a species.

Essential gene prediction using limited gene essentiality information-An integrative semi-supervised machine learning strategy

Essential gene prediction helps to find minimal genes indispensable for the survival of any organism. Machine learning (ML) algorithms have been useful for the prediction of gene essentiality. However, currently available ML pipelines perform poorly for organisms with limited experimental data. The objective is the development of a new ML pipeline to help in the annotation of essential genes of less explored disease-causing organisms for which minimal experimental data is available. The proposed strategy combines unsupervised feature selection technique, dimension reduction using the Kamada-Kawai algorithm, and semi-supervised ML algorithm employing Laplacian Support Vector Machine (LapSVM) for prediction of essential and non-essential genes from genome-scale metabolic networks using very limited labeled dataset. A novel scoring technique, Semi-Supervised Model Selection Score, equivalent to area under the ROC curve (auROC), has been proposed for the selection of the best model when supervised performance metrics calculation is difficult due to lack of data. The unsupervised feature selection followed by dimension reduction helped to observe a distinct circular pattern in the clustering of essential and non-essential genes. LapSVM then created a curve that dissected this circle for the classification and prediction of essential genes with high accuracy (auROC > 0.85) even with 1% labeled data for model training. After successful validation of this ML pipeline on both Eukaryotes and Prokaryotes that show high accuracy even when the labeled dataset is very limited, this strategy is used for the prediction of essential genes of organisms with inadequate experimentally known data, such as Leishmania sp. Using a graph-based semi-supervised machine learning scheme, a novel integrative approach has been proposed for essential gene prediction that shows universality in application to both Prokaryotes and Eukaryotes with limited labeled data. The essential genes predicted using the pipeline provide an important lead for the prediction of gene essentiality and identification of novel therapeutic targets for antibiotic and vaccine development against disease-causing parasites.

Challenges in gene-oriented approaches for pangenome content discovery

Given a group of genomes, represented as the sets of genes that belong to them, the discovery of the pangenomic content is based on the search of genetic homology among the genes for clustering them into families. Thus, pangenomic analyses investigate the membership of the families to the given genomes. This approach is referred to as the gene-oriented approach in contrast to other definitions of the problem that takes into account different genomic features. In the past years, several tools have been developed to discover and analyse pangenomic contents. Because of the hardness of the problem, each tool applies a different strategy for discovering the pangenomic content. This results in a differentiation of the performance of each tool that depends on the composition of the input genomes. This review reports the main analysis instruments provided by the current state of the art tools for the discovery of pangenomic contents. Moreover, unlike previous works, the presented study compares pangenomic tools from a methodological perspective, analysing the causes that lead a given methodology to outperform other tools. The analysis is performed by taking into account different bacterial populations, which are synthetically generated by changing evolutionary parameters. The benchmarks used to compare the pangenomic tools, in addition to the computational pipeline developed for this purpose, are available at https://github.com/InfOmics/pangenes-review. Contact: V. Bonnici, R. Giugno Supplementary information: Supplementary data are available at Briefings in Bioinformatics online.

It's a small, small world: unravelling the role and evolution of small RNAs in organelle and endosymbiont genomes

Organelles and host-restricted bacterial symbionts are characterized by having highly reduced genomes that lack many key regulatory genes and elements. Thus, it has been hypothesized that the eukaryotic nuclear genome is primarily responsible for regulating these symbioses. However, with the discovery of organelle- and symbiont-expressed small RNAs (sRNAs) there is emerging evidence that these sRNAs may play a role in gene regulation as well. Here, we compare the diversity of organelle and bacterial symbiont sRNAs recently identified using genome-enabled '-omic' technologies and discuss their potential role in gene regulation. We also discuss how the genome architecture of small genomes may influence the evolution of these sRNAs and their potential function. Additionally, these new studies suggest that some sRNAs are conserved within organelle and symbiont taxa and respond to changes in the environment and/or their hosts. In summary, these results suggest that organelle and symbiont sRNAs may play a role in gene regulation in addition to nuclear-encoded host mechanisms.

Large-Scale Analysis of the Mycoplasma bovis Genome Identified Non-essential, Adhesion- and Virulence-Related Genes

Mycoplasma bovis is an important pathogen of cattle causing bovine mycoplasmosis. Clinical manifestations are numerous, but pneumonia, mastitis, and arthritis cases are mainly reported. Currently, no efficient vaccine is available and antibiotic treatments are not always satisfactory. The design of new, efficient prophylactic and therapeutic approaches requires a better understanding of the molecular mechanisms responsible for M. bovis pathogenicity. Random transposon mutagenesis has been widely used in Mycoplasma species to identify potential gene functions. Such an approach can also be used to screen genomes and search for essential and non-essential genes for growth. Here, we generated a random transposon mutant library of M. bovis strain JF4278 containing approximately 4000 independent insertion sites. We then coupled high-throughput screening of this mutant library to transposon sequencing and bioinformatic analysis to identify M. bovis non-essential, adhesion- and virulence-related genes. Three hundred and fifty-two genes of M. bovis were assigned as essential for growth in rich medium. Among the remaining non-essential genes, putative virulence-related factors were subsequently identified. The complete mutant library was screened for adhesion using primary bovine mammary gland epithelial cells. Data from this assay resulted in a list of conditional-essential genes with putative adhesion-related functions by identifying non-essential genes for growth that are essential for host cell-adhesion. By individually assessing the adhesion capacity of six selected mutants, two previously unknown factors and the adhesin TrmFO were associated with a reduced adhesion phenotype. Overall, our study (i) uncovers new, putative virulence-related genes; (ii) offers a list of putative adhesion-related factors; and (iii) provides valuable information for vaccine design and for exploring M. bovis biology, pathogenesis, and host-interaction.

Detection of Mycoplasma anatis, M. anseris, M. cloacale and Mycoplasma sp. 1220 in waterfowl using species-specific PCR assays

Mycoplasma anatis, M. anseris, M. cloacale and M. sp. 1220 colonise geese and ducks, and could be associated with infections of avian respiratory and nervous systems, cause mild to severe inflammation of cloaca and genital tracts, and embryo lethality. Co-occurrence of these Mycoplasma species in waterfowl is frequently detected and the identification of these mycoplasmas to the species level at a regular microbiology laboratory is difficult due to their similar morphological, cultural and biochemical properties. Moreover, species differentiation is only possible based on the sequence analysis of the product of a genus-specific PCR assay. Therefore, the aim of the current study was to develop an effective and robust method for the identification of these species in avian clinical specimens. Polymerase chain reaction (PCR) assays using species-specific primers, which target housekeeping genes in order to identify these species, were designed in the present study. The developed PCR assays can precisely identify these four mycoplasmas to the species level directly from DNA samples extracted from clinical specimens, and no cross-amplification was observed among these species and with other well-known avian mycoplasmas. The average sensitivity of the assays was 10¹−10² genomic equivalents per reaction. These conventional PCR assays can be run simultaneously at the same PCR cycling program, and the species can be differentiated directly (without sequence analysis) by gel electrophoresis due to the specific sizes of the amplicons. In conclusion, the presented species-specific assays were found to be suitable for routine use at regular veterinary diagnostic laboratories and promote the rapid, simple and cost-effective differentiation of these waterfowl Mycoplasma species.

Genome-driven cell engineering review: in vivo and in silico metabolic and genome engineering

Producing 'designer cells' with specific functions is potentially feasible in the near future. Recent developments, including whole-cell models, genome design algorithms and gene editing tools, have advanced the possibility of combining biological research and mathematical modelling to further understand and better design cellular processes. In this review, we will explore computational and experimental approaches used for metabolic and genome design. We will highlight the relevance of modelling in this process, and challenges associated with the generation of quantitative predictions about cell behaviour as a whole: although many cellular processes are well understood at the subsystem level, it has proved a hugely complex task to integrate separate components together to model and study an entire cell. We explore these developments, highlighting where computational design algorithms compensate for missing cellular information and underlining where computational models can complement and reduce lab experimentation. We will examine issues and illuminate the next steps for genome engineering.

Comparative genome-scale modelling of the pathogenic Flavobacteriaceae species Riemerella anatipestifer in China

Riemerella anatipestifer (RA) is a gram-negative bacterium that has a high potential to infect waterfowl. Although more and more genomes of RA have been generated comparaed to genomic analysis of RA still remains at the level of individual species. In this study, we analysed the pan-genome of 27 RA virulent isolates to reveal the intraspecies genomic diversity from various aspects. The multi-locus sequence typing (MLST) analysis suggests that the geographic origin of R. anatipestifer is Guangdong province, China. Results of pan-genome analysis revealed an open pan-genome for all 27 species with the sizes of 2967 genes. We identified 387 genes among 555 unique genes originated by horizontal gene transfer. Further studies showed 204 strain-specific HGT genes were predicted as virulent proteins. Screening the 1113 core genes in RA through subtractive genomic approach, 70 putative vaccine targets out of 125 non-cytoplasmic proteins have been predicted. Further analysis of these non A. platyrhynchos homologous proteins predicted that 56 essential proteins as drug target with more interaction partners were involved in unique metabolic pathways of RA. In conclusion, the present study indicated the essence and the diversity of RA and also provides useful information for identification of vaccine and drugs candidates in future.

An outbreak of Mycoplasma mycoides subspecies capri arthritis in young goats: a case study

Mycoplasmosis is a well-known cause of morbidity and mortality in small ruminants. Previously recognized outbreaks have involved arthritis, and pneumonia or pleuropneumonia. Modern bacteriology procedures rely less on isolation techniques that require special media for mollicutes given that these species are notoriously difficult to isolate, and rely more on PCR tests. We report an outbreak of arthritis, pleuropneumonia, and mild meningitis affecting dairy goat kids, spanning a period of 3 y, which had unusual epidemiologic characteristics related to husbandry practices. Lesions were characterized by polyarthritis of the appendicular joints, with copious joint fluid and extension of arthritic exudate beyond the joint itself. The cause remained unknown until serendipitous isolation of a mycoplasma on blood agar. Mycoplasmosis was not detected from synovial samples by a general mycoplasma PCR, despite multiple attempts. Isolated colonies were also negative by this general PCR assay. The isolate was identified as Mycoplasma mycoides subspecies capri, using universal 16S primers and amplicon sequencing. Testing of additional isolates from other diseased goats in the herd confirmed that this was the cause of illness. A failure to recognize the distinct nature of organisms of the M. mycoides group of mycoplasmas meant that a PCR test that cannot detect this group of organisms was utilized at first, and the etiology of the illness was overlooked for a period of time. Veterinary pathologists and microbiologists must be aware of the limitations of some PCR assays when confronted with joint disease and pleuropneumonia in small ruminants.

Genome Analyses of a New Mycoplasma Species from the Scorpion Centruroides vittatus

Arthropod Mycoplasma are little known endosymbionts in insects, primarily known as plant disease vectors. Mycoplasma in other arthropods such as arachnids are unknown. We report the first complete Mycoplasma genome sequenced, identified, and annotated from a scorpion, Centruroides vittatus, and designate it as Mycoplasma vittatus. We find the genome is at least a 683,827 bp single circular chromosome with a GC content of 42.7% and with 987 protein-coding genes. The putative virulence determinants include 11 genes associated with the virulence operon associated with protein synthesis or DNA transcription and ten genes with antibiotic and toxic compound resistance. Comparative analysis revealed that the M. vittatus genome is smaller than other Mycoplasma genomes and exhibits a higher GC content. Phylogenetic analysis shows M. vittatus as part of the Hominis group of Mycoplasma. As arthropod genomes accumulate, further novel Mycoplasma genomes may be identified and characterized.

Extracellular DNA release from the genome-reduced pathogen Mycoplasma hyopneumoniae is essential for biofilm formation on abiotic surfaces

Mycoplasma hyopneumoniae is an economically devastating, globally disseminated pathogen that can maintain a chronic infectious state within its host, swine. Here, we depict the events underpinning M. hyopneumoniae biofilm formation on an abiotic surface and demonstrate for the first time, biofilms forming on porcine epithelial cell monolayers and in the lungs of pigs, experimentally infected with M. hyopneumoniae. Nuclease treatment prevents biofilms forming on glass but not on porcine epithelial cells indicating that extracellular DNA (eDNA), which localises at the base of biofilms, is critical in the formation of these structures on abiotic surfaces. Subpopulations of M. hyopneumoniae cells, denoted by their ability to take up the dye TOTO-1 and release eDNA, were identified. A visually distinct sub-population of pleomorphic cells, that we refer to here as large cell variants (LCVs), rapidly transition from phase dark to translucent "ghost" cells. The translucent cells accumulate the membrane-impermeable dye TOTO-1, forming readily discernible membrane breaches immediately prior to lysis and the possible release of eDNA and other intracellular content (public goods) into the extracellular environment. Our novel observations expand knowledge of the lifestyles adopted by this wall-less, genome-reduced pathogen and provide further insights to its survival within farm environments and swine.

The Genomic Basis of Intrinsic and Acquired Antibiotic Resistance in the Genus Serratia

Serratia marcescens, a member of the Enterobacteriaceae family, was long thought to be a non-pathogenic bacterium prevalent in environmental habitats. Together with other members of this genus, it has emerged in recent years as an opportunistic nosocomial pathogen causing various types of infections. One important feature of pathogens belonging to this genus is their intrinsic and acquired resistance to a variety of antibiotic families, including β-lactam, aminoglycosides, quinolones and polypeptide antibiotics. The aim of this study was to elucidate which genes participate in the intrinsic and acquired antibiotic resistance of this genus in order to determine the Serratia genus resistome. We performed phylogenomic and comparative genomic analyses using 32 Serratia spp. genomes deposited in the NCBI GenBank from strains isolated from different ecological niches and different lifestyles. S. marcescens strain SmUNAM836, which was previously isolated from a Mexican adult with obstructive pulmonary disease, was included in this study. The results show that most of the antibiotic resistance genes (ARGs) were found on the chromosome, and to a lesser degree, on plasmids and transposons acquired through horizontal gene transfer. Four strains contained the gyrA point mutation in codon Ser83 that confers quinolone resistance. Pathogenic and environmental isolates presented a high number of ARGs, especially genes associated with efflux systems. Pathogenic strains, specifically nosocomial strains, presented more acquired resistance genes than environmental isolates. We may conclude that the environment provides a natural reservoir for antibiotic resistance, which has been underestimated in the medical field.

Minimal Cells-Real and Imagined

A minimal cell is one whose genome only encodes the minimal set of genes necessary for the cell to survive. Scientific reductionism postulates the best way to learn the first principles of cellular biology would be to use a minimal cell in which the functions of all genes and components are understood. The genes in a minimal cell are, by definition, essential. In 2016, synthesis of a genome comprised of only the set of essential and quasi-essential genes encoded by the bacterium Mycoplasma mycoides created a near-minimal bacterial cell. This organism performs the cellular functions common to all organisms. It replicates DNA, transcribes RNA, translates proteins, undergoes cell division, and little else. In this review, we examine this organism and contrast it with other bacteria that have been used as surrogates for a minimal cell.

A Comprehensive Overview of Online Resources to Identify and Predict Bacterial Essential Genes

Genes critical for the survival or reproduction of an organism in certain circumstances are classified as essential genes. Essential genes play a significant role in deciphering the survival mechanism of life. They may be greatly applied to pharmaceutics and synthetic biology. The continuous progress of experimental method for essential gene identification has accelerated the accumulation of gene essentiality data which facilitates the study of essential genes in silico. In this article, we present some available online resources related to gene essentiality, including bioinformatic software tools for transposon sequencing (Tn-seq) analysis, essential gene databases and online services to predict bacterial essential genes. We review several computational approaches that have been used to predict essential genes, and summarize the features used for gene essentiality prediction. In addition, we evaluate the available online bacterial essential gene prediction servers based on the experimentally validated essential gene sets of 30 bacteria from DEG. This article is intended to be a quick reference guide for the microbiologists interested in the essential genes.

Coevolution of the translational machinery optimizes initiation with unusual initiator tRNAs and initiation codons in mycoplasmas

Initiator tRNAs (i-tRNAs) are characterized by the presence of three consecutive GC base pairs (GC/GC/GC) in their anticodon stems in all domains of life. However, many mycoplasmas possess unconventional i-tRNAs wherein the highly conserved sequence of GC/GC/GC is represented by AU/GC/GC, GC/GC/GU or AU/GC/GU. These mycoplasmas also tend to preferentially utilize non-AUG initiation codons. To investigate if initiation with the unconventional i-tRNAs and non-AUG codons in mycoplasmas correlated with the changes in the other components of the translation machinery, we carried out multiple sequence alignments of genes encoding initiation factors (IF), 16S rRNAs, and the ribosomal proteins such as uS9, uS12 and uS13. In addition, the occurrence of Shine-Dalgarno sequences in mRNAs was analyzed. We observed that in the mycoplasmas harboring AU/GC/GU i-tRNAs, a highly conserved position of R131 in IF3, is represented by P, F or Y and, the conserved C-terminal tail (SKR) of uS9 is represented by the TKR sequence. Using the Escherichia coli model, we show that the change of R131 in IF3 optimizes initiation with the AU/GC/GU i-tRNAs. Also, the SKR to TKR change in uS9 was compatible with the R131P variation in IF3 for initiation with the AU/GC/GU i-tRNA variant. Interestingly, the mycoplasmas harboring AU/GC/GU i-tRNAs are also human pathogens. We propose that these mycoplasmas might have evolved a relaxed translational apparatus to adapt to the environment they encounter in the host.

The essential genomic landscape of the commensal Bifidobacterium breve UCC2003

Bifidobacteria are common gut commensals with purported health-promoting effects. This has encouraged scientific research into bifidobacteria, though recalcitrance to genetic manipulation and scarcity of molecular tools has hampered our knowledge on the precise molecular determinants of their health-promoting attributes and gut adaptation. To overcome this problem and facilitate functional genomic analyses in bifidobacteria, we created a large Tn5 transposon mutant library of the commensal Bifidobacterium breve UCC2003 that was further characterized by means of a Transposon Directed Insertion Sequencing (TraDIS) approach. Statistical analysis of transposon insertion distribution revealed a set of 453 genes that are essential for or markedly contribute to growth of this strain under laboratory conditions. These essential genes encode functions involved in the so-called bifid-shunt, most enzymes related to nucleotide biosynthesis and a range of housekeeping functions. Comparison to the Bifidobacterium and B. breve core genomes highlights a high degree of conservation of essential genes at the species and genus level, while comparison to essential gene datasets from other gut bacteria identified essential genes that appear specific to bifidobacteria. This work establishes a useful molecular tool for scientific discovery of bifidobacteria and identifies targets for further studies aimed at characterizing essential functions not previously examined in bifidobacteria.

Detecting the Diversity of Mycoplasma and Ureaplasma Endosymbionts Hosted by Trichomonas vaginalis Isolates

Objectives: The symbiosis of Trichomonas vaginalis and Mycoplasma hominis is the first described association between two obligate human parasites. Trichomonas is the niche and the vector for the transmission of M. hominis infection. This clinically significant symbiosis may affect T. vaginalis virulence and susceptibility to treatment. The aims of this study were to investigate the intracellularly present Mycoplasma and Ureaplasma species in T. vaginalis strains isolated from the vaginal discharge of infected women as well as to trace the diversity pattern among the species detected in the isolated strains. Methods: Hundred pure T. vaginalis cultures were isolated from ~7,500 patient specimens presented with clinical purulent vaginitis. PCR and sequencing for Mycoplasma/Ureaplasma spp. were performed in DNA extracted from the pure cultures. In addition, vaginal discharge samples were cultured for the presence of M. hominis and U. urealyticum. Phylogenetic analysis assisted the identification of interspecies relationships between the Mycoplasma and Ureaplasma isolates. Results: Fifty four percentage of T. vaginalis isolates were harboring Mycoplasma spp. Phylogenetic analysis revealed three distinct clusters, two with already characterized M. hominis and Ureaplasma spp. (37% of total Mycoplasma spp.), whereas one group formed a distinct cluster matched with the newly identified species Candidatus Mycoplasma girerdii (59.3%) and one or more unknown Mycoplasma spp. (3.7%). Conclusions:T. vaginalis strains associated with vaginal infection might host intracellular mycoplasmas or ureaplasmas. Intracellular Mollicutes that remain undetected in the extracellular environment when conventional diagnostic methods are implemented may comprise either novel species, such as Candidatus M. giredii, or unknown species with yet unexplored clinical significance.

Comprehensive bioinformatics analysis of Mycoplasma pneumoniae genomes to investigate underlying population structure and type-specific determinants

Mycoplasma pneumoniae is a significant cause of respiratory illness worldwide. Despite a minimal and highly conserved genome, genetic diversity within the species may impact disease. We performed whole genome sequencing (WGS) analysis of 107 M. pneumoniae isolates, including 67 newly sequenced using the Pacific BioSciences RS II and/or Illumina MiSeq sequencing platforms. Comparative genomic analysis of 107 genomes revealed >3,000 single nucleotide polymorphisms (SNPs) in total, including 520 type-specific SNPs. Population structure analysis supported the existence of six distinct subgroups, three within each type. We developed a predictive model to classify an isolate based on whole genome SNPs called against the reference genome into the identified subtypes, obviating the need for genome assembly. This study is the most comprehensive WGS analysis for M. pneumoniae to date, underscoring the power of combining complementary sequencing technologies to overcome difficult-to-sequence regions and highlighting potential differential genomic signatures in M. pneumoniae.

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.

Repetitive Elements in Mycoplasma hyopneumoniae Transcriptional Regulation

Transcriptional regulation, a multiple-step process, is still poorly understood in the important pig pathogen Mycoplasma hyopneumoniae. Basic motifs like promoters and terminators have already been described, but no other cis-regulatory elements have been found. DNA repeat sequences have been shown to be an interesting potential source of cis-regulatory elements. In this work, a genome-wide search for tandem and palindromic repetitive elements was performed in the intergenic regions of all coding sequences from M. hyopneumoniae strain 7448. Computational analysis demonstrated the presence of 144 tandem repeats and 1,171 palindromic elements. The DNA repeat sequences were distributed within the 5' upstream regions of 86% of transcriptional units of M. hyopneumoniae strain 7448. Comparative analysis between distinct repetitive sequences found in related mycoplasma genomes demonstrated different percentages of conservation among pathogenic and nonpathogenic strains. qPCR assays revealed differential expression among genes showing variable numbers of repetitive elements. In addition, repeats found in 206 genes already described to be differentially regulated under different culture conditions of M. hyopneumoniae strain 232 showed almost 80% conservation in relation to M. hyopneumoniae strain 7448 repeats. Altogether, these findings suggest a potential regulatory role of tandem and palindromic DNA repeats in the M. hyopneumoniae transcriptional profile.

The Evolution of Advanced Molecular Diagnostics for the Detection and Characterization of Mycoplasma pneumoniae

Over the past decade there have been significant advancements in the methods used for detecting and characterizing Mycoplasma pneumoniae, a common cause of respiratory illness and community-acquired pneumonia worldwide. The repertoire of available molecular diagnostics has greatly expanded from nucleic acid amplification techniques (NAATs) that encompass a variety of chemistries used for detection, to more sophisticated characterizing methods such as multi-locus variable-number tandem-repeat analysis (MLVA), Multi-locus sequence typing (MLST), matrix-assisted laser desorption ionization-time-of-flight mass spectrometry (MALDI-TOF MS), single nucleotide polymorphism typing, and numerous macrolide susceptibility profiling methods, among others. These many molecular-based approaches have been developed and employed to continually increase the level of discrimination and characterization in order to better understand the epidemiology and biology of M. pneumoniae. This review will summarize recent molecular techniques and procedures and lend perspective to how each has enhanced the current understanding of this organism and will emphasize how Next Generation Sequencing may serve as a resource for researchers to gain a more comprehensive understanding of the genomic complexities of this insidious pathogen.

The bacterial pangenome as a new tool for analysing pathogenic bacteria

The bacterial pangenome was introduced in 2005 and, in recent years, has been the subject of many studies. Thanks to progress in next-generation sequencing methods, the pangenome can be divided into two parts, the core (common to the studied strains) and the accessory genome, offering a large panel of uses. In this review, we have presented the analysis methods, the pangenome composition and its application as a study of lifestyle. We have also shown that the pangenome may be used as a new tool for redefining the pathogenic species. We applied this to the Escherichia coli and Shigella species, which have been a subject of controversy regarding their taxonomic and pathogenic position.

•

Pangenome is a new way of studying pathogenic bacteria.

•

Pangenome can be used as a taxonomic tool.

•

This review describes pangenome in the world of pathogenic bacteria.

Draft genome sequences of three Holospora species (Holospora obtusa, Holospora undulata, and Holospora elegans), endonuclear symbiotic bacteria of the ciliate Paramecium caudatum

We present draft genome sequences of three Holospora species, hosted by the ciliate Paramecium caudatum; that is, the macronucleus-specific H. obtusa and the micronucleus-specific H. undulata and H. elegans. We investigate functions of orthologous core genes conserved across the three Holospora species, which may be essential for the infection and survival in the host nucleus.

Mycoplasma Contamination of Cell Cultures: Vesicular Traffic in Bacteria and Control over Infectious Agents

Cell cultures are subject to contamination either with cells of other cultures or with microorganisms, including fungi, viruses, and bacteria. Mycoplasma contamination of cell cultures is of particular importance. Since cell cultures are used for the production of vaccines and physiologically active compounds, designing a system for controlling contaminants becomes topical for fundamental science and biotechnological production. The discovery of extracellular membrane vesicles in mycoplasmas makes it necessary to take into consideration the bacterial vesicular traffic in systems designed for controlling infectious agents. The extracellular vesicles of bacteria mediate the traffic of proteins and genes, participate in cell-to-cell interactions, as well as in the pathogenesis and development of resistance to antibiotics. The present review discusses the features of mycoplasmas, their extracellular vesicles, and the interaction between contaminants and eukaryotic cells. Furthermore, it provides an analysis of the problems associated with modern methods of diagnosis and eradication of mycoplasma contamination from cell cultures and prospects for their solution.

Genes found essential in other mycoplasmas are dispensable in Mycoplasma bovis

Mycoplasmas are regarded to be useful models for studying the minimum genetic complement required for independent survival of an organism. Mycoplasma bovis is a globally distributed pathogen causing pneumonia, mastitis, arthritis, otitis media and reproductive tract disease, and genome sequences of three strains, the type strain PG45 and two strains isolated in China, have been published. In this study, several Tn4001 based transposon constructs were generated and used to create a M. bovis PG45 insertional mutant library. Direct genome sequencing of 319 independent insertions detected disruptions in 129 genes in M. bovis, 48 of which had homologues in Mycoplasma mycoides subspecies mycoides SC and 99 of which had homologues in Mycoplasma agalactiae. Sixteen genes found to be essential in previous studies on other mycoplasma species were found to be dispensable. Five of these genes have previously been predicted to be part of the core set of 153 essential genes in mycoplasmas. Thus this study has extended the list of non-essential genes of mycoplasmas from that previously generated by studies in other species.

Comparative genomics and phylogenomics of hemotrophic mycoplasmas

Hemotrophic mycoplasmas (hemoplasmas) are a group of animal pathogens of the Mollicutes class. Recently, the genomes of 8 hemoplasmas have been completely sequenced. The aim of this study was to gain a better understanding of their genomic features and relationship to other Mycoplasma species. The genome structure and dynamics of hemoplasmas were analyzed by evaluating gene synteny, adaptive evolution of paralogous gene families (PGF) and horizontal gene transfer (HGT). The Mollicutes class was then phylogenetically analyzed by constructing a distance matrix of the 16S rRNA genes and a phylogenetic tree with 32 conserved, concatenated proteins. Our results suggest that the hemoplasmas have dynamic genomes. The genome size variation (from 547 to 1,545 genes) indicates substantial gene gain/loss throughout evolution. Poorly conserved gene syntenies among hemoplasmas, positional shuffling of paralogous genes between strains, HGT, and codons under positive selection in PGFs were also observed. When compared to other Mollicutes species, the hemoplasmas experienced further metabolic reduction, and the 16S rRNA gene distance matrix of the available mollicutes suggests that these organisms presently constitute the most divergent clade within its class. Our phylogenetic tree of concatenated proteins showed some differences when compared to the 16S rRNA gene tree, but non-mycoplasma organisms, such as Ureaplasma spp. and Mesoplasma spp., continue to branch within Mycoplasma clades. In conclusion, while the hemoplasmas experienced further metabolic shrinkage through gene loss, PGFs with positively selected codons are likely beneficial to these species. Phylogeny of the mollicutes based on 16S rRNA genes or concatenated proteins do not obey the current taxonomy. The metabolism and genetic diversity of the mollicutes, the presence of HGT, and lack of standard for genus circumscription are likely to hinder attempts to classify these organisms based on phylogenetic analyses.

The transcriptome landscape of Prochlorococcus MED4 and the factors for stabilizing the core genome

BACKGROUND

Gene gain and loss frequently occurs in the cyanobacterium Prochlorococcus, a phototroph that numerically dominates tropical and subtropical open oceans. However, little is known about the stabilization of its core genome, which contains approximately 1250 genes, in the context of genome streamlining. Using Prochlorococcus MED4 as a model organism, we investigated the constraints on core genome stabilization using transcriptome profiling.

RESULTS

RNA-Seq technique was used to obtain the transcriptome map of Prochlorococcus MED4, including operons, untranslated regions, non-coding RNAs, and novel genes. Genome-wide expression profiles revealed that three factors contribute to core genome stabilization. First, a negative correlation between gene expression levels and protein evolutionary rates was observed. Highly expressed genes were overrepresented in the core genome but not in the flexible genome. Gene necessity was determined as a second powerful constraint on genome evolution through functional enrichment analysis. Third, quick mRNA turnover may increase corresponding proteins' fidelity among genes that were abundantly expressed. Together, these factors influence core genome stabilization during MED4 genome evolution.

CONCLUSIONS

Gene expression, gene necessity, and mRNA turnover contribute to core genome maintenance during cyanobacterium Prochlorococcus genus evolution.

Short and long-term genome stability analysis of prokaryotic genomes

Background

Gene organization dynamics is actively studied because it provides useful evolutionary information, makes functional annotation easier and often enables to characterize pathogens. There is therefore a strong interest in understanding the variability of this trait and the possible correlations with life-style. Two kinds of events affect genome organization: on one hand translocations and recombinations change the relative position of genes shared by two genomes (i.e. the backbone gene order); on the other, insertions and deletions leave the backbone gene order unchanged but they alter the gene neighborhoods by breaking the syntenic regions. A complete picture about genome organization evolution therefore requires to account for both kinds of events.

Results

We developed an approach where we model chromosomes as graphs on which we compute different stability estimators; we consider genome rearrangements as well as the effect of gene insertions and deletions. In a first part of the paper, we fit a measure of backbone gene order conservation (hereinafter called backbone stability) against phylogenetic distance for over 3000 genome comparisons, improving existing models for the divergence in time of backbone stability. Intra- and inter-specific comparisons were treated separately to focus on different time-scales. The use of multiple genomes of a same species allowed to identify genomes with diverging gene order with respect to their conspecific. The inter-species analysis indicates that pathogens are more often unstable with respect to non-pathogens. In a second part of the text, we show that in pathogens, gene content dynamics (insertions and deletions) have a much more dramatic effect on genome organization stability than backbone rearrangements.

Conclusion

In this work, we studied genome organization divergence taking into account the contribution of both genome order rearrangements and genome content dynamics. By studying species with multiple sequenced genomes available, we were able to explore genome organization stability at different time-scales and to find significant differences for pathogen and non-pathogen species. The output of our framework also allows to identify the conserved gene clusters and/or partial occurrences thereof, making possible to explore how gene clusters assembled during evolution.

Identification of the Treponema pallidum subsp. pallidum TP0092 (RpoE) regulon and its implications for pathogen persistence in the host and syphilis pathogenesis

Bacteria often respond to harmful environmental stimuli with the induction of extracytoplasmic function (ECF) sigma (σ) factors that in turn direct RNA polymerase to transcribe specific groups of response genes (or regulons) to minimize cellular damage and favor adaptation to the changed extracellular milieu. In Treponema pallidum subsp. pallidum, the agent of syphilis, the TP0092 gene is predicted to code for the pathogen's only annotated ECF σ factor, homologous to RpoE, known in Escherichia coli to control a key transduction pathway for maintenance of envelope homeostasis in response to external stress and cell growth. Here we have shown that TP0092 is highly transcribed during experimental syphilis. Furthermore, TP0092 transcription levels significantly increase as infection progresses toward immune clearance of the pathogen, suggesting a role for TP0092 in helping T. pallidum respond to harmful stimuli in the host environment. To investigate this hypothesis, we determined the TP0092 regulon at two different time points during infection using chromatin immunoprecipitation followed by high-throughput sequencing. A total of 22 chromosomal regions, all containing putative TP0092-binding sites and corresponding to as many T. pallidum genes, were identified. Noteworthy among them are the genes encoding desulfoferrodoxin and thioredoxin, involved in detoxification of reactive oxygen species (ROS). Because T. pallidum does not possess other enzymes for ROS detoxification, such as superoxide dismutase, catalase, or glutathione peroxidase, our results suggest that the TP0092 regulon is important in protecting the syphilis spirochete from damage caused by ROS produced at the site of infection during the inflammatory response.

Coping with complexity in metabolic engineering

In the past decade, systems biology has revealed great metabolic and regulatory complexity even in seemingly simple microbial systems. Metabolic engineering aims to control this complexity in order to establish sustainable and economically viable production routes for valuable chemicals. Recent advances in systems-level data generation and modeling of cellular metabolism and regulation together with tremendous progress in synthetic biology will provide the tools to put biotechnologists on the fast track for implementing novel production processes. Great potential lies in the reduction of cellular complexity by orthogonalization of metabolic modules. Here, we review recent advances that will eventually enable metabolic engineers to predict, design, and build streamlined microbial cell factories with reduced time and effort.

Streamlining and core genome conservation among highly divergent members of the SAR11 clade

SAR11 is an ancient and diverse clade of heterotrophic bacteria that are abundant throughout the world’s oceans, where they play a major role in the ocean carbon cycle. Correlations between the phylogenetic branching order and spatiotemporal patterns in cell distributions from planktonic ocean environments indicate that SAR11 has evolved into perhaps a dozen or more specialized ecotypes that span evolutionary distances equivalent to a bacterial order. We isolated and sequenced genomes from diverse SAR11 cultures that represent three major lineages and encompass the full breadth of the clade. The new data expand observations about genome evolution and gene content that previously had been restricted to the SAR11 Ia subclade, providing a much broader perspective on the clade’s origins, evolution, and ecology. We found small genomes throughout the clade and a very high proportion of core genome genes (48 to 56%), indicating that small genome size is probably an ancestral characteristic. In their level of core genome conservation, the members of SAR11 are outliers, the most conserved free-living bacteria known. Shared features of the clade include low GC content, high gene synteny, a large hypervariable region bounded by rRNA genes, and low numbers of paralogs. Variation among the genomes included genes for phosphorus metabolism, glycolysis, and C1 metabolism, suggesting that adaptive specialization in nutrient resource utilization is important to niche partitioning and ecotype divergence within the clade. These data provide support for the conclusion that streamlining selection for efficient cell replication in the planktonic habitat has occurred throughout the evolution and diversification of this clade.

The SAR11 clade is the most abundant group of marine microorganisms worldwide, making them key players in the global carbon cycle. Growing knowledge about their biochemistry and metabolism is leading to a more mechanistic understanding of organic carbon oxidation and sequestration in the oceans. The discovery of small genomes in SAR11 provided crucial support for the theory that streamlining selection can drive genome reduction in low-nutrient environments. Study of isolates in culture revealed atypical organic nutrient requirements that can be attributed to genome reduction, such as conditional auxotrophy for glycine and its precursors, a requirement for reduced sulfur compounds, and evidence for widespread cycling of C1 compounds in marine environments. However, understanding the genetic variation and distribution of such pathways and characteristics like streamlining throughout the group has required the isolation and genome sequencing of diverse SAR11 representatives, an analysis of which we provide here.