Phylogenetic Diversity and Cosymbiosis in the Bioluminescent Symbioses of “Photobacterium mandapamensis”

“Photobacterium mandapamensis” (proposed name) and Photobacterium leiognathi are closely related, phenotypically similar marine bacteria that form bioluminescent symbioses with marine animals. Despite their similarity, however, these bacteria can be distinguished phylogenetically by sequence divergence of their luminescence genes, luxCDAB(F)E, by the presence (P. mandapamensis) or the absence (P. leiognathi) of luxF and, as shown here, by the sequence divergence of genes involved in the synthesis of riboflavin, ribBHA. To gain insight into the possibility that P. mandapamensis and P. leiognathi are ecologically distinct, we used these phylogenetic criteria to determine the incidence of P. mandapamensis as a bioluminescent symbiont of marine animals. Five fish species, Acropoma japonicum (Perciformes, Acropomatidae), Photopectoralis panayensis and Photopectoralis bindus (Perciformes, Leiognathidae), Siphamia versicolor (Perciformes, Apogonidae), and Gadella jordani (Gadiformes, Moridae), were found to harbor P. mandapamensis in their light organs. Specimens of A. japonicus, P. panayensis, and P. bindus harbored P. mandapamensis and P. leiognathi together as cosymbionts of the same light organ. Regardless of cosymbiosis, P. mandapamensis was the predominant symbiont of A. japonicum, and it was the apparently exclusive symbiont of S. versicolor and G. jordani. In contrast, P. leiognathi was found to be the predominant symbiont of P. panayensis and P. bindus, and it appears to be the exclusive symbiont of other leiognathid fishes and a loliginid squid. A phylogenetic test for cospeciation revealed no evidence of codivergence between P. mandapamensis and its host fishes, indicating that coevolution apparently is not the basis for this bacterium's host preferences. These results, which are the first report of bacterial cosymbiosis in fish light organs and the first demonstration that P. leiognathi is not the exclusive light organ symbiont of leiognathid fishes, demonstrate that the host species ranges of P. mandapamensis and P. leiognathi are substantially distinct. The host range difference underscores possible differences in the environmental distributions and physiologies of these two bacterial species.     

The Genome Sequence of ‘Mycobacterium massiliense’ Strain CIP 108297 Suggests the Independent Taxonomic Status of the Mycobacterium abscessus Complex at the Subspecies Level

Members of the Mycobacterium abscessus complex are rapidly growing mycobacteria that are emerging as human pathogens. The M. abscessus complex was previously composed of three species, namely M. abscessus sensu stricto, ‘M. massiliense’, and ‘M. bolletii’. In 2011, ‘M. massiliense’ and ‘M. bolletii’ were united and reclassified as a single subspecies within M. abscessus: M. abscessus subsp. bolletii. However, the placement of ‘M. massiliense’ within the boundary of M. abscessus subsp. bolletii remains highly controversial with regard to clinical aspects. In this study, we revisited the taxonomic status of members of the M. abscessus complex based on comparative analysis of the whole-genome sequences of 53 strains. The genome sequence of the previous type strain of ‘Mycobacterium massiliense’ (CIP 108297) was determined using next-generation sequencing. The genome tree based on average nucleotide identity (ANI) values supported the differentiation of ‘M. bolletii’ and ‘M. massiliense’ at the subspecies level. The genome tree also clearly illustrated that ‘M. bolletii’ and ‘M. massiliense’ form a distinct phylogenetic clade within the radiation of the M. abscessus complex. The genomic distances observed in this study suggest that the current M. abscessus subsp. bolletii taxon should be divided into two subspecies, M. abscessus subsp. massiliense subsp. nov. and M. abscessus subsp. bolletii, to correspondingly accommodate the previously known ‘M. massiliense’ and ‘M. bolletii’ strains.     

Unique Features of a Japanese ‘Candidatus Liberibacter asiaticus’ Strain Revealed by Whole Genome Sequencing

