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.

Functional genes and gene array analysis as tools for monitoring hydrocarbon biodegradation

Bioremediation, which is the exploitation of the intrinsic ability of environmental microbes to degrade and remove harmful compounds from nature, is considered to be an environmentally sustainable and cost-effective means for environmental clean-up. However, a comprehensive understanding of the biodegradation potential of microbial communities and their response to decontamination measures is required for the effective management of bioremediation processes. In this thesis, the potential to use hydrocarbon-degradative genes as indicators of aerobic hydrocarbon biodegradation was investigated. Small-scale functional gene macro- and microarrays targeting aliphatic, monoaromatic and low molecular weight polyaromatic hydrocarbon biodegradation were developed in order to simultaneously monitor the biodegradation of mixtures of hydrocarbons. The validity of the array analysis in monitoring hydrocarbon biodegradation was evaluated in microcosm studies and field-scale bioremediation processes by comparing the hybridization signal intensities to hydrocarbon mineralization, real-time polymerase chain reaction (PCR), dot blot hybridization and both chemical and microbiological monitoring data. The results obtained by real-time PCR, dot blot hybridization and gene array analysis were in good agreement with hydrocarbon biodegradation in laboratory-scale microcosms. Mineralization of several hydrocarbons could be monitored simultaneously using gene array analysis. In the field-scale bioremediation processes, the detection and enumeration of hydrocarbon-degradative genes provided important additional information for process optimization and design. In creosote-contaminated groundwater, gene array analysis demonstrated that the aerobic biodegradation potential that was present at the site, but restrained under the oxygen-limited conditions, could be successfully stimulated with aeration and nutrient infiltration. During ex situ bioremediation of diesel oil- and lubrication oil-contaminated soil, the functional gene array analysis revealed inefficient hydrocarbon biodegradation, caused by poor aeration during composting. The functional gene array specifically detected upper and lower biodegradation pathways required for complete mineralization of hydrocarbons. Bacteria representing 1 % of the microbial community could be detected without prior PCR amplification. Molecular biological monitoring methods based on functional genes provide powerful tools for the development of more efficient remediation processes. The parallel detection of several functional genes using functional gene array analysis is an especially promising tool for monitoring the biodegradation of mixtures of hydrocarbons.

Rrp1B gene polymorphism (1307T > C) in metastatic progression of breast cancer

Rrp1B (ribosomal RNA processing1 homolog B) is a novel candidate metastasis modifier gene in breast cancer. Functional gene assays demonstrated that a physical and functional interaction existing between Rrp1b and metastasis modifier gene SIPA1 causes reduction in the tumor growth and metastatic potential. Ectopic expression of Rrp1B modulates various metastasis predictive extra cellular matrix (ECM) genes associated with tumor suppression. The aim of this study is to determine the functional significance of single nucleotide polymorphism (SNP) in human Rrp1B gene (1307 T > C; rs9306160) with breast cancer development and progression. The study consists of 493 breast cancer cases recruited from Nizam's Institute of Medical Sciences, Hyderabad, and 558 age-matched healthy female controls from rural and urban areas. Genomic DNA was isolated by non-enzymatic method. Genotyping was done by amplification refractory mutation system (ARMS-PCR) method. Genotypes were reconfirmed by sequencing and results were analyzed statistically. We have performed Insilco analysis to know the RNA secondary structure by using online tool m fold. The TT genotype and T allele frequencies of Rrp1B1307 T > C polymorphism were significantly elevated in breast cancer (chi (2); p = < 0.008) cases compared to controls under different genetic models. The presence of T allele had conferred 1.75-fold risk for breast cancer development (OR = 1.75; 95 % CI = 1.15-2.67). The frequency of TT genotype of Rrp1b 1307T > C polymorphism was significantly elevated in obese patients (chi (2); p = 0.008) and patients with advanced disease (chi (2); p = 0.01) and with increased tumor size (chi (2); p = 0.01). Moreover, elevated frequency of T allele was also associated with positive lymph node status (chi (2); p = 0.04) and Her2 negative receptor status (chi (2); p = 0.006). Presence of Rrp1b1307TT genotype and T allele confer strong risk for breast cancer development and progression.

A de novo originated gene depresses budding yeast mating pathway and is repressed by the protein encoded by its antisense strand

Recent transcription profiling studies have revealed an unexpectedly large proportion of antisense transcripts in eukaryotic genomes. These antisense genes seem to regulate gene expression by interacting with sense genes. Previous studies have focused on the non-coding antisense genes, but the possible regulatory role of the antisense protein is poorly understood. In this study, we found that a protein encoded by the antisense gene ADF1 acts as a transcription suppressor, regulating the expression of sense gene MDF1 in Saccharomyces cerevisiae. Based on the evolutionary, genetic, cytological and biochemical evidence, we show that the protein-coding sense gene MDF1 most likely originated de novo from a previously non-coding sequence and can significantly suppress the mating efficiency of baker's yeast in rich medium by binding MAT alpha 2 and thus promote vegetative growth. These results shed new light on several important issues, including a new sense-antisense interaction mechanism, the de novo origination of a functional gene, and the regulation of yeast mating pathway.

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.

