BMC Genomics

Unified translation repression mechanism for microRNAs and upstream AUGs.

Publication Date: 2010 Mar 5 PMID: 20205738
Authors: Ajay, S. S. - Athey, B. D. - Lee, I.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: MicroRNAs (miRNAs) are endogenous small RNAs that modulate gene expression at the post-transcriptional level by binding complementary sites in the 3'-UTR. In a recent genome-wide study reporting a new miRNA target class (miBridge), we identified and validated interactions between 5'-UTRs and miRNAs. Separately, upstream AUGs (uAUGs) in 5'-UTRs are known to regulate genes translationally without affecting mRNA levels, one of the mechanisms for miRNA-mediated repression. RESULTS: Using sequence data from whole-genome cDNA alignments we identified 1418 uAUG sequences on the 5'-UTR that specifically interact with 3'-ends of conserved miRNAs. We computationally identified miRNAs that can target six genes through their uAUGs that were previously reported to suppress translation. We extended this meta-analysis by confirming expression of these miRNAs in cell-lines used in the uAUG studies. Similarly, seven members of the KLF family of genes containing uAUGs were computationally identified as interacting with several miRNAs. Using KLF9 as an example (whose protein expression is limited to brain tissue despite the mRNA being expressed ubiquitously), we show computationally that miRNAs expressed only in HeLa cells and not in neuroblastoma (N2A) cells can bind the uAUGs responsible for translation inhibition. Our computed results show that tissue- or cell-line specific repression of protein translation by uAUGs can be explained by the presence or absence of miRNAs that target these uAUG sequences. We propose that these uAUGs represent a subset of miRNA interaction sites on 5'-UTRs in miBridge, whereby a miRNA binding an uAUG hinders the progression of ribosome scanning the mRNA before it reaches the open reading frame (ORF). CONCLUSIONS: While both miRNAs and uAUGs are separately known to down-regulate protein expression, we show that they may be functionally related by identifying potential interactions through a sequence-specific binding mechanism. Using prior experimental evidence that shows uAUG effects on translation repression together with miRNA expression data specific to cell lines, we demonstrate through computational analysis that cell-specific down-regulation of protein expression (while maintaining mRNA levels) correlates well with the simultaneous presence of miRNA and target uAUG sequences in one cell type and not others, suggesting tissue-specific translation repression by miRNAs through uAUGs.

post to: CiteULike

Comprehensive analysis of MHC class I genes from the U-, S-, and Z-lineages in Atlantic salmon.

Publication Date: 2010 Mar 5 PMID: 20205726
Authors: Lukacs, M. F. - Harstad, H. - Bakke, H. G. - Beetz-Sargent, M. - McKinnel, L. - Lubieniecki, K. P. - Koop, B. F. - Grimholt, U.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: We have previously sequenced more than 500kb of the duplicated MHC class I regions in Atlantic salmon. In the IA region we identified the loci for the MHC class I gene Sasa-UBA in addition to a soluble MHC class I molecule, Sasa-ULA. A pseudolocus for Sasa-UCA was identified in the nonclassical IB region. Both regions contained genes for antigen presentation, as wells as orthologues to other genes residing in the human MHC region. RESULTS: The genomic localisation of two MHC class I lineages (Z and S) has been resolved. 7 BACs were sequenced using a combination of standard Sanger and 454 sequencing. The new sequence data extended the IA region with 150kb identifying the location of one Z-lineage locus, ZAA. The IB region was extended with 350kb including three new Z-lineage loci, ZBA, ZCA and ZDA in addition to a UGA locus. An allelic version of the IB region contained a functional UDA locus in addition to the UCA pseudolocus. Additionally a BAC harbouring two MHC class I genes (UHA) was placed on linkage group 14, while a BAC containing the S-lineage locus SAA (previously known as UAA) was placed on LG10. Gene expression studies showed limited expression range for all class I genes with exception of UBA being dominantly expressed in gut, spleen and gills, and ZAA with high expression in blood. CONCLUSION: Here we describe the genomic organization of MHC class I loci from the U-, Z-, and S-lineages in Atlantic salmon. Nine of the described class I genes are located in the extension of the duplicated IA and IB regions, while three class I genes are found on two separate linkage groups. The gene organization of the two regions indicates that the IB region is evolving at a different pace than the IA region. Expression profiling, polymorphic content, peptide binding properties and phylogenetic relationship show that Atlantic salmon has only one MHC class Ia gene (UBA), in addition to a multitude of nonclassical MHC class I genes from the U-, S- and Z-lineages.

post to: CiteULike

Chemogenomic and transcriptome analysis identifies mode of action of the chemosensitizing agent CTBT (7-chlorotetrazolo[5,1-c]benzo[1,2,4]triazine).

Publication Date: 2010 Mar 4 PMID: 20202201
Authors: Batova, M. - Klobucnikova, V. - Oblasova, Z. - Gregan, J. - Zahradnik, P. - Hapala, I. - Subik, J. - Schuller, C.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: CTBT (7-chlorotetrazolo[5,1-c]benzo[1,2,4]triazine) increases efficacy of commonly used antifungal agents by an unknown mechanism. It increases the susceptibility of Saccharomyces cerevisiae, Candida albicans and Candida glabrata cells to cycloheximide, 5-fluorocytosine and azole antimycotic drugs. Here we elucidate CTBT mode of action with a combination of systematic genetic and transcriptome analysis. RESULTS: To identify the cellular processes affected by CTBT, we screened the systematic haploid deletion mutant collection for CTBT sensitive mutants. We identified 169 hypersensitive deletion mutants. The deleted genes encode proteins mainly involved in mitochondrial functions, DNA repair, transcription and chromatin remodeling, and oxidative stress response. We found that the susceptibility of yeast cells to CTBT depends on molecular oxygen. Transcriptome analysis of the immediate early response to CTBT revealed rapid induction of oxidant and stress response defense genes. Many of these genes depend on the transcription factors Yap1 and Cin5. Yap1 accumulates rapidly in the nucleus in CTBT treated cells suggesting acute oxidative stress. Moreover, molecular calculations supported a superoxide generating activity of CTBT. Superoxide production in vivo by CTBT was found associated to mitochondria as indicated by oxidation of MitoSOX Red. CONCLUSION: We conclude that CTBT causes intracellular superoxide production and oxidative stress in fungal cells and is thus enhancing antimycotic drug effects drug effects by a secondary stress.

post to: CiteULike

Genome dynamics of Bartonella grahamii in micro-populations of woodland rodents.

Publication Date: 2010 Mar 4 PMID: 20202191
Authors: Berglund, E. C. - Ehrenborg, C. - Vinnere Pettersson, O. - Granberg, F. - Naslund, K. - Holmberg, M. - Andersson, S. G.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: Rodents represent a high-risk reservoir for the emergence of new human pathogens. The recent completion of the 2.3 Mb genome of Bartonella grahamii, one of the most prevalent blood-borne bacteria in wild rodents, revealed a higher abundance of genes for host-cell interaction systems than in the genomes of closely related human pathogens. The sequence variability within the global B. grahamii population was recently investigated by multi locus sequence typing, but no study on the variability of putative host-cell interaction systems has been performed. RESULTS: To study the population dynamics of B. grahamii, we analyzed the genomic diversity on a whole-genome scale of 27 B. grahamii strains isolated from four different species of wild rodents in three geographic locations separated by less than 30 km. Even using highly variable spacer regions, only 3 sequence types were identified. This low sequence diversity contrasted with a high variability in genome content. Microarray comparative genome hybridizations identified genes for outer surface proteins, including a repeated region containing the fha gene for filamentous hemaggluttinin and a plasmid that encodes a type IV secretion system, as the most variable. The estimated generation times in liquid culture medium for a subset of strains ranged from 5 to 22 hours, but did not correlate with sequence type or presence/absence patterns of the fha gene or the plasmid. CONCLUSION: Our study has revealed a geographic microstructure of B. grahamii in wild rodents. Despite near-identity in nucleotide sequence, major differences were observed in gene presence/absence patterns that did not segregate with host species. This suggests that genetically similar strains can infect a range of different hosts.

post to: CiteULike

Promiscuity of enhancer, coding and non-coding transcription functions in ultraconserved elements.

Publication Date: 2010 Mar 4 PMID: 20202189
Authors: Licastro, D. - Gennarino, V. A. - Petrera, F. - Sanges, R. - Banfi, S. - Stupka, E.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: Ultraconserved elements (UCEs) are highly constrained elements of mammalian genomes, whose functional role has not been completely elucidated yet. Previous studies have shown that some of them act as enhancers in mouse, while some others are expressed in both normal and cancer-derived human tissues. Only one UCE element so far was shown to present these two functions concomitantly, as had been observed in other isolated instances of single, non ultraconserved enhancer elements. RESULTS: We used a custom microarray to assess the levels of UCE transcription during mouse development and integrated these data with published microarray and next-generation sequencing datasets as well as with newly produced PCR validation experiments. We show that a large fraction of non-exonic UCEs is transcribed across all developmental stages examined from only one DNA strand. Although the nature of these transcripts remains a mistery, our meta-analysis of RNA-Seq datasets indicates that they are unlikely to be short RNAs and that some of them might encode nuclear transcripts. In the majority of cases this function overlaps with the already established enhancer function of these elements during mouse development. Utilizing several next-generation sequencing datasets, we were further able to show that the level of expression observed in non-exonic UCEs is significantly higher than in random regions of the genome and that this is also seen in other regions which act as enhancers. CONCLUSION: Our data shows that the concurrent presence of enhancer and transcript function in non-exonic UCE elements is more widespread than previously shown. Moreover through our own experiments as well as the use of next-generation sequencing datasets, we were able to show that the RNAs encoded by non-exonic UCEs are likely to be long RNAs transcribed from only one DNA strand.

post to: CiteULike

Genome-wide transcriptome analysis of the transition from primary to secondary stem development in Populus trichocarpa.

Publication Date: 2010 Mar 4 PMID: 20199690
Authors: Dharmawardhana, P. - Brunner, A. M. - Strauss, S. H.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: With its genome sequence and other experimental attributes, Populus trichocarpa has become the model species for genomic studies of wood development. Wood is derived from secondary growth of tree stems, and begins with the development of a ring of vascular cambium in the young developing stem. The terminal region of the developing shoot provides a steep developmental gradient from primary to secondary growth that facilitates identification of genes that play specialized functions during each of these phases of growth. RESULTS: Using a genomic microarray representing the majority of the transcriptome, we profiled gene expression in stem segments that spanned primary to secondary growth. We found 3,016 genes that were differentially expressed during stem development (Q-value < 0.05; >2-fold expression variation), and 15% of these genes encode proteins with no significant identities to known genes. We identified all gene family members putatively involved in secondary growth for carbohydrate active enzymes, tubulins, actins, actin depolymerizing factors, fasciclin-like AGPs, and vascular development-associated transcription factors. Almost 70% of expressed transcription factors were upregulated during the transition to secondary growth. The primary shoot elongation region of the stem contained specific carbohydrate active enzyme and expansin family members that are likely to function in primary cell wall synthesis and modification. Genes involved in plant defense and protective functions were also dominant in the primary growth region. CONCLUSION: Our results describe the global patterns of gene expression that occur during the transition from primary to secondary stem growth. We were able to identify three major patterns of gene expression and over-represented gene ontology categories during stem development. The new regulatory factors and cell wall biogenesis genes that we identified provide candidate genes for further functional characterization, as well as new tools for molecular breeding and biotechnology aimed at improvement of tree growth rate, crown form, and wood quality.

post to: CiteULike

Transcribed-ultra conserved region expression profiling from low-input total RNA.

Publication Date: 2010 Mar 3 PMID: 20199688
Authors: Scaruffi, P. - Stigliani, S. - Coco, S. - Valdora, F. - De Vecchi, C. - Bonassi, S. - Tonini, G. P.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: Ultra Conserved Regions (UCRs) are a class of 481 noncoding sequences located in both intra- and inter-genic regions of the genome. The recent findings that they are significantly altered in adult chronic lymphocytic leukemias, carcinomas, and pediatric neuroblastomas lead to the hypothesis that UCRs may play a role in tumorigenesis. RESULTS: We present a novel application of Ribo-SPIA isothermal linear amplification of minute RNA quantities for quantifying Transcribed-UCR (T-UCR) expression by quantitative PCR. Direct comparison of non-amplified with amplified cDNA in two neuroblastoma cell lines showed that the amplification approach increases sensitivity and repeatability in T-UCR quantification. It is noteworthy that the Ribo-SPIA step allowed us to analyze all 481 T-UCRs by using 150 ng of RNA, while introducing a minimal bias and preserving the magnitude of relative expression. Only the less abundant T-UCRs have high intra-assay variability, consistently with the Poisson distribution statistics and stochastic effects on PCR repeatability. CONCLUSIONS: We demonstrated that the quantification procedure shown here is an accurate and reliable technique for genome-wide non coding gene (i.e., UCRs) profiling using small amounts of RNA. This issue is particularly important because studies of transcription regulation are increasingly conducted in small homogeneous samples, such as laser capture microdissected or sorted cell populations.

