Segregate or cooperate- a study of the interaction between two species of Dictyostelium

Background

Plastids have inherited their own genomes from a single cyanobacterial ancestor, but the majority of cyanobacterial genes, once retained in the ancestral plastid genome, have been lost or transferred into the eukaryotic host nuclear genome via endosymbiotic gene transfer. Although previous studies showed that cyanobacterial gnd genes, which encode 6-phosphogluconate dehydrogenase, are present in several plastid-lacking protists as well as primary and secondary plastid-containing phototrophic eukaryotes, the evolutionary paths of these genes remain elusive.

Results

Here we show an extended phylogenetic analysis including novel gnd gene sequences from Excavata and Glaucophyta. Our analysis demonstrated the patchy distribution of the excavate genes in the gnd gene phylogeny. The Diplonema gene was related to cytosol-type genes in red algae and Opisthokonta, while heterolobosean genes occupied basal phylogenetic positions with plastid-type red algal genes within the monophyletic eukaryotic group that is sister to cyanobacterial genes. Statistical tests based on exhaustive maximum likelihood analyses strongly rejected that heterolobosean gnd genes were derived from a secondary plastid of green lineage. In addition, the cyanobacterial gnd genes from phototrophic and phagotrophic species in Euglenida were robustly monophyletic with Stramenopiles, and this monophyletic clade was moderately separated from those of red algae. These data suggest that these secondary phototrophic groups might have acquired the cyanobacterial genes independently of secondary endosymbioses.

Conclusion

We propose an evolutionary scenario in which plastid-lacking Excavata acquired cyanobacterial gnd genes via eukaryote-to-eukaryote lateral gene transfer or primary endosymbiotic gene transfer early in eukaryotic evolution, and then lost either their pre-existing or cyanobacterial gene.     

Phylogenomic evidence supports past endosymbiosis, intracellular and horizontal gene transfer in Cryptosporidium parvum

Background

The apicomplexan parasite Cryptosporidium parvum is an emerging pathogen capable of causing illness in humans and other animals and death in immunocompromised individuals. No effective treatment is available and the genome sequence has recently been completed. This parasite differs from other apicomplexans in its lack of a plastid organelle, the apicoplast. Gene transfer, either intracellular from an endosymbiont/donor organelle or horizontal from another organism, can provide evidence of a previous endosymbiotic relationship and/or alter the genetic repertoire of the host organism. Given the importance of gene transfers in eukaryotic evolution and the potential implications for chemotherapy, it is important to identify the complement of transferred genes in Cryptosporidium.

Results

We have identified 31 genes of likely plastid/endosymbiont (n = 7) or prokaryotic (n = 24) origin using a phylogenomic approach. The findings support the hypothesis that Cryptosporidium evolved from a plastid-containing lineage and subsequently lost its apicoplast during evolution. Expression analyses of candidate genes of algal and eubacterial origin show that these genes are expressed and developmentally regulated during the life cycle of C. parvum.

Conclusions

Cryptosporidium is the recipient of a large number of transferred genes, many of which are not shared by other apicomplexan parasites. Genes transferred from distant phylogenetic sources, such as eubacteria, may be potential targets for therapeutic drugs owing to their phylogenetic distance or the lack of homologs in the host. The successful integration and expression of the transferred genes in this genome has changed the genetic and metabolic repertoire of the parasite.     

Virulence determinants, drug resistance and mobile genetic elements of Laribacter hongkongensis: a genome-wide analysis

Background

Laribacter hongkongensis is associated with community-acquired gastroenteritis and traveler's diarrhea. In this study, we performed an in-depth annotation of the genes in its genome related to the various steps in the infective process, drug resistance and mobile genetic elements.

Results

