Identification of transcripts with enriched expression in the developing and adult pancreas

Genome-wide association study on overweight in Brazilian children with asthma: Old stories and new discoveries

INTRODUCTION

Overweight and obesity are chronic and multifactorial diseases with a strong genetic component contributing to weight gain across all age groups. This study aimed to conduct a Genome-wide Association Study (GWAS) on a cohort of 1,004 Brazilian children (5-11 years old) to identify specific DNA regions associated with susceptibility to overweight.

METHODS

The GWAS was performed on children participating in the SCAALA (Asthma and Allergy Social Changes in Latin America) program, with participants classified as either overweight or non-overweight. Genotyping was carried out using the Illumina 2.5 Human Omni bead chip. Using ELISA, cytokine levels (IL-5, IL-13, IL-10, and IFN) were measured in the blood culture supernatant. Furthermore, pathway analyses were conducted utilizing the Gene Ontology tool.

RESULTS

Our analysis revealed eight significant signals distributed across the genome. The most prominent single nucleotide variant (SNV) was identified in the IL1R1 gene, followed by three variants in the LOC105377841 region (located between the ADH5P4 and EYS genes), as well as variants in the KNTC1, RAPTOR, and DSCAM genes. Among the identified variants, three (IL1R1, RAPTOR, and DSCAM) are associated with immune mechanisms, one (ST18) is linked to the death pathway, and one (KNTC1) is associated with mitotic spindle assembly. The genetic risk score analysis demonstrated that having one or more variants among the six analyzed significantly increased the risk of being overweight by eightfold.

CONCLUSIONS

Our study uncovered genetic loci within pathways with strong biological plausibility, including identifying a novel region (LOC105377841) not previously associated with overweight. Understanding the genetic variants involved in overweight and obesity is crucial for advancing our knowledge of these diseases, particularly within mixed populations such as the Brazilian population.

A single-cell time-lapse of mouse prenatal development from gastrula to birth

The house mouse (Mus musculus) is an exceptional model system, combining genetic tractability with close evolutionary affinity to humans1,2. Mouse gestation lasts only 3 weeks, during which the genome orchestrates the astonishing transformation of a single-cell zygote into a free-living pup composed of more than 500 million cells. Here, to establish a global framework for exploring mammalian development, we applied optimized single-cell combinatorial indexing3 to profile the transcriptional states of 12.4 million nuclei from 83 embryos, precisely staged at 2- to 6-hour intervals spanning late gastrulation (embryonic day 8) to birth (postnatal day 0). From these data, we annotate hundreds of cell types and explore the ontogenesis of the posterior embryo during somitogenesis and of kidney, mesenchyme, retina and early neurons. We leverage the temporal resolution and sampling depth of these whole-embryo snapshots, together with published data4-8 from earlier timepoints, to construct a rooted tree of cell-type relationships that spans the entirety of prenatal development, from zygote to birth. Throughout this tree, we systematically nominate genes encoding transcription factors and other proteins as candidate drivers of the in vivo differentiation of hundreds of cell types. Remarkably, the most marked temporal shifts in cell states are observed within one hour of birth and presumably underlie the massive physiological adaptations that must accompany the successful transition of a mammalian fetus to life outside the womb.

Single-cell transcriptomic and spatial landscapes of the developing human pancreas

Current differentiation protocols have not been successful in reproducibly generating fully functional human beta cells in vitro, partly due to incomplete understanding of human pancreas development. Here, we present detailed transcriptomic analysis of the various cell types of the developing human pancreas, including their spatial gene patterns. We integrated single-cell RNA sequencing with spatial transcriptomics at multiple developmental time points and revealed distinct temporal-spatial gene cascades. Cell trajectory inference identified endocrine progenitor populations and branch-specific genes as the progenitors differentiate toward alpha or beta cells. Spatial differentiation trajectories indicated that Schwann cells are spatially co-located with endocrine progenitors, and cell-cell connectivity analysis predicted that they may interact via L1CAM-EPHB2 signaling. Our integrated approach enabled us to identify heterogeneity and multiple lineage dynamics within the mesenchyme, showing that it contributed to the exocrine acinar cell state. Finally, we have generated an interactive web resource for investigating human pancreas development for the research community.

Colorectal Cancer-Associated Genes Are Associated with Tooth Agenesis and May Have a Role in Tooth Development

Previously reported co-occurrence of colorectal cancer (CRC) and tooth agenesis (TA) and the overlap in disease-associated gene variants suggest involvement of similar molecular pathways. Here, we took an unbiased approach and tested genome-wide significant CRC-associated variants for association with isolated TA. Thirty single nucleotide variants (SNVs) in CRC-predisposing genes/loci were genotyped in a discovery dataset composed of 440 individuals with and without isolated TA. Genome-wide significant associations were found between TA and ATF1 rs11169552 (P = 4.36 × 10-10) and DUSP10 rs6687758 (P = 1.25 × 10-9), and positive association found with CASC8 rs10505477 (P = 8.2 × 10-5). Additional CRC marker haplotypes were also significantly associated with TA. Genotyping an independent dataset consisting of 52 cases with TA and 427 controls confirmed the association with CASC8. Atf1 and Dusp10 expression was detected in the mouse developing teeth from early bud stages to the formation of the complete tooth, suggesting a potential role for these genes and their encoded proteins in tooth development. While their individual contributions in tooth development remain to be elucidated, these genes may be considered candidates to be tested in additional populations.

Postimplantation Mga expression and embryonic lethality of two gene-trap alleles

BACKGROUND

The dual-specificity T-box/basic helix-loop-helix leucine zipper transcription factor MGA is part of the MAX-interacting network of proteins. In the mouse, MGA is necessary for the survival of the pluripotent epiblast cells of the peri-implantation embryo and a null, gene-trap allele MgaGt results in embryonic lethality shortly after implantation. We have used this allele to document expression of Mga in postimplantation embryos and also investigated a second, hypomorphic gene-trap allele, MgaInv.

