cAMP prevents the glucose-mediated stimulation of GLUT2 gene transcription in hepatocytes.

Glucose homoeostasis necessitates the presence in the liver of the high Km glucose transporter GLUT2. In hepatocytes, we and others have demonstrated that glucose stimulates GLUT2 gene expression in vivo and in vitro. This effect is transcriptionally regulated and requires glucose metabolism within the hepatocytes. In this report, we further characterized the cis-elements of the murine GLUT2 promoter, which confers glucose responsiveness on a reporter gene coding the chloramphenicol acetyl transferase (CAT) gene. 5'-Deletions of the murine GLUT2 promoter linked to the CAT reporter gene were transfected into a GLUT2 expressing hepatoma cell line (mhAT3F) and into primary cultured rat hepatocytes, and subsequently incubated at low and high glucose concentrations. Glucose stimulates gene transcription in a manner similar to that observed for the endogenous GLUT2 mRNA in both cell types; the -1308 to -212 bp region of the promoter contains the glucose-responsive elements. Furthermore, the -1308 to -338 bp region of the promoter contains repressor elements when tested in an heterologous thymidine kinase promoter. The glucose-induced GLUT2 mRNA accumulation was decreased by dibutyryl-cAMP both in mhAT3F cells and in primary hepatocytes. A putative cAMP-responsive element (CRE) is localized at the -1074/-1068 bp region of the promoter. The inhibitory effect of cAMP on GLUT2 gene expression was observed in hepatocytes transfected with constructs containing this CRE (-1308/+49 bp fragment), as well as with constructs not containing the consensus CRE (-312/+49 bp fragment). This suggests that the inhibitory effect of cAMP is not mediated by the putative binding site located in the repressor fragment of the GLUT2 promoter. Taken together, these data demonstrate that the elements conferring glucose and cAMP responsiveness on the GLUT2 gene are located within the -312/+49 region of the promoter.

Biased gene conversion and GC-content evolution in the coding sequences of reptiles and vertebrates.

Mammalian and avian genomes are characterized by a substantial spatial heterogeneity of GC-content, which is often interpreted as reflecting the effect of local GC-biased gene conversion (gBGC), a meiotic repair bias that favors G and C over A and T alleles in high-recombining genomic regions. Surprisingly, the first fully sequenced nonavian sauropsid (i.e., reptile), the green anole Anolis carolinensis, revealed a highly homogeneous genomic GC-content landscape, suggesting the possibility that gBGC might not be at work in this lineage. Here, we analyze GC-content evolution at third-codon positions (GC3) in 44 vertebrates species, including eight newly sequenced transcriptomes, with a specific focus on nonavian sauropsids. We report that reptiles, including the green anole, have a genome-wide distribution of GC3 similar to that of mammals and birds, and we infer a strong GC3-heterogeneity to be already present in the tetrapod ancestor. We further show that the dynamic of coding sequence GC-content is largely governed by karyotypic features in vertebrates, notably in the green anole, in agreement with the gBGC hypothesis. The discrepancy between third-codon positions and noncoding DNA regarding GC-content dynamics in the green anole could not be explained by the activity of transposable elements or selection on codon usage. This analysis highlights the unique value of third-codon positions as an insertion/deletion-free marker of nucleotide substitution biases that ultimately affect the evolution of proteins.

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.

Differential gene expression in response to mechanical wounding and insect feeding in Arabidopsis.

Wounding in multicellular eukaryotes results in marked changes in gene expression that contribute to tissue defense and repair. Using a cDNA microarray technique, we analyzed the timing, dynamics, and regulation of the expression of 150 genes in mechanically wounded leaves of Arabidopsis. Temporal accumulation of a group of transcripts was correlated with the appearance of oxylipin signals of the jasmonate family. Analysis of the coronatine-insensitive coi1-1 Arabidopsis mutant that is also insensitive to jasmonate allowed us to identify a large number of COI1-dependent and COI1-independent wound-inducible genes. Water stress was found to contribute to the regulation of an unexpectedly large fraction of these genes. Comparing the results of mechanical wounding with damage by feeding larvae of the cabbage butterfly (Pieris rapae) resulted in very different transcript profiles. One gene was specifically induced by insect feeding but not by wounding; moreover, there was a relative lack of water stress-induced gene expression during insect feeding. These results help reveal a feeding strategy of P. rapae that may minimize the activation of a subset of water stress-inducible, defense-related genes.

