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 .

Mutational activation of niche-specific genes provides insight into regulatory networks and bacterial function in a complex environment

The genome of the plant-colonizing bacterium Pseudomonas fluorescens SBW25 harbors a subset of genes that are expressed specifically on plant surfaces. The function of these genes is central to the ecological success of SBW25, but their study poses significant challenges because no phenotype is discernable in vitro. Here, we describe a genetic strategy with general utility that combines suppressor analysis with IVET (SPyVET) and provides a means of identifying regulators of niche-specific genes. Central to this strategy are strains carrying operon fusions between plant environment-induced loci (EIL) and promoterless 'dapB. These strains are prototrophic in the plant environment but auxotrophic on laboratory minimal medium. Regulatory elements were identified by transposon mutagenesis and selection for prototrophs on minimal medium. Approximately 106 mutants were screened for each of 27 strains carrying 'dapB fusions to plant EIL and the insertion point for the transposon determined in approximately 2,000 putative regulator mutants. Regulators were functionally characterized and used to provide insight into EIL phenotypes. For one strain carrying a fusion to the cellulose-encoding wss operon, five different regulators were identified including a diguanylate cyclase, the flagella activator, FleQ, and alginate activator, AmrZ (AlgZ). Further rounds of suppressor analysis, possible by virtue of the SPyVET strategy, revealed an additional two regulators including the activator AlgR, and allowed the regulatory connections to be determined.

In vivo selection for metastasis promoting genes in the mouse

Here, we report the identification of a metastasis promoting factor by a forward genetic screen in mice. A retroviral cDNA library was introduced into the nonmetastatic cancer cell line 168FARN, which was then orthotopically transplanted into mouse mammary fat pads, followed by selection for cells that metastasize to the lung. The genes encoding the disulfide isomerase ERp5 and beta-catenin were found to promote breast cancer invasion and metastasis. Disulfide isomerases (thiol isomerases), which catalyze disulfide bond formation, reduction, and isomerization, have not previously been implicated in cancer cell signaling and tumor metastasis. Overexpression of ERp5 promotes both in vitro migration and invasion and in vivo metastasis of breast cancer cells. These effects were shown to involve activation of ErbB2 and phosphoinositicle 3-kinase (PI3K) pathways through dimerization of ErbB2. Activation of ErbB2 and PI3K subsequently stimulates RhoA and beta-catenin, which mediate the migration and invasion of tumor cells. Inhibition of ErbB2 and PI3K reverses the phenotypes induced by ERp5. Finally, ERp5 was shown to be up-regulated in human surgical samples of invasive breast cancers. These data identify a link between disulfide isomerases and tumor development, and provide a mechanism that modulates ErbB2 and PI3K signaling in the promotion of cancer progression.

Expression of metallothionein genes I and II in bank vole clethrionomys glareolus populations chronically exposed in situ to heavy metals

The expression of two metallothionein genes (Mt-I and Mt-II) in the liver, kidney, and gonad of bank voles collected at four metal-contaminated sites (Cd, Zn, Pb, and Fe) were measured using the quantitative real-time PCR method (QPCR). Relative Mt gene expression was calculated by applying a normalization factor (NF) using the expression of two housekeeping genes, ribosomal 18S and beta-actin. Relative Mt expression in tissues of animals from contaminated sites was up to 54.8-fold higher than those from the reference site for Mt-I and up to 91.6-fold higher for Mt-II. Mt-II gene expression in the livers of bank voles from contaminated sites was higher than Mt-I gene expression. Inversely, Mt-II expression in the kidneys of voles was lower than Mt-I expression. Positive correlations between cadmium levels in the tissues and Mt-I were obtained in all studied tissues. Zinc, which undergoes homeostatic regulation, correlated positively with both Mt-I and Mt-II gene expression only in the kidney. Results showed that animals living in chronically contaminated environments intensively activate detoxifying mechanisms such as metallothionein expression. This is the first time that QPCR techniques to measure MT gene expression have been applied to assess the impact of environmental metal pollution on field collected bank voles.

Expression of target and reference genes in Daphnia magna exposed to ibuprofen

