Abstract

Background

Methylotrophic yeast species (e.g. Hansenula polymorpha, Pichia pastoris) can grow on methanol as sole source of carbon and energy. These organisms are important cell factories for the production of recombinant proteins, but are also used in fundamental research as model organisms to study peroxisome biology. During exponential growth on glucose, cells of H. polymorpha typically contain a single, small peroxisome that is redundant for growth while on methanol multiple, enlarged peroxisomes are present. These organelles are crucial to support growth on methanol, as they contain key enzymes of methanol metabolism.

In this study, changes in the transcriptional profiles during adaptation of H. polymorpha cells from glucose- to methanol-containing media were investigated using DNA-microarray analyses.

Results

Two hours after the shift of cells from glucose to methanol nearly 20% (1184 genes) of the approximately 6000 annotated H. polymorpha genes were significantly upregulated with at least a two-fold differential expression. Highest upregulation (> 300-fold) was observed for the genes encoding the transcription factor Mpp1 and formate dehydrogenase, an enzyme of the methanol dissimilation pathway. Upregulated genes also included genes encoding other enzymes of methanol metabolism as well as of peroxisomal β-oxidation.

A moderate increase in transcriptional levels (up to 4-fold) was observed for several PEX genes, which are involved in peroxisome biogenesis. Only PEX11 and PEX32 were higher upregulated. In addition, an increase was observed in expression of the several ATG genes, which encode proteins involved in autophagy and autophagy processes. The strongest upregulation was observed for ATG8 and ATG11.

Approximately 20% (1246 genes) of the genes were downregulated. These included glycolytic genes as well as genes involved in transcription and translation.

Conclusion

Transcriptional profiling of H. polymorpha cells shifted from glucose to methanol showed the expected downregulation of glycolytic genes together with upregulation of the methanol utilisation pathway. This serves as a confirmation and validation of the array data obtained. Consistent with this, also various PEX genes were upregulated. The strong upregulation of ATG genes is possibly due to induction of autophagy processes related to remodeling of the cell architecture required to support growth on methanol. These processes may also be responsible for the enhanced peroxisomal β-oxidation, as autophagy leads to recycling of membrane lipids. The prominent downregulation of transcription and translation may be explained by the reduced growth rate on methanol (t_d glucose 1 h vs t_d methanol 4.5 h).

Abstract

Background

In most species of mammals, the TRB locus has the common feature of a library of TRBV genes positioned at the 5'- end of two in tandem aligned D-J-C gene clusters, each composed of a single TRBD gene, 6-7 TRBJ genes and one TRBC gene. An enhancer located at the 3'end of the last TRBC and a well-defined promoter situated at the 5'end of the TRBD gene and/or a undefined promoter situated at the 5'end of the TRBD2 are sufficient to generate the full recombinase accessibility at the locus. In ruminant species, the 3'end of the TRB locus is characterized by the presence of three D-J-C clusters, each constituted by a single TRBD, 5-7 TRBJ and one TRBC genes with the center cluster showing a structure combined with the clusters upstream and downstream, suggesting that a unequal crossover occurred in the duplication. An enhancer downstream the last TRBC, and a promoter at the 5'-end of each TRBD gene are also present.

Results

In this paper we focused our attention on the analysis of a large number of sheep TR β-chain transcripts derived from four different lymphoid tissues of three diverse sheep breed animals to certify the use and frequency of the three gene clusters in the β-chain repertoire. As the sheep TRB locus genomic organization is known, the exact interpretation of the V-D-J rearrangements was fully determined. Our results clearly demonstrate that sheep β-chain constitutes a level of variability that is substantially larger than that described in other mammalian species. This is due not only to the increase of the number of D and J genes available to the somatic recombination, but also to the presence of the trans-rearrangement process. Moreover, the functional complexity of β-chain repertoire is resolved by other mechanisms such as alternative cis- and trans-splicing and recombinational diversification that seems to affect the variety of the constant region.

Conclusion

All together our data demonstrate that a disparate set of molecular mechanisms operate to perform a diversified repertoire in the sheep β-chain and this could confer some special biological properties to the corresponding αβ T cells in the ruminant lineage.

Abstract

Background

Studies using microarray analysis of colorectal cancer have been generally beleaguered by the lack of a normal cell population of the same lineage as the tumor cell. One of the main objectives of this study was to generate a reference gene expression data set for normal colonic epithelium which can be used in comparisons with diseased tissues, as well as to provide a dataset that could be used as a baseline for studies in alternative splicing.

Results

We present a dependable expression reference data set for non-neoplastic colonic epithelial cells. An enriched population of fresh colon epithelial cells were obtained from non-neoplastic, colectomy specimens and analyzed using Affymetrix GeneChip EXON 1.0 ST arrays. For demonstration purposes, we have compared the data derived from these cells to a publically available set of tumor and matched normal colon data. This analysis allowed an assessment of global gene expression alterations and demonstrated that adjacent normal tissues, with a high degree of cellular heterogeneity, are not always representative of normal cells for comparison to tumors which arise from the colon epithelium. We also examined alternative splicing events in tumors compared to normal colon epithelial cells.

Conclusions

The findings from this study represent the first comprehensive expression profile for non-neoplastic colonic epithelial cells reported. Our analysis of splice variants illustrate that this is a very labor intensive procedure, requiring vigilant examination of the data. It is projected that the contribution of this set of data derived from pure colonic epithelial cells will enhance studies in colon-related disease and offer a vital baseline for studies aimed at elucidating the mechanisms of alternative splicing.

Abstract

Background

RNAi screens via pooled short hairpin RNAs (shRNAs) have recently become a powerful tool for the identification of essential genes in mammalian cells. In the past years, several pooled large-scale shRNA screens have identified a variety of genes involved in cancer cell proliferation. All of those studies employed microarray analysis, utilizing either the shRNA's half hairpin sequence or an additional shRNA-associated 60 nt barcode sequence as a molecular tag. Here we describe a novel method to decode pooled RNAi screens, namely barcode tiling array analysis, and demonstrate how this approach can be used to precisely quantify the abundance of individual shRNAs from a pool.

Results

We synthesized DNA microarrays with six overlapping 25 nt long tiling probes complementary to each unique 60 nt molecular barcode sequence associated with every shRNA expression construct. By analyzing dilution series of expression constructs we show how our approach allows quantification of shRNA abundance from a pool and how it clearly outperforms the commonly used analysis via the shRNA's half hairpin sequences. We further demonstrate how barcode tiling arrays can be used to predict anti-proliferative effects of individual shRNAs from pooled negative selection screens. Out of a pool of 305 shRNAs, we identified 28 candidate shRNAs to fully or partially impair the viability of the breast carcinoma cell line MDA-MB-231. Individual validation of a subset of eleven shRNA expression constructs with potential inhibitory, as well as non-inhibitory, effects on the cell line proliferation provides further evidence for the accuracy of the barcode tiling approach.

Conclusions

In summary, we present an improved method for the rapid, quantitative and statistically robust analysis of pooled RNAi screens. Our experimental approach, coupled with commercially available lentiviral vector shRNA libraries, has the potential to greatly facilitate the discovery of putative targets for cancer therapy as well as sensitizers of drug toxicity.

Abstract

Background

Non-genotoxic carcinogens are notoriously difficult to identify as they do not damage DNA directly and have diverse modes of action, necessitating long term in vivo studies. The early effects of the classic rodent non-genotoxic hepatocarcinogen phenobarbital have been investigated in the Fisher rat using a combination of metabolomics and transcriptomics, to investige early stage mechanistic changes that are predictive of longer term pathology.

Results

Liver and blood plasma were profiled across 14 days, and multivariate statistics used to identify perturbed pathways. Both metabolomics and transcriptomics detected changes in the liver which were dose dependent, even after one day of exposure. Integration of the two datasets associated perturbations with specific pathways. Hepatic glycogen was decreased due to a decrease in synthesis, and plasma triglycerides were decreased due to an increase in fatty acid uptake by the liver. Hepatic succinate was increased and this was associated with increased heme biosynthesis. Glutathione synthesis was also increased, presumably in response to oxidative stress. Liquid Chromatography Mass Spectrometry demonstrated a remodeling of lipid species, possibly resulting from proliferation of the smooth endoplasmic reticulum.

Conclusions

The data fusion of metabolomic and transcriptomic changes proved to be a highly sensitive approach for monitoring early stage changes in altered hepatic metabolism, oxidative stress and cytochrome P450 induction simultaneously. This approach is particularly useful in interpreting changes in metabolites such as succinate which are hubs of metabolism.

Abstract

Background

Biological nitrogen fixation is highly controlled at the transcriptional level by regulatory networks that respond to the availability of fixed nitrogen. In many diazotrophs, addition of excess ammonium in the growth medium results in immediate repression of nif gene transcription. Although the regulatory cascades that control the transcription of the nif genes in proteobacteria have been well investigated, there are limited data on the kinetics of ammonium-dependent repression of nitrogen fixation.

Results

Here we report a global transcriptional profiling analysis of nitrogen fixation and ammonium repression in Pseudomonas stutzeri A1501, a root-associated and nitrogen-fixing bacterium. A total of 166 genes, including those coding for the global nitrogen regulation (Ntr) and Nif-specific regulatory proteins, were upregulated under nitrogen fixation conditions but rapidly downregulated as early as 10 min after ammonium shock. Among these nitrogen fixation-inducible genes, 95 have orthologs in each of Azoarcus sp. BH72 and Azotobacter vinelandii AvoP. In particular, a 49-kb expression island containing nif and other associated genes was markedly downregulated by ammonium shock. Further functional characterization of pnfA, a new NifA-σ⁵⁴-dependent gene chromosomally linked to nifHDK, is reported. This gene encodes a protein product with an amino acid sequence similar to that of five hypothetical proteins found only in diazotrophic strains. No noticeable differences in the transcription of nifHDK were detected between the wild type strain and pnfA mutant. However, the mutant strain exhibited a significant decrease in nitrogenase activity under microaerobic conditions and lost its ability to use nitrate as a terminal electron acceptor for the support of nitrogen fixation under anaerobic conditions.

Conclusions

Based on our results, we conclude that transcriptional regulation of nif gene expression in A1501 is mediated by the nif-specific and ntr gene regulatory systems. Furthermore, microarray and mutational analyses revealed that many genes of unknown function may play some essential roles in controlling the expression or activity of nitrogenase. The findings presented here establish the foundation for further studies on the physiological function of nitrogen fixation-inducible genes.

Abstract

Background

Considerable progress has been made in our understanding of sex determination and dosage compensation mechanisms in model organisms such as C. elegans, Drosophila and M. musculus. Strikingly, the mechanism involved in sex determination and dosage compensation are very different among these three model organisms. Birds present yet another situation where the heterogametic sex is the female. Sex determination is still poorly understood in birds and few key determinants have so far been identified. In contrast to most other species, dosage compensation of bird sex chromosomal genes appears rather ineffective.

Results

By comparing microarrays from microdissected primitive streak from single chicken embryos, we identified a large number of genes differentially expressed between male and female embryos at a very early stage (Hamburger and Hamilton stage 4), long before any sexual differentiation occurs. Most of these genes are located on the Z chromosome, which indicates that dosage compensation is ineffective in early chicken embryos. Gene ontology analyses, using an enhanced annotation tool for Affymetrix probesets of the chicken genome developed in our laboratory (called Manteia), show that among these male-biased genes found on the Z chromosome, more than 20 genes play a role in sex differentiation.

Conclusions

These results corroborate previous studies demonstrating the rather inefficient dosage compensation for Z chromosome in birds and show that this sexual dimorphism in gene regulation is observed long before the onset of sexual differentiation. These data also suggest a potential role of non-compensated Z-linked genes in somatic sex differentiation in birds.

Abstract

Background

Microarray technologies produced large amount of data. In a previous study, we have shown the interest of k-Nearest Neighbour approach for restoring the missing gene expression values, and its positive impact of the gene clustering by hierarchical algorithm. Since, numerous replacement methods have been proposed to impute missing values (MVs) for microarray data. In this study, we have evaluated twelve different usable methods, and their influence on the quality of gene clustering. Interestingly we have used several datasets, both kinetic and non kinetic experiments from yeast and human.

Results

We underline the excellent efficiency of approaches proposed and implemented by Bo and co-workers and especially one based on expected maximization (EM_array). These improvements have been observed also on the imputation of extreme values, the most difficult predictable values. We showed that the imputed MVs have still important effects on the stability of the gene clusters. The improvement on the clustering obtained by hierarchical clustering remains limited and, not sufficient to restore completely the correct gene associations. However, a common tendency can be found between the quality of the imputation method and the gene cluster stability. Even if the comparison between clustering algorithms is a complex task, we observed that k-means approach is more efficient to conserve gene associations.

Conclusions

More than 6.000.000 independent simulations have assessed the quality of 12 imputation methods on five very different biological datasets. Important improvements have so been done since our last study. The EM_array approach constitutes one efficient method for restoring the missing expression gene values, with a lower estimation error level. Nonetheless, the presence of MVs even at a low rate is a major factor of gene cluster instability. Our study highlights the need for a systematic assessment of imputation methods and so of dedicated benchmarks. A noticeable point is the specific influence of some biological dataset.

Abstract

Background

The SLC17 family of transporters transports the amino acids: glutamate and aspartate, and, as shown recently, also nucleotides. Vesicular glutamate transporters are found in distinct species, such as C. elegans, but the evolutionary origin of most of the genes in this family has been obscure.