Citrus greening (huanglongbing) is the most destructive disease of citrus worldwide. It is spread by citrus psyllids and is associated with phloem-limited bacteria of three species of α-Proteobacteria, namely, ‘Candidatus Liberibacter asiaticus’, ‘Ca. L. americanus’, and ‘Ca. L. africanus’. Recent findings suggested that some Japanese strains lack the bacteriophage-type DNA polymerase region (DNA pol), in contrast to the Floridian psy62 strain. The whole genome sequence of the pol-negative ‘Ca. L. asiaticus’ Japanese isolate Ishi-1 was determined by metagenomic analysis of DNA extracted from ‘Ca. L. asiaticus’-infected psyllids and leaf midribs. The 1.19-Mb genome has an average 36.32% GC content. Annotation revealed 13 operons encoding rRNA and 44 tRNA genes, but no typical bacterial pathogenesis-related genes were located within the genome, similar to the Floridian psy62 and Chinese gxpsy. In contrast to other ‘Ca. L. asiaticus’ strains, the genome of the Japanese Ishi-1 strain lacks a prophage-related region.     

Recognizing the pseudogenes in bacterial genomes

Pseudogenes are now known to be a regular feature of bacterial genomes and are found in particularly high numbers within the genomes of recently emerged bacterial pathogens. As most pseudogenes are recognized by sequence alignments, we use newly available genomic sequences to identify the pseudogenes in 11 genomes from 4 bacterial genera, each of which contains at least 1 human pathogen. The numbers of pseudogenes range from 27 in Staphylococcus aureus MW2 to 337 in Yersinia pestis CO92 (e.g. 1–8% of the annotated genes in the genome). Most pseudogenes are formed by small frameshifting indels, but because stop codons are A + T-rich, the two low-G + C Gram-positive taxa (Streptococcus and Staphylococcus) have relatively high fractions of pseudogenes generated by nonsense mutations when compared with more G + C-rich genomes. Over half of the pseudogenes are produced from genes whose original functions were annotated as ‘hypothetical’ or ‘unknown’; however, several broadly distributed genes involved in nucleotide processing, repair or replication have become pseudogenes in one of the sequenced Vibrio vulnificus genomes. Although many of our comparisons involved closely related strains with broadly overlapping gene inventories, each genome contains a largely unique set of pseudogenes, suggesting that pseudogenes are formed and eliminated relatively rapidly from most bacterial genomes.     

Quantifying bacterial evolution in the wild: A birthday problem for Campylobacter lineages

Measuring molecular evolution in bacteria typically requires estimation of the rate at which nucleotide changes accumulate in strains sampled at different times that share a common ancestor. This approach has been useful for dating ecological and evolutionary events that coincide with the emergence of important lineages, such as outbreak strains and obligate human pathogens. However, in multi-host (niche) transmission scenarios, where the pathogen is essentially an opportunistic environmental organism, sampling is often sporadic and rarely reflects the overall population, particularly when concentrated on clinical isolates. This means that approaches that assume recent common ancestry are not applicable. Here we present a new approach to estimate the molecular clock rate in Campylobacter that draws on the popular probability conundrum known as the ‘birthday problem’. Using large genomic datasets and comparative genomic approaches, we use isolate pairs that share recent common ancestry to estimate the rate of nucleotide change for the population. Identifying synonymous and non-synonymous nucleotide changes, both within and outside of recombined regions of the genome, we quantify clock-like diversification to estimate synonymous rates of nucleotide change for the common pathogenic bacteria Campylobacter coli (2.4 x 10⁻⁶ s/s/y) and Campylobacter jejuni (3.4 x 10⁻⁶ s/s/y). Finally, using estimated total rates of nucleotide change, we infer the number of effective lineages within the sample time frame–analogous to a shared birthday–and assess the rate of turnover of lineages in our sample set over short evolutionary timescales. This provides a generalizable approach to calibrating rates in populations of environmental bacteria and shows that multiple lineages are maintained, implying that large-scale clonal sweeps may take hundreds of years or more in these species.     

