Comparative Genometrics (CG): a database dedicated to biometric comparisons of whole genomes

The ever increasing rate at which whole genome sequences are becoming accessible to the scientific community has created an urgent need for tools enabling comparison of chromosomes of different species. We have applied biometric methods to available chromosome sequences and posted the results on our Comparative Genometrics (CG) web site. By genometrics, a term coined by Elston and Wilson [Genet. Epidemiol. (1990), 7, 17–19], we understand a biometric analysis of chromosomes. During the initial phase, our web site displays, for all completely sequenced prokaryotic genomes, three genometric analyses: the DNA walk [Lobry (1999) Microbiology Today, 26, 164–165] and two complementary representations, i.e. the cumulative GC- and TA-skew analyses, capable of identifying, at the level of whole genomes, features inherent to chromosome organization and functioning. It appears that the latter features are taxon-specific. Although primarily focused on prokaryotic chromosomes, the CG web site contains genometric information on paradigm plasmids, phages, viruses and eukaryotic organelles. Relevant data and methods can be readily used by the scientific community for further analyses as well as for tutorial purposes. Our data posted at the CG web site are freely available on the World Wide Web at http://www.unil.ch/comparativegenometrics.     

Over 50,000 Metagenomically Assembled Draft Genomes for the Human Oral Microbiome Reveal New Taxa

The oral cavity of each person is home to hundreds of bacterial species. While taxa for oral diseases have been studied using culture-based characterization as well as amplicon sequencing,metagenomic and genomic information remains scarce compared to the fecal microbiome. Here,using metagenomic shotgun data for 3346 oral metagenomic samples together with 808 published samples, we obtain 56,213 metagenome-assembled genomes(MAGs), and more than 64% of the3589 species-level genome bins(SGBs) contai...     

Nucleomorph and plastid genome sequences of the chlorarachniophyte Lotharella oceanica: convergent reductive evolution and frequent recombination in nucleomorph-bearing algae

Background

Nucleomorphs are residual nuclei derived from eukaryotic endosymbionts in chlorarachniophyte and cryptophyte algae. The endosymbionts that gave rise to nucleomorphs and plastids in these two algal groups were green and red algae, respectively. Despite their independent origin, the chlorarachniophyte and cryptophyte nucleomorph genomes share similar genomic features such as extreme size reduction and a three-chromosome architecture. This suggests that similar reductive evolutionary forces have acted to shape the nucleomorph genomes in the two groups. Thus far, however, only a single chlorarachniophyte nucleomorph and plastid genome has been sequenced, making broad evolutionary inferences within the chlorarachniophytes and between chlorarachniophytes and cryptophytes difficult. We have sequenced the nucleomorph and plastid genomes of the chlorarachniophyte Lotharella oceanica in order to gain insight into nucleomorph and plastid genome diversity and evolution.

Results

The L. oceanica nucleomorph genome was found to consist of three linear chromosomes totaling ~610 kilobase pairs (kbp), much larger than the 373 kbp nucleomorph genome of the model chlorarachniophyte Bigelowiella natans. The L. oceanica plastid genome is 71 kbp in size, similar to that of B. natans. Unexpectedly long (~35 kbp) sub-telomeric repeat regions were identified in the L. oceanica nucleomorph genome; internal multi-copy regions were also detected. Gene content analyses revealed that nucleomorph house-keeping genes and spliceosomal intron positions are well conserved between the L. oceanica and B. natans nucleomorph genomes. More broadly, gene retention patterns were found to be similar between nucleomorph genomes in chlorarachniophytes and cryptophytes. Chlorarachniophyte plastid genomes showed near identical protein coding gene complements as well as a high level of synteny.

Conclusions

We have provided insight into the process of nucleomorph genome evolution by elucidating the fine-scale dynamics of sub-telomeric repeat regions. Homologous recombination at the chromosome ends appears to be frequent, serving to expand and contract nucleomorph genome size. The main factor influencing nucleomorph genome size variation between different chlorarachniophyte species appears to be expansion-contraction of these telomere-associated repeats rather than changes in the number of unique protein coding genes. The dynamic nature of chlorarachniophyte nucleomorph genomes lies in stark contrast to their plastid genomes, which appear to be highly stable in terms of gene content and synteny.

