An ensemble biclustering approach for querying gene expression compendia with experimental lists

MOTIVATION: Query-based biclustering techniques allow interrogating a gene expression compendium with a given gene or gene list. They do so by searching for genes in the compendium that have a profile close to the average expression profile of the genes in this query-list. As it can often not be guaranteed that the genes in a long query-list will all be mutually coexpressed, it is advisable to use each gene separately as a query. This approach, however, leaves the user with a tedious post-processing of partially redundant biclustering results. The fact that for each query-gene multiple parameter settings need to be tested in order to detect the 'most optimal bicluster size' adds to the redundancy problem. RESULTS: To aid with this post-processing, we developed an ensemble approach to be used in combination with query-based biclustering. The method relies on a specifically designed consensus matrix in which the biclustering outcomes for multiple query-genes and for different possible parameter settings are merged in a statistically robust way. Clustering of this matrix results in distinct, non-redundant consensus biclusters that maximally reflect the information contained within the original query-based biclustering results. The usefulness of the developed approach is illustrated on a biological case study in Escherichia coli. Availability and implementation: Compiled Matlab code is available from http://homes.esat.kuleuven.be/~kmarchal/Supplementary_Information_DeSmet_2011/.     

An alignment-free approach for eukaryotic ITS2 annotation and phylogenetic inference

The ITS2 gene class shows a high sequence divergence among its members that have complicated its annotation and its use for reconstructing phylogenies at a higher taxonomical level (beyond species and genus). Several alignment strategies have been implemented to improve the ITS2 annotation quality and its use for phylogenetic inferences. Although, alignment based methods have been exploited to the top of its complexity to tackle both issues, no alignment-free approaches have been able to successfully address both topics. By contrast, the use of simple alignment-free classifiers, like the topological indices (TIs) containing information about the sequence and structure of ITS2, may reveal to be a useful approach for the gene prediction and for assessing the phylogenetic relationships of the ITS2 class in eukaryotes. Thus, we used the TI2BioP (Topological Indices to BioPolymers) methodology [1], [2], freely available at http://ti2biop.sourceforge.net/ to calculate two different TIs. One class was derived from the ITS2 artificial 2D structures generated from DNA strings and the other from the secondary structure inferred from RNA folding algorithms. Two alignment-free models based on Artificial Neural Networks were developed for the ITS2 class prediction using the two classes of TIs referred above. Both models showed similar performances on the training and the test sets reaching values above 95% in the overall classification. Due to the importance of the ITS2 region for fungi identification, a novel ITS2 genomic sequence was isolated from Petrakia sp. This sequence and the test set were used to comparatively evaluate the conventional classification models based on multiple sequence alignments like Hidden Markov based approaches, revealing the success of our models to identify novel ITS2 members. The isolated sequence was assessed using traditional and alignment-free based techniques applied to phylogenetic inference to complement the taxonomy of the Petrakia sp. fungal isolate.     

Profiling treatment-specific post-translational modifications in a complex proteome with subtractive substrate phage display

Proteolytic activation of zymogens or controlled degradation of inhibitory factors is part of a major regulatory system on the post-translational level to regulate treatment induced cellular stress responses. The identification of differential activity based substrates is thus of high interest to prioritize and validate candidate targets for drug discovery. Here we present a novel subtractive substrate phage display screening method for the selection of treatment induced post-translational peptide modifications in complex proteomes. We investigated this approach with tumor cells in response to a protease activating anticancer treatment modality using subtractive and iterative screening of cellular extracts derived from control and treated cells. Specific phage were identified that served as substrates for proteolytic activities in response to treatment related activity changes and could be distinguished from substrates for unspecific proteolytic background activities. Novel, selected peptide substrates were investigated in vitro and in vivo and showed high substrate specificity and functional biological significance.     

Sulfato/thiosulfato reducing bacteria characterization by FT-IR spectroscopy: a new approach to biocorrosion control

Sulfato and Thiosulfato Reducing Bacteria (SRB, TRB) can induce corrosion process on steel immersed in seawater. This phenomenon, called biocorrosion, costs approximatively 5 billion euros in France each year. We provide the first evidence that Fourier Transformed InfraRed (FTIR) spectroscopy is a competitive technique to evaluate the sulfurogen flora involved in biocorrosion in comparison with time consuming classical identification methods or PCR analyses. A great discrimination was obtained between SRB, TRB and some contamination bacteria known to be present in seawater and seem to be able to reduce sulfate under particular conditions. Moreover, this preliminary study demonstrates that FTIR spectroscopic and genotypic results present a good correlation (these results are confirmed by other data obtained before or later, data not shown here). The advantages gained by FTIR spectroscopy are to give information on strain phenotype and bacterial metabolism which are of great importance in corrosion processes.     

Deletion mutagenesis of p22phox subunit of flavocytochrome b558: identification of regions critical for gp91phox maturation and NADPH oxidase activity