Results

Our phylogenetic analysis shows that the SLC17 family consists of four main phylogenetic clades which were all present before the divergence of the insect lineage. One of these clades has not been previously described and it is not found in vertebrates. The clade containing Slc17a9 had the most restricted evolutionary history with only one member in most species. We detected expression of Slc17a1-17a4 only in the peripheral tissues but not in the CNS, while Slc17a5- Slc17a9 are highly expressed in both the CNS and periphery.

Conclusions

The in situ hybridization studies on vesicular nucleotide transporter revealed high expression throughout the cerebral cortex, certain areas in the hippocampus and in specific nuclei of the hypothalamus and thalamus. Some of the regions with high expression, such as the medial habenula and the dentate gyrus of the hippocampus, are important sites for purinergic neurotransmission. Noteworthy, other areas relying on purine-mediated signaling, such as the molecular layer of the dentate gyrus and the periaqueductal gray, lack or have a very low expression of Slc17a9, suggesting that there could be another nucleotide transporter in these regions.

Abstract

Background

Transmembrane receptor kinases play critical roles in both animal and plant signaling pathways regulating growth, development, differentiation, cell death, and pathogenic defense responses. In Arabidopsis thaliana, there are at least 223 Leucine-rich repeat receptor-like kinases (LRR-RLKs), representing one of the largest protein families. Although functional roles for a handful of LRR-RLKs have been revealed, the functions of the majority of members in this protein family have not been elucidated.

Results

As a resource for the in-depth analysis of this important protein family, the complementary DNA sequences (cDNAs) of 194 LRR-RLKs were cloned into the Gateway^R donor vector pDONR/Zeo^R and analyzed by DNA sequencing. Among them, 157 clones showed sequences identical to the predictions in the Arabidopsis sequence resource, TAIR8. The other 37 cDNAs showed gene structures distinct from the predictions of TAIR8, which was mainly caused by alternative splicing of pre-mRNA. Most of the genes have been further cloned into Gateway^R destination vectors with GFP or FLAG epitope tags and have been transformed into Arabidopsis for in planta functional analysis. All clones from this study have been submitted to the Arabidopsis Biological Resource Center (ABRC) at Ohio State University for full accessibility by the Arabidopsis research community.

Conclusions

Most of the Arabidopsis LRR-RLK genes have been isolated and the sequence analysis showed a number of alternatively spliced variants. The generated resources, including cDNA entry clones, expression constructs and transgenic plants, will facilitate further functional analysis of the members of this important gene family.

Abstract

Background

The approaches for shotgun-based sequencing of vertebrate genomes are now well-established, and have resulted in the generation of numerous draft whole-genome sequence assemblies. In contrast, the process of refining those assemblies to improve contiguity and increase accuracy (known as 'sequence finishing') remains tedious, labor-intensive, and expensive. As a result, the vast majority of vertebrate genome sequences generated to date remain at a draft stage.

Results

To date, our genome sequencing efforts have focused on comparative studies of targeted genomic regions, requiring sequence finishing of large blocks of orthologous sequence (average size 0.5-2 Mb) from various subsets of 75 vertebrates. This experience has provided a unique opportunity to compare the relative effort required to finish shotgun-generated genome sequence assemblies from different species, which we report here. Importantly, we found that the sequence assemblies generated for the same orthologous regions from various vertebrates show substantial variation with respect to misassemblies and, in particular, the frequency and characteristics of sequence gaps. As a consequence, the work required to finish different species' sequences varied greatly. Application of the same standardized methods for finishing provided a novel opportunity to "assay" characteristics of genome sequences among many vertebrate species. It is important to note that many of the problems we have encountered during sequence finishing reflect unique architectural features of a particular vertebrate's genome, which in some cases may have important functional and/or evolutionary implications. Finally, based on our analyses, we have been able to improve our procedures to overcome some of these problems and to increase the overall efficiency of the sequence-finishing process, although significant challenges still remain.

Conclusion

Our findings have important implications for the eventual finishing of the draft whole-genome sequences that have now been generated for a large number of vertebrates.

Abstract

Background

Early embryonic loss is a large contributor to infertility in cattle. Although genetic factors are known to affect early embryonic development, the discovery of such factors has been a serious challenge. The objective of this study was to identify genes differentially expressed between blastocysts and degenerative embryos at early stages of development.

Results

Using microarrays, genome-wide RNA expression was profiled and compared for in vitro fertilization (IVF) - derived blastocysts and embryos undergoing degenerative development up to the same time point. Surprisingly similar transcriptomic profiles were found in degenerative embryos and blastocysts. Nonetheless, we identified 67 transcripts that significantly differed between these two groups of embryos at a 15% false discovery rate, including 33 transcripts showing at least a two-fold difference. Several signaling and metabolic pathways were found to be associated with the developmental status of embryos, among which were previously known important steroid biosynthesis and cell communication pathways in early embryonic development.

Conclusions

This study presents the first direct and comprehensive comparison of transcriptomes between IVF blastocysts and degenerative embryos, providing important information for potential genes and pathways associated with early embryonic development.

Abstract

Background

Neurofibrillary tangles (NFT), a cardinal neuropathological feature of Alzheimer's disease (AD) that is highly correlated with synaptic loss and dementia severity, appear to be partly attributable to increased phosphorylation of the microtubule stabilizing protein tau at certain AD-related residues. Identifying the kinases involved in the pathologic phosphorylation of tau may provide targets at which to aim new AD-modifying treatments.

Results

We report results from a screen of 572 kinases in the human genome for effects on tau hyperphosphorylation using a loss of function, high-throughput RNAi approach. We confirm effects of three kinases from this screen, the eukaryotic translation initiation factor 2 α kinase 2 (EIF2AK2), the dual-specificity tyrosine-(Y)-phosphorylation regulated kinase 1A (DYRK1A), and the A-kinase anchor protein 13 (AKAP13) on tau phosphorylation at the 12E8 epitope (serine 262/serine 356). We provide evidence that EIF2AK2 effects may result from effects on tau protein expression, whereas DYRK1A and AKAP13 are likely more specifically involved in tau phosphorylation pathways.

Conclusions

These findings identify novel kinases that phosphorylate tau protein and provide a valuable reference data set describing the kinases involved in phosphorylating tau at an AD-relevant epitope.

Abstract

Background

A large amount of experimental data generated by modern high-throughput technologies is available through various public repositories. Our knowledge about molecular interaction networks, functional biological pathways and transcriptional regulatory modules is rapidly expanding, and is being organized in lists of functionally related genes. Jointly, these two sources of information hold a tremendous potential for gaining new insights into functioning of living systems.

Results

Genomics Portals platform integrates access to an extensive knowledge base and a large database of human, mouse, and rat genomics data with basic analytical visualization tools. It provides the context for analyzing and interpreting new experimental data and the tool for effective mining of a large number of publicly available genomics datasets stored in the back-end databases. The uniqueness of this platform lies in the volume and the diversity of genomics data that can be accessed and analyzed (gene expression, ChIP-chip, ChIP-seq, epigenomics, computationally predicted binding sites, etc), and the integration with an extensive knowledge base that can be used in such analysis.

Conclusion

The integrated access to primary genomics data, functional knowledge and analytical tools makes Genomics Portals platform a unique tool for interpreting results of new genomics experiments and for mining the vast amount of data stored in the Genomics Portals backend databases. Genomics Portals can be accessed and used freely at http://GenomicsPortals.org.

Abstract

Background

The chemical components of sex pheromones have been determined for more than a thousand moth species, but so far only a handful of genes encoding enzymes responsible for the biosynthesis of these compounds have been identified. For understanding the evolution of moth sexual communication, it is essential to know which genes are involved in the production of specific pheromone components and what controls the variation in their relative frequencies in the pheromone blend. We used a transcriptomic approach to characterize the pheromone gland of the Noctuid moth Heliothis virescens, an important agricultural pest, in order to obtain substantial general sequence information and to identify a range of candidate genes involved in the pheromone biosynthetic pathway.

Results

To facilitate identifying sets of genes involved in a broad range of processes and to capture rare transcripts, we developed our majority of ESTs from a normalized cDNA library of Heliothis virescens pheromone glands (PG). Combining these with a non-normalized library yielded a total of 17,233 ESTs, which assembled into 2,082 contigs and 6,228 singletons. Using BLAST searches of the NR and Swissprot databases we were able to identify a large number of putative unique gene elements (unigenes), which we compared to those derived from previous transcriptomic surveys of the larval stage of Heliothis virescens. The distribution of unigenes among GO Biological Process functional groups shows an overall similarity between PG and larval transcriptomes, but with distinct enrichment of specific pathways in the PG. In addition, we identified a large number of candidate genes in the pheromone biosynthetic pathways.

Conclusion

These data constitute one of the first large-scale EST-projects for Noctuidae, a much-needed resource for exploring these pest species. Our analysis shows a surprisingly complex transcriptome and we identified a large number of potential pheromone biosynthetic pathway and immune-related genes that can be applied to population and systematic studies of Heliothis virescens and other Noctuidae.

Abstract

Background

The extracellular promastigote and the intracellular amastigote stages alternate in the digenetic life cycle of the trypanosomatid parasite Leishmania. Amastigotes develop inside parasitophorous vacuoles of mammalian phagocytes, where they tolerate extreme environmental conditions. Temperature increase and pH decrease are crucial factors in the multifactorial differentiation process of promastigotes to amastigotes. Although expression profiling approaches for axenic, cell culture- and lesion-derived amastigotes have already been reported, the specific influence of temperature increase and acidification of the environment on developmental regulation of genes has not been previously studied. For the first time, we have used custom L. infantum genomic DNA microarrays to compare the isolated and the combined effects of both factors on the transcriptome.

Results

Immunofluorescence analysis of promastigote-specific glycoprotein gp46 and expression modulation analysis of the amastigote-specific A2 gene have revealed that concomitant exposure to temperature increase and acidification leads to amastigote-like forms. The temperature-induced gene expression profile in the absence of pH variation resembles the profile obtained under combined exposure to both factors unlike that obtained for exposure to acidification alone. In fact, the subsequent fold change-based global iterative hierarchical clustering analysis supports these findings.

Conclusions

The specific influence of temperature and pH on the differential regulation of genes described in this study and the evidence provided by clustering analysis is consistent with the predominant role of temperature increase over extracellular pH decrease in the amastigote differentiation process, which provides new insights into Leishmania physiology.

Abstract

Background

Recently, a number of high-resolution genome-wide maps of nucleosome locations in S. cerevisiae have been derived experimentally. However, nucleosome positions are determined in vivo by the combined effects of numerous factors. Consequently, nucleosomes are not simple static units, which may explain the discrepancies in reported nucleosome positions as measured by different experiments. In order to more accurately depict the genome-wide nucleosome distribution, we integrated multiple nucleosomal positioning datasets using a multi-angle analysis strategy.

Results

To evaluate the contribution of chromatin structure to transcription, we used the vast amount of available nucleosome analyzed data. Analysis of this data allowed for the comprehensive identification of the connections between promoter nucleosome positioning patterns and various transcription-dependent properties. Further, we characterised the function of nucleosome destabilisation in the context of transcription regulation. Our results indicate that genes with similar nucleosome occupancy patterns share general transcription attributes. We identified the local regulatory correlation (LRC) regions for two distinct types of nucleosomes and we assessed their regulatory properties. We also estimated the nucleosome reproducibility and measurement accuracy for high-confidence transcripts. We found that by maintaining a distance of ~13 bp between the upstream border of the +1 nucleosome and the transcription start sites (TSSs), the stable +1 nucleosome may form a barrier against the accessibility of the TSS and shape an optimum chromatin conformation for gene regulation. An in-depth analysis of nucleosome positioning in normally growing and heat shock cells suggested that the extent and patterns of nucleosome sliding are associated with gene activation.

Conclusions

Our results, which combine different types of data, suggest that cross-platform information, including discrepancy and consistency, reflects the mechanisms of nucleosome packaging in vivo more faithfully than individual studies. Furthermore, nucleosomes can be divided into two classes according to their stable and dynamic characteristics. We found that two different nucleosome-positioning characteristics may significantly impact transcription programs. Besides, some positioned-nucleosomes are involved in the transition from stable state to dynamic state in response to abrupt environmental changes.

Abstract

Background

Plasmodium falciparum is the main causative agent of malaria. Of the 5 484 predicted genes of P. falciparum, about 57% do not have sufficient sequence similarity to characterized genes in other species to warrant functional assignments. Non-homology methods are thus needed to obtain functional clues for these uncharacterized genes. Gene expression data have been widely used in the recent years to help functional annotation in an intra-species way via the so-called Guilt By Association (GBA) principle.

Results

We propose a new method that uses gene expression data to assess inter-species annotation transfers. Our approach starts from a set of likely orthologs between a reference species (here S. cerevisiae and D. melanogaster) and a query species (P. falciparum). It aims at identifying clusters of coexpressed genes in the query species whose coexpression has been conserved in the reference species. These conserved clusters of coexpressed genes are then used to assess annotation transfers between genes with low sequence similarity, enabling reliable transfers of annotations from the reference to the query species. The approach was used with transcriptomic data sets of P. falciparum, S. cerevisiae and D. melanogaster, and enabled us to propose with high confidence new/refined annotations for several dozens hypothetical/putative P. falciparum genes. Notably, we revised the annotation of genes involved in ribosomal proteins and ribosome biogenesis and assembly, thus highlighting several potential drug targets.

Conclusions

