Functional analysis of pyochelin-/enantiopyochelin-related genes from a pathogenicity island of Pseudomonas aeruginosa strain PA14.

Genomic islands are foreign DNA blocks inserted in so-called regions of genomic plasticity (RGP). Depending on their gene content, they are classified as pathogenicity, symbiosis, metabolic, fitness or resistance islands, although a detailed functional analysis is often lacking. Here we focused on a 34-kb pathogenicity island of Pseudomonas aeruginosa PA14 (PA14GI-6), which is inserted at RGP5 and carries genes related to those for pyochelin/enantiopyochelin biosynthesis. These enantiomeric siderophores of P. aeruginosa and certain strains of Pseudomonas protegens are assembled by a thiotemplate mechanism from salicylate and two molecules of cysteine. The biochemical function of several proteins encoded by PA14GI-6 was investigated by a series of complementation analyses using mutants affected in potential homologs. We found that PA14_54940 codes for a bifunctional salicylate synthase/salicyl-AMP ligase (for generation and activation of salicylate), that PA14_54930 specifies a dihydroaeruginoic acid (Dha) synthetase (for coupling salicylate with a cysteine-derived thiazoline ring), that PA14_54910 produces a type II thioesterase (for quality control), and that PA14_54880 encodes a serine O-acetyltransferase (for increased cysteine availability). The structure of the PA14GI-6-specified metabolite was determined by mass spectrometry, thin-layer chromatography, and HPLC as (R)-Dha, an iron chelator with antibacterial, antifungal and antitumor activity. The conservation of this genomic island in many clinical and environmental P. aeruginosa isolates of different geographical origin suggests that the ability for Dha production may confer a selective advantage to its host.

A non-circadian role for clock-genes in sleep homeostasis: a strain comparison.

BACKGROUND: We have previously reported that the expression of circadian clock-genes increases in the cerebral cortex after sleep deprivation (SD) and that the sleep rebound following SD is attenuated in mice deficient for one or more clock-genes. We hypothesized that besides generating circadian rhythms, clock-genes also play a role in the homeostatic regulation of sleep. Here we follow the time course of the forebrain changes in the expression of the clock-genes period (per)-1, per2, and of the clock-controlled gene albumin D-binding protein (dbp) during a 6 h SD and subsequent recovery sleep in three inbred strains of mice for which the homeostatic sleep rebound following SD differs. We reasoned that if clock genes are functionally implicated in sleep homeostasis then the SD-induced changes in gene expression should vary according to the genotypic differences in the sleep rebound. RESULTS: In all three strains per expression was increased when animals were kept awake but the rate of increase during the SD as well as the relative increase in per after 6 h SD were highest in the strain for which the sleep rebound was smallest; i.e., DBA/2J (D2). Moreover, whereas in the other two strains per1 and per2 reverted to control levels with recovery sleep, per2 expression specifically, remained elevated in D2 mice. dbp expression increased during the light period both during baseline and during SD although levels were reduced during the latter condition compared to baseline. In contrast to per2, dbp expression reverted to control levels with recovery sleep in D2 only, whereas in the two other strains expression remained decreased. CONCLUSION: These findings support and extend our previous findings that clock genes in the forebrain are implicated in the homeostatic regulation of sleep and suggest that sustained, high levels of per2 expression may negatively impact recovery sleep.

Analysis of stop-gain and frameshift variants in human innate immunity genes.

