Comparative phenomics and targeted use of genomics reveals variation in carbon and nitrogen assimilation among different Brettanomyces bruxellensis strains.

Recent studies have suggested a correlation between genotype groups of Brettanomyces bruxellensis and their source of isolation. To further explore this relationship, the objective of this study was to assess metabolic differences in carbon and nitrogen assimilation between different B. bruxellensis strains from three beverages, including beer, wine, and soft drink, using Biolog Phenotype Microarrays. While some similarities of physiology were noted, many traits were variable among strains. Interestingly, some phenotypes were found that could be linked to strain origin, especially for the assimilation of particular α- and ß-glycosides as well as α- and ß-substituted monosaccharides. Based upon gene presence or absence, an α-glucosidase and ß-glucosidase were found explaining the observed phenotypes. Further, using a PCR screen on a large number of isolates, we have been able to specifically link a genomic deletion to the beer strains, suggesting that this region may have a fitness cost for B. bruxellensis in certain fermentation systems such as brewing. More specifically, none of the beer strains were found to contain a ß-glucosidase, which may have direct impacts on the ability for these strains to compete with other microbes or on flavor production.

Exoproteome analysis of Starmerella bombicola results in the discovery of an esterase required for lactonization of sophorolipids.

The yeast Starmerella bombicola secretes sophorolipids, a family of biosurfactants that find applications in green household products and cosmetics. Over the past years, a gene cluster was discovered that is responsible for the entire synthesis of the open (acidic) form of these molecules from glucose, fatty acids and acetyl-CoA building blocks. However, a significant fraction of the natural product is obtained as a ring closed form (lactonic). Both genetic and proteomic approaches hitherto failed to discover an enzyme responsible for the esterification reaction required for the ring closure step. We hypothesized that this enzyme is extracellularly secreted. Therefore, we characterized the composition of the S. bombicola exoproteome at different time points of the growth and compared it with known yeast exoproteomes. We identified 44 proteins, many of them commonly found in other fungi. Curiously, we discovered an enzyme with homology to Pseudozyma antarctica lipase A. A deletion mutation of its gene resulted in complete abolishment of the sophorolipid lactonization providing evidence that this might be the missing enzyme in the sophorolipid biosynthetic pathway. BIOLOGICAL SIGNIFICANCE: Growing concern about the impact of chemical processes on the environment increases consumers' demand for bio-based products. Lately, the household care and personal care sectors show increasing interest in naturally occurring biosurfactants, which constitute environment-friendly alternatives for chemical surfactants, typically derived from mineral oils. A particular group of biosurfactants, sophorolipids, already found their way to the market, being used in a range of household detergent products and in cosmetics. This work describes how proteomic approaches have led to the completion of our knowledge on the biosynthetic pathway of sophorolipids as performed by Starmerella bombicola, a fungus used in the industrial production of these biosurfactants. Moreover, we proved that by creating a deletion mutant in the lactone esterase discovered in this study, we can shape the biosynthesis towards custom-made sophorolipids with desired functions. Herewith, we demonstrate the potential of proteomics in industrial biotechnology.

Overview of tomato (Solanum lycopersicum) candidate pathogen recognition genes reveals important Solanum R locus dynamics.

To investigate the genome-wide spatial arrangement of R loci, a complete catalogue of tomato (Solanum lycopersicum) and potato (Solanum tuberosum) nucleotide-binding site (NBS) NBS, receptor-like protein (RLP) and receptor-like kinase (RLK) gene repertories was generated. Candidate pathogen recognition genes were characterized with respect to structural diversity, phylogenetic relationships and chromosomal distribution. NBS genes frequently occur in clusters of related gene copies that also include RLP or RLK genes. This scenario is compatible with the existence of selective pressures optimizing coordinated transcription. A number of duplication events associated with lineage-specific evolution were discovered. These findings suggest that different evolutionary mechanisms shaped pathogen recognition gene cluster architecture to expand and to modulate the defence repertoire. Analysis of pathogen recognition gene clusters associated with documented resistance function allowed the identification of adaptive divergence events and the reconstruction of the evolution history of these loci. Differences in candidate pathogen recognition gene number and organization were found between tomato and potato. Most candidate pathogen recognition gene orthologues were distributed at less than perfectly matching positions, suggesting an ongoing lineage-specific rearrangement. Indeed, a local expansion of Toll/Interleukin-1 receptor (TIR)-NBS-leucine-rich repeat (LRR) (TNL) genes in the potato genome was evident. Taken together, these findings have implications for improved understanding of the mechanisms of molecular adaptive selection at Solanum R loci.

Divergence of regulatory sequences in duplicated fish genes.