Genome sequencing and analysis of two early-flowering cherry (Cerasus × kanzakura) varieties, ‘Kawazu-zakura’ and ‘Atami-zakura’

To gain genetic insights into the early-flowering phenotype of ornamental cherry, also known as sakura, we determined the genome sequences of two early-flowering cherry (Cerasus × kanzakura) varieties, ‘Kawazu-zakura’ and ‘Atami-zakura’. Because the two varieties are interspecific hybrids, likely derived from crosses between Cerasus campanulata (early-flowering species) and Cerasus speciosa, we employed the haplotype-resolved sequence assembly strategy. Genome sequence reads obtained from each variety by single-molecule real-time sequencing (SMRT) were split into two subsets, based on the genome sequence information of the two probable ancestors, and assembled to obtain haplotype-phased genome sequences. The resultant genome assembly of ‘Kawazu-zakura’ spanned 519.8 Mb with 1,544 contigs and an N50 value of 1,220.5 kb, while that of ‘Atami-zakura’ totalled 509.6 Mb with 2,180 contigs and an N50 value of 709.1 kb. A total of 72,702 and 69,528 potential protein-coding genes were predicted in the genome assemblies of ‘Kawazu-zakura’ and ‘Atami-zakura’, respectively. Gene clustering analysis identified 2,634 clusters uniquely presented in the C. campanulata haplotype sequences, which might contribute to its early-flowering phenotype. Genome sequences determined in this study provide fundamental information for elucidating the molecular and genetic mechanisms underlying the early-flowering phenotype of ornamental cherry tree varieties and their relatives.     

Phylogenetic analysis of the <Emphasis Type="Italic">lux</Emphasis> operon distinguishes two evolutionarily distinct clades of <Emphasis Type="Italic">Photobacterium leiognathi</Emphasis>

The luminous marine bacterium Photobacterium mandapamensis was synonymized several years ago with Photobacterium leiognathi based on a high degree of phenotypic and genetic similarity. To test the possibility that P. leiognathi as now formulated, however, actually contains two distinct bacterial groups reflecting the earlier identification of P. mandapamensis and P. leiognathi as separate species, we compared P. leiognathi strains isolated from light-organ symbiosis with leiognathid fishes (i.e., ATCC 25521^T, ATCC 25587, lequu.1.1 and lleuc.1.1) with strains from seawater originally described as P. mandapamensis and later synonymized as P. leiognathi (i.e., ATCC 27561^T and ATCC 33981) and certain strains initially identified as P. leiognathi (i.e., PL-721, PL-741, 554). Analysis of the 16S rRNA and gyrB genes did not resolve distinct clades, affirming a close relationship among these strains. However, strains ATCC 27561^T, ATCC 33981, PL-721, PL-741 and 554 were found to bear a luxF gene in the lux operon (luxABFE), whereas ATCC 25521^T, ATCC 25587, lequu.1.1 and lleuc.1.1 lack this gene (luxABE). Phylogenetic analysis of the luxAB(F)E region confirmed this distinction. Furthermore, ATCC 27561^T, ATCC 33981, PL-721, PL-741 and 554 all produced a higher level of luminescence on high-salt medium, as previously described for PL-721, whereas ATCC 25521^T, ATCC 25587, lequu.1.1 and lleuc.1.1 all produced a higher level of luminescence on low-salt medium, a characteristic of P. leiognathi from leiognathid fish light organs. These results demonstrate that P. leiognathi contains two evolutionarily and phenotypically distinct clades, P. leiognathi subsp. leiognathi (strains ATCC 25521^T, ATCC 25587, lequu.1.1 and lleuc.1.1), and P. leiognathi subsp. mandapamensis (strains ATCC 27561^T, ATCC 33981, PL-721, PL-741 and 554).Electronic Supplementary Material Supplementary material is available for this article at .     