Electronic supplementary material

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

Investigations of oligonucleotide usage variance within and between prokaryotes

Oligonucleotide usage in archaeal and bacterial genomes can be linked to a number of properties, including codon usage (trinucleotides), DNA base-stacking energy (dinucleotides), and DNA structural conformation (di- to tetranucleotides). We wanted to assess the statistical information potential of different DNA ‘word-sizes’ and explore how oligonucleotide frequencies differ in coding and non-coding regions. In addition, we used oligonucleotide frequencies to investigate DNA composition and how DNA sequence patterns change within and between prokaryotic organisms. Among the results found was that prokaryotic chromosomes can be described by hexanucleotide frequencies, suggesting that prokaryotic DNA is predominantly short range correlated, i.e., information in prokaryotic genomes is encoded in short oligonucleotides. Oligonucleotide usage varied more within AT-rich and host-associated genomes than in GC-rich and free-living genomes, and this variation was mainly located in non-coding regions. Bias (selectional pressure) in tetranucleotide usage correlated with GC content, and coding regions were more biased than non-coding regions. Non-coding regions were also found to be approximately 5.5% more AT-rich than coding regions, on average, in the 402 chromosomes examined. Pronounced DNA compositional differences were found both within and between AT-rich and GC-rich genomes. GC-rich genomes were more similar and biased in terms of tetranucleotide usage in non-coding regions than AT-rich genomes. The differences found between AT-rich and GC-rich genomes may possibly be attributed to lifestyle, since tetranucleotide usage within host-associated bacteria was, on average, more dissimilar and less biased than free-living archaea and bacteria.     

Genomics and phylogeny of the proposed phylum ‘Candidatus Poribacteria’ associated with the excavating sponge Thoosa mismalolli

Members of the proposed phylum ‘Candidatus Poribacteria’ are among the most abundant microorganisms in the highly diverse microbiome of the sponge mesohyl. Genomic and phylogenetic characteristics of this proposed phylum are barely known. In this study, we analyzed metagenome-assembled genomes (MAGs) obtained from the coral reef excavating sponge Thoosa mismalolli from the Mexican Pacific Ocean. Two MAGs were extracted and analyzed together with 32 MAGs and single-amplified genomes (SAGs) obtained from NCBI. The phylogenetic tree based on the sequences of 139 single-copy genes (SCG) showed two clades. Clade A (23 genomes) represented 67.7% of the total of the genomes, while clade B (11 genomes) comprised 32.3% of the genomes. The Average Nucleotide Identity (ANI) showed values between 66 and 99% for the genomes of the proposed phylum, and the pangenome of genomes revealed a total of 37,234 genes that included 1722 core gene. The number of genes used in the phylogenetic analysis increased from 28 (previous studies) to 139 (this study), which allowed a better resolution of the phylogeny of the proposed phylum. The results supported the two previously described classes, ‘Candidatus Entoporibacteria’ and ‘Candidatus Pelagiporibacteria’, and the genomes SB0101 and SB0202 obtained in this study belong to two new species of the class ‘Candidatus Entoporibacteria’. This is the first comparative study that includes MAGs from a non-sponge host (Porites lutea) to elucidate the taxonomy of the poorly known Candidatus phylum in a polyphasic approach. Finally, our study also contributes to the sponge microbiome project by reporting the first MAGs of the proposed phylum ‘Candidatus Poribacteria’ isolated from the excavating sponge T. mismalolli.

     

Chromosome dynamics in multichromosome bacteria