RESULTS

Compound heterozygotes, MgaGt/MgaInv, die prior to midgestation. The extraembryonic portion of the embryos appears to develop relatively normally while the embryonic portion, including the pluripotent cells of the epiblast, is severely retarded by E7.5. Mga expression is initially limited to the pluripotent inner cell mass of the blastocyst and epiblast, but during organogenesis it is widely expressed, notably in the central nervous system and sensory organs, reproductive and excretory systems, heart, somites and limbs.

CONCLUSIONS

Widespread yet specific areas of expression of Mga during organogenesis raise the possibility that the transcription factor may play roles in controlling proliferation and potency in the progenitor cell populations of different organ systems. Documentation of these patterns sets the stage for the investigation of specific progenitor cell types.

Modeling coexistence of oscillation and Delta/Notch-mediated lateral inhibition in pancreas development and neurogenesis

During pancreas development, Neurog3 positive endocrine progenitors are specified by Delta/Notch (D/N) mediated lateral inhibition in the growing ducts. During neurogenesis, genes that determine the transition from the proneural state to neuronal or glial lineages are oscillating before their expression is sustained. Although the basic gene regulatory network is very similar, cycling gene expression in pancreatic development was not investigated yet, and previous simulations of lateral inhibition in pancreas development excluded by design the possibility of oscillations. To explore this possibility, we developed a dynamic model of a growing duct that results in an oscillatory phase before the determination of endocrine progenitors by lateral inhibition. The basic network (D/N + Hes1 + Neurog3) shows scattered, stable Neurog3 expression after displaying transient expression. Furthermore, we included the Hes1 negative feedback as previously discussed in neurogenesis and show the consequences for Neurog3 expression in pancreatic duct development. Interestingly, a weakened HES1 action on the Hes1 promoter allows the coexistence of stable patterning and oscillations. In conclusion, cycling gene expression and lateral inhibition are not mutually exclusive. In this way, we argue for a unified mode of D/N mediated lateral inhibition in neurogenic and pancreatic progenitor specification.

Single cell transcriptomic profiling of mouse pancreatic progenitors

The heterogeneity of the developing pancreatic epithelium and low abundance of endocrine progenitors limit the information derived from traditional expression studies. To identify genes that characterize early developmental tissues composed of multiple progenitor lineages, we applied single-cell RNA-Seq to embryonic day (e)13.5 mouse pancreata and performed integrative analysis with single cell data from mature pancreas. We identified subpopulations expressing macrophage or endothelial markers and new pancreatic progenitor markers. We also identified potential α-cell precursors expressing glucagon (Gcg) among the e13.5 pancreatic cells. Despite their high Gcg expression levels, these cells shared greater transcriptomic similarity with other e13.5 cells than with adult α-cells, indicating their immaturity. Comparative analysis identified the sodium-dependent neutral amino acid transporter, Slc38a5, as a characteristic gene expressed in α-cell precursors but not mature cells. By immunofluorescence analysis, we observed SLC38A5 expression in pancreatic progenitors, including in a subset of NEUROG3+ endocrine progenitors and MAFB+ cells and in all GCG+ cells. Expression declined in α-cells during late gestation and was absent in the adult islet. Our results suggest SLC38A5 as an early marker of α-cell lineage commitment.

Glucose tolerance female-specific QTL mapped in collaborative cross mice

Type-2 diabetes (T2D) is a complex metabolic disease characterized by impaired glucose tolerance. Despite environmental high risk factors, host genetic background is a strong component of T2D development. Herein, novel highly genetically diverse strains of collaborative cross (CC) lines from mice were assessed to map quantitative trait loci (QTL) associated with variations of glucose-tolerance response. In total, 501 mice of 58 CC lines were maintained on high-fat (42 % fat) diet for 12 weeks. Thereafter, an intraperitoneal glucose tolerance test (IPGTT) was performed for 180 min. Subsequently, the values of Area under curve for the glucose at zero and 180 min (AUC_0-180), were measured, and used for QTL mapping. Heritability and coefficient of variations in glucose tolerance (CVg) were calculated. One-way analysis of variation was significant (P < 0.001) for AUC_0-180 between the CC lines as well between both sexes. Despite Significant variations for both sexes, QTL analysis was significant, only for females, reporting a significant female-sex-dependent QTL (~2.5 Mbp) associated with IPGTT AUC_0-180 trait, located on Chromosome 8 (32-34.5 Mbp, containing 51 genes). Gene browse revealed QTL for body weight/size, genes involved in immune system, and two main protein-coding genes involved in the Glucose homeostasis, Mboat4 and Leprotl1. Heritability and coefficient of genetic variance (CVg) were 0.49 and 0.31 for females, while for males, these values 0.34 and 0.22, respectively. Our findings demonstrate the roles of genetic factors controlling glucose tolerance, which significantly differ between sexes requiring independent studies for females and males toward T2D prevention and therapy.

Myt3 suppression sensitizes islet cells to high glucose-induced cell death via Bim induction

