Predicting functional gene links from phylogenetic-statistical analyses of whole genomes

An important element of the developing field of proteomics is to understand protein-protein interactions and other functional links amongst genes. Across-species correlation methods for detecting functional links work on the premise that functionally linked proteins will tend to show a common pattern of presence and absence across a range of genomes. We describe a maximum likelihood statistical model for predicting functional gene linkages. The method detects independent instances of the correlated gain or loss of pairs of proteins on phylogenetic trees, reducing the high rates of false positives observed in conventional across-species methods that do not explicitly incorporate a phylogeny. We show, in a dataset of 10,551 protein pairs, that the phylogenetic method improves by up to 35% on across-species analyses at identifying known functionally linked proteins. The method shows that protein pairs with at least two to three correlated events of gain or loss are almost certainly functionally linked. Contingent evolution, in which one gene's presence or absence depends upon the presence of another, can also be detected phylogenetically, and may identify genes whose functional significance depends upon its interaction with other genes. Incorporating phylogenetic information improves the prediction of functional linkages. The improvement derives from having a lower rate of false positives and from detecting trends that across-species analyses miss. Phylogenetic methods can easily be incorporated into the screening of large-scale bioinformatics datasets to identify sets of protein links and to characterise gene networks.

Report of a chimeric origin of transposable elements in a bovine-coding gene

Despite the wide distribution of transposable elements (TEs) in mammalian genomes, part of their evolutionary significance remains to be discovered. Today there is a substantial amount of evidence showing that TEs are involved in the generation of new exons in different species. In the present study, we searched 22,805 genes and reported the occurrence of TE-cassettes in coding sequences of 542 cow genes using the RepeatMasker program. Despite the significant number (542) of genes with TE insertions in exons only 14 (2.6%) of them were translated into protein, which we characterized as chimeric genes. From these chimeric genes, only the FAST kinase domains 3 (FASTKD3) gene, present on chromosome BTA 20, is a functional gene and showed evidence of the exaptation event. The genome sequence analysis showed that the last exon coding sequence of bovine FASTKD3 is similar to 85% similar to the ART2A retrotransposon sequence. In addition, comparison among FASTKD3 proteins shows that the last exon is very divergent from those of Homo sapiens, Pan troglodytes and Canis familiares. We suggest that the gene structure of bovine FASTKD3 gene could have originated by several ectopic recombinations between TE copies. Additionally, the absence of TE sequences in all other species analyzed suggests that the TE insertion is clade-specific, mainly in the ruminant lineage.

Identification and characterisation of a cytochrome P450 gene and processed pseudogene from an arachnid: the cattle tick, Boophilus microplus

We isolated and sequenced the first known cytochrome P450 gene and pseudogene from an arachnid, the cattle tick, Boophilus microplus. Bath the gene and pseudogene belong to the family CYP4, but a new subfamily, CYP4W, had to be created for these genes because they are substantially different to other CYP4 genes. The gene, CPP4W1, has greatest homology with CYP4C1 from a cockroach, Blaberus discoidalis. The predicted molecular weight of the protein encoded by CYP4W1 (63 KDa) is greater than that of the other CYP4 genes. The pseudogene, CYP4W1P, is probably a processed pseudogene derived from the functional gene CYP4W1. This is only the third CYP processed pseudogene to be identified. The pseudogene is 98% identical to the functional gene, CYP4W1, therefore we hypothesise that this pseudogene evolved recently from the functional gene. The CYP4 genes from arthropods have diverged from each other more than those of mammals; consequently the phylogeny of the arthropod genes could not be resolved. (C) 1999 Elsevier Science Ltd. All rights reserved.

The murine chaperonin 10 gene family contains an intronless, putative gene for early pregnancy factor, Cpn10-rs1

Early pregnancy factor (EPF) is a secreted protein with growth regulatory and immunomodulatory properties. Human platelet-derived EPF shares amino acid sequence identity with chaperonin 10 (Cpn10), a mitochondrial matrix protein which functions as a molecular chaperone. The striking differences in cellular localization and function of the two proteins suggest differential regulation of production reflecting either alternative transcription of the same gene or transcription from different genes. In mammals and more distantly related genera, there is a large gene family with homology to CPN 10 cDNA, which includes intronless copies of the coding sequence. To determine whether this could represent the gene for EPF, we have screened a mouse genomic library and sequenced representative Cpn10 family members, looking for a functional gene distinct from that of Cpn 10, which could encode EPF. Eight distinct genes were identified. Cpn10 contains introns, while other members are intronless. Six of these appear to be pseudogenes, and the remaining member, Cpn10-rs1, would encode a full-length protein. The 309-bp open reading frame (ORF) is identical to that of mouse Cpn10 cDNA with the exception of three single-base changes, two resulting in amino acid changes. Only one further single nucleotide difference between the Cpn10-rs1 and Cpn10 cDNAs is observed, located in the 3' UTR. Single nucleotide primer extension was applied to discriminate between Cpn10-rs1 and Cpn10 expression. Cpn10, which is ubiquitous, was detected in all tissue samples tested, whereas Cpn10-rs1 was expressed selectively. The pattern was completely coincident with known patterns of EPF activity, strongly suggesting that Cpn10-rs1 does encode EPF. The complete ORF of Cpn10-rs1 was expressed in E. coli. The purified recombinant protein was found to be equipotent with native human platelet-derived EPF in the bioassay for EPF, the rosette inhibition test.

