REEF: searching REgionally Enriched Features in genomes

Maximum-scoring path sets on pangenome graphs of constant treewidth

We generalize a problem of finding maximum-scoring segment sets, previously studied by Csűrös (IEEE/ACM Transactions on Computational Biology and Bioinformatics, 2004, 1, 139-150), from sequences to graphs. Namely, given a vertex-weighted graph G and a non-negative startup penalty c, we can find a set of vertex-disjoint paths in G with maximum total score when each path's score is its vertices' total weight minus c. We call this new problem maximum-scoring path sets (MSPS). We present an algorithm that has a linear-time complexity for graphs with a constant treewidth. Generalization from sequences to graphs allows the algorithm to be used on pangenome graphs representing several related genomes and can be seen as a common abstraction for several biological problems on pangenomes, including searching for CpG islands, ChIP-seq data analysis, analysis of region enrichment for functional elements, or simple chaining problems.

Comparative analysis of the transcriptomes of two rice subspecies during domestication

Two subspecies of rice, Oryza sativa ssp. indica and O. sativa ssp. japonica, with reproductive isolation and differences in morphology and phenotypic differences, were established during the process of rice domestication. To understand how domestication has changed the transcriptomes of the two rice subspecies and given rise to the phenotypic differences, we obtained approximately 700 Gb RNA-Seq data from 26 indica and 25 japonica accessions, and identified 97,005 transcribed fragments and 4579 novel transcriptionally active regions. The two rice subspecies had significantly different gene expression profiles, we identified 1,357 (3.3% in all genes) differentially expressed genes (DEGs) between indica and japonica rice. Combining existing gene function studies, it is found that some of these differential genes are related to the differentiation of the two subspecies, such as grain shape and cold tolerance, etc. Functional annotation of these DEGs indicates that they are involved in cell wall biosynthesis and reproductive processes. Furthermore, compared with the non-DEGs, the DEGs from both subspecies had more 5'flanking regions with low polymorphism to divergence ratios, indicating a stronger positive selection pressure on the regulation of the DEGs. This study improves our understanding of the rice genome by comparatively analyzing the transcriptomes of indica and japonica rice and identifies DEGs those may be responsible for the reproductive isolation and phenotypic differences between the two rice subspecies.

Genome-wide analyses reveal the role of noncoding variation in complex traits during rice domestication

Genomes carry millions of noncoding variants, and identifying the tiny fraction with functional consequences is a major challenge for genomics. We assessed the role of selection on long noncoding RNAs (lncRNAs) for domestication-related changes in rice grains. Among 3363 lncRNA transcripts identified in early developing panicles, 95% of those with differential expression (329 lncRNAs) between Oryza sativa ssp. japonica and wild rice were significantly down-regulated in the domestication event. Joint genome and transcriptome analyses reveal that directional selection on lncRNAs altered the expression of energy metabolism genes during domestication. Transgenic experiments and population analyses with three focal lncRNAs illustrate that selection on these loci led to increased starch content and grain weight. Together, our findings indicate that genome-wide selection for lncRNA down-regulation was an important mechanism for the emergence of rice domestication traits.

The role of PKAc1 in gene regulation and trichodimerol production in Trichoderma reesei

Background

Trichoderma reesei represents a model system for investigation of plant cell wall degradation and its connection to light response. The cyclic adenosine monophosphate pathway (cAMP pathway) plays an important role in both physiological outputs, being crucial for regulation of photoreceptor function as well as for cellulase regulation on different carbon sources. Phosphorylation of photoreceptors and of the carbon catabolite repressor CRE1 was shown in ascomycetes, indicating a relevance of protein kinase A in regulation of the target genes of these transcription factors as well as an impact on regulation of induction specific genes. Moreover, the cAMP pathway impacts growth and development.

Results

Here, we investigated gene regulation by the catalytic subunit of protein kinase A (PKAc1) upon growth on cellulose. We found distinct gene sets for regulation upon growth in light and darkness with an overlap of only 13 genes. PKAc1 regulates metabolic genes as well as transport and defense functions. The overlap of gene regulation by PKAc1 with the genes representing the cAMP dependent regulatory output of the photoreceptor ENV1 indicates an involvement of PKA in this pathway, which counteracts its effects by contrasting regulation. Moreover, we found considerable overlap with the gene sets regulated under cellulase inducing conditions and by the carbon catabolite repressor CRE1. Our analysis also showed that PKAc1 regulates the genes of the SOR cluster associated with the biosynthesis of sorbicillinoids. The homologue of gin4, encoding a CAMK type kinase, which is regulated by PKAc1, CRE1 and YPR2 showed a moderate impact on trichodimerol production. We isolated trichodimerol as representative sorbicillin compound and established a method for its quantification in large sample sets using high performance thin layer chromatography (HPTLC), which can be broadly applied for secondary metabolite screening of mutants or different growth conditions. Due to the high expression levels of the SOR cluster under conditions of sexual development we crosschecked the relevance of PKAc1 under these conditions. We could show that PKAc1 impacts biosynthesis of trichodimerol in axenic growth and upon mating.