For acid and bile resistance, L. hongkongensis possessed a urease gene cassette, two arc gene clusters and bile salt efflux systems. For intestinal colonization, it possessed a putative adhesin of the autotransporter family homologous to those of diffusely adherent Escherichia coli (E. coli) and enterotoxigenic E. coli. To evade from host defense, it possessed superoxide dismutase and catalases. For lipopolysaccharide biosynthesis, it possessed the same set of genes that encode enzymes for synthesizing lipid A, two Kdo units and heptose units as E. coli, but different genes for its symmetrical acylation pattern, and nine genes for polysaccharide side chains biosynthesis. It contained a number of CDSs that encode putative cell surface acting (RTX toxin and hemolysins) and intracellular cytotoxins (patatin-like proteins) and enzymes for invasion (outer membrane phospholipase A). It contained a broad variety of antibiotic resistance-related genes, including genes related to β-lactam (n = 10) and multidrug efflux (n = 54). It also contained eight prophages, 17 other phage-related CDSs and 26 CDSs for transposases.

Conclusions

The L. hongkongensis genome possessed genes for acid and bile resistance, intestinal mucosa colonization, evasion of host defense and cytotoxicity and invasion. A broad variety of antibiotic resistance or multidrug resistance genes, a high number of prophages, other phage-related CDSs and CDSs for transposases, were also identified.     

The complexity of gene expression dynamics revealed by permutation entropy

Background

High complexity is considered a hallmark of living systems. Here we investigate the complexity of temporal gene expression patterns using the concept of Permutation Entropy (PE) first introduced in dynamical systems theory. The analysis of gene expression data has so far focused primarily on the identification of differentially expressed genes, or on the elucidation of pathway and regulatory relationships. We aim to study gene expression time series data from the viewpoint of complexity.

Results

Applying the PE complexity metric to abiotic stress response time series data in Arabidopsis thaliana, genes involved in stress response and signaling were found to be associated with the highest complexity not only under stress, but surprisingly, also under reference, non-stress conditions. Genes with house-keeping functions exhibited lower PE complexity. Compared to reference conditions, the PE of temporal gene expression patterns generally increased upon stress exposure. High-complexity genes were found to have longer upstream intergenic regions and more cis-regulatory motifs in their promoter regions indicative of a more complex regulatory apparatus needed to orchestrate their expression, and to be associated with higher correlation network connectivity degree. Arabidopsis genes also present in other plant species were observed to exhibit decreased PE complexity compared to Arabidopsis specific genes.

Conclusions

We show that Permutation Entropy is a simple yet robust and powerful approach to identify temporal gene expression profiles of varying complexity that is equally applicable to other types of molecular profile data.     

Functional diversity of human protein kinase splice variants marks significant expansion of human kinome

Background

The low pH environment of the human stomach is lethal for most microorganisms; but not Escherichia coli, which can tolerate extreme acid stress. Acid resistance in E. coli is hierarchically controlled by numerous regulators among which are small noncoding RNAs (sncRNA).

Results

In this study, we individually deleted seventy-nine sncRNA genes from the E. coli K12-MG1655 chromosome, and established a single-sncRNA gene knockout library. By systematically screening the sncRNA mutant library, we show that the sncRNA GcvB is a novel regulator of acid resistance in E. coli. We demonstrate that GcvB enhances the ability of E. coli to survive low pH by upregulating the levels of the alternate sigma factor RpoS.

Conclusion

GcvB positively regulates acid resistance by affecting RpoS expression. These data advance our understanding of the sncRNA regulatory network involved in modulating acid resistance in E. coli.     

A set of novel multiplex Taqman real-time PCRs for the detection of diarrhoeagenic Escherichia coli and its use in determining the prevalence of EPEC and EAEC in a university hospital

Background

Escherichia coli isolates of equine faecal origin were investigated for antibiotic resistance, resistance genes and their ability to perform horizontal transfer.

Methods

In total, 264 faecal samples were collected from 138 horses in hospital and community livery premises in northwest England, yielding 296 resistant E. coli isolates. Isolates were tested for susceptibility to antimicrobial drugs by disc diffusion and agar dilution methods in order to determine minimum inhibitory concentrations (MIC). PCR amplification was used to detect genes conferring resistance to: ampicillin (TEM and SHV beta-lactamase), chloramphenicol (catI, catII, catIII and cml), tetracycline (tetA, tetB, tetC, tetD, tet E and tetG), and trimethoprim (dfrA1, dfrA9, dfrA12, dfrA13, dfr7, and dfr17).

Results

