High-resolution Transcript Profiling Of The Atypical Biotrophic Interaction Between Theobroma Cacao And The Fungal Pathogen Moniliophthora Perniciosa.

Witches' broom disease (WBD), caused by the hemibiotrophic fungus Moniliophthora perniciosa, is one of the most devastating diseases of Theobroma cacao, the chocolate tree. In contrast to other hemibiotrophic interactions, the WBD biotrophic stage lasts for months and is responsible for the most distinctive symptoms of the disease, which comprise drastic morphological changes in the infected shoots. Here, we used the dual RNA-seq approach to simultaneously assess the transcriptomes of cacao and M. perniciosa during their peculiar biotrophic interaction. Infection with M. perniciosa triggers massive metabolic reprogramming in the diseased tissues. Although apparently vigorous, the infected shoots are energetically expensive structures characterized by the induction of ineffective defense responses and by a clear carbon deprivation signature. Remarkably, the infection culminates in the establishment of a senescence process in the host, which signals the end of the WBD biotrophic stage. We analyzed the pathogen's transcriptome in unprecedented detail and thereby characterized the fungal nutritional and infection strategies during WBD and identified putative virulence effectors. Interestingly, M. perniciosa biotrophic mycelia develop as long-term parasites that orchestrate changes in plant metabolism to increase the availability of soluble nutrients before plant death. Collectively, our results provide unique insight into an intriguing tropical disease and advance our understanding of the development of (hemi)biotrophic plant-pathogen interactions.

SAGE transcript profiling of the juvenile cambial region of Eucalyptus grandis

Despite the importance of Eucalyptus spp. in the pulp and paper industry, functional genomic approaches have only recently been applied to understand wood formation in this genus. We attempted to establish a global view of gene expression in the juvenile cambial region of Eucalyptus grandis Hill ex Maiden. The expression profile was obtained from serial analysis of gene expression (SAGE) library data produced from 3- and 6-year-old trees. Fourteen-base expressed sequence tags (ESTs) were searched against public Eucalyptus ESTs and annotated with GenBank. Altogether 43,304 tags were generated producing 3066 unigenes with three or more copies each, 445 with a putative identity, 215 with unknown function and 2406 without an EST match. The expression profile of the juvenile cambial region revealed the presence of highly frequent transcripts related to general metabolism and energy metabolism, cellular processes, transport, structural components and information pathways. We made a quantitative analysis of a large number of genes involved in the biosynthesis of cellulose, pectin, hemicellulose and lignin. Our findings provide insight into the expression of functionally related genes involved in juvenile wood formation in young fast-growing E. grandis trees.

Transcript profiling of papaya fruit reveals differentially expressed genes associated with fruit ripening

Papaya (Carica papaya L) fruit has a short shelf life due to fast ripening induced by ethylene, but little is known about the genetic control of ripening and attributes of fruit quality. Therefore, we identified ripening-related genes affected by ethylene using cDNA-AFLP (Amplified Fragment Length Polymorphism of cDNA). Transcript profiling of non-induced and ethylene-induced fruit samples was performed, and 71 differentially expressed genes were identified. Among those genes some involved in ethylene biosynthesis, regulation of transcription, and stress responses or plant defence were found (heat shock proteins, polygalacturonase-inhibiting protein, and acyl-CoA oxidases). Several transcription factors were isolated, and except for a 14-3-3 protein, an AP2 domain-containing factor, a salt-tolerant zinc finger protein, and a suppressor of PhyA-105 1, most of them were negatively affected by ethylene, including fragments of transcripts similar to VRN1, and ethylene responsive factors (ERF). With respect to fruit quality, genes related to cell wall structure or metabolism, volatiles or pigment precursors, and vitamin biosynthesis were also found. (C) 2010 Elsevier Ireland Ltd. All rights reserved.

Versatile gene-specific sequence tags for Arabidopsis functional genomics: transcript profiling and reverse genetics applications.

Microarray transcript profiling and RNA interference are two new technologies crucial for large-scale gene function studies in multicellular eukaryotes. Both rely on sequence-specific hybridization between complementary nucleic acid strands, inciting us to create a collection of gene-specific sequence tags (GSTs) representing at least 21,500 Arabidopsis genes and which are compatible with both approaches. The GSTs were carefully selected to ensure that each of them shared no significant similarity with any other region in the Arabidopsis genome. They were synthesized by PCR amplification from genomic DNA. Spotted microarrays fabricated from the GSTs show good dynamic range, specificity, and sensitivity in transcript profiling experiments. The GSTs have also been transferred to bacterial plasmid vectors via recombinational cloning protocols. These cloned GSTs constitute the ideal starting point for a variety of functional approaches, including reverse genetics. We have subcloned GSTs on a large scale into vectors designed for gene silencing in plant cells. We show that in planta expression of GST hairpin RNA results in the expected phenotypes in silenced Arabidopsis lines. These versatile GST resources provide novel and powerful tools for functional genomics.