Conclusions

We conclude that PKAc1 is involved in light dependent regulation of plant cell wall degradation, including carbon catabolite repression as well as secondary metabolism and development in T. reesei.

Dataset of differential gene expression between total normal human thyroid and histologically normal thyroid adjacent to papillary thyroid carcinoma

This article contains further data and information from our published manuscript [1]. We aim to identify significant transcriptome alterations of total normal human thyroid vs. histologically normal thyroid adjacent to papillary thyroid carcinoma. We performed a systematic meta-analysis of all the available gene expression profiles for the whole organ also collecting gene expression data for the normal thyroid adjacent to papillary thyroid carcinoma. A differential quantitative transcriptome reference map was generated by using TRAM (Transcriptome Mapper) software able to combine, normalize and integrate a total of 35 datasets from total normal thyroid and 40 datasets from histologically normal thyroid adjacent to papillary thyroid carcinoma from different sources. This analysis identified genes and genome segments that significantly discriminated the two groups of samples. Differentially expressed genes were grouped and enrichment function analyses were performed identifying the main features of the differentially expressed genes between total normal thyroid and histologically normal thyroid adjacent to papillary thyroid carcinoma. The search for housekeeping genes retrieved 414 loci.

Characterization of Imprinted Genes in Rice Reveals Conservation of Regulation and Imprinting with Other Plant Species

Genomic imprinting is an epigenetic phenomenon by which certain genes display differential expression in a parent-of-origin-dependent manner. Hundreds of imprinted genes have been identified from several plant species. Here, we identified, with a high level of confidence, 208 imprinted gene candidates from rice (Oryza sativa). Imprinted genes of rice showed limited association with the transposable elements, which contrasts with findings from Arabidopsis (Arabidopsis thaliana). Generally, imprinting in rice is conserved within a species, but intraspecific variation also was detected. The imprinted rice genes do not show signatures of selection, which suggests that domestication has had a limited evolutionary consequence on genomic imprinting. Although conservation of imprinting in plants is limited, we show that some loci are imprinted in several different species. Moreover, our results suggest that different types of epigenetic regulation can be established either before or after fertilization. Imprinted 24-nucleotide small RNAs and their neighboring genes tend to express alleles from different parents. This association was not observed between 21-nucleotide small RNAs and their neighboring genes. Together, our findings suggest that the regulation of imprinting can be diverse, and genomic imprinting has evolutionary and biological significance.

A molecular view of the normal human thyroid structure and function reconstructed from its reference transcriptome map

BACKGROUND

The thyroid is the earliest endocrine structure to appear during human development, and thyroid hormones are necessary for proper organism development, in particular for the nervous system and heart, normal growth and skeletal maturation. To date a quantitative, validated transcriptional atlas of the whole normal human thyroid does not exist and the availability of a detailed expression map might be an excellent occasion to investigate the many features of the thyroid transcriptome.

RESULTS

We present a view at the molecular level of the normal human thyroid histology and physiology obtained by a systematic meta-analysis of all the available gene expression profiles for the whole organ. A quantitative transcriptome reference map was generated by using the TRAM (Transcriptome Mapper) software able to combine, normalize and integrate a total of 35 suitable datasets from different sources thus providing a typical reference expression value for each of the 27,275 known, mapped transcripts obtained. The experimental in vitro validation of data was performed by "Real-Time" reverse transcription polymerase chain reaction showing an excellent correlation coefficient (r = 0.93) with data obtained in silico.

CONCLUSIONS

Our study provides a quantitative global reference portrait of gene expression in the normal human thyroid and highlights differential expression between normal human thyroid and a pool of non-thyroid tissues useful for modeling correlations between thyroidal gene expression and specific thyroid functions and diseases. The experimental in vitro validation supports the possible usefulness of the human thyroid transcriptome map as a reference for molecular studies of the physiology and pathology of this organ.

A CRE1- regulated cluster is responsible for light dependent production of dihydrotrichotetronin in Trichoderma reesei

