Comparative genome analysis of Mycoplasma pneumoniae

Background

Mycoplasma pneumoniae is a common pathogen that causes upper and lower respiratory tract infections in people of all ages, responsible for up to 40 % of community-acquired pneumonias. It also causes a wide array of extrapulmonary infections and autoimmune phenomena. Phylogenetic studies of the organism have been generally restricted to specific genes or regions of the genome, because whole genome sequencing has been completed for only 4 strains. To better understand the physiology and pathogenicity of this important human pathogen, we performed comparative genomic analysis of 15 strains of M. pneumoniae that were isolated between the 1940s to 2009 from respiratory specimens and cerebrospinal fluid originating from the USA, China and England.

Results

Illumina MiSeq whole genome sequencing was performed on the 15 strains and all genome sequences were completed. Results from the comparative genomic analysis indicate that although about 1500 SNP and indel variants exist between type1 and type 2 strains, there is an overall high degree of sequence similarity among the strains (>99 % identical to each other). Within the two subtypes, conservation of most genes, including the CARDS toxin gene and arginine deiminase genes, was observed. The major variation occurs in the P1 and ORF6 genes associated with the adhesin complex. Multiple hsdS genes (encodes S subunit of type I restriction enzyme) with variable tandem repeat copy numbers were found in all 15 genomes.

Conclusions

These data indicate that despite conclusions drawn from 16S rRNA sequences suggesting rapid evolution, the M. pneumoniae genome is extraordinarily stable over time and geographic distance across the globe with a striking lack of evidence of horizontal gene transfer.

Electronic supplementary material

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

人类基因组上CpG岛的定位和系统分析

为了研究CpG岛产生和消失机制以及位于基因启动子区域外的CpG岛保守性等问题,我们通过序列比对和进化保守性分析等方法,分析在人类和小鼠中保守的基因上的CpG岛。结果显示已有保守序列的突变以及序列插入删除是CpG岛产生和消失的主要原因,进一步分析发现52%的在小鼠基因组上保守序列完全缺失的CpG岛位于两个转座子之间,提示转座子所介导的序列插入是CpG岛形成和消失的重要原因。人类基因组上在启动子区域外的CpG岛中约有79%为新产生的CpG岛,显著高于启动子区域内新产生的CpG岛比例(41%)。GO分析表明与这些CpG岛相关的部分基因与神经系统发育显著相关,提示新产生的CpG岛参与神经发育过程。     

The Mycoplasma genome

A review is presented of the available data on the nature of chromosomal and extrachromosomal DNA of Mollicutes (mycoplasmas)--the smallest and simplest procaryotic organisms.     

DNA methylation and CpG suppression

Cytosine methylation in vertebrate genomes occurs predominantly at the dinucleotide CpG. This dinucleotide is deficient in vertebrate DNA, an observation which has hitherto been explained by passive deamination of S-methylcytosine to thymidine. Since the frequency and distribution of CpG may prove to be a useful indirect means to study the function of DNA methylation, it is of interest that the observed 'CpG suppression' is less apparent within and around coding sequences. A variety of different mechanisms now appear to be responsible for maintaining a relatively high CpG level in these regions despite the apparent attendant disadvantage of mutation.     

Long-range autocorrelations of CpG islands in the human genome

In this paper, we use a statistical estimator developed in astrophysics to study the distribution and organization of features of the human genome. Using the human reference sequence we quantify the global distribution of CpG islands (CGI) in each chromosome and demonstrate that the organization of the CGI across a chromosome is non-random, exhibits surprisingly long range correlations (10 Mb) and varies significantly among chromosomes. These correlations of CGI summarize functional properties of the genome that are not captured when considering variation in any particular separate (and local) feature. The demonstration of the proposed methods to quantify the organization of CGI in the human genome forms the basis of future studies. The most illuminating of these will assess the potential impact on phenotypic variation of inter-individual variation in the organization of the functional features of the genome within and among chromosomes, and among individuals for particular chromosomes.     

Physical mapping of the Mycoplasma capricolum genome