The closest lineage of Archaeplastida is revealed by phylogenomics analyses that include Microheliella maris

By clarifying the phylogenetic positions of ‘orphan’ protists (unicellular micro-eukaryotes with no affinity to extant lineages), we may uncover the novel affiliation between two (or more) major lineages in eukaryotes. Microheliella maris was an orphan protist, which failed to be placed within the previously described lineages by pioneering phylogenetic analyses. In this study, we analysed a 319-gene alignment and demonstrated that M. maris represents a basal lineage of one of the major eukaryotic lineages, Cryptista. We here propose a new clade name ‘Pancryptista’ for Cryptista plus M. maris. The 319-gene analyses also indicated that M. maris is a key taxon to recover the monophyly of Archaeplastida and the sister relationship between Archaeplastida and Pancryptista, which is collectively called ‘CAM clade’ here. Significantly, Cryptophyceae tend to be attracted to Rhodophyta depending on the taxon sampling (ex., in the absence of M. maris and Rhodelphidia) and the particular phylogenetic ‘signal’ most likely hindered the stable recovery of the monophyly of Archaeplastida in previous studies.     

Ecological fitness, genomic islands and bacterial pathogenicity. A Darwinian view of the evolution of microbes

The compositions of bacterial genomes can be changed rapidly and dramatically through a variety of processes including horizontal gene transfer. This form of change is key to bacterial evolution, as it leads to ‘evolution in quantum leaps’. Horizontal gene transfer entails the incorporation of genetic elements transferred from another organism—perhaps in an earlier generation—directly into the genome, where they form ‘genomic islands’, i.e. blocks of DNA with signatures of mobile genetic elements. Genomic islands whose functions increase bacterial fitness, either directly or indirectly, have most likely been positively selected and can be termed ‘fitness islands’. Fitness islands can be divided into several subtypes: ‘ecological islands’ in environmental bacteria and ‘saprophytic islands’, ‘symbiosis islands’ or ‘pathogenicity islands’ (PAIs) in microorganisms that interact with living hosts. Here we discuss ways in which PAIs contribute to the pathogenic potency of bacteria, and the idea that genetic entities similar to genomic islands may also be present in the genomes of eukaryotes.     

galign: A Tool for Rapid Genome Polymorphism Discovery

Background

Highly parallel sequencing technologies have become important tools in the analysis of sequence polymorphisms on a genomic scale. However, the development of customized software to analyze data produced by these methods has lagged behind.

Methods/Principal Findings

Here I describe a tool, ‘galign’, designed to identify polymorphisms between sequence reads obtained using Illumina/Solexa technology and a reference genome. The ‘galign’ alignment tool does not use Smith-Waterman matrices for sequence comparisons. Instead, a simple algorithm comparing parsed sequence reads to parsed reference genome sequences is used. ‘galign’ output is geared towards immediate user application, displaying polymorphism locations, nucleotide changes, and relevant predicted amino-acid changes for ease of information processing. To do so, ‘galign’ requires several accessory files easily derived from an annotated reference genome. Direct sequencing as well as in silico studies demonstrate that ‘galign’ provides lesion predictions comparable in accuracy to available prediction programs, accompanied by greater processing speed and more user-friendly output. We demonstrate the use of ‘galign’ to identify mutations leading to phenotypic consequences in C. elegans.

Conclusion/Significance

Our studies suggest that ‘galign’ is a useful tool for polymorphism discovery, and is of immediate utility for sequence mining in C. elegans.     

Paleovirology of ‘syncytins’, retroviral env genes exapted for a role in placentation