Changing light conditions, caused by the rotation of earth resulting in day and night or growth on the surface or within a substrate, result in considerably altered physiological processes in fungi. For the biotechnological workhorse Trichoderma reesei, regulation of glycoside hydrolase gene expression, especially cellulase expression was shown to be a target of light dependent gene regulation. Analysis of regulatory targets of the carbon catabolite repressor CRE1 under cellulase inducing conditions revealed a secondary metabolite cluster to be differentially regulated in light and darkness and by photoreceptors. We found that this cluster is involved in production of trichodimerol and that the two polyketide synthases of the cluster are essential for biosynthesis of dihydrotrichotetronine (syn. bislongiquinolide or bisorbibutenolide). Additionally, an indirect influence on production of the peptaibol antibiotic paracelsin was observed. The two polyketide synthetase genes as well as the monooxygenase gene of the cluster were found to be connected at the level of transcription in a positive feedback cycle in darkness, but negative feedback in light, indicating a cellular sensing and response mechanism for the products of these enzymes. The transcription factor TR_102497/YPR2 residing within the cluster regulates the cluster genes in a light dependent manner. Additionally, an interrelationship of this cluster with regulation of cellulase gene expression was detected. Hence the regulatory connection between primary and secondary metabolism appears more widespread than previously assumed, indicating a sophisticated distribution of resources either to degradation of substrate (feed) or to antagonism of competitors (fight), which is influenced by light.

Investigation of orthologous pathogen recognition gene-rich regions in solanaceous species

Pathogen receptor proteins such as receptor-like protein (RLP), receptor-like kinase (RLK), and nucleotide-binding leucine-rich repeat (NLR) play a leading role in plant immunity activation. The genome architecture of such genes has been extensively investigated in several plant species. However, we still know little about their elaborate reorganization that arose during the plant speciation process. Using recently released pepper and eggplant genome sequences, we were able to identify 1097 pathogen recognition genes (PRGs) in the cultivated pepper Zunla-1 and 775 in the eggplant line Nakate-Shinkuro. The retrieved genes were analysed for their tendency to cluster, using different methods to infer the means of grouping. Orthologous relationships among clustering loci were found, and interesting reshuffling within given loci was observed for each analysed species. The information obtained was integrated into a comparative map to highlight the evolutionary dynamics in which the PRG loci were involved. Diversification of 14 selected PRG-rich regions was also explored using a DNA target-enrichment approach. A large number of gene variants were found as well as rearrangements of sequences encoding single protein domain and changes in chromosome gene order among species. Gene duplication and transposition activity have clearly influenced plant genome R-gene architecture and diversification. Our findings contribute to addressing several biological questions concerning the parallel evolution that occurred between genomes of the family Solanaceae. Moreover, the integration of different methods proved a powerful approach to reconstruct the evolutionary history in plant families and to transfer important biology findings among plant genomes.

SegCorr a statistical procedure for the detection of genomic regions of correlated expression

BACKGROUND

Detecting local correlations in expression between neighboring genes along the genome has proved to be an effective strategy to identify possible causes of transcriptional deregulation in cancer. It has been successfully used to illustrate the role of mechanisms such as copy number variation (CNV) or epigenetic alterations as factors that may significantly alter expression in large chromosomal regions (gene silencing or gene activation).

RESULTS

The identification of correlated regions requires segmenting the gene expression correlation matrix into regions of homogeneously correlated genes and assessing whether the observed local correlation is significantly higher than the background chromosomal correlation. A unified statistical framework is proposed to achieve these two tasks, where optimal segmentation is efficiently performed using dynamic programming algorithm, and detection of highly correlated regions is then achieved using an exact test procedure. We also propose a simple and efficient procedure to correct the expression signal for mechanisms already known to impact expression correlation. The performance and robustness of the proposed procedure, called SegCorr, are evaluated on simulated data. The procedure is illustrated on cancer data, where the signal is corrected for correlations caused by copy number variation. It permitted the detection of regions with high correlations linked to epigenetic marks like DNA methylation.

CONCLUSIONS

SegCorr is a novel method that performs correlation matrix segmentation and applies a test procedure in order to detect highly correlated regions in gene expression.

Analysis of Light- and Carbon-Specific Transcriptomes Implicates a Class of G-Protein-Coupled Receptors in Cellulose Sensing

In fungi, most metabolic processes are subject to regulation by light. For Trichoderma reesei, light-dependent regulation of cellulase gene expression is specifically shown. Therefore, we intended to unravel the relationship between regulation of enzymes by the carbon source and regulation of enzymes by light. Our two-dimensional analysis included inducing and repressing carbon sources which we used to compare light-specific regulation to dark-specific regulation and to rule out effects specific for a single carbon source. We found close connections with respect to gene regulation as well as significant differences in dealing with carbon in the environment in light and darkness. Moreover, our analyses showed an intricate regulation mechanism for substrate degradation potentially involving surface sensing and provide a basis for knowledge-based screening for strain improvement.