Loss-of-function variants in innate immunity genes are associated with Mendelian disorders in the form of primary immunodeficiencies. Recent resequencing projects report that stop-gains and frameshifts are collectively prevalent in humans and could be responsible for some of the inter-individual variability in innate immune response. Current computational approaches evaluating loss-of-function in genes carrying these variants rely on gene-level characteristics such as evolutionary conservation and functional redundancy across the genome. However, innate immunity genes represent a particular case because they are more likely to be under positive selection and duplicated. To create a ranking of severity that would be applicable to innate immunity genes we evaluated 17,764 stop-gain and 13,915 frameshift variants from the NHLBI Exome Sequencing Project and 1,000 Genomes Project. Sequence-based features such as loss of functional domains, isoform-specific truncation and nonsense-mediated decay were found to correlate with variant allele frequency and validated with gene expression data. We integrated these features in a Bayesian classification scheme and benchmarked its use in predicting pathogenic variants against Online Mendelian Inheritance in Man (OMIM) disease stop-gains and frameshifts. The classification scheme was applied in the assessment of 335 stop-gains and 236 frameshifts affecting 227 interferon-stimulated genes. The sequence-based score ranks variants in innate immunity genes according to their potential to cause disease, and complements existing gene-based pathogenicity scores. Specifically, the sequence-based score improves measurement of functional gene impairment, discriminates across different variants in a given gene and appears particularly useful for analysis of less conserved genes.

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.

Invasion and fixation of sex-reversal genes.

We simulated a meta-population with random dispersal among demes but local mating within demes to investigate conditions under which a dominant female-determining gene W, with no individual selection advantage, can invade and become fixed in females, changing the population from male to female heterogamety. Starting with one mutant W in a single deme, the interaction of sex ratio selection and random genetic drift causes W to be fixed among females more often than a comparable neutral mutation with no influence on sex determination, even when YY males have slightly reduced viability. Meta-population structure and interdeme selection can also favour the fixation of W. The reverse transition from female to male heterogamety can also occur with higher probability than for a comparable neutral mutation. These results help to explain the involvement of sex-determining genes in the evolution of sex chromosomes and in sexual selection and speciation.

CYP3A5 and ABCB1 genes influence blood pressure and response to treatment, and their effect is modified by salt.

The permeability-glycoprotein efflux-transporter encoded by the multidrug resistance 1 (ABCB1) gene and the cytochromes P450 3A4/5 encoded by the CYP3A4/5 genes are known to interact in the transport and metabolism of many drugs. Recent data have shown that the CYP3A5 genotypes influence blood pressure and that permeability-glycoprotein activity might influence the activity of the renin-angiotensin system. Hence, these 2 genes may contribute to blood pressure regulation in humans. We analyzed the association of variants of the ABCB1 and CYP3A5 genes with ambulatory blood pressure, plasma renin activity, plasma aldosterone, endogenous lithium clearance, and blood pressure response to treatment in 72 families (373 individuals; 55% women; mean age: 46 years) of East African descent. The ABCB1 and CYP3A5 genes interact with urinary sodium excretion in their effect on ambulatory blood pressure (daytime systolic: P=0.05; nighttime systolic and diastolic: P<0.01), suggesting a gene-gene-environment interaction. The combined action of these genes is also associated with postproximal tubular sodium reabsorption, plasma renin activity, plasma aldosterone, and with an altered blood pressure response to the angiotensin-converting enzyme inhibitor lisinopril (P<0.05). This is the first reported association of the ABCB1 gene with blood pressure in humans and demonstration that genes encoding for proteins metabolizing and transporting drugs and endogenous substrates contribute to blood pressure regulation.

Essential Bacillus subtilis genes.

To estimate the minimal gene set required to sustain bacterial life in nutritious conditions, we carried out a systematic inactivation of Bacillus subtilis genes. Among approximately 4,100 genes of the organism, only 192 were shown to be indispensable by this or previous work. Another 79 genes were predicted to be essential. The vast majority of essential genes were categorized in relatively few domains of cell metabolism, with about half involved in information processing, one-fifth involved in the synthesis of cell envelope and the determination of cell shape and division, and one-tenth related to cell energetics. Only 4% of essential genes encode unknown functions. Most essential genes are present throughout a wide range of Bacteria, and almost 70% can also be found in Archaea and Eucarya. However, essential genes related to cell envelope, shape, division, and respiration tend to be lost from bacteria with small genomes. Unexpectedly, most genes involved in the Embden-Meyerhof-Parnas pathway are essential. Identification of unknown and unexpected essential genes opens research avenues to better understanding of processes that sustain bacterial life.