The development of the emerging field of ‘paleovirology’ allows biologists to reconstruct the evolutionary history of fossil endogenous retroviral sequences integrated within the genome of living organisms and has led to the retrieval of conserved, ancient retroviral genes ‘exapted’ by ancestral hosts to fulfil essential physiological roles, syncytin genes being undoubtedly among the most remarkable examples of such a phenomenon. Indeed, syncytins are ‘new’ genes encoding proteins derived from the envelope protein of endogenous retroviral elements that have been captured and domesticated on multiple occasions and independently in diverse mammalian species, through a process of convergent evolution. Knockout of syncytin genes in mice provided evidence for their absolute requirement for placenta development and embryo survival, via formation by cell–cell fusion of syncytial cell layers at the fetal–maternal interface. These genes of exogenous origin, acquired ‘by chance’ and yet still ‘necessary’ to carry out a basic function in placental mammals, may have been pivotal in the emergence of mammalian ancestors with a placenta from egg-laying animals via the capture of a founding retroviral env gene, subsequently replaced in the diverse mammalian lineages by new env-derived syncytin genes, each providing its host with a positive selective advantage.     

Optimization of linkage mapping strategy and construction of a high-density American lotus linkage map

Background

Lotus is a diploid plant with agricultural, medicinal, and ecological significance. Genetic linkage maps are fundamental resources for genome and genetic study, and also provide molecular markers for breeding in agriculturally important species. Genotyping by sequencing revolutionized genetic mapping, the restriction-site associated DNA sequencing (RADseq) allowed rapid discovery of thousands of SNPs markers, and a crucial aspect of the sequence based mapping strategy is the reference sequences used for marker identification.

Results

We assessed the effectiveness of linkage mapping using three types of references for scoring markers: the unmasked genome, repeat masked genome, and gene models. Overall, the repeat masked genome produced the optimal genetic maps. A high-density genetic map of American lotus was constructed using an F₁ population derived from a cross between Nelumbo nucifera ‘China Antique’ and N. lutea ‘AL1’. A total of 4,098 RADseq markers were used to construct the American lotus ‘AL1’ genetic map, and 147 markers were used to construct the Chinese lotus ‘China Antique’ genetic map. The American lotus map has 9 linkage groups, and spans 494.3 cM, with an average distance of 0.7 cM between adjacent markers. The American lotus map was used to anchor scaffold sequences in the N. nucifera ‘China Antique’ draft genome. 3,603 RADseq markers anchored 234 individual scaffold sequences into 9 megascaffolds spanning 67% of the 804 Mb draft genome.

Conclusions

Among the unmasked genome, repeat masked genome and gene models, the optimal reference sequences to call RADseq markers for map construction is repeat masked genome. This high density genetic map is a valuable resource for genomic research and crop improvement in lotus.     

Microarray Comparative Genomic Hybridisation Analysis Incorporating Genomic Organisation,and Application to Enterobacterial Plant Pathogens

Microarray comparative genomic hybridisation (aCGH) provides an estimate of the relative abundance of genomic DNA (gDNA) taken from comparator and reference organisms by hybridisation to a microarray containing probes that represent sequences from the reference organism. The experimental method is used in a number of biological applications, including the detection of human chromosomal aberrations, and in comparative genomic analysis of bacterial strains, but optimisation of the analysis is desirable in each problem domain.We present a method for analysis of bacterial aCGH data that encodes spatial information from the reference genome in a hidden Markov model. This technique is the first such method to be validated in comparisons of sequenced bacteria that diverge at the strain and at the genus level: Pectobacterium atrosepticum SCRI1043 (Pba1043) and Dickeya dadantii 3937 (Dda3937); and Lactococcus lactis subsp. lactis IL1403 and L. lactis subsp. cremoris MG1363. In all cases our method is found to outperform common and widely used aCGH analysis methods that do not incorporate spatial information. This analysis is applied to comparisons between commercially important plant pathogenic soft-rotting enterobacteria (SRE) Pba1043, P. atrosepticum SCRI1039, P. carotovorum 193, and Dda3937.Our analysis indicates that it should not be assumed that hybridisation strength is a reliable proxy for sequence identity in aCGH experiments, and robustly extends the applicability of aCGH to bacterial comparisons at the genus level. Our results in the SRE further provide evidence for a dynamic, plastic ‘accessory’ genome, revealing major genomic islands encoding gene products that provide insight into, and may play a direct role in determining, variation amongst the SRE in terms of their environmental survival, host range and aetiology, such as phytotoxin synthesis, multidrug resistance, and nitrogen fixation.     