In fungi, most metabolic processes are subject to regulation by light. Trichoderma reesei is adapted to degradation of plant cell walls and regulates production of the required enzymes in a manner dependent on the nutrient source and the light status. Here we investigated the interrelated relevance of two regulation levels of the transcriptome of T. reesei: light regulation and carbon source-dependent control. We show that the carbon source (cellulose, lactose, sophorose, glucose, or glycerol) is the major source of variation, with light having a modulating effect on transcript regulation. A total of 907 genes were regulated under cellulase-inducing conditions in light, and 947 genes were regulated in darkness, with 530 genes overlapping (1,324 in total). Only 218 of the 1,324 induction-specific genes were independent of light and not regulated by the BLR1, BLR2, and ENV1 photoreceptors. Analysis of the genomic distribution of genes regulated by light upon growth on cellulose revealed considerable overlap of light-regulated clusters with induction-specific clusters and carbohydrate-active enzyme (CAZyme) clusters. Further, we found evidence for the operation of a sensing mechanism for solid cellulosic substrates, with regulation of genes such as swo1, cip1, and cip2 or of genes encoding hydrophobins which is related to the cyclic AMP (cAMP)-dependent regulatory output of ENV1. We identified class XIII G-protein-coupled receptors (GPCRs) CSG1 and CSG2 in T. reesei as putative cellulose/glucose-sensing GPCRs. Our data indicate that the cellulase regulation pathway is bipartite, comprising a section corresponding to transcriptional regulation and one corresponding to posttranscriptional regulation, with the two connected by the function of CSG1.

IMPORTANCE In fungi, most metabolic processes are subject to regulation by light. For Trichoderma reesei, light-dependent regulation of cellulase gene expression is specifically shown. Therefore, we intended to unravel the relationship between regulation of enzymes by the carbon source and regulation of enzymes by light. Our two-dimensional analysis included inducing and repressing carbon sources which we used to compare light-specific regulation to dark-specific regulation and to rule out effects specific for a single carbon source. We found close connections with respect to gene regulation as well as significant differences in dealing with carbon in the environment in light and darkness. Moreover, our analyses showed an intricate regulation mechanism for substrate degradation potentially involving surface sensing and provide a basis for knowledge-based screening for strain improvement.

Study of Gene Trafficking between Acanthamoeba and Giant Viruses Suggests an Undiscovered Family of Amoeba-Infecting Viruses

The nucleocytoplasmic large DNA viruses (NCLDV) are a group of extremely complex double-stranded DNA viruses, which are major parasites of a variety of eukaryotes. Recent studies showed that certain unicellular eukaryotes contain fragments of NCLDV DNA integrated in their genome, when surprisingly many of these organisms were not previously shown to be infected by NCLDVs. These findings prompted us to search the genome of Acanthamoeba castellanii strain Neff (Neff), one of the most prolific hosts in the discovery of giant NCLDVs, for possible DNA inserts of viral origin. We report the identification of 267 markers of lateral gene transfer with viruses, approximately half of which are clustered in Neff genome regions of viral origins, transcriptionally inactive or exhibit nucleotide-composition signatures suggestive of a foreign origin. The integrated viral genes had diverse origin among relatives of viruses that infect Neff, including Mollivirus, Pandoravirus, Marseillevirus, Pithovirus, and Mimivirus However, phylogenetic analysis suggests the existence of a yet-undiscovered family of amoeba-infecting NCLDV in addition to the five already characterized. The active transcription of some apparently anciently integrated virus-like genes suggests that some viral genes might have been domesticated during the amoeba evolution. These insights confirm that genomic insertion of NCLDV DNA is a common theme in eukaryotes. This gene flow contributed fertilizing the eukaryotic gene repertoire and participated in the occurrence of orphan genes, a long standing issue in genomics. Search for viral inserts in eukaryotic genomes followed by environmental screening of the original viruses should be used to isolate radically new NCLDVs.

Transposable Elements versus the Fungal Genome: Impact on Whole-Genome Architecture and Transcriptional Profiles

Transposable elements (TEs) are exceptional contributors to eukaryotic genome diversity. Their ubiquitous presence impacts the genomes of nearly all species and mediates genome evolution by causing mutations and chromosomal rearrangements and by modulating gene expression. We performed an exhaustive analysis of the TE content in 18 fungal genomes, including strains of the same species and species of the same genera. Our results depicted a scenario of exceptional variability, with species having 0.02 to 29.8% of their genome consisting of transposable elements. A detailed analysis performed on two strains of Pleurotus ostreatus uncovered a genome that is populated mainly by Class I elements, especially LTR-retrotransposons amplified in recent bursts from 0 to 2 million years (My) ago. The preferential accumulation of TEs in clusters led to the presence of genomic regions that lacked intra- and inter-specific conservation. In addition, we investigated the effect of TE insertions on the expression of their nearby upstream and downstream genes. Our results showed that an important number of genes under TE influence are significantly repressed, with stronger repression when genes are localized within transposon clusters. Our transcriptional analysis performed in four additional fungal models revealed that this TE-mediated silencing was present only in species with active cytosine methylation machinery. We hypothesize that this phenomenon is related to epigenetic defense mechanisms that are aimed to suppress TE expression and control their proliferation.