The heterodimeric flavocytochrome b558, comprised of the two integral membrane proteins p22phox and gp91phox, mediates the transfer of electrons from NADPH to molecular oxygen in the phagocyte NADPH oxidase to generate the superoxide precursor of microbicidal oxidants. This study uses deletion mutagenesis to identify regions of p22phox required for maturation of gp91phox and for NADPH oxidase activity. N-terminal, C-terminal, or internal deletions of human p22phox were generated and expressed in Chinese hamster ovary cells with transgenes for gp91phox and two other NADPH oxidase subunits, p47phox, and p67phox. The results demonstrate that p22phox-dependent maturation of gp91phox carbohydrate, cell surface expression of gp91phox, and the enzymatic function of flavocytochrome b558 are closely correlated. Whereas the 5 N-terminal and 25 C-terminal amino acids are dispensable for these functions, the N-terminal 11 amino acids of p22phox are required, as is a hydrophilic region between amino acids 65 and 90. Upon deletion of 54 residues at the C terminus of p22phox (amino acids 142-195), maturation and cell surface expression of gp91phox was still preserved, although NADPH oxidase activity was absent, as expected, due to removal of a proline-rich domain between amino acids 151-160 that is required for recruitment of p47phox. Antibody binding studies indicate that the extreme N terminus of p22phox is inaccessible in the absence of cell permeabilization, supporting a model in which both the N- and C-terminal domains of p22phox extend into the cytoplasm, anchored by two membrane-embedded regions.     

Isolectins I-A and I-B of Griffonia (Bandeiraea) simplicifolia. Crystal structure of metal-free GS I-B(4) and molecular basis for metal binding and monosaccharide specificity.

Seeds from the African legume shrub Griffonia simplicifolia contain several lectins. Among them the tetrameric lectin GS I-B(4) has strict specificity for terminal alpha Gal residues, whereas the closely related lectin GS I-A(4) can also bind to alpha GalNAc. These two lectins are commonly used as markers in histology or for research in xenotransplantation. To elucidate the basis for the fine difference in specificity, the amino acid sequences of both lectins have been determined and show 89% identity. The crystal structure of GS I-B(4), determined at 2.5-A resolution, reveals a new quaternary structure that has never been observed in other legume lectins. An unexpected loss of both Ca(2+) and Mn(2+) ions, which are necessary for carbohydrate binding in legume lectins, may be related to a particular amino acid sequence Pro-Glu-Pro in the metal binding loop. Comparison with demetallized concanavalin A reveals a different process for the loss of metal ions and for the subsequent loss of carbohydrate binding activity. The GS I-A x alpha GalNAc and GS I-B x alpha Gal complexes were constructed using homology modeling and docking approaches. The unusual presence of an aromatic amino acid at position 47 (Tyr in I-A and Trp in I-B) explains the strong preference for alpha-anomeric sugars in both isolectins. Alteration at one amino acid position, Ala(106) in I-A versus Glu(106) in I-B, is the basis for the observed specificities toward alpha GalNAc and alpha Gal.     

Disuniting uniformity: a pied cladistic canvas of mtDNA haplogroup H in Eurasia

It has been often stated that the overall pattern of human maternal lineages in Europe is largely uniform. Yet this uniformity may also result from an insufficient depth and width of the phylogenetic analysis, in particular of the predominant western Eurasian haplogroup (Hg) H that comprises nearly a half of the European mitochondrial DNA (mtDNA) pool. Making use of the coding sequence information from 267 mtDNA Hg H sequences, we have analyzed 830 mtDNA genomes, from 11 European, Near and Middle Eastern, Central Asian, and Altaian populations. In addition to the seven previously specified subhaplogroups, we define fifteen novel subclades of Hg H present in the extant human populations of western Eurasia. The refinement of the phylogenetic resolution has allowed us to resolve a large number of homoplasies in phylogenetic trees of Hg H based on the first hypervariable segment (HVS-I) of mtDNA. As many as 50 out of 125 polymorphic positions in HVS-I were found to be mutated in more than one subcluster of Hg H. The phylogeographic analysis revealed that sub-Hgs H1*, H1b, H1f, H2a, H3, H6a, H6b, and H8 demonstrate distinct phylogeographic patterns. The monophyletic subhaplogroups of Hg H provide means for further progress in the understanding of the (pre)historic movements of women in Eurasia and for the understanding of the present-day genetic diversity of western Eurasians in general.     

Distinct unfolding and refolding pathways of ribonuclease a revealed by heating and cooling temperature jumps

Heating and cooling temperature jumps (T-jumps) were performed using a newly developed technique to trigger unfolding and refolding of wild-type ribonuclease A and a tryptophan-containing variant (Y115W). From the linear Arrhenius plots of the microscopic folding and unfolding rate constants, activation enthalpy (ΔH^#), and activation entropy (ΔS^#) were determined to characterize the kinetic transition states (TS) for the unfolding and refolding reactions. The single TS of the wild-type protein was split into three for the Y115W variant. Two of these transition states, TS1 and TS2, characterize a slow kinetic phase, and one, TS3, a fast phase. Heating T-jumps induced protein unfolding via TS2 and TS3; cooling T-jumps induced refolding via TS1 and TS3. The observed speed of the fast phase increased at lower temperature, due to a strongly negative ΔH^# of the folding-rate constant. The results are consistent with a path-dependent protein folding/unfolding mechanism. TS1 and TS2 are likely to reflect X-Pro¹¹⁴ isomerization in the folded and unfolded protein, respectively, and TS3 the local conformational change of the β-hairpin comprising Trp¹¹⁵. A very fast protein folding/unfolding phase appears to precede both processes. The path dependence of the observed kinetics is suggestive of a rugged energy protein folding funnel.