Our approach uses both sequence similarity and gene expression data to help inter-species gene annotation transfers. Experiments show that this strategy improves the accuracy achieved when using solely sequence similarity and outperforms the accuracy of the GBA approach. In addition, our experiments with P. falciparum show that it can infer a function for numerous hypothetical genes.

Abstract

Background

Synaptotagmin proteins were first identified in nervous tissue, residing in synaptic vesicles. Synaptotagmins were subsequently found to form a large family, some members of which play important roles in calcium triggered exocytic events. These members have been investigated intensively, but other family members are not well understood, making it difficult to grasp the meaning of family membership in functional terms. Further difficulty arises as families are defined quite legitimately in different ways: by common descent or by common possession of distinguishing features. One definition does not necessarily imply the other. The evolutionary range of genome sequences now available, can shed more light on synaptotagmin gene phylogeny and clarify family relationships. The aim of compiling this open access collection of synaptotagmin and synaptotagmin-like sequences, is that its use may lead to greater understanding of the biological function of these proteins in an evolutionary context.

Results

46 metazoan genomes were examined and their complement of Syt, Esyt, Rph3a, Rph3al, Doc2 and Dblc2 genes identified. All of the sequences were compared, named, then examined in detail. Esyt genes were formerly named Fam62. The species in this collection are Trichoplax, Nematostella, Capitella, Helobdella, Lottia, Ciona, Strongylocentrotus, Branchiostoma, Ixodes, Daphnia, Acyrthosiphon, Tribolium, Nasonia, Apis, Anopheles, Drosophila, Caenorhabditis, Takifugu, Tetraodon, Gasterosteus, Oryzias, Danio, Xenopus, Anolis, Gallus, Taeniopygia,Ornithorhynchus, Monodelphis, Mus and Homo. All of the data described in this paper is available as additional files.

Conclusions

Only a subset of synaptotagmin proteins appear able to function as calcium triggers. Syt1, Syt7 and Syt9 are ancient conserved synaptotagmins of this type. Some animals carry extensive repertoires of synaptotagmin genes. Other animals of no less complexity, carry only a small repertoire. Current understanding does not explain why this is so. The biological roles of many synaptotagmins remain to be understood. This collection of genes offers prospects for fruitful speculation about the functional roles of the synaptotagmin repertoires of different animals and includes a great range of biological complexity. With reference to this gene collection, functional relationships among Syt, Esyt, Rph3a, Rph3al, Doc2 and Dblc2 genes, which encode similar proteins, can better be assessed in future.

Abstract

Background

Excessive fat deposition is one of the largest problems faced by salmon aquaculture industries, leading to production losses due to high volume of adipose tissue offal. In addition, increased lipid accumulation may impose considerable stress on adipocytes leading to adipocyte activation and production and secretion of inflammatory mediators, as observed in mammals.

Results

Microarray and qPCR analyses were performed to follow transcriptome changes during adipogenesis in the primary culture of adipose stromo-vascular fraction (aSVF) of Atlantic salmon. Cellular heterogeneity decreased by confluence as evidenced by the down-regulation of markers of osteo/chondrogenic, myogenic, immune and vasculature lineages. Transgelin (TAGLN), a marker of the multipotent pericyte, was prominently expressed around confluence while adipogenic PPARγ was up-regulated already in subconfluent cells. Proliferative activity and subsequent cell cycle arrest were reflected in the fluctuations of pro- and anti-mitotic regulators. Marked regulation of genes involved in lipid and glucose metabolism and pathways producing NADPH and glycerol-3-phosphate (G3P) was seen during the terminal differentiation, also characterised by diverse stress responses. Activation of the glutathione and thioredoxin antioxidant systems and changes in the iron metabolism suggested the need for protection against oxidative stress. Signs of endoplasmic reticulum (ER) stress and unfolded protein response (UPR) occured in parallel with the increased lipid droplet (LD) formation and production of secretory proteins (adipsin, visfatin). The UPR markers XBP1 and ATF6 were induced together with genes involved in ubiquitin-proteasome and lysosomal proteolysis. Concurrently, translation was suppressed as evidenced by the down-regulation of genes encoding elongation factors and components of the ribosomal machinery. Notably, expression changes of a panel of genes that belong to different immune pathways were seen throughout adipogenesis. The induction of AP1 (Jun, Fos), which is a master regulator of stress responses, culminated by the end of adipogenesis, concurrent with the maximal observed lipid deposition.

Conclusions

Our data point to an intimate relationship between metabolic regulation and immune responses in white adipocytes of a cold-blooded vertebrate. Stress imposed on adipocytes by LD formation and expansion is prominently reflected in the ER compartment and the activated UPR response could have an important role at visceral obesity in fish.

Abstract

Background

Chronic myelogenous leukemia (CML) results from the neoplastic transformation of a haematopoietic stem cell. The hallmark genetic abnormality of CML is a chimeric BCR/ABL1 fusion gene resulting from the Philadelphia chromosome rearrangement t(9;22)(q34;q11). Clinical and laboratory studies indicate that the BCR/ABL1 fusion protein is essential for initiation, maintenance and progression of CML, yet the event(s) driving the transformation from chronic phase to blast phase are poorly understood.

Results

Here we report multiple genome aberrations in a collection of 78 CML and 14 control samples by oligonucleotide array comparative genomic hybridization. We found a unique signature of genome deletions within the immunoglobulin heavy chain (IGH) and T cell receptor regions (TCR), frequently accompanied by concomitant loss of sequences within the short arm regions of chromosomes 7 and 9, including IKZF1, HOXA7, CDKN2A/2B, MLLT3, IFNA/B, RNF38, PAX5, JMJD2C and PDCD1LG2 genes.

Conclusions

None of these genome losses were detected in any of the CML samples with myeloid transformation, chronic phase or controls, indicating that their presence is obligatory for the development of a malignant clone with a lymphoid phenotype. Notably, the coincidental deletions at IGH and TCR regions appear to precede the loss of IKZF1 and/or p16 genes in CML indicating a possible involvement of RAG in these deletions.

Abstract

Background

Peach fruit undergoes a rapid softening process that involves a number of metabolic changes. Storing fruit at low temperatures has been widely used to extend its postharvest life. However, this leads to undesired changes, such as mealiness and browning, which affect the quality of the fruit. In this study, a 2-D DIGE approach was designed to screen for differentially accumulated proteins in peach fruit during normal softening as well as under conditions that led to fruit chilling injury.

Results

The analysis allowed us to identify 43 spots -representing about 18% of the total number analyzed- that show statistically significant changes. Thirty-nine of the proteins could be identified by mass spectrometry. Some of the proteins that changed during postharvest had been related to peach fruit ripening and cold stress in the past. However, we identified other proteins that had not been linked to these processes. A graphical display of the relationship between the differentially accumulated proteins was obtained using pairwise average-linkage cluster analysis and principal component analysis. Proteins such as endopolygalacturonase, catalase, NADP-dependent isocitrate dehydrogenase, pectin methylesterase and dehydrins were found to be very important for distinguishing between healthy and chill injured fruit. A categorization of the differentially accumulated proteins was performed using Gene Ontology annotation. The results showed that the 'response to stress', 'cellular homeostasis', 'metabolism of carbohydrates' and 'amino acid metabolism' biological processes were affected the most during the postharvest.

Conclusions

Using a comparative proteomic approach with 2-D DIGE allowed us to identify proteins that showed stage-specific changes in their accumulation pattern. Several proteins that are related to response to stress, cellular homeostasis, cellular component organization and carbohydrate metabolism were detected as being differentially accumulated. Finally, a significant proportion of the proteins identified had not been associated with softening, cold storage or chilling injury-altered fruit before; thus, comparative proteomics has proven to be a valuable tool for understanding fruit softening and postharvest.

Abstract

Background

Chemical communication in mammals involves globular lipocalins that protect and transport pheromones during their passage out of the body. Efficient communication via this protein - pheromone complex is essential for triggering multiple responses including aggression, mate choice, copulatory behaviour, and onset and synchronization of oestrus. The roles of lipocalins in communication were studied in many organisms and especially in mice (i.e. Mus musculus domesticus) which excrete Major Urinary Proteins (Mup) in excessive amounts in saliva and urine. Other mammals, however, often lack the genes for Mups or their expression is very low. Therefore, we aimed at characterization of candidate lipocalins in Myodes glareolus which are potentially linked to chemical communication. One of them is Aphrodisin which is a unique lipocalin that was previously described from hamster vaginal discharge and is known to carry pheromones stimulating copulatory behaviour in males.

Results

Here we show that Aphrodisin-like proteins exist in other species, belong to a group of Odorant Binding Proteins (Obp), and contrary to the expression of Aphrodisin only in hamster genital tract and parotid glands of females, we have detected these transcripts in both sexes of M. glareolus with the expression confirmed in various tissues including prostate, prepucial and salivary glands, liver and uterus. On the level of mRNA, we have detected three different gene variants. To assess their relevance for chemical communication we investigated the occurrence of particular proteins in saliva, urine and vaginal discharge. On the protein level we confirmed the presence of Obp2 and Obp3 in both saliva and urine. Appropriate bands in the range of 17-20 kDa from vaginal discharge were, however, beyond the MS detection limits.

Conclusion

Our results demonstrate that three novel Obps (Obp1, Obp2, and Obp3) are predominant lipocalins in Myodes urine and saliva. On the protein level we have detected further variants and thus we assume that similarly as Major Urinary Proteins in mice, these proteins may be important in chemical communication in this Cricetid rodent.

Abstract

Background

The antenno-maxilary complex (AMC) forms the chemosensory system of the Drosophila larva and is involved in gustatory and olfactory perception. We have previously shown that a mutant allele of the homeodomain transcription factor Prospero (prosVoila1, V1), presents several developmental defects including abnormal growth and altered taste responses. In addition, many neural tracts connecting the AMC to the central nervous system (CNS) were affected. Our earlier reports on larval AMC did not argue in favour of a role of pros in cell fate decision, but strongly suggested that pros could be involved in the control of other aspect of neuronal development. In order to identify these functions, we used microarray analysis of larval AMC and CNS tissue isolated from the wild type, and three other previously characterised prospero alleles, including the V1 mutant, considered as a null allele for the AMC.

Results

A total of 17 samples were first analysed with hierarchical clustering. To determine those genes affected by loss of pros function, we calculated a discriminating score reflecting the differential expression between V1 mutant and other pros alleles. We identified a total of 64 genes in the AMC. Additional manual annotation using all the computed information on the attributed role of these genes in the Drosophila larvae nervous system, enabled us to identify one functional category of potential Prospero target genes known to be involved in neurite outgrowth, synaptic transmission and more specifically in neuronal connectivity remodelling. The second category of genes found to be differentially expressed between the null mutant AMC and the other alleles concerned the development of the sensory organs and more particularly the larval olfactory system. Surprisingly, a third category emerged from our analyses and suggests an association of pros with the genes that regulate autophagy, growth and insulin pathways. Interestingly, EGFR and Notch pathways were represented in all of these three functional categories. We now propose that Pros could perform all of these different functions through the modulation of these two antagonistic and synergic pathways.

Conclusions

The current data contribute to the clarification of the prospero function in the larval AMC and show that pros regulates different function in larvae as compared to those controlled by this gene in embryos. In the future, the possible mechanism by which Pros could achieve its function in the AMC will be explored in detail.

Abstract

Background

Tissue gene expression is generally regulated by multiple transcription factors (TFs). A major first step toward understanding how tissues achieve their specificity is to identify, at the genome scale, interacting TFs regulating gene expression in different tissues. Despite previous discoveries, the mechanisms that control tissue gene expression are not fully understood.

Results

We have integrated a function conservation approach, which is based on evolutionary conservation of biological function, and genes with highest expression level in human tissues to predict TF pairs controlling tissue gene expression. To this end, we have identified 2549 TF pairs associated with a certain tissue. To find interacting TFs controlling tissue gene expression in a broad spatial and temporal manner, we looked for TF pairs common to the same type of tissues and identified 379 such TF pairs, based on which TF-TF interaction networks were further built. We also found that tissue-specific TFs may play an important role in recruiting non-tissue-specific TFs to the TF-TF interaction network, offering the potential for coordinating and controlling tissue gene expression across a variety of conditions.

Conclusion

The findings from this study indicate that tissue gene expression is regulated by large sets of interacting TFs either on the same promoter of a gene or through TF-TF interaction networks.

Abstract

Background

Saliva of adult female mosquitoes help sugar and blood feeding by providing enzymes and polypeptides that help sugar digestion, control microbial growth and counteract their vertebrate host hemostasis and inflammation. Mosquito saliva also potentiates the transmission of vector borne pathogens, including arboviruses. Culex tarsalis is a bird feeding mosquito vector of West Nile Virus closely related to C. quinquefasciatus, a mosquito relatively recently adapted to feed on humans, and the only mosquito of the genus Culex to have its sialotranscriptome so far described.

Results

A total of 1,753 clones randomly selected from an adult female C. tarsalis salivary glands (SG) cDNA library were sequenced and used to assemble a database that yielded 809 clusters of related sequences, 675 of which were singletons. Primer extension experiments were performed in selected clones to further extend sequence coverage, allowing for the identification of 283 protein sequences, 80 of which code for putative secreted proteins.

Conclusion