CATMA, a comprehensive genome-scale resource for silencing and transcript profiling of Arabidopsis genes.

BACKGROUND: The Complete Arabidopsis Transcript MicroArray (CATMA) initiative combines the efforts of laboratories in eight European countries 1 to deliver gene-specific sequence tags (GSTs) for the Arabidopsis research community. The CATMA initiative offers the power and flexibility to regularly update the GST collection according to evolving knowledge about the gene repertoire. These GST amplicons can easily be reamplified and shared, subsets can be picked at will to print dedicated arrays, and the GSTs can be cloned and used for other functional studies. This ongoing initiative has already produced approximately 24,000 GSTs that have been made publicly available for spotted microarray printing and RNA interference. RESULTS: GSTs from the CATMA version 2 repertoire (CATMAv2, created in 2002) were mapped onto the gene models from two independent Arabidopsis nuclear genome annotation efforts, TIGR5 and PSB-EuGène, to consolidate a list of genes that were targeted by previously designed CATMA tags. A total of 9,027 gene models were not tagged by any amplified CATMAv2 GST, and 2,533 amplified GSTs were no longer predicted to tag an updated gene model. To validate the efficacy of GST mapping criteria and design rules, the predicted and experimentally observed hybridization characteristics associated to GST features were correlated in transcript profiling datasets obtained with the CATMAv2 microarray, confirming the reliability of this platform. To complete the CATMA repertoire, all 9,027 gene models for which no GST had yet been designed were processed with an adjusted version of the Specific Primer and Amplicon Design Software (SPADS). A total of 5,756 novel GSTs were designed and amplified by PCR from genomic DNA. Together with the pre-existing GST collection, this new addition constitutes the CATMAv3 repertoire. It comprises 30,343 unique amplified sequences that tag 24,202 and 23,009 protein-encoding nuclear gene models in the TAIR6 and EuGène genome annotations, respectively. To cover the remaining untagged genes, we identified 543 additional GSTs using less stringent design criteria and designed 990 sequence tags matching multiple members of gene families (Gene Family Tags or GFTs) to cover any remaining untagged genes. These latter 1,533 features constitute the CATMAv4 addition. CONCLUSION: To update the CATMA GST repertoire, we designed 7,289 additional sequence tags, bringing the total number of tagged TAIR6-annotated Arabidopsis nuclear protein-coding genes to 26,173. This resource is used both for the production of spotted microarrays and the large-scale cloning of hairpin RNA silencing vectors. All information about the resulting updated CATMA repertoire is available through the CATMA database http://www.catma.org.

Transcript profiling suggests that differential metabolic adaptation of mice to a high fat diet is associated with changes in liver to muscle lipid fluxes.

Genetically homogenous C57Bl/6 mice display differential metabolic adaptation when fed a high fat diet for 9 months. Most become obese and diabetic, but a significant fraction remains lean and diabetic or lean and non-diabetic. Here, we performed microarray analysis of "metabolic" transcripts expressed in liver and hindlimb muscles to evaluate: (i) whether expressed transcript patterns could indicate changes in metabolic pathways associated with the different phenotypes, (ii) how these changes differed from the early metabolic adaptation to short term high fat feeding, and (iii) whether gene classifiers could be established that were characteristic of each metabolic phenotype. Our data indicate that obesity/diabetes was associated with preserved hepatic lipogenic gene expression and increased plasma levels of very low density lipoprotein and, in muscle, with an increase in lipoprotein lipase gene expression. This suggests increased muscle fatty acid uptake, which may favor insulin resistance. In contrast, the lean mice showed a strong reduction in the expression of hepatic lipogenic genes, in particular of Scd-1, a gene linked to sensitivity to diet-induced obesity; the lean and non-diabetic mice presented an additional increased expression of eNos in liver. After 1 week of high fat feeding the liver gene expression pattern was distinct from that seen at 9 months in any of the three mouse groups, thus indicating progressive establishment of the different phenotypes. Strikingly, development of the obese phenotype involved re-expression of Scd-1 and other lipogenic genes. Finally, gene classifiers could be established that were characteristic of each metabolic phenotype. Together, these data suggest that epigenetic mechanisms influence gene expression patterns and metabolic fates.