Modular analysis of gene expression data with R.

SUMMARY: Large sets of data, such as expression profiles from many samples, require analytic tools to reduce their complexity. The Iterative Signature Algorithm (ISA) is a biclustering algorithm. It was designed to decompose a large set of data into so-called 'modules'. In the context of gene expression data, these modules consist of subsets of genes that exhibit a coherent expression profile only over a subset of microarray experiments. Genes and arrays may be attributed to multiple modules and the level of required coherence can be varied resulting in different 'resolutions' of the modular mapping. In this short note, we introduce two BioConductor software packages written in GNU R: The isa2 package includes an optimized implementation of the ISA and the eisa package provides a convenient interface to run the ISA, visualize its output and put the biclusters into biological context. Potential users of these packages are all R and BioConductor users dealing with tabular (e.g. gene expression) data. AVAILABILITY: http://www.unil.ch/cbg/ISA CONTACT: sven.bergmann@unil.ch

Variant within the promoter region of the CHRNA3 gene associated with FTN dependence is not related to self-reported willingness to quit smoking.

INTRODUCTION: Common variation in the CHRNA5-CHRNA3-CHRNB4 gene region is robustly associated with smoking quantity. Conversely, the association between one of the most significant single nucleotide polymorphisms (SNPs; rs1051730 within the CHRNA3 gene) with perceived difficulty or willingness to quit smoking among current smokers is unknown. METHODS: Cross-sectional study including current smokers, 502 women, and 552 men. Heaviness of smoking index (HSI), difficulty, attempting, and intention to quit smoking were assessed by questionnaire. RESULTS: The rs1051730 SNP was associated with increased HSI (age, gender, and education-adjusted mean ± SE: 2.6 ± 0.1, 2.2 ± 0.1, and 2.0 ± 0.1 for AA, AG, and GG genotypes, respectively, p < .01). Multivariate logistic regression adjusting for gender, age, education, leisure-time physical activity, and personal history of cardiovascular or lung disease showed rs1051730 to be associated with higher smoking dependence (odds ratio [OR] and 95% CI for each additional A-allele: 1.38 [1.11-1.72] for smoking more than 20 cigarette equivalents/day; 1.31 [1.00-1.71] for an HSI ≥5 and 1.32 [1.05-1.65] for smoking 5 min after waking up) and borderline associated with difficulty to quit (OR = 1.29 [0.98-1.70]), but this relationship was no longer significant after adjusting for nicotine dependence. Also, no relationship was found with willingness (OR = 1.03 [0.85-1.26]), attempt (OR = 1.00 [0.83-1.20]), or preparation (OR = 0.95 [0.38-2.38]) to quit. Similar findings were obtained for other SNPs, but their effect on nicotine dependence was no longer significant after adjusting for rs1051730. Conclusions: These data confirm the effect of rs1051730 on nicotine dependence but failed to find any relationship with difficulty, willingness, and motivation to quit.

Chromosomal gene movements reflect the recent origin and biology of therian sex chromosomes

Mammalian sex chromosomes stem from ancestral autosomes and have substantially differentiated. It was shown that X-linked genes have generated duplicate intronless gene copies (retrogenes) on autosomes due to this differentiation. However, the precise driving forces for this out-of-X gene "movement" and its evolutionary onset are not known. Based on expression analyses of male germ-cell populations, we here substantiate and extend the hypothesis that autosomal retrogenes functionally compensate for the silencing of their X-linked housekeeping parental genes during, but also after, male meiotic sex chromosome inactivation (MSCI). Thus, sexually antagonistic forces have not played a major role for the selective fixation of X-derived gene copies in mammals. Our dating analyses reveal that although retrogenes were produced ever since the common mammalian ancestor, selectively driven retrogene export from the X only started later, on the placental mammal (eutherian) and marsupial (metatherian) lineages, respectively. Together, these observations suggest that chromosome-wide MSCI emerged close to the eutherian-marsupial split approximately 180 million years ago. Given that MSCI probably reflects the spread of the recombination barrier between the X and Y, crucial for their differentiation, our data imply that these chromosomes became more widely differentiated only late in the therian ancestor, well after the divergence of the monotreme lineage. Thus, our study also provides strong independent support for the recent notion that our sex chromosomes emerged, not in the common ancestor of all mammals, but rather in the therian ancestor, and therefore are much younger than previously thought

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.