Comparison of the C. tarsalis sialotranscriptome with that of C. quinquefasciatus reveals accelerated evolution of salivary proteins as compared to housekeeping proteins. The average amino acid identity among salivary proteins is 70.1%, while that for housekeeping proteins is 91.2% (P < 0.05), and the codon volatility of secreted proteins is significantly higher than those of housekeeping proteins. Several protein families previously found exclusive of mosquitoes, including only in the Aedes genus have been identified in C. tarsalis. Interestingly, a protein family so far unique to C. quinquefasciatus, with 30 genes, is also found in C. tarsalis, indicating it was not a specific C. quinquefasciatus acquisition in its evolution to optimize mammal blood feeding.

Abstract

Background

Non-coding RNAs (ncRNAs) play a crucial role in the intricate regulation of bacterial gene expression, allowing bacteria to quickly adapt to changing environments. In the past few years, a growing number of regulatory RNA elements have been predicted by computational methods, mostly in well-studied γ-proteobacteria but lately in several α-proteobacteria as well. Here, we have compared an extensive compilation of these non-coding RNA predictions to intergenic expression data of a whole-genome high-resolution tiling array in the soil-dwelling α-proteobacterium Rhizobium etli.

Results

Expression of 89 candidate ncRNAs was detected, both on the chromosome and on the six megaplasmids encompassing the R. etli genome. Of these, 11 correspond to functionally well characterized ncRNAs, 12 were previously identified in other α-proteobacteria but are as yet uncharacterized and 66 were computationally predicted earlier but had not been experimentally identified and were therefore classified as novel ncRNAs. The latter comprise 17 putative sRNAs and 49 putative cis-regulatory ncRNAs. A selection of these candidate ncRNAs was validated by RT-qPCR, Northern blotting and 5' RACE, confirming the existence of 4 ncRNAs. Interestingly, individual transcript levels of numerous ncRNAs varied during free-living growth and during interaction with the eukaryotic host plant, pointing to possible ncRNA-dependent regulation of these specialized processes.

Conclusions

Our data support the practical value of previous ncRNA prediction algorithms and significantly expand the list of candidate ncRNAs encoded in the intergenic regions of R. etli and, by extension, of α-proteobacteria. Moreover, we show high-resolution tiling arrays to be suitable tools for studying intergenic ncRNA transcription profiles across the genome. The differential expression levels of some of these ncRNAs may indicate a role in adaptation to changing environmental conditions.

Abstract

Background

Small endogenous non-coding RNAs (sncRNAs) such as small interfering RNA (siRNA), microRNA and other small RNA transcripts are derived from distinct loci in the genome and play critical roles in RNA-mediated gene silencing mechanisms in plants and metazoa. They are approximately 22 nucleotides long; regulate mRNA stability through perfect or imperfect match to the targets. The biological activities of sncRNAs have been related to many biological events, from resistance to microbe infections to cellular differentiation. The development of the zoonotic parasite Schistosoma japonicum parasite includes multiple steps of morphological alterations and biological differentiations, which provide a unique model for studies on the functions of small RNAs. Characterization of the genome-wide transcription of the sncRNAs will be a major step in understanding of the parasite biology. The objective of this study is to investigate the transcriptional profile and potential function of the small non-coding RNAs in the development of S. japanicum.

Results

The endogenous siRNAs were found mainly derived from transposable elements (TE) or transposons and the natural antisense transcripts (NAT). In contrast to other organisms, the TE-derived siRNAs in S. japonicum were more predominant than other sncRNAs including microRNAs (miRNAs). Further, there were distinct length and 3'end variations in the sncRNAs, which were associated with the developmental differentiation of the parasite. Among the identified miRNA transcripts, there were 38 unique to S. japonicum and 16 that belonged to 13 miRNA families are common to other metazoan lineages. These miRNAs were either ubiquitously expressed, or they exhibited specific expression patterns related to the developmental stages or sex. Genes that encoded miRNAs are mainly located in clusters within the genome of S. japonicum. However, genes within one cluster could be differentially transcribed, which suggested that individual genes might be regulated by distinct mechanisms during parasite development.

Conclusions

Many miRNA and endogenous siRNA transcripts were identified in S. japonicum and the amount of siRNA was at least 4.4 and 1.6 times more than that of miRNA in both schistosomulum and adult worm stages respectively. SiRNAs are mainly derived from transposable elements (or transposons); while natural antisense transcripts (NAT)-derived siRNAs were much less. A majority of miRNA transcripts identified in the parasite were species-specific and the expression of certain miRNAs was found developmentally regulated. Both miRNA and siRNAs are potentially important regulators in the development of schistosomal parasites.

Abstract

Background

As uricoletic animals, chickens produce cleidoic eggs, which are self-contained bacteria-resistant biological packages for extra-uterine development of the chick embryo. The eggshell constitutes a natural physical barrier against bacterial penetration if it forms correctly and remains intact. The eggshell's remarkable mechanical properties are due to interactions among mineral components and the organic matrix proteins. The purpose of our study was to identify novel eggshell proteins by examining the transcriptome of the uterus during calcification of the eggshell. An extensive bioinformatic analysis on genes over-expressed in the uterus allowed us to identify novel eggshell proteins that contribute to the egg's natural defenses.

Results

Our 14 K Del-Mar Chicken Integrated Systems microarray was used for transcriptional profiling in the hen's uterus during eggshell deposition. A total of 605 transcripts were over-expressed in the uterus compared with the magnum or white isthmus across a wide range of abundance (1.1- to 79.4-fold difference). The 605 highly-expressed uterine transcripts correspond to 469 unique genes, which encode 437 different proteins. Gene Ontology (GO) analysis was used for interpretation of protein function. The most over-represented GO terms are related to genes encoding ion transport proteins, which provide eggshell mineral precursors. Signal peptide sequence was found for 54 putative proteins secreted by the uterus during eggshell formation. Many functional proteins are involved in calcium binding or biomineralization--prerequisites for interacting with the mineral phase during eggshell fabrication. While another large group of proteins could be involved in proper folding of the eggshell matrix. Many secreted uterine proteins possess antibacterial properties, which would protect the egg against microbial invasion. A final group includes proteases and protease inhibitors that regulate protein activity in the acellular uterine fluid where eggshell formation takes place.

Conclusions

Our original study provides the first detailed description of the chicken uterus transcriptome during formation of the eggshell. We have discovered a cache of about 600 functional genes and identified a large number of encoded proteins secreted into uterine fluid for fabrication of the eggshell and chemical protection of the egg. Some of these uterine genes could prove useful as biological markers for genetic improvement of phenotypic traits (i.e., egg and eggshell quality).

Abstract

Background

Histone acetyltransferase enzymes (HATs) are implicated in regulation of transcription. HATs from different families may overlap in target and substrate specificity.

Results

We isolated the elp3⁺ gene encoding the histone acetyltransferase subunit of the Elongator complex in fission yeast and characterized the phenotype of an Δelp3 mutant. We examined genetic interactions between Δelp3 and two other HAT mutants, Δmst2 and Δgcn5 and used whole genome microarray analysis to analyze their effects on gene expression.

Conclusions

Comparison of phenotypes and expression profiles in single, double and triple mutants indicate that these HAT enzymes have overlapping functions. Consistent with this, overlapping specificity in histone H3 acetylation is observed. However, there is no evidence for overlap with another HAT enzyme, encoded by the essential mst1⁺ gene.

Abstract

Background

Recent studies have demonstrated that non-protein-coding RNAs (npcRNAs/ncRNAs) play important roles during eukaryotic development, species evolution, and in the etiology of disease. Rhesus macaques are the most widely used primate model in both biomedical research and primate evolutionary studies. However, most reports on these animals focus on the functional roles of protein-coding sequences, whereas very little is known about macaque ncRNAs.

Results

In the present study, we performed the first systematic profiling of intermediate-size ncRNAs (50 to 500 nt) from the rhesus monkey by constructing a cDNA library. We identified 117 rhesus monkey ncRNAs, including 80 small nucleolar RNAs (snoRNAs), 29 other types of known RNAs (snRNAs, Y RNA, and others), and eight unclassified ncRNAs. Comparative genomic analysis and northern blot hybridizations demonstrated that some snoRNAs were lineage- or species-specific. Paralogous sequences were found for most rhesus monkey snoRNAs, the expression of which might be attributable to extensive duplication within the rhesus monkey genome. Further investigation of snoRNA flanking sequences showed that some rhesus monkey snoRNAs are retrogenes derived from L1-mediated integration. Finally, phylogenetic analysis demonstrated that birds and primates share some snoRNAs and host genes thereof, suggesting that both the relevant host genes and the snoRNAs contained therein may be inherited from a common ancestor. However, some rhesus monkey snoRNAs hosted by non-ribosome-related genes appeared after the evolutionary divergence between birds and mammals.

Conclusions

We provide the first experimentally-derived catalog of rhesus monkey ncRNAs and uncover some interesting genomic and evolutionary features. These findings provide important information for future functional characterization of snoRNAs during primate evolution.

Abstract

Background

The availability of two genetically very similar cell lines (A and B) derived from the laboratory isolate Entamoeba histolytica HM-1:IMSS, which differ in their virulence properties, provides a powerful tool for identifying pathogenicity factors of the causative agent of human amoebiasis. Cell line A is incapable inducing liver abscesses in gerbils, whereas interaction with cell line B leads to considerable abscess formation. Phenotypic characterization of both cell lines revealed that trophozoites from the pathogenic cell line B have a larger cell size, an increased growth rate in vitro, an increased cysteine peptidase activity and higher resistance to nitric oxide stress. To find proteins that may serve as virulence factors, the proteomes of both cell lines were previously studied, resulting in the identification of a limited number of differentially synthesized proteins. This study aims to identify additional genes, serving as virulence factors, or virulence markers.

Results

To obtain a comprehensive picture of the differences between the cell lines, we compared their transcriptomes using an oligonucleotide-based microarray and confirmed findings with quantitative real-time PCR. Out of 6242 genes represented on the array, 87 are differentially transcribed (≥two-fold) in the two cell lines. Approximately 50% code for hypothetical proteins. Interestingly, only 19 genes show a five-fold or higher differential expression. These include three rab7 GTPases, which were found with a higher abundance in the non-pathogenic cell line A. The aig1-like GTPasesare of special interest because the majority of them show higher levels of transcription in the pathogenic cell line B. Only two molecules were found to be differentially expressed between the two cell lines in both this study and our previous proteomic approach.

Conclusions

In this study we have identified a defined set of genes that are differentially transcribed between the non-pathogenic cell line A and the pathogenic cell line B of E. histolytica. The identification of transcription profiles unique for amoebic cell lines with pathogenic phenotypes may help to elucidate the transcriptional framework of E. histolytica pathogenicity and serve as a basis for identifying transcriptional markers and virulence factors.

Abstract

Background

Tungsten carbide (WC) and tungsten carbide cobalt (WC-Co) nanoparticles are of occupational health relevance because of the increasing usage in hard metal industries. Earlier studies showed an enhanced toxic potential for WC-Co compared to WC or cobalt ions alone. Therefore, we investigated the impact of these particles, compared to cobalt ions applied as CoCl₂, on the global gene expression level in human keratinocytes (HaCaT) in vitro.

Results

WC nanoparticles exerted very little effects on the transcriptomic level after 3 hours and 3 days of exposure. In contrast, WC-Co nanoparticles caused significant transcriptional changes that were similar to those provoked by CoCl₂. However, CoCl₂ exerted even more pronounced changes in the transcription patterns. Gene set enrichment analyses revealed that the differentially expressed genes were related to hypoxia response, carbohydrate metabolism, endocrine pathways, and targets of several transcription factors. The role of the transcription factor HIF1 (hypoxia inducible factor 1) is particularly highlighted and aspects of downstream events as well as the role of other transcription factors related to cobalt toxicity are considered.

Conclusions

This study provides extensive data useful for the understanding of nanoparticle and cobalt toxicity. It shows that WC nanoparticles caused low transcriptional responses while WC-Co nanoparticles are able to exert responses similar to that of free cobalt ions, particularly the induction of hypoxia-like effects via interactions with HIF1α in human keratinocytes. However, the enhanced toxicity of WC-Co particles compared to CoCl₂ could not be explained by differences in gene transcription.

Abstract

Background

Helicobacter canadensis is an emerging human pathogen and zoonotic agent. The genome of H. canadensis was sequenced previously and determined to contain 29 annotated coding regions associated with homopolymeric tracts.

Results

Twenty-one of the repeat-associated coding regions were determined to be potentially transcriptionally or translationally phase variable. In each case the homopolymeric tract was within the predicted promoter region or at the 5' end of the coding region, respectively. However, eight coding sequences were identified with simple sequence repeats toward the 3' end of the open reading frame. In these cases, the repeat tract would be too far into the coding region to be mediating translational phase variation. All of the 29 coding region-associated homopolymeric tracts display variability in tract length in the sequencing read data.

Conclusions

Twenty-nine coding regions have been identified in the genome sequence of Helicobacter canadensis strain NCTC13241 that show variations in homopolymeric tract length in the bacterial population, indicative of phase variation. Five of these are potentially associated with promoter regions, which would lead to transcriptional phase variation. Translational phase variation usually switches expression of a gene ON and OFF due to the repeat region being located sufficiently close to the initiation codon for the resulting frame-shift to lead to a premature termination codon and stop the translation of the protein. Sixteen of the 29 coding regions have homopolymeric tracts characteristic of translational phase variation. For eight coding sequences with repeats located later in the reading frame, changes in the repeat tract length would alter the protein sequence at the C-terminus but not stop the expression of the protein. This mechanism of C-terminal phase variation has implications for stochastic switching of protein sequence in bacterial species that already undergo transcriptional and translational phase variation.