Diabetes is a chronic disease that results from the body's inability to properly control circulating blood glucose levels. The loss of glucose homoeostasis can arise from a loss of β-cell mass because of immune-cell-mediated attack, as in type 1 diabetes, and/or from dysfunction of individual β-cells (in conjunction with target organ insulin resistance), as in type 2 diabetes. A better understanding of the transcriptional pathways regulating islet-cell survival is of great importance for the development of therapeutic strategies that target β-cells for diabetes. To this end, we previously identified the transcription factor Myt3 as a pro-survival factor in islets following acute suppression of Myt3 in vitro. To determine the effects of Myt3 suppression on islet-cell survival in vivo, we used an adenovirus to express an shRNA targeting Myt3 in syngeneic optimal and marginal mass islet transplants, and demonstrate that suppression of Myt3 impairs the function of marginal mass grafts. Analysis of grafts 5 weeks post-transplant revealed that grafts transduced with the shMyt3 adenovirus contained ~20% the number of transduced cells as grafts transduced with a control adenovirus. In fact, increased apoptosis and significant cell loss in the shMyt3-transduced grafts was evident after only 5 days, suggesting that Myt3 suppression sensitizes islet cells to stresses present in the early post-transplant period. Specifically, we find that Myt3 suppression sensitizes islet cells to high glucose-induced cell death via upregulation of the pro-apoptotic Bcl2 family member Bim. Taken together these data suggest that Myt3 may be an important link between glucotoxic and immune signalling pathways.

β-Cell-Specific Mafk Overexpression Impairs Pancreatic Endocrine Cell Development

The MAF family transcription factors are homologs of v-Maf, the oncogenic component of the avian retrovirus AS42. They are subdivided into 2 groups, small and large MAF proteins, according to their structure, function, and molecular size. MAFK is a member of the small MAF family and acts as a dominant negative form of large MAFs. In previous research we generated transgenic mice that overexpress MAFK in order to suppress the function of large MAF proteins in pancreatic β-cells. These mice developed hyperglycemia in adulthood due to impairment of glucose-stimulated insulin secretion. The aim of the current study is to examine the effects of β-cell-specific Mafk overexpression in endocrine cell development. The developing islets of Mafk-transgenic embryos appeared to be disorganized with an inversion of total numbers of insulin+ and glucagon+ cells due to reduced β-cell proliferation. Gene expression analysis by quantitative RT-PCR revealed decreased levels of β-cell-related genes whose expressions are known to be controlled by large MAF proteins. Additionally, these changes were accompanied with a significant increase in key β-cell transcription factors likely due to compensatory mechanisms that might have been activated in response to the β-cell loss. Finally, microarray comparison of gene expression profiles between wild-type and transgenic pancreata revealed alteration of some uncharacterized genes including Pcbd1, Fam132a, Cryba2, and Npy, which might play important roles during pancreatic endocrine development. Taken together, these results suggest that Mafk overexpression impairs endocrine development through a regulation of numerous β-cell-related genes. The microarray analysis provided a unique data set of differentially expressed genes that might contribute to a better understanding of the molecular basis that governs the development and function of endocrine pancreas.

PDX1 binds and represses hepatic genes to ensure robust pancreatic commitment in differentiating human embryonic stem cells

Inactivation of the Pancreatic and Duodenal Homeobox 1 (PDX1) gene causes pancreatic agenesis, which places PDX1 high atop the regulatory network controlling development of this indispensable organ. However, little is known about the identity of PDX1 transcriptional targets. We simulated pancreatic development by differentiating human embryonic stem cells (hESCs) into early pancreatic progenitors and subjected this cell population to PDX1 chromatin immunoprecipitation sequencing (ChIP-seq). We identified more than 350 genes bound by PDX1, whose expression was upregulated on day 17 of differentiation. This group included known PDX1 targets and many genes not previously linked to pancreatic development. ChIP-seq also revealed PDX1 occupancy at hepatic genes. We hypothesized that simultaneous PDX1-driven activation of pancreatic and repression of hepatic programs underlie early divergence between pancreas and liver. In HepG2 cells and differentiating hESCs, we found that PDX1 binds and suppresses expression of endogenous liver genes. These findings rebrand PDX1 as a context-dependent transcriptional repressor and activator within the same cell type.

•

Early pancreatic progenitor (ePP) cells are efficiently derived from hESCs

•

High levels of the homeobox transcription factor PDX1 label ePP cells

•

PDX1 binds a battery of foregut/midgut and early pancreatic genes in ePP cells

•

PDX1 binds and represses hepatic genes

In vivo genome-wide analysis of multiple tissues identifies gene regulatory networks, novel functions and downstream regulatory genes for Bapx1 and its co-regulation with Sox9 in the mammalian vertebral column

Background

Vertebrate organogenesis is a highly complex process involving sequential cascades of transcription factor activation or repression. Interestingly a single developmental control gene can occasionally be essential for the morphogenesis and differentiation of tissues and organs arising from vastly disparate embryological lineages.

Results

Here we elucidated the role of the mammalian homeobox gene Bapx1 during the embryogenesis of five distinct organs at E12.5 - vertebral column, spleen, gut, forelimb and hindlimb - using expression profiling of sorted wildtype and mutant cells combined with genome wide binding site analysis. Furthermore we analyzed the development of the vertebral column at the molecular level by combining transcriptional profiling and genome wide binding data for Bapx1 with similarly generated data sets for Sox9 to assemble a detailed gene regulatory network revealing genes previously not reported to be controlled by either of these two transcription factors.

Conclusions

The gene regulatory network appears to control cell fate decisions and morphogenesis in the vertebral column along with the prevention of premature chondrocyte differentiation thus providing a detailed molecular view of vertebral column development.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1072) contains supplementary material, which is available to authorized users.

An integrated cell purification and genomics strategy reveals multiple regulators of pancreas development