SNCA, MAPT, and GSK3B in Parkinson disease: a gene-gene interaction study.

Background and purpose:  Recent evidence suggests that variation in the SNCA, MAPT, and GSK3B genes interacts in affecting risk for Parkinson disease (PD). In the current study, we attempt to validate previously published findings, evaluating gene-gene interactions between SNCA, MAPT, and GSK3B in association with PD. Methods:  Three Caucasian PD patient-control series from the United States, Ireland, and Norway (combined n = 1020 patients and 1095 controls) were genotyped for SNCA rs356219, MAPT H1/H2-discriminating SNP rs1052553, and GSK3B rs334558 and rs6438552. Results:  Our findings indicate that as previously reported, the SNCA rs356219-G allele and MAPT rs1052553 (H1 haplotype) were both associated with an increased risk of PD, whilst contrary to previous reports, GSK3B variants were not. No pair-wise interaction was observed between SNCA, MAPT, and GSK3B; the risk effects of SNCA rs356219-G and MAPT rs1052553-H1 were seen in a similar manner across genotypes of other variants, with no evidence suggesting synergistic, antagonistic, or deferential effects. Conclusions:  In the Caucasian patient-control series examined, risk for PD was influenced by variation in SNCA and MAPT but not GSK3B. Additionally, those three genes did not interact in determining disease risk.

Fine-scale spatial genetic structure and gene dispersal in Silene latifolia.

Plants are sessile organisms, often characterized by limited dispersal. Seeds and pollen are the critical stages for gene flow. Here we investigate spatial genetic structure, gene dispersal and the relative contribution of pollen vs seed in the movement of genes in a stable metapopulation of the white campion Silene latifolia within its native range. This short-lived perennial plant is dioecious, has gravity-dispersed seeds and moth-mediated pollination. Direct measures of pollen dispersal suggested that large populations receive more pollen than small isolated populations and that most gene flow occurs within tens of meters. However, these studies were performed in the newly colonized range (North America) where the specialist pollinator is absent. In the native range (Europe), gene dispersal could fall on a different spatial scale. We genotyped 258 individuals from large and small (15) subpopulations along a 60 km, elongated metapopulation in Europe using six highly variable microsatellite markers, two X-linked and four autosomal. We found substantial genetic differentiation among subpopulations (global F(ST)=0.11) and a general pattern of isolation by distance over the whole sampled area. Spatial autocorrelation revealed high relatedness among neighboring individuals over hundreds of meters. Estimates of gene dispersal revealed gene flow at the scale of tens of meters (5-30 m), similar to the newly colonized range. Contrary to expectations, estimates of dispersal based on X and autosomal markers showed very similar ranges, suggesting similar levels of pollen and seed dispersal. This may be explained by stochastic events of extensive seed dispersal in this area and limited pollen dispersal.

Implication of DNA methylation and PAX5 factor in the transcriptional regulation of hTERT, the human telomerase reverse transcriptase gene

