Comparative genomic indexing reveals the phylogenomics of Escherichia coli pathogens

The Escherichia coli O26 serogroup includes important food-borne pathogens associated with human and animal diarrheal disease. Current typing methods have revealed great genetic heterogeneity within the O26 group; the data are often inconsistent and focus only on verotoxin (VT)-positive O26 isolates. To improve current understanding of diversity within this serogroup, the genomic relatedness of VT-positive and -negative O26 strains was assessed by comparative genomic indexing. Our results clearly demonstrate that irrespective of virulence characteristics and pathotype designation, the O26 strains show greater genomic similarity to each other than to any other strain included in this study. Our data suggest that enteropathogenic and VT-expressing E. coli O26 strains represent the same clonal lineage and that W-expressing E. coli O26 strains have gained additional virulence characteristics. Using this approach, we established the core genes which are central to the E. coli species and identified regions of variation from the E. coli K-12 chromosomal backbone.

Epidemic multidrug-resistant (MDR-AmpC) Salmonella enterica serovar Newport strains contain three phage regions and a MDR resistance plasmid

Multidrug-resistant (MDR-AmpC) Salmonella enterica serovar Newport has caused serious disease in animals and humans in North America, whereas in the UK S. enterica serovar Newport is not associated with severe disease and usually sensitive to antibiotics; MDR S. Newport (not AmpC) strains have only been isolated from poultry. We found that UK poultry strains belonged to MLST type ST166 and were distinct from cattle isolates for being able to utilize D-tagotose and when compared by pulsed-field gel electrophoresis (PFGE), comparative genomic hybridization (CGH) and diversity arrays technology (DArT). Cattle strains belonged to the ST45 complex differing from ST166 at all seven loci. PFGE showed that 19 out of 27 cattle isolates were more than 85% similar to each other and some UK and US strains were indistinguishable. Both CGH and DArT identified genes (including phage-related ones) that were uniquely present in the US isolates and two such genes identified by DArT showed sequence similarities with the pertussis-like (artAB) toxin. This work demonstrates that MDR-AmpC S. Newport from the USA are genetically closely related to pan-susceptible strains from the UK, but contained three extra phage regions and a MDR plasmid.

Identification of genetic and phenotypic differences associated with prevalent and non-prevalent Salmonella Enteritidis phage types: analysis of variation in amino acid transport

In this study, differences at the genetic level of 37 Salmonella Enteritidis strains from five phage types (PTs) were compared using comparative genomic hybridization (CGH) to assess differences between PTs. There were approximately 400 genes that differentiated prevalent (4, 6, 8 and 13a) and sporadic (11) PTs, of which 35 were unique to prevalent PTs, including six plasmid-borne genes, pefA, B, C, D, srgC and rck, and four chromosomal genes encoding putative amino acid transporters. Phenotype array studies also demonstrated that strains from prevalent PTs were less susceptible to urea stress and utilized L-histidine, L-glutamine, L-proline, L-aspartic acid, gly-asn and gly-gln more efficiently than PT11 strains. Complementation of a PT11 strain with the transporter genes from PT4 resulted in a significant increase in utilization of the amino acids and reduced susceptibility to urea stress. In epithelial cell association assays, PT11 strains were less invasive than other prevalent PTs. Most strains from prevalent PTs were better biofilm formers at 37 degrees C than at 28 degrees C, whilst the converse was true for PT11 strains. Collectively, the results indicate that genetic and corresponding phenotypic differences exist between strains of the prevalent PTs 4, 6, 8 and 13a and non-prevalent PT11 strains that are likely to provide a selective advantage for strains from the former PTs and could help them to enter the food chain and cause salmonellosis.

Comparative genome analysis provides insights into the evolution and adaptation of Pseudomonas syringae pv. aesculi on Aesculus hippocastanum.