On the basis of limited information, bacteria were once assumed to have no more than one chromosome. In the era of genomics, it has become clear that some, like eukaryotes, have more than one chromosome. Multichromosome bacteria provide opportunities to investigate how split genomes emerged, whether the individual chromosomes communicate to coordinate their replication and segregation, and what selective advantages split genomes might provide. Our current knowledge of these topics comes mostly from studies in Vibrio cholerae, which has two chromosomes, chr1 and chr2. Chr1 carries out most of the house-keeping functions and is considered the main chromosome, whereas chr2 appears to have originated from a plasmid and has acquired genes of mostly unknown origin and function. Nevertheless, unlike plasmids, chr2 replicates once and only once per cell cycle, like a bona fide chromosome. The two chromosomes replicate and segregate using separate programs, unlike eukaryotic chromosomes. They terminate replication synchronously, suggesting that there might be communication between them. Replication of the chromosomes is affected by segregation genes but in a chromosome specific fashion, a new development in the field of DNA replication control. The split genome allows genome duplication to complete in less time and with fewer replication forks, which could be beneficial for genome maintenance during rapid growth, which is the norm for V. cholerae in broth cultures and in the human host. In the latter, the expression of chr2 genes increases preferentially. Studies of chromosome maintenance in multichromosomal bacteria, although in their infancy, are already broadening our view of chromosome biology. This article is part of a Special Issue entitled: Chromatin in time and space.     

Classification of Bacteria and Archaea: Past,present and future

The late 19th century was the beginning of bacterial taxonomy and bacteria were classified on the basis of phenotypic markers. The distinction of prokaryotes and eukaryotes was introduced in the 1960s. Numerical taxonomy improved phenotypic identification but provided little information on the phylogenetic relationships of prokaryotes. Later on, chemotaxonomic and genotypic methods were widely used for a more satisfactory classification. Archaea were first classified as a separate group of prokaryotes in 1977. The current classification of Bacteria and Archaea is based on an operational-based model, the so-called polyphasic approach, comprised of phenotypic, chemotaxonomic and genotypic data, as well as phylogenetic information. The provisional status Candidatus has been established for describing uncultured prokaryotic cells for which their phylogenetic relationship has been determined and their authenticity revealed by in situ probing.     

Gene identification in novel eukaryotic genomes by self-training algorithm

Finding new protein-coding genes is one of the most important goals of eukaryotic genome sequencing projects. However, genomic organization of novel eukaryotic genomes is diverse and ab initio gene finding tools tuned up for previously studied species are rarely suitable for efficacious gene hunting in DNA sequences of a new genome. Gene identification methods based on cDNA and expressed sequence tag (EST) mapping to genomic DNA or those using alignments to closely related genomes rely either on existence of abundant cDNA and EST data and/or availability on reference genomes. Conventional statistical ab initio methods require large training sets of validated genes for estimating gene model parameters. In practice, neither one of these types of data may be available in sufficient amount until rather late stages of the novel genome sequencing. Nevertheless, we have shown that gene finding in eukaryotic genomes could be carried out in parallel with statistical models estimation directly from yet anonymous genomic DNA. The suggested method of parallelization of gene prediction with the model parameters estimation follows the path of the iterative Viterbi training. Rounds of genomic sequence labeling into coding and non-coding regions are followed by the rounds of model parameters estimation. Several dynamically changing restrictions on the possible range of model parameters are added to filter out fluctuations in the initial steps of the algorithm that could redirect the iteration process away from the biologically relevant point in parameter space. Tests on well-studied eukaryotic genomes have shown that the new method performs comparably or better than conventional methods where the supervised model training precedes the gene prediction step. Several novel genomes have been analyzed and biologically interesting findings are discussed. Thus, a self-training algorithm that had been assumed feasible only for prokaryotic genomes has now been developed for ab initio eukaryotic gene identification.     

真核生物DNA非编码区的组分分析