A physical map of Mycoplasma capricolum ATCC 27343 genome was constructed, based on estimation of the restriction fragment sizes by pulse-field electrophoresis. The linkage order of restriction fragments was determined by two-dimensional electrophoresis of partial and complete single digests and complete double digests and by Southern hybridization analysis. The genome size was established at 1155.5 kb, and 26 cleavage sites for 7 endonucleases were assigned to the map.     

Sequence context analysis in the mouse genome: single nucleotide polymorphisms and CpG island sequences

A genome-wide view of sequence mutability in mice is still limited, although biologists usually assume the same scenario for mice as for humans. In this study, we examined the sequence context in the local environment of 482,528 mouse single nucleotide polymorphisms (SNPs). We found that CpG-containing short sequences, in general, had more representation in the local sequences of SNPs compared to the genome sequences. The extent of this overrepresentation was stronger in mice than in humans, which is inconsistent with previous observations of the weaker neighboring-nucleotide biases on mouse SNPs. To exclude the CpG effect, we compared the distribution patterns of short sequences among the six categories of SNPs. The results revealed an even stronger pattern in the CpG-containing group for C/G substitution compared to for A/G or C/T substitutions. We next performed the first genome-wide sequence context analysis of SNPs in the mouse CpG islands. SNPs occurring at CpG sites were 3.14-fold less prevalent than expected, suggesting the suppression of methylation-dependent deamination in the CpG islands. The extent of this suppression was less in mice than in humans. Finally, compared with humans, the observations of a greater deficit of CpG dinucleotides, a stronger overrepresentation of CpG-containing n-mers surrounding the polymorphic sites, and a higher SNP/genome ratio of CpG dinucleotides in the mouse genome support the "loss of CpG islands" model in the mouse lineage.     

Cloning of the complete Mycoplasma pneumoniae genome.

The complete genome of Mycoplasma pneumoniae was cloned in an ordered library consisting of 34 overlapping or adjacent cosmids, one plasmid and two lambda phages. The genome size was determined by adding up the sizes of either the individual unique EcoRI restriction fragments of the gene bank or of the XhoI fragments of genomic M. pneumoniae DNA. The values from these calculations, 835 and 849 kbp, are in good agreement. An XhoI restriction map was constructed by identifying adjacent DNA fragments by probing with selected cosmid clones.     

A physical map of the Mycoplasma genitalium genome

We report the construction of a physical map of the genome of the human pathogen Mycoplasma genitalium through the use of pulse-field gel electrophoresis. The small size and relative simplicity of this genome permit the arrangement of restriction fragments without having to construct linking clones. The size of the genome has been calculated to be approximately 600 kb and several important genetic determinants have been assigned specific loci on the map.     

Physical mapping of the Mycoplasma pneumoniae genome.

In order to study the genome organization of Mycoplasma pneumoniae a cosmid library of M. pneumoniae DNA was established using a newly designed cosmid vector (pcosRW2). From this library 32 overlapping clones were isolated covering a contiguous 720 kbp DNA segment representing about 90% of the genome assuming a genome size of about 800 kbp.     

DNA repair in reduced genome: the Mycoplasma model

The occurrence of bacteria with a reduced genome, such as that found in Mycoplasmas, raises the question as to which genes should be enough to guarantee the genomic stability indispensable for the maintenance of life. The aim of this work was to compare nine Mycoplasma genomes in regard to DNA repair genes. An in silico analysis was done using six Mycoplasma species, whose genomes are accessible at GenBank, and M. synoviae, and two strains of M. hyopneumoniae, whose genomes were recently sequenced by The Brazilian National Genome Project Consortium and Southern Genome Investigation Program (Brazil) respectively. Considering this reduced genome model, our comparative analysis suggests that the DNA integrity necessary for life can be primarily maintained by nucleotide excision repair (NER), which is the only complete repair pathway. Furthermore, some enzymes involved with base excision repair (BER) and recombination are also present and can complement the NER activity. The absence of RecR and RecO-like ORFs was observed only in M. genitalium and M. pneumoniae, which can be involved with the conservation of gene order observed between these two species. We also obtained phylogenetic evidence for the recent acquisition of the ogt gene in M. pulmonis and M. penetrans by a lateral transference event. In general, the presence or nonexistence of repair genes is shared by all species analyzed, suggesting that the loss of the majority of repair genes was an ancestral event, which occurred before the divergence of the Mycoplasma species.     