post to: CiteULike

MicroRNAs of Bombyx mori identified by Solexa sequencing.

Publication Date: 2010 Mar 3 PMID: 20199675
Authors: Liu, S. - Li, D. - Li, Q. - Zhao, P. - Xiang, Z. - Xia, Q.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: MicroRNA (miRNA) and other small regulatory RNAs contribute to the modulation of a large number of cellular processes. We sequenced three small RNA libraries prepared from the whole body, and the anterior-middle and posterior silk glands of Bombyx mori, with a view to expanding the repertoire of silkworm miRNAs and exploring transcriptional differences in miRNAs between segments of the silk gland. RESULTS: With the aid of large-scale Solexa sequencing technology, we validated 257 unique miRNA genes, including 202 novel and 55 previously reported genes, corresponding to 324 loci in the silkworm genome. Over 30 known silkworm miRNAs were further corrected in their sequence constitutes and length. A number of reads originated from the loop regions of the precursors of two previously reported miRNAs (bmo-miR-1920 and miR-1921). Interestingly, the majority of the newly identified miRNAs were silkworm-specific, 23 unique miRNAs were widely conserved from invertebrates to vertebrates, 13 unique miRNAs were limited to invertebrates, and 32 were confined to insects. We identified 24 closely positioned clusters and 45 paralogs of miRNAs in the silkworm genome. However, sequence tags showed that paralogs or clusters were not prerequisites for coordinated transcription and accumulation. The majority of silkworm-specific miRNAs were located in transposable elements, and displayed significant differences in abundance between the anterior-middle and posterior silk gland. CONCLUSIONS: Conservative analysis revealed that miRNAs can serve as phylogenetic markers and function in evolutionary signaling. The newly identified miRNAs greatly enrich the repertoire of insect miRNAs, and provide insights into miRNA evolution, biogenesis, and expression in insects. The differential expression of miRNAs in the anterior-middle and posterior silk glands supports their involvement as new levels in the regulation of the silkworm silk gland.

post to: CiteULike

Conservation and divergence of known apicomplexan transcriptional regulons.

Publication Date: 2010 Mar 3 PMID: 20199665
Authors: Essien, K. - Stoeckert, C. J. Jr
Journal: BMC Genomics

ABSTRACT: BACKGROUND: The apicomplexans are a diverse phylum of parasites causing an assortment of diseases including malaria in a wide variety of animals and lymphoproliferation in cattle. Little is known about how these varied parasites regulate their transcriptional regulons. Even less is known about how regulon systems, consisting of transcription factors and target genes together with their associated biological process, evolve in these diverse parasites. RESULTS: In order to obtain insights into the differences in transcriptional regulation between these parasites we compared the orthology profiles of putative malaria transcription factors across species and examined the enrichment patterns of four binding sites across eleven apicomplexans. About three-fifths of the factors are broadly conserved in several phylogenetic orders of sequenced apicomplexans. This observation suggests the existence of regulons whose regulation is conserved across this ancient phylum. Transcription factors not broadly conserved across the phylum are possibly involved in regulon systems that have diverged between species. Examining binding site enrichment patterns in light of transcription factor conservation patterns suggests a second mode via which regulon systems may diverge - rewiring of existing transcription factors and their associated binding sites in specific ways. Integrating binding sites with transcription factor conservation patterns also facilitated prediction of putative regulators for one of the binding sites. CONCLUSIONS: Even though transcription factors are underrepresented in apicomplexans, the distribution of these factors and their associated regulons reflect common and family-specific transcriptional regulatory processes.

post to: CiteULike

Distribution of candidate genes for experimentally induced arthritis in rats.

Publication Date: 2010 Mar 2 PMID: 20196835
Authors: Andersson, L. - Stahl, F.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: Rat models are frequently used to link genomic regions to experimentally induced arthritis in quantitative trait locus (QTL) analyses. To facilitate the search for candidate genes within such regions, we have previously developed an application (CGC) that uses weighted keywords to rank genes based on their descriptive text. In this study, CGC is used for analyzing the localization of candidate genes from two viewpoints: distribution over the rat genome and functional connections between arthritis QTLs. METHODS: To investigate if candidate genes identified by CGC are more likely to be found inside QTLs, we ranked 2403 genes genome wide in rat. The number of genes within different ranges of CGC scores localized inside and outside QTLs was then calculated. Furthermore, we investigated if candidate genes within certain QTLs share similar functions, and if these functions could be connected to genes within other QTLs. Based on references between genes in OMIM, we created connections between genes in QTLs identified in two distinct rat crosses. In this way, QTL pairs with one QTL from each cross that share an unexpectedly high number of gene connections were identified. The genes that were found to connect a pair of QTLs were then functionally analysed using a publicly available classification tool. RESULTS: Out of the 2403 genes ranked by the CGC application, 1160 were localized within QTL regions. No difference was observed between highly and lowly rated genes. Hence, highly rated candidate genes for arthritis seem to be distributed randomly inside and outside QTLs. Furthermore, we found five pairs of QTLs that shared a significantly high number of interconnected genes. When functionally analyzed, most genes connecting two QTLs could be included in a single functional cluster. Thus, the functional connections between these genes could very well be involved in the development of an arthritis phenotype. CONCLUSIONS: From the genome wide CGC search, we conclude that candidate genes for arthritis in rat are randomly distributed between QTL and non-QTL regions. We do however find certain pairs of QTLs that share a large number of functionally connected candidate genes, suggesting that these QTLs contain a number of genes involved in similar functions contributing to the arthritis phenotype.

post to: CiteULike

Genome-wide transcription factor binding site/promoter databases for the analysis of gene sets and co-occurrence of transcription factor binding motifs.

Publication Date: 2010 Mar 1 PMID: 20193056
Authors: Veerla, S. - Ringner, M. - Hoglund, M.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: The use of global gene expression profiling is a well established approach to understand biological processes. One of the major goals of these investigations is to identify sets of genes with similar expression patterns. Such gene signatures may be very informative and reveal new aspects of particular biological processes. A logical and systematic next step is to reduce the identified gene signatures to the regulatory components that induce the relevant gene expression changes. A central issue in this context is to identify transcription factors, or transcription factor binding sites (TFBS), likely to be of importance for the expression of the gene signatures. RESULTS: We develop a strategy that efficiently produces TFBS/promoter databases based on user-defined criteria. The resulting databases constitute all genes in the Santa Cruz database and the positions for all TFBS provided by the user as PWMs. These databases are then used for two purposes, to identify significant TFBS in the promoters in sets of genes and to identify clusters of co-occurring TFBS. We use two criteria for significance, significantly enriched TFBS in terms of total number of binding sites for the promoters, and significantly present TFBS in terms of the fraction of promoters with binding sites. Significant TFBS are identified by a re-sampling procedure in which the query gene set is compared with typically 100 000 gene lists of similar size randomly drawn from the TFBS/promoter database. We apply this strategy to a large number of published ChIP-Chip data sets and show that the proposed approach faithfully reproduces ChIP-Chip results. The strategy also identifies relevant TFBS when analyzing gene signatures obtained from the MSigDB database. In addition, we show that several TFBS are highly correlated and that co-occurring TFBS define functionally related sets of genes. CONCLUSIONS: The presented approach of promoter analysis faithfully reproduces the results from several ChIP-Chip and MSigDB derived gene sets and hence may prove to be an important method in the analysis of gene signatures obtained through ChIP-Chip or global gene expression experiments. We show that TFBS are organized in clusters of co-occurring TFBS that together define highly coherent sets of genes.

post to: CiteULike

Complex sense-antisense architecture of TNFAIP1/POLDIP2 on 17q11.2 represents a novel transcriptional structural-functional gene module involved in breast cancer progression.

Publication Date: 2010 PMID: 20158880
Authors: Grinchuk, O. V. - Motakis, E. - Kuznetsov, V. A.
Journal: BMC Genomics

BACKGROUND: A sense-antisense gene pair (SAGP) is a gene pair where two oppositely transcribed genes share a common nucleotide sequence region. In eukaryotic genomes, SAGPs can be organized in complex sense-antisense architectures (CSAGAs) in which at least one sense gene shares loci with two or more antisense partners. As shown in several case studies, SAGPs may be involved in cancers, neurological diseases and complex syndromes. However, CSAGAs have not yet been characterized in the context of human disease or cancer. RESULTS: We characterize five genes (TMEM97, IFT20, TNFAIP1, POLDIP2 and TMEM199) organized in a CSAGA on 17q11.2 (we term this the TNFAIP1/POLDIP2 CSAGA) and demonstrate their strong and reproducible co-regulatory transcription pattern in breast cancer tumours. Genes of the TNFAIP1/POLDIP2 CSAGA are located inside the smallest region of recurrent amplification on 17q11.2 and their expression profile correlates with the DNA copy number of the region. Survival analysis of a group of 410 breast cancer patients revealed significant survival-associated individual genes and gene pairs in the TNFAIP1/POLDIP2 CSAGA. Moreover, several of the gene pairs associated with survival, demonstrated synergistic effects. Expression of genes-members of the TNFAIP1/POLDIP2 CSAGA also strongly correlated with expression of genes of ERBB2 core region of recurrent amplification on 17q12. We clearly demonstrate that the observed co-regulatory transcription profile of the TNFAIP1/POLDIP2 CSAGA is maintained not only by a DNA amplification mechanism, but also by chromatin remodelling and local transcription activation. CONCLUSION: We have identified a novel TNFAIP1/POLDIP2 CSAGA and characterized its co-regulatory transcription profile in cancerous breast tissues. We suggest that the TNFAIP1/POLDIP2 CSAGA represents a clinically significant transcriptional structural-functional gene module associated with amplification of the genomic region on 17q11.2 and correlated with expression ERBB2 amplicon core genes in breast cancer. Co-expression pattern of this module correlates with histological grades and a poor prognosis in breast cancer when over-expressed. TNFAIP1/POLDIP2 CSAGA maps the risks of breast cancer relapse onto the complex genomic locus on 17q11.2.

post to: CiteULike

Bimodal gene expression patterns in breast cancer.

Publication Date: 2010 PMID: 20158879
Authors: Bessarabova, M. - Kirillov, E. - Shi, W. - Bugrim, A. - Nikolsky, Y. - Nikolskaya, T.
Journal: BMC Genomics

We identified a set of genes with an unexpected bimodal distribution among breast cancer patients in multiple studies. The property of bimodality seems to be common, as these genes were found on multiple microarray platforms and in studies with different end-points and patient cohorts. Bimodal genes tend to cluster into small groups of four to six genes with synchronised expression within the group (but not between the groups), which makes them good candidates for robust conditional descriptors. The groups tend to form concise network modules underlying their function in cancerogenesis of breast neoplasms.

post to: CiteULike

Data integration and exploration for the identification of molecular mechanisms in tumor-immune cells interaction.

Publication Date: 2010 PMID: 20158878
Authors: Mlecnik, B. - Sanchez-Cabo, F. - Charoentong, P. - Bindea, G. - Pages, F. - Berger, A. - Galon, J. - Trajanoski, Z.
Journal: BMC Genomics

Cancer progression is a complex process involving host-tumor interactions by multiple molecular and cellular factors of the tumor microenvironment. Tumor cells that challenge immune activity may be vulnerable to immune destruction. To address this question we have directed major efforts towards data integration and developed and installed a database for cancer immunology with more than 1700 patients and associated clinical data and biomolecular data. Mining of the database revealed novel insights into the molecular mechanisms of tumor-immune cell interaction. In this paper we present the computational tools used to analyze integrated clinical and biomolecular data. Specifically, we describe a database for heterogeneous data types, the interfacing bioinformatics and statistical tools including clustering methods, survival analysis, as well as visualization methods. Additionally, we discuss generic issues relevant to the integration of clinical and biomolecular data, as well as recent developments in integrative data analyses including biomolecular network reconstruction and mathematical modeling.