The regulatory logic underlying global transcriptional programs controlling development of visceral organs like the pancreas remains undiscovered. Here, we profiled gene expression in 12 purified populations of fetal and adult pancreatic epithelial cells representing crucial progenitor cell subsets, and their endocrine or exocrine progeny. Using probabilistic models to decode the general programs organizing gene expression, we identified co-expressed gene sets in cell subsets that revealed patterns and processes governing progenitor cell development, lineage specification, and endocrine cell maturation. Purification of Neurog3 mutant cells and module network analysis linked established regulators such as Neurog3 to unrecognized gene targets and roles in pancreas development. Iterative module network analysis nominated and prioritized transcriptional regulators, including diabetes risk genes. Functional validation of a subset of candidate regulators with corresponding mutant mice revealed that the transcription factors Etv1, Prdm16, Runx1t1 and Bcl11a are essential for pancreas development. Our integrated approach provides a unique framework for identifying regulatory genes and functional gene sets underlying pancreas development and associated diseases such as diabetes mellitus.

Discovery of specific pancreas developmental regulators has accelerated in recent years. In contrast, the global regulatory programs controlling pancreas development are poorly understood compared to other organs or tissues like heart or blood. Decoding this regulatory logic may accelerate development of replacement organs from renewable sources like stem cells, but this goal requires identification of regulators and assessment of their functions on a global scale. To address this important challenge for pancreas biology, we combined purification of normal and mutant cells with genome-scale methods to generate and analyze expression profiles from developing pancreas cells. Our work revealed regulatory gene sets governing development of pancreas progenitor cells and their progeny. Our integrative approach nominated multiple pancreas developmental regulators, including suspected risk genes for human diabetes, which we validated by phenotyping mutant mice on a scale not previously reported. Selection of these candidate regulators was unbiased; thus it is remarkable that all were essential for pancreatic islet development. Thus, our studies provide a new heuristic resource for identifying genetic functions underlying pancreas development and diseases like diabetes mellitus.

New discoveries of old SON: a link between RNA splicing and cancer

The SON protein is a ubiquitously expressed DNA- and RNA-binding protein primarily localized to nuclear speckles. Although several early studies implicated SON in DNA-binding, tumorigenesis and apoptosis, functional significance of this protein had not been recognized until recent studies discovered SON as a novel RNA splicing co-factor. During constitutive RNA splicing, SON ensures efficient intron removal from the transcripts containing suboptimal splice sites. Importantly, SON-mediated splicing is required for proper processing of selective transcripts related to cell cycle, microtubules, centrosome maintenance, and genome stability. Moreover, SON regulates alternative splicing of RNAs from the genes involved in apoptosis and epigenetic modification. In addition to the role in RNA splicing, SON has an ability to suppress transcriptional activation at certain promoter/enhancer DNA sequences. Considering the multiple SON target genes which are directly involved in cell proliferation, genome stability and chromatin modifications, SON is an emerging player in gene regulation during cancer development and progression. Here, we summarize available information from several early studies on SON, and highlight recent discoveries describing molecular mechanisms of SON-mediated gene regulation. We propose that our future effort on better understanding of diverse SON functions would reveal novel targets for cancer therapy.

Profiling of embryonic stem cell differentiation

Embryonic stem (ES) cells have been shown to recapitulate normal developmental stages. They are therefore a highly useful tool in the study of developmental biology. Profiling of ES cell-derived cells has yielded important information about the characteristics of differentiated cells, and allowed the identification of novel marker genes and pathways of differentiation. In this review, we focus on recent results from profiling studies of mouse embryos, human islets, and human ES cell-derived differentiated cells from several research groups. Global gene expression data from mouse embryos have been used to identify novel genes or pathways involved in the developmental process, and to search for transcription factors that regulate direct reprogramming. We introduce gene expression databases of human pancreas cells (Beta Cell Gene Atlas, EuroDia database), and summarize profiling studies of islet- or human ES cell-derived pancreatic cells, with a focus on gene expression, microRNAs, epigenetics, and protein expression. Then, we describe our gene expression profile analyses and our search for novel endoderm, or pancreatic, progenitor marker genes. We differentiated mouse ES cells into mesendoderm, definitive endoderm (DE), mesoderm, ectoderm, and Pdx1-expressing pancreatic lineages, and performed DNA microarray analyses. Genes specifically expressed in DE, and/or in Pdx1-expressing cells, were extracted and their expression patterns in normal embryonic development were studied by in situ hybridization. Out of 54 genes examined, 27 were expressed in the DE of E8.5 mouse embryos, and 15 genes were expressed in distinct domains in the pancreatic buds of E14.5 mouse embryos. Akr1c19, Aebp2, Pbxip1, and Creb3l1 were all novel, and none has been described as being expressed, either in the DE, or in the pancreas. By introducing the profiling results of ES cell-derived cells, the benefits of using ES cells to study early embryonic development will be discussed.

Characterization of the transcriptional activity of the basic helix-loop-helix (bHLH) transcription factor Atoh8

The atonal-related Neurogenin/NeuroD family of basic helix-loop-helix (bHLH) transcription factors comprises potent inducers of neuronal and endocrine differentiation programs in the nervous and digestive system. Atonal homolog 8 (Atoh8) displays high similarity in the bHLH domain with NeuroD proteins. Yet, available evidences indicate that Atoh8 has distinctive features including a ubiquitous expression pattern in embryonic tissues and the ability to inhibit differentiation. To gain insights into Atoh8 function, we aimed at identifying Atoh8 targets and investigated the effects of Atoh8 on global gene expression patterns in pancreatic mPAC cells, a model of bHLH-dependent endocrine differentiation. Our data reveal that Atoh8 is a weak transcriptional activator and does not exhibit proendocrine activity. Conversely, it blocks the induction of a reduced group of gene targets of the atonal-related proendocrine factor Neurogenin3. We show that Atoh8 lacks a transactivation domain and possesses intrinsic repressor activity that depends on a conserved Proline-rich domain. Atoh8 binds the ubiquitous E protein E47 and its ability to repress transcription may partly result from its ability to inhibit E47/E47 and Neurogenin3/E47 dimer activities. These results reveal distinctive transcriptional properties of Atoh8 within the atonal-related bHLH family that may be associated with the acquisition of new biological functions.