Transposable elements (TEs) are enigmatic genetic units that have played important roles in the evolution of eukaryotic genomes. Since their discovery in the 1950s, they have gained increasing attention and are known today as active genome modelers in multiple species. Although these elements have been widely studied in plants, much less is known about their occurrence and impact on the fungal kingdom. Using a diverse set of basidiomycete and ascomycete fungi, we quantified and characterized a huge diversity of DNA and RNA transposable elements, and we identified species that had 0.02 to 29.8% of their genomes occupied by transposable elements. In addition, using our basidiomycete model Pleurotus ostreatus, we demonstrated how TE insertions produced detrimental effects on the expression of upstream and downstream genes, which were downregulated compared with the control groups. This silencing mechanism was present in the basidiomycetes tested but exhibited a patchy distribution in ascomycetes, and might be related to specific genome defense mechanisms that control transposon proliferation. This finding reveals the broader impact of transposable elements in fungi. In addition to their importance as long-term evolutionary forces, they play major roles in the more dynamic transcriptome regulation of certain species.

Widespread and Adaptive Alterations in Genome-Wide Gene Expression Associated with Ecological Divergence of Two Oryza Species

Ecological speciation is a common mechanism by which new species arise. Despite great efforts, the role of gene expression in ecological divergence and speciation is poorly understood. Here, we conducted a genome-wide gene expression investigation of two Oryza species that are evolutionarily young and distinct in ecology and morphology. Using digital gene expression technology and the paired-end RNA sequencing method, we obtained 21,415 expressed genes across three reproduction-related tissues. Of them, approximately 8% (1,717) differed significantly in expression levels between the two species and these differentially expressed genes are randomly distributed across the genome. Moreover, 62% (1,064) of the differentially expressed genes exhibited a signature of directional selection in at least one species. Importantly, the genes with differential expression between species evolved more rapidly at the 5' flanking sequences than the genes without differential expression relative to coding sequences, suggesting that cis-regulatory changes are likely adaptive and play an important role in the ecological divergence of the two species. Finally, we showed evidence of significant differentiation between species in phenotype traits and observed that genes with differential expression were overrepresented with functional terms involving phenotypic and ecological differentiation between the two species, including reproduction- and stress-related characteristics. Our findings demonstrate that ecological speciation is associated with widespread and adaptive alterations in genome-wide gene expression and provide new insights into the importance of regulatory evolution in ecological speciation in plants.

The transcriptome of lae1 mutants of Trichoderma reesei cultivated at constant growth rates reveals new targets of LAE1 function

Background

The putative methyltransferase LaeA is a global regulator that affects the expression of multiple secondary metabolite gene clusters in several fungi. In Trichoderma reesei, its ortholog LAE1 appears to predominantly regulate genes involved in increasing competitive fitness in its environment, including expression of cellulases and polysaccharide hydrolases. A drawback in all studies related to LaeA/LAE1 function so far, however, is that the respective loss-of-function and overexpressing mutants display different growth rates. Thus some of the properties attributed to LaeA/LAE1 could be simply due to changes of the growth rate.

Results

We cultivated T. reesei, a Δlae1 mutant and a lae1-overexpressing strain in chemostats on glucose at two different growth rates (0.075 and 0.020 h^-1) which resemble growth rates at repressing and derepressing conditions, respectively. Under these conditions, the effect of modulating LAE1 expression was mainly visible in the Δlae1 mutant, whereas the overexpressing strain showed little differences to the parent strain. The effect on the expression of some gene categories identified earlier (polyketide synthases, heterokaryon incompatibility proteins, PTH11-receptors) was confirmed, but in addition GCN5-N-acetyltransferases, amino acid permeases and flavin monooxygenases were identified as so far unknown major targets of LAE1 action. LAE1 was also shown to interfere with the regulation of expression of several genes by the growth rate. About a tenth of the genes differentially expressed in the Δlae1 mutant under either growth condition were found to be clustered in the genome, but no specific gene group was associated with this phenomenon.

Conclusions

Our data show that – using T. reesei LAE1 as a model - the investigation of transcriptome in regulatory mutants at constant growth rates leads to new insights into the physiological roles of the respective regulator.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-447) contains supplementary material, which is available to authorized users.

Targets of light signalling in Trichoderma reesei

Background

The tropical ascomycete Trichoderma reesei (Hypocrea jecorina) represents one of the most efficient plant cell wall degraders. Regulation of the enzymes required for this process is affected by nutritional signals as well as other environmental signals including light.

Results