Mapping of a Leishmania major gene/locus that confers pentamidine resistance by deletion and insertion of transposable element

Pentamidine (PEN) is an alternative compound to treat antimony-resistant leishmaniasis patients, which cellular target remains unclear. One approach to the identification of prospective targets is to identify genes able to mediate PEN resistance following overexpression. Starting from a genomic library of transfected parasites bearing a multicopy episomal cosmid vector containing wild-type Leishmania major DNA, we isolated one locus capable to render PEN resistance to wild type cells after DNA transfection. In order to map this Leishmania locus, cosmid insert was deleted by two successive sets of partial digestion with restriction enzymes, followed by transfection into wild type cells, overexpression, induction and functional tests in the presence of PEN. To determine the Leishmania gene related to PEN resistance, nucleotide sequencing experiments were done through insertion of the transposon Mariner element of Drosophila melanogaster (mosK) into the deleted insert to work as primer island. Using general molecular techniques, we described here this method that permits a quickly identification of a functional gene facilitating nucleotide sequence experiments from large DNA fragments. Followed experiments revealed the presence of a P-Glycoprotein gene in this locus which role in Leishmania metabolism has now been analyzed.

A comprehensive analysis of gene expression profiles in distal parts of the mouse renal tubule.

The distal parts of the renal tubule play a critical role in maintaining homeostasis of extracellular fluids. In this review, we present an in-depth analysis of microarray-based gene expression profiles available for microdissected mouse distal nephron segments, i.e., the distal convoluted tubule (DCT) and the connecting tubule (CNT), and for the cortical portion of the collecting duct (CCD; Zuber et al., Proc Natl Acad Sci USA 106:16523-16528, 2009). Classification of expressed transcripts in 14 major functional gene categories demonstrated that all principal proteins involved in maintaining the salt and water balance are represented by highly abundant transcripts. However, a significant number of transcripts belonging, for instance, to categories of G-protein-coupled receptors or serine/threonine kinases exhibit high expression levels but remain unassigned to a specific renal function. We also established a list of genes differentially expressed between the DCT/CNT and the CCD. This list is enriched by genes related to segment-specific transport functions and by transcription factors directing the development of the distal nephron or collecting ducts. Collectively, this in silico analysis provides comprehensive information about relative abundance and tissue specificity of the DCT/CNT and the CCD expressed transcripts and identifies new candidate genes for renal homeostasis.

The role of gene duplication in the origin and evolution of new biological functions

Abstract : Gene duplication is an essential source of material for the origin of genetic novelty and the evolution of lineage- or species-specific phenotypic traits. The reverse transcription of source gene mRNA followed by the genomic insertion of the resulting cDNA - retroposition - has provided the human genome with a significant number of gene copies during the last ~63 million years (MYA) of primate evolution. We estimated that at least 1 new functional gene (retrogene) per MYA emerged by retroposition in the primate lineage leading to humans. Using a combination of comparative sequencing and evolutionary simulations, we obtained strong evidence of functionality for 7 primate specific retrogenes. Most of these genes are specifically expressed in testis suggesting that retroposition has contributed with genetic raw material necessary for the evolution ofmale-specific functions in primates. We characterized CDC14Bretro (identified in the previous survey) that originated from the retroposition of a cell cycle gene - CDC14B - in the common ancestor of humans and apes. We demonstrate that CDC14Bretro experienced a period of intense positive selection in the African ape ancestor. By virtue of the amino acid substitutions that occurred during this period CDC 14Bretro adapted to a new subcellular compartment in African apes. Further analyses indicate that this subcellular shift reflects the evolution of anew functional role of CDC 14Bretro. Prompted by this result, we used yeast (Saccharomyces cerevisiae) to investigate on a global scale the extent of functional diversification of duplicate genes through the subcellular adaptation of their encoded proteins. We found that duplicate proteins frequently evolved new cellular localization patterns, either by partitioning of ancestral localizations ("sublocalization"), or more frequently by relocalization to previously unoccupied compartments ("neolocalization"). Interestingly, proteins involved in processes with a wider subcellular distribution more frequently evolved new localization patterns suggesting that subcellular localization changes are dependent on progenitor gene functions. Relocated proteins adapted to their new subcellular environments and evolved new functional roles through changes of their physio-chemical properties, expression levels, and interaction partners. Our work suggests an important role of subcellular adaptation for the emergence of new gene functions.