The proportion of antibiotic resistant isolates, and multidrug resistant isolates (MDR) was significantly higher in hospital samples compared to livery samples (MDR: 48% of hospital isolates; 12% of livery isolates, p < 0.001). Resistance to ciprofloxacin and florfenicol were identified mostly within the MDR phenotypes. Resistance genes included dfr, TEM beta-lactamase, tet and cat, conferring resistance to trimethoprim, ampicillin, tetracycline and chloramphenicol, respectively. Within each antimicrobial resistance group, these genes occurred at frequencies of 93% (260/279), 91%, 86.8% and 73.5%, respectively; with 115/296 (38.8%) found to be MDR isolates. Conjugation experiments were performed on selected isolates and MDR phenotypes were readily transferred.

Conclusions

Our findings demonstrate that E. coli of equine faecal origin are commonly resistant to antibiotics used in human and veterinary medicine. Furthermore, our results suggest that most antibiotic resistance observed in equine E. coli is encoded by well-known and well-characterized resistant genes common to E. coli from man and domestic animals. These data support the ongoing concern about antimicrobial resistance, MDR, antimicrobial use in veterinary medicine and the zoonotic risk that horses could potentially pose to public health.     

Protein Connectivity and Protein Complexity Promotes Human Gene Duplicability in a Mutually Exclusive Manner

It has previously been reported that protein complexity (i.e. number of subunits in a protein complex) is negatively correlated to gene duplicability in yeast as well as in humans. However, unlike in yeast, protein connectivity in a protein–protein interaction network has a positive correlation with gene duplicability in human genes. In the present study, we have analyzed 1732 human and 1269 yeast proteins that are present both in a protein–protein interaction network as well as in a protein complex network. In the human case, we observed that both protein connectivity and protein complexity complement each other in a mutually exclusive manner over gene duplicability in a positive direction. Analysis of human haploinsufficient proteins and large protein complexes (complex size >10) shows that when protein connectivity does not have any direct association with gene duplicability, there exists a positive correlation between gene duplicability and protein complexity. The same trend, however, is not found in case of yeast, where both protein connectivity and protein complexity independently guide gene duplicability in the negative direction. We conclude that the higher rate of duplication of human genes may be attributed to organismal complexity either by increasing connectivity in the protein–protein interaction network or by increasing protein complexity.     

Identification and characterization of insect-specific proteins by genome data analysis

Background

Previous studies have sought to identify a link between the distribution of variable genes amongst isolates of Campylobacter jejuni and particular host preferences. The genomic sequence data available currently was obtained using only isolates from human or chicken hosts. In order to identify variable genes present in isolates from alternative host species, five subtractions between C. jejuni isolates from different sources (rabbit, cattle, wild bird) were carried out, designed to assess genomic variability within and between common multilocus sequence type (MLST) clonal complexes (ST-21, ST-42, ST-45 and ST-61).

Results

The vast majority (97%) of the 195 subtracted sequences identified had a best BLASTX match with a Campylobacter protein. However, there was considerable variation within and between the four clonal complexes included in the subtractions. The distributions of eight variable sequences, including four with putative roles in the use of alternative terminal electron acceptors, amongst a panel of C. jejuni isolates representing diverse sources and STs, were determined.

Conclusion

There was a clear correlation between clonal complex and the distribution of the metabolic genes. In contrast, there was no evidence to support the hypothesis that the distribution of such genes may be related to host preference. The other variable genes studied were also generally distributed according to MLST type. Thus, we found little evidence for widespread horizontal gene transfer between clonal complexes involving these genes.     

Prevalence of antimicrobial resistance in bacteria isolated from central nervous system specimens as reported by U.S. hospital laboratories from 2000 to 2002

Background

Class 1 integrons contain genetic elements for site-specific recombination, capture and mobilization of resistance genes. Studies investigating the prevalence, distribution and types of integron located resistance genes are important for surveillance of antimicrobial resistance and to understand resistance development at the molecular level.

Methods

We determined the prevalence and genetic content of class 1 integrons in Enterobacteriaceae (strain collection 1, n = 192) and E. coli (strain collection 2, n = 53) from bloodstream infections in patients from six Norwegian hospitals by molecular techniques. Class 1 integrons were also characterized in 54 randomly selected multiresistant E. coli isolates from gastrointestinal human infections (strain collection 3).