Fundamental structural units of the Escherichia coli nucleoid revealed by atomic force microscopy

A small container of several to a few hundred µm³ (i.e. bacterial cells and eukaryotic nuclei) contains extremely long genomic DNA (i.e. mm and m long, respectively) in a highly organized fashion. To understand how such genomic architecture could be achieved, Escherichia coli nucleoids were subjected to structural analyses under atomic force microscopy, and found to change their structure dynamically during cell growth, i.e. the nucleoid structure in the stationary phase was more tightly compacted than in the log phase. However, in both log and stationary phases, a fundamental fibrous structure with a diameter of ~80 nm was found. In addition to this ‘80 nm fiber’, a thinner ‘40 nm fiber’ and a higher order ‘loop’ structure were identified in the log phase nucleoid. In the later growth phases, the nucleoid turned into a ‘coral reef structure’ that also possessed the 80 nm fiber units, and, finally, into a ‘tightly compacted nucleoid’ that was stable in a mild lysis buffer. Mutant analysis demonstrated that these tight compactions of the nucleoid required a protein, Dps. From these results and previously available information, we propose a structural model of the E.coli nucleoid.     

Protein interaction mapping on a functional shotgun sequence of Rickettsia sibirica

Protein interaction maps can reveal novel pathways and functional complexes, allowing ‘guilt by association’ annotation of uncharacterized proteins. To address the need for large-scale protein interaction analyses, a bacterial two-hybrid system was coupled with a whole genome shotgun sequencing approach for microbial genome analysis. We report the first large-scale proteomics study using this system, integrating de novo genome sequencing with functional interaction mapping and annotation in a high-throughput format. We apply the approach by shotgun sequencing and annotating the genome of Rickettsia sibirica strain 246, an obligate intracellular human pathogen among the Spotted Fever Group rickettsiae. The bacteria invade endothelial cells and cause lysis after large amounts of progeny have accumulated. Little is known about specific Rickettsial virulence factors and their mode of pathogenicity. Analysis of the combined genomic sequence and protein–protein interaction data for a set of virulence related Type IV secretion system (T4SS) proteins revealed over 250 interactions and will provide insight into the mechanism of Rickettsial pathogenicity.     

Cassette-like variation of restriction enzyme genes in Escherichia coli C and relatives

A surprising result of comparative bacterial genomics has been the large amount of DNA found to be present in one strain but not in another of the same species. We examine in detail one location where gene content varies extensively, the restriction cluster in Escherichia coli. This region is designated the Immigration Control Region (ICR) for the density and variability of restriction functions found there. To better define the boundaries of this variable locus, we determined the sequence of the region from a restrictionless strain, E.coli C. Here we compare the 13.7 kb E.coli C sequence spanning the site of the ICR with corresponding sequences from five E.coli strains and Salmonella typhimurium LT2. To discuss this variation, we adopt the term ‘framework’ to refer to genes that are stable components of genomes within related lineages, while ‘migratory’ genes are transient inhabitants of the genome. Strikingly, seven different migratory DNA segments, encoding different sets of genes and gene fragments, alternatively occupy a single well-defined location in the seven strains examined. The flanking framework genes, yjiS and yjiA, display approximately normal patterns of conservation. The patterns observed are consistent with the action of a site-specific recombinase. Since no nearby gene codes for a likely recombinase of known families, such a recombinase must be of a new family or unlinked.