Duplicated genes can undergo different fates, from nonfunctionalization to subfunctionalization and neofunctionalization. In particular, changes in regulatory sequences affecting the expression domain of genes seem to be responsible for the latter two fates. In this study we used in silico motif detection to show how alterations in the composition of regulatory motifs between paralogous genes in zebrafish and Tetraodon might reflect the functional divergence of duplicates.

Deletions involving long-range conserved nongenic sequences upstream and downstream of FOXL2 as a novel disease-causing mechanism in blepharophimosis syndrome.

The expression of a gene requires not only a normal coding sequence but also intact regulatory regions, which can be located at large distances from the target genes, as demonstrated for an increasing number of developmental genes. In previous mutation studies of the role of FOXL2 in blepharophimosis syndrome (BPES), we identified intragenic mutations in 70% of our patients. Three translocation breakpoints upstream of FOXL2 in patients with BPES suggested a position effect. Here, we identified novel microdeletions outside of FOXL2 in cases of sporadic and familial BPES. Specifically, four rearrangements, with an overlap of 126 kb, are located 230 kb upstream of FOXL2, telomeric to the reported translocation breakpoints. Moreover, the shortest region of deletion overlap (SRO) contains several conserved nongenic sequences (CNGs) harboring putative transcription-factor binding sites and representing potential long-range cis-regulatory elements. Interestingly, the human region orthologous to the 12-kb sequence deleted in the polled intersex syndrome in goat, which is an animal model for BPES, is contained in this SRO, providing evidence of human-goat conservation of FOXL2 expression and of the mutational mechanism. Surprisingly, in a fifth family with BPES, one rearrangement was found downstream of FOXL2. In addition, we report nine novel rearrangements encompassing FOXL2 that range from partial gene deletions to submicroscopic deletions. Overall, genomic rearrangements encompassing or outside of FOXL2 account for 16% of all molecular defects found in our families with BPES. In summary, this is the first report of extragenic deletions in BPES, providing further evidence of potential long-range cis-regulatory elements regulating FOXL2 expression. It contributes to the enlarging group of developmental diseases caused by defective distant regulation of gene expression. Finally, we demonstrate that CNGs are candidate regions for genomic rearrangements in developmental genes.

The first green lineage cdc25 dual-specificity phosphatase.

The Cdc25 protein phosphatase is a key enzyme involved in the regulation of the G(2)/M transition in metazoans and yeast. However, no Cdc25 ortholog has so far been identified in plants, although functional studies have shown that an activating dephosphorylation of the CDK-cyclin complex regulates the G(2)/M transition. In this paper, the first green lineage Cdc25 ortholog is described in the unicellular alga Ostreococcus tauri. It encodes a protein which is able to rescue the yeast S. pombe cdc25-22 conditional mutant. Furthermore, microinjection of GST-tagged O. tauri Cdc25 specifically activates prophase-arrested starfish oocytes. In vitro histone H1 kinase assays and anti-phosphotyrosine Western Blotting confirmed the in vivo activating dephosphorylation of starfish CDK1-cyclinB by recombinant O. tauri Cdc25. We propose that there has been coevolution of the regulatory proteins involved in the control of M-phase entry in the metazoan, yeast and green lineages.

Analysis of the evolutionary relationships of HIV-1 and SIVcpz sequences using bayesian inference: implications for the origin of HIV-1.

The most plausible origin of HIV-1 group M is an SIV lineage currently represented by SIVcpz isolated from the chimpanzee subspecies Pan troglodytes troglodytes. The origin of HIV-1 group O is less clear. Putative recombination between any of the HIV-1 and SIVcpz sequences was tested using bootscanning and Bayesian-scanning plots, as well as a new method using a Bayesian multiple change-point (BMCP) model to infer parental sequences and crossing-over points. We found that in the case of highly divergent sequences, such as HIV-1/SIVcpz, Bayesian scanning and BMCP methods are more appropriate than bootscanning analysis to investigate spatial phylogenetic variation, including estimating the boundaries of the regions with discordant evolutionary relationships and the levels of support of the phylogenetic clusters under study. According to the Bayesian scanning plots and BMCP method, there was strong evidence for discordant phylogenetic clustering throughout the genome: (1) HIV-1 group O clustered with SIVcpzANT/TAN in middle pol, and partial vif/env; (2) SIVcpzGab1 clustered with SIVcpzANT/TAN in 3'pol/vif, and middle env; (3) HIV-1 group O grouped with SIVcpzCamUS and SIVcpzGab1 in p17/p24; (4) HIV-1 group M was more closely related to SIVcpzCamUS in 3'gag/pol and in middle pol, whereas in partial gp120 group M clustered with group O. Conditionally independent phylogenetic analysis inferred by maximum likelihood (ML) and Bayesian methods further confirmed these findings. The discordant phylogenetic relationships between the HIV-1/SIVcpz sequences may have been caused by ancient recombination events, but they are also due, at least in part, to altered rates of evolution between parental SIVcpz lineages.