Our transcriptome analysis of strains lacking the photoreceptors BLR1 and BLR2 as well as ENV1 revealed a considerable increase in the number of genes showing significantly different transcript levels in light and darkness compared to wild-type. We show that members of all glycoside hydrolase families can be subject to light dependent regulation, hence confirming nutrient utilization including plant cell wall degradation as a major output pathway of light signalling. In contrast to N. crassa, photoreceptor mediated regulation of carbon metabolism in T. reesei occurs primarily by BLR1 and BLR2 via their positive effect on induction of env1 transcription, rather than by a presumed negative effect of ENV1 on the function of the BLR complex. Nevertheless, genes consistently regulated by photoreceptors in N. crassa and T. reesei are significantly enriched in carbon metabolic functions. Hence, different regulatory mechanisms are operative in these two fungi, while the light dependent regulation of plant cell wall degradation appears to be conserved.

Analysis of growth on different carbon sources revealed that the oxidoreductive D-galactose and pentose catabolism is influenced by light and ENV1. Transcriptional regulation of the target enzymes in these pathways is enhanced by light and influenced by ENV1, BLR1 and/or BLR2. Additionally we detected an ENV1-regulated genomic cluster of 9 genes including the D-mannitol dehydrogenase gene lxr1, with two genes of this cluster showing consistent regulation in N. crassa.

Conclusions

We show that one major output pathway of light signalling in Trichoderma reesei is regulation of glycoside hydrolase genes and the degradation of hemicellulose building blocks. Targets of ENV1 and BLR1/BLR2 are for the most part distinct and indicate individual functions for ENV1 and the BLR complex besides their postulated regulatory interrelationship.

Systems biological approaches towards understanding cellulase production by Trichoderma reesei

Recent progress and improvement in "-omics" technologies has made it possible to study the physiology of organisms by integrated and genome-wide approaches. This bears the advantage that the global response, rather than isolated pathways and circuits within an organism, can be investigated ("systems biology"). The sequencing of the genome of Trichoderma reesei (teleomorph Hypocrea jecorina), a fungus that serves as a major producer of biomass-degrading enzymes for the use of renewable lignocellulosic material towards production of biofuels and biorefineries, has offered the possibility to study this organism and its enzyme production on a genome wide scale. In this review, I will highlight the use of genomics, transcriptomics, proteomics and metabolomics towards an improved and novel understanding of the biochemical processes that involve in the massive overproduction of secreted proteins.

Expression of biomass-degrading enzymes is a major event during conidium development in Trichoderma reesei

The conidium plays a critical role in the life cycle of many filamentous fungi, being the primary means for survival under unfavorable conditions. To investigate the transcriptional changes taking place during the transition from growing hyphae to conidia in Trichoderma reesei, microarray experiments were performed. A total of 900 distinct genes were classified as differentially expressed, relative to their expression at time zero of conidiation, at least at one of the time points analyzed. The main functional categories (FunCat) overrepresented among the upregulated genes were those involving solute transport, metabolism, transcriptional regulation, secondary metabolite synthesis, lipases, proteases, and, particularly, cellulases and hemicellulases. Categories overrepresented among the downregulated genes were especially those associated with ribosomal and mitochondrial functions. The upregulation of cellulase and hemicellulase genes was dependent on the function of the positive transcriptional regulator XYR1, but XYR1 exerted no influence on conidiation itself. At least 20% of the significantly regulated genes were nonrandomly distributed within the T. reesei genome, suggesting an epigenetic component in the regulation of conidiation. The significant upregulation of cellulases and hemicellulases during this process, and thus cellulase and hemicellulase content in the spores of T. reesei, contributes to the hypothesis that the ability to hydrolyze plant biomass is a major trait of this fungus enabling it to break dormancy and reinitiate vegetative growth after a period of facing unfavorable conditions.

TRAM (Transcriptome Mapper): database-driven creation and analysis of transcriptome maps from multiple sources

Background

Several tools have been developed to perform global gene expression profile data analysis, to search for specific chromosomal regions whose features meet defined criteria as well as to study neighbouring gene expression. However, most of these tools are tailored for a specific use in a particular context (e.g. they are species-specific, or limited to a particular data format) and they typically accept only gene lists as input.

Results

TRAM (Transcriptome Mapper) is a new general tool that allows the simple generation and analysis of quantitative transcriptome maps, starting from any source listing gene expression values for a given gene set (e.g. expression microarrays), implemented as a relational database. It includes a parser able to assign univocal and updated gene symbols to gene identifiers from different data sources. Moreover, TRAM is able to perform intra-sample and inter-sample data normalization, including an original variant of quantile normalization (scaled quantile), useful to normalize data from platforms with highly different numbers of investigated genes. When in 'Map' mode, the software generates a quantitative representation of the transcriptome of a sample (or of a pool of samples) and identifies if segments of defined lengths are over/under-expressed compared to the desired threshold. When in 'Cluster' mode, the software searches for a set of over/under-expressed consecutive genes. Statistical significance for all results is calculated with respect to genes localized on the same chromosome or to all genome genes. Transcriptome maps, showing differential expression between two sample groups, relative to two different biological conditions, may be easily generated. We present the results of a biological model test, based on a meta-analysis comparison between a sample pool of human CD34+ hematopoietic progenitor cells and a sample pool of megakaryocytic cells. Biologically relevant chromosomal segments and gene clusters with differential expression during the differentiation toward megakaryocyte were identified.