RESUME La télomérase est une enzyme dite "d'immortalité" qui permet aux cellules de maintenir la longueur de leurs télomères, ce qui confère une capacité de réplication illimitée aux cellules reproductrices et cancéreuses. A l'inverse, les cellules somatiques normales, qui n'expriment pas la télomérase, ont une capacité de réplication limitée. La sous-unité catalytique de la télomérase, hTERT, est définie comme le facteur limitant l'activité télomérasique. Entre activateurs et répresseurs, le rôle de la méthylation de l'ADN et de l'acétylation des histones, de nombreux modèles ont été suggérés. La découverte de l'implication de CTCF dans la régulation transcriptionnelle de hTERT explique en partie le mécanisme de répression de la télomérase dans la plupart des cellules somatiques et sa réactivation dans les cellules tumorales. Dans les cellules télomérase-positives, l'activité inhibitrice de CTCF est bloquée par un mécanisme dépendent ou non de la méthylation. Dans la plupart des carcinomes, une hyperméthylation de la région 5' de hTERT bloque l'effet inhibiteur de CTCF, alors qu'une petite région hypométhylée permet un faible niveau de transcription du gène. Nous avons démontré que la protéine MBD2 se lie spécifiquement sur la région 5' méthylée de hTERT dans différentes lignées cellulaires et qu'elle est impliquée dans la répression partielle de la transcription de hTERT dans les cellules tumorales méthylées. Par contre, nous avons montré que dans les lymphocytes B normaux et néoplasiques, la régulation de hTERT est indépendante de la méthylation. Dans ces cellules, le facteur PAX5 se lie sur la région 5' de hTERT en aval du site d'initiation de la traduction (ATG). L'expression exogène de PAX5 dans les cellules télomérase-négatives active la transcription de hTERT, alors que la répression de PAX5 dans les cellules lymphomateuses inhibe la transcription du gène. PAX5 est donc directement impliqué dans l'activation de l'expression de hTERT dans les lymphocytes B exprimant la télomérase. Ces résultats révèlent des différences entre les niveaux de méthylation de hTERT dans les cellules de carcinomes et les lymphocytes B exprimant la télomérase. La méthylation de hTERT en tant que biomarqueur de cancer a été évaluée, puis appliquée à la détection de métastases. Nous avons ainsi montré que la méthylation de hTERT est positivement corrélée au diagnostic cytologique dans les liquides céphalorachidiens. Nos résultats conduisent à un modèle de régulation de hTERT, qui aide à comprendre comment la transcription de ce gène est régulée par CTCF, avec un mécanisme lié ou non à la méthylation du gène hTERT. La méthylation de hTERT s'est aussi révélée être un nouveau et prometteur biomarqueur de cancer. SUMMARY Human telomerase is an "immortalizing" enzyme that enables cells to maintain telomere length, allowing unlimited replicative capacity to reproductive and cancer cells. Conversely, normal somatic cells that do not express telomerase have a finite replicative capacity. The catalytic subunit of telomerase, hTERT, is defined as the limiting factor for telomerase activity. Between activators and repressors, and the role of DNA methylation and histone acetylation, an abundance of hTERT regulatory models have been suggested. The discovery of the implication of CTCF in the transcriptional regulation of hTERT in part explained the mechanism of silencing of telomerase in most somatic cells and its reactivation in neoplastic cells. In telomerase-positive cells, the inhibitory activity of CTCF is blocked by methylation-dependent and -independent mechanisms. In most carcinoma cells, hypermethylation of the hTERT 5' region has been shown to block the inhibitory effect of CTCF, while a short hypomethylated region allows a low transcription level of the gene. We have demonstrated that MBD2 protein specifically binds the methylated 5' region of hTERT in different cell lines and is therefore involved in the partial repression of hTERT transcription in methylated tumor cells. In contrast, we have shown that in normal and neoplastic B cells, hTERT regulation is methylation-independent. The PAX5 factor has been shown to bind to the hTERT 5'region downstream of the ATG translational start site. Ectopic expression of PAX5 in telomerase-negative cells or repression of PAX5 expression in B lymphoma cells respectively activated and repressed hTERT transcription. Thus, PAX5 is strongly implicated in hTERT expression activation in telomerase-positive B cells. These results reveal differences between the hTERT methylation patterns in telomerase-positive carcinoma cells and telomerase-positive normal B cells. The potential of hTERT methylation as a cancer biomarker was evaluated and applied to the detection of metastasis. We have shown that hTERT methylation correlates with the cytological diagnosis in cerebrospinal fluids. Our results suggest a model of hTERT gene regulation, which helps us to better understand how hTERT transcription is regulated by CTCF in methylation-dependant and independent mechanisms. Our data also indicate that hTERT methylation is a promising new cancer biomarker.

Post-genomic update on a classical candidate gene for coronary artery disease: ESR1.

BACKGROUND: After age, sex is the most important risk factor for coronary artery disease (CAD). The mechanism through which women are protected from CAD is still largely unknown, but the observed sex difference suggests the involvement of the reproductive steroid hormone signaling system. Genetic association studies of the gene-encoding Estrogen Receptor α (ESR1) have shown conflicting results, although only a limited range of variation in the gene has been investigated. METHODS AND RESULTS: We exploited information made available by advanced new methods and resources in complex disease genetics to revisit the question of ESR1's role in risk of CAD. We performed a meta-analysis of 14 genome-wide association studies (CARDIoGRAM discovery analysis, N=≈87,000) to search for population-wide and sex-specific associations between CAD risk and common genetic variants throughout the coding, noncoding, and flanking regions of ESR1. In addition to samples from the MIGen (N=≈6000), WTCCC (N=≈7400), and Framingham (N=≈3700) studies, we extended this search to a larger number of common and uncommon variants by imputation into a panel of haplotypes constructed using data from the 1000 Genomes Project. Despite the widespread expression of ERα in vascular tissues, we found no evidence for involvement of common or low-frequency genetic variation throughout the ESR1 gene in modifying risk of CAD, either in the general population or as a function of sex. CONCLUSIONS: We suggest that future research on the genetic basis of sex-related differences in CAD risk should initially prioritize other genes in the reproductive steroid hormone biosynthesis system.

