Algorithm to find gene expression profiles of deregulation and identify families of disease-altered genes

MOTIVATION: Alteration of gene expression often results in up- or down-regulated genes and the most common analysis strategies look for such differentially expressed genes. However, molecular disease mechanisms typically constitute abnormalities in the regulation of genes producing strong alterations in the expression levels. The search for such deregulation states in the genomic expression profiles will help to identify disease-altered genes better. RESULTS: We have developed an algorithm that searches for the genes which present a significant alteration in the variability of their expression profiles, by comparing an altered state with a control state. The algorithm provides groups of genes and assigns a statistical measure of significance to each group of genes selected. The method also includes a prefilter tool to select genes with a threshold of differential expression that can be set by the user ad casum. The method is evaluated using an experimental set of microarrays of human control and cancer samples from patients with acute promyelocytic leukemia.     

Identification of the human methylmalonyl-CoA racemase gene based on the analysis of prokaryotic gene arrangements. Implications for decoding the human genome

In this report, we identify the human DL-methylmalonyl-CoA racemase gene by analyzing prokaryotic gene arrangements and extrapolating the information obtained to human genes by homology searches. Sequence similarity searches were used to identify two groups of homologues that were frequently arranged with prokaryotic methylmalonyl-CoA mutase genes, and that were of unknown function. Both gene groups had homologues in the human genome. Because methylmalonyl-CoA mutases are involved in the metabolism of propionyl-CoA, we inferred that conserved neighbors of methylmalonyl-CoA mutase genes and their human homologues were also involved in this process. Subsequent biochemical studies confirmed this inference by showing that the prokaryotic gene PH0272 and its human homologue both encode DL-methylmalonyl-CoA racemases. To our knowledge this is the first report in which the function of a eukaryotic gene was determined based on the analysis of prokaryotic gene arrangements. Importantly, such analyses are rapid and may be generally applicable for the identification of human genes that lack homologues of known function or that have been misidentified on the basis of sequence similarity searches.     

Coalescent-based species tree inference from gene tree topologies under incomplete lineage sorting by maximum likelihood

Incomplete lineage sorting can cause incongruence between the phylogenetic history of genes (the gene tree) and that of the species (the species tree), which can complicate the inference of phylogenies. In this article, I present a new coalescent-based algorithm for species tree inference with maximum likelihood. I first describe an improved method for computing the probability of a gene tree topology given a species tree, which is much faster than an existing algorithm by Degnan and Salter (2005). Based on this method, I develop a practical algorithm that takes a set of gene tree topologies and infers species trees with maximum likelihood. This algorithm searches for the best species tree by starting from initial species trees and performing heuristic search to obtain better trees with higher likelihood. This algorithm, called STELLS (which stands for Species Tree InfErence with Likelihood for Lineage Sorting), has been implemented in a program that is downloadable from the author's web page. The simulation results show that the STELLS algorithm is more accurate than an existing maximum likelihood method for many datasets, especially when there is noise in gene trees. I also show that the STELLS algorithm is efficient and can be applied to real biological datasets.     

Trait-to-gene: a computational method for predicting the function of uncharacterized genes

The function of unknown genes is often inferred from comparisons to well-characterized homologs. In this paper, we show that, even if all of the homologs of a gene are unannotated, its function may be deduced through phylogenetic profiling. We have designed a series of algorithms that make functional predictions of genes based on orthology and set theory, but our approach to predicting gene function requires no previous knowledge of homolog function. With this technique, we successfully identified 94% of the clusters of orthologous groups that are known to be involved in flagella development or function. As a test, we removed the function of three putative flagellar genes that had been previously uncharacterized in Bacillus subtilis. We observed a motility phenotype for two of these three genes. Thus, these algorithms allow for high-throughput functional prediction of genes beyond that provided by simple orthology-based annotation endeavors.     

Resistance gene analogues of Arabidopsis thaliana: recognition by structure

Following completion of Arabidopsis thaliana sequencing projects, multiple resistance gene analogues (RGAs) have been identified. In this work a review of the current state of knowledge available in protein databases and scientific articles is presented. Putative resistance genes were identified by using BLAST searches as well as HMM fingerprints (the latter to infer existence of characteristic domains). The representation of all five classes of putative resistance genes in Col-0 ecotype was examined, along with the statistics on RGAs present on all five chromosomes of Arabidopsis thaliana.     

Condition-specific coregulation with cis-regulatory motifs and modules in the mouse genome

Deciphering genetic regulatory codes remains a challenge. Here, we present an effective approach to identifying in vivo condition-specific coregulation with cis-regulatory motifs and modules in the mouse genome. A resampling-based algorithm was adopted to cluster our microarray data of a stress response, which generated 35 tight clusters with unique expression patterns containing 811 genes of 5652 genes significantly altered. Database searches identified many known motifs within the 3-kb regulatory regions of 40 genes from 3 clusters and modules with six to nine motifs that were commonly shared by 60-100% of these genes. The upstream regulatory region contained the highest frequency of these common motifs. CisModule program predictions were comparable with the results from database searches and found four potentially novel motifs. This result indicates that these motifs and modules could be responsible for gene coregulation of the stress response in the lacrimal gland.     