Abstract

Background

Plants engineered for abiotic stress tolerance may soon be commercialized. The engineering of these plants typically involves the manipulation of complex multigene networks and may therefore have a greater potential to introduce pleiotropic effects than the simple monogenic traits that currently dominate the plant biotechnology market. While research on unintended effects in transgenic plant systems has been instrumental in demonstrating the substantial equivalence of many transgenic plant systems, it is essential that such analyses be extended to transgenic plants engineered for stress tolerance. Drought-tolerant Arabidopsis thaliana were engineered through overexpression of the transcription factor ABF3 in order to investigate unintended pleiotropic effects. In order to eliminate position effects, the Cre/lox recombination system was used to create control plant lines that contain identical T-DNA insertion sites but with the ABF3 transgene excised. This additionally allowed us to determine if Cre recombinase can cause unintended effects that impact the transcriptome.

Results

Microarray analysis of control plant lines that underwent Cre-mediated excision of the ABF3 transgene revealed only two genes that were differentially expressed in more than one plant line, suggesting that the impact of Cre recombinase on the transcriptome was minimal. In the absence of drought stress, overexpression of ABF3 had no effect on the transcriptome, but following drought stress, differences were observed in the gene expression patterns of plants overexpressing ABF3 relative to control plants. Examination of the functional distribution of the differentially expressed genes revealed strong similarity indicating that unintended pathways were not activated.

Conclusions

The action of ABF3 is tightly controlled in Arabidopsis. In the absence of drought stress, ectopic activation of drought response pathways does not occur. In response to drought stress, overexpression of ABF3 results in a reprogramming of the drought response, which is characterized by changes in the timing or strength of expression of some drought response genes, without activating any unexpected gene networks. These results illustrate that important gene networks are highly regulated in Arabidopsis and that engineering stress tolerance may not necessarily cause extensive changes to the transcriptome.

Abstract

Background

Genes in bacteria may be organised into operons, leading to strict co-expression of the genes that participate in the same operon. However, comparisons between different bacterial genomes have shown that much of the operon structure is dynamic on an evolutionary time scale. This indicates that there are opposing effects influencing the tendency for operon formation, and these effects may be reflected in properties like evolutionary rate, complex formation, metabolic pathways and gene fusion.

Results

We have used multi-species protein-protein comparisons to generate a high-quality set of genes that are persistent in bacterial genomes (i.e. they have close to universal distribution). We have analysed these genes with respect to operon participation and important functional properties, including evolutionary rate and protein-protein interactions.

Conclusions

Genes for ribosomal proteins show a very slow rate of evolution. This is consistent with a strong tendency for the genes to participate in operons and for their proteins to be involved in essential and well defined complexes. Persistent genes for non-ribosomal proteins can be separated into two classes according to tendency to participate in operons. Those with a strong tendency for operon participation make proteins with fewer interaction partners that seem to participate in relatively static complexes and possibly linear pathways. Genes with a weak tendency for operon participation tend to produce proteins with more interaction partners, but possibly in more dynamic complexes and convergent pathways. Genes that are not regulated through operons are therefore more evolutionary constrained than the corresponding operon-associated genes and will on average evolve more slowly.

Abstract

Background

Glutathione S-transferases (GSTs) are the ubiquitous enzymes that play a key role in cellular detoxification. Although several GSTs have been identified and characterized in various plant species, the knowledge about their role in developmental processes and response to various stimuli is still very limited. In this study, we report genome-wide identification, characterization and comprehensive expression analysis of members of GST gene family in crop plant rice, to reveal their function(s).

Results

A systematic analysis revealed the presence of at least 79 GST genes in the rice genome. Phylogenetic analysis grouped GST proteins into seven classes. Sequence analysis together with the organization of putative motifs indicated the potential diverse functions of GST gene family members in rice. The tandem gene duplications have contributed a major role in expansion of this gene family. Microarray data analysis revealed tissue-/organ- and developmental stage-specific expression patterns of several rice GST genes. At least 31 GST genes showed response to plant hormones auxin and cytokinin. Furthermore, expression analysis showed the differential expression of quite a large number of GST genes during various abiotic stress (20), arsenate stress (32) and biotic stress (48) conditions. Many of the GST genes were commonly regulated by developmental processes, hormones, abiotic and biotic stresses.

Conclusion

The transcript profiling suggests overlapping and specific role(s) of GSTs during various stages of development in rice. Further, the study provides evidence for the role of GSTs in mediating crosstalk between various stress and hormone response pathways and represents a very useful resource for functional analysis of selected members of this family in rice.

Abstract

Background

Gene regulation at transcript level can provide a good indication of the complex signaling mechanisms underlying physiological and pathological processes. Transcriptomic methods such as microarray and quantitative real-time PCR require stable reference genes for accurate normalization of gene expression. Some but not all studies have shown that housekeeping genes (HGKs), β-actin (ACTB) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH), which are routinely used for normalization, may vary significantly depending on the cell/tissue type and experimental conditions. It is currently unclear if these genes are stably expressed in cells undergoing drastic morphological changes during neuronal differentiation. Recent meta-analysis of microarray datasets showed that some but not all of the ribosomal protein genes are stably expressed. To test the hypothesis that some ribosomal protein genes can serve as reference genes for neuronal differentiation, a genome-wide analysis was performed and putative reference genes were identified based on stability of expressions. The stabilities of these potential reference genes were then analyzed by reverse transcription quantitative real-time PCR in six differentiation conditions.

Results

Twenty stably expressed genes, including thirteen ribosomal protein genes, were selected from microarray analysis of the gene expression profiles of GDNF and NGF induced differentiation of PC12 cells. The expression levels of these candidate genes as well as ACTB and GAPDH were further analyzed by reverse transcription quantitative real-time PCR in PC12 cells differentiated with a variety of stimuli including NGF, GDNF, Forskolin, KCl and ROCK inhibitor, Y27632. The performances of these candidate genes as stable reference genes were evaluated with two independent statistical approaches, geNorm and NormFinder.

Conclusions

The ribosomal protein genes, RPL19 and RPL29, were identified as suitable reference genes during neuronal differentiation of PC12 cells, regardless of the type of differentiation conditions. The combination of these two novel reference genes, but not the commonly used HKG, GAPDH, allows robust and accurate normalization of differentially expressed genes during PC12 differentiation.

Abstract

Background

The increasing interest in small non-coding RNAs (ncRNAs) such as microRNAs (miRNAs), small interfering RNAs (siRNAs) and Piwi-interacting RNAs (piRNAs) and recent advances in sequencing technology have yielded large numbers of short (18-32 nt) RNA sequences from different organisms, some of which are derived from small nucleolar RNAs (snoRNAs) and transfer RNAs (tRNAs). We observed that these short ncRNAs frequently cover the entire length of annotated snoRNAs or tRNAs, which suggests that other loci specifying similar ncRNAs can be identified by clusters of short RNA sequences.

Results

We combined publicly available datasets of tens of millions of short RNA sequence tags from Drosophila melanogaster, and mapped them to the Drosophila genome. Approximately 6 million perfectly mapping sequence tags were then assembled into 521,302 tag-contigs (TCs) based on tag overlap. Most transposon-derived sequences, exons and annotated miRNAs, tRNAs and snoRNAs are detected by TCs, which show distinct patterns of length and tag-depth for different categories. The typical length and tag-depth of snoRNA-derived TCs was used to predict 7 previously unrecognized box H/ACA and 26 box C/D snoRNA candidates. We also identified one snRNA candidate and 86 loci with a high number of tags that are yet to be annotated, 7 of which have a particular 18mer motif and are located in introns of genes involved in development. A subset of new snoRNA candidates and putative ncRNA candidates was verified by Northern blot.

Conclusions

In this study, we have introduced a new approach to identify new members of known classes of ncRNAs based on the features of TCs corresponding to known ncRNAs. A large number of the identified TCs are yet to be examined experimentally suggesting that many more novel ncRNAs remain to be discovered.

Abstract

Background

Proteins do not act in isolation; they frequently act together in protein complexes to carry out concerted cellular functions. The evolution of complexes is poorly understood, especially in organisms other than yeast, where little experimental data has been available.

Results

We generated accurate, high coverage datasets of protein complexes for E. coli and yeast in order to study differences in the evolution of complexes between these two species. We show that substantial differences exist in how complexes have evolved between these organisms. A previously proposed model of complex evolution identified complexes with cores of interacting homologues. We support findings of the relative importance of this mode of evolution in yeast, but find that it is much less common in E. coli. Additionally it is shown that those homologues which do cluster in complexes are involved in eukaryote-specific functions. Furthermore we identify correlated pairs of non-homologous domains which occur in multiple protein complexes. These were identified in both yeast and E. coli and we present evidence that these too may represent complex cores in yeast but not those of E. coli.

Conclusions

Our results suggest that there are differences in the way protein complexes have evolved in E. coli and yeast. Whereas some yeast complexes have evolved by recruiting paralogues, this is not apparent in E. coli. Furthermore, such complexes are involved in eukaryotic-specific functions. This implies that the increase in gene family sizes seen in eukaryotes in part reflects multiple family members being used within complexes. However, in general, in both E. coli and yeast, homologous domains are used in different complexes.

Abstract

Background

More and more completely sequenced fungal genomes are becoming available and many more sequencing projects are in progress. This deluge of data should improve our knowledge of the various primary and secondary metabolisms of Fungi, including their synthesis of useful compounds such as antibiotics or toxic molecules such as mycotoxins. Functional annotation of many fungal genomes is imperfect, especially of genes encoding enzymes, so we need dedicated tools to analyze their metabolic pathways in depth.

Description

FUNGIpath is a new tool built using a two-stage approach. Groups of orthologous proteins predicted using complementary methods of detection were collected in a relational database. Each group was further mapped on to steps in the metabolic pathways published in the public databases KEGG and MetaCyc. As a result, FUNGIpath allows the primary and secondary metabolisms of the different fungal species represented in the database to be compared easily, making it possible to assess the level of specificity of various pathways at different taxonomic distances. It is freely accessible at http://www.fungipath.u-psud.fr.

Conclusions

As more and more fungal genomes are expected to be sequenced during the coming years, FUNGIpath should help progressively to reconstruct the ancestral primary and secondary metabolisms of the main branches of the fungal tree of life and to elucidate the evolution of these ancestral fungal metabolisms to various specific derived metabolisms.

Abstract

Background

Phosphatidylinositol 3-phosphate is involved in regulation of several key cellular processes, mainly endocytosis, signaling, nuclear processes, cytoskeletal remodelling, cell survival, membrane trafficking, phagosome maturation and autophagy. In most cases effector proteins bind to this lipid, using either FYVE or PX domain. These two domains are distributed amongst varied life forms such as virus, protists, fungi, viridiplantae and metazoa. As the binding ligand is identical for both domains, the goal of this study was to understand if there is any selectivity for either of these domains in different taxa. Further, to understand the different cellular functions that these domains may be involved in, we analyzed the taxonomic distribution of additional domains that associate with FYVE and PX.

Results

There is selectivity for either FYVE or PX in individual genomes where both domains are present. Fungi and metazoa encode more PX, whereas streptophytes in viridiplantae encode more FYVE. Excess of FYVE in streptophytes results from proteins containing RCC1and DZC domains and FYVE domains in these proteins have a non-canonical ligand-binding site. Within a taxonomic group the selected domain associates with a higher number of other domains and is thus expected to discharge a larger number of cellular functions. Also, while certain associated domains are present in all taxonomic groups, most of them are unique to a specific group indicating that while certain common functions are discharged by these domains in all taxonomic groups, some functions appear to be group specific.

Conclusions

Although both FYVE and PX bind to PtdIns(3)P, genomes of different taxa show distinct selectivity of encoding either of the two. Higher numbers of taxonomic group specific domains co-occur with the more abundant domain (FYVE/PX) indicating that group-specific rare domain architectures might have emerged to accomplish certain group-specific functions.

Abstract

Background

MicroRNAs (miRNAs) repress target genes at the post-transcriptional level, and function in the development and cell-lineage pathways of host species. Tissue-specific expression of miRNAs is highly relevant to their physiological roles in the corresponding tissues. However, to date, few miRNAs have been spatially identified in the silkworm.

Results

We establish for the first time the spatial expression patterns of nearly 100 miRNAs in multiple normal tissues (organs) of Bombyx mori females and males using microarray and Northern-blotting analyses. In all, only 10 miRNAs were universally distributed (including bmo-let-7 and bmo-bantam), while the majority were expressed exclusively or preferentially in specific tissue types (e.g., bmo-miR-275 and bmo-miR-1). Additionally, we examined the developmental patterns of miRNA expression during metamorphosis of the body wall, silk glands, midgut and fat body. In total, 63 miRNAs displayed significant alterations in abundance in at least 1 tissue during the developmental transition from larvae to pupae (e.g., bmo-miR-263b and bmo-miR-124). Expression patterns of five miRNAs were significantly increased during metamorphosis in all four tissues (e.g., bmo-miR-275 and bmo-miR-305), and two miRNA pairs, bmo-miR-10b-3p/5p and bmo-miR-281-3p/5p, showed coordinate expression.

Conclusions

In this study, we conducted preliminary spatial measurements of several miRNAs in the silkworm. Periods of rapid morphological change were associated with alterations in miRNA expression patterns in the body wall, silk glands, midgut and fat body during metamorphosis. Accordingly, we propose that corresponding ubiquitous or tissue-specific expression of miRNAs supports their critical roles in tissue specification. These results should facilitate future functional analyses.