post to: CiteULike

Analysis of deep sequencing microRNA expression profile from human embryonic stem cells derived mesenchymal stem cells reveals possible role of let-7 microRNA family in downstream targeting of hepatic nuclear factor 4 alpha.

Publication Date: 2010 PMID: 20158877
Authors: Koh, W. - Sheng, C. T. - Tan, B. - Lee, Q. Y. - Kuznetsov, V. - Kiang, L. S. - Tanavde, V.
Journal: BMC Genomics

BACKGROUND: Recent literature has revealed that genetic exchange of microRNA between cells can be essential for cell-cell communication, tissue-specificity and developmental processes. In stem cells, as in other cells, this can be accomplished through microvesicles or exosome mediated transfer. However, molecular profiles and functions of microRNAs within the cells and in their exosomes are poorly studied. Next generation sequencing technologies could provide a broad-spectrum of microRNAs and their expression and identify possible microRNA targets. In this work, we performed deep sequencing of microRNAs to understand the profile and expression of the microRNAs in microvesicles and intracellular environment of human embryonic stem cells derived mesenchymal stem cells (hES-MSC). We outline a workflow pertaining to visualizing, statistical analysis and interpreting deep sequencing data of known intracellular and extracellular microRNAs from hES-MSC). We utilized these results of which directed our attention towards establishing hepatic nuclear factor 4 alpha (HNF4A) as a downstream target of let-7 family of microRNAs. RESULTS: In our study, significant differences in expression profile of microRNAs were found in the intracellular and extracellular environment of hES-MSC. However, a high level of let-7 family of microRNAs is predominant in both intra- and extra- cellular samples of hES-MSC. Further results derived from visualization of our alignment data and network analysis showed that let-7 family microRNAs could affect the downstream target HNF4A, which is a known endodermal differentiation marker. The elevated presence of let-7 microRNA in both intracellular and extra cellular environment further suggests a possible intercellular signalling mechanism through microvesicles transfer. We suggest that let-7 family microRNAs might play a signalling role via such a mechanism amongst populations of stem cells in maintaining self renewal property by suppressing HNF4A expression. This is in line with recent paradigm where microRNAs regulate self-renewal and differentiation pathways of embryonic stem cells by forming an integral biological network with transcription factors. CONCLUSION: In summary, our study using a combination of alignment, statistical and network analysis tools to examine deep sequencing data of microRNAs in hES-MSC has led to a result that (i) identifies intracellular and exosome microRNA expression profiles of hES-MSC with a possible mechanism of miRNA mediated intercellular regulation by these cells and (ii) placed HNF4A within the cross roads of regulation by the let-7 family of microRNAs.

post to: CiteULike

Differences in the transactivation domains of p53 family members: a computational study.

Publication Date: 2010 PMID: 20158876
Authors: Mavinahalli, J. N. - Madhumalar, A. - Beuerman, R. W. - Lane, D. P. - Verma, C.
Journal: BMC Genomics

The N terminal transactivation domain of p53 is regulated by ligases and coactivator proteins. The functional conformation of this region appears to be an alpha helix which is necessary for its appropriate interactions with several proteins including MDM2 and p300. Folding simulation studies have been carried out to examine the propensity and stability of this region and are used to understand the differences between the family members with the ease of helix formation following the order p53 > p73 > p63. It is clear that hydrophobic clusters control the kinetics of helix formation, while electrostatic interactions control the thermodynamic stability of the helix. Differences in these interactions between the family members may partially account for the differential binding to, and regulation by, MDM2 (and MDMX). Phosphorylations of the peptides further modulate the stability of the helix and control associations with partner proteins.

post to: CiteULike

Projection of gene-protein networks to the functional space of the proteome and its application to analysis of organism complexity.

Publication Date: 2010 PMID: 20158875
Authors: Kanapin, A. A. - Mulder, N. - Kuznetsov, V. A.
Journal: BMC Genomics

We consider the problem of biological complexity via a projection of protein-coding genes of complex organisms onto the functional space of the proteome. The latter can be defined as a set of all functions committed by proteins of an organism. Alternative splicing (AS) allows an organism to generate diverse mature RNA transcripts from a single mRNA strand and thus it could be one of the key mechanisms of increasing of functional complexity of the organism's proteome and a driving force of biological evolution. Thus, the projection of transcription units (TU) and alternative splice-variant (SV) forms onto proteome functional space could generate new types of relational networks (e.g. SV-protein function networks, SFN) and lead to discoveries of novel evolutionarily conservative functional modules. Such types of networks might provide new reliable characteristics of organism complexity and a better understanding of the evolutionary integration and plasticity of interconnection of genome-transcriptome-proteome functions. RESULTS: We use the InterPro and UniProt databases to attribute descriptive features (keywords) to protein sequences. UniProt database includes a controlled and curated vocabulary of specific descriptors or keywords. The keywords have been assigned to a protein sequence via conserved domains or via similarity with annotated sequences. Then we consider the unique combinations of keywords as the protein functional labels (FL), which characterize the biological functions of the given protein and construct the contingency tables and graphs providing the projections of transcription units (TU) and alternative splice-variants (SV) onto all FL of the proteome of a given organism. We constructed SFNs for organisms with different evolutionary history and levels of complexity, and performed detailed statistical parameterization of the networks. CONCLUSIONS: The application of the algorithm to organisms with different evolutionary history and level of biological complexity (nematode, fruit fly, vertebrata) reveals that the parameters describing SFN correlate with the complexity of a given organism. Using statistical analysis of the links of the functional networks, we propose new features of evolution of protein function acquisition. We reveal a group of genes and corresponding functions, which could be attributed to an early conservative part of the cellular machinery essential for cell viability and survival. We identify and provide characteristics of functional switches in the polyform group of TUs in different organisms. Based on comparison of mouse and human SFNs, a role of alternative splicing as a necessary source of evolution towards more complex organisms is demonstrated. The entire set of FL across many organisms could be used as a draft of the catalogue of the functional space of the proteome world.

post to: CiteULike

Computational approaches for detecting protein complexes from protein interaction networks: a survey.

Publication Date: 2010 PMID: 20158874
Authors: Li, X. - Wu, M. - Kwoh, C. K. - Ng, S. K.
Journal: BMC Genomics

BACKGROUND: Most proteins form macromolecular complexes to perform their biological functions. However, experimentally determined protein complex data, especially of those involving more than two protein partners, are relatively limited in the current state-of-the-art high-throughput experimental techniques. Nevertheless, many techniques (such as yeast-two-hybrid) have enabled systematic screening of pairwise protein-protein interactions en masse. Thus computational approaches for detecting protein complexes from protein interaction data are useful complements to the limited experimental methods. They can be used together with the experimental methods for mapping the interactions of proteins to understand how different proteins are organized into higher-level substructures to perform various cellular functions. RESULTS: Given the abundance of pairwise protein interaction data from high-throughput genome-wide experimental screenings, a protein interaction network can be constructed from protein interaction data by considering individual proteins as the nodes, and the existence of a physical interaction between a pair of proteins as a link. This binary protein interaction graph can then be used for detecting protein complexes using graph clustering techniques. In this paper, we review and evaluate the state-of-the-art techniques for computational detection of protein complexes, and discuss some promising research directions in this field. CONCLUSIONS: Experimental results with yeast protein interaction data show that the interaction subgraphs discovered by various computational methods matched well with actual protein complexes. In addition, the computational approaches have also improved in performance over the years. Further improvements could be achieved if the quality of the underlying protein interaction data can be considered adequately to minimize the undesirable effects from the irrelevant and noisy sources, and the various biological evidences can be better incorporated into the detection process to maximize the exploitation of the increasing wealth of biological knowledge available.

post to: CiteULike

Collective motions and specific effectors: a statistical mechanics perspective on biological regulation.

Publication Date: 2010 PMID: 20158873
Authors: Giuliani, A.
Journal: BMC Genomics

BACKGROUND: The interaction of a multiplicity of scales in both time and space is a fundamental feature of biological systems. The complementation of macroscopic (entire organism) and microscopic (molecular biology) views with a mesoscopic level of analysis able to connect the different planes of investigation is urgently needed. This will allow to both obtain a general frame of reference for rationalizing the burden of data coming from high throughput technologies and to derive effective operational views on biological systems. RESULTS: The network paradigm in which microscopic level elements (nodes) are each other related by functional links so giving rise to both global (entire network) and local (specific) behavior is a promising metaphor to try and develop a statistical mechanics inspired approach for biological systems. Here we show the application of this paradigm to different systems going from yeast metabolism to murine macrophages response to immune stimulation. CONCLUSIONS: The need to complement the purely molecular view with mesoscopic approaches is evident in all the studied examples that in turn demonstrate the untenability of the simple ergodic approach dominant in molecular biology in which the data coming from huge ensemble of cells are considered as relative to a single 'average' cell.

post to: CiteULike

Parameterization of disorder predictors for large-scale applications requiring high specificity by using an extended benchmark dataset.

Publication Date: 2010 PMID: 20158872
Authors: Sirota, F. L. - Ooi, H. S. - Gattermayer, T. - Schneider, G. - Eisenhaber, F. - Maurer-Stroh, S.
Journal: BMC Genomics

BACKGROUND: Algorithms designed to predict protein disorder play an important role in structural and functional genomics, as disordered regions have been reported to participate in important cellular processes. Consequently, several methods with different underlying principles for disorder prediction have been independently developed by various groups. For assessing their usability in automated workflows, we are interested in identifying parameter settings and threshold selections, under which the performance of these predictors becomes directly comparable. RESULTS: First, we derived a new benchmark set that accounts for different flavours of disorder complemented with a similar amount of order annotation derived for the same protein set. We show that, using the recommended default parameters, the programs tested are producing a wide range of predictions at different levels of specificity and sensitivity. We identify settings, in which the different predictors have the same false positive rate. We assess conditions when sets of predictors can be run together to derive consensus or complementary predictions. This is useful in the framework of proteome-wide applications where high specificity is required such as in our in-house sequence analysis pipeline and the ANNIE webserver. CONCLUSIONS: This work identifies parameter settings and thresholds for a selection of disorder predictors to produce comparable results at a desired level of specificity over a newly derived benchmark dataset that accounts equally for ordered and disordered regions of different lengths.

post to: CiteULike

Bioinformatic search of plant microtubule-and cell cycle related serine-threonine protein kinases.

Publication Date: 2010 PMID: 20158871
Authors: Karpov, P. A. - Nadezhdina, E. S. - Yemets, A. I. - Matusov, V. G. - Nyporko, A. Y. - Shashina, N. Y. - Blume, Y. B.
Journal: BMC Genomics

A bioinformatic search was carried for plant homologues of human serine-threonine protein kinases involved in regulation of cell division and microtubule protein phosphorylation (SLK, PAK6, PAK7, MARK1, MAST2, TTBK1, TTBK2, AURKA, PLK1, PLK4 and PASK). A number of SLK, MAST2 and AURKA plant homologues were identified. The closest identified homologue of human AURKA kinase was a protein of unknown function, A7PY12/GSVIVT00026259001 from Vitis vinifera (herein named as "STALK", Serine-Threonine Aurora-Like Kinase). Analysis of STALK's three-dimensional structure confirmed its relationship to the subgroup of AURKA-like protein kinases.

post to: CiteULike

Reducing the haystack to find the needle: improved protein identification after fast elimination of non-interpretable peptide MS/MS spectra and noise reduction.

Publication Date: 2010 PMID: 20158870
Authors: Mujezinovic, N. - Schneider, G. - Wildpaner, M. - Mechtler, K. - Eisenhaber, F.
Journal: BMC Genomics

BACKGROUND: Tandem mass spectrometry (MS/MS) has become a standard method for identification of proteins extracted from biological samples but the huge number and the noise contamination of MS/MS spectra obstruct swift and reliable computer-aided interpretation. Typically, a minor fraction of the spectra per sample (most often, only a few %) and about 10% of the peaks per spectrum contribute to the final result if protein identification is not prevented by the noise at all. RESULTS: Two fast preprocessing screens can substantially reduce the haystack of MS/MS data. (1) Simple sequence ladder rules remove spectra non-interpretable in peptide sequences. (2) Modified Fourier-transform-based criteria clear background in the remaining data. In average, only a remainder of 35% of the MS/MS spectra (each reduced in size by about one quarter) has to be handed over to the interpretation software for reliable protein identification essentially without loss of information, with a trend to improved sequence coverage and with proportional decrease of computer resource consumption. CONCLUSIONS: The search for sequence ladders in tandem MS/MS spectra with subsequent noise suppression is a promising strategy to reduce the number of MS/MS spectra from electro-spray instruments and to enhance the reliability of protein matches. Supplementary material and the software are available from an accompanying WWW-site with the URL http://mendel.bii.a-star.edu.sg/mass-spectrometry/MSCleaner-2.0/.