Genomic exploration of the hemiascomycetous yeasts: 21. Comparative functional classification of genes

We explored the biological diversity of hemiascomycetous yeasts using a set of 22000 newly identified genes in 13 species through BLASTX searches. Genes without clear homologue in Saccharomyces cerevisiae appeared to be conserved in several species, suggesting that they were recently lost by S. cerevisiae. They often identified well-known species-specific traits. Cases of gene acquisition through horizontal transfer appeared to occur very rarely if at all. All identified genes were ascribed to functional classes. Functional classes were differently represented among species. Species classification by functional clustering roughly paralleled rDNA phylogeny. Unequal distribution of rapidly evolving, ascomycete-specific, genes among species and functions was shown to contribute strongly to this clustering. A few cases of gene family amplification were documented, but no general correlation could be observed between functional differentiation of yeast species and variations of gene family sizes. Yeast biological diversity seems thus to result from limited species-specific gene losses or duplications, and for a large part from rapid evolution of genes and regulatory factors dedicated to specific functions.     

Functional characterization of SR and SR-related genes in Caenorhabditis elegans

The SR proteins constitute a family of nuclear phosphoproteins, which are required for constitutive splicing and also influence alternative splicing regulation. Initially, it was suggested that SR proteins were functionally redundant in constitutive splicing. However, differences have been observed in alternative splicing regulation, suggesting unique functions for individual SR proteins. Homology searches of the Caenorhabditis elegans genome identified seven genes encoding putative orthologues of the human factors SF2/ASF, SRp20, SC35, SRp40, SRp75 and p54, and also several SR-related genes. To address the issue of functional redundancy, we used dsRNA interference (RNAi) to inhibit specific SR protein function during C.elegans development. RNAi with CeSF2/ASF caused late embryonic lethality, suggesting that this gene has an essential function during C.elegans development. RNAi with other SR genes resulted in no obvious phenotype, which is indicative of gene redundancy. Simultaneous interference of two or more SR proteins in certain combinations caused lethality or other developmental defects. RNAi with CeSRPK, an SR protein kinase, resulted in early embryonic lethality, suggesting an essential role for SR protein phosphorylation during development.     

Genome-wide analysis of a land plant-specific acyl:coenzyme A synthetase (ACS) gene family in Arabidopsis, poplar, rice and Physcomitrella

The plant enzyme 4-coumarate:coenzyme A ligase (4CL) is part of a family of adenylate-forming enzymes present in all organisms. Analysis of genome sequences shows the presence of '4CL-like' enzymes in plants and other organisms, but their evolutionary relationships and functions remain largely unknown. 4CL and 4CL-like genes were identified by BLAST searches in Arabidopsis, Populus, rice, Physcomitrella, Chlamydomonas and microbial genomes. Evolutionary relationships were inferred by phylogenetic analysis of aligned amino acid sequences. Expression patterns of a conserved set of Arabidopsis and poplar 4CL-like acyl-CoA synthetase (ACS) genes were assayed. The conserved ACS genes form a land plant-specific class. Angiosperm ACS genes grouped into five clades, each of which contained representatives in three fully sequenced genomes. Expression analysis revealed conserved developmental and stress-induced expression patterns of Arabidopsis and poplar genes in some clades. Evolution of plant ACS enzymes occurred early in land plants. Differential gene expansion of angiosperm ACS clades has occurred in some lineages. Evolutionary and gene expression data, combined with in vitro and limited in vivo protein function data, suggest that angiosperm ACS enzymes play conserved roles in octadecanoid and fatty acid metabolism, and play roles in organ development, for example in anthers.     

Gene and enhancer trap tagging of vascular-expressed genes in poplar trees

We report a gene discovery system for poplar trees based on gene and enhancer traps. Gene and enhancer trap vectors carrying the beta-glucuronidase (GUS) reporter gene were inserted into the poplar genome via Agrobacterium tumefaciens transformation, where they reveal the expression pattern of genes at or near the insertion sites. Because GUS expression phenotypes are dominant and are scored in primary transformants, this system does not require rounds of sexual recombination, a typical barrier to developmental genetic studies in trees. Gene and enhancer trap lines defining genes expressed during primary and secondary vascular development were identified and characterized. Collectively, the vascular gene expression patterns revealed that approximately 40% of genes expressed in leaves were expressed exclusively in the veins, indicating that a large set of genes is required for vascular development and function. Also, significant overlap was found between the sets of genes responsible for development and function of secondary vascular tissues of stems and primary vascular tissues in other organs of the plant, likely reflecting the common evolutionary origin of these tissues. Chromosomal DNA flanking insertion sites was amplified by thermal asymmetric interlaced PCR and sequenced and used to identify insertion sites by reference to the nascent Populus trichocarpa genome sequence. Extension of the system was demonstrated through isolation of full-length cDNAs for five genes of interest, including a new class of vascular-expressed gene tagged by enhancer trap line cET-1-pop1-145. Poplar gene and enhancer traps provide a new resource that allows plant biologists to directly reference the poplar genome sequence and identify novel genes of interest in forest biology.