Dataset Information

Comparative genomic analysis of Mycoplasma anatis strains.

ABSTRACT:

Background

The Gram-negative intracellular bacterium Mycoplasma anatis is a pathogen of respiratory infectious diseases in ducks and has caused significant economic losses in the poultry industry.

Objective

This study, as the first report of the structure and function of the pan-genome of Mycoplasma anatis, may provide a valuable genetic basis for many aspects of future research on the pathogens of waterfowl.

Methods

We sequenced the whole genomes of 15 Mycoplasma anatis isolated from ducks in China. Draft genome sequencing was carried out and whole-genome sequencing was performed by the sequencers of the PacBio Sequel and an IonTorrent Personal Genome Machine (PGM). Then the common genic elements of protein-coding genes, tRNAs, and rRNAs of Mycoplasma anatis genomes were predicted by using the pipeline Prokka v1.13.7. To investigate homologous protein clusters across Mycoplasma anatis genomes, we adopted Roary v3.13.0 to cluster orthologous genes (OGs) based on the following criteria.

Results

We obtained one complete genome and 14 genome sketches. Microbial mobile genetic element analysis revealed the distribution of insertion sequences (IS30, IS3, and IS1634), prophage regions, and CRISPR arrays in the genome of Mycoplasma anatis. Comparative genomic analysis decoded the genetic components and functional classification of the pan-genome of Mycoplasma anatis that comprised 646 core genes, 231 dispensable genes and among them 110 was strain-specific. Virulence-related gene profiles of Mycoplasma anatis were systematically identified, and the products of these genes included bacterial ABC transporter systems, iron transport proteins, toxins, and secretion systems.

Conclusion

A complete virulence-related gene profile of Mycoplasma anatis has been identified, most of the genes are highly conserved in all strains. Sequencing results are relevant to the molecular mechanisms of drug resistance, adaptive evolution of pathogens, population structure, and vaccine development.

SUBMITTER: Zhou Q

PROVIDER: S-EPMC8237044 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Mycoplasma hyopneumoniae is the causative agent of porcine enzootic pneumonia and a major factor in the porcine respiratory disease complex. A clear understanding of the mechanisms of pathogenesis does not exist although it is clear that M. hyopneumoniae adheres to porcine ciliated epithelium by action of a protein called P97. Previous studies have shown variation in the gene encoding the P97cilium adhesin within different strains of M. hyopneumoniae, but the extent of genetic variation among field strains across the genome is not known. Since M. hyopneumoniae is a worldwide problem, it is reasonable to expect that a wide range of genetic variability may exist given all of the different breed and housing conditions. This variation may impact the overall virulence of a single strain. Using microarray technology, this study examined potential variation of fourteen field strains in comparison to strain 232 on which the array was based. Genomic DNA was obtained, amplified with TempliPhiâ¢, and labeled indirectly with Alexa dyes. Post genomic hybridization, the arrays were scanned and data analyzed using a linear statistical model. Results indicate that genetic variation could be detected in all fourteen field strains but across different loci, suggesting that variation occurs throughout the genome. Fifty-nine percent of the variable loci were hypothetical genes. Twenty-two percent of the lipoprotein genes showed variation in at least one field strain. A permutation test identified a location in M. hyopneumoniae genome where spatial clustering of variability between the field strains and strain 232 exists. TempliPhiâ¢ amplified DNA samples from field isolates were compared to control strain 232 using a two-color experimental microarray design. Independent samples from one isolate labeled with one dye were paired with control samples labeled with the alternate dye; the samples were mixed and hybridized to the microarray. For nine of the fourteen isolates (95MP1501, 95MP1502, 95MP1503, 95MP1504, 95MP1508, 95MP1509, 97MP0001, 00MP1301, and 05MP2301), four independent field isolate DNA samples were paired with four independent DNA samples from control 232. In two of the four arrays, the control sample was labeled with Alexa 555 dye and compared to the field isolate sample labeled with Alexa 647 dye (Molecular Probes, Inc., Eugene, Ore.). The dye assignment to control and treated samples was reversed for the other two arrays (dye swap). The arrays were hybridized under identical conditions as described below. This procedure was repeated for isolates 95MP1510 for a total of four arrays, where the control sample was labeled with Alexa 647 dye in three of the arrays and Alexa 555 dye for the fourth array. For isolates 95MP1505, 95MP1506 and 95MP1507, a total of five arrays each, including two dye swaps, were done; and for isolate 00MP1502, a total of six arrays were done with the control labeled with Alexa 555 dye for four of the arrays and Alexa 647 for the other two arrays.

Project description:BACKGROUND: While the genomic era is accumulating a tremendous amount of data, the question of how genomics can describe a bacterial species remains to be fully addressed. The recent sequencing of the genome of the Mycoplasma agalactiae type strain has challenged our general view on mycoplasmas by suggesting that these simple bacteria are able to exchange significant amount of genetic material via horizontal gene transfer. Yet, events that are shaping mycoplasma genomes and that are underlining diversity within this species have to be fully evaluated. For this purpose, we compared two strains that are representative of the genetic spectrum encountered in this species: the type strain PG2 which genome is already available and a field strain, 5632, which was fully sequenced and annotated in this study. RESULTS: The two genomes differ by ca. 130 kbp with that of 5632 being the largest (1006 kbp). The make up of this additional genetic material mainly corresponds (i) to mobile genetic elements and (ii) to expanded repertoire of gene families that encode putative surface proteins and display features of highly-variable systems. More specifically, three entire copies of a previously described integrative conjugative element are found in 5632 that accounts for ca. 80 kbp. Other mobile genetic elements, found in 5632 but not in PG2, are the more classical insertion sequences which are related to those found in two other ruminant pathogens, M. bovis and M. mycoides subsp. mycoides SC. In 5632, repertoires of gene families encoding surface proteins are larger due to gene duplication. Comparative proteomic analyses of the two strains indicate that the additional coding capacity of 5632 affects the overall architecture of the surface and suggests the occurrence of new phase variable systems based on single nucleotide polymorphisms. CONCLUSION: Overall, comparative analyses of two M. agalactiae strains revealed a very dynamic genome which structure has been shaped by gene flow among ruminant mycoplasmas and expansion-reduction of gene repertoires encoding surface proteins, the expression of which is driven by localized genetic micro-events.