post to: CiteULike

Statistics of protein-DNA binding and the total number of binding sites for a transcription factor in the mammalian genome.

Publication Date: 2010 PMID: 20158869
Authors: Kuznetsov, V. A. - Singh, O. - Jenjaroenpun, P.
Journal: BMC Genomics

BACKGROUND: Transcription factor (TF)-DNA binding loci are explored by analyzing massive datasets generated with application of Chromatin Immuno-Precipitation (ChIP)-based high-throughput sequencing technologies. These datasets suffer from a bias in the information about binding loci availability, sample incompleteness and diverse sources of technical and biological noises. Therefore adequate mathematical models of ChIP-based high-throughput assay(s) and statistical tools are required for a robust identification of specific and reliable TF binding sites (TFBS), a precise characterization of TFBS avidity distribution and a plausible estimation the total number of specific TFBS for a given TF in the genome for a given cell type. RESULTS: We developed an exploratory mixture probabilistic model for a specific and non-specific transcription factor-DNA (TF-DNA) binding. Within ChiP-seq data sets, the statistics of specific and non-specific DNA-protein binding is defined by a mixture of sample size-dependent skewed functions described by Kolmogorov-Waring (K-W) function (Kuznetsov, 2003) and exponential function, respectively. Using available Chip-seq data for eleven TFs, essential for self-maintenance and differentiation of mouse embryonic stem cells (SC) (Nanog, Oct4, sox2, KLf4, STAT3, E2F1, Tcfcp211, ZFX, n-Myc, c-Myc and Essrb) reported in Chen et al (2008), we estimated (i) the specificity and the sensitivity of the ChiP-seq binding assays and (ii) the number of specific but not identified in the current experiments binding sites (BSs) in the genome of mouse embryonic stem cells. Motif finding analysis applied to the identified c-Myc TFBSs supports our results and allowed us to predict many novel c-Myc target genes. CONCLUSION: We provide a novel methodology of estimating the specificity and the sensitivity of TF-DNA binding in massively paralleled ChIP sequencing (ChIP-seq) binding assay. Goodness-of fit analysis of K-W functions suggests that a large fraction of low- and moderate- avidity TFBSs cannot be identified by the ChIP-based methods. Thus the task to identify the binding sensitivity of a TF cannot be technically resolved yet by current ChIP-seq, compared to former experimental techniques. Considering our improvement in measuring the sensitivity and the specificity of the TFs obtained from the ChIP-seq data, the models of transcriptional regulatory networks in embryonic cells and other cell types derived from the given ChIp-seq data should be carefully revised.

post to: CiteULike

Modeling neutral evolution of Alu elements using a branching process.

Publication Date: 2010 PMID: 20158868
Authors: Kimmel, M. - Mathaes, M.
Journal: BMC Genomics

BACKGROUND: Alu elements occupy about eleven percent of the human genome and are still growing in copy numbers. Since Alu elements substantially impact the shape of our genome, there is a need for modeling the amplification, mutation and selection forces of these elements. METHODS: Our proposed theoretical neutral model follows a discrete-time branching process described by Griffiths and Pakes. From this model, we derive a limit frequency spectrum of the Alu element distribution, which serves as the theoretical, neutral frequency to which real Alu insertion data can be compared through statistical goodness of fit tests. Departures from the neutral frequency spectrum may indicate selection. RESULTS: A comparison of the Alu sequence data, obtained by courtesy of Dr. Jerzy Jurka, with our model shows that the distributions of Alu sequences in the AluY family systematically deviate from the expected distribution derived from the branching process. CONCLUSIONS: This observation suggests that Alu sequences do not evolve neutrally and might be under selection.

post to: CiteULike

Transposon identification using profile HMMs.

Publication Date: 2010 PMID: 20158867
Authors: Edlefsen, P. T. - Liu, J. S.
Journal: BMC Genomics

BACKGROUND: Transposons are "jumping genes" that account for large quantities of repetitive content in genomes. They are known to affect transcriptional regulation in several different ways, and are implicated in many human diseases. Transposons are related to microRNAs and viruses, and many genes, pseudogenes, and gene promoters are derived from transposons or have origins in transposon-induced duplication. Modeling transposon-derived genomic content is difficult because they are poorly conserved. Profile hidden Markov models (profile HMMs), widely used for protein sequence family modeling, are rarely used for modeling DNA sequence families. The algorithm commonly used to estimate the parameters of profile HMMs, Baum-Welch, is prone to prematurely converge to local optima. The DNA domain is especially problematic for the Baum-Welch algorithm, since it has only four letters as opposed to the twenty residues of the amino acid alphabet. RESULTS: We demonstrate with a simulation study and with an application to modeling the MIR family of transposons that two recently introduced methods, Conditional Baum-Welch and Dynamic Model Surgery, achieve better estimates of the parameters of profile HMMs across a range of conditions. CONCLUSIONS: We argue that these new algorithms expand the range of potential applications of profile HMMs to many important DNA sequence family modeling problems, including that of searching for and modeling the virus-like transposons that are found in all known genomes.

post to: CiteULike

Mechanisms of chromosomal rearrangement in the human genome.

Publication Date: 2010 PMID: 20158866
Authors: Tsai, A. G. - Lieber, M. R.
Journal: BMC Genomics

Many human cancers are associated with characteristic chromosomal rearrangements, especially hematopoietic cancers such as leukemias and lymphomas. The first and most critical step in the rearrangement process is the induction of two DNA double-strand breaks (DSB). In all cases, at least one of the two DSBs is generated by a pathologic process, such as (1) randomly-positioned breaks due to ionizing radiation, free radical oxidative damage, or spontaneous hydrolysis; (2) breaks associated with topoisomerase inhibitor treatment; or (3) breaks at direct or inverted repeat sequences, mediated by unidentified strand breakage mechanisms. In lymphoid cells, one of the two requisite DSBs is often physiologic, the result of V(D)J recombination or class switch recombination (CSR) at the lymphoid antigen receptor loci. The RAG complex, which causes the DSBs in V(D)J recombination, can cause (4) sequence-specific, pathologic DSBs at sites that fit the consensus of their normal V(D)J recombination signal targets; or (5) structure-specific, pathologic DSBs at regions of single- to double-strand transition. CSR occurs specifically in the B-cell lineage, and requires (6) activation-induced cytidine deaminase (AID) action at sites of single-stranded DNA, which may occur pathologically outside of the normal target loci of class switch recombination regions and somatic hypermutation (SHM) zones. Recent work proposes a seventh mechanism: the sequential action of AID and the RAG complex at CpG sites provides a coherent model for the pathologic DSBs at some of the most common sites of translocation in human lymphoma - the bcl-2 gene in follicular lymphoma and diffuse large B-cell lymphoma, and the bcl-1 gene in mantle cell lymphoma.

post to: CiteULike

Whole genome analysis of p38 SAPK-mediated gene expression upon stress.

Publication Date: 2010 Mar 1 PMID: 20187982
Authors: Ferreiro, I. - Joaquin, M. - Islam, A. - Gomez-Lopez, G. - Barragan, M. - Lombardia, L. - Dominguez, O. - Pisano, D. G. - Lopez-Bigas, N. - Nebreda, A. R. - Posas, F.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: Cells have the ability to respond and adapt to environmental changes through the activation of stress-activated protein kinases (SAPKs). Although p38 SAPK signalling is known to participate in the regulation of gene expression little is known on the molecular mechanisms used by this SAPK to regulate stress-responsive genes and the overall set of genes regulated by p38 in response to different stimuli is unknown. RESULTS: Here, we report a whole genome expression analyses on mouse embryonic fibroblasts (MEFs) treated with three different p38 SAPK activating-stresses, namely osmostress, the cytokine TNFa and the protein synthesis inhibitor anisomycin. In addition, we have analysed the contribution of p38a, the major p38 family member present in MEFs, to the overall stress-induced transcriptional response by using both a chemical inhibitor (SB203580) and p38a deficient (p38a-/-) MEFs. We have also analysed the implication of p38 SAPK in the kinetics of the gene expression response to osmostress. CONCLUSIONS: Our genome-wide transcriptional analyses shows a selective response to specific stimuli and a restricted common response of up to 20% of the stress up-regulated early genes that involves an important set of transcription factors, which might be critical for either cell adaptation or preparation for continuous extracellular changes. Interestingly, up to 85% of the up-regulated genes are under the transcriptional control of p38 SAPK. Thus, activation of p38 SAPK is critical to elicit the early gene expression program required for cell adaptation to stress.

post to: CiteULike

Chloroplast genome sequence of the moss Tortula ruralis: gene content, polymorphism, and structural arrangement relative to other green plant chloroplast genomes.

Publication Date: 2010 Feb 27 PMID: 20187961
Authors: Oliver, M. J. - Murdock, A. G. - Mishler, B. D. - Kuehl, J. V. - Boore, J. L. - Mandoli, D. F. - Everett, K. D. - Wolf, P. G. - Duffy, A. M. - Karol, K. G.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: Tortula ruralis, a widely distributed species in the moss family Pottiaceae, is increasingly used as a model organism for the study of desiccation tolerance and mechanisms of cellular repair. In this paper, we present the chloroplast genome sequence of T. ruralis, only the second published chloroplast genome for a moss, and the first for a vegetatively desiccation-tolerant plant. RESULTS: The Tortula chloroplast genome is ~123,500 bp, and differs in a number of ways from that of Physcomitrella patens, the first published moss chloroplast genome. For example, Tortula lacks the ~71 kb inversion found in the large single copy region of the Physcomitrella genome and other members of the Funariales. Also, the Tortula chloroplast genome lacks petN, a gene found in all known land plant plastid genomes. In addition, an unusual case of nucleotide polymorphism was discovered. CONCLUSIONS: In this report we have presented the first sequence of the chloroplast genome of a vegetative desiccation-tolerant plant: the moss Tortula ruralis. Although the genome differs from that of the only other sequenced moss, Physcomitrella patens, we have yet to determine the biological significance of the differences. The polymorphisms we have uncovered in the sequencing of the genome offer a rare possibility (for mosses) of the generation of DNA markers for fine-level phylogenetic studies, or to investigate individual variation within populations.

post to: CiteULike

TobEA: an atlas of tobacco gene expression from seed to senescence.

Publication Date: 2010 Feb 26 PMID: 20187945
Authors: Edwards, K. D. - Bombarely, A. - Story, G. W. - Allen, F. - Mueller, L. A. - Coates, S. A. - Jones, L.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: Transcriptomics has resulted in the development of large data sets and tools for the progression of functional genomics and systems biology in many model organisms. Currently there is no commercially available microarray to allow such expression studies in Nicotiana tabacum (tobacco). RESULTS: A custom designed Affymetrix tobacco expression microarray was generated from a set of over 40k unigenes and used to measure gene expression in 19 different tobacco samples to produce the Tobacco Expression Atlas (TobEA). TobEA provides a snap shot of the transcriptional activity for thousands of tobacco genes in different tissues throughout the lifecycle of the plant and enables the identification of the biological processes occurring in these different tissues. 772 of 2513 transcription factors previously identified in tobacco were mapped to the array, with 87% of them being expressed in at least one tissue in the atlas. Putative transcriptional networks were identified based on the co-expression of these transcription factors. Several interactions in a floral identity transcription factor network were consistent with previous results from other plant species. To broaden access and maximise the benefit of TobEA a set of tools were developed to provide researchers with expression information on their genes of interest via the Solanaceae Genomics Network (SGN) web site. The array has also been made available for public use via the Nottingham Arabidopsis Stock Centre microarray service. CONCLUSIONS: The generation of a tobacco expression microarray is an important development for research in this model plant. The data provided by TobEA represents a valuable resource for plant functional genomics and systems biology research and can be used to identify gene targets for both fundamental and applied scientific applications in tobacco.

post to: CiteULike