Reconstituting pancreas development from purified progenitor cells reveals genes essential for islet differentiation

Developmental biology is challenged to reveal the function of numerous candidate genes implicated by recent genome-scale studies as regulators of organ development and diseases. Recapitulating organogenesis from purified progenitor cells that can be genetically manipulated would provide powerful opportunities to dissect such gene functions. Here we describe systems for reconstructing pancreas development, including islet β-cell and α-cell differentiation, from single fetal progenitor cells. A strict requirement for native genetic regulators of in vivo pancreas development, such as Ngn3, Arx, and Pax4, revealed the authenticity of differentiation programs in vitro. Efficient genetic screens permitted by this system revealed that Prdm16 is required for pancreatic islet development in vivo. Discovering the function of genes regulating pancreas development with our system should enrich strategies for regenerating islets for treating diabetes mellitus.

The pancreatic β cell and type 1 diabetes: innocent bystander or active participant?

Type 1 diabetes mellitus (T1DM) is a chronic disease resulting from destruction of insulin-producing pancreatic β cells. Genetic and environmental factors contribute to T1DM onset. Use of high-throughput DNA sequencing has allowed geneticists to perform genome-wide association studies (GWAS) to identify novel gene loci associated with T1DM. Interestingly, >50% of these genes encode products that are expressed in β cells. These studies, coupled with emerging molecular evidence that β cells are impaired by gain-of-function or loss-of-function of these loci, suggest an active role for the β cell in eliciting its own demise. Although immune dysregulation plays a vital role in T1DM pathogenesis, understanding the mechanisms contributing to β cell failure may lead to new strategies to preserve or improve β cell function in patients with T1DM.

Microarray analysis of novel candidate genes responsible for glucose-stimulated insulin secretion in mouse pancreatic β cell line MIN6

Elucidating the regulation of glucose-stimulated insulin secretion (GSIS) in pancreatic islet β cells is important for understanding and treating diabetes. MIN6 cells, a transformed β-cell line derived from a mouse insulinoma, retain GSIS and are a popular in vitro model for insulin secretion. However, in long-term culture, MIN6 cells' GSIS capacity is lost. We previously isolated a subclone, MIN6 clone 4, from the parental MIN6 cells, that shows well-regulated insulin secretion in response to glucose, glybenclamide, and KCl, even after prolonged culture. To investigate the molecular mechanisms responsible for preserving GSIS in this subclone, we compared four groups of MIN6 cells: Pr-LP (parental MIN6, low passage number), Pr-HP (parental MIN6, high passage number), C4-LP (MIN6 clone 4, low passage number), and C4-HP (MIN6 clone 4, high passage number). Based on their capacity for GSIS, we designated the Pr-LP, C4-LP, and C4-HP cells as “responder cells.” In a DNA microarray analysis, we identified a group of genes with high expression in responder cells (“responder genes”), but extremely low expression in the Pr-HP cells. Another group of genes (“non-responder genes”) was expressed at high levels in the Pr-HP cells, but at extremely low levels in the responder cells. Some of the responder genes were involved in secretory machinery or glucose metabolism, including Chrebp, Scgn, and Syt7. Among the non-responder genes were Car2, Maf, and Gcg, which are not normally expressed in islet β cells. Interestingly, we found a disproportionate number of known imprinted genes among the responder genes. Our findings suggest that the global expression profiling of GSIS-competent and GSIS-incompetent MIN6 cells will help delineate the gene regulatory networks for insulin secretion.

Identification and analysis of murine pancreatic islet enhancers

AIMS/HYPOTHESIS

The paucity of information on the epigenetic barriers that are blocking reprogramming protocols, and on what makes a beta cell unique, has hampered efforts to develop novel beta cell sources. Here, we aimed to identify enhancers in pancreatic islets, to understand their developmental ontologies, and to identify enhancers unique to islets to increase our understanding of islet-specific gene expression.

METHODS

We combined H3K4me1-based nucleosome predictions with pancreatic and duodenal homeobox 1 (PDX1), neurogenic differentiation 1 (NEUROD1), v-Maf musculoaponeurotic fibrosarcoma oncogene family, protein A (MAFA) and forkhead box A2 (FOXA2) occupancy data to identify enhancers in mouse islets.

RESULTS

We identified 22,223 putative enhancer loci in in vivo mouse islets. Our validation experiments suggest that nearly half of these loci are active in regulating islet gene expression, with the remaining regions probably poised for activity. We showed that these loci have at least nine developmental ontologies, and that islet enhancers predominately acquire H3K4me1 during differentiation. We next discriminated 1,799 enhancers unique to islets and showed that these islet-specific enhancers have reduced association with annotated genes, and identified a subset that are instead associated with novel islet-specific long non-coding RNAs (lncRNAs).

CONCLUSIONS/INTERPRETATIONS

Our results indicate that genes with islet-specific expression and function tend to have enhancers devoid of histone methylation marks or, less often, that are bivalent or repressed, in embryonic stem cells and liver. Further, we identify a subset of enhancers unique to islets that are associated with novel islet-specific genes and lncRNAs. We anticipate that these data will facilitate the development of novel sources of functional beta cell mass.

Transcriptome analysis of epigenetically modulated genome indicates signature genes in manifestation of type 1 diabetes and its prevention in NOD mice