The functional role of the tumor suppressor MMAC /PTEN in glioblastoma multiforme and prostate cancer

Investigations into the molecular basis of glioblastoma multiforme led to the identification of a putative tumor suppressor gene, MMAC/ PTEN. Initial studies implicated MMAC/PTEN in many different tumor types, and identified a protein phosphatase motif in its sequence. This project aimed to identify the biological and biochemical functions of MMAC/PTEN by transiently expressing the gene in cancer cells that lack a functional gene product. ^ Expression of MMAC/PTEN mildly suppressed the growth of U251 human glioma cells and abrogated the growth advantage mediated by overexpression of the epidermal growth factor receptor (EGFR). Immunoblotting demonstrated that MMAC/PTEN expression did not affect the phosphorylation of the EGFR itself, or the intermediates of several downstream signaling pathways. However, MMAC/PTEN expression significantly reduced the phosphorylation and catalytic activity of the proto-oncogene Akt/PKB. While Akt/PKB regulates the survival of many cell types, expression of MMAC/PTEN did not induce apoptosis in adherent U251 cells. Instead, MMAC/PTEN expression sensitized the cells to apoptosis when maintained in suspension (anoikis). As the survival of suspended cells is one of the hallmarks leading to metastasis, MMAC/PTEN expression was examined in a system in which metastasis is more clinically relevant, prostate cancer. ^ Expression of MMAC/PTEN in both LNCaP and PC3-P human prostate cancer cells specifically inhibited Akt/PKB phosphorylation. MMAC/PTEN expression in LNCaP cells resulted in a profound inhibition of growth that was significantly greater than that achieved with expression of p53. Expression of MMAC/PTEN in PC3-P cells resulted in greater growth inhibition than was observed in U251 glioma cells, but less than was observed in LNCaP cells, or upon p53 expression. To determine if MMAC/PTEN could function as a tumor suppressor in vivo, the effects of MMAC/PTEN expression on PC3-P cells implanted orthotopically in nude mice were examined. The ex-vivo expression of MMAC/PTEN did not decrease tumor incidence, but it did significantly decrease tumor size and metastasis. In-vivo expression of MMAC/PTEN in pre-established PC3-P tumors did not significantly inhibit tumor incidence or size, but did inhibit metastasis formation. ^ These studies demonstrate that MMAC/PTEN is a novel and important tumor suppressor gene, which functions to downregulate an important cell survival signaling pathway. ^

Insecticide resistance resulting from an absence of target-site gene product

Genetic changes in insects that lead to insecticide resistance include point mutations and up-regulation/amplification of detoxification genes. Here, we report a third mechanism, resistance caused by an absence of gene product. Mutations of the Methoprene-tolerant (Met) gene of Drosophila melanogaster result in resistance to both methoprene, a juvenile hormone (JH) agonist insecticide, and JH. Previous results have demonstrated a mechanism of resistance involving an intracellular JH binding protein that has reduced ligand affinity in Met flies. We show that a γ-ray induced allele, Met27, completely lacks Met transcript during the insecticide-sensitive period in development. Although Met27 homozygotes have reduced oogenesis, they are viable, demonstrating that Met is not a vital gene. Most target-site resistance genes encode vital proteins and thus have few mutational changes that permit both resistance and viability. In contrast, resistance genes such as Met that encode nonvital insecticide target proteins can have a variety of mutational changes that result in an absence of functional gene product and thus should show higher rates of resistance evolution.

Scrambling of the actin I gene in two Oxytricha species.

The DNA in a germ-line nucleus (a micronucleus) undergoes extensive processing when it develops into a somatic nucleus (a macronucleus) after cell mating in hypotrichous ciliates. Processing includes destruction of a large amount of spacer DNA between genes and excision of gene-sized molecules from chromosomes. Before processing, micronuclear genes are interrupted by numerous noncoding segments called internal eliminated sequences (IESs). The IESs are excised and destroyed, and the retained macro-nuclear-destined sequences (MDSs) are spliced. MDSs in some micronuclear genes are not in proper order and must be reordered during processing to create functional gene-sized molecules for the macronucleus. Here we report that the micronuclear actin I gene in Oxytricha trifallax WR consists of 10 MDSs and 9 IESs compared to the previously reported 9 MDSs and 8 IESs in the micronuclear actin I gene of Oxytricha nova. The MDSs in the actin I gene are scrambled in a similar pattern in the two species, but the positions of MDS-IES junctions are shifted by up to 14 bp for scrambled and 138 bp for the nonscrambled MDSs. The shifts in MDS-IES junctions create differences in the repeat sequences that are believed to guide MDS splicing. Also, the sizes and sequences of IESs in the micronuclear actin I genes are different in the two Oxytricha species. These observations give insight about the possible origins of IES insertion and MDS scrambling in evolution and show the extraordinary malleability of the germ-line DNA in hypotrichs.