Hyper-expansion of large DNA segments in the genome of kuruma shrimp, Marsupenaeus japonicus.

Publication Date: 2010 Feb 26 PMID: 20187930
Authors: Koyama, T. - Asakawa, S. - Katagiri, T. - Shimizu, A. - Fagutao, F. F. - Mavichak, R. - Santos, M. D. - Fuji, K. - Sakamoto, T. - Kitakado, T. - Kondo, H. - Shimizu, N. - Aoki, T. - Hirono, I.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: Higher crustaceans (class Malacostraca) represent the most species-rich and morphologically diverse group of non-insect arthropods and many of its members are commercially important. Although the crustacean DNA sequence information is growing exponentially, little is known about the genome organization of Malacostraca. Here, we constructed a bacterial artificial chromosome (BAC) library and performed BAC-end sequencing to provide genomic information for kuruma shrimp (Marsupenaeus japonicus), one of the most widely cultured species among crustaceans, and found the presence of a redundant sequence in the BAC library. We examined the BAC clone that includes the redundant sequence to further analyze its length, copy number and location in the kuruma shrimp genome. RESULTS: Mj024A04 BAC clone, which includes one redundant sequence, contained 27 putative genes and seemed to display a normal genomic DNA structure. Notably, of the putative genes, 3 genes encode homologous proteins to the inhibitor of apoptosis protein and 7 genes encode homologous proteins to white spot syndrome virus, a virulent pathogen known to affect crustaceans. Colony hybridization and PCR analysis of 381 BAC clones showed that almost half of the BAC clones maintain DNA segments whose sequences are homologous to the representative BAC clone Mj024A04. The Mj024A04 partial sequence was detected multiple times in the kuruma shrimp nuclear genome with a calculated copy number of at least 100. Microsatellites based BAC genotyping clearly showed that Mj024A04 homologous sequences were cloned from at least 48 different chromosomal loci. The absence of micro-syntenic relationships with the available genomic sequences of Daphnia and Drosophila suggests the uniqueness of these fragments in kuruma shrimp from current arthropod genome sequences. CONCLUSIONS: Our results demonstrate that hyper-expansion of large DNA segments took place in the kuruma shrimp genome. Although we analyzed only a part of the duplicated DNA segments, our result suggested that it is difficult to analyze the shrimp genome following normal analytical methodology. Hence, it is necessary to avoid repetitive sequence (such as segmental duplications) when studying the other unique structures in the shrimp genome.

post to: CiteULike

Gene-specific FACS sorting method for target selection in high-throughput amplicon sequencing.

Publication Date: 2010 Feb 26 PMID: 20184782
Authors: Sandberg, J. - Neiman, M. G. - Ahmadian, A. - Lundeberg, J.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: In addition to shotgun sequencing, next generation sequencing has been shown to be suitable for deep sequencing of many specific PCR-amplified target genes in parallel. However, unspecific product formation is a common problem in amplicon sequencing since these fragments are difficult to fully remove by gel purification, and their presence inevitably reduces the number of target sequence reads that can be obtained in each sequencing run. RESULTS: We have used a novel flow cytometric sorting approach to specifically enrich Roche/454 DNA Capture beads carrying target DNA sequences on their surface, and reject beads carrying unspecific sequences. This procedure gives a nearly three-fold increase in the fraction of informative sequences obtained. Presented results also show that there are no significant differences in the distribution or presence of different genotypes between a FACS-enriched sample and a standard-enriched control sample. CONCLUSIONS: Target-specific FACS enrichment prior to Roche/454 sequencing provides a quick, inexpensive way of increasing the amount of high quality data obtained in a single sequencing run, without introducing any sequence bias.

post to: CiteULike

Comparative assessment of methods for estimating individual genome-wide homozygosity-by-descent from human genomic data.

Publication Date: 2010 Feb 25 PMID: 20184767
Authors: Polasek, O. - Hayward, C. - Bellenguez, C. - Vitart, V. - Kolcic, I. - McQuillan, R. - Saftic, V. - Gyllensten, U. - Wilson, J. F. - Rudan, I. - Wright, A. F. - Campbell, H. - Leutenegger, A. L.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: Genome-wide homozygosity estimation from genomic data is becoming an increasingly interesting research topic. The aim of this study was to compare different methods for estimating individual homozygosity-by-descent based on the information from human genome-wide scans rather than genealogies. We considered the four most commonly used methods and investigated their applicability to single-nucleotide polymorphism (SNP) data in both a simulation study and by using the human genotyped data. A total of 986 inhabitants from the isolated Island of Vis, Croatia (where inbreeding is present, but no pedigree-based inbreeding was observed at the level of F>0.0625) were included in this study. All individuals were genotyped with the Illumina HumanHap300 array with 317,503 SNP markers. RESULTS: Simulation data suggested that multi-point FEstim is the method most strongly correlated to true homozygosity-by-descent. Correlation coefficients between the homozygosity-by-descent estimates were high but only for inbred individuals, with nearly absolute correlation between single-point measures. CONCLUSIONS: Deciding who is really inbred is a methodological challenge where multi-point approaches can be very helpful once the set of SNP markers is filtered to remove linkage disequilibrium. The use of several different methodological approaches and hence different homozygosity measures can help to distinguish between homozygosity-by-state and homozygosity-by-descent in studies investigating the effects of genomic autozygosity on human health.

post to: CiteULike

Escherichia coli infection induces distinct local and systemic transcriptome responses in the mammary gland.

Publication Date: 2010 Feb 25 PMID: 20184744
Authors: Mitterhuemer, S. - Petzl, W. - Krebs, S. - Mehne, D. - Klanner, A. - Wolf, E. - Zerbe, H. - Blum, H.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: Coliform bacteria are the most common etiologic agents in severe mastitis of cows. Escherichia coli infections are mostly restricted to a single udder quarter. Neighboring quarters stay clinically inapparent, implicating the presence of a systemic defense reaction. To address its underlying mechanism, we performed a transcriptome study of mammary tissue from udder quarters inoculated with E. coli (6 h and 24 h post infection), from neighboring quarters of the same animals, and from untreated control animals. RESULTS: After 6 h 13 differentially expressed genes (DEG) were detected in infected quarters versus control animals. Eighteen hours later 2103 and 453 DEG were found in infected and in neighboring quarters vs. control animals. Cluster analysis revealed DEG found only in infected quarters (local response) and DEG detected in both infected and neighboring quarters (systemic response). The first group includes genes mainly involved in immune response and inflammation, while the systemic reaction comprises antigen processing and presentation, cytokines, chaperones, the p38 pathway, protein degradation and apoptosis. Enhanced expression of antimicrobial genes (S100A8, S100A9, S100A12, CXCL2, GNLY), acute phase genes (LBP, SAA1, CFB, C6, IF) and indicators of oxidative stress (GPX3, MT2A, SOD2) point to an active defense reaction in infected and neighboring healthy quarters. Its early onset is indicated by increased transcription of NFIL3 at 6 h. NFIL3 is a predicted regulator of many genes of the systemic response at 24 h. The significance of our transcriptome study was evidenced by some recent findings with candidate gene based approaches. CONCLUSIONS: The discovery and holistic analysis of an extensive systemic reaction in the mammary gland significantly expands the knowledge of host-pathogen interactions in mastitis which may be relevant for the development of novel therapies and for genetic selection towards mastitis resistance.

post to: CiteULike

Genome-wide analysis of aberrant methylation in human breast cancer cells using methyl-DNA immunoprecipitation combined with high-throughput sequencing.

Publication Date: 2010 Feb 25 PMID: 20181289
Authors: Ruike, Y. - Imanaka, Y. - Sato, F. - Shimizu, K. - Tsujimoto, G.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: Cancer cells undergo massive alterations to their DNA methylation patterns that result in aberrant gene expression and malignant phenotypes. However, the mechanisms that underlie methylome changes are not well understood nor is the genomic distribution of DNA methylation changes well characterized. RESULTS: Here, we performed methylated DNA immunoprecipitation combined with high-throughput sequencing (MeDIP-seq) to obtain whole-genome DNA methylation profiles for eight human breast cancer cell (BCC) lines and for normal human mammary epithelial cells (HMEC). The MeDIP-seq analysis generated non-biased DNA methylation maps by covering almost the entire genome with sufficient depth and resolution. The most prominent feature of the BCC lines compared to HMEC was a massively reduced methylation level particularly in CpG-poor regions. While hypomethylation did not appear to be associated with particular genomic features, hypermethylation preferentially occurred at CpG-rich gene-related regions independently of the distance from transcription start sites. We also investigated methylome alterations during epithelial-to-mesenchymal transition (EMT) in MCF7 cells. EMT induction was associated with specific alterations to the methylation patterns of gene-related CpG-rich regions, although overall methylation levels were not significantly altered. Moreover, approximately 40% of the epithelial cell-specific methylation patterns in gene-related regions were altered to those typical of mesenchymal cells, suggesting a cell-type specific regulation of DNA methylation. CONCLUSIONS: This study provides the most comprehensive analysis to date of the methylome of human mammary cell lines and has produced novel insights into the mechanisms of methylome alteration during tumorigenesis and the interdependence between DNA methylome alterations and morphological changes.

post to: CiteULike

Flux of transcript patterns during soybean seed development.

Publication Date: 2010 Feb 24 PMID: 20181280
Authors: Jones, S. I. - Gonzalez, D. O. - Vodkin, L. O.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: To understand gene expression networks leading to functional properties of the soybean seed, we have undertaken a detailed examination of soybean seed development during the stages of major accumulation of oils, proteins, and starches, as well as the desiccating and mature stages, using microarrays consisting of up to 27,000 soybean cDNAs. A subset of these genes on a highly-repetitive 70-mer oligonucleotide microarray was also used to support the results. RESULTS: It was discovered that genes related to cell growth and maintenance processes, as well as energy processes like photosynthesis, decreased in expression levels as the cotyledons approached the mature, dry stage. Genes involved with some storage proteins had their highest expression levels at the stage of highest fresh weight. However, genes encoding many transcription factors and DNA binding proteins showed higher expression levels in the desiccating and dry seeds than in most of the green stages. CONCLUSIONS: Data on 27,000 cDNAs have been obtained over five stages of soybean development, including the stages of major accumulation of agronomically-important products, using two different types of microarrays. Of particular interest are the genes found to peak in expression at the desiccating and dry seed stages, such as those annotated as transcription factors, which may indicate the preparation of pathways that will be needed later in the early stages of imbibition and germination.

post to: CiteULike

A consensus linkage map of the grass carp (Ctenopharyngodon idella) based on microsatellites and SNPs.

Publication Date: 2010 Feb 24 PMID: 20181260
Authors: Xia, J. H. - Liu, F. - Zhu, Z. Y. - Fu, J. - Feng, J. - Li, J. - Yue, G. H.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: Grass carp (Ctenopharyngodon idella) belongs to the family Cyprinidae which includes more than 2000 fish species. It is one of the most important freshwater food fish species in world aquaculture. A linkage map is an essential framework for mapping traits of interest and is often the first step towards understanding genome evolution. The aim of this study is to construct a first generation genetic map of grass carp using microsatellites and SNPs to generate a new resource for mapping QTL for economically important traits and to conduct a comparative mapping analysis to shed new insights into the evolution of fish genomes. RESULTS: We constructed a first generation linkage map of grass carp with a mapping panel containing two F1 families including 192 progenies. Sixteen SNPs in genes and 263 microsatellite markers were mapped to twenty-four linkage groups (LGs). The number of LGs was corresponding to the haploid chromosome number of grass carp. The sex-specific map was 1149.4 and 888.8 cM long in females and males respectively whereas the sex-averaged map spanned 1176.1 cM. The average resolution of the map was 4.2 cM/locus. BLAST searches of sequences of mapped markers of grass carp against the whole genome sequence of zebrafish revealed substantial macrosynteny relationship and extensive colinearity of markers between grass carp and zebrafish. CONCLUSIONS: The linkage map of grass carp presented here is the first linkage map of a food fish species based on co-dominant markers in the family Cyprinidae. This map provides a valuable resource for mapping phenotypic variations and serves as a reference to approach comparative genomics and understand the evolution of fish genomes and could be complementary to grass carp genome sequencing project.

post to: CiteULike

Correcting for intra-experiment variation in Illumina BeadChip data is necessary to generate robust gene-expression profiles.