Abstract

Background

Analysis of fungal genome sequence assemblies reveals that telomeres are poorly represented even though telomeric reads tend to be superabundant. We surmised that the problem might lie in the DNA shearing conditions used to create clone libraries for genome sequencing.

Results

A shotgun strategy was used to sequence and assemble circular and linear cosmid DNAs sheared using conditions typical for a genome project. The DNA sheared in circular form assembled into a single sequence contig. However, the linearized cosmid produced an incomplete assembly because the two DNA termini, though greatly overrepresented in the clone library used for sequencing, were separated from neighboring sequences by gaps of ~1.4 and 1.8 kb. These gap sizes were reduced, but not eliminated, by shearing the linear cosmid into smaller fragments. Mapping of shearing breakpoints revealed a paucity of breaks in the subterminal regions of the linearized cosmid and also near chromosome ends of the fungus Neurospora crassa.

Conclusion

Together, our data indicate that the ends of linear DNA molecules are recalcitrant to hydrodynamic shearing. We propose that this causes DNA termini to be overrepresented in the resulting fragment population but ultimately prevents their incorporation into sequence assemblies.

Abstract

Background

Retrotransposons are key players in the evolution of eukaryotic genomes. Moreover, it is now known that some retrotransposon classes, like the abundant and plant-specific Sireviruses, have intriguingly distinctive host preferences. Yet, it is largely unknown if this bias is supported by different genome structures.

Results

We performed sensitive comparative analysis of the genomes of a large set of Ty1/copia retrotransposons. We discovered that Sireviruses are unique among Pseudoviridae in that they constitute an ancient genus characterized by vastly divergent members, which however contain highly conserved motifs in key non-coding regions: multiple polypurine tract (PPT) copies cluster upstream of the 3' long terminal repeat (3'LTR), of which the terminal PPT tethers to a distinctive attachment site and is flanked by a precisely positioned inverted repeat. Their LTRs possess a novel type of repeated motif (RM) defined by its exceptionally high copy number, symmetry and core CGG-CCG signature. These RM boxes form CpG islands and lie a short distance upstream of a conserved promoter region thus hinting towards regulatory functions. Intriguingly, in the envelope-containing Sireviruses additional boxes cluster at the 5' vicinity of the envelope. The 5'LTR/internal domain junction and a polyC-rich integrase signal are also highly conserved domains of the Sirevirus genome.

Conclusions

Our comparative analysis of retrotransposon genomes using advanced in silico methods highlighted the unique genome organization of Sireviruses. Their structure may dictate a life cycle with different regulation and transmission strategy compared to other Pseudoviridae, which may contribute towards their pattern of distribution within and across plants.

Abstract

Background

Corynebacterium aurimucosum is a slightly yellowish, non-lipophilic, facultative anaerobic member of the genus Corynebacterium and predominantly isolated from human clinical specimens. Unusual black-pigmented variants of C. aurimucosum (originally named as C. nigricans) continue to be recovered from the female urogenital tract and they are associated with complications during pregnancy. C. aurimucosum ATCC 700975 (C. nigricans CN-1) was originally isolated from a vaginal swab of a 34-year-old woman who experienced a spontaneous abortion during month six of pregnancy. For a better understanding of the physiology and lifestyle of this potential urogenital pathogen, the complete genome sequence of C. aurimucosum ATCC 700975 was determined.

Results

Sequencing and assembly of the C. aurimucosum ATCC 700975 genome yielded a circular chromosome of 2,790,189 bp in size and the 29,037-bp plasmid pET44827. Specific gene sets associated with the central metabolism of C. aurimucosum apparently provide enhanced metabolic flexibility and adaptability in aerobic, anaerobic and low-pH environments, including gene clusters for the uptake and degradation of aromatic amines, L-histidine and L-tartrate as well as a gene region for the formation of selenocysteine and its incorporation into formate dehydrogenase. Plasmid pET44827 codes for a non-ribosomal peptide synthetase that plays the pivotal role in the synthesis of the characteristic black pigment of C. aurimucosum ATCC 700975.

Conclusions

The data obtained by the genome project suggest that C. aurimucosum could be both a resident of the human gut and possibly a pathogen in the female genital tract causing complications during pregnancy. Since hitherto all black-pigmented C. aurimucosum strains have been recovered from female genital source, biosynthesis of the pigment is apparently required for colonization by protecting the bacterial cells against the high hydrogen peroxide concentration in the vaginal environment. The location of the corresponding genes on plasmid pET44827 explains why black-pigmented (formerly C. nigricans) and non-pigmented C. aurimucosum strains were isolated from clinical specimens.

Abstract

Background

Even though the process of potato tuber starch biosynthesis is well understood, mechanisms regulating biosynthesis are still unclear. Transcriptome analysis provides valuable information as to how genes are regulated. Therefore, this work aimed at investigating transcriptional regulation of starch biosynthetic genes in leaves and tubers of potato plants under various conditions. More specifically we looked at gene expression diurnally in leaves and tubers, during tuber induction and in tubers growing at different velocities. To determine velocity of potato tuber growth a new method based on X-ray Computed Tomography (X-ray CT) was established.

Results

Comparative transcriptome analysis between leaves and tubers revealed striking similarities with the same genes being differentially expressed in both tissues. In tubers, oscillation of granule bound starch synthase (GBSS) expression) was observed which could be linked to sucrose supply from source leaves. X-ray CT was used to determine time-dependent changes in tuber volume and the growth velocity was calculated. Although there is not a linear correlation between growth velocity and expression of starch biosynthetic genes, there are significant differences between growing and non-growing tubers. Co-expression analysis was used to identify transcription factors positively correlating with starch biosynthetic genes possibly regulating starch biosynthesis.

Conclusion

Most starch biosynthetic enzymes are encoded by gene families. Co-expression analysis revealed that the same members of these gene families are co-regulated in leaves and tubers. This suggests that regulation of transitory and storage starch biosynthesis in leaves and tubers, respectively, is surprisingly similar. X-ray CT can be used to monitor growth and development of belowground organs and allows to link tuber growth to changes in gene expression. Comparative transcriptome analysis provides a useful tool to identify transcription factors possibly involved in the regulation of starch biosynthesis.

Abstract

Background

ClpB-cyt/HSP100 protein acts as chaperone, mediating disaggregation of denatured proteins. Previous studies have shown that ClpB-cyt/HSP100 gene belongs to the group class I Clp ATPase proteins and ClpB-cyt/HSP100 transcript is regulated by heat stress and developmental cues.

Results

Nine ORFs were noted to constitute rice class I Clp ATPases in the following manner: 3 ClpB proteins (ClpB-cyt, Os05g44340; ClpB-m, Os02g08490; ClpB-c, Os03g31300), 4 ClpC proteins (ClpC1, Os04g32560; ClpC2, Os12g12580; ClpC3, Os11g16590; ClpC4, Os11g16770) and 2 ClpD proteins (ClpD1, Os02g32520; ClpD2, Os04g33210). Using the respective signal sequences cloned upstream to GFP/CFP reporter proteins and transient expression studies with onion epidermal cells, evidence is provided that rice ClpB-m and Clp-c proteins are indeed localized to their respective cell locations mitochondria and chloroplasts, respectively. Associated with their diverse cell locations, domain structures of OsClpB-c, OsClpB-m and OsClpB-cyt proteins are noted to possess a high-level conservation. OsClpB-cyt transcript is shown to be enriched at milk and dough stages of seed development. While expression of OsClpB-m was significantly less as compared to its cytoplasmic and chloroplastic counterparts in different tissues, this transcript showed highest heat-induced expression amongst the 3 ClpB proteins. OsClpC1 and OsClpC2 are predicted to be chloroplast-localized as is the case with all known plant ClpC proteins. However, the fact that OsClpC3 protein appears mitochondrial/chloroplastic with equal probability and OsClpC4 a plasma membrane protein reflects functional diversity of this class. Different class I Clp ATPase transcripts were noted to be cross-induced by a host of different abiotic stress conditions. Complementation assays of Δhsp104 mutant yeast cells showed that OsClpB-cyt, OsClpB-m, OsClpC1 and OsClpD1 have significantly positive effects. Remarkably, OsClpD1 gene imparted appreciably high level tolerance to the mutant yeast cells.

Conclusions

Rice class I Clp ATPase gene family is constituted of 9 members. Of these 9, only 3 belonging to ClpB group are heat stress regulated. Distribution of ClpB proteins to different cell organelles indicates that their functioning might be critical in different cell locations. From the complementation assays, OsClpD1 appears to be more effective than OsClpB-cyt protein in rescuing the thermosensitive defect of the yeast ScΔhsp104 mutant cells.

Abstract

Background

Although polyploidy has long been recognized as a major force in the evolution of plants, most of what we know about the genetic consequences of polyploidy comes from the study of crops and model systems. Furthermore, although many polyploid species have formed repeatedly, patterns of genome evolution and gene expression are largely unknown for natural polyploid populations of independent origin. We therefore examined patterns of loss and expression in duplicate gene pairs (homeologs) in multiple individuals from seven natural populations of independent origin of Tragopogon mirus (Asteraceae), an allopolyploid that formed repeatedly within the last 80 years from the diploids T. dubius and T. porrifolius.

Results

Using cDNA-AFLPs, we found differential band patterns that could be attributable to gene silencing, novel expression, and/or maternal/paternal effects between T. mirus and its diploid parents. Subsequent cleaved amplified polymorphic sequence (CAPS) analyses of genomic DNA and cDNA revealed that 20 of the 30 genes identified through cDNA-AFLP analysis showed additivity, whereas nine of the 30 exhibited the loss of one parental homeolog in at least one individual. Homeolog loss (versus loss of a restriction site) was confirmed via sequencing. The remaining gene (ADENINE-DNA GLYCOSYLASE) showed ambiguous patterns in T. mirus because of polymorphism in the diploid parent T. dubius. Most (63.6%) of the homeolog loss events were of the T. dubius parental copy. Two genes, NUCLEAR RIBOSOMAL DNA and GLYCERALDEHYDE-3-PHOSPHATE DEHYDROGENASE, showed differential expression of the parental homeologs, with the T. dubius copy silenced in some individuals of T. mirus.

Conclusions

Genomic and cDNA CAPS analyses indicated that plants representing multiple populations of this young natural allopolyploid have experienced frequent and preferential elimination of homeologous loci. Comparable analyses of synthetic F₁ hybrids showed only additivity. These results suggest that loss of homeologs and changes in gene expression are not the immediate result of hybridization, but are processes that occur following polyploidization, occurring during the early (<40) generations of the young polyploid. Both T. mirus and a second recently formed allopolyploid, T. miscellus, exhibit more homeolog losses than gene silencing events. Furthermore, both allotetraploids undergo biased loss of homeologs contributed by their shared diploid parent, T. dubius. Further studies are required to assess whether the results for the 30 genes so far examined are representative of the entire genome.

Abstract

Background

Trichomonas vaginalis is the most common non-viral human sexually transmitted pathogen and importantly, contributes to facilitating the spread of HIV. Yet very little is known about its surface and secreted proteins mediating interactions with, and permitting the invasion and colonisation of, the host mucosa. Initial annotations of T. vaginalis genome identified a plethora of candidate extracellular proteins.

Results

Data mining of the T. vaginalis genome identified 911 BspA-like entries (TvBspA) sharing TpLRR-like leucine-rich repeats, which represent the largest gene family encoding potential extracellular proteins for the pathogen. A broad range of microorganisms encoding BspA-like proteins was identified and these are mainly known to live on mucosal surfaces, among these T. vaginalis is endowed with the largest gene family. Over 190 TvBspA proteins with inferred transmembrane domains were characterised by a considerable structural diversity between their TpLRR and other types of repetitive sequences and two subfamilies possessed distinct classic sorting signal motifs for endocytosis. One TvBspA subfamily also shared a glycine-rich protein domain with proteins from Clostridium difficile pathogenic strains and C. difficile phages. Consistent with the hypothesis that TvBspA protein structural diversity implies diverse roles, we demonstrated for several TvBspA genes differential expression at the transcript level in different growth conditions. Identified variants of repetitive segments between several TvBspA paralogues and orthologues from two clinical isolates were also consistent with TpLRR and other repetitive sequences to be functionally important. For one TvBspA protein cell surface expression and antibody responses by both female and male T. vaginalis infected patients were also demonstrated.

Conclusions

The biased mucosal habitat for microbial species encoding BspA-like proteins, the characterisation of a vast structural diversity for the TvBspA proteins, differential expression of a subset of TvBspA genes and the cellular localisation and immunological data for one TvBspA; all point to the importance of the TvBspA proteins to various aspects of T. vaginalis pathobiology at the host-pathogen interface.

Abstract

Background

In many eukaryotes, microRNAs (miRNAs) bind to complementary sites in the 3'-untranslated regions (3'-UTRs) of target messenger RNAs (mRNAs) and regulate their expression at the stage of translation. Recent studies have revealed that many miRNAs are evolutionarily conserved; however, the evolution of their target genes has yet to be systematically characterized. We sought to elucidate a set of conserved miRNA/target-gene pairs and to analyse the mechanism underlying miRNA-mediated gene regulation in the early stage of bilaterian evolution.

Results