Characterization of the hcnABC gene cluster encoding hydrogen cyanide synthase and anaerobic regulation by ANR in the strictly aerobic biocontrol agent Pseudomonas fluorescens CHA0.

The secondary metabolite hydrogen cyanide (HCN) is produced by Pseudomonas fluorescens from glycine, essentially under microaerophilic conditions. The genetic basis of HCN synthesis in P. fluorescens CHA0 was investigated. The contiguous structural genes hcnABC encoding HCN synthase were expressed from the T7 promoter in Escherichia coli, resulting in HCN production in this bacterium. Analysis of the nucleotide sequence of the hcnABC genes showed that each HCN synthase subunit was similar to known enzymes involved in hydrogen transfer, i.e., to formate dehydrogenase (for HcnA) or amino acid oxidases (for HcnB and HcnC). These similarities and the presence of flavin adenine dinucleotide- or NAD(P)-binding motifs in HcnB and HcnC suggest that HCN synthase may act as a dehydrogenase in the reaction leading from glycine to HCN and CO2. The hcnA promoter was mapped by primer extension; the -40 sequence (TTGGC ... ATCAA) resembled the consensus FNR (fumarate and nitrate reductase regulator) binding sequence (TTGAT ... ATCAA). The gene encoding the FNR-like protein ANR (anaerobic regulator) was cloned from P. fluorescens CHA0 and sequenced. ANR of strain CHA0 was most similar to ANR of P. aeruginosa and CydR of Azotobacter vinelandii. An anr mutant of P. fluorescens (CHA21) produced little HCN and was unable to express an hcnA-lacZ translational fusion, whereas in wild-type strain CHA0, microaerophilic conditions strongly favored the expression of the hcnA-lacZ fusion. Mutant CHA21 as well as an hcn deletion mutant were impaired in their capacity to suppress black root rot of tobacco, a disease caused by Thielaviopsis basicola, under gnotobiotic conditions. This effect was most pronounced in water-saturated artificial soil, where the anr mutant had lost about 30% of disease suppression ability, compared with wild-type strain CHA0. These results show that the anaerobic regulator ANR is required for cyanide synthesis in the strictly aerobic strain CHA0 and suggest that ANR-mediated cyanogenesis contributes to the suppression of black root rot.

The fou2 gain-of-function allele and the wild-type allele of Two Pore Channel 1 contribute to different extents or by different mechanisms to defense gene expression in Arabidopsis.

The fatty acid oxygenation up-regulated 2 (fou2) mutant in Arabidopsis thaliana creates a gain-of-function allele in a non-selective cation channel encoded by the Two Pore Channel 1 (TPC1) gene. This mutant genetically implicates cation fluxes in the control of the positive feedback loop whereby jasmonic acid (JA) stimulates its own synthesis. In this study we observed extensive transcriptome reprogramming in healthy fou2 leaves closely resembling that induced by treatment with methyl jasmonate, biotic stresses and the potassium starvation response. Proteomic analysis of fou2 leaves identified increased levels of seven biotic stress- and JA-inducible proteins. In agreement with these analyses, epistasis studies performed by crossing fou2 with aos indicated that elevated levels of JA in fou2 are the major determinant of the mutant phenotype. In addition, generation of fou2 aba1-5, fou2 etr1-1 and fou2 npr1-1 double mutants showed that the fou2 phenotype was only weakly affected by ABA levels and unaffected by mutations in NPR1 and ETR1. The results now suggest possible mechanisms whereby fou2 could induce JA synthesis/signaling early in the wound response. In contrast to fou2, transcriptome analysis of a loss-of-function allele of TPC1, tpc1-2, revealed no differential expression of JA biosynthesis genes in resting leaves. However, the analysis disclosed reduced mRNA levels of the pathogenesis-related genes PDF1.2a and THI2.1 in healthy and diseased tpc1-2 leaves. The results suggest that wild-type TPC1 contributes to their expression by mechanisms somewhat different from those affecting their expression in fou2.