Promoter and transcription analysis of penicillin-binding protein genes in Streptococcus gordonii.

An optimally cross-linked peptidoglycan requires both transglycosylation and transpeptidation, provided by class A and class B penicillin-binding proteins (PBPs). Streptococcus gordonii possesses three class A PBPs (PBPs 1A, 1B, and 2A) and two class B PBPs (PBPs 2B and 2X) that are important for penicillin resistance. High-level resistance (MIC, > or =2 microg/ml) requires mutations in class B PBPs. However, although unmutated, class A PBPs are critical to facilitate resistance development (M. Haenni and P. Moreillon, Antimicrob. Agents Chemother. 50:4053-4061, 2006). Thus, their overexpression might be important to sustain the drug. Here, we determined the promoter regions of the S. gordonii PBPs and compared them to those of other streptococci. The extended -10 box was highly conserved and complied with a sigma(A)-type promoter consensus sequence. In contrast, the -35 box was poorly conserved, leaving the possibility of differential PBP regulation. Gene expression in a penicillin-susceptible parent (MIC, 0.008 microg/ml) and a high-level-resistant mutant (MIC, 2 microg/ml) was monitored using luciferase fusions. In the absence of penicillin, all PBPs were constitutively expressed, but their expression was globally increased (1.5 to 2 times) in the resistant mutant. In the presence of penicillin, class A PBPs were specifically overexpressed both in the parent (PBP 2A) and in the resistant mutant (PBPs 1A and 2A). By increasing transglycosylation, class A PBPs could promote peptidoglycan stability when transpeptidase is inhibited by penicillin. Since penicillin-related induction of class A PBPs occurred in both susceptible and resistant cells, such a mutation-independent facilitating mechanism could be operative at each step of resistance development.

Messenger RNA expression of transporter and ion channel genes in undifferentiated and differentiated Caco-2 cells compared to human intestines.

PURPOSE: The purpose of this work was to study the influence of cell differentiation on the mRNA expression of transporters and channels in Caco-2 cells and to assess Caco-2 cells as a model for carrier-mediated drug transport in the intestines. METHOD: Gene mRNA expression was measured using a custom-designed microarray chip with 750 deoxyoligonucleotide probes (70mers). Each oligomer was printed four times on poly-lysine-coated glass slides. Expression profiles were expressed as ratio values between fluorescence intensities of Cy3 and Cy5 dye-labeled cDNA derived from poly(A) + RNA samples of Caco-2 cells and total RNA of human intestines. RESULTS: Significant differences in the mRNA expression profile of transporters and channels were observed upon differentiation of Caco-2 cells from 5 days to 2 weeks in culture, including changes for MAT8, S-protein, and Nramp2. Comparing Caco-2 cells of different passage number revealed few changes in mRNAs except for GLUT3, which was down-regulated 2.4-fold within 13 passage numbers. Caco-2 cells had a similar expression profile when either cultured in flasks or on filters but differed more strongly from human small and large intestine, regardless of the differentiation state of Caco-2 cells. Expression of several genes highly transcribed in small or large intestines differed fourfold or more in Caco-2 cells. CONCLUSIONS: Although Caco-2 cells have proven a suitable model for studying carrier-mediated transport in human intestines, the expression of specific transporter and ion channel genes may differ substantially.

Classification of human astrocytic gliomas on the basis of gene expression: a correlated group of genes with angiogenic activity emerges as a strong predictor of subtypes.