Initially, we extracted five evolutionarily conserved miRNAs (let-7, miR-1, miR-124, miR-125/lin-4, and miR-34) among five diverse bilaterian animals. Subsequently, we designed a procedure to predict evolutionarily conserved miRNA/target-gene pairs by introducing orthologous gene information. As a result, we extracted 31 orthologous miRNA/target-gene pairs that were conserved among at least four diverse bilaterian animals; the prediction set showed prominent enrichment of orthologous miRNA/target-gene pairs that were verified experimentally. Approximately 84% of the target genes were regulated by three miRNAs (let-7, miR-1, and miR-124) and their function was classified mainly into the following categories: development, muscle formation, cell adhesion, and gene regulation. We used a reporter gene assay to experimentally verify the downregulation of six candidate pairs (out of six tested pairs) in HeLa cells.

Conclusions

The application of our new method enables the identification of 31 miRNA/target-gene pairs that were expected to have been regulated from the era of the common bilaterian ancestor. The downregulation of all six candidate pairs suggests that orthologous information contributed to the elucidation of the primordial set of genes that has been regulated by miRNAs; it was also an efficient tool for the elimination of false positives from the predicted candidates. In conclusion, our study identified potentially important miRNA-target pairs that were evolutionarily conserved throughout diverse bilaterian animals and that may provide new insights into early-stage miRNA functions.

Abstract

Background

The marine cyanobacterium Prochlorococcus marinus, having multiple ecotypes of distinct genotypic/phenotypic traits and being the first documented example of genome shrinkage in free-living organisms, offers an ideal system for studying niche-driven molecular micro-diversity in closely related microbes. The present study, through an extensive comparative analysis of various genomic/proteomic features of 6 high light (HL) and 6 low light (LL) adapted strains, makes an attempt to identify molecular determinants associated with their vertical niche partitioning.

Results

Pronounced strand-specific asymmetry in synonymous codon usage is observed exclusively in LL strains. Distinct dinucleotide abundance profiles are exhibited by 2 LL strains with larger genomes and G+C-content ≈ 50% (group LLa), 4 LL strains having reduced genomes and G+C-content ≈ 35-37% (group LLb), and 6 HL strains. Taking into account the emergence of LLa, LLb and HL strains (based on 16S rRNA phylogeny), a gradual increase in average aromaticity, pI values and beta- & coil-forming propensities and a decrease in mean hydrophobicity, instability indices and helix-forming propensities of core proteins are observed. Greater variations in orthologous gene repertoire are found between LLa and LLb strains, while higher number of positively selected genes exist between LL and HL strains.

Conclusion

Strains of different Prochlorococcus groups are characterized by distinct compositional, physicochemical and structural traits that are not mere remnants of a continuous genetic drift, but are potential outcomes of a grand scheme of niche-oriented stepwise diversification, that might have driven them chronologically towards greater stability/fidelity and invoked upon them a special ability to inhabit diverse oceanic environments.

Abstract

Background

Fungi secrete various proteins that have diverse functions. Prediction of secretory proteins using only one program is unsatisfactory. To enhance prediction accuracy, we constructed Fungal Secretome Database (FSD).

Description

A three-layer hierarchical identification rule based on nine prediction programs was used to identify putative secretory proteins in 158 fungal/oomycete genomes (208,883 proteins, 15.21% of the total proteome). The presence of putative effectors containing known host targeting signals such as RXLX [EDQ] and RXLR was investigated, presenting the degree of bias along with the species. The FSD's user-friendly interface provides summaries of prediction results and diverse web-based analysis functions through Favorite, a personalized repository.

Conclusions

The FSD can serve as an integrated platform supporting researches on secretory proteins in the fungal kingdom. All data and functions described in this study can be accessed on the FSD web site at http://fsd.snu.ac.kr/.

Abstract

Background

Delineating the genetic basis of body composition is important to agriculture and medicine. In addition, the incorporation of gene-gene interactions in the statistical model provides further insight into the genetic factors that underlie body composition traits. We used Bayesian model selection to comprehensively map main, epistatic and sex-specific QTL in an F₂ reciprocal intercross between two chicken lines divergently selected for high or low growth rate.

Results

We identified 17 QTL with main effects across 13 chromosomes and several sex-specific and sex-antagonistic QTL for breast meat yield, thigh + drumstick yield and abdominal fatness. Different sets of QTL were found for both breast muscles [Pectoralis (P) major and P. minor], which suggests that they could be controlled by different regulatory mechanisms. Significant interactions of QTL by sex allowed detection of sex-specific and sex-antagonistic QTL for body composition and abdominal fat. We found several female-specific P. major QTL and sex-antagonistic P. minor and abdominal fatness QTL. Also, several QTL on different chromosomes interact with each other to affect body composition and abdominal fatness.

Conclusions

The detection of main effects, epistasis and sex-dimorphic QTL suggest complex genetic regulation of somatic growth. An understanding of such regulatory mechanisms is key to mapping specific genes that underlie QTL controlling somatic growth in an avian model.

Abstract

The version of this article published in BMC Genomics 2009, 10:558, contains data in Table 1 which are now known to be unreliable, and an illustration, in Figure 1, of unusual miRNA processing events predicted by these unreliable data. In this full-length correction, new data replace those found to be unreliable, leading to a more straightforward interpretation without altering the principle conclusions of the study. Table 1 and associated methods have been corrected, Figure 1 deleted, supplementary file 1 added, and modifications made to the sections "Deep sequencing of small RNAs from grapevine leaf tissue" and "Microarray analysis of miRNA expression". The editors and authors regret the inconvenience caused to readers by premature publication of the original paper.

Background

MicroRNAs are short (~21 base) single stranded RNAs that, in plants, are generally coded by specific genes and cleaved specifically from hairpin precursors. MicroRNAs are critical for the regulation of multiple developmental, stress related and other physiological processes in plants. The recent annotation of the genome of the grapevine (Vitis vinifera L.) allowed the identification of many putative conserved microRNA precursors, grouped into multiple gene families.

Results

Here we use oligonucleotide arrays to provide the first indication that many of these microRNAs show differential expression patterns between tissues and during the maturation of fruit in the grapevine. Furthermore we demonstrate that whole transcriptome sequencing and deep-sequencing of small RNA fractions can be used both to identify which microRNA precursors are expressed in different tissues and to estimate genomic coordinates and patterns of splicing and alternative splicing for many primary miRNA transcripts.

Conclusions

Our results show that many microRNAs are differentially expressed in different tissues and during fruit maturation in the grapevine. Furthermore, the demonstration that whole transcriptome sequencing can be used to identify candidate splicing events and approximate primary microRNA transcript coordinates represents a significant step towards the large-scale elucidation of mechanisms regulating the expression of microRNAs at the transcriptional and post-transcriptional levels.

Abstract

Background

Cyclotides are a family of circular peptides that exhibit a range of biological activities, including anti-bacterial, cytotoxic, anti-HIV activities, and are proposed to function in plant defence. Their high stability has motivated their development as scaffolds for the stabilisation of peptide drugs. Oldenlandia affinis is a member of the Rubiaceae (coffee) family from which 18 cyclotides have been sequenced to date, but the details of their processing from precursor proteins have only begun to be elucidated. To increase the speed at which genes involved in cyclotide biosynthesis and processing are being discovered, an expressed sequence tag (EST) project was initiated to survey the transcript profile of O. affinis and to propose some future directions of research on in vivo protein cyclisation.

Results

Using flow cytometry the holoploid genome size (1C-value) of O. affinis was estimated to be 4,210 - 4,284 Mbp, one of the largest genomes of the Rubiaceae family. High-quality ESTs were identified, 1,117 in total, from leaf cDNAs and assembled into 502 contigs, comprising 202 consensus sequences and 300 singletons. ESTs encoding the cyclotide precursors for kalata B1 (Oak1) and kalata B2 (Oak4) were among the 20 most abundant ESTs. In total, 31 ESTs encoded cyclotide precursors, representing a distinct commitment of 2.8% of the O. affinis transcriptome to cyclotide biosynthesis. The high expression levels of cyclotide precursor transcripts are consistent with the abundance of mature cyclic peptides in O. affinis. A new cyclotide precursor named Oak5 was isolated and represents the first cDNA for the bracelet class of cyclotides in O. affinis. Clones encoding enzymes potentially involved in processing cyclotides were also identified and include enzymes involved in oxidative folding and proteolytic processing.

Conclusion

The EST library generated in this study provides a valuable resource for the study of the cyclisation of plant peptides. Further analysis of the candidates for cyclotide processing discovered in this work will increase our understanding and aid in reconstructing cyclotide production using transgenic systems and will benefit their development in pharmaceutical applications and insect-resistant crop plants.

Abstract

Background

Transposable elements are the most abundant components of all characterized genomes of higher eukaryotes. It has been documented that these elements not only contribute to the shaping and reshaping of their host genomes, but also play significant roles in regulating gene expression, altering gene function, and creating new genes. Thus, complete identification of transposable elements in sequenced genomes and construction of comprehensive transposable element databases are essential for accurate annotation of genes and other genomic components, for investigation of potential functional interaction between transposable elements and genes, and for study of genome evolution. The recent availability of the soybean genome sequence has provided an unprecedented opportunity for discovery, and structural and functional characterization of transposable elements in this economically important legume crop.

Description

Using a combination of structure-based and homology-based approaches, a total of 32,552 retrotransposons (Class I) and 6,029 DNA transposons (Class II) with clear boundaries and insertion sites were structurally annotated and clearly categorized, and a soybean transposable element database, SoyTEdb, was established. These transposable elements have been anchored in and integrated with the soybean physical map and genetic map, and are browsable and visualizable at any scale along the 20 soybean chromosomes, along with predicted genes and other sequence annotations. BLAST search and other infrastracture tools were implemented to facilitate annotation of transposable elements or fragments from soybean and other related legume species. The majority (> 95%) of these elements (particularly a few hundred low-copy-number families) are first described in this study.

Conclusion

SoyTEdb provides resources and information related to transposable elements in the soybean genome, representing the most comprehensive and the largest manually curated transposable element database for any individual plant genome completely sequenced to date. Transposable elements previously identified in legumes, the third largest family of flowering plants, are relatively scarce. Thus this database will facilitate structural, evolutionary, functional, and epigenetic analyses of transposable elements in soybean and other legume species.

Abstract

Background

Physarum polycephalum is a free-living amoebozoan protist displaying a complex life cycle, including alternation between single- and multinucleate stages through sporulation, a simple form of cell differentiation. Sporulation in Physarum can be experimentally induced by several external factors, and Physarum displays many biochemical features typical for metazoan cells, including metazoan-type signaling pathways, which makes this organism a model to study cell cycle, cell differentiation and cellular reprogramming.

Results

In order to identify the genes associated to the light-induced sporulation in Physarum, especially those related to signal transduction, we isolated RNA before and after photoinduction from sporulation- competent cells, and used these RNAs to synthesize cDNAs, which were then analyzed using the 454 sequencing technology. We obtained 16,669 cDNAs that were annotated at every computational level. 13,169 transcripts included hit count data, from which 2,772 displayed significant differential expression (upregulated: 1,623; downregulated: 1,149). Transcripts with valid annotations and significant differential expression were later integrated into putative networks using interaction information from orthologs.

Conclusions

Gene ontology analysis suggested that most significantly downregulated genes are linked to DNA repair, cell division, inhibition of cell migration, and calcium release, while highly upregulated genes were involved in cell death, cell polarization, maintenance of integrity, and differentiation. In addition, cell death- associated transcripts were overrepresented between the upregulated transcripts. These changes are associated to a network of actin-binding proteins encoded by genes that are differentially regulated before and after light induction.

Abstract

Background

Downy mildew is a destructive grapevine disease caused by Plasmopara viticola (Berk. and Curt.) Berl. and de Toni, which can only be controlled by intensive fungicide treatments. Natural sources of resistance from wild grapevine (Vitis) species are used in conventional breeding approaches, but the signals and effectors involved in resistance in this important crop species are not well understood.

Results

Early transcriptional changes associated with P. viticola infection in susceptible V. vinifera and resistant V. riparia plants were analyzed using the Combimatrix microarray platform. Transcript levels were measured 12 and 24 h post-inoculation, reflecting the time points immediately preceding the onset of resistance in V. riparia, as determined by microscopic analysis. Our data indicate that resistance in V. riparia is induced after infection, and is not based on differences in basal gene expression between the two species. The strong and rapid transcriptional reprogramming involves the induction of pathogenesis-related proteins and enzymes required for the synthesis of phenylpropanoid-derived compounds, many of which are also induced, albeit to a lesser extent, in V. vinifera. More interestingly, resistance in V. riparia also involves the specific modulation of numerous transcripts encoding components of signal transduction cascades, hypersensitive reaction markers and genes involved in jasmonate biosynthesis. The limited transcriptional modulation in V. vinifera represents a weak attempted defense response rather than the activation of compatibility-specific pathways.

Conclusions

Several candidate resistance genes were identified that could be exploited in future biotechnological approaches to increase disease resistance in susceptible grapevine species. Measurements of jasmonic acid and methyl jasmonate in infected leaves suggest that this hormone may also be involved in V. riparia resistance to P. viticola.

Abstract

Background