A recently emerging bleeding canker disease, caused by Pseudomonas syringae pathovar aesculi (Pae), is threatening European horse chestnut in northwest Europe. Very little is known about the origin and biology of this new disease. We used the nucleotide sequences of seven commonly used marker genes to investigate the phylogeny of three strains isolated recently from bleeding stem cankers on European horse chestnut in Britain (E-Pae). On the basis of these sequences alone, the E-Pae strains were identical to the Pae type-strain (I-Pae), isolated from leaf spots on Indian horse chestnut in India in 1969. The phylogenetic analyses also showed that Pae belongs to a distinct clade of P. syringae pathovars adapted to woody hosts. We generated genome-wide Illumina sequence data from the three E-Pae strains and one strain of I-Pae. Comparative genomic analyses revealed pathovar-specific genomic regions in Pae potentially implicated in virulence on a tree host, including genes for the catabolism of plant-derived aromatic compounds and enterobactin synthesis. Several gene clusters displayed intra-pathovar variation, including those encoding type IV secretion, a novel fatty acid biosynthesis pathway and a sucrose uptake pathway. Rates of single nucleotide polymorphisms in the four Pae genomes indicate that the three E-Pae strains diverged from each other much more recently than they diverged from I-Pae. The very low genetic diversity among the three geographically distinct E-Pae strains suggests that they originate from a single, recent introduction into Britain, thus highlighting the serious environmental risks posed by the spread of an exotic plant pathogenic bacterium to a new geographic location. The genomic regions in Pae that are absent from other P. syringae pathovars that infect herbaceous hosts may represent candidate genetic adaptations to infection of the woody parts of the tree.

Genetic diversity among Escherichia coli O157 : H7 isolates and identification of genes linked to human infections

An Escherichia coli oligonucleotide microarray based on three sequenced genomes was validated for comparative genomic microarray hybridization and used to study the diversity of E. coli O157 isolates from human infections and food and animal sources. Among 26 test strains, 24 (including both Shiga toxin [Stx]-positive and -negative strains) were found to be related to the two sequenced E. coli O157:117 strains, EDL933 and Sakai. However, these strains showed much greater genetic diversity than those reported previously, and most of them could not be categorized as either lineage I or H. Some genes were found more often in isolates from human than from nonhuman sources; e.g., ECs1202 and ECs2976, associated with stx2AB and stx1AB, were in all isolates from human sources but in only 40% of those from nonhuman sources. Some (but not all) lineage I-specific or -dominant genes were also more frequently associated with isolates from human. The results suggested that it might be more effective to concentrate our efforts on finding markers that are directly related to infection rather than those specific to certain lineages. In addition, two Stx-negative O157 cattle isolates (one confirmed to be 117) were significantly different from other Stx-positive and -negative E. coli O157:117 strains and were more similar to MG1655 in their gene content. This work demonstrates that not all E. coli O157:117 strains belong to the same clonal group, and those that were similar to E. coli K-12 might be less virulent.

A process for analysis of microarray comparative genomics hybridisation studies for bacterial genomes

Background: Microarray based comparative genomic hybridisation (CGH) experiments have been used to study numerous biological problems including understanding genome plasticity in pathogenic bacteria. Typically such experiments produce large data sets that are difficult for biologists to handle. Although there are some programmes available for interpretation of bacterial transcriptomics data and CGH microarray data for looking at genetic stability in oncogenes, there are none specifically to understand the mosaic nature of bacterial genomes. Consequently a bottle neck still persists in accurate processing and mathematical analysis of these data. To address this shortfall we have produced a simple and robust CGH microarray data analysis process that may be automated in the future to understand bacterial genomic diversity. Results: The process involves five steps: cleaning, normalisation, estimating gene presence and absence or divergence, validation, and analysis of data from test against three reference strains simultaneously. Each stage of the process is described and we have compared a number of methods available for characterising bacterial genomic diversity, for calculating the cut-off between gene presence and absence or divergence, and shown that a simple dynamic approach using a kernel density estimator performed better than both established, as well as a more sophisticated mixture modelling technique. We have also shown that current methods commonly used for CGH microarray analysis in tumour and cancer cell lines are not appropriate for analysing our data. Conclusion: After carrying out the analysis and validation for three sequenced Escherichia coli strains, CGH microarray data from 19 E. coli O157 pathogenic test strains were used to demonstrate the benefits of applying this simple and robust process to CGH microarray studies using bacterial genomes.

Comparative genomics suggests that an ancestral polyploidy event leads to enhanced root nodule symbiosis in the Papilionoideae

