Genome-wide and single-base resolution DNA methylomes of the Pacific oyster Crassostrea gigas provide insight into the evolution of invertebrate CpG methylation

Background

Studies of DNA methylomes in a wide range of eukaryotes have revealed both conserved and divergent characteristics of DNA methylation among phylogenetic groups. However, data on invertebrates particularly molluscs are limited, which hinders our understanding of the evolution of DNA methylation in metazoa. The sequencing of the Pacific oyster Crassostrea gigas genome provides an opportunity for genome-wide profiling of DNA methylation in this model mollusc.

Results

Homologous searches against the C. gigas genome identified functional orthologs for key genes involved in DNA methylation: DNMT1, DNMT2, DNMT3, MBD2/3 and UHRF1. Whole-genome bisulfite sequencing (BS-seq) of the oyster’s mantle tissues revealed that more than 99% methylation modification was restricted to cytosines in CpG context and methylated CpGs accumulated in the bodies of genes that were moderately expressed. Young repeat elements were another major targets of CpG methylation in oysters. Comparison with other invertebrate methylomes suggested that the 5’-end bias of gene body methylation and the negative correlation between gene body methylation and gene length were the derived features probably limited to the insect lineage. Interestingly, phylostratigraphic analysis showed that CpG methylation preferentially targeted genes originating in the common ancestor of eukaryotes rather than the oldest genes originating in the common ancestor of cellular organisms.

Conclusions

Comparative analysis of the oyster DNA methylomes and that of other animal species revealed that the characteristics of DNA methylation were generally conserved during invertebrate evolution, while some unique features were derived in the insect lineage. The preference of methylation modification on genes originating in the eukaryotic ancestor rather than the oldest genes is unexpected, probably implying that the emergence of methylation regulation in these ''relatively young’ genes was critical for the origin and radiation of eukaryotes.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1119) contains supplementary material, which is available to authorized users.     

Genome-wide and molecular evolution analysis of the subtilase gene family in Vitis vinifera

Background

Vitis vinifera (grape) is one of the most economically significant fruit crops in the world. The availability of the recently released grape genome sequence offers an opportunity to identify and analyze some important gene families in this species. Subtilases are a group of subtilisin-like serine proteases that are involved in many biological processes in plants. However, no comprehensive study incorporating phylogeny, chromosomal location and gene duplication, gene organization, functional divergence, selective pressure and expression profiling has been reported so far for the grape.

Results

In the present study, a comprehensive analysis of the subtilase gene family in V. vinifera was performed. Eighty subtilase genes were identified. Phylogenetic analyses indicated that these subtilase genes comprised eight groups. The gene organization is considerably conserved among the groups. Distribution of the subtilase genes is non-random across the chromosomes. A high proportion of these genes are preferentially clustered, indicating that tandem duplications may have contributed significantly to the expansion of the subtilase gene family. Analyses of divergence and adaptive evolution show that while purifying selection may have been the main force driving the evolution of grape subtilases, some of the critical sites responsible for the divergence may have been under positive selection. Further analyses of real-time PCR data suggested that many subtilase genes might be important in the stress response and functional development of plants.

Conclusions

Tandem duplications as well as purifying and positive selections have contributed to the functional divergence of subtilase genes in V. vinifera. The data may contribute to a better understanding of the grape subtilase gene family.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1116) contains supplementary material, which is available to authorized users.     

Clostridium sordellii genome analysis reveals plasmid localized toxin genes encoded within pathogenicity loci

Background

Clostridium sordellii can cause severe infections in animals and humans, the latter associated with trauma, toxic shock and often-fatal gynaecological infections. Strains can produce two large clostridial cytotoxins (LCCs), TcsL and TcsH, related to those produced by Clostridium difficile, Clostridium novyi and Clostridium perfringens, but the genetic basis of toxin production remains uncharacterised.

Results

Phylogenetic analysis of the genome sequences of 44 strains isolated from human and animal infections in the UK, US and Australia placed the species into four clades. Although all strains originated from animal or clinical disease, only 5 strains contained LCC genes: 4 strains contain tcsL alone and one strain contains tcsL and tcsH. Four toxin-positive strains were found within one clade. Where present, tcsL and tcsH were localised in a pathogenicity locus, similar to but distinct from that present in C. difficile. In contrast to C. difficile, where the LCCs are chromosomally localised, the C. sordellii tcsL and tcsH genes are localised on plasmids. Our data suggest gain and loss of entire toxigenic plasmids in addition to horizontal transfer of the pathogenicity locus. A high quality, annotated sequence of ATCC9714 reveals many putative virulence factors including neuraminidase, phospholipase C and the cholesterol-dependent cytolysin sordellilysin that are highly conserved between all strains studied.