The ghost of selection past: rates of evolution and functional divergence of anciently duplicated genes.

The duplication of genes and even complete genomes may be a prerequisite for major evolutionary transitions and the origin of evolutionary novelties. However, the evolutionary mechanisms of gene evolution and the origin of novel gene functions after gene duplication have been a subject of many debates. Recently, we compiled 26 groups of orthologous genes, which included one gene from human, mouse, and chicken, one or two genes from the tetraploid Xenopus and two genes from zebrafish. Comparative analysis and mapping data showed that these pairs of zebrafish genes were probably produced during a fish-specific genome duplication that occurred between 300 and 450 Mya, before the teleost radiation (Taylor et al. 2001). As discussed here, many of these retained duplicated genes code for DNA binding proteins. Different models have been developed to explain the retention of duplicated genes and in particular the subfunctionalization model of Force et al. (1999) could explain why so many developmental control genes have been retained. Other models are harder to reconcile with this particular set of duplicated genes. Most genes seem to have been subjected to strong purifying selection, keeping properties such as charge and polarity the same in both duplicates, although some evidence was found for positive Darwinian selection, in particular for Hox genes. However, since only the cumulative pattern of nucleotide substitutions can be studied, clear indications of positive Darwinian selection or neutrality may be hard to find for such anciently duplicated genes. Nevertheless, an increase in evolutionary rate in about half of the duplicated genes seems to suggest that either positive Darwinian selection has occurred or that functional constraints have been relaxed at one point in time during functional divergence.

Comparative genomics provides evidence for an ancient genome duplication event in fish.

There are approximately 25 000 species in the division Teleostei and most are believed to have arisen during a relatively short period of time ca. 200 Myr ago. The discovery of 'extra' Hox gene clusters in zebrafish (Danio rerio), medaka (Oryzias latipes), and pufferfish (Fugu rubripes), has led to the hypothesis that genome duplication provided the genetic raw material necessary for the teleost radiation. We identified 27 groups of orthologous genes which included one gene from man, mouse and chicken, one or two genes from tetraploid Xenopus and two genes from zebrafish. A genome duplication in the ancestor of teleost fishes is the most parsimonious explanation for the observations that for 15 of these genes, the two zebrafish orthologues are sister sequences in phylogenies that otherwise match the expected organismal tree, the zebrafish gene pairs appear to have been formed at approximately the same time, and are unlinked. Phylogenies of nine genes differ a little from the tree predicted by the fish-specific genome duplication hypothesis: one tree shows a sister sequence relationship for the zebrafish genes but differs slightly from the expected organismal tree and in eight trees, one zebrafish gene is the sister sequence to a clade which includes the second zebrafish gene and orthologues from Xenopus, chicken, mouse and man. For these nine gene trees, deviations from the predictions of the fish-specific genome duplication hypothesis are poorly supported. The two zebrafish orthologues for each of the three remaining genes are tightly linked and are, therefore, unlikely to have been formed during a genome duplication event. We estimated that the unlinked duplicated zebrafish genes are between 300 and 450 Myr. Thus, genome duplication could have provided the genetic raw material for teleost radiation. Alternatively, the loss of different duplicates in different populations (i.e. 'divergent resolution') may have promoted speciation in ancient teleost populations.

Phylogenetic relationships among algae based on complete large-subunit rRNA sequences.

The complete or nearly complete large-subunit rRNA (LSU rRNA) sequences were determined for representatives of several algal groups such as the chlorarachniophytes, cryptomonads, haptophytes, bacillariophytes, dictyochophytes and pelagophytes. Our aim was to study the phylogenetic position and relationships of the different groups of algae, and in particular to study the relationships among the different classes of heterokont algae. In LSU rRNA phylogenies, the chlorarachniophytes, cryptomonads and haptophytes seem to form independent evolutionary lineages, for which a specific relationship with any of the other eukaryotic taxa cannot be demonstrated. This is in accordance with phylogenies inferred on the basis of the small-subunit rRNA (SSU rRNA). Regarding the heterokont algae, which form a well-supported monophyletic lineage on the basis of LSU rRNA, resolution between the different classes could be improved by combining the SSU and LSU rRNA data. Based on a concatenated alignment of both molecules, the phaeophytes and the xanthophytes are sister taxa, as well as the pelagophytes and the dictyochophytes, and the chrysophytes and the eustigmatophytes. All these sister group relationships are highly supported by bootstrap analysis and by different methods of tree construction.