Publication Date: 2010 Feb 24 PMID: 20181233
Authors: Kitchen, R. R. - Sabine, V. S. - Sims, A. H. - Macaskill, E. J. - Renshaw, L. - Thomas, J. S. - van Hemert, J. I. - Dixon, J. M. - Bartlett, J. M.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: Microarray technology is a popular means of producing whole genome transcriptional profiles, however high cost and scarcity of mRNA has led many studies to be conducted based on the analysis of single samples. We exploit the design of the Illumina platform, specifically multiple arrays on each chip, to evaluate intra-experiment technical variation using repeated hybridisations of universal human reference RNA (UHRR) and duplicate hybridisations of primary breast tumour samples from a clinical study. RESULTS: A clear batch-specific bias was detected in the measured expressions of both the UHRR and clinical samples. This bias was found to persist following standard microarray normalisation techniques. However, when mean-centering or empirical Bayes batch-correction methods (ComBat) were applied to the data, inter-batch variation in the UHRR and clinical samples were greatly reduced. Correlation between replicate UHRR samples improved by two orders of magnitude following batch-correction using ComBat (ranging from 0.9833-0.9991 to 0.9997-0.9999) and increased the consistency of the gene-lists from the duplicate clinical samples, from 11.6% in quantile normalised data to 66.4% in batch-corrected data. The use of UHRR as an inter-batch calibrator provided a small additional benefit when used in conjunction with ComBat, further increasing the agreement between the two gene-lists, up to 74.1%. CONCLUSION: In the interests of practicalities and cost, these results suggest that single samples can generate reliable data, but only after careful compensation for technical bias in the experiment. We recommend that investigators appreciate the propensity for such variation in the design stages of a microarray experiment and that the use of suitable correction methods become routine during the statistical analysis of the data.

post to: CiteULike

Comparative analysis of fungal protein kinases and associated domains.

Publication Date: 2010 Feb 24 PMID: 20178650
Authors: Kosti, I. - Mandel-Gutfreund, Y. - Glaser, F. - Horwitz, B. A.
Journal: BMC Genomics

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.

post to: CiteULike

Pediatric primary central nervous system germ cell tumors of different prognosis groups show characteristic miRNome traits and chromosome copy number variations.

Publication Date: 2010 Feb 24 PMID: 20178649
Authors: Wang, H. W. - Wu, Y. H. - Hsieh, J. Y. - Liang, M. L. - Chao, M. E. - Liu, D. J. - Hsu, M. T. - Wong, T. T.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: Intracranial pediatric germ cell tumors (GCTs) are rare and heterogeneous neoplasms and vary in histological differentiation, prognosis and clinical behavior. Germinoma and mature teratoma are GCTs that have a good prognosis, while other types of GCTs, termed nongerminomatous malignant germ cell tumors (NGMGCTs), are tumors with an intermediate or poor prognosis. The second group of tumors requires more extensive drug and irradiation treatment regimens. The mechanisms underlying the differences in incidence and prognosis of the various GCT subgroups are unclear. RESULTS: We identified a distinct mRNA profile correlating with GCT histological differentiation and prognosis, and also present in this study the first miRNA profile of pediatric primary intracranial GCTs. Most of the differentially expressed miRNAs were downregulated in germinomas, but miR-142-5p and miR-146a were upregulated. Genes responsible for self-renewal (such as POU5F1 (OCT4), NANOG and KLF4) and the immune response were abundant in germinomas, while genes associated with neuron differentiation, Wnt/beta-catenin pathway, invasiveness and epithelial-mesenchymal transition (including SNAI2 (SLUG) and TWIST2) were abundant in NGMGCTs. Clear transcriptome segregation based on patient survival was observed, with malignant NGMGCTs being closest to embryonic stem cells. Chromosome copy number variations (CNVs) at cytobands 4q13.3-4q28.3 and 9p11.2-9q13 correlated with GCT malignancy and clinical risk. Six genes (BANK1, CXCL9, CXCL11, DDIT4L, ELOVL6 and HERC5) within 4q13.3-4q28.3 were more abundant in germinomas. CONCLUSIONS: Our results integrate molecular profiles with clinical observations and provide insights into the underlying mechanisms causing GCT malignancy. The genes, pathways and microRNAs identified have the potential to be novel therapeutic targets.

post to: CiteULike

Genomic sequence of a mutant strain of Caenorhabditis elegans with an altered recombination pattern.

Publication Date: 2010 Feb 23 PMID: 20178641
Authors: Rose, A. M. - O'Neil, N. J. - Bilenky, M. - Butterfield, Y. S. - Malhis, N. - Flibotte, S. - Jones, M. R. - Marra, M. - Baillie, D. L. - Jones, S. J.
Journal: BMC Genomics

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.

post to: CiteULike

Identification and characterization of Tc1/mariner-like DNA transposons in genomes of the pathogenic fungi of the Paracoccidioides species complex.

Publication Date: 2010 Feb 23 PMID: 20178623
Authors: Marini, M. M. - Zanforlin, T. - Santos, P. C. - Barros, R. R. - Guerra, A. C. - Puccia, R. - Felipe, M. S. - Brigido, M. - Soares, C. M. - Ruiz, J. C. - Silveira, J. F. - Cisalpino, P. S.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: Paracoccidioides brasiliensis (Eukaryota, Fungi, Ascomycota) is a thermodimorphic fungus, the etiological agent of paracoccidioidomycosis, the most important systemic mycoses in Latin America. Three isolates corresponding to distinct phylogenetic lineages of the Paracoccidioides species complex had their genomes sequenced. In this study the identification and characterization of class II transposable elements in the genomes of these fungi was carried out. RESULTS: A genomic survey for DNA transposons in the sequence assemblies of Paracoccidioides, a genus recently proposed to encompass species P. brasiliensis (harboring phylogenetic lineages S1, PS2, PS3) and P. lutzii (Pb01-like isolates), has been completed. Eight new Tc1/mariner families, referred to as Trem (Transposable element mariner), labeled A through H were identified. Elements from each family have 65-80% sequence similarity with other Tc1/mariner elements. They are flanked by 2-bp TA target site duplications and different termini. Encoded DDD-transposases, some of which have complete ORFs, indicated that they could be functionally active. The distribution of Trem elements varied between the genomic sequences characterized as belonging to P. brasiliensis (S1 and PS2) and P. lutzii. TremC and H elements would have been present in a hypothetical ancestor common to P. brasiliensis and P. lutzii, while TremA, B and F elements were either acquired by P. brasiliensis or lost by P. lutzii after speciation. Although TremD and TremE share about 70% similarity, they are specific to P. brasiliensis and P. lutzii, respectively. This suggests that these elements could either have been present in a hypothetical common ancestor and have evolved divergently after the split between P. brasiliensis and P. Lutzii, or have been independently acquired by horizontal transfer. CONCLUSIONS: New families of Tc1/mariner DNA transposons in the genomic assemblies of the Paracoccidioides species complex are described. Families were distinguished based on significant BLAST identities between transposases and/or TIRs. The expansion of Trem in a putative ancestor common to the species P. brasiliensis and P. lutzii would have given origin to TremC and TremH, while other elements could have been acquired or lost after speciation had occurred. The results may contribute to our understanding of the organization and architecture of genomes in the genus Paracoccidioides.

post to: CiteULike

High-density linkage mapping and evolution of paralogs and orthologs in Salix and Populus.

Publication Date: 2010 PMID: 20178595
Authors: Berlin, S. - Lagercrantz, U. - von Arnold, S. - Ost, T. - Ronnberg-Wastljung, A. C.
Journal: BMC Genomics

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, S1, contains 495 markers with a total genetic distance of 2477 cM and an average distance of 5.0 cM between the markers. The S3 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 x 10-9 and 1.68 x 10-9 per site and year, which are close to rates found in other perennials but much lower than rates in annuals.

post to: CiteULike

Defining species specific genome differences in malaria parasites.

Publication Date: 2010 Feb 23 PMID: 20175934
Authors: Liew, K. J. - Hu, G. - Bozdech, Z. - Preiser, P. R.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: In recent years a number of genome sequences for different plasmodium species have become available. This has allowed the identification of numerous conserved genes across the different species and has significantly enhanced our understanding of parasite biology. In contrast little is known about species specific differences between the different genomes partly due to the lower sequence coverage and therefore relatively poor annotation of some of the draft genomes particularly the rodent malarias parasite species. RESULTS: To improve the current annotation and gene identification status of the draft genomes of P. berghei, P. chabaudi and P. yoelii, we performed genome-wide comparisons between these three species. Through analyses via comparative genome hybridizations using a newly designed pan-rodent array as well as in depth bioinformatic analysis, we were able to improve on the coverage of the draft rodent parasite genomes by detecting orthologous genes between these related rodent parasite species. More than 1,000 orthologs for P. yoelii were now newly associated with a P. falciparum gene. In addition to extending the current core gene set for all plasmodium species this analysis also for the first time identifies a relatively small number of genes that are unique to the primate malaria parasites while a larger gene set is uniquely conserved amongst the rodent malaria parasites. CONCLUSIONS: These findings allow a more thorough investigation of the genes that are important for host specificity in malaria.

post to: CiteULike

Expression profiling of mouse embryonic fibroblasts with a deletion in the helicase domain of the Werner Syndrome gene homologue treated with hydrogen peroxide.

Publication Date: 2010 Feb 22 PMID: 20175907
Authors: Labbe, A. - Turaga, R. V. - Paquet, E. R. - Garand, C. - Lebel, M.
Journal: BMC Genomics

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 Wrndeltahel/deltahel). RESULTS: Wrndeltahel/deltahel 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 Wrndeltahel/deltahel cells. Importantly, untreated Wrndeltahel/deltahel 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 Wrndeltahel/deltahel 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 Wrndeltahel/deltahel mutant cells already exhibited a modest chronic intracellular oxidative stress.

post to: CiteULike

Transcriptome screen for fast evolving genes by Inter-Specific Selective Hybridization (ISSH).

Publication Date: 2010 Feb 22 PMID: 20175901
Authors: Montoya-Burgos, J. I. - Foulon, A. - Bahechar, I.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: Fast evolving genes are targets of an increasing panel of biological studies, from cancer research to population genetics and species specific adaptations. Yet, their identification and isolation are still laborious, particularly for non-model organisms. We developed a method, named the Inter-Specific Selective Hybridization (ISSH) method, for generating cDNA libraries enriched in fast evolving genes. It utilizes transcripts of homologous tissues of distinct yet related species. Experimental hybridization conditions are monitored in order to discard transcripts that do not find their homologous counterparts in the two species sets as well as transcripts that display a strong complementarity between the two species. Only heteroduplexes that disanneal at low stringency are used for constructing the resulting cDNA library. RESULTS: We demonstrate the efficiency of the ISSH method by generating a brain cDNA library enriched in fast evolving transcripts of a non-model catfish species as well as a control, non-enriched library. Our results indicate that the enriched library contains effectively more fast evolving sequences than the control library. Gene annotation analyses also indicate enrichment in genes with low expression levels and non-ubiquitously expressed genes, both categories encompassing the majority of fast evolving genes. Furthermore, most of the identified transcripts show higher sequence divergence between two closely related catfish species as compared to recognized fast evolving DNA markers. CONCLUSIONS: The ISSH method offers a simple, inexpensive and efficient way to screen the transcriptome for isolating fast evolving genes. This method opens new opportunities in the investigation of biological mechanisms that include fast evolving genes, such as the evolution of lineage specific processes and traits responsible for species adaptation to their environment.

post to: CiteULike

Comparative gene expression profiling between human cultured myotubes and skeletal muscle tissue.

Publication Date: 2010 Feb 22 PMID: 20175888
Authors: Raymond, F. - Metairon, S. - Kussmann, M. - Colomer, J. - Nascimento, A. - Mormeneo, E. - Garcia-Martinez, C. - Gomez-Foix, A. M.
Journal: BMC Genomics

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.

post to: CiteULike

Digital gene expression analysis of two life cycle stages of the human-infective parasite, Trypanosoma brucei gambiense reveals differentially expressed clusters of co-regulated genes.

