BMC Genomics - Latest Articles

Relative entropy differences in bacterial chromosomes, plasmids, phages and genomic islands

Jon Bohlin — 2012-02-10T00:00:00Z

Background: We sought to assess whether the concept of relative entropy (information capacity), could aid our understanding of the process of horizontal gene transfer in microbes. We analyzed the differences in information capacity between prokaryotic chromosomes, genomic islands (GI), phages, and plasmids. Relative entropy was estimated using the Kullback-Leibler measure. Results: Relative entropy was highest in bacterial chromosomes and had the sequence chromosomes >GI>phage>plasmid. There was an association between relative entropy and AT content in chromosomes, phages, plasmids and GIs with the strongest association being in phages. Relative entropy was also found to be lower in the obligate intracellular Mycobacterium leprae than in the related M. tuberculosis when measured on a shared set of highly conserved genes. Conclusions: We argue that relative entropy differences reflect how plasmids, phages and GIs interact with microbial host chromosomes and that all these biological entities are, or have been, subjected to different selective pressures. The rate at which amelioration of horizontally acquired DNA occurs within the chromosome is likely to account for the small differences between chromosomes and stably incorporated GIs compared to the transient or independent replicons such as phages and plasmids.

Peptide Markers of Aminoacyl tRNA Synthetases Facilitate Taxa Counting in Metagenomic Data

Erez Persi — 2012-02-10T00:00:00Z

Background: Taxa counting is a major problem faced by analysis of metagenomic data. The most popular method relies on analysis of 16S rRNA sequences, but some studies employ also protein based analyses. It would be advantageous to have a method that is applicable directly to short sequences, of the kind extracted from samples in modern metagenomic research. This is achieved by the technique proposed here. Results: We employ specific peptides, deduced from aminoacyl tRNA synthetases, as markers for the occurrence of single genes in data. Sequences carrying these markers are aligned and compared with each other to provide a lower limit for taxa counts in metagenomic data. The method is compared with 16S rRNA searches on a set of known genomes. The taxa counting problem is analyzed mathematically and a heuristic algorithm is proposed. When applied to genomic contigs of a recent human gut microbiome study, the taxa counting method provides information on numbers of different species and strains. We then apply our method to short read data and demonstrate how it can be calibrated to cope with errors. Comparison to known databases leads to estimates of the percentage of novelties, and the type of phyla involved. Conclusions: A major advantage of our method is its simplicity: it relies on searching sequences for the occurrence of just 4000 specific peptides belonging to the S61 subgroup of aaRS enzymes. When compared to other methods, it provides additional insight into the taxonomic contents of metagenomic data. Furthermore, it can be directly applied to short read data, avoiding the need for genomic contig reconstruction, and taking into account short reads that are otherwise discarded as singletons. Hence it is very suitable for a fast analysis of next generation sequencing data.

Transcriptome profiling of low temperature-treated cassava apical shoots showed dynamic responses of tropical plant to cold stress

Dong An — 2012-02-10T00:00:00Z

Background: Cassava is an important tropical root crop adapted to a wide range of environmental stimuli such as drought and acid soils. Nevertheless, it is an extremely cold-sensitive tropical species. Thus far, there is limited information about gene regulation and signalling pathways related to the cold stress response in cassava. The development of microarray technology has accelerated the study of global transcription profiling under certain conditions. Results: A 60-mer oligonucleotide microarray representing 20,840 genes was used to perform transcriptome profiling in apical shoots of cassava subjected to cold at 7degreesC for 0, 4 and 9 h. A total of 508 transcripts were identified as early cold-responsive genes in which 319 sequences had functional descriptions when aligned with Arabidopsis proteins. Gene ontology annotation analysis identified many cold-relevant categories, including 'Response to abiotic and biotic stimulus', 'Response to stress', 'Transcription factor activity', and 'Chloroplast'. Various stress-associated genes with a wide range of biological functions were found, such as signal transduction components (e.g., MAP kinase 4), transcription factors (TFs, e.g., RAP2.11), and reactive oxygen species (ROS) scavenging enzymes (e.g., catalase 2), as well as photosynthesis-related genes (e.g., PsaL). Seventeen major TF families including many well-studied members (e.g., AP2-EREBP) were also involved in the early response to cold stress. Meanwhile, KEGG pathway analysis uncovered many important pathways, such as 'Plant hormone signal transduction' and 'Starch and sucrose metabolism'. Furthermore, the expression changes of 32 genes under cold and other abiotic stress conditions were validated by real-time RT-PCR. Importantly, most of the tested stress-responsive genes were primarily expressed in mature leaves, stem cambia, and fibrous roots rather than apical buds and young leaves. As a response to cold stress in cassava, an increase in transcripts and enzyme activities of ROS scavenging genes and the accumulation of total soluble sugars (including sucrose and glucose) were also detected. Conclusions: The dynamic expression changes reflect the integrative controlling and transcriptome regulation of the networks in the cold stress response of cassava. The biological processes involved in the signal perception and physiological response might shed light on the molecular mechanisms related to cold tolerance in tropical plants and provide useful candidate genes for genetic improvement.