Functional diversity of bacterial genes associated with aromatic hydrocarbon degradation in anthropogenic dark earth of Amazonia

The objective of this work was to evaluate the catabolic gene diversity for the bacterial degradation of aromatic hydrocarbons in anthropogenic dark earth of Amazonia (ADE) and their biochar (BC). Functional diversity analyses in ADE soils can provide information on how adaptive microorganisms may influence the fertility of soils and what is their involvement in biogeochemical cycles. For this, clone libraries containing the gene encoding for the alpha subunit of aromatic ring-hydroxylating dioxygenases (α-ARHD bacterial gene) were constructed, totaling 800 clones. These libraries were prepared from samples of an ADE soil under two different land uses, located at the Caldeirão Experimental Station - secondary forest (SF) and agriculture (AG) -, and the biochar (SF_BC and AG_BC, respectively). Heterogeneity estimates indicated greater diversity in BC libraries; and Venn diagrams showed more unique operational protein clusters (OPC) in the SF_BC library than the ADE soil, which indicates that specific metabolic processes may occur in biochar. Phylogenetic analysis showed unidentified dioxygenases in ADE soils. Libraries containing functional gene encoding for the alpha subunit of the aromatic ring-hydroxylating dioxygenases (ARHD) gene from biochar show higher diversity indices than those of ADE under secondary forest and agriculture.

RNA interference is ineffective as a routine method for gene silencing in chick embryos as monitored by fgf8 silencing.

The in vivo accessibility of the chick embryo makes it a favoured model system for experimental developmental biology. Although the range of available techniques now extends to miss-expression of genes through in ovo electroporation, it remains difficult to knock out individual gene expression. Recently, the possibility of silencing gene expression by RNAi in chick embryos has been reported. However, published studies show only discrete quantitative differences in the expression of the endogenous targeted genes and unclear morphological alterations. To elucidate whether the tools currently available are adequate to silence gene expression sufficiently to produce a clear and specific null-like mutant phenotype, we have performed several experiments with different molecules that trigger RNAi: dsRNA, siRNA, and shRNA produced from a plasmid coexpressing green fluorescent protein as an internal marker. Focussing on fgf8 expression in the developing isthmus, we show that no morphological defects are observed, and that fgf8 expression is neither silenced in embryos microinjected with dsRNA nor in embryos microinjected and electroporated with a pool of siRNAs. Moreover, fgf8 expression was not significantly silenced in most isthmic cells transformed with a plasmid producing engineered shRNAs to fgf8. We also show that siRNA molecules do not spread significantly from cell to cell as reported for invertebrates, suggesting the existence of molecular differences between different model systems that may explain the different responses to RNAi. Although our results are basically in agreement with previously reported studies, we suggest, in contrast to them, that with currently available tools and techniques the number of cells in which fgf8 gene expression is decreased, if any, is not sufficient to generate a detectable mutant phenotype, thus making RNAi useless as a routine method for functional gene analysis in chick embryos.

Constrained models of evolution lead to improved prediction of functional linkage from correlated gain and loss of genes

Motivation: We compare phylogenetic approaches for inferring functional gene links. The approaches detect independent instances of the correlated gain and loss of pairs of genes from species' genomes. We investigate the effect on results of basing evidence of correlations on two phylogenetic approaches, Dollo parsminony and maximum likelihood (ML). We further examine the effect of constraining the ML model by fixing the rate of gene gain at a low value, rather than estimating it from the data. Results: We detect correlated evolution among a test set of pairs of yeast (Saccharomyces cerevisiae) genes, with a case study of 21 eukaryotic genomes and test data derived from known yeast protein complexes. If the rate at which genes are gained is constrained to be low, ML achieves by far the best results at detecting known functional links. The model then has fewer parameters but it is more realistic by preventing genes from being gained more than once. Availability: BayesTraits by M. Pagel and A. Meade, and a script to configure and repeatedly launch it by D. Barker and M. Pagel, are available at http://www.evolution.reading.ac.uk .