Spatio-temporal transcript profiling of rice roots and shoots in response to phosphate starvation and recovery.

Using rice (Oryza sativa) as a model crop species, we performed an in-depth temporal transcriptome analysis, covering the early and late stages of Pi deprivation as well as Pi recovery in roots and shoots, using next-generation sequencing. Analyses of 126 paired-end RNA sequencing libraries, spanning nine time points, provided a comprehensive overview of the dynamic responses of rice to Pi stress. Differentially expressed genes were grouped into eight sets based on their responses to Pi starvation and recovery, enabling the complex signaling pathways involved in Pi homeostasis to be untangled. A reference annotation-based transcript assembly was also generated, identifying 438 unannotated loci that were differentially expressed under Pi starvation. Several genes also showed induction of unannotated splice isoforms under Pi starvation. Among these, PHOSPHATE2 (PHO2), a key regulator of Pi homeostasis, displayed a Pi starvation-induced isoform, which was associated with increased translation activity. In addition, microRNA (miRNA) expression profiles after long-term Pi starvation in roots and shoots were assessed, identifying 20 miRNA families that were not previously associated with Pi starvation, such as miR6250. In this article, we present a comprehensive spatio-temporal transcriptome analysis of plant responses to Pi stress, revealing a large number of potential key regulators of Pi homeostasis in plants.

What can digital transcript profiling reveal about human cancers?

Important biological and clinical features of malignancy are reflected in its transcript pattern. Recent advances in gene expression technology and informatics have provided a powerful new means to obtain and interpret these expression patterns. A comprehensive approach to expression profiling is serial analysis of gene expression (SAGE), which provides digital information on transcript levels. SAGE works by counting transcripts and storing these digital values electronically, providing absolute gene expression levels that make historical comparisons possible. SAGE produces a comprehensive profile of gene expression and can be used to search for candidate tumor markers or antigens in a limited number of samples. The Cancer Genome Anatomy Project has created a SAGE database of human gene expression levels for many different tumors and normal reference tissues and provides online tools for viewing, comparing, and downloading expression profiles. Digital expression profiling using SAGE and informatics have been useful for identifying genes that have a role in tumor invasion and other aspects of tumor progression.

Ensemble transcript interaction networks: A case study on Alzheimer's disease

Systems biology techniques are a topic of recent interest within the neurological field. Computational intelligence (CI) addresses this holistic perspective by means of consensus or ensemble techniques ultimately capable of uncovering new and relevant findings. In this paper, we propose the application of a CI approach based on ensemble Bayesian network classifiers and multivariate feature subset selection to induce probabilistic dependences that could match or unveil biological relationships. The research focuses on the analysis of high-throughput Alzheimer's disease (AD) transcript profiling. The analysis is conducted from two perspectives. First, we compare the expression profiles of hippocampus subregion entorhinal cortex (EC) samples of AD patients and controls. Second, we use the ensemble approach to study four types of samples: EC and dentate gyrus (DG) samples from both patients and controls. Results disclose transcript interaction networks with remarkable structures and genes not directly related to AD by previous studies. The ensemble is able to identify a variety of transcripts that play key roles in other neurological pathologies. Classical statistical assessment by means of non-parametric tests confirms the relevance of the majority of the transcripts. The ensemble approach pinpoints key metabolic mechanisms that could lead to new findings in the pathogenesis and development of AD

Construction and validation of a Bovine Innate Immune Microarray

