Expression of FoxC, FoxF, FoxL1, and FoxQ1 genes in the dogfish scyliorhinus canicula defines ancient and derived roles for fox genes in vertebrate development

In the human genome, members of the FoxC, FoxF, FoxL1, and FoxQ1 gene families are found in two paralagous clusters. Here we characterize all four gene families in the dogfish Seyliorhinus canicula, a member of the cartilaginous fish lineage that diverged before the radiation of osteichthyan vertebrates. We identify two FoxC genes, two FoxF genes, and single FoxQ1 and FoxL1 genes, demonstrating cluster duplication preceded the radiation of gnathostomes. The expression of all six genes was analyzed by in situ hybridization. The results show conserved expression of FoxL1, FoxF, and FoxC genes in different compartments of the mesoderm and of FoxQ1 in pharyngeal endoderm and its derivatives, confirming these as ancient sites of Fox gene expression, and also illustrate multiple cases of lineage-specific expression domains. Comparison to invertebrate chordates shows that the majority of conserved vertebrate expression domains mark tissues that are part of the primitive chordate body plan.

The Fox genes of Branchiostoma floridae

The Fox genes are united by encoding a fork head domain, a deoxyribonucleic acid (DNA)-binding domain of the winged-helix type that marks these genes as encoding transcription factors. Vertebrate Fox genes are classified into 23 subclasses named from FoxA to FoxS. We have surveyed the genome of the amphioxus Branchiostoma floridae, identifying 32 distinct Fox genes representing 21 of these 23 subclasses. The missing subclasses, FoxR and FoxS, are specific to vertebrates, and in addition, B. floridae has one further group, FoxAB, that is not found in vertebrates. Hence, we conclude B. floridae has maintained a high level of Fox gene diversity. Expressed sequence tag and complementary DNA sequence data support the expression of 23 genes. Several linkages between B. floridae Fox genes were noted, including some that have evolved relatively recently via tandem duplication in the amphioxus lineage and others that are more ancient.

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.

Phylogenetic analysis of the pPT23A plasmid family of Pseudomonas syringae

The pPT23A plasmid family of Pseudomonas syringae contains members that contribute to the ecological and pathogenic fitness of their P. syringae hosts. In an effort to understand the evolution of these plasmids and their hosts, we undertook a comparative analysis of the phylogeny of plasmid genes and that of conserved chromosomal genes from P. syringae. In total, comparative sequence and phylogenetic analyses were done utilizing 47 pPT23A family plasmids (PFPs) from 16 pathovars belonging to six genomospecies. Our results showed that the plasmid replication gene (repA), the only gene currently known to be distributed among all the PFPs, had a phylogeny that was distinct from that of the P. syringae hosts of these plasmids and from those of other individual genes on PFPs. The phylogenies of two housekeeping chromosomal genes, those for DNA gyrase B subunit (gyrB) and primary sigma factor (rpoD), however, were strongly associated with genomospecies of P. syringae. Based on the results from this study, we conclude that the pPT23A plasmid family represents a dynamic genome that is mobile among P. syringae pathovars.

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.