Genome duplication, divergent resolution and speciation.

What are the evolutionary consequences of gene duplication? One answer is speciation, according to a model initially called Reciprocal Silencing and recently expanded and renamed Divergent Resolution. This model shows how the loss of different copies of a duplicated gene in allopatric populations (divergent resolution) can promote speciation by genetically isolating these populations should they become reunited. Genome duplication events produce thousands of duplicated genes. Therefore, lineages with a history of genome duplication might have been especially prone to speciation via divergent resolution.

The European Large Subunit Ribosomal RNA Database.

The European Large Subunit Ribosomal RNA Database compiles all complete or nearly complete large subunit ribosomal RNA sequences available from public sequence databases. These are provided in aligned format and the secondary structure, as derived by comparative sequence analysis, is included. Additional information about the sequences such as literature references and taxonomic information is also included. The database is available from our WWW server at http://rrna.uia.ac.be/lsu/.

Distribution of substitution rates and location of insertion sites in the tertiary structure of ribosomal RNA.

The relative substitution rate of each nucleotide site in bacterial small subunit rRNA, large subunit rRNA and 5S rRNA was calculated from sequence alignments for each molecule. Two-dimensional and three-dimensional variability maps of the rRNAs were obtained by plotting the substitution rates on secondary structure models and on the tertiary structure of the rRNAs available from X-ray diffraction results. This showed that the substitution rates are generally low near the centre of the ribosome, where the nucleotides essential for its function are situated, and that they increase towards the surface. An inventory was made of insertions characteristic of the Archaea, Bacteria and Eucarya domains, and for additional insertions present in specific eukaryotic taxa. All these insertions occur at the ribosome surface. The taxon-specific insertions seem to arise randomly in the eukaryotic evolutionary tree, without any phylogenetic relatedness between the taxa possessing them.

An updated and comprehensive rRNA phylogeny of (crown) eukaryotes based on rate-calibrated evolutionary distances.

Recent experience with molecular phylogeny has shown that all molecular markers have strengths and weaknesses. Nonetheless, despite several notable discrepancies with phylogenies obtained from protein data, the merits of the small subunit ribosomal RNA (SSU rRNA) as a molecular phylogenetic marker remain indisputable. Over the last 10 to 15 years a massive SSU rRNA database has been gathered, including more then 3000 complete sequences from eukaryotes. This creates a huge computational challenge, which is exacerbated by phenomena such as extensive rate variation among sites in the molecule. A few years ago, a fast phylogenetic method was developed that takes into account among-site rate variation in the estimation of evolutionary distances. This "substitution rate calibration" (SRC) method not only corrects for a major source of artifacts in phylogeny reconstruction but, because it is based on a distance approach, allows comprehensive trees including thousands of sequences to be constructed in a reasonable amount of time. In this study, a nucleotide variability map and a phylogenetic tree were constructed, using the SRC method, based on all available (January 2000) complete SSU rRNA sequences (2551) for species belonging to the so-called eukaryotic crown. The resulting phylogeny constitutes the most complete description of overall eukaryote diversity and relationships to date. Furthermore, branch lengths estimated with the SRC method better reflect the huge differences in evolutionary rates among and within eukaryotic lineages. The ribosomal RNA tree is compared with a recent protein phylogeny obtained from concatenated actin, alpha-tubulin, beta-tubulin, and elongation factor 1-alpha amino acid sequences. A consensus phylogeny of the eukaryotic crown based on currently available molecular data is discussed, as well as specific problems encountered in analyzing sequences when large differences in substitution rate are present, either between different sequences (rate variation among lineages) or between different positions within the same sequence (among-site rate variation).

Comparative analysis of more than 3000 sequences reveals the existence of two pseudoknots in area V4 of eukaryotic small subunit ribosomal RNA.

The secondary structure of V4, the largest variable area of eukaryotic small subunit ribosomal RNA, was re-examined by comparative analysis of 3253 nucleotide sequences distributed over the animal, plant and fungal kingdoms and a diverse set of protist taxa. An extensive search for compensating base pair substitutions and for base covariation revealed that in most eukaryotes the secondary structure of the area consists of 11 helices and includes two pseudoknots. In one of the pseudoknots, exchange of base pairs between the two stems seems to occur, and covariation analysis points to the presence of a base triple. The area also contains three potential insertion points where additional hairpins or branched structures are present in a number of taxa scattered throughout the eukaryotic domain.