Evaluation of single CpG sites as proxies of CpG island methylation states at the genome scale

Methylation of a CpG island is a faithful marker of silencing of its associated gene. Different approaches report the methylation status of a CpG island based on the determination of one or a few CpG sites by assuming the homogeneity of methylation along the element. This strategy is frequently applied in both locus-specific and genome-wide studies, but often without a validation of the representativeness of the interrogated CpG site compared with the whole element. We have evaluated the predictive informativeness of the HpaII sites located in CpG islands using data from high-resolution methylome maps, which offer the possibility to assess the methylation homogeneity of each CpG island and to determine the reporter accuracy of single sites as surrogate markers. An excellent correlation was observed between the HpaII and CpG island methylation levels (r > 0.93). At the qualitative level, the predictive sensitivity of HpaII was >95% with >92% specificity for methylated CpG islands and >90% sensitivity with >95% specificity for unmethylated CpG islands. This analysis provides a global validation framework for strategies based on the use of the methylation-sensitive HpaII restriction enzyme.     

Mapping of replication initiation site in Mycoplasma capricolum genome by two-dimensional gel-electrophoretic analysis.

The homolog of the dnaA gene, which has been reported to be present in the vicinity of the initiation site of replication in the genome of Mycoplasma capricolum (M.Miyata, L.Wang, and T.Fukumura, J. Bacteriol. 175: 655-660, 1993) was mapped precisely. A 9540-bp region containing the dnaA gene was cloned and the entire region was sequenced with the exception of a previously reported region of 2517 bp (Fujita, M.Q., Yoshikawa, H. and Ogasawara, N. Gene 93: 73-78, 1992). The organization of the 9540-bp region was compared with that of corresponding regions in other bacteria. The arrangement and directions of rnpA, rpmH, dnaA, dnaN were conserved, but no other open reading frames were found that were homologous to those that are commonly found around dnaA genes in other bacteria. The directions of movement of the replication fork around the dnaA gene were analyzed by neutral/alkaline two-dimensional gel electrophoresis. The forks developed in a 1569-bp region that consisted of the dnaA structural gene and its downstream non-coding region, and then they proceeded bidirectionally.     

Heterologous CpG island becomes extensively methylated in the genome of transgenic mice

It is demonstrated that a heterologous (chicken) CpG island containing five Sp1 canonical recognition sequences becomes highly methylated in the genome of transgenic mice bearing one or several copies of the transgene. Similar levels of methylation of the chicken CpG island were observed in different tissues of transgenic mice except the brain where the level of methylation of this chicken CpG-rich fragment was significantly lower than in other tissues. Analysis of susceptibility of the "transgenic" CpG island to Hpa II and Msp I restriction nucleases revealed an unusual methylation pattern interfering with the action of both of these enzymes. A conclusion has been drawn that heterologous CpG island per se does not contain all necessary signals permitting to maintain its own non-methylated status in the genome of transgenic animals.     

CpG suppression in vertebrate genomes does not account for the rarity of (CpG)n microsatellite repeats.

Simple microsatellite repetitive sequences are widely distributed in eukaryotic genomes. Using the GCG Find program, the distribution of each type of mono- and dinucleotide repetitive sequence has been examined in GenBank sequences. Examples of each type of simple satellite sequence could be found, although the frequency of (CpG)n greater than or equal to 8 repeats was extremely low. The suppression of CpG dinucleotides in vertebrates does not adequately explain the rarity of this repeat since (CpG)n repeats are also extremely infrequent in species genomes where CpG dinucleotides are not suppressed. Instead, it is proposed that (CpG)n repeats must possess a DNA conformation that has a deleterious structural effect.