Results

Class 1 integrons were present in 10.9% of the Enterobacteriaceae blood culture isolates of collection 1, all but one (S. Typhi) being E. coli. Data indicated variations in class 1 integron prevalence between hospitals. Class 1 integrons were present in 37% and 34% of the resistant blood culture isolates (collection 1 and 2, respectively) and in 42% of the resistant gastrointestinal E. coli. We detected a total of 10 distinct integron cassette PCR amplicons that varied in size between 0.15 kb and 2.2 kb and contained between zero and three resistance genes. Cassettes encoding resistance to trimethoprim and aminoglycosides were most common. We identified and characterized a novel plasmid-located integron with a cassette-bound novel gene (linF) located downstream of an aadA2 gene cassette. The linF gene encoded a putative 273 aa lincosamide nucleotidyltransferase resistance protein and conferred resistance to lincomycin and clindamycin. The deduced LinF amino acid sequence displayed approximately 35% identity to the Enterococcus faecium and Enterococcus faecalis nucleotidyl transferases encoded by linB and linB'

Conclusions

The present study demonstrated an overall low and stable prevalence of class 1 integron gene cassettes in clinical Enterobacteriaceae and E. coli isolates in Norway. Characterization of the novel lincosamide resistance gene extends the growing list of class 1 integron gene cassettes that confer resistance to an increasing number of antibiotics.     

Metabolic flux balance analysis and the in silicoanalysis of Escherichia coliK-12 gene deletions

Background

Genome sequencing and bioinformatics are producing detailed lists of the molecular components contained in many prokaryotic organisms. From this 'parts catalogue' of a microbial cell, in silicorepresentations of integrated metabolic functions can be constructed and analyzed using flux balance analysis (FBA). FBA is particularly well-suited to study metabolic networks based on genomic, biochemical, and strain specific information.

Results

Herein, we have utilized FBA to interpret and analyze the metabolic capabilities of Escherichia coli. We have computationally mapped the metabolic capabilities of E. coliusing FBA and examined the optimal utilization of the E. colimetabolic pathways as a function of environmental variables. We have used an in silicoanalysis to identify seven gene products of central metabolism (glycolysis, pentose phosphate pathway, TCA cycle, electron transport system) essential for aerobic growth of E. colion glucose minimal media, and 15 gene products essential for anaerobic growth on glucose minimal media. The in silico tpi ^-, zwf, and pta ^-mutant strains were examined in more detail by mapping the capabilities of these in silicoisogenic strains.

Conclusions

We found that computational models of E. colimetabolism based on physicochemical constraints can be used to interpret mutant behavior. These in silicaresults lead to a further understanding of the complex genotype-phenotype relation. Supplementary information: 10.1186/1471-2105-1-1     

Relationship between gene duplicability and diversifiability in the topology of biochemical networks

Background

Selective gene duplicability, the extensive expansion of a small number of gene families, is universal. Quantitatively, the number of genes (P_(K)) with K duplicates in a genome decreases precipitously as K increases, and often follows a power law (P_(k)∝k^-α). Functional diversification, either neo- or sub-functionalization, is a major evolution route for duplicate genes.

Results

Using three lines of genomic datasets, we studied the relationship between gene duplicability and diversifiability in the topology of biochemical networks. First, we explored scenario where two pathways in the biochemical networks antagonize each other. Synthetic knockout of respective genes for the two pathways rescues the phenotypic defects of each individual knockout. We identified duplicate gene pairs with sufficient divergences that represent this antagonism relationship in the yeast S. cerevisiae. Such pairs overwhelmingly belong to large gene families, thus tend to have high duplicability. Second, we used distances between proteins of duplicate genes in the protein interaction network as a metric of their diversification. The higher a gene’s duplicate count, the further the proteins of this gene and its duplicates drift away from one another in the networks, which is especially true for genetically antagonizing duplicate genes. Third, we computed a sequence-homology-based clustering coefficient to quantify sequence diversifiability among duplicate genes – the lower the coefficient, the more the sequences have diverged. Duplicate count (K) of a gene is negatively correlated to the clustering coefficient of its duplicates, suggesting that gene duplicability is related to the extent of sequence divergence within the duplicate gene family.