Conclusions

Genome analysis of C. sordellii reveals that the LCCs, the major virulence factors, are localised on plasmids. Many strains do not contain the LCC genes; it is probable that in several of these cases the plasmid has been lost upon laboratory subculture. Our data are consistent with LCCs being the primary virulence factors in the majority of infections, but LCC-negative strains may precipitate certain categories of infection. A high quality genome sequence reveals putative virulence factors whose role in virulence can be investigated.

Electronic supplementary material

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

Proteogenomic mapping of Mycoplasma hyopneumoniae virulent strain 232

Background

Mycoplasma hyopneumoniae causes respiratory disease in swine and contributes to the porcine respiratory disease complex, a major disease problem in the swine industry. The M. hyopneumoniae strain 232 genome is one of the smallest and best annotated microbial genomes, containing only 728 annotated genes and 691 known proteins. Standard protein databases for mass spectrometry only allow for the identification of known and predicted proteins, which if incorrect can limit our understanding of the biological processes at work. Proteogenomic mapping is a methodology which allows the entire 6-frame genome translation of an organism to be used as a mass spectrometry database to help identify unknown proteins as well as correct and confirm existing annotations. This methodology will be employed to perform an in-depth analysis of the M. hyopneumoniae proteome.

Results

Proteomic analysis indicates 483 of 691 (70%) known M. hyopneumoniae strain 232 proteins are expressed under the culture conditions given in this study. Furthermore, 171 of 328 (52%) hypothetical proteins have been confirmed. Proteogenomic mapping resulted in the identification of previously unannotated genes gatC and rpmF and 5-prime extensions to genes mhp063, mhp073, and mhp451, all conserved and annotated in other M. hyopneumoniae strains and Mycoplasma species. Gene prediction with Prodigal, a prokaryotic gene predicting program, completely supports the new genomic coordinates calculated using proteogenomic mapping.

Conclusions

Proteogenomic mapping showed that the protein coding genes of the M. hyopneumoniae strain 232 identified in this study are well annotated. Only 1.8% of mapped peptides did not correspond to genes defined by the current genome annotation. This study also illustrates how proteogenomic mapping can be an important tool to help confirm, correct and append known gene models when using a genome sequence as search space for peptide mass spectra. Using a gene prediction program which scans for a wide variety of promoters can help ensure genes are accurately predicted or not missed completely. Furthermore, protein extraction using differential detergent fractionation effectively increases the number of membrane and cytoplasmic proteins identifiable my mass spectrometry.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-576) contains supplementary material, which is available to authorized users.     

Complete mitochondrial genomes of the human follicle mites Demodex brevis and D. folliculorum: novel gene arrangement,truncated tRNA genes,and ancient divergence between species

Background

Follicle mites of the genus Demodex are found on a wide diversity of mammals, including humans; surprisingly little is known, however, about the evolution of this association. Additional sequence information promises to facilitate studies of Demodex variation within and between host species. Here we report the complete mitochondrial genome sequences of two species of Demodex known to live on humans—Demodex brevis and D. folliculorum—which are the first such genomes available for any member of the genus. We analyzed these sequences to gain insight into the evolution of mitochondrial genomes within the Acariformes. We also used relaxed molecular clock analyses, based on alignments of mitochondrial proteins, to estimate the time of divergence between these two species.

Results

Both Demodex genomes shared a novel gene order that differs substantially from the ancestral chelicerate pattern, with transfer RNA (tRNA) genes apparently having moved much more often than other genes. Mitochondrial tRNA genes of both species were unusually short, with most of them unable to encode tRNAs that could fold into the canonical cloverleaf structure; indeed, several examples lacked both D- and T-arms. Finally, the high level of sequence divergence observed between these species suggests that these two lineages last shared a common ancestor no more recently than about 87 mya.

Conclusions

Among Acariformes, rearrangements involving tRNA genes tend to occur much more often than those involving other genes. The truncated tRNA genes observed in both Demodex species would seem to require the evolution of extensive tRNA editing capabilities and/or coevolved interacting factors. The molecular machinery necessary for these unusual tRNAs to function might provide an avenue for developing treatments of skin disorders caused by Demodex. The deep divergence time estimated between these two species sets a lower bound on the time that Demodex have been coevolving with their mammalian hosts, and supports the hypothesis that there was an early split within the genus Demodex into species that dwell in different skin microhabitats.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1124) contains supplementary material, which is available to authorized users.