Background: Transcriptomic techniques are now being applied in ecotoxicology and toxicology to measure the impact of stressors and develop understanding of mechanisms of toxicity. Microarray technology in particular offers the potential to measure thousands of gene responses simultaneously. However, it is important that microarrays responses should be validated, at least initially, using real-time quantitative polymerase chain reaction (QPCR). The accurate measurement of target gene expression requires normalisation to an invariant internal control e. g., total RNA or reference genes. Reference genes are preferable, as they control for variation inherent in the cDNA synthesis and PCR. However, reference gene expression can vary between tissues and experimental conditions, which makes it crucial to validate them prior to application. Results: We evaluated 10 candidate reference genes for QPCR in Daphnia magna following a 24 h exposure to the non-steroidal anti-inflammatory drug (NSAID) ibuprofen (IB) at 0, 20, 40 and 80 mg IB l(-1). Six of the 10 candidates appeared suitable for use as reference genes. As a robust approach, we used a combination normalisation factor (NF), calculated using the geNorm application, based on the geometric mean of three selected reference genes: glyceraldehyde-3-phosphate dehydrogenase, ubiquitin conjugating enzyme and actin. The effects of normalisation are illustrated using as target gene leukotriene B4 12-hydroxydehydrogenase (Ltb4dh), which was upregulated following 24 h exposure to 63-81 mg IB l(-1). Conclusions: As anticipated, use of the NF clarified the response of Ltb4dh in daphnids exposed to sublethal levels of ibuprofen. Our findings emphasise the importance in toxicogenomics of finding and applying invariant internal QPCR control(s) relevant to the study conditions.

Comparative expression analysis of transcription factor genes in the endostyle of invertebrate chordates

The endostyle of invertebrate chordates is a pharyngeal organ that is thought to be homologous with the follicular thyroid of vertebrates. Although thyroid-like features such as iodine-concentrating and peroxidase activities are located in the dorsolateral part of both ascidian and amphioxus endostyles, the structural organization and numbers of functional units are different. To estimate phylogenetic relationships of each functional zone with special reference to the evolution of the thyroid, we have investigated, in ascidian and amphioxus, the expression patterns of thyroid-related transcription factors such as TTF-2/MoxE4 and Pax2/5/8, as well as the forkhead transcription factors FoxQ1 and FoxA. Comparative gene expression analyses depicted an overall similarity between ascidians and amphioxus endostyles, while differences in expression patterns of these genes might be specifically related to the addition or elimination of a pair of glandular zones. Expressions of Ci-FoxE and BbFoxE4 suggest that the ancestral FoxE class might have been recruited for the formation of thyroid-like region in a possible common ancestor of chordates. Furthermore, coexpression of FoxE4, Pax2/5/8, and TPO in the dorsolateral part of both ascidian and amphioxus endostyles suggests that genetic basis of the thyroid function was already in place before the vertebrate lineage. (c) 2005 Wiley-Liss, Inc.

Chromosomal mapping of ANTP class homeobox genes in amphioxus: piecing together ancestral genomes

Homeobox genes encode DNA-binding proteins, many of which are implicated in the control of embryonic development. Evolutionarily, most homeobox genes fall into two related clades: the ANTP and the PRD classes. Some genes in ANTP class, notably Hox, ParaHox, and NK genes, have an intriguing arrangement into physical clusters. To investigate the evolutionary history of these gene clusters, we examined homeobox gene chromosomal locations in the cephalochordate amphioxus, Branchiostoma floridae. We deduce that 22 amphioxus ANTP class homeobox genes localize in just three chromosomes. One contains the Hox cluster plus AmphiEn, AmphiMnx, and AmphiDll. The ParaHox cluster resides in another chromosome, whereas a third chromosome contains the NK type homeobox genes, including AmphiMsx and ArnphiTlx. By comparative analysis we infer that clustering of ANTP class homeobox genes evolved just once, during a series of extensive cis-duplication events of genes early in animal evolution. A trans-duplication event occurred later to yield the Hox and ParaHox gene clusters on different chromosomes. The results obtained have implications for understanding the origin of homeobox gene clustering, the diversification of the ANTP class of homeobox genes, and the evolution of animal genomes.

Protochordate Zic genes define primitive somite compartments and highlight molecular changes underlying neural crest evolution

The vertebrate Zic gene family encodes C2H2 zinc finger transcription factors closely related to the Gli proteins. Zic genes are expressed in multiple areas of developing vertebrate embryos, including the dorsal neural tube where they act as potent neural crest inducers. Here we describe the characterization of a Zic ortholog from the amphioxus Branchiostoma floridae and further describe the expression of a Zic ortholog from the ascidian Ciona intestinalis. Molecular phylogenetic analysis and sequence comparisons suggest the gene duplications that formed the vertebrate Zic family were specific to the vertebrate lineage. In Ciona maternal CiZic/Ci-macho1 transcripts are localized during cleavage stages by asymmetric cell division, whereas zygotic expression by neural plate cells commences during neurulation. The amphioxus Zic ortholog AmphiZic is expressed in dorsal mesoderm and ectoderm during gastrulation, before being eliminated first from midline cells and then from all neurectoderm during neurulation. After neurulation, expression is reactivated in the dorsal neural tube and dorsolateral somite. Comparison of CiZic and AmphiZic expression with vertebrate Zic expression leads to two main conclusions. First, Zic expression allows us to define homologous compartments between vertebrate and amphioxus somites, showing primitive subdivision of vertebrate segmented mesoderm. Second, we show that neural Zic expression is a chordate synapomorphy, whereas the precise pattern of neural expression has evolved differently on the different chordate lineages. Based on these observations we suggest that a change in Zic regulation, specifically the evolution of a dorsal neural expression domain in vertebrate neurulae, was an important step in the evolution of the neural crest.