Conclusion

Thus, a positive correlation exists between gene diversifiability and duplicability in the context of biochemical networks – an improvement of our understanding of gene duplicability.     

Unraveling the genomic mosaic of a ubiquitous genus of marine cyanobacteria

Background

The picocyanobacterial genus Synechococcus occurs over wide oceanic expanses, having colonized most available niches in the photic zone. Large scale distribution patterns of the different Synechococcus clades (based on 16S rRNA gene markers) suggest the occurrence of two major lifestyles ('opportunists'/'specialists'), corresponding to two distinct broad habitats ('coastal'/'open ocean'). Yet, the genetic basis of niche partitioning is still poorly understood in this ecologically important group.

Results

Here, we compare the genomes of 11 marine Synechococcus isolates, representing 10 distinct lineages. Phylogenies inferred from the core genome allowed us to refine the taxonomic relationships between clades by revealing a clear dichotomy within the main subcluster, reminiscent of the two aforementioned lifestyles. Genome size is strongly correlated with the cumulative lengths of hypervariable regions (or 'islands'). One of these, encompassing most genes encoding the light-harvesting phycobilisome rod complexes, is involved in adaptation to changes in light quality and has clearly been transferred between members of different Synechococcus lineages. Furthermore, we observed that two strains (RS9917 and WH5701) that have similar pigmentation and physiology have an unusually high number of genes in common, given their phylogenetic distance.

Conclusion

We propose that while members of a given marine Synechococcus lineage may have the same broad geographical distribution, local niche occupancy is facilitated by lateral gene transfers, a process in which genomic islands play a key role as a repository for transferred genes. Our work also highlights the need for developing picocyanobacterial systematics based on genome-derived parameters combined with ecological and physiological data.     

Validation of suitable house keeping genes for hypoxia-cultured human chondrocytes

Background

Sleep is a restorative process and is essential for maintenance of mental and physical health. In an attempt to understand the complexity of sleep, multidisciplinary strategies, including genetic approaches, have been applied to sleep research. Although quantitative real time PCR has been used in previous sleep-related gene expression studies, proper validation of reference genes is currently lacking. Thus, we examined the effect of total or paradoxical sleep deprivation (TSD or PSD) on the expression stability of the following frequently used reference genes in brain and blood: beta-actin (b-actin), beta-2-microglobulin (B2M), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), and hypoxanthine guanine phosphoribosyl transferase (HPRT).

Results

Neither TSD nor PSD affected the expression stability of all tested genes in both tissues indicating that b-actin, B2M, GAPDH and HPRT are appropriate reference genes for the sleep-related gene expression studies. In order to further verify these results, the relative expression of brain derived neurotrophic factor (BDNF) and glycerol-3-phosphate dehydrogenase1 (GPD1) was evaluated in brain and blood, respectively. The normalization with each of four reference genes produced similar pattern of expression in control and sleep deprived rats, but subtle differences in the magnitude of expression fold change were observed which might affect the statistical significance.

Conclusion

This study demonstrated that sleep deprivation does not alter the expression stability of commonly used reference genes in brain and blood. Nonetheless, the use of multiple reference genes in quantitative RT-PCR is required for the accurate results.     

The influence of gene-environment interactions on GHR and IGF-1 expression and their association with growth in brook charr, Salvelinus fontinalis (Mitchill)

Background

Information transfer systems in Archaea, including many components of the DNA replication machinery, are similar to those found in eukaryotes. Functional assignments of archaeal DNA replication genes have been primarily based upon sequence homology and biochemical studies of replisome components, but few genetic studies have been conducted thus far. We have developed a tractable genetic system for knockout analysis of genes in the model halophilic archaeon, Halobacterium sp. NRC-1, and used it to determine which DNA replication genes are essential.

Results

Using a directed in-frame gene knockout method in Halobacterium sp. NRC-1, we examined nineteen genes predicted to be involved in DNA replication. Preliminary bioinformatic analysis of the large haloarchaeal Orc/Cdc6 family, related to eukaryotic Orc1 and Cdc6, showed five distinct clades of Orc/Cdc6 proteins conserved in all sequenced haloarchaea. Of ten orc/cdc6 genes in Halobacterium sp. NRC-1, only two were found to be essential, orc10, on the large chromosome, and orc2, on the minichromosome, pNRC200. Of the three replicative-type DNA polymerase genes, two were essential: the chromosomally encoded B family, polB1, and the chromosomally encoded euryarchaeal-specific D family, polD1/D2 (formerly called polA1/polA2 in the Halobacterium sp. NRC-1 genome sequence). The pNRC200-encoded B family polymerase, polB2, was non-essential. Accessory genes for DNA replication initiation and elongation factors, including the putative replicative helicase, mcm, the eukaryotic-type DNA primase, pri1/pri2, the DNA polymerase sliding clamp, pcn, and the flap endonuclease, rad2, were all essential. Targeted genes were classified as non-essential if knockouts were obtained and essential based on statistical analysis and/or by demonstrating the inability to isolate chromosomal knockouts except in the presence of a complementing plasmid copy of the gene.

Conclusion

The results showed that ten out of nineteen eukaryotic-type DNA replication genes are essential for Halobacterium sp. NRC-1, consistent with their requirement for DNA replication. The essential genes code for two of ten Orc/Cdc6 proteins, two out of three DNA polymerases, the MCM helicase, two DNA primase subunits, the DNA polymerase sliding clamp, and the flap endonuclease.     

Extensive genomic diversity and selective conservation of virulence-determinants in enterohemorrhagic Escherichia coli strains of O157 and non-O157 serotypes

Background

Enterohemorrhagic Escherichia coli (EHEC) O157 causes severe food-borne illness in humans. The chromosome of O157 consists of 4.1 Mb backbone sequences shared by benign E. coli K-12, and 1.4 Mb O157-specific sequences encoding many virulence determinants, such as Shiga toxin genes (stx genes) and the locus of enterocyte effacement (LEE). Non-O157 EHECs belonging to distinct clonal lineages from O157 also cause similar illness in humans. According to the 'parallel' evolution model, they have independently acquired the major virulence determinants, the stx genes and LEE. However, the genomic differences between O157 and non-O157 EHECs have not yet been systematically analyzed.

Results

Using microarray and whole genome PCR scanning analyses, we performed a whole genome comparison of 20 EHEC strains of O26, O111, and O103 serotypes with O157. In non-O157 EHEC strains, although genome sizes were similar with or rather larger than O157 and the backbone regions were well conserved, O157-specific regions were very poorly conserved. Around only 20% of the O157-specific genes were fully conserved in each non-O157 serotype. However, the non-O157 EHECs contained a significant number of virulence genes that are found on prophages and plasmids in O157, and also multiple prophages similar to, but significantly divergent from, those in O157.

Conclusion

Although O157 and non-O157 EHECs have independently acquired a huge amount of serotype- or strain-specific genes by lateral gene transfer, they share an unexpectedly large number of virulence genes. Independent infections of similar but distinct bacteriophages carrying these virulence determinants are deeply involved in the evolution of O157 and non-O157 EHECs.     

A positive feedback-based gene circuit to increase the production of a membrane protein

Background

Membrane proteins are an important class of proteins, playing a key role in many biological processes, and are a promising target in pharmaceutical development. However, membrane proteins are often difficult to produce in large quantities for the purpose of crystallographic or biochemical analyses.

Results

In this paper, we demonstrate that synthetic gene circuits designed specifically to overexpress certain genes can be applied to manipulate the expression kinetics of a model membrane protein, cytochrome bd quinol oxidase in E. coli, resulting in increased expression rates. The synthetic circuit involved is an engineered, autoinducer-independent variant of the lux operon activator LuxR from V. fischeri in an autoregulatory, positive feedback configuration.

Conclusions

Our proof-of-concept experiments indicate a statistically significant increase in the rate of production of the bd oxidase membrane protein. Synthetic gene networks provide a feasible solution for the problem of membrane protein production.