Background: Microarray transcript profiling has the potential to illuminate the molecular processes that are involved in the responses of cattle to disease challenges. This knowledge may allow the development of strategies that exploit these genes to enhance resistance to disease in an individual or animal population. Results: The Bovine Innate Immune Microarray developed in this study consists of 1480 characterised genes identified by literature searches, 31 positive and negative control elements and 5376 cDNAs derived from subtracted and normalised libraries. The cDNA libraries were produced from 'challenged' bovine epithelial and leukocyte cells. The microarray was found to have a limit of detection of 1 pg/mu g of total RNA and a mean slide-to-slide correlation co-efficient of 0.88. The profiles of differentially expressed genes from Concanavalin A ( ConA) stimulated bovine peripheral blood lymphocytes were determined. Three distinct profiles highlighted 19 genes that were rapidly up-regulated within 30 minutes and returned to basal levels by 24 h; 76 genes that were upregulated between 2 - 8 hours and sustained high levels of expression until 24 h and 10 genes that were down-regulated. Quantitative real-time RT-PCR on selected genes was used to confirm the results from the microarray analysis. The results indicate that there is a dynamic process involving gene activation and regulatory mechanisms re-establishing homeostasis in the ConA activated lymphocytes. The Bovine Innate Immune Microarray was also used to determine the cross-species hybridisation capabilities of an ovine PBL sample. Conclusion: The Bovine Innate Immune Microarray has been developed which contains a set of well-characterised genes and anonymous cDNAs from a number of different bovine cell types. The microarray can be used to determine the gene expression profiles underlying innate immune responses in cattle and sheep.

Differential display and suppression subtractive hybridization analysis of the pulp of ripening banana

The fruit of banana undergoes several important physico-chemical changes during ripening. Analysis of gene expression would permit identification of important genes and regulatory elements involved in this process. Therefore, transcript profiling of preclimacteric and climacteric fruit was performed using differential display and Suppression subtractive hybridization. Our analyses resulted in the isolation of 12 differentially expressed cDNAs, which were confirmed by dot-blots and northern blots. Among the sequences identified were sequences homologous to plant aquaporins, adenine nucleotide translocator, immunophilin, legumin-like proteins, deoxyguanosine kinase and omega-3 fatty acid desaturase. Some of these cDNAs correspond to newly isolated genes involved in changes related to the respiratory climacteric, or stress-defense responses. Functional characterization of ripening-associated genes could provide information useful in controlling biochemical pathways that would have an impact on banana quality and shelf life. (C) 2009 Elsevier B.V. All rights reserved.

The transcriptome of the arbuscular mycorrhizal fungus Glomus intraradices (DAOM 197198) reveals functional tradeoffs in an obligate symbiont.

? The arbuscular mycorrhizal symbiosis is arguably the most ecologically important eukaryotic symbiosis, yet it is poorly understood at the molecular level. To provide novel insights into the molecular basis of symbiosis-associated traits, we report the first genome-wide analysis of the transcriptome from Glomus intraradices DAOM 197198. ? We generated a set of 25,906 nonredundant virtual transcripts (NRVTs) transcribed in germinated spores, extraradical mycelium and symbiotic roots using Sanger and 454 sequencing. NRVTs were used to construct an oligoarray for investigating gene expression. ? We identified transcripts coding for the meiotic recombination machinery, as well as meiosis-specific proteins, suggesting that the lack of a known sexual cycle in G. intraradices is not a result of major deletions of genes essential for sexual reproduction and meiosis. Induced expression of genes encoding membrane transporters and small secreted proteins in intraradical mycelium, together with the lack of expression of hydrolytic enzymes acting on plant cell wall polysaccharides, are all features of G. intraradices that are shared with ectomycorrhizal symbionts and obligate biotrophic pathogens. ? Our results illuminate the genetic basis of symbiosis-related traits of the most ancient lineage of plant biotrophs, advancing future research on these agriculturally and ecologically important symbionts.

Transcriptional regulation of MDR1, a gene involved in antifungal drug resistance in Candida albicans