A genomewide survey of developmentally relevant genes in Ciona intestinalis - III. Genes for Fox, ETS, nuclear receptors and NF kappa B

A survey against the draft genome sequence and the cDNA/EST database of Ciona intestinalis identified a number of genes encoding transcription factors regulating a variety of processes including development. In the present study, we describe almost complete sets of genes for Fox, ETS-domain transcription factors, nuclear receptors, and NFkappaB as well as other factors regulating NFkappaB activity, with their phylogenetic nature. Vertebrate Fox transcription factors are currently delineated into 17 subfamilies: FoxA to FoxQ. The present survey yielded 29 genes of this family in the Ciona genome, 24 of which were Ciona orthologues of known Fox genes. In addition, we found 15 ETS aenes, 17 nuclear receptor genes, and several NFkappaB signaling pathway genes in the Ciona genome. The number of Ciona genes in each family is much smaller than that of vertebrates, which represents a simplified feature of the ascidian genome. For example, humans have two NFkappaB genes, three Rel genes, and five NFAT genes, while Ciona has one gene for each family. The Ciona genome also contains smaller numbers of genes for the NFkappaB regulatory system, i.e. after the split of ascidians/vertebrates, vertebrates evolved a more complex NFkappaB system. The present results therefore provide molecular information for the investigation of complex developmental processes, and an insight into chordate evolution.

Characterization of the 23S and 5S rRNA genes and 23S-5S intergenic spacer region (ITS-2) of Photobacterium damselae

The 23S ribosomal RNA (rRNA) gene has been sequenced in strains of the fish pathogens Photobacterium damselae subsp. damselae (ATCC 33539) and subsp. piscicida (ATCC 29690), showing that 3 nucleotide positions are clearly different between subspecies. In addition, the 5S rRNA gene plus the intergenic spacer region between the 23S and 5S rRNA genes (ITS-2) were amplified, cloned and sequenced for the 2 reference strains as well as the field isolates RG91 (subsp. damselae) and DI21 (subsp. piscicida). A 100% similarity was found for the consensus 5S rRNA gene sequence in the 2 subspecies, although some microheterogeneity was detected as inter-cistronic variability within the same chromosome. Sequence analysis of the spacer region between the 23S and 5S rRNA genes revealed 2 conserved and 3 variable nucleotide sequence blocks, and 4 different modular organizations were found. The ITS-2 spacer region exhibited both inter-subspecies and inter-cistronic polymorphism, with a mosaic-like structure. The EMBL accession numbers for the 23S, 5S and ITS-2 sequences are: P. damselae subsp. piscicida 5S gene (AJ274379), P. damselae subsp. damselae 23S gene (Y18520), subsp. piscicida 23S gene (Y17901), R damselae subsp. piscicida ITS-2 (AJ250695, AJ250696), P. damselae subsp. damselae ITS-2 (AJ250697, AJ250698).

Molecular cloning and comparative analysis of four beta-galactosidase genes from Bifidobacterium bifidum NCIMB41171

Bifidobacterium bifidum NCIMB41171 carries four genes encoding different beta-galactosidases. One of them, named bbgIII, consisted of an open reading frame of 1,935 amino acid (a.a.) residues encoding a protein with a multidomain structure, commonly identified on cell wall bound enzymes, having a signal peptide, a membrane anchor, FIVAR domains, immunoglobulin Ig-like and discoidin-like domains. The other three genes, termed bbgI, bbgII and bbgIV, encoded proteins of 1,291, 689 and 1,052 a.a. residues, respectively, which were most probably intracellularly located. Two cases of protein evolution between strains of the same species were identified when the a.a. sequences of the BbgI and BbgIII were compared with homologous proteins from B. bifidum DSM20215. The homologous proteins were found to be differentiated at the C-terminal a.a. part either due to a single nucleotide insertion or to a whole DNA sequence insertion, respectively. The bbgIV gene was located in a gene organisation surrounded by divergently transcribed genes putatively for sugar transport (galactoside-symporter) and gene regulation (LacI-transcriptional regulator), a structure that was found to be highly conserved in B. longum, B. adolescentis and B. infantis, suggesting optimal organisation shared amongst those species.