Root nodule symbiosis (RNS) is one of the most efficient biological systems for nitrogen fixation and it occurs in 90% of genera in the Papilionoideae, the largest subfamily of legumes. Most papilionoid species show evidence of a polyploidy event occurred approximately 58 million years ago. Although polyploidy is considered to be an important evolutionary force in plants, the role of this papilionoid polyploidy event, especially its association with RNS, is not understood. In this study, we explored this role using an integrated comparative genomic approach and conducted gene expression comparisons and gene ontology enrichment analyses. The results show the following: (1) approximately a quarter of the papilionoid-polyploidy-derived duplicate genes are retained; (2) there is a striking divergence in the level of expression of gene duplicate pairs derived from the polyploidy event; and (3) the retained duplicates are frequently involved in the processes crucial for RNS establishment, such as symbiotic signalling, nodule organogenesis, rhizobial infection and nutrient exchange and transport. Thus, we conclude that the papilionoid polyploidy event might have further refined RNS and induced a more robust and enhanced symbiotic system. This conclusion partly explains the widespread occurrence of the Papilionoideae.

Origin of avian genome size and structure in non-avian dinosaurs

Avian genomes are small and streamlined compared with those of other amniotes by virtue of having fewer repetitive elements and less non-coding DNA(1,2). This condition has been suggested to represent a key adaptation for flight in birds, by reducing the metabolic costs associated with having large genome and cell sizes(3,4). However, the evolution of genome architecture in birds, or any other lineage, is difficult to study because genomic information is often absent for long-extinct relatives. Here we use a novel bayesian comparative method to show that bone-cell size correlates well with genome size in extant vertebrates, and hence use this relationship to estimate the genome sizes of 31 species of extinct dinosaur, including several species of extinct birds. Our results indicate that the small genomes typically associated with avian flight evolved in the saurischian dinosaur lineage between 230 and 250 million years ago, long before this lineage gave rise to the first birds. By comparison, ornithischian dinosaurs are inferred to have had much larger genomes, which were probably typical for ancestral Dinosauria. Using comparative genomic data, we estimate that genome-wide interspersed mobile elements, a class of repetitive DNA, comprised 5 - 12% of the total genome size in the saurischian dinosaur lineage, but was 7 - 19% of total genome size in ornithischian dinosaurs, suggesting that repetitive elements became less active in the saurischian lineage. These genomic characteristics should be added to the list of attributes previously considered avian but now thought to have arisen in non-avian dinosaurs, such as feathers(5), pulmonary innovations 6, and parental care and nesting

The chromosomal assessment of salt tolerant substituted tritipyrum using genomic fluorescent in situ hybridiization (FISH)

Wheat, although moderately tolerant to salt, can not be cultivated in many areas. However, in the triticeae tribe, some of the wild wheat relatives are highly tolerant, e.g. Thinopyrum bessarabicum, which grows on the sea shore. Eight primary hexaploid tritipyrum lines, amphiploids between Triticum durum and Thinopyrum bessarabicum have been produced which can set seed in at least 250 mM NaCl. These tritipyrums (2n=6x=42, AABBEbEb) due to reasons such as brittle rachis, continuous production of tillers, late maturity, tall stature and meiotic instability will not fulfill the requirements of a successful commercial salt tolerant crop. To overcome such problems the substituted tritipyrum, in which selected Eb chromosomes are replaced by D genome chromosomes of 6x wheat, was produced from 6x tritipyrum x 6x wheat hybrids (F1: 2n=6x=42, AABBDEb) followed by selfing and backcrossing with 6x tritipyrum. The fertile plants among the above progenies were screened by the genomic fluorescent in situ hybridization technique to identify their Eb and D chromosome constitution. This study showed that producing tritiprum with variable numbers of Eb and D genome chromosomes is feasible and that FISH is a useful technique for determining the number of Eb chromosomes present.

Comparative genomics of Brachyspira pilosicoli strains: genome rearrangements, reductions and correlation of genetic compliment with phenotypic diversity.

Background. The anaerobic spirochaete Brachyspira pilosicoli causes enteric disease in avian, porcine and human hosts, amongst others. To date, the only available genome sequence of B. pilosicoli is that of strain 95/1000, a porcine isolate. In the first intra-species genome comparison within the Brachyspira genus, we report the whole genome sequence of B. pilosicoli B2904, an avian isolate, the incomplete genome sequence of B. pilosicoli WesB, a human isolate, and the comparisons with B. pilosicoli 95/1000. We also draw on incomplete genome sequences from three other Brachyspira species. Finally we report the first application of the high-throughput Biolog phenotype screening tool on the B. pilosicoli strains for detailed comparisons between genotype and phenotype. Results. Feature and sequence genome comparisons revealed a high degree of similarity between the three B. pilosicoli strains, although the genomes of B2904 and WesB were larger than that of 95/1000 (~2,765, 2.890 and 2.596 Mb, respectively). Genome rearrangements were observed which correlated largely with the positions of mobile genetic elements. Through comparison of the B2904 and WesB genomes with the 95/1000 genome, features that we propose are non-essential due to their absence from 95/1000 include a peptidase, glycine reductase complex components and transposases. Novel bacteriophages were detected in the newly-sequenced genomes, which appeared to have involvement in intra- and inter-species horizontal gene transfer. Phenotypic differences predicted from genome analysis, such as the lack of genes for glucuronate catabolism in 95/1000, were confirmed by phenotyping. Conclusions. The availability of multiple B. pilosicoli genome sequences has allowed us to demonstrate the substantial genomic variation that exists between these strains, and provides an insight into genetic events that are shaping the species. In addition, phenotype screening allowed determination of how genotypic differences translated to phenotype. Further application of such comparisons will improve understanding of the metabolic capabilities of Brachyspira species.