ABSTRACT Upregulation of the Major Facilitator transporter gene MDR1 (Multi_drug Resistance 1) is one of the mechanisms observed in Candida albicans clinical isolates developing resistance to azole antifungal agents. To better understand this phenomenon, the cis-acting regulatory elements present in a modulatable reporter system under the control of the MDR1 promoter were characterized. In an azole-susceptible strain, transcription of this reporter is transiently upregulated in response to either benomyl or H2O2, whereas its expression is constitutively high in an azole-resistant strain (FR2). Two cis-acting regulatory elements, that are necessary and sufficient to convey the same transcriptional responses to a heterologous promoter (CDR2), were identified within the MDR1promoter. The first element, called BRE (for Benomyl Response Element, -296 to -260 with respect to the ATG start codon), is required for benomyl-dependent MDR1 upregulation and for constitutive high expression of MDR1 in FR2. The second element, termed HRE (for H2O2 Response Element, -561 to -520), is required for H2O2-dependent MDR1 upregulation, but is dispensable for constitutive high expression. Two potential binding sites (TTAG/CTAA) for the blip transcription factor Cap1p lie within the HRE. Moreover, inactivation of CAP1 abolished the transient response to H2O2 and diminished significantly the transient response to benomyl. Cap1p, which has been previously implicated in cellular responses to oxidative stress, may thus play a transacting and positive regulatory role in benomyl- and H2O2-dependent transcription of MDR1. However, it is not the only transcription factor involved in the response of MDR1 to benomyl. A minimal BRE element (-290 to -273) that is sufficient to detect in vitro sequence-specific binding of protein complexes in crude extracts prepared from C. albicans was also delimited. Genome-wide transcript profiling analyses undertaken with a matched pair of clinical isolates, one of which being azole-resistant and upregulating MDR1, and with an azole-susceptible strain exposed to benomyl, revealed that genes specifically upregulated by benomyl harbour in their promoters Cap1p binding site(s). This strengthened the idea that Cap1p plays a role in benomyl-dependent upregulation of MDR1. BRE-like sequences were also identified in several genes co-regulated with MDR1 in both conditions, which was consistent with the involvement of the BRE in both processes. A set of 147 mutants lacking a single transcription factor gene was next screened for loss of MDR1response to benomyl. Unfortunately, none of the tested mutants showed a loss of benomyl-dependent MDR1 upregulation. Nevertheless, a significant diminution of the response was observed in the mutants in which the MADS-box transcription factor Mcm1p and the C2H2 zinc finger transcription factor orf19.13374p were inactivated, suggesting that Mcm1p and orf19.13374p are involved in MDR1response to benomyl. Interestingly, the BRE contains a perfect match to the binding consensus of Mcm1p, raising the possibility that MDR1may be a direct target of this transcriptional activator. In conclusion, while the identity of the trans-acting factors that bind to the BRE and HRE remains to be confirmed, the tools we have developed during characterization of the cis-acting elements of the MDR1promoter should now serve to elucidate the nature of the components that modulate its activity. RESUME La surexpression du gène MDR1 (pour Résistance Multidrogue 1), qui code pour un transporteur de la famille des Major Facilitators, est l'un des mécanismes observés dans les isolats cliniques de la levure Candida albicans développant une résistance aux agents antifongiques appelés azoles. Pour mieux comprendre ce phénomène, les éléments de régulation agissant en cis dans un système rapporteur modulable sous le contrôle du promoteur MDR1 ont été caractérisés. Dans une souche sensible aux azoles, la transcription de ce rapporteur est transitoirement surélevée en réponse soit au bénomyl soit à l'agent oxydant H2O2, alors que son expression est constitutivement élevée dans une souche résistante aux azoles (souche FR2). Deux éléments de régulation agissant en cis, nécessaires et suffisants pour transmettre les mêmes réponses transcriptionnelles à un promoteur hétérologue (CDR2), ont été identifiés dans le promoteur MDR1. Le premier élément, appelé BRE (pour Elément de Réponse au Bénomyl, de -296 à -260 par rapport au codon d'initiation ATG) est requis pour la surexpression de MDR1dépendante du bénomyl et pour l'expression constitutive de MDR1 dans FR2. Le deuxième élément, appelé HRE (pour Elément de Réponse à l'H2O2, de -561 à -520), est requis pour la surexpression de MDR1 dépendante de l'H2O2, mais n'est pas impliqué dans l'expression constitutive du gène MDR1. Deux sites de fixation potentiels (TTAG/CTAA) pour le facteur de transcription Cap1p ont été identifiés dans l'élément HRE. De plus, l'inactivation de CAP1 abolit la réponse transitoire à l'H2O2 et diminua significativement la réponse transitoire au bénomyl. Cap1p, qui est impliqué dans les réponses de la cellule au stress oxydatif, doit donc jouer un rôle positif en trans dans la surexpression de MDR1 dépendante du bénomyl et de l'H2O2. Cependant, ce n'est pas le seul facteur de transcription impliqué dans la réponse au bénomyl. Un élément BRE d'une longueur minimale (de -290 à -273) a également été défini et est suffisant pour détecter une interaction spécifique in vitro avec des protéines provenant d'extraits bruts de C. albicans. L'analyse du profil de transcription d'une paire d'isolats cliniques comprenant une souche résistante aux azoles surexprimant MDR1, et d'une souche sensible aux azoles exposée au bénomyl, a révélé que les gènes spécifiquement surexprimés par le bénomyl contiennent dans leurs promoteurs un ou plusieurs sites de fixation pour Cap1p. Ceci renforce l'idée que Cap1p joue un rôle dans la surexpression de MDR1dépendante du bénomyl. Une ou deux séquences ressemblant à l'élément BRE ont également été identifiées dans la plupart des gènes corégulés avec MDR1 dans ces deux conditions, ce qui était attendu compte-tenu du rôle joué par cet élément dans les deux processus. Une collection de 147 mutants dans lesquels un seul facteur de transcription est inactivé a été testée pour la perte de réponse au bénomyl de MDR1. Malheureusement, la surexpression de MDR1 dépendante du bénomyl n'a été perdue dans aucun des mutants testés. Néanmoins, une diminution significative de la réponse a été observée chez des mutants dans lesquels le facteur de transcription à MADS-box Mcm1p et le facteur de transcription à doigts de zinc de type C2H2 orf19.13374p ont été inactivés, suggérant que Mcm1p et orf19.13374p sont impliqués dans la réponse de MDR1au bénomyl. Il est intéressant de noter que la BRE contient une séquence qui s'aligne parfaitement avec la séquence consensus du site de fixation de Mcm1p, ce qui soulève la possibilité que MDR1 pourrait être une cible directe de cet activateur transcriptionnel. En conclusion, alors que l'identité des facteurs agissant en trans en se fixant à la BRE et à la HRE reste à être confirmée, les outils que nous avons développés au cours de la caractérisation des éléments agissant en cis sur le promoteur MDR1 peut maintenant servir à élucider la nature des composants modulant son activité.

Contribution of CgPDR1-Regulated Genes in Enhanced Virulence of Azole-Resistant Candida glabrata.

In Candida glabrata, the transcription factor CgPdr1 is involved in resistance to azole antifungals via upregulation of ATP binding cassette (ABC)-transporter genes including at least CgCDR1, CgCDR2 and CgSNQ2. A high diversity of GOF (gain-of-function) mutations in CgPDR1 exists for the upregulation of ABC-transporters. These mutations enhance C. glabrata virulence in animal models, thus indicating that CgPDR1 might regulate the expression of yet unidentified virulence factors. We hypothesized that CgPdr1-dependent virulence factor(s) should be commonly regulated by all GOF mutations in CgPDR1. As deduced from transcript profiling with microarrays, a high number of genes (up to 385) were differentially regulated by a selected number (7) of GOF mutations expressed in the same genetic background. Surprisingly, the transcriptional profiles resulting from expression of GOF mutations showed minimal overlap in co-regulated genes. Only two genes, CgCDR1 and PUP1 (for PDR1upregulated and encoding a mitochondrial protein), were commonly upregulated by all tested GOFs. While both genes mediated azole resistance, although to different extents, their deletions in an azole-resistant isolate led to a reduction of virulence and decreased tissue burden as compared to clinical parents. As expected from their role in C. glabrata virulence, the two genes were expressed as well in vitro and in vivo. The individual overexpression of these two genes in a CgPDR1-independent manner could partially restore phenotypes obtained in clinical isolates. These data therefore demonstrate that at least these two CgPDR1-dependent and -upregulated genes contribute to the enhanced virulence of C. glabrata that acquired azole resistance.

Hepatic glucose sensing is required to preserve β cell glucose competence.

Liver glucose metabolism plays a central role in glucose homeostasis and may also regulate feeding and energy expenditure. Here we assessed the impact of glucose transporter 2 (Glut2) gene inactivation in adult mouse liver (LG2KO mice). Loss of Glut2 suppressed hepatic glucose uptake but not glucose output. In the fasted state, expression of carbohydrate-responsive element-binding protein (ChREBP) and its glycolytic and lipogenic target genes was abnormally elevated. Feeding, energy expenditure, and insulin sensitivity were identical in LG2KO and control mice. Glucose tolerance was initially normal after Glut2 inactivation, but LG2KO mice exhibited progressive impairment of glucose-stimulated insulin secretion even though β cell mass and insulin content remained normal. Liver transcript profiling revealed a coordinated downregulation of cholesterol biosynthesis genes in LG2KO mice that was associated with reduced hepatic cholesterol in fasted mice and reduced bile acids (BAs) in feces, with a similar trend in plasma. We showed that chronic BAs or farnesoid X receptor (FXR) agonist treatment of primary islets increases glucose-stimulated insulin secretion, an effect not seen in islets from Fxr-/- mice. Collectively, our data show that glucose sensing by the liver controls β cell glucose competence and suggest BAs as a potential mechanistic link.