在全基因组水平上,用直方图、混沌表示灰度图、距离差异度和信息熵差异度四种方法,研究了拟南芥、线虫、果蝇的DNA内含子、基因间隔区DNA、外显子三种区域的核苷酸短序列组分及组分复杂度.结果表明：a.不同基因组之间,不管基因数目多少,用4种方法得到的外显子部分其组分复杂度都比较接近,而非编码区部分的组分复杂度却很大.这一点定量地说明了物种之间的复杂程度,主要不体现在编码区部分,而体现在非编码区部分.b.同一基因组中,内含子的核苷酸短序列组分复杂度都是相似的,外显子和intergenic DNA部分的组分复杂度也是相似的.c.内含子和intergenic DNA在转录、剪切、二级结构等方面有很大的不同,但它们在核苷酸短序列组分上的差异却很小,说明内含子和intergenic DNA在转录、剪切、二级结构上的不同并不通过核苷酸短序列组分来进行限制.     

A Novel Prokaryotic Promoter Identified in the Genome of Some Monopartite Begomoviruses

Geminiviruses are known to exhibit both prokaryotic and eukaryotic features in their genomes, with the ability to express their genes and even replicate in bacterial cells. We have demonstrated previously the existence of unit-length single-stranded circular DNAs of Ageratum yellow vein virus (AYVV, a species in the genus Begomovirus, family Geminiviridae) in Escherichia coli cells, which prompted our search for unknown prokaryotic functions in the begomovirus genomes. By using a promoter trapping strategy, we identified a novel prokaryotic promoter, designated AV3 promoter, in nts 762-831 of the AYVV genome. Activity assays revealed that the AV3 promoter is strong, unidirectional, and constitutive, with an endogenous downstream ribosome binding site and a translatable short open reading frame of eight amino acids. Sequence analyses suggested that the AV3 promoter might be a remnant of prokaryotic ancestors that could be related to certain promoters of bacteria from marine or freshwater environments. The discovery of the prokaryotic AV3 promoter provided further evidence for the prokaryotic origin in the evolutionary history of geminiviruses.     

‘Candidatus ferrigenium straubiae’ sp. nov., ‘Candidatus ferrigenium bremense’ sp. nov., ‘Candidatus ferrigenium altingense’ sp. nov., are autotrophic Fe(II)-oxidizing bacteria of the family Gallionellaceae

Iron(II) [Fe(II)] oxidation coupled to denitrification is recognized as an environmentally important process in many ecosystems. However, the Fe(II)-oxidizing bacteria (FeOB) dominating autotrophic nitrate-reducing Fe(II)-oxidizing enrichment cultures, affiliated with the family Gallionellaceae, remain poorly taxonomically defined due to lack of representative isolates. We describe the taxonomic classification of three novel FeOB based on metagenome-assembled genomes (MAGs) acquired from the autotrophic nitrate-reducing enrichment cultures KS, BP and AG. Phylogenetic analysis of nearly full-length 16S rRNA gene sequences demonstrated that these three FeOB were most closely affiliated to the genera Ferrigenium, Sideroxydans and Gallionella, with up to 96.5%, 95.4% and 96.2% 16S rRNA gene sequence identities to representative isolates of these genera, respectively. In addition, average amino acid identities (AAI) of the genomes compared to the most closely related genera revealed highest AAI with Ferrigenium kumadai An22 (76.35–76.74%), suggesting that the three FeOB are members of this genus. Phylogenetic analysis of conserved functional genes further supported that these FeOB represent three novel species of the genus Ferrigenium. Moreover, the three novel FeOB likely have characteristic features, performing partial denitrification coupled to Fe(II) oxidation and carbon fixation. Scanning electron microscopy of the enrichment cultures showed slightly curved rod-shaped cells, ranging from 0.2-0.7 μm in width and 0.5–2.3 μm in length. Based on the phylogenetic, genomic and physiological characteristics, we propose that these FeOB represent three novel species, ‘Candidatus Ferrigenium straubiae’ sp. nov., ‘Candidatus Ferrigenium bremense’ sp. nov. and ‘Candidatus Ferrigenium altingense’ sp. nov. that might have unique metabolic features among the genus Ferrigenium.     

Environmental genomics, the big picture?