Classic genetic studies implicated several genes including immune response genes in the risk of developing type 1 diabetes in humans. However, recent evidence including discordant diabetes incidence among monozygotic twins suggested a role for epigenetics in disease manifestation. NOD mice spontaneously develop type 1 diabetes like humans and serve as an excellent model system to study the mechanisms of type 1 diabetes as well as the efficacy of maneuvers to manipulate the disease. Using this preclinical model, we have recently demonstrated that pharmacological inhibition of histone deacetylases can lead to histone hyperacetylation, selective up-regulation of interferon-γ and its transactivator Tbx21/Tbet, and amelioration of autoimmune diabetes. In the current study, we show that chromatin remodeling can render splenocytes incapable of transferring diabetes into immunodeficient NOD.scid mice. To elucidate the underlying mechanisms of drug-mediated protection against type 1 diabetes, we performed global gene expression profiling of splenocytes using high throughput microarray technology. This unbiased transcriptome analysis unraveled the exaggerated expression of a novel set of closely related inflammatory genes in splenocytes of acutely diabetic mice and their repression in mice cured of diabetes by chromatin remodeling. Analysis of gene expression by qRT-PCR using RNA derived from spleens and pancreata of cured mice validated the suppression of most of these genes, indicating an inverse correlation between the high levels of these inflammatory genes and protection against diabetes in NOD mice. In addition, higher-level expression of genes involved in insulin sensitivity, erythropoiesis, hemangioblast generation, and cellular redox control was evident in spleens of cured mice, indicating their possible contribution to protection against type 1 diabetes. Taken together, these results are consistent with the involvement of epistatic mechanisms in the manifestation of autoimmune diabetes and further indicate the utility of chromatin remodeling in curing this complex autoimmune disorder.

The transcription factor Myt3 acts as a pro-survival factor in β-cells

Aims/Hypothesis

We previously identified the transcription factor Myt3 as specifically expressed in pancreatic islets. Here, we sought to determine the expression and regulation of Myt3 in islets and to determine its significance in regulating islet function and survival.

Methods

Myt3 expression was determined in embryonic pancreas and adult islets by qPCR and immunohistochemistry. ChIP-seq, ChIP-qPCR and luciferase assays were used to evaluate regulation of Myt3 expression. Suppression of Myt3 was used to evaluate gene expression, insulin secretion and apoptosis in islets.

Results

We show that Myt3 is the most abundant MYT family member in adult islets and that it is expressed in all the major endocrine cell types in the pancreas after E18.5. We demonstrate that Myt3 expression is directly regulated by Foxa2, Pdx1, and Neurod1, which are critical to normal β-cell development and function, and that Ngn3 induces Myt3 expression through alterations in the Myt3 promoter chromatin state. Further, we show that Myt3 expression is sensitive to both glucose and cytokine exposure. Of specific interest, suppressing Myt3 expression reduces insulin content and increases β-cell apoptosis, at least in part, due to reduced Pdx1, Mafa, Il-6, Bcl-xl, c-Iap2 and Igfr1 levels, while over-expression of Myt3 protects islets from cytokine induced apoptosis.

Conclusion/Interpretation

We have identified Myt3 as a novel transcriptional regulator with a critical role in β-cell survival. These data are an important step in clarifying the regulatory networks responsible for β-cell survival, and point to Myt3 as a potential therapeutic target for improving functional β-cell mass.

Identification of four genes required for mammalian blastocyst formation

RNA transcription, processing and translation are fundamental molecular processes required for development, growth and cell viability. Towards the functional annotation of the genome, we are engaged in a reverse genetic screen using mammalian preimplantation embryos as a model system. Here we report the essential function of four RNA processing/splicing factors (Sf3b14, Sf3b1, Rpl7l1, and Rrp7a) and show that each of these genes is required for blastocyst formation in the mouse. As very little information is known about these genes, we characterized their normal expression and localization in mouse embryos as well as phenotypic analysis of loss of function during preimplantation development. Functional knockdown of each gene product results in normal morula development but there is failure to form a blastocoel cavity or morphologically differentiated trophectoderm. We show that zygotic genome activation does occur as well as initial lineage specification in the absence of each factor. Consistent with a role in RNA splicing, we demonstrate that the absence of Sf3b14 and Sf3b1 in 8-cell and morula-stage embryos results in a specific reduction of intron containing transcripts, but no reduction single-exon genes. Taken together, we show critical developmental and molecular requirements of Sf3b14, Sf3b1, Rpl7l1, and Rrp7a during mammalian preimplantation.

Mouse strain specific gene expression differences for illumina microarray expression profiling in embryos

Background

In the field of mouse genetics the advent of technologies like microarray based expression profiling dramatically increased data availability and sensitivity, yet these advanced methods are often vulnerable to the unavoidable heterogeneity of in vivo material and might therefore reflect differentially expressed genes between mouse strains of no relevance to a targeted experiment. The aim of this study was not to elaborate on the usefulness of microarray analysis in general, but to expand our knowledge regarding this potential “background noise” for the widely used Illumina microarray platform surpassing existing data which focused primarily on the adult sensory and nervous system, by analyzing patterns of gene expression at different embryonic stages using wild type strains and modern transgenic models of often non-isogenic backgrounds.

Results

Wild type embryos of 11 mouse strains commonly used in transgenic and molecular genetic studies at three developmental time points were subjected to Illumina microarray expression profiling in a strain-by-strain comparison. Our data robustly reflects known gene expression patterns during mid-gestation development. Decreasing diversity of the input tissue and/or increasing strain diversity raised the sensitivity of the array towards the genetic background. Consistent strain sensitivity of some probes was attributed to genetic polymorphisms or probe design related artifacts.

Conclusion

Our study provides an extensive reference list of gene expression profiling background noise of value to anyone in the field of developmental biology and transgenic research performing microarray expression profiling with the widely used Illumina microarray platform. Probes identified as strain specific background noise further allow for microarray expression profiling on its own to be a valuable tool for establishing genealogies of mouse inbred strains.

Islet organogenesis, angiogenesis and innervation

