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.

Reverse transcription quantitative real-time polymerase chain reaction reference genes in the spared nerve injury model of neuropathic pain : validation and literature search

La douleur neuropathique est définie comme une douleur causée par une lésion du système nerveux somato-sensoriel. Elle se caractérise par des douleurs exagérées, spontanées, ou déclenchées par des stimuli normalement non douloureux (allodynie) ou douloureux (hyperalgésie). Bien qu'elle concerne 7% de la population, ses mécanismes biologiques ne sont pas encore élucidés. L'étude des variations d'expressions géniques dans les tissus-clés des voies sensorielles (notamment le ganglion spinal et la corne dorsale de la moelle épinière) à différents moments après une lésion nerveuse périphérique permettrait de mettre en évidence de nouvelles cibles thérapeutiques. Elles se détectent de manière sensible par reverse transcription quantitative real-time polymerase chain reaction (RT- qPCR). Pour garantir des résultats fiables, des guidelines ont récemment recommandé la validation des gènes de référence utilisés pour la normalisation des données ("Minimum information for publication of quantitative real-time PCR experiments", Bustin et al 2009). Après recherche dans la littérature des gènes de référence fréquemment utilisés dans notre modèle de douleur neuropathique périphérique SNI (spared nerve injury) et dans le tissu nerveux en général, nous avons établi une liste de potentiels bons candidats: Actin beta (Actb), Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), ribosomal proteins 18S (18S), L13a (RPL13a) et L29 (RPL29), hypoxanthine phosphoribosyltransferase 1 (HPRT1) et hydroxymethyl-bilane synthase (HMBS). Nous avons évalué la stabilité d'expression de ces gènes dans le ganglion spinal et dans la corne dorsale à différents moments après la lésion nerveuse (SNI) en calculant des coefficients de variation et utilisant l'algorithme geNorm qui compare les niveaux d'expression entre les différents candidats et détermine la paire de gènes restante la plus stable. Il a aussi été possible de classer les gènes selon leur stabilité et d'identifier le nombre de gènes nécessaires pour une normalisation la plus précise. Les gènes les plus cités comme référence dans le modèle SNI ont été GAPDH, HMBS, Actb, HPRT1 et 18S. Seuls HPRT1 and 18S ont été précédemment validés dans des arrays de RT-qPCR. Dans notre étude, tous les gènes testés dans le ganglion spinal et dans la corne dorsale satisfont au critère de stabilité exprimé par une M-value inférieure à 1. Par contre avec un coefficient de variation (CV) supérieur à 50% dans le ganglion spinal, 18S ne peut être retenu. La paire de gènes la plus stable dans le ganglion spinal est HPRT1 et Actb et dans la corne dorsale il s'agit de RPL29 et RPL13a. L'utilisation de 2 gènes de référence stables suffit pour une normalisation fiable. Nous avons donc classé et validé Actb, RPL29, RPL13a, HMBS, GAPDH, HPRT1 et 18S comme gènes de référence utilisables dans la corne dorsale pour le modèle SNI chez le rat. Dans le ganglion spinal 18S n'a pas rempli nos critères. Nous avons aussi déterminé que la combinaison de deux gènes de référence stables suffit pour une normalisation précise. Les variations d'expression génique de potentiels gènes d'intérêts dans des conditions expérimentales identiques (SNI, tissu et timepoints post SNI) vont pouvoir se mesurer sur la base d'une normalisation fiable. Non seulement il sera possible d'identifier des régulations potentiellement importantes dans la genèse de la douleur neuropathique mais aussi d'observer les différents phénotypes évoluant au cours du temps après lésion nerveuse.

Sleep loss reduces the DNA-binding of BMAL1, CLOCK, and NPAS2 to specific clock genes in the mouse cerebral cortex.

We have previously demonstrated that clock genes contribute to the homeostatic aspect of sleep regulation. Indeed, mutations in some clock genes modify the markers of sleep homeostasis and an increase in homeostatic sleep drive alters clock gene expression in the forebrain. Here, we investigate a possible mechanism by which sleep deprivation (SD) could alter clock gene expression by quantifying DNA-binding of the core-clock transcription factors CLOCK, NPAS2, and BMAL1 to the cis-regulatory sequences of target clock genes in mice. Using chromatin immunoprecipitation (ChIP), we first showed that, as reported for the liver, DNA-binding of CLOCK and BMAL1 to target clock genes changes in function of time-of-day in the cerebral cortex. Tissue extracts were collected at ZT0 (light onset), -6, -12, and -18, and DNA enrichment of E-box or E'-box containing sequences was measured by qPCR. CLOCK and BMAL1 binding to Cry1, Dbp, Per1, and Per2 depended on time-of-day, with maximum values reached at around ZT6. We then observed that SD, performed between ZT0 and -6, significantly decreased DNA-binding of CLOCK and BMAL1 to Dbp, consistent with the observed decrease in Dbp mRNA levels after SD. The DNA-binding of NPAS2 and BMAL1 to Per2 was also decreased by SD, although SD is known to increase Per2 expression in the cortex. DNA-binding to Per1 and Cry1 was not affected by SD. Our results show that the sleep-wake history can affect the clock molecular machinery directly at the level of chromatin binding thereby altering the cortical expression of Dbp and Per2 and likely other targets. Although the precise dynamics of the relationship between DNA-binding and mRNA expression, especially for Per2, remains elusive, the results also suggest that part of the reported circadian changes in DNA-binding of core clock components in tissues peripheral to the suprachiasmatic nuclei could, in fact, be sleep-wake driven.

Identification of corticosteroid-regulated genes in cardiomyocytes by serial analysis of gene expression.

Corticosteroids (aldosterone, cortisol/corticosterone) exert direct functional effects on cardiomyocytes. However, gene networks activated by corticosteroids in cardiomyocytes, as well as the involvement of the mineralocorticoid receptor (MR) vs the glucocorticoid receptor (GR) in these effects, remain largely unknown. Here we characterized the corticosteroid-dependent transcriptome in primary culture of neonatal mouse cardiomyocytes treated with 10(-6) M aldosterone, a concentration predicted to occupy both MR and GR. Serial analysis of gene expression revealed 101 aldosterone-regulated genes. The MR/GR specificity was characterized for one regulated transcript, namely ecto-ADP-ribosyltransferase-3 (Art3). Using cardiomyocytes from GR(null/null) or MR(null/null) mice we demonstrate that in GR(null/null) cardiomyocytes the response is abrogated, but it is fully maintained in MR(null/null) cardiomyocytes. We conclude that Art3 expression is regulated exclusively via the GR. Our study identifies a new set of corticosteroid-regulated genes in cardiomyocytes and demonstrates a new approach to studying the selectivity of MR- vs GR-dependent effects.