The enormous sequencing capabilities of our times might be reaching the point of overflowing the possibilities to analyse data and allow for a feedback on where to focus the available resources. We have now a foreseeable future in which most bacterial species will have an annotated genome. However, we know also that most prokaryotic diversity would not be included there. On the one hand, there is the problem of many groups not being easily amenable to culture and hence not represented in culture-centred microbial taxonomy. On the other hand, the gene pools present in one species can be orders of magnitude larger than the genome of one strain (selected for genome sequencing). Contrasting with eukaryotic genomes, the repertoire of genes present in one prokaryotic cell genome does not correlate stringently with its taxonomic identity. Hence gene catalogues from one environment might provide more meaningful information than the classical species catalogues. Metagenomics or microbial environmental genomics provide a different tool that gravitates around the habitat rather than the species. Such a tool could be just the right way to complement "organismal genomics". Its potential to advance our understanding of microbial ecology and prokaryotic diversity and evolution is discussed.     

Patterns of prokaryotic lateral gene transfers affecting parasitic microbial eukaryotes

Background

The influence of lateral gene transfer on gene origins and biology in eukaryotes is poorly understood compared with those of prokaryotes. A number of independent investigations focusing on specific genes, individual genomes, or specific functional categories from various eukaryotes have indicated that lateral gene transfer does indeed affect eukaryotic genomes. However, the lack of common methodology and criteria in these studies makes it difficult to assess the general importance and influence of lateral gene transfer on eukaryotic genome evolution.

Results

We used a phylogenomic approach to systematically investigate lateral gene transfer affecting the proteomes of thirteen, mainly parasitic, microbial eukaryotes, representing four of the six eukaryotic super-groups. All of the genomes investigated have been significantly affected by prokaryote-to-eukaryote lateral gene transfers, dramatically affecting the enzymes of core pathways, particularly amino acid and sugar metabolism, but also providing new genes of potential adaptive significance in the life of parasites. A broad range of prokaryotic donors is involved in such transfers, but there is clear and significant enrichment for bacterial groups that share the same habitats, including the human microbiota, as the parasites investigated.

Conclusions

Our data show that ecology and lifestyle strongly influence gene origins and opportunities for gene transfer and reveal that, although the outlines of the core eukaryotic metabolism are conserved among lineages, the genes making up those pathways can have very different origins in different eukaryotes. Thus, from the perspective of the effects of lateral gene transfer on individual gene ancestries in different lineages, eukaryotic metabolism appears to be chimeric.     

Knowledge-based discovery for designing CRISPR-CAS systems against invading mobilomes in thermophiles

Clustered regularly interspaced short palindromic repeats (CRISPRs) are direct features of the prokaryotic genomes involved in resistance to their bacterial viruses and phages. Herein, we have identified CRISPR loci together with CRISPR-associated sequences (CAS) genes to reveal their immunity against genome invaders in the thermophilic archaea and bacteria. Genomic survey of this study implied that genomic distribution of CRISPR-CAS systems was varied from strain to strain, which was determined by the degree of invading mobiloms. Direct repeats found to be equal in some extent in many thermopiles, but their spacers were differed in each strain. Phylogenetic analyses of CAS superfamily revealed that genes cmr, csh, csx11, HD domain, devR were belonged to the subtypes of cas gene family. The members in cas gene family of thermophiles were functionally diverged within closely related genomes and may contribute to develop several defense strategies. Nevertheless, genome dynamics, geological variation and host defense mechanism were contributed to share their molecular functions across the thermophiles. A thermophilic archaean, Thermococcus gammotolerans and thermophilic bacteria, Petrotoga mobilis and Thermotoga lettingae have shown superoperons-like appearance to cluster cas genes, which were typically evolved for their defense pathways. A cmr operon was identified with a specific promoter in a thermophilic archaean, Caldivirga maquilingensis. Overall, we concluded that knowledge-based genomic survey and phylogeny-based functional assignment have suggested for designing a reliable genetic regulatory circuit naturally from CRISPR-CAS systems, acquired defense pathways, to thermophiles in future synthetic biology.

Electronic supplementary material

The online version of this article (doi:10.1007/s11693-015-9176-8) contains supplementary material, which is available to authorized users.     