The pancreas is characterized by a major component, an exocrine and ductal system involved in digestion, and a minor component, the endocrine islets represented by islet micro-organs that tightly regulate glucose homoeostasis. Pancreatic organogenesis is strictly co-ordinated by transcription factors that are expressed sequentially to yield functional islets capable of maintaining glucose homoeostasis. Angiogenesis and innervation complete islet development, equipping islets to respond to metabolic demands. Proper regulation of this triad of processes during development is critical for establishing functional islets.

Purines as potential morphogens during embryonic development

Components of purinergic signalling are expressed in the early embryo raising the possibility that ATP, ADP and adenosine may contribute to the mechanisms of embryonic development. We summarize the available data from four developmental models—mouse, chick, Xenopus and zebrafish. While there are some notable examples where purinergic signalling is indeed important during development, e.g. development of the eye in the frog, it is puzzling that deletion of single components of purinergic signalling often results in rather minor developmental phenotypes. We suggest that a key step in further analysis is to perform combinatorial alterations of expression of purinergic signalling components to uncover their roles in development. We introduce the concept that purinergic signalling could create novel morphogenetic fields to encode spatial location via the concentration of ATP, ADP and adenosine. We show that using minimal assumptions and the known properties of the ectonucleotidases, complex spatial patterns of ATP and adenosine can be set up. These patterns may provide a new way to assess the potential of purinergic signalling in developmental processes.

BioTapestry: a tool to visualize the dynamic properties of gene regulatory networks

BioTapestry is an open source, freely available software tool that has been developed to handle the -challenges of modeling genetic regulatory networks (GRNs). Using BioTapestry, a researcher can -construct a network model and use it to visualize and understand the dynamic behavior of a complex, spatially and temporally distributed GRN. Here we provide a step-by-step example of a way to use BioTapestry to build a GRN model and discuss some common issues that can arise during this process.

Targeting SOX17 in human embryonic stem cells creates unique strategies for isolating and analyzing developing endoderm

Human embryonic stem cells (hESCs) can provide insights into development of inaccessible human tissues such as embryonic endoderm. Progress in this area has been hindered by a lack of methods for isolating endodermal cells and tracing fates of their differentiated progeny. By using homologous recombination in human ESCs, we inserted an enhanced green fluorescent protein (eGFP) transgene into the SOX17 locus, a postulated marker of human endoderm. FACS purification and gene expression profiling confirmed that SOX17(+)-hESC progeny expressed endodermal markers and unveiled specific cell surface protein combinations that permitted FACS-based isolation of primitive gut tube endodermal cells produced from unmodified human ESCs and from induced pluripotent stem cells (iPSC). Differentiating SOX17(+) endodermal cells expressed markers of liver, pancreas, and intestinal epithelium in vitro and gave rise to endodermal progeny in vivo. Thus, prospective isolation, lineage tracing, and developmental studies of SOX17(+) hESC progeny have revealed fundamental aspects of human endodermal biology.

An expression profile analysis of ES cell-derived definitive endodermal cells and Pdx1-expressing cells

Background

We developed an efficient in vitro method to differentiate mouse ES cells into the definitive endoderm (DE) and then Pdx1-expressing pancreatic lineages using mesodermal-derived supporting cells, M15. Using this method, resulting ES cell-derived DE and Pdx1-expressing cells were isolated by cell sorting, and their gene expression profiles were investigated with DNA microarray. Genes that were specifically expressed in DE and/or in Pdx1-expressing cells were extracted and their expression patterns in normal embryonic development were studied.

Results

Genes whose expression increased in DE and Pdx1 positive cells compared to the undifferentiated ES cells were chosen and in situ hybridizations were performed. Out of 54 genes examined, 27 were expressed in the DE of E8.5 mouse embryos and 15 genes were expressed in distinct domains in the pancreatic buds of E14.5 embryos. Among those genes expressed were Foxq1, CpM, Foxp4, Pcdh1, and Zmiz1, which were previously reported in other endodermal tissues. Genes, such as Parm1, Tmem184a, Hipk2 and Sox4 were reported to be expressed during early pancreatic development. Nptx2, C2cd4b, Tcf7l2 and Kiss1r were reported to be associated with beta cell or pancreatic functions in the adult. Akr1c19, Aebp2, Pbxip1 and Creb3l1, were novel and have not been described as being expressed either in DE or the pancreas.

Conclusions

We identified 27 genes, including 4 novel genes expressed in DE and pancreatic progenitor cells during normal development using an ES cell in vitro differentiation system. These results showed that DE cells and Pdx1/GFP-expressing cells obtained from our M15 based differentiation method mimic cells during the normal developmental processes. Additionally, ES cells are an excellent model for studies of early developmental processes.

Gill transcriptome response to changes in environmental calcium in the green spotted puffer fish

Background

Calcium ion is tightly regulated in body fluids and for euryhaline fish, which are exposed to rapid changes in environmental [Ca²⁺], homeostasis is especially challenging. The gill is the main organ of active calcium uptake and therefore plays a crucial role in the maintenance of calcium ion homeostasis. To study the molecular basis of the short-term responses to changing calcium availability, the whole gill transcriptome obtained by Super Serial Analysis of Gene Expression (SuperSAGE) of the euryhaline teleost green spotted puffer fish, Tetraodon nigroviridis, exposed to water with altered [Ca²⁺] was analysed.

Results

Transfer of T. nigroviridis from 10 ppt water salinity containing 2.9 mM Ca²⁺ to high (10 mM Ca²⁺ ) and low (0.01 mM Ca²⁺) calcium water of similar salinity for 2-12 h resulted in 1,339 differentially expressed SuperSAGE tags (26-bp transcript identifiers) in gills. Of these 869 tags (65%) were mapped to T. nigroviridis cDNAs or genomic DNA and 497 (57%) were assigned to known proteins. Thirteen percent of the genes matched multiple tags indicating alternative RNA transcripts. The main enriched gene ontology groups belong to Ca²⁺ signaling/homeostasis but also muscle contraction, cytoskeleton, energy production/homeostasis and tissue remodeling. K-means clustering identified co-expressed transcripts with distinct patterns in response to water [Ca²⁺] and exposure time.