Conclusions

TRAM is designed to create, and statistically analyze, quantitative transcriptome maps, based on gene expression data from multiple sources. The release includes FileMaker Pro database management runtime application and it is freely available at http://apollo11.isto.unibo.it/software/, along with preconfigured implementations for mapping of human, mouse and zebrafish transcriptomes.

Differential gene expression in a rat model of depression based on persistent differences in exploratory activity

Affective disorders are often accompanied by changes in motivation and anxiety. We investigated the genome-wide gene expression patterns in an animal model of depression that separates Wistar rats belonging into clusters of persistently high anxiety/low motivation to explore and low anxiety/high motivation to explore (low explorers and high explorers, LE and HE, respectively), in three brain regions previously implicated in mood disorders (raphe, hippocampus and the frontal cortex). Several serotonin-, GABA-, and glutamatergic genes were differentially expressed in LE- and HE-rats. The analysis of Gene Ontology biological process terms associated with the differentially regulated genes identified a significant overrepresentation of genes involved in the neuron development, morphogenesis, and differentiation; the most enriched pathways from the Kyoto Encyclopedia of Genes and Genomes were the Wnt signalling, MAPK signalling, long-term potentiation, and long-term depression pathways. These findings corroborate some expression data from other models of depression, and suggest additional targets.

Liverbase: a comprehensive view of human liver biology