Characterization of para-Nitrophenol-degrading bacterial communities in river water by using functional markers and stable isotope probing

Microbial degradation is a major determinant of the fate of pollutants in the environment. para-Nitrophenol (PNP) is an EPA listed priority pollutant with a wide environmental distribution, but little is known about the microorganisms that degrade it in the environment. We studied the diversity of active PNP-degrading bacterial populations in river water using a novel functional marker approach coupled with [13C6]PNP stable isotope probing (SIP). Culturing together with culture-independent terminal restriction fragment length polymorphism analysis of 16S rRNA gene amplicons identified Pseudomonas syringae to be the major driver of PNP degradation in river water microcosms. This was confirmed by SIP-pyrosequencing of amplified 16S rRNA. Similarly, functional gene analysis showed that degradation followed the Gram-negative bacterial pathway and involved pnpA from Pseudomonas spp. However, analysis of maleylacetate reductase (encoded by mar), an enzyme common to late stages of both Gram-negative and Gram-positive bacterial PNP degradation pathways, identified a diverse assemblage of bacteria associated with PNP degradation, suggesting that mar has limited use as a specific marker of PNP biodegradation. Both the pnpA and mar genes were detected in a PNP-degrading isolate, P. syringae AKHD2, which was isolated from river water. Our results suggest that PNP-degrading cultures of Pseudomonas spp. are representative of environmental PNP-degrading populations.

Targeted Nucleotide Exchange in Bovine Myostatin Gene

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Report of a chimeric origin of transposable elements in a bovine-coding gene

Despite the wide distribution of transposable elements (TEs) in mammalian genomes, part of their evolutionary significance remains to be discovered. Today there is a substantial amount of evidence showing that TEs are involved in the generation of new exons in different species. In the present study, we searched 22,805 genes and reported the occurrence of TE-cassettes in coding sequences of 542 cow genes using the RepeatMasker program. Despite the significant number (542) of genes with TE insertions in exons only 14 (2.6%) of them were translated into protein, which we characterized as chimeric genes. From these chimeric genes, only the FAST kinase domains 3 (FASTKD3) gene, present on chromosome BTA 20, is a functional gene and showed evidence of the exaptation event. The genome sequence analysis showed that the last exon coding sequence of bovine FASTKD3 is ∼85% similar to the ART2A retrotransposon sequence. In addition, comparison among FASTKD3 proteins shows that the last exon is very divergent from those of Homo sapiens, Pan troglodytes and Canis familiares. We suggest that the gene structure of bovine FASTKD3 gene could have originated by several ectopic recombinations between TE copies. Additionally, the absence of TE sequences in all other species analyzed suggests that the TE insertion is clade-specific, mainly in the ruminant lineage. ©FUNPEC-RP.

Functional diversity of bacterial genes associated with aromatic hydrocarbon degradation in anthropogenic dark earth of Amazonia

The objective of this work was to evaluate the catabolic gene diversity for the bacterial degradation of aromatic hydrocarbons in anthropogenic dark earth of Amazonia (ADE) and their biochar (BC). Functional diversity analyses in ADE soils can provide information on how adaptive microorganisms may influence the fertility of soils and what is their involvement in biogeochemical cycles. For this, clone libraries containing the gene encoding for the alpha subunit of aromatic ring-hydroxylating dioxygenases (alpha-A RH D bacterial gene) were constructed, totaling 800 clones. These libraries were prepared from samples of an ADE soil under two different land uses, located at the Caldeirao Experimental Station secondary forest (SF) and agriculture (AG)-, and the biochar (SF_BC and AG_BC, respectively). Heterogeneity estimates indicated greater diversity in BC libraries; and Venn diagrams showed more unique operational protein clusters (OPC) in the SF_BC library than the ADE soil, which indicates that specific metabolic processes may occur in biochar. Phylogenetic analysis showed unidentified dioxygenases in ADE soils. Libraries containing functional gene encoding for the alpha subunit of the aromatic ring-hydroxylating dioxygenases (ARHD) gene from biochar show higher diversity indices than those of ADE under secondary forest and agriculture.

Isolation of an active gene and of two pseudogenes for mouse U7 small nuclear RNA

Three U7 RNA-related sequences were isolated from mouse genomic DNA libraries. Only one of the sequences completely matches the published mouse U7 RNA sequence, whereas the other two apparently represent pseudogenes. The matching sequence represents a functional gene, as it is expressed after microinjection into Xenopus laevis oocytes. Sequence variations of the conserved cis-acting 5' and 3' elements of U RNA genes may partly explain the low abundance of U7 RNA.