The development of targeted treatment strategies adapted to individual patients requires identification of the different tumor classes according to their biology and prognosis. We focus here on the molecular aspects underlying these differences, in terms of sets of genes that control pathogenesis of the different subtypes of astrocytic glioma. By performing cDNA-array analysis of 53 patient biopsies, comprising low-grade astrocytoma, secondary glioblastoma (respective recurrent high-grade tumors), and newly diagnosed primary glioblastoma, we demonstrate that human gliomas can be differentiated according to their gene expression. We found that low-grade astrocytoma have the most specific and similar expression profiles, whereas primary glioblastoma exhibit much larger variation between tumors. Secondary glioblastoma display features of both other groups. We identified several sets of genes with relatively highly correlated expression within groups that: (a). can be associated with specific biological functions; and (b). effectively differentiate tumor class. One prominent gene cluster discriminating primary versus nonprimary glioblastoma comprises mostly genes involved in angiogenesis, including VEGF fms-related tyrosine kinase 1 but also IGFBP2, that has not yet been directly linked to angiogenesis. In situ hybridization demonstrating coexpression of IGFBP2 and VEGF in pseudopalisading cells surrounding tumor necrosis provided further evidence for a possible involvement of IGFBP2 in angiogenesis. The separating groups of genes were found by the unsupervised coupled two-way clustering method, and their classification power was validated by a supervised construction of a nearly perfect glioma classifier.

Substitution rates at neutral genes depend on population size under fluctuating demography and overlapping generations.

It is widely accepted that the rate of evolution (substitution rate) at neutral genes is unaffected by population size fluctuations. This result has implications for the analysis of genetic data in population genetics and phylogenetics, and provides, in particular, a justification for the concept of the molecular clock. Here, we show that the substitution rate at neutral genes does depend on population size fluctuations in the presence of overlapping generations. As both population size fluctuations and overlapping generations are expected to be the norm rather than the exception in natural populations, this observation may be relevant for understanding variation in substitution rates within and between lineages.

Selecting control genes for RT-QPCR using public microarray data.

Background: Gene expression analysis has emerged as a major biological research area, with real-time quantitative reverse transcription PCR (RT-QPCR) being one of the most accurate and widely used techniques for expression profiling of selected genes. In order to obtain results that are comparable across assays, a stable normalization strategy is required. In general, the normalization of PCR measurements between different samples uses one to several control genes (e. g. housekeeping genes), from which a baseline reference level is constructed. Thus, the choice of the control genes is of utmost importance, yet there is not a generally accepted standard technique for screening a large number of candidates and identifying the best ones. Results: We propose a novel approach for scoring and ranking candidate genes for their suitability as control genes. Our approach relies on publicly available microarray data and allows the combination of multiple data sets originating from different platforms and/or representing different pathologies. The use of microarray data allows the screening of tens of thousands of genes, producing very comprehensive lists of candidates. We also provide two lists of candidate control genes: one which is breast cancer-specific and one with more general applicability. Two genes from the breast cancer list which had not been previously used as control genes are identified and validated by RT-QPCR. Open source R functions are available at http://www.isrec.isb-sib.ch/similar to vpopovic/research/ Conclusion: We proposed a new method for identifying candidate control genes for RT-QPCR which was able to rank thousands of genes according to some predefined suitability criteria and we applied it to the case of breast cancer. We also empirically showed that translating the results from microarray to PCR platform was achievable.

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.

Evolution of vitellogenin genes: comparative analysis of the nucleotide sequences downstream of the transcription initiation site of four Xenopus laevis and one chicken gene.

Electron microscopic analysis of heteroduplexes between the most distantly related Xenopus vitellogenin genes (A genes X B genes) has revealed the distribution of homologous regions that have been preferentially conserved after the duplication events that gave rise to the multigene family in Xenopus laevis. DNA sequence analysis was limited to the region downstream of the transcription initiation site of the Xenopus genes A1, B1 and B2 and a comparison with the Xenopus A2 and the major chicken vitellogenin gene is presented. Within the coding regions of the first three exons, nucleotide substitutions resulting in amino acid changes accumulate at a rate similar to that observed in globin genes. This suggests that the duplication event which led to the formation of the A and B ancestral genes in Xenopus laevis occurred about 150 million years ago. Homologous exons of the A1-A2 and B1-B2 gene pairs, which formed about 30 million years ago, show a quite similar sequence divergence. In contrast, A1-A2 homologous introns seem to have evolved much faster than their B1-B2 counterparts.