Sequencing and identification of expressed Schistosoma mansoni genes by random selection of cDNA clones from a directional library

We have initiated a gene discovery program in Schistosoma mansoni based on the technique of Expressed Sequence Tags (ESTs), i.e. partial sequences of cDNAs obtained from single passes in automatic DNA sequencers. ESTs can be used to identify genese onf the basis of their homology whith sequences from other species deposited in DNA or protein databases. Trasncripts with sequences without matches in teh databases may represent novel parasite-specific genes. This approach has shown to be very efficient and in less than two years a broad range of novel genes has already been ascertained, more than doubling the number of known S. mansoni genes.

trEST, trGEN and Hits: access to databases of predicted protein sequences.

High throughput genome (HTG) and expressed sequence tag (EST) sequences are currently the most abundant nucleotide sequence classes in the public database. The large volume, high degree of fragmentation and lack of gene structure annotations prevent efficient and effective searches of HTG and EST data for protein sequence homologies by standard search methods. Here, we briefly describe three newly developed resources that should make discovery of interesting genes in these sequence classes easier in the future, especially to biologists not having access to a powerful local bioinformatics environment. trEST and trGEN are regularly regenerated databases of hypothetical protein sequences predicted from EST and HTG sequences, respectively. Hits is a web-based data retrieval and analysis system providing access to precomputed matches between protein sequences (including sequences from trEST and trGEN) and patterns and profiles from Prosite and Pfam. The three resources can be accessed via the Hits home page (http://hits. isb-sib.ch).

Social insect genomes exhibit dramatic evolution in gene composition and regulation while preserving regulatory features linked to sociality.

Genomes of eusocial insects code for dramatic examples of phenotypic plasticity and social organization. We compared the genomes of seven ants, the honeybee, and various solitary insects to examine whether eusocial lineages share distinct features of genomic organization. Each ant lineage contains ∼4000 novel genes, but only 64 of these genes are conserved among all seven ants. Many gene families have been expanded in ants, notably those involved in chemical communication (e.g., desaturases and odorant receptors). Alignment of the ant genomes revealed reduced purifying selection compared with Drosophila without significantly reduced synteny. Correspondingly, ant genomes exhibit dramatic divergence of noncoding regulatory elements; however, extant conserved regions are enriched for novel noncoding RNAs and transcription factor-binding sites. Comparison of orthologous gene promoters between eusocial and solitary species revealed significant regulatory evolution in both cis (e.g., Creb) and trans (e.g., fork head) for nearly 2000 genes, many of which exhibit phenotypic plasticity. Our results emphasize that genomic changes can occur remarkably fast in ants, because two recently diverged leaf-cutter ant species exhibit faster accumulation of species-specific genes and greater divergence in regulatory elements compared with other ants or Drosophila. Thus, while the "socio-genomes" of ants and the honeybee are broadly characterized by a pervasive pattern of divergence in gene composition and regulation, they preserve lineage-specific regulatory features linked to eusociality. We propose that changes in gene regulation played a key role in the origins of insect eusociality, whereas changes in gene composition were more relevant for lineage-specific eusocial adaptations.

Identification of Transcribed Sequences (ESTs) in the Trypanosoma cruzi Genome Project

Random single pass sequencing of cDNA fragments, also known as generation of Expressed Sequence Tags (ESTs), has been highly successful in the study of the gene content of higher organisms, and forms an integral part of most genome projects, with the objective to identify new genes and targets for disease control and prevention and to generate mapping probes. In the Trypanosoma cruzi genome project, EST sequencing has also been a starting point, and here we report data on the first 797 sequences obtained, partly from a CL Brener epimastigote non-normalized library, partly on a normalized library. Only around 30% of the sequences obtained showed similarity with Genbank and dbEST databases, half of which with sequences already reported for T. cruzi.

Postinfarction heart failure in rats is associated with upregulation of GLUT-1 and downregulation of genes of fatty acid metabolism.

OBJECTIVES: Increasing evidence suggests that left ventricular remodeling is associated with a shift from fatty acid to glucose metabolism for energy production. The aim of this study was to determine whether left ventricular remodeling with and without late-onset heart failure after myocardial infarction is associated with regional changes in the expression of regulatory proteins of glucose or fatty acid metabolism. METHODS: Myocardial infarction was induced in rats by ligation of the left anterior descending coronary artery (LAD). In infarcted and sham-operated hearts the peri-infarction region (5-mm zone surrounding the region at risk), the interventricular septum and the right ventricular free wall were separated for analysis. RESULTS: At 8 and 20 weeks after LAD ligation, the peri-infarction region and the septum exhibited marked re-expression of atrial natriuretic factor [+252+/-37 and +1093+/-279%, respectively, in the septum (P<0.05)] and of alpha-smooth muscle actin [+34+/-10 and +43+/-14%, respectively, in the septum (P<0.05)]. At 8 weeks, when left ventricular hypertrophy was present without signs of heart failure, myocardial mRNA expression of glucose transporters (GLUT-1 and GLUT-4) was not altered, whereas mRNA expression of medium-chain acyl-CoA dehydrogenase (MCAD) was significantly reduced in the peri-infarction region (-25+/-7%; P<0.05). In hearts exhibiting heart failure 20 weeks after infarct-induction there was a change in all three ventricular regions of both mRNA and protein content of GLUT-1 [+72+/-28 and +121+/-15%, respectively, in the peri-infarction region (P<0.05)] and MCAD [-29+/-9 and -56+/-4%, respectively, in the peri-infarction region (P<0.05)]. CONCLUSION: In rats with large myocardial infarction, progression from compensated remodeling to overt heart failure is associated with upregulation of GLUT-1 and downregulation of MCAD in both the peri-infarction region and the septum.

Novel H3K4me3 marks are enriched at human- and chimpanzee-specific cytogenetic structures.

Human and chimpanzee genomes are 98.8% identical within comparable sequences. However, they differ structurally in nine pericentric inversions, one fusion that originated human chromosome 2, and content and localization of heterochromatin and lineage-specific segmental duplications. The possible functional consequences of these cytogenetic and structural differences are not fully understood and their possible involvement in speciation remains unclear. We show that subtelomeric regions-regions that have a species-specific organization, are more divergent in sequence, and are enriched in genes and recombination hotspots-are significantly enriched for species-specific histone modifications that decorate transcription start sites in different tissues in both human and chimpanzee. The human lineage-specific chromosome 2 fusion point and ancestral centromere locus as well as chromosome 1 and 18 pericentric inversion breakpoints showed enrichment of human-specific H3K4me3 peaks in the prefrontal cortex. Our results reveal an association between plastic regions and potential novel regulatory elements.

Hard-wired heterogeneity in blood stem cells revealed using a dynamic regulatory network model.

MOTIVATION: Combinatorial interactions of transcription factors with cis-regulatory elements control the dynamic progression through successive cellular states and thus underpin all metazoan development. The construction of network models of cis-regulatory elements, therefore, has the potential to generate fundamental insights into cellular fate and differentiation. Haematopoiesis has long served as a model system to study mammalian differentiation, yet modelling based on experimentally informed cis-regulatory interactions has so far been restricted to pairs of interacting factors. Here, we have generated a Boolean network model based on detailed cis-regulatory functional data connecting 11 haematopoietic stem/progenitor cell (HSPC) regulator genes. RESULTS: Despite its apparent simplicity, the model exhibits surprisingly complex behaviour that we charted using strongly connected components and shortest-path analysis in its Boolean state space. This analysis of our model predicts that HSPCs display heterogeneous expression patterns and possess many intermediate states that can act as 'stepping stones' for the HSPC to achieve a final differentiated state. Importantly, an external perturbation or 'trigger' is required to exit the stem cell state, with distinct triggers characterizing maturation into the various different lineages. By focusing on intermediate states occurring during erythrocyte differentiation, from our model we predicted a novel negative regulation of Fli1 by Gata1, which we confirmed experimentally thus validating our model. In conclusion, we demonstrate that an advanced mammalian regulatory network model based on experimentally validated cis-regulatory interactions has allowed us to make novel, experimentally testable hypotheses about transcriptional mechanisms that control differentiation of mammalian stem cells. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Teleost fish larvae adapt to dietary arachidonic acid supply through modulation of the expression of lipid metabolism and stress response genes

Dietary fatty acid supply can affect stress response in fish during early development. Although knowledge on the mechanisms involved in fatty acid regulation of stress tolerance is scarce, it has often been hypothesised that eicosanoid profiles can influence cortisol production. Genomic cortisol actions are mediated by cytosolic receptors which may respond to cellular fatty acid signalling. An experiment was designed to test the effects of feeding gilthead sea-bream larvae with four microdiets, containing graded arachidonic acid (ARA) levels (0·4, 0·8, 1·5 and 3·0 %), on the expression of genes involved in stress response (steroidogenic acute regulatory protein, glucocorticoid receptor and phosphoenolpyruvate carboxykinase), lipid and, particularly, eicosanoid metabolism (hormone-sensitive lipase, PPARα, phospholipase A2, cyclo-oxygenase-2 and 5-lipoxygenase), as determined by real-time quantitative PCR. Fish fatty acid phenotypes reflected dietary fatty acid profiles. Growth performance, survival after acute stress and similar whole-body basal cortisol levels suggested that sea-bream larvae could tolerate a wide range of dietary ARA levels. Transcription of all genes analysed was significantly reduced at dietary ARA levels above 0·4 %. Nonetheless, despite practical suppression of phospholipase A2 transcription, higher leukotriene B4 levels were detected in larvae fed 3·0 % ARA, whereas a similar trend was observed regarding PGE2 production. The present study demonstrates that adaptation to a wide range of dietary ARA levels in gilthead sea-bream larvae involves the modulation of the expression of genes related to eicosanoid synthesis, lipid metabolism and stress response. The roles of ARA, other polyunsaturates and eicosanoids as signals in this process are discussed.

Mouse interleukin-2 receptor alpha gene expression. Delimitation of cis-acting regulatory elements in transgenic mice and by mapping of DNase-I hypersensitive sites.

The alpha chain of the interleukin-2 receptor (IL-2R alpha) is a key regulator of lymphocyte proliferation. To analyze the mechanisms controlling its expression in normal cells, we used the 5'-flanking region (base pairs -2539/+93) of the mouse gene to drive chloramphenicol acetyltransferase expression in four transgenic mouse lines. Constitutive transgene activity was restricted to lymphoid organs. In mature T lymphocytes, transgene and endogenous IL-2R alpha gene expression was stimulated by concanavalin A and up-regulated by IL-2 with very similar kinetics. In thymic T cell precursors, IL-1 and IL-2 cooperatively induced transgene and IL-2R alpha gene expression. These results show that regulation of the endogenous IL-2R alpha gene occurs mainly at the transcriptional level. They demonstrate that cis-acting elements in the 5'-flanking region present in the transgene confer correct tissue specificity and inducible expression in mature T cells and their precursors in response to antigen, IL-1, and IL-2. In a complementary approach, we screened the 5' end of the endogenous IL-2R alpha gene for DNase-I hypersensitive sites. We found three lymphocyte specific DNase-I hypersensitive sites. Two, at -0.05 and -5.3 kilobase pairs, are present in resting T cells. A third site appears at -1.35 kilobase pairs in activated T cells. It co-localizes with IL-2-responsive elements identified by transient transfection experiments.

Promoter architecture of mouse olfactory receptor genes.

Odorous chemicals are detected by the mouse main olfactory epithelium (MOE) by about 1100 types of olfactory receptors (OR) expressed by olfactory sensory neurons (OSNs). Each mature OSN is thought to express only one allele of a single OR gene. Major impediments to understand the transcriptional control of OR gene expression are the lack of a proper characterization of OR transcription start sites (TSSs) and promoters, and of regulatory transcripts at OR loci. We have applied the nanoCAGE technology to profile the transcriptome and the active promoters in the MOE. nanoCAGE analysis revealed the map and architecture of promoters for 87.5% of the mouse OR genes, as well as the expression of many novel noncoding RNAs including antisense transcripts. We identified candidate transcription factors for OR gene expression and among them confirmed by chromatin immunoprecipitation the binding of TBP, EBF1 (OLF1), and MEF2A to OR promoters. Finally, we showed that a short genomic fragment flanking the major TSS of the OR gene Olfr160 (M72) can drive OSN-specific expression in transgenic mice.

Genes as leaders and followers in evolution.

A major question for the study of phenotypic evolution is whether intra- and interspecific diversity originates directly from genetic variation, or instead, as plastic responses to environmental influences initially, followed later by genetic change. In species with discrete alternative phenotypes, evolutionary sequences can be inferred from transitions between environmental and genetic phenotype control, and from losses of phenotypic alternatives. From the available evidence, sequences appear equally probable to start with genetic polymorphism as with polyphenism, with a possible dominance of one or the other for specific trait types. We argue in this review that to evaluate the prevalence of each route, an investigation of both genetic and environmental cues for phenotype determination in several related rather than in isolated species is required.

Detection of hilA gene sequences in serovars of Salmonella enterica sufigbspecies enterica

hilA gene promoter, component of the Salmonella Pathogenicity Island 1, has been found in Salmonella serovar Typhimurium, being important for the regulation of type III secretion apparatus genes. We detected hilA gene sequences in Salmonella serovars Typhi, Enteritidis, Choleraesuis, Paratyphi A and B, and Pullorum, by polymerase chain reaction (PCR) and hybridization techniques. The primers to carry out PCR were designed according to hilA sequence. A low stringency hybridization with the probe pVV441 (hilA open-reading-frame plasmid) was carried out. To find hilA gene sequences in other Salmonella sp. suggest that these serovars could have similar sequences of this kind of virulence genes.

Genetic analysis of phenoxyalkanoic acid degradation in Sphingomonas herbicidovorans MH.

Phenoxyalkanoic acid degradation is well studied in Beta- and Gammaproteobacteria, but the genetic background has not been elucidated so far in Alphaproteobacteria. We report the isolation of several genes involved in dichlor- and mecoprop degradation from the alphaproteobacterium Sphingomonas herbicidovorans MH and propose that the degradation proceeds analogously to that previously reported for 2,4-dichlorophenoxyacetic acid (2,4-D). Two genes for alpha-ketoglutarate-dependent dioxygenases, sdpA(MH) and rdpA(MH), were found, both of which were adjacent to sequences with potential insertion elements. Furthermore, a gene for a dichlorophenol hydroxylase (tfdB), a putative regulatory gene (cadR), two genes for dichlorocatechol 1,2-dioxygenases (dccA(I/II)), two for dienelactone hydrolases (dccD(I/II)), part of a gene for maleylacetate reductase (dccE), and one gene for a potential phenoxyalkanoic acid permease were isolated. In contrast to other 2,4-D degraders, the sdp, rdp, and dcc genes were scattered over the genome and their expression was not tightly regulated. No coherent pattern was derived on the possible origin of the sdp, rdp, and dcc pathway genes. rdpA(MH) was 99% identical to rdpA(MC1), an (R)-dichlorprop/alpha-ketoglutarate dioxygenase from Delftia acidovorans MC1, which is evidence for a recent gene exchange between Alpha- and Betaproteobacteria. Conversely, DccA(I) and DccA(II) did not group within the known chlorocatechol 1,2-dioxygenases, but formed a separate branch in clustering analysis. This suggests a different reservoir and reduced transfer for the genes of the modified ortho-cleavage pathway in Alphaproteobacteria compared with the ones in Beta- and Gammaproteobacteria.

Gene expression variation and expression quantitative trait mapping of human chromosome 21 genes.

Inter-individual differences in gene expression are likely to account for an important fraction of phenotypic differences, including susceptibility to common disorders. Recent studies have shown extensive variation in gene expression levels in humans and other organisms, and that a fraction of this variation is under genetic control. We investigated the patterns of gene expression variation in a 25 Mb region of human chromosome 21, which has been associated with many Down syndrome (DS) phenotypes. Taqman real-time PCR was used to measure expression variation of 41 genes in lymphoblastoid cells of 40 unrelated individuals. For 25 genes found to be differentially expressed, additional analysis was performed in 10 CEPH families to determine heritabilities and map loci harboring regulatory variation. Seventy-six percent of the differentially expressed genes had significant heritabilities, and genomewide linkage analysis led to the identification of significant eQTLs for nine genes. Most eQTLs were in trans, with the best result (P=7.46 x 10(-8)) obtained for TMEM1 on chromosome 12q24.33. A cis-eQTL identified for CCT8 was validated by performing an association study in 60 individuals from the HapMap project. SNP rs965951 located within CCT8 was found to be significantly associated with its expression levels (P=2.5 x 10(-5)) confirming cis-regulatory variation. The results of our study provide a representative view of expression variation of chromosome 21 genes, identify loci involved in their regulation and suggest that genes, for which expression differences are significantly larger than 1.5-fold in control samples, are unlikely to be involved in DS-phenotypes present in all affected individuals.

Identification of estrogen-responsive DNA sequences by transient expression experiments in a human breast cancer cell line.

The expression of a hybrid gene formed by the promoter region of the Xenopus laevis vitellogenin gene B1 and the CAT coding region is regulated by estrogen when the gene is transfected into hormone-responsive MCF-7 cells. Furthermore, the 5' flanking region of the gene B1 alone can confer inducibility to heterologous promoters, although to a varying extent depending on the promoter used. Deletion mapping of he vitellogenin hormone-responsive sequences revealed that a 13 bp element 5'-AGTCACTGTGACC-3' at position -334 is essential for estrogen inducibility. We have shown previously that this 13 bp element is present upstream of several liver-specific estrogen-inducible genes.