Genistein affects the expression of genes involved in blood pressure regulation and angiogenesis in primary human endothelial cells

Background: Several lines of evidence suggest that the dietary isoflavone genistein (Gen) has beneficial effects with regard to cardiovascular disease and in particular on aspects related to blood pressure and angiogenesis. The biological action of Gen may be, at Least in part, attributed to its ability to affect cell signalling and response. However, so far, most of the molecular mechanisms underlying the activity of Gen in the endothelium are unknown. Methods and results: To examine the transcriptional response to 2.5 mu M Gen on primary human endothelial cells (HUVEC), we applied cDNA array technology both under baseline condition and after treatment with the pro-atherogenic stimulus, copper-oxidized LDL. The alteration of the expression patterns of individual transcripts was substantiated using either RT-PCR or Northern blotting. Gen significantly affected the expression of genes encoding for proteins centrally involved in the vascular tone such as endothelin-converting enzyme-1, endothetin-2, estrogen related receptor a and atria[ natriuretic peptide receptor A precursor. Furthermore, Gen countered the effect of oxLDL on mRNA levels encoding for vascular endothelial growth factor receptor 165, types 1 and 2. Conclusions: Our data indicate that physiologically achievable levels of Gen change the expression of mRNA encoding for proteins involved in the control of blood pressure under baseline conditions and reduce the angiogenic response to oxLDL in the endothelium. (c) 2005 Elsevier B.V. All rights reserved.

Polymorphisms in dopamine system genes are associated with individual differences in attention in infancy

Knowledge about the functional status of the frontal cortex in infancy is limited. This study investigated the effects of polymorphisms in four dopamine system genes on performance in a task developed to assess such functioning, the Freeze-Frame task, at 9 months of age. Polymorphisms in the catechol-O-methyltransferase (COMT) and the dopamine D4 receptor (DRD4) genes are likely to impact directly on the functioning of the frontal cortex, whereas polymorphisms in the dopamine D2 receptor (DRD2) and dopamine transporter (DAT1) genes might influence frontal cortex functioning indirectly via strong frontostriatal connections. A significant effect of the COMT valine158methionine (Val158Met) polymorphism was found. Infants with the Met/Met genotype were significantly less distractible than infants with the Val/Val genotype in Freeze-Frame trials presenting an engaging central stimulus. In addition, there was an interaction with the DAT1 3′ variable number of tandem repeats polymorphism; the COMT effect was present only in infants who did not have two copies of the DAT1 10-repeat allele. These findings indicate that dopaminergic polymorphisms affect selective aspects of attention as early as infancy and further validate the Freeze-Frame task as a frontal cortex task.

Sleeping with half open eyes: dreams and realities in The Cry of the Owl

This article is a close analysis of The Cry of the Owl (Thraves, 2009). It is also part of larger project to bring together traditions of detailed criticism with those of production history, which culminates in second article on the film due to be published in 2011. The detail of the argument concerns analysing a range of the film’s key signifying systems, with a particular interest in the way the film explores the gap between images / impressions and characters’ realities; engages in a complex way with generic traditions and modes of address; establishes complex patterns of connection and contrast through blocking, camera strategies and narrative structure.

A specific group of genes respond to cold dehydration stress in cut Alstroemeria flowers whereas ambient dehydration stress accelerates developmental senescence expression patterns

Petal development and senescence entails a normally irreversible process. It starts with petal expansion and pigment production, and ends with nutrient remobilization and ultimately cell death. In many species this is accompanied by petal abscission. Post-harvest stress is an important factor in limiting petal longevity in cut flowers and accelerates some of the processes of senescence such as petal wilting and abscission. However, some of the effects of moderate stress in young flowers are reversible with appropriate treatments. Transcriptomic studies have shown that distinct gene sets are expressed during petal development and senescence. Despite this, the overlap in gene expression between developmental and stress-induced senescence in petals has not been fully investigated in any species. Here a custom-made cDNA microarray from Alstroemeria petals was used to investigate the overlap in gene expression between developmental changes (bud to first sign of senescence) and typical post-harvest stress treatments. Young flowers were stressed by cold or ambient temperatures without water followed by a recovery and rehydration period. Stressed flowers were still at the bud stage after stress treatments. Microarray analysis showed that ambient dehydration stress accelerates many of the changes in gene expression patterns that would normally occur during developmental senescence. However, a higher proportion of gene expression changes in response to cold stress were specific to this stimulus and not senescence related. The expression of 21 transcription factors was characterized, showing that overlapping sets of regulatory genes are activated during developmental senescence and by different stresses.