Publication Date: 2010 Feb 22 PMID: 20175885
Authors: Veitch, N. J. - Johnson, P. C. - Trivedi, U. - Terry, S. - Wildridge, D. - Macleod, A.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: The evolutionarily ancient parasite, Trypanosoma brucei, is unusual in that the majority of its genes are regulated post-transcriptionally, leading to the suggestion that transcript abundance of most genes does not vary significantly between different life cycle stages despite the fact that the parasites undergo substantial cellular remodelling and metabolic changes throughout its complex life cycle. To investigate this in the clinically relevant sub-species, Trypanosoma brucei gambiense, which is the causative agent of the fatal human disease African sleeping sickness, we have compared the transcriptome of two different life cycle stages, the potentially human-infective bloodstream forms with the non-human-infective procyclic stage using digital gene expression (DGE) analysis. RESULTS: Over eleven million unique tags were generated, producing expression data for 7360 genes, covering 81% of the genes in the genome. Compared to microarray analysis of the related T. b. brucei parasite, approximately 10 times more genes with a 2.5 fold change in expression levels were detected. The transcriptome analysis revealed the existence of several differentially expressed gene clusters within the genome, indicating that contiguous genes presumably from the same polycistronic unit are co-regulated either at the level of transcription or transcript stability. CONCLUSIONS: DGE analysis is extremely sensitive for detecting gene expression differences, revealing firstly that a far greater number of genes are stage-regulated than had previously been identified and secondly and more importantly, this analysis has revealed the existence of several differentially expressed clusters of genes present on what appears to be the same polycistronic units, a phenomenon which had not previously been observed in microarray studies. These differentially regulated clusters of genes are in addition to the previously identified RNA polymerase I polycistronic units of variant surface glycoproteins and procyclin expression sites, which encode the major surface proteins of the parasite. This raises a number of questions regarding the function and regulation of the gene clusters that clearly warrant further study.

post to: CiteULike

Structural classification by the Lipase Engineering Database: a case study of Candida antarctica lipase A.

Publication Date: 2010 Feb 19 PMID: 20170513
Authors: Widmann, M. - Juhl, P. B. - Pleiss, J.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: The Lipase Engineering Database (LED) integrates information on sequence, structure and function of lipases, esterases and related proteins with the alpha/beta 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 alpha/beta 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 alpha/beta 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 alpha/beta hydrolase fold. By specifically comparing the substrate binding pocket of CALA to other binding pockets of alpha/beta 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 alpha/beta 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.

post to: CiteULike

Feasibility of physical map construction from fingerprinted bacterial artificial chromosome libraries of polyploid plant species.

Publication Date: 2010 Feb 19 PMID: 20170511
Authors: Luo, M. C. - Ma, Y. - You, F. M. - Anderson, O. D. - Kopecky, D. - Simkova, H. - Safar, J. - Dolezel, J. - Gill, B. S. - McGuire, P. E. - Dvorak, J.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: The presence of closely related genomes in polyploid species makes the assembly of total genomic sequence from shotgun sequence reads produced by the current sequencing platforms exceedingly difficult, if not impossible. Genomes of polyploid species could be sequenced following the ordered-clone sequencing approach employing contigs of bacterial artificial chromosome (BAC) clones and BAC-based physical maps. Although BAC contigs can currently be constructed for virtually any diploid organism with the SNaPshot high-information-content-fingerprinting (HICF) technology, it is currently unknown if this is also true for polyploid species. It is possible that BAC clones from orthologous regions of homoeologous chromosomes would share numerous restriction fragments and be therefore included into common contigs. Because of this and other concerns, physical mapping utilizing the SNaPshot HICF of BAC libraries of polyploid species has not been pursued and the possibility of doing so has not been assessed. The sole exception has been in common wheat, an allohexaploid in which it is possible to construct single-chromosome or single-chromosome-arm BAC libraries from DNA of flow-sorted chromosomes and bypass the obstacles created by polyploidy. RESULTS: The potential of the SNaPshot HICF technology for physical mapping of polyploid plants utilizing global BAC libraries was evaluated by assembling contigs of fingerprinted clones in an in silico merged BAC library composed of single-chromosome libraries of two wheat homoeologous chromosome arms, 3AS and 3DS, and complete chromosome 3B. Because the chromosome arm origin of each clone was known, it was possible to estimate the fidelity of contig assembly. On average 97.78% or more clones, depending on the library, were from a single chromosome arm. A large portion of the remaining clones was shown to be library contamination from other chromosomes, a feature that is unavoidable during the construction of single-chromosome BAC libraries. CONCLUSIONS: The negligibly low level of incorporation of clones from homoeologous chromosome arms into a contig during contig assembly suggested that it is feasible to construct contigs and physical maps using global BAC libraries of wheat and almost certainly also of other plant polyploid species with genome sizes comparable to that of wheat. Because of the high purity of the resulting assembled contigs, they can be directly used for genome sequencing. It is currently unknown but possible that equally good BAC contigs can be also constructed for polyploid species containing smaller, more gene-rich genomes.

post to: CiteULike

Massive gene losses in Asian cultivated rice unveiled by comparative genome analysis.

Publication Date: 2010 PMID: 20167122
Authors: Sakai, H. - Itoh, T.
Journal: BMC Genomics

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.

post to: CiteULike

High-throughput genome sequencing of two Listeria monocytogenes clinical isolates during a large foodborne outbreak.

Publication Date: 2010 PMID: 20167121
Authors: Gilmour, M. W. - Graham, M. - Van Domselaar, G. - Tyler, S. - Kent, H. - Trout-Yakel, K. M. - Larios, O. - Allen, V. - Lee, B. - Nadon, C.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: A large, multi-province outbreak of listeriosis associated with ready-to-eat meat products contaminated with Listeria monocytogenes serotype 1/2a occurred in Canada in 2008. Subtyping of outbreak-associated isolates using pulsed-field gel electrophoresis (PFGE) revealed two similar but distinct AscI PFGE patterns. High-throughput pyrosequencing of two L. monocytogenes isolates was used to rapidly provide the genome sequence of the primary outbreak strain and to investigate the extent of genetic diversity associated with a change of a single restriction enzyme fragment during PFGE. RESULTS: The chromosomes were collinear, but differences included 28 single nucleotide polymorphisms (SNPs) and three indels, including a 33 kbp prophage that accounted for the observed difference in AscI PFGE patterns. The distribution of these traits was assessed within further clinical, environmental and food isolates associated with the outbreak, and this comparison indicated that three distinct, but highly related strains may have been involved in this nationwide outbreak. Notably, these two isolates were found to harbor a 50 kbp putative mobile genomic island encoding translocation and efflux functions that has not been observed in other Listeria genomes. CONCLUSIONS: High-throughput genome sequencing provided a more detailed real-time assessment of genetic traits characteristic of the outbreak strains than could be achieved with routine subtyping methods. This study confirms that the latest generation of DNA sequencing technologies can be applied during high priority public health events, and laboratories need to prepare for this inevitability and assess how to properly analyze and interpret whole genome sequences in the context of molecular epidemiology.

post to: CiteULike

Identification of microRNAs expressed in two mosquito vectors, Aedes albopictus and Culex quinquefasciatus.

Publication Date: 2010 PMID: 20167119
Authors: Skalsky, R. L. - Vanlandingham, D. L. - Scholle, F. - Higgs, S. - Cullen, B. R.
Journal: BMC Genomics

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.

post to: CiteULike

Increased incidence of rare codon clusters at 5' and 3' gene termini:implications for function.

Publication Date: 2010 PMID: 20167116
Authors: Clarke, T. F. 4th - Clark, P. L.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: The process of translation can be affected by the use of rare versus common codons within the mRNA transcript. RESULTS: Here, we show that rare codons are enriched at the 5' and 3' termini of genes from E. coli and other prokaryotes. Genes predicted to be secreted show significant enrichment in 5' rare codon clusters, but not 3' rare codon clusters. Surprisingly, no correlation between 5' mRNA structure and rare codon usage was observed. CONCLUSIONS: Potential functional roles for the enrichment of rare codons at terminal positions are explored.

post to: CiteULike

General and species-specific transcriptional responses to downy mildew infection in a susceptible (Vitis vinifera) and a resistant (V. riparia) grapevine species.

Publication Date: 2010 PMID: 20167053
Authors: Polesani, M. - Bortesi, L. - Ferrarini, A. - Zamboni, A. - Fasoli, M. - Zadra, C. - Lovato, A. - Pezzotti, M. - Delledonne, M. - Polverari, A.
Journal: BMC Genomics

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.

post to: CiteULike

Genomic sequencing and analyses of Lymantria xylina multiple nucleopolyhedrovirus.

Publication Date: 2010 PMID: 20167051
Authors: Nai, Y. S. - Wu, C. Y. - Wang, T. C. - Chen, Y. R. - Lau, W. H. - Lo, C. F. - Tsai, M. F. - Wang, C. H.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: Outbreaks of the casuarina moth, Lymantria xylina Swinehoe (Lepidoptera: Lymantriidae), which is a very important forest pest in Taiwan, have occurred every five to 10 years. This moth has expanded its range of host plants to include more than 65 species of broadleaf trees. LyxyMNPV (L. xylina multiple nucleopolyhedrovirus) is highly virulent to the casuarina moth and has been investigated as a possible biopesticide for controlling this moth. LdMNPV-like virus has also been isolated from Lymantria xylina larvae but LyxyMNPV was more virulent than LdMNPV-like virus both in NTU-LY and IPLB-LD-652Y cell lines. To better understand LyxyMNPV, the nucleotide sequence of the LyxyMNPV DNA genome was determined and analysed. RESULTS: The genome of LyxyMNPV consists of 156,344 bases, has a G+C content of 53.4% and contains 157 putative open reading frames (ORFs). The gene content and gene order of LyxyMNPV were similar to those of LdMNPV, with 151 ORFs identified as homologous to those reported in the LdMNPV genome. Two genes (Lyxy49 and Lyxy123) were homologous to other baculoviruses, and four unique LyxyMNPV ORFs (Lyxy11, Lyxy19, Lyxy130 and Lyxy131) were identified in the LyxyMNPV genome, including a gag-like gene that was not reported in baculoviruses. LdMNPV contains 23 ORFs that are absent in LyxyMNPV. Readily identifiable homologues of the gene host range factor-1 (hrf-1), which appears to be involved in the susceptibility of L. dispar to NPV infection, were not present in LyxyMNPV. Additionally, two putative odv-e27 homologues were identified in LyxyMNPV. The LyxyMNPV genome encoded 14 bro genes compared with 16 in LdMNPV, which occupied more than 8% of the LyxyMNPV genome. Thirteen homologous regions (hrs) were identified containing 48 repeated sequences composed of 30-bp imperfect palindromes. However, they differed in the relative positions, number of repeats and orientation in the genome compared to LdMNPV. CONCLUSION: The gene parity plot analysis, percent identity of the gene homologues and a phylogenetic analysis suggested that LyxyMNPV is a Group II NPV that is most closely related to LdMNPV but with a highly distinct genomic organisation.

post to: CiteULike

Transcriptomic changes arising during light-induced sporulation in Physarum polycephalum.

Publication Date: 2010 Feb 17 PMID: 20163733
Authors: Barrantes, I. - Glockner, G. - Meyer, S. - Marwan, W.
Journal: BMC Genomics

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.

post to: CiteULike

Genome-wide analysis of alternative splicing in Chlamydomonas reinhardtii.

Publication Date: 2010 PMID: 20163725
Authors: Labadorf, A. - Link, A. - Rogers, M. F. - Thomas, J. - Reddy, A. S. - Ben-Hur, A.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: Genome-wide computational analysis of alternative splicing (AS) in several flowering plants has revealed that pre-mRNAs from about 30% of genes undergo AS. Chlamydomonas, a simple unicellular green alga, is part of the lineage that includes land plants. However, it diverged from land plants about one billion years ago. Hence, it serves as a good model system to study alternative splicing in early photosynthetic eukaryotes, to obtain insights into the evolution of this process in plants, and to compare splicing in simple unicellular photosynthetic and non-photosynthetic eukaryotes. We performed a global analysis of alternative splicing in Chlamydomonas reinhardtii using its recently completed genome sequence and all available ESTs and cDNAs. RESULTS: Our analysis of AS using BLAT and a modified version of the Sircah tool revealed AS of 498 transcriptional units with 611 events, representing about 3% of the total number of genes. As in land plants, intron retention is the most prevalent form of AS. Retained introns and skipped exons tend to be shorter than their counterparts in constitutively spliced genes. The splice site signals in all types of AS events are weaker than those in constitutively spliced genes. Furthermore, in alternatively spliced genes, the prevalent splice form has a stronger splice site signal than the non-prevalent form. Analysis of constitutively spliced introns revealed an over-abundance of motifs with simple repetitive elements in comparison to introns involved in intron retention. In almost all cases, AS results in a truncated ORF, leading to a coding sequence that is around 50% shorter than the prevalent splice form. Using RT-PCR we verified AS of two genes and show that they produce more isoforms than indicated by EST data. All cDNA/EST alignments and splice graphs are provided in a website at http://combi.cs.colostate.edu/as/chlamy. CONCLUSIONS: The extent of AS in Chlamydomonas that we observed is much smaller than observed in land plants, but is much higher than in simple unicellular heterotrophic eukaryotes. The percentage of different alternative splicing events is similar to flowering plants. Prevalence of constitutive and alternative splicing in Chlamydomonas, together with its simplicity, many available public resources, and well developed genetic and molecular tools for this organism make it an excellent model system to elucidate the mechanisms involved in regulated splicing in photosynthetic eukaryotes.