‘Genomemark’: Detecting Word Periodicity in Biological Sequences

Abstract

Identifying and predicting the structural characteristics of novel repeats throughout the genome can lend insight into biological function. Specific repeats are believed to have biological significance as a function of their distribution patterns. We have developed ‘GenomeMark,’ a computer program that detects and statistically analyzes candidate repeats. Specifically, ‘GenomeMark’ identifies the periodic distribution of unique words, calculating their χ² and Z-score values. Using ‘GenomeMark,’ we identified novel sequence words present in tandem throughout genomes. We found that these sequences have remarkable spacer sequence distributions and many were genome specific, validating the genome signature theory. Further analysis confirmed that many of these sequences have a specific biological function. The program is available from the authors upon request and is freely available for non-commercial and academic entities.     

Distribution of centromere-like parS sites in bacteria: insights from comparative genomics

Partitioning of low-copy-number plasmids to daughter cells often depends on ParA and ParB proteins acting on centromere-like parS sites. Similar chromosome-encoded par loci likely also contribute to chromosome segregation. Here, we used bioinformatic approaches to search for chromosomal parS sites in 400 prokaryotic genomes. Although the consensus sequence matrix used to search for parS sites was derived from two gram-positive species, putative parS sites were identified on the chromosomes of 69% of strains from all branches of bacteria. Strains that were not found to contain parS sites clustered among relatively few branches of the prokaryotic evolutionary tree. In the vast majority of cases, parS sites were identified in origin-proximal regions of chromosomes. The widespread conservation of parS sites across diverse bacteria suggests that par loci evolved very early in the evolution of bacterial chromosomes and that the absence of parS, parA, and/or parB in certain strains likely reflects the loss of one of more of these loci much later in evolution. Moreover, the highly conserved origin-proximal position of parS suggests par loci are primarily devoted to regulating processes that involve the origin region of bacterial chromosomes. In species containing multiple chromosomes, the parS sites found on secondary chromosomes diverge significantly from those found on their primary chromosomes, suggesting that chromosome segregation of multipartite genomes requires distinct replicon-specific par loci. Furthermore, parS sites on secondary chromosomes are not well conserved among different species, suggesting that the evolutionary histories of secondary chromosomes are more diverse than those of primary chromosomes.     

Weeding out the genes: the Arabidopsis genome project

The Arabidopsis genome sequence is scheduled for completion at the end of this year (December 2000). It will be the first higher plant genome to be sequenced, and will allow a detailed comparison with bacterial, yeast and animal genomes. Already, two of the five chromosomes have been sequenced, and we have had our first glimpse of higher eukaryotic centromeres, and the structure of heterochromatin. The implications for understanding plant gene function, genome structure and genome organization are profound. In this review, the lessons learned for future genome projects are reviewed as well as a summary of the initial findings in Arabidopsis. Electronic Publication     

Candidatus Abditibacter,a novel genus within the Cryomorphaceae,thriving in the North Sea

Coastal phytoplankton blooms are frequently followed by successive blooms of heterotrophic bacterial clades. The class Flavobacteriia within the Bacteroidetes has been shown to play an important role in the degradation of high molecular weight substrates that become available in the later stages of such blooms. One of the flavobacterial clades repeatedly observed over the course of several years during phytoplankton blooms off the coast of Helgoland, North Sea, is Vis6. This genus-level clade belongs to the family Cryomorphaceae and has been resistant to cultivation to date. Based on metagenome assembled genomes, comparative 16S rRNA gene sequence analyses and fluorescence in situ hybridization, we here propose a novel candidate genus Abditibacter, comprising three novel species Candidatus Abditibacter vernus, Candidatus Abditibacter forsetii and Candidatus Abditibacter autumni. While the small genomes of the three novel photoheterotrophic species encode highly similar gene repertoires, including genes for degradation of proteins and algal storage polysaccharides such as laminarin, two of them – Ca. A. vernus and Ca. A. forsetii – seem to have a preference for spring blooms, while Ca. A. autumni almost exclusively occurs in late summer and autumn.