The Liverbase ( http://liverbase.hupo.org.cn ) integrates information on the human liver proteome, including the function, abundance, and subcellular localization of proteins as well as associated disease information. The overall objective of the Liverbase is to provide a unique public resource for the liver community by providing comprehensive functional annotation of proteins implicated in liver development and disease. The central database features are manually annotated proteins localized in or functionally associated with human liver. In this first version of Liverbase, the associated data includes the human liver proteome (6788 proteins) and transcriptome (11205 significantly expressed genes: 10224 from CHIP and 5422 from MPSS, respecively) from the Chinese human liver proteome project (CNHLPP). As a database made publicly available through the Web site, Liverbase provides browsing and searching capabilities and a compilation of external links to other databases and homepages. Liverbase enables (i) the establishment of liver GO slim with 51 nonredundant items; (ii) systematic searches for proteins within specific functional or metabolic pathways; (iii) systematic searches that aim to find the proteins that underlie common and rare liver diseases; and (iv) the integration of detailed protein annotations derived from the literature. Liverbase also contains an external links page with links to other biological databases and homepages, including GO, KEGG, pfam, SWISS-PROT, and GNF databases. Liverbase users can utilize all these information to conduct systems biology research on liver.

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

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.

AIRE activated tissue specific genes have histone modifications associated with inactive chromatin

The Autoimmune Regulator (AIRE) protein is expressed in thymic medullary epithelial cells, where it promotes the ectopic expression of tissue-restricted antigens needed for efficient negative selection of developing thymocytes. Mutations in AIRE cause APECED syndrome, which is characterized by a breakdown of self-tolerance. The molecular mechanism by which AIRE increases the expression of a variety of different genes remains unknown. Here, we studied AIRE-regulated genes using whole genome expression analysis and chromatin immunoprecipitation. We show that AIRE preferentially activates genes that are tissue-specific and characterized by low levels of initial expression in stably transfected HEK293 cell model and mouse thymic medullary epithelial cells. In addition, the AIRE-regulated genes lack active chromatin marks, such as histone H3 trimethylation (H3K4me3) and acetylation (AcH3), on their promoters. We also show that during activation by AIRE, the target genes acquire histone H3 modifications associated with transcription and RNA polymerase II. In conclusion, our data show that AIRE is able to promote ectopic gene expression from chromatin associated with histone modifications characteristic to inactive genes.

Filling gaps in PPAR-alpha signaling through comparative nutrigenomics analysis

Background

The application of high-throughput genomic tools in nutrition research is a widespread practice. However, it is becoming increasingly clear that the outcome of individual expression studies is insufficient for the comprehensive understanding of such a complex field. Currently, the availability of the large amounts of expression data in public repositories has opened up new challenges on microarray data analyses. We have focused on PPARα, a ligand-activated transcription factor functioning as fatty acid sensor controlling the gene expression regulation of a large set of genes in various metabolic organs such as liver, small intestine or heart. The function of PPARα is strictly connected to the function of its target genes and, although many of these have already been identified, major elements of its physiological function remain to be uncovered. To further investigate the function of PPARα, we have applied a cross-species meta-analysis approach to integrate sixteen microarray datasets studying high fat diet and PPARα signal perturbations in different organisms.

Results

We identified 164 genes (MDEGs) that were differentially expressed in a constant way in response to a high fat diet or to perturbations in PPARs signalling. In particular, we found five genes in yeast which were highly conserved and homologous of PPARα targets in mammals, potential candidates to be used as models for the equivalent mammalian genes. Moreover, a screening of the MDEGs for all known transcription factor binding sites and the comparison with a human genome-wide screening of Peroxisome Proliferating Response Elements (PPRE), enabled us to identify, 20 new potential candidate genes that show, both binding site, both change in expression in the condition studied. Lastly, we found a non random localization of the differentially expressed genes in the genome.

Conclusion

The results presented are potentially of great interest to resume the currently available expression data, exploiting the power of in silico analysis filtered by evolutionary conservation. The analysis enabled us to indicate potential gene candidates that could fill in the gaps with regards to the signalling of PPARα and, moreover, the non-random localization of the differentially expressed genes in the genome, suggest that epigenetic mechanisms are of importance in the regulation of the transcription operated by PPARα.

Motif discovery in promoters of genes co-localized and co-expressed during myeloid cells differentiation

Genes co-expressed may be under similar promoter-based and/or position-based regulation. Although data on expression, position and function of human genes are available, their true integration still represents a challenge for computational biology, hampering the identification of regulatory mechanisms. We carried out an integrative analysis of genomic position, functional annotation and promoters of genes expressed in myeloid cells. Promoter analysis was conducted by a novel multi-step method for discovering putative regulatory elements, i.e. over-represented motifs, in a selected set of promoters, as compared with a background model. The combination of transcriptional, structural and functional data allowed the identification of sets of promoters pertaining to groups of genes co-expressed and co-localized in regions of the human genome. The application of motif discovery to 26 groups of genes co-expressed in myeloid cells differentiation and co-localized in the genome showed that there are more over-represented motifs in promoters of co-expressed and co-localized genes than in promoters of simply co-expressed genes (CEG). Motifs, which are similar to the binding sequences of known transcription factors, non-uniformly distributed along promoter sequences and/or occurring in highly co-expressed subset of genes were identified. Co-expressed and co-localized gene sets were grouped in two co-expressed genomic meta-regions, putatively representing functional domains of a high-level expression regulation.

Conserved chromosomal clustering of genes governed by chromatin regulators in Drosophila

BACKGROUND

The trithorax group (trxG) and Polycomb group (PcG) proteins are responsible for the maintenance of stable transcriptional patterns of many developmental regulators. They bind to specific regions of DNA and direct the post-translational modifications of histones, playing a role in the dynamics of chromatin structure.

RESULTS

We have performed genome-wide expression studies of trx and ash2 mutants in Drosophila melanogaster. Using computational analysis of our microarray data, we have identified 25 clusters of genes potentially regulated by TRX. Most of these clusters consist of genes that encode structural proteins involved in cuticle formation. This organization appears to be a distinctive feature of the regulatory networks of TRX and other chromatin regulators, since we have observed the same arrangement in clusters after experiments performed with ASH2, as well as in experiments performed by others with NURF, dMyc, and ASH1. We have also found many of these clusters to be significantly conserved in D. simulans, D. yakuba, D. pseudoobscura and partially in Anopheles gambiae.

CONCLUSION

The analysis of genes governed by chromatin regulators has led to the identification of clusters of functionally related genes conserved in other insect species, suggesting this chromosomal organization is biologically important. Moreover, our results indicate that TRX and other chromatin regulators may act globally on chromatin domains that contain transcriptionally co-regulated genes.

Genomic expression during human myelopoiesis

Background

Human myelopoiesis is an exciting biological model for cellular differentiation since it represents a plastic process where multipotent stem cells gradually limit their differentiation potential, generating different precursor cells which finally evolve into distinct terminally differentiated cells. This study aimed at investigating the genomic expression during myeloid differentiation through a computational approach that integrates gene expression profiles with functional information and genome organization.

Results

Gene expression data from 24 experiments for 8 different cell types of the human myelopoietic lineage were used to generate an integrated myelopoiesis dataset of 9,425 genes, each reliably associated to a unique genomic position and chromosomal coordinate. Lists of genes constitutively expressed or silent during myelopoiesis and of genes differentially expressed in commitment phase of myelopoiesis were first identified using a classical data analysis procedure. Then, the genomic distribution of myelopoiesis genes was investigated integrating transcriptional and functional characteristics of genes. This approach allowed identifying specific chromosomal regions significantly highly or weakly expressed, and clusters of differentially expressed genes and of transcripts related to specific functional modules.

Conclusion

The analysis of genomic expression during human myelopoiesis using an integrative computational approach allowed discovering important relationships between genomic position, biological function and expression patterns and highlighting chromatin domains, including genes with coordinated expression and lineage-specific functions.