Conclusions

The generated transcript expression patterns provide a framework of novel water calcium-responsive genes in the gill during the initial response after transfer to different [Ca²⁺]. This molecular response entails initial perception of alterations, activation of signaling networks and effectors and suggests active remodeling of cytoskeletal proteins during the initial acclimation process. Genes related to energy production and energy homeostasis are also up-regulated, probably reflecting the increased energetic needs of the acclimation response. This study is the first genome-wide transcriptome analysis of fish gills and is an important resource for future research on the short-term mechanisms involved in the gill acclimation responses to environmental Ca²⁺ changes and osmoregulation.

Locus co-occupancy, nucleosome positioning, and H3K4me1 regulate the functionality of FOXA2-, HNF4A-, and PDX1-bound loci in islets and liver

The liver and pancreas share a common origin and coexpress several transcription factors. To gain insight into the transcriptional networks regulating the function of these tissues, we globally identify binding sites for FOXA2 in adult mouse islets and liver, PDX1 in islets, and HNF4A in liver. Because most eukaryotic transcription factors bind thousands of loci, many of which are thought to be inactive, methods that can discriminate functionally active binding events are essential for the interpretation of genome-wide transcription factor binding data. To develop such a method, we also generated genome-wide H3K4me1 and H3K4me3 localization data in these tissues. By analyzing our binding and histone methylation data in combination with comprehensive gene expression data, we show that H3K4me1 enrichment profiles discriminate transcription factor occupied loci into three classes: those that are functionally active, those that are poised for activation, and those that reflect pioneer-like transcription factor activity. Furthermore, we demonstrate that the regulated presence of H3K4me1-marked nucleosomes at transcription factor occupied promoters and enhancers controls their activity, implicating both tissue-specific transcription factor binding and nucleosome remodeling complex recruitment in determining tissue-specific gene expression. Finally, we apply these approaches to generate novel insights into how FOXA2, PDX1, and HNF4A cooperate to drive islet- and liver-specific gene expression.

A seriation approach for visualization-driven discovery of co-expression patterns in Serial Analysis of Gene Expression (SAGE) data

Background

Serial Analysis of Gene Expression (SAGE) is a DNA sequencing-based method for large-scale gene expression profiling that provides an alternative to microarray analysis. Most analyses of SAGE data aimed at identifying co-expressed genes have been accomplished using various versions of clustering approaches that often result in a number of false positives.

Principal Findings

Here we explore the use of seriation, a statistical approach for ordering sets of objects based on their similarity, for large-scale expression pattern discovery in SAGE data. For this specific task we implement a seriation heuristic we term ‘progressive construction of contigs’ that constructs local chains of related elements by sequentially rearranging margins of the correlation matrix. We apply the heuristic to the analysis of simulated and experimental SAGE data and compare our results to those obtained with a clustering algorithm developed specifically for SAGE data. We show using simulations that the performance of seriation compares favorably to that of the clustering algorithm on noisy SAGE data.

Conclusions

We explore the use of a seriation approach for visualization-based pattern discovery in SAGE data. Using both simulations and experimental data, we demonstrate that seriation is able to identify groups of co-expressed genes more accurately than a clustering algorithm developed specifically for SAGE data. Our results suggest that seriation is a useful method for the analysis of gene expression data whose applicability should be further pursued.

Expression of Groucho/TLE proteins during pancreas development

Background

The full-length mammalian homologs of groucho, Tle1, 2, 3, and 4, act as transcriptional corepressors and are recruited by transcription factors containing an eh1 or WRPW/Y domain. Many transcription factors critical to pancreas development contain a Gro/TLE interaction domain and several have been shown to require Gro/TLE interactions for proper function during neuronal development. However, a detailed analysis of the expression patterns of the Gro/TLE proteins in pancreas development has not been performed. Moreover, little is known about the ability of Gro/TLE proteins to interact with transcription factors in the pancreas.

Results

We describe the expression of Gro/TLE family members, and of 34 different transcription factors that contain a Gro/TLE interaction motif, in the pancreas utilizing nine SAGE libraries created from the developing and adult pancreas, as well as the GenePaint database. Next, we show the dynamic expression of Tle1, 2, 3, 4, 5 and 6 during pancreas development by qRT-PCR. To further define the cell-type specificity of the expression of these proteins we use immunofluorescence to co-localize them with Pdx1 at embryonic day 12.5 (E12.5), Ngn3 at E14.5, Pdx1, Nkx2-2, Insulin, Glucagon, Pancreatic polypeptide and Somatostatin at E18.5, as well as Insulin and Glucagon in the adult. We then show that Tle2 can interact with Nkx2-2, Hes1, Arx, and Nkx6-1 which are all critical factors in pancreas development. Finally, we demonstrate that Tle2 modulates the repressive abilities of Arx in a β-cell line.

Conclusion

Although Tle1, 2, 3, and 4 show overlapping expression in pancreatic progenitors and in the adult islet, the expression of these factors is restricted to different cell types during endocrine cell maturation. Of note, Tle2 and Tle3 are co-expressed with Gro/TLE interaction domain containing transcription factors that are essential for endocrine pancreas development. We further demonstrate that Tle2 can interact with several of these factors and that Tle2 modulate Arx's repressive activity. Taken together our studies suggest that Gro/TLE proteins play a role in the repression of target genes during endocrine cell specification.