Comparative genomics of Salmonella enterica serovars Derby and Mbandaka, two prevalent serovars associated with different livestock species in the UK

Background Despite the frequent isolation of Salmonella enterica sub. enterica serovars Derby and Mbandaka from livestock in the UK and USA little is known about the biological processes maintaining their prevalence. Statistics for Salmonella isolations from livestock production in the UK show that S. Derby is most commonly associated with pigs and turkeys and S. Mbandaka with cattle and chickens. Here we compare the first sequenced genomes of S. Derby and S. Mbandaka as a basis for further analysis of the potential host adaptations that contribute to their distinct host species distributions. Results Comparative functional genomics using the RAST annotation system showed that predominantly mechanisms that relate to metabolite utilisation, in vivo and ex vivo persistence and pathogenesis distinguish S. Derby from S. Mbandaka. Alignment of the genome nucleotide sequences of S. Derby D1 and D2 and S. Mbandaka M1 and M2 with Salmonella pathogenicity islands (SPI) identified unique complements of genes associated with host adaptation. We also describe a new genomic island with a putative role in pathogenesis, SPI-23. SPI-23 is present in several S. enterica serovars, including S. Agona, S. Dublin and S. Gallinarum, it is absent in its entirety from S. Mbandaka. Conclusions We discovered a new 37 Kb genomic island, SPI-23, in the chromosome sequence of S. Derby, encoding 42 ORFS, ten of which are putative TTSS effector proteins. We infer from full-genome synonymous SNP analysis that these two serovars diverged, between 182kya and 625kya coinciding with the divergence of domestic pigs. The differences between the genomes of these serovars suggest they have been exposed to different stresses including, phage, transposons and prolonged externalisation. The two serovars possess distinct complements of metabolic genes; many of which cluster into pathways for catabolism of carbon sources.

Comparative analysis of performance and stability among composite cross populations, variety mixtures and pure lines of winter wheat in organic and conventional cropping systems

This study investigated the effects of increased genetic diversity in winter wheat (Triticum aestivum L.), either from hybridization across genotypes or from physical mixing of lines, on grain yield, grain quality, and yield stability in different cropping environments. Sets of pure lines (no diversity), chosen for high yielding ability or high quality, were compared with line mixtures (intermediate level of diversity), and lines crossed with each other in composite cross populations (CCPn, high diversity). Additional populations containing male sterility genes (CCPms) to increase outcrossing rates were also tested. Grain yield, grain protein content, and protein yield were measured at four sites (two organically-managed and two conventionally-managed) over three years, using seed harvested locally in each preceding year. CCPn and mixtures out-yielded the mean of the parents by 2.4% and 3.6%, respectively. These yield differences were consistent across genetic backgrounds but partly inconsistent across cropping environments and years. Yield stability measured by environmental variance was higher in CCPn and CCPms than the mean of the parents. An index of yield reliability tended to be higher in CCPn, CCPms and mixtures than the mean of the parents. Lin and Binns’ superiority values of yield and protein yield were consistently and significantly lower (i.e. better) in the CCPs than in the mean of the parents, but not different between CCPs and mixtures. However, CCPs showed greater early ground cover and plant height than mixtures. When compared with the (locally non-predictable) best-yielding pure line, CCPs and mixtures exhibited lower mean yield and somewhat lower yield reliability but comparable superiority values. Thus, establishing CCPs from smaller sets of high-performing parent lines might optimize their yielding ability. On the whole, the results demonstrate that using increased within-crop genetic diversity can produce wheat crops with improved yield stability and good yield reliability across variable and unpredictable cropping environments.