Global changes in gene expression by the opportunistic pathogen Burkholderia cenocepacia in response to internalization by murine macrophages

Jennifer Tolman — 2012-02-09T00:00:00Z

Background: Burkholderia cenocepacia is an opportunistic pathogen causing life-threatening infections in patients with cystic fibrosis. The bacterium survives within macrophages by interfering with endocytic trafficking and delaying the maturation of the B. cenocepacia-containing phagosome. We hypothesize that B. cenocepacia undergoes changes in gene expression after internalization by macrophages, inducing genes involved in intracellular survival and host adaptation. Results: We examined gene expression by intracellular B. cenocepacia using selective capture of transcribed sequences (SCOTS) combined with microarray analysis. We identified 767 genes with significantly different levels of expression by intracellular bacteria, of which 330 showed increased expression and 437 showed decreased expression. Affected genes represented all aspects of cellular life including information storage and processing, cellular processes and signaling, and metabolism. In general, intracellular gene expression demonstrated a pattern of environmental sensing, bacterial response, and metabolic adaptation to the phagosomal environment. Deletion of various SCOTS-identified genes affected bacterial entry into macrophages and intracellular replication. We also show that intracellular B. cenocepacia is cytotoxic towards the macrophage host, and capable of spread to neighboring cells, a role dependent on SCOTS-identified genes. In particular, genes involved in bacterial motility, cobalamin biosynthesis, the type VI secretion system, and membrane modification contributed greatly to macrophage entry and subsequent intracellular behavior of B. cenocepacia. Conclusions: B. cenocepacia enters macrophages, adapts to the phagosomal environment, replicates within a modified phagosome, and exhibits cytotoxicity towards the host cells. The analysis of the transcriptomic response of intracellular B. cenocepacia reveals that metabolic adaptation to a new niche plays a major role in the survival of B. cenocepacia in macrophages. This adaptive response does not require the expression of any specific virulence-associated factor, which is consistent with the opportunistic nature of this microorganism. Further investigation into the remaining SCOTS-identified genes will provide a more complete picture of the adaptive response of B. cenocepacia to the host cell environment.

Transcriptomic and proteomic analyses of the Aspergillus fumigatus hypoxia response using an oxygen-controlled fermenter

Bridget Barker — 2012-02-06T00:00:00Z

Background: Aspergillus fumigatus is a mold responsible for the majority of cases of aspergillosis in humans. To survive in the human body, A. fumigatus must adapt to microenvironments that are often characterized by low nutrient and oxygen availability. Recent research suggests that the ability of A. fumigatus and other pathogenic fungi to adapt to hypoxia contributes to their virulence. However, molecular mechanisms of A. fumigatus hypoxia adaptation are poorly understood. Thus, to better understand how A. fumigatus adapts to hypoxic microenvironments found in vivo during human fungal pathogenesis, the dynamic changes of the fungal transcriptome and proteome in hypoxia were investigated over a period of 24 hours utilizing an oxygen-controlled fermenter system. Results: Significant increases in transcripts associated with iron and sterol metabolism, the cell wall, the GABA shunt, and transcriptional regulators were observed in response to hypoxia. A concomitant reduction in transcripts was observed with ribosome and terpenoid backbone biosynthesis, TCA cycle, amino acid metabolism and RNA degradation. Analysis of changes in transcription factor mRNA abundance shows that hypoxia induces significant positive and negative changes that may be important for regulating the hypoxia response in this pathogenic mold. Growth in hypoxia resulted in changes in the protein levels of several glycolytic enzymes, but these changes were not always reflected by the corresponding transcriptional profiling data. However, a good correlation overall (R2 = 0.2, p < 0.05) existed between the transcriptomic and proteomics datasets for all time points. The lack of correlation between some transcript levels and their subsequent protein levels suggests another regulatory layer of the hypoxia response in A. fumigatus. Conclusions: Taken together, our data suggest a robust cellular response that is likely regulated both at the transcriptional and post-transcriptional level in response to hypoxia by the human pathogenic mold A. fumigatus. As with other pathogenic fungi, the induction of glycolysis and transcriptional down-regulation of the TCA cycle and oxidative phosphorylation appear to major components of the hypoxia response in this pathogenic mold. In addition, a significant induction of the transcripts involved in ergosterol biosynthesis is consistent with previous observations in the pathogenic yeasts Candida albicans and Cryptococcus neoformans indicating conservation of this response to hypoxia in pathogenic fungi. Because ergosterol biosynthesis enzymes also require iron as a co-factor, the increase in iron uptake transcripts is consistent with an increased need for iron under hypoxia. However, unlike C. albicans and C. neoformans, the GABA shunt appears to play an important role in reducing NADH levels in response to hypoxia in A. fumigatus and it will be intriguing to determine whether this is critical for fungal virulence. Overall, regulatory mechanisms of the A. fumigatus hypoxia response appear to involve both transcriptional and post-transcriptional control of transcript and protein levels and thus provide candidate genes for future analysis of their role in hypoxia adaptation and fungal virulence.