MicroRNAs (miRNAs) are small non-coding RNAs that post-transcriptionally regulate gene expression in a variety of organisms, including insects, vertebrates, and plants. miRNAs play important roles in cell development and differentiation as well as in the cellular response to stress and infection. To date, there are limited reports of miRNA identification in mosquitoes, insects that act as essential vectors for the transmission of many human pathogens, including flaviviruses. West Nile virus (WNV) and dengue virus, members of the Flaviviridae family, are primarily transmitted by Aedes and Culex mosquitoes. Using high-throughput deep sequencing, we examined the miRNA repertoire in Ae. albopictus cells and Cx. quinquefasciatus mosquitoes.

Results

We identified a total of 65 miRNAs in the Ae. albopictus C7/10 cell line and 77 miRNAs in Cx. quinquefasciatus mosquitoes, the majority of which are conserved in other insects such as Drosophila melanogaster and Anopheles gambiae. The most highly expressed miRNA in both mosquito species was miR-184, a miRNA conserved from insects to vertebrates. Several previously reported Anopheles miRNAs, including miR-1890 and miR-1891, were also found in Culex and Aedes, and appear to be restricted to mosquitoes. We identified seven novel miRNAs, arising from nine different precursors, in C7/10 cells and Cx. quinquefasciatus mosquitoes, two of which have predicted orthologs in An. gambiae. Several of these novel miRNAs reside within a ~350 nt long cluster present in both Aedes and Culex. miRNA expression was confirmed by primer extension analysis. To determine whether flavivirus infection affects miRNA expression, we infected female Culex mosquitoes with WNV. Two miRNAs, miR-92 and miR-989, showed significant changes in expression levels following WNV infection.

Conclusions

Aedes and Culex mosquitoes are important flavivirus vectors. Recent advances in both mosquito genomics and high-throughput sequencing technologies enabled us to interrogate the miRNA profile in these two species. Here, we provide evidence for over 60 conserved and seven novel mosquito miRNAs, expanding upon our current understanding of insect miRNAs. Undoubtedly, some of the miRNAs identified will have roles not only in mosquito development, but also in mediating viral infection in the mosquito host.

Abstract

Background

Rice is one of the most important food crops in the world. With increasing world demand for food crops, there is an urgent need to develop new cultivars that have enhanced performance with regard to yield, disease resistance, and so on. Wild rice is expected to provide useful genetic resources that could improve the present cultivated species. However, the quantity and quality of these unexplored resources remain unclear. Recent accumulation of the genomic information of both cultivated and wild rice species allows for their comparison at the molecular level. Here, we compared the genome sequence of Oryza sativa ssp. japonica with sets of bacterial artificial chromosome end sequences (BESs) from two wild rice species, O. rufipogon and O. nivara, and an African rice species, O. glaberrima.

Results

We found that about four to five percent of the BESs of the two wild rice species and about seven percent of the African rice could not be mapped to the japonica genome, suggesting that a substantial number of genes have been lost in the japonica rice lineage; however, their close relatives still possess their counterpart genes. We estimated that during evolution, O. sativa has lost at least one thousand genes that are still preserved in the genomes of the other species. In addition, our BLASTX searches against the non-redundant protein sequence database showed that disease resistance-related proteins were significantly overrepresented in the close relative-specific genomic portions. In total, 235 unmapped BESs of the three relatives matched 83 non-redundant proteins that contained a disease resistance protein domain, most of which corresponded to an NBS-LRR domain.

Conclusion

We found that the O. sativa lineage appears to have recently experienced massive gene losses following divergence from its wild ancestor. Our results imply that the domestication process accelerated large-scale genomic deletions in the lineage of Asian cultivated rice and that the close relatives of cultivated rice have the potential to restore the lost traits.

Abstract

Background

The Lipase Engineering Database (LED) integrates information on sequence, structure and function of lipases, esterases and related proteins with the α/β hydrolase fold. A new superfamily for Candida antarctica lipase A (CALA) was introduced including the recently published crystal structure of CALA. Since CALA has a highly divergent sequence in comparison to other α/β hydrolases, the Lipase Engineering Database was used to classify CALA in the frame of the already established classification system. This involved the comparison of CALA to similar structures as well as sequence-based comparisons against the content of the LED.

Results

The new release 3.0 (December 2009) of the Lipase Engineering Database contains 24783 sequence entries for 18585 proteins as well as 656 experimentally determined protein structures, including the structure of CALA. In comparison to the previous release 1 with 4322 protein and 167 structure entries this update represents a significant increase in data volume. By comparing CALA to representative structures from all superfamilies, a structure from the deacetylase superfamily was found to be most similar to the structure of CALA. While the α/β hydrolase fold is conserved in both proteins, the major difference is found in the cap region. Sequence alignments between both proteins show a sequence similarity of only 15%. A multisequence alignment of both protein families was used to create hidden Markov models for the cap region of CALA and showed that the cap region of CALA is unique among all other proteins of the α/β hydrolase fold. By specifically comparing the substrate binding pocket of CALA to other binding pockets of α/β hydrolases, the binding pocket of Candida rugosa lipase was identified as being highly similar. This similarity also applied to the lid of Candida rugosa lipase in comparison to the potential lid of CALA.

Conclusion

The LED serves as a valuable tool for the systematic analysis of single proteins or protein families. The updated release 3.0 was used for the evaluation of α/β hydrolases. The HTML version of the database with new features is available at http://www.led.uni-stuttgart.de and provides sequences, structures and a set of analysis tools including phylogenetic trees and HMM profiles

Abstract

Background

A high-sensitivity DNA microarray platform requiring nanograms of RNA input facilitates the application of transcriptome analysis to individual skeletal muscle (SM) tissue samples. Culturing myotubes from SM-biopsies enables investigating transcriptional defects and assaying therapeutic strategies. This study compares the transcriptome of aneurally cultured human SM cells versus that of tissue biopsies.

Results

We used the Illumina expression BeadChips to determine the transcriptomic differences between tissue and cultured SM samples from five individuals. Changes in the expression of several genes were confirmed by QuantiGene Plex assay or reverse transcription real-time PCR. In cultured myotubes compared to the tissue, 1216 genes were regulated: 583 down and 633 up. Gene ontology analysis showed that downregulated genes were mainly associated with cytoplasm, particularly mitochondria, and involved in metabolism and the muscle-system/contraction process. Upregulated genes were predominantly related to cytoplasm, endoplasmic reticulum, and extracellular matrix. The most significantly regulated pathway was mitochondrial dysfunction. Apoptosis genes were also modulated. Among the most downregulated genes detected in this study were genes encoding metabolic proteins AMPD1, PYGM, CPT1B and UCP3, muscle-system proteins TMOD4, MYBPC1, MYOZ1 and XIRP2, the proteolytic CAPN3 and the myogenic regulator MYF6. Coordinated reduced expression of five members of the GIMAP gene family, which form a cluster on chromosome 7, was shown, and the GIMAP4-reduction was validated. Within the most upregulated group were genes encoding senescence/apoptosis-related proteins CDKN1A and KIAA1199 and potential regulatory factors HIF1A, TOP2A and CCDC80.

Conclusions

Cultured muscle cells display reductive metabolic and muscle-system transcriptome adaptations as observed in muscle atrophy and they activate tissue-remodeling and senescence/apoptosis processes.

Abstract

Background

Werner Syndrome (WS) is a rare disorder characterized by the premature onset of a number of age-related diseases. The gene responsible for WS encodes a DNA helicase/exonuclease protein believed to affect different aspects of transcription, replication, and/or DNA repair. In addition to genomic instability, human WS cells exhibit oxidative stress. In this report, we have examined the impact of exogenous hydrogen peroxide on the expression profile of mouse embryonic fibroblasts lacking part of the helicase domain of the WRN homologue (here referred to as Wrn^Δhel/Δhel).

Results

Wrn^Δhel/Δhel mutant mouse embryonic fibroblasts exhibit increased oxidative stress. This was reflected by increased intracellular reactive oxygen species (ROS), increased oxidative damage in genomic DNA, changes in ATP/ADP ratios, and a disruption of the inner mitochondrial transmembrane potential when compared to wild type mouse embryonic fibroblasts. Expression profile analyses of hydrogen peroxide-treated wild type cells have indicated significant decreases in the expression of genes involved in mitosis, glycolysis, fatty acid metabolism, nucleic acid metabolism, and cell cycle control, as well as protein modification and stability. Such decreases in these biological processes were not observed in hydrogen peroxide-treated Wrn^Δhel/Δhel cells. Importantly, untreated Wrn^Δhel/Δhel cells already exhibited down regulation of several biological processes decreased in wild type cells that had been treated with hydrogen peroxide.

Conclusion

Expression profiling of Wrn^Δhel/Δhel mutant cells revealed a very different response to exogenous addition of hydrogen peroxide in culture compared to wild type cells. This is due in part to the fact that Wrn^Δhel/Δhel mutant cells already exhibited a modest chronic intracellular oxidative stress.

Abstract

Background

Salix (willow) and Populus (poplar) are members of the Salicaceae family and they share many ecological as well as genetic and genomic characteristics. The interest of using willow for biomass production is growing, which has resulted in increased pressure on breeding of high yielding and resistant clones adapted to different environments. The main purpose of this work was to develop dense genetic linkage maps for mapping of traits related to yield and resistance in willow. We used the Populus trichocarpa genome to extract evenly spaced markers and mapped the orthologous loci in the willow genome. The marker positions in the two genomes were used to study genome evolution since the divergence of the two lineages some 45 mya.

Results

We constructed two linkage maps covering the 19 linkage groups in willow. The most detailed consensus map, S₁, contains 495 markers with a total genetic distance of 2477 cM and an average distance of 5.0 cM between the markers. The S₃ consensus map contains 221 markers and has a total genetic distance of 1793 cM and an average distance of 8.1 cM between the markers. We found high degree of synteny and gene order conservation between willow and poplar. There is however evidence for two major interchromosomal rearrangements involving poplar LG I and XVI and willow LG Ib, suggesting a fission or a fusion in one of the lineages, as well as five intrachromosomal inversions. The number of silent substitutions were three times lower (median: 0.12) between orthologs than between paralogs (median: 0.37 - 0.41).

Conclusions

The relatively slow rates of genomic change between willow and poplar mean that the genomic resources in poplar will be most useful in genomic research in willow, such as identifying genes underlying QTLs of important traits. Our data suggest that the whole-genome duplication occurred long before the divergence of the two genera, events which have until now been regarded as contemporary. Estimated silent substitution rates were 1.28 × 10^-9 and 1.68 × 10^-9 per site and year, which are close to rates found in other perennials but much lower than rates in annuals.

Abstract

Background

The original sequencing and annotation of the Caenorhabditis elegans genome along with recent advances in sequencing technology provide an exceptional opportunity for the genomic analysis of wild-type and mutant strains. Using the Illumina Genome Analyzer, we sequenced the entire genome of Rec-1, a strain that alters the distribution of meiotic crossovers without changing the overall frequency. Rec-1 was derived from ethylmethane sulfonate (EMS)-treated strains, one of which had a high level of transposable element mobility. Sequencing of this strain provides an opportunity to examine the consequences on the genome of altering the distribution of meiotic recombination events.

Results

Using Illumina sequencing and MAQ software, 83% of the base pair sequence reads were aligned to the reference genome available at Wormbase, providing a 21-fold coverage of the genome. Using the software programs MAQ and Slider, we observed 1124 base pair differences between Rec-1 and the reference genome in Wormbase (WS190), and 441 between the mutagenized Rec-1 (BC313) and the wild-type N2 strain (VC2010). The most frequent base-substitution was G:C to A:T, 141 for the entire genome most of which were on chromosomes I or X, 55 and 31 respectively. With this data removed, no obvious pattern in the distribution of the base differences along the chromosomes was apparent. No major chromosomal rearrangements were observed, but additional insertions of transposable elements were detected. There are 11 extra copies of Tc1, and 8 of Tc2 in the Rec-1 genome, most likely the remains of past high-hopper activity in a progenitor strain.

Conclusion

Our analysis of high-throughput sequencing was able to detect regions of direct repeat sequences, deletions, insertions of transposable elements, and base pair differences. A subset of sequence alterations affecting coding regions were confirmed by an independent approach using oligo array comparative genome hybridization. The major phenotype of the Rec-1 strain is an alteration in the preferred position of the meiotic recombination event with no other significant phenotypic consequences. In this study, we observed no evidence of a mutator effect at the nucleotide level attributable to the Rec-1 mutation.

Abstract

Background

Protein phosphorylation is responsible for a large portion of the regulatory functions of eukaryotic cells. Although the list of sequenced genomes of filamentous fungi has grown rapidly, the kinomes of recently sequenced species have not yet been studied in detail. The objective of this study is to apply a comparative analysis of the kinase distribution in different fungal phyla, and to explore its relevance to understanding the evolution of fungi and their taxonomic classification. We have analyzed in detail 12 subgroups of kinases and their distribution over 30 species, as well as their potential use as a classifier for members of the fungal kingdom.

Results

Our findings show that despite the similarity of the kinase distribution in all fungi, their domain distributions and kinome density can potentially be used to classify them and give insight into their evolutionary origin. In general, we found that the overall representation of kinase groups is similar across fungal genomes, the only exception being a large number of tyrosine kinase-like (TKL) kinases predicted in Laccaria bicolor. This unexpected finding underscores the need to continue to sequence fungal genomes, since many species or lineage-specific properties may remain to be discovered. Furthermore, we found that the domain organization significantly varies between the fungal species. Our results suggest that protein kinases and their functional domains strongly reflect fungal taxonomy.

Conclusions

Comparison of the predicted kinomes of sequenced fungi suggests essential signaling functions common to all species, but also specific adaptations of the signal transduction networks to particular species.