post to: CiteULike

SoyTEdb: a comprehensive database of transposable elements in the soybean genome.

Publication Date: 2010 PMID: 20163715
Authors: Du, J. - Grant, D. - Tian, Z. - Nelson, R. T. - Zhu, L. - Shoemaker, R. C. - Ma, J.
Journal: BMC Genomics

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.

post to: CiteULike

Intensity-based analysis of dual-color gene expression data as an alternative to ratio-based analysis to enhance reproducibility.

Publication Date: 2010 Feb 17 PMID: 20163706
Authors: Bossers, K. - Ylstra, B. - Brakenhoff, R. H. - Smeets, S. J. - Verhaagen, J. - van de Wiel, M. A.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: Ratio-based analysis is the current standard for the analysis of dual-color microarray data. Indeed, this method provides a powerful means to account for potential technical variations such as differences in background signal, spot size and spot concentration. However, current high density dual-color array platforms are of very high quality, and inter-array variance has become much less pronounced. We therefore raised the question whether it is feasible to use an intensity-based analysis rather than ratio-based analysis of dual-color microarray datasets. Furthermore, we compared performance of both ratio- and intensity-based analyses in terms of reproducibility and sensitivity for differential gene expression. RESULTS: By analyzing three distinct and technically replicated datasets with either ratio- or intensity-based models, we determined that, when applied to the same dataset, intensity-based analysis of dual-color gene expression experiments yields 1) more reproducible results, and 2) is more sensitive in the detection of differentially expressed genes. These effects were most pronounced in experiments with large biological variation and complex hybridization designs. Furthermore, a power analysis revealed that for direct two-group comparisons above a certain sample size, ratio-based models have higher power, although the difference with intensity-based models is very small. CONCLUSIONS: Intensity-based analysis of dual-color datasets results in more reproducible results and increased sensitivity in the detection of differential gene expression than the analysis of the same dataset with ratio-based analysis. Complex dual-color setups such as interwoven loop designs benefit most from ignoring the array factor. The applicability of our approach to array platforms other than dual-color needs to be further investigated.

post to: CiteULike

Identification of candidates for cyclotide biosynthesis and cyclisation by expressed sequence tag analysis of Oldenlandia affinis.

Publication Date: 2010 Feb 16 PMID: 20158917
Authors: Qin, Q. - McCallum, E. J. - Kaas, Q. - Suda, J. - Saska, I. - Craik, D. J. - Mylne, J. S.
Journal: BMC Genomics

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.

post to: CiteULike

Analysis of a normalised expressed sequence tag (EST) library from a key pollinator, the bumblebee Bombus terrestris.

Publication Date: 2010 Feb 15 PMID: 20156341
Authors: Sadd, B. M. - Kube, M. - Klages, S. - Reinhardt, R. - Schmid-Hempel, P.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: The bumblebee, Bombus terrestris (Order Hymenoptera), is of widespread importance. This species is widely used for commercial pollination in Europe, and along with other Bombus spp. is a key member of natural pollinator assemblages. Furthermore, the species is studied in a wide variety of biological fields. The objective of this project was to create a B. terrestris EST resource that will prove to be valuable in obtaining a deeper understanding of this significant social insect. RESULTS: A normalised cDNA library was constructed from the thorax and abdomen of B. terrestris workers in order to enhance the discovery of rare genes. A total of 29'428 ESTs were sequenced. Subsequent clustering resulted in 13'333 unique sequences. Of these, 58.8 percent had significant similarities to known proteins, with 54.5 percent having a "best-hit" to existing Hymenoptera sequences. Comparisons with the honeybee and other insects allowed the identification of potential candidates for gene loss, pseudogene evolution and possible incomplete annotation in the honeybee genome. Further, given the focus of much basic research and the perceived threat of disease to natural and commercial populations, the immune system of bumblebees is a particularly relevant component. Although the library is derived from unchallenged bees we still uncover transcription of a number of immune genes, spanning the principally described insect immune pathways. Additionally, the EST library provides a resource for the discovery of genetic markers that can be used in population level studies. Indeed, initial screens identified 589 simple sequence repeats and 854 potential single nucleotide polymorphisms. CONCLUSION: The resource that these B. terrestris ESTs is valuable for ongoing work. The ESTs provide direct evidence of transcriptionally active regions, but they will also facilitate further functional genomics work, gene discovery and future genome annotation. These are important aspects in obtaining a greater understanding of this key pollinator species.

post to: CiteULike

Correction: High throughput approaches reveal splicing of primary microRNA transcripts and tissue specific expression of mature microRNAs in Vitis vinifera.

Publication Date: 2010 PMID: 20152027
Authors: Mica, E. - Piccolo, V. - Delledonne, M. - Ferrarini, A. - Pezzotti, M. - Casati, C. - Del Fabbro, C. - Valle, G. - Policriti, A. - Morgante, M. - Pesole, G. - Pe, M. E. - Horner, D. S.
Journal: BMC Genomics

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.

post to: CiteULike

Genomics based analysis of interactions between developing B-lymphocytes and stromal cells reveal complex interactions and two-way communication.

Publication Date: 2010 PMID: 20152025
Authors: Zetterblad, J. - Qian, H. - Zandi, S. - Mansson, R. - Lagergren, A. - Hansson, F. - Bryder, D. - Paulsson, N. - Sigvardsson, M.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: The use of functional genomics has largely increased our understanding of cell biology and promises to help the development of systems biology needed to understand the complex order of events that regulates cellular differentiation in vivo. One model system clearly dependent on the integration of extra and intra cellular signals is the development of B-lymphocytes from hematopoietic stem cells in the bone marrow. This developmental pathway involves several defined differentiation stages associated with specific expression of genes including surface markers that can be used for the prospective isolation of the progenitor cells directly from the bone marrow to allow for ex vivo gene expression analysis. The developmental process can be simulated in vitro making it possible to dissect information about cell/cell communication as well as to address the relevance of communication pathways in a rather direct manner. Thus we believe that B-lymphocyte development represents a useful model system to take the first steps towards systems biology investigations in the bone marrow. RESULTS: In order to identify extra cellular signals that promote B lymphocyte development we created a database with approximately 400 receptor ligand pairs and software matching gene expression data from two cell populations to obtain information about possible communication pathways. Using this database and gene expression data from NIH3T3 cells (unable to support B cell development), OP-9 cells (strongly supportive of B cell development), pro-B and pre-B cells as well as mature peripheral B-lineage cells, we were able to identify a set of potential stage and stromal cell restricted communication pathways. Functional analysis of some of these potential ways of communication allowed us to identify BMP-4 as a potent stimulator of B-cell development in vitro. Further, the analysis suggested that there existed possibilities for progenitor B cells to send signals to the stroma. The functional consequences of this were investigated by co-culture experiments revealing that the co-incubation of stromal cells with B cell progenitors altered both the morphology and the gene expression pattern in the stromal cells. CONCLUSIONS: We believe that this gene expression data analysis method allows for the identification of functionally relevant interactions and therefore could be applied to other data sets to unravel novel communication pathways.

post to: CiteULike

Mapping main, epistatic and sex-specific QTL for body composition in a chicken population divergently selected for low or high growth rate.

Publication Date: 2010 PMID: 20149241
Authors: Ankra-Badu, G. A. - Shriner, D. - Le Bihan-Duval, E. - Mignon-Grasteau, S. - Pitel, F. - Beaumont, C. - Duclos, M. J. - Simon, J. - Porter, T. E. - Vignal, A. - Cogburn, L. A. - Allison, D. B. - Yi, N. - Aggrey, S. E.
Journal: BMC Genomics

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 F2 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.

post to: CiteULike

U12 type introns were lost at multiple occasions during evolution.

Publication Date: 2010 Feb 11 PMID: 20149226
Authors: Bartschat, S. - Samuelsson, T.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: Two categories of introns are known, a common U2 type and a rare U12 type. These two types of introns are removed by distinct spliceosomes. The phylogenetic distribution of spliceosomal RNAs that are characteristic of the U12 spliceosome, i.e. the U11, U12, U4atac and U6atac RNAs, suggest that U12 spliceosomes were lost in many phylogenetic groups. We have now examined the distribution of U2 and U12 introns in many of these groups. RESULTS: U2 and U12 introns were predicted by making use of available EST and genomic sequences. The results show that in species or branches where U12 spliceosomal components are missing, also U12 type of introns are lacking. Examples are the choanoflagellate Monosiga brevicollis, Entamoeba histolytica, green algae, diatoms, and the fungal lineage Basidiomycota. Furthermore, whereas U12 splicing does not occur in Caenorhabditis elegans, U12 introns as well as U12 snRNAs are present in Trichinella spiralis, which is deeply branching in the nematode tree. A comparison of homologous genes in T. spiralis and C. elegans revealed different mechanisms whereby U12 introns were lost. CONCLUSIONS: The phylogenetic distribution of U12 introns and spliceosomal RNAs give further support to an early origin of U12 dependent splicing. In addition, this distribution identifies a large number of instances during eukaryotic evolution where such splicing was lost.

post to: CiteULike

Fungal Secretome Database: Integrated platform for annotation of fungal secretomes.

Publication Date: 2010 Feb 11 PMID: 20146824
Authors: Choi, J. - Park, J. - Kim, D. - Jung, K. - Kang, S. - Lee, Y. H.
Journal: BMC Genomics

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/.

post to: CiteULike

Identification and analysis of in planta expressed genes of Magnaporthe oryzae.

Publication Date: 2010 PMID: 20146797
Authors: Kim, S. - Park, J. - Park, S. Y. - Mitchell, T. K. - Lee, Y. H.
Journal: BMC Genomics

ABSTRACT: BACKGROUND: Infection of plants by pathogens and the subsequent disease development involves substantial changes in the biochemistry and physiology of both partners. Analysis of genes that are expressed during these interactions represents a powerful strategy to obtain insights into the molecular events underlying these changes. We have employed expressed sequence tag (EST) analysis to identify rice genes involved in defense responses against infection by the blast fungus Magnaporthe oryzae and fungal genes involved in infectious growth within the host during a compatible interaction. RESULTS: A cDNA library was constructed with RNA from rice leaves (Oryza sativa cv. Hwacheong) infected with M. oryzae strain KJ201. To enrich for fungal genes, subtraction library using PCR-based suppression subtractive hybridization was constructed with RNA from infected rice leaves as a tester and that from uninfected rice leaves as the driver. A total of 4,148 clones from two libraries were sequenced to generate 2,302 non-redundant ESTs. Of these, 712 and 1,562 ESTs could be identified to encode fungal and rice genes, respectively. To predict gene function, Gene Ontology (GO) analysis was applied, with 31% and 32% of rice and fungal ESTs being assigned to GO terms, respectively. One hundred uniESTs were found to be specific to fungal infection EST. More than 80 full-length fungal cDNA sequences were used to validate ab initio annotated gene model of M. oryzae genome sequence. CONCLUSION: This study shows the power of ESTs to refine genome annotation and functional characterization. Results of this work have advanced our understanding of the molecular mechanisms underpinning fungal-plant interactions and formed the basis for new hypothesis.

post to: CiteULike

Distinct, ecotype-specific genome and proteome signatures in the marine cyanobacteria Prochlorococcus.

Publication Date: 2010 Feb 10 PMID: 20146791
Authors: Paul, S. - Dutta, A. - Bag, S. K. - Das, S. - Dutta, C.
Journal: BMC Genomics

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 [almost equal to] 50% (group LLa), 4 LL strains having reduced genomes and G+C-content [almost equal to] 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.

post to: CiteULike