Ac/Ds-transposon activation tagging in poplar: a powerful tool for gene discovery

Matthias Fladung — 2012-02-06T00:00:00Z

Background: Rapid improvements in the development of new sequencing technologies have led to the availability of genome sequences of more than 300 organisms today. Thanks to bioinformatic analyses, prediction of gene models and protein-coding transcripts has become feasible. Various reverse and forward genetics strategies have been followed to determine the functions of these gene models and regulatory sequences. Using T-DNA or transposons as tags, significant progress has been made by using "Knock-in" approaches ("gain-of-function" or "activation tagging") in different plant species but not in perennial plants species, e.g. long-lived trees. Here, large scale gene tagging resources are still lacking. Results: We describe the first application of an inducible transposon-based activation tagging system for a perennial plant species, as example a poplar hybrid (P. tremula L. x P. tremuloides Michx). Four activation-tagged populations comprising a total of 12,083 individuals derived from 23 independent "Activation Tagging Ds" (ATDs) transgenic lines were produced and phenotyped. To date, 29 putative variants have been isolated and new ATDs genomic positions were successfully determined for 24 of those. Sequences obtained were blasted against the publicly available genome sequence of P. trichocarpa v2.0 (Phytozome v7.0; http://www.phytozome.net/poplar) revealing possible transcripts for 17 variants.In a second approach, 300 randomly selected individuals without any obvious phenotypic alterations were screened for ATDs excision. For one third of those transposition of ATDs was confirmed and in about 5% of these cases genes were tagged. Conclusions: The novel strategy of first genotyping and then phenotyping a tagging population as proposed here is, in particular, applicable for long-lived, difficult to transform plant species. We could demonstrate the power of the ATDs transposon approach and the simplicity to induce ATDs transposition in vitro. Since a transposon is able to pass chromosomal boundaries, only very few primary transposon-carrying transgenic lines are required for the establishment of large transposon tagging populations. In contrast to T-DNA-based activation tagging, which is plagued by a lack of transformation efficiency and its time consuming nature, this for the first time, makes it feasible one day to tag (similarly to Arabidopsis) every gene within a perennial plant genome.

Medicago truncatula transporter database: a comprehensive database resource for M. truncatula transporters

Zhenyan Miao — 2012-02-06T00:00:00Z

Background: Medicago truncatula has been chosen as a model species for genomic studies. It is closely related to an important legume, alfalfa. Transporters are a large group of membrane-spanning proteins. They deliver essential nutrients, eject waste products, and assist the cell in sensing environmental conditions by forming a complex system of pumps and channels. Although studies have effectively characterized individual M. truncatula transporters in several databases, until now there has been no available systematic database that includes all transporters in M. truncatula.DescriptionThe M. truncatula transporter database (MTDB) contains comprehensive information on the transporters in M. truncatula. Based on the TransportTP method, we have presented a novel prediction pipeline. A total of 3,665 putative transporters have been annotated based on International Medicago Genome Annotated Group (IMGAG) V3.5 V3 and the M. truncatula Gene Index (MTGI) V10.0 releases and assigned to 162 families according to the transporter classification system. These families were further classified into seven types according to their transport mode and energy coupling mechanism. Extensive annotations referring to each protein were generated, including basic protein function, expressed sequence tag (EST) mapping, genome locus, three-dimensional template prediction, transmembrane segment, and domain annotation. A chromosome distribution map and text-based Basic Local Alignment Search Tools were also created. In addition, we have provided a way to explore the expression of putative M. truncatula transporter genes under stress treatments. Conclusions: In summary, the MTDB enables the exploration and comparative analysis of putative transporters in M. truncatula. A user-friendly web interface and regular updates make MTDB valuable to researchers in related fields. The MTDB is freely available now to all users at http://bioinformatics.cau.edu.cn/MtTransporter/.

Heritable genome-wide variation of gene expression and promoter methylation between wild and domesticated chickens

Daniel Natt — 2012-02-04T00:00:00Z

Background: Variations in gene expression, mediated by epigenetic mechanisms, may cause broad phenotypic effects in animals. However, it has been debated to what extent expression variation and epigenetic modifications, such as patterns of DNA methylation, are transferred across generations, and therefore it is uncertain what role epigenetic variation may play in adaptation. Results: In Red Junglefowl, ancestor of domestic chickens, gene expression and methylation profiles in thalamus/hypothalamus differed substantially from that of a domesticated egg laying breed. Expression as well as methylation differences were largely maintained in the offspring, demonstrating reliable inheritance of epigenetic variation. Some of the inherited methylation differences were tissue-specific, and the differential methylation at specific loci were little changed after eight generations of intercrossing between Red Junglefowl and domesticated laying hens. There was an over-representation of differentially expressed and methylated genes in selective sweep regions associated with chicken domestication. Conclusions: Our results show that epigenetic variation is inherited in chickens, and we suggest that selection of favourable epigenomes, either by selection of genotypes affecting epigenetic states, or by selection of methylation states which are inherited independently of sequence differences, may have been an important aspect of chicken domestication.

Analysis of the SOS response of Vibrio and other bacteria with multiple chromosomes

Neus Sanchez-Alberola — 2012-02-03T00:00:00Z

Background: The SOS response is a well-known regulatory network present in most bacteria and aimed at addressing DNA damage. It has also been linked extensively to stress-induced mutagenesis, virulence and the emergence and dissemination of antibiotic resistance determinants. Recently, the SOS response has been shown to regulate the activity of integrases in the chromosomal superintegrons of the Vibrionaceae, which encompasses a wide range of pathogenic species harboring multiple chromosomes. Here we combine in silico and in vitro techniques to perform a comparative genomics analysis of the SOS regulon in the Vibrionaceae, and we extend the methodology to map this transcriptional network in other bacterial species harboring multiple chromosomes. Results: Our analysis provides the first comprehensive description of the SOS response in a family (Vibrionaceae) that includes major human pathogens. It also identifies several previously unreported members of the SOS transcriptional network, including two proteins of unknown function. The analysis of the SOS response in other bacterial species with multiple chromosomes uncovers additional regulon members and reveals that there is a conserved core of SOS genes, and that specialized additions to this basic network take place in different phylogenetic groups. Our results also indicate that across all groups the main elements of the SOS response are always found in the large chromosome, whereas specialized additions are found in the smaller chromosomes and plasmids. Conclusions: Our findings confirm that the SOS response of the Vibrionaceae is strongly linked with pathogenicity and suggest that the characterization of the newly identified members of this regulon could provide key insights into the pathogenesis of Vibrio. The persistent location of key SOS genes in the large chromosome across several bacterial groups confirms that the SOS response plays an essential role in these organisms and sheds light into the mechanisms of evolution of global transcriptional networks involved in adaptability and rapid response to environmental changes, suggesting that small chromosomes may act as evolutionary test beds for the rewiring of transcriptional networks.

Post-genomic analyses of fungal lignocellulosic biomass degradation reveal the unexpected potential of the plant pathogen Ustilago maydis

Marie Couturier — 2012-02-02T00:00:00Z

Background: Filamentous fungi are potent biomass degraders due to their ability to thrive in ligno(hemi)cellulose-rich environments. During the last decade, fungal genome sequencing initiatives have yielded abundant information on the genes that are putatively involved in lignocellulose degradation. At present, additional experimental studies are essential to provide insights into the fungal secreted enzymatic pools involved in lignocellulose degradation. Results: In this study, we performed a wide analysis of 20 filamentous fungi for which genomic data are available to investigate their biomass-hydrolysis potential. A comparison of fungal genomes and secretomes using enzyme activity profiling revealed discrepancies in carbohydrate active enzymes (CAZymes) sets dedicated to plant cell wall. Investigation of the contribution made by each secretome to the saccharification of wheat straw demonstrated that most of them individually supplemented the industrial Trichoderma reesei CL847 enzymatic cocktail. Unexpectedly, the most striking effect was obtained with the phytopathogen Ustilago maydis that improved the release of total sugars by 57% and of glucose by 22%. Proteomic analyses of the best-performing secretomes indicated a specific enzymatic mechanism of U. maydis that is likely to involve oxido-reductases and hemicellulases. Conclusion: This study provides insight into the lignocellulose-degradation mechanisms by filamentous fungi and allows for the identification of a number of enzymes that are potentially useful to further improve the industrial lignocellulose bioconversion process.