<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">JDM</journal-id><journal-title-group><journal-title>Journal of Diabetes Mellitus</journal-title></journal-title-group><issn pub-type="epub">2160-5831</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jdm.2021.111001</article-id><article-id pub-id-type="publisher-id">JDM-107018</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Medicine&amp;Healthcare</subject></subj-group></article-categories><title-group><article-title>
 
 
  Gene Expression Profiling of the Mouse Pancreas during the Secondary Transition in the Organogenesis of the Pancreatic Gland*
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Stefanie</surname><given-names>J. Willmann</given-names></name><xref ref-type="aff" rid="aff1"><sub>1</sub></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><label>1</label><addr-line>TUM, Munich, Germany</addr-line></aff><pub-date pub-type="epub"><day>03</day><month>02</month><year>2021</year></pub-date><volume>11</volume><issue>01</issue><fpage>1</fpage><lpage>9</lpage><history><date date-type="received"><day>23,</day>	<month>November</month>	<year>2020</year></date><date date-type="rev-recd"><day>31,</day>	<month>January</month>	<year>2021</year>	</date><date date-type="accepted"><day>3,</day>	<month>February</month>	<year>2021</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  Diabetes mellitus
   is a chronic disease that impacts the homeostasis of blood sugar levels caused by loss or defect of insulin-producing β-cells in the Islets of Langerhans. 
  Type 1 diabetes
   (
  T1D
  )
   
  is caused by auto-immune mediated destruction of β-cells, 
  whereas in T2D
  , insulin is produced but used inefficiently.
   T2D
   accounts for 90% of people with diabetes worldwide (
  WHO
   
  1999)
   and is 
  the fastest increasing disease worldwide (
  https://diabetesatlas.org/en/
  )
  . For an improved understanding of the pathomechanism of diabetes, profound knowledge of pancreas organogenesis and the associated gene regulatory networks is required. Therefore, we dissected and profiled the pancreatic endodermal and non-endodermal compartment between the embryonic stages (E) 12.5 and E 15.5 when progenitor cells commit to their different pancreatic lineages. Our associated study mined the global mRNA expression profile to increase the understanding of the secondary transition, endodermal-non-endodermal tissue interaction, and diabetic-related gene regulation. Furthermore, we validated 635 regulated pancreatic genes using the publicly available GenePaint.org, respective gp3.mpg.de
   to evaluate genes associated with genetic variants in Single-nucleotide polymorphism (SNP) related to 
  T2D
  .
 
</p></abstract><kwd-group><kwd>Gene Expression Profile</kwd><kwd> Pancreas Organogenesis</kwd><kwd> Single-Nucleotide  Polymorphism</kwd><kwd> Type 1 Diabetes</kwd><kwd> Type 2 Diabetes</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Experimental Design, Materials and Methods</title><sec id="s1_1"><title>1.1. Dissection, Isolation, and Processing</title><p>Pancreatic tissue was isolated from FVF+ embryos at E 12.5 to E 15.5 in triplicates [<xref ref-type="bibr" rid="scirp.107018-ref1">1</xref>]. A single-cell suspension was obtained by trypsinization with 0.25% Trypsin and the specific incubated pancreatic tissue resuspension. Fluorescent-Activated cell sorting (FACs) was used to separate the FVF+ endodermal and FVF− non-endodermal cells. Total RNA was isolated using Trizol (Qiagen) and purified by using the RNeasy Mini kit (Qiagen) according to the manufacturer’s protocol. The Agilent 2100 Bioanalyzer had assessed RNA quality. The total RNA of three biological replicates for each embryonic stage was used to generate cRNA by two rounds of transcription. In general, a poly (dT) primer was utilized to synthesize cDNA from total RNA. In the next step, the cDNA amplified cRNA through T7 polymerase, using the T7 promoter sequence of the poly (dT) primer in the first round by random priming. Next, the cDNA was employed for the transcription into biotinylated cRNA, which hybridized to the 169 Affymetrix&#174; array format according to the manufacturer protocol (GeneChip&#174; Gene mouse 1.0 ST Array Card). Subsequently, processing of the Fluidics script FS450_0007 and scanning (GeneChip&#174; Scanner 3000 7G Whole-Genome Association System) further lead to the probe set outcome.</p><p>The Affymetrix&#174; Expression console normalized the robust multichip-analyses (RMA) on gene-level (ratio &gt; 200 compared endodermal and non-endodermal). Comprehensive R based microarray analyses web frontend (CARMAweb) inducted gene-wise testing by Limma t-test and Benjamini-Hochberg multiple testing correction with a false discovery rate (FDR &lt; 5%) offered statistically significant differential expression of 2921 probe sets [<xref ref-type="bibr" rid="scirp.107018-ref2">2</xref>]. Data were submitted to the Gene Expression Omnibus (GSE66856).</p></sec><sec id="s1_2"><title>1.2. Gene ontology (GO) Term Analysis and Validation of Gene Expression by Gene Paint in Silico in Situ</title><p>Subsequently, GO term enrichment, annotated through association of InterPro records with GO term, grouped 635 genes in the following clusters: extracellular/located at plasma membrane (304 genes), signaling molecules (222 genes), transcription factors (TF; 178 genes), cilium related (40 genes) and previously not described (42 genes). Afterwards, pancreatic-specific gene expression and localization were analyzed using the database GenePaint.org, respective gp3.mpg.de, which offers in silico in situs of whole mount embryosat E1 4.5. Pancreatic genes had been arranged according to their expression pattern in mesenchyme (non-epithelial), tip (exocrine), trunk (endocrine) and epithelium (pancreatic multipotent cells).</p></sec><sec id="s1_3"><title>1.3. Analysis of Mouse Orthologous Genes with Associated Single-Nucleotide Polymorphism (SNP) in T1D and T2D GWAS Hits</title><p>SNP of genes associated with T1D and T2D had been analyzed for their expression level during secondary transition and depicted in a heatmap (<xref ref-type="fig" rid="fig1">Figure 1</xref>). The heatmap represents the relative fold change in gene expression during the secondary transition stage between the non-endodermal and endodermal tissue compartments. The clustering of the genes was accomplished through extraction of gene symbols by Marullo et al. and Hakonarson et al. [<xref ref-type="bibr" rid="scirp.107018-ref3">3</xref>] [<xref ref-type="bibr" rid="scirp.107018-ref4">4</xref>]. Based on gene expression levels in the two tissue compartments we identified different clusters of the selected genes: transcriptionally regulated in both tissue populations or elevated either in the endodermal or the non-endodermal tissue compartment.</p><p>Heatmap shows clustering of 153 T1D- and T2D-associated genes according to their expression level during lineage commitment between E12.5 - 15.5. Predominantly, maturity-onset diabetes of the young (MODY) genes are upregulated in the endodermal compartment during pancreas organogenesis (Gck, Pdx1, NeuroD, Hnf1α, Hnf1&#223;). The black arrow marks the gene Ins2.</p><p>A subset of these T2D-associated genes which are regulated in our mRNA expression profile is illustrated in <xref ref-type="table" rid="table1">Table 1</xref>. Additional information of the mRNA localization (endodermal/non-endodermal), gene locus, chromosomal localization, the associated SNP and the disease as well as the related publications is provided as follows [<xref ref-type="bibr" rid="scirp.107018-ref4">4</xref>] - [<xref ref-type="bibr" rid="scirp.107018-ref12">12</xref>].</p><p>The GenePaint.org, respective gp3.mpg.dein silico in situ platform confirmed expression of mouse orthologues with associated SNPs in T2D shown in <xref ref-type="table" rid="table1">Table 1</xref>. Additional information about mRNA localization of the selected genes is obtained by inspection of the in situ whole mount embryos (see <xref ref-type="fig" rid="fig2">Figure 2</xref>).</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Genes associated with SNP in GWAS studies and related to T1D or T2D</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  colspan="6"  >Genetic variants associated with T2D or glycemic traits</th></tr></thead><tr><td align="center" valign="middle" >Localization</td><td align="center" valign="middle" >Locus name</td><td align="center" valign="middle" >Chromosome</td><td align="center" valign="middle" >associated SNP</td><td align="center" valign="middle" >established loci</td><td align="center" valign="middle" >Literature</td></tr><tr><td align="center" valign="middle" >non- and endodermal</td><td align="center" valign="middle" >TCF7L2</td><td align="center" valign="middle" >10</td><td align="center" valign="middle" >rs4506565, rs7903146, rs12243326</td><td align="center" valign="middle" >T2DM</td><td align="center" valign="middle" >Dupuis et al., 2012; Saxena et al., 2010</td></tr><tr><td align="center" valign="middle" >endodermal</td><td align="center" valign="middle" >KCNQ1</td><td align="center" valign="middle" >11</td><td align="center" valign="middle" >rs2237892, rs231362</td><td align="center" valign="middle" >T2DM</td><td align="center" valign="middle" >Yasuda et al., 2008; Voight et al., 2010</td></tr><tr><td align="center" valign="middle" >endodermal</td><td align="center" valign="middle" >PROX1</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >rs340874</td><td align="center" valign="middle" >T2DM</td><td align="center" valign="middle" >Dupuis et al., 2012</td></tr><tr><td align="center" valign="middle" >endodermal</td><td align="center" valign="middle" >HNF1β</td><td align="center" valign="middle" >17</td><td align="center" valign="middle" >rs757210, rs4430796, rs7501939</td><td align="center" valign="middle" >T2DM</td><td align="center" valign="middle" >Pearson et al., 2007; Winckler et al., 2007; Wang et al., 2014</td></tr><tr><td align="center" valign="middle" >endodermal</td><td align="center" valign="middle" >KCNJ11</td><td align="center" valign="middle" >11</td><td align="center" valign="middle" >rs5219, rs757110</td><td align="center" valign="middle" >T2DM</td><td align="center" valign="middle" >Gloyn et al., 2003</td></tr><tr><td align="center" valign="middle" >endodermal</td><td align="center" valign="middle" >GCK</td><td align="center" valign="middle" >7</td><td align="center" valign="middle" >rs1799881, rs6975024, rs4607517, rs3757840</td><td align="center" valign="middle" >T2DM</td><td align="center" valign="middle" >Dupuis et al., 2012; Soranzo et al., 2010; Scott et al., 2012</td></tr><tr><td align="center" valign="middle" >endodermal</td><td align="center" valign="middle" >FOXA2</td><td align="center" valign="middle" >20</td><td align="center" valign="middle" >rs6113722, rs6048205</td><td align="center" valign="middle" >T2DM</td><td align="center" valign="middle" >Scott et al., 2012; Manning et al., 2012</td></tr><tr><td align="center" valign="middle" >endodermal</td><td align="center" valign="middle" >PDX1</td><td align="center" valign="middle" >13</td><td align="center" valign="middle" >rs11619319, rs2293941</td><td align="center" valign="middle" >T2DM</td><td align="center" valign="middle" >Scott et al., 2012; Manning et al., 2012</td></tr></tbody></table></table-wrap><p>Sagittal sections of whole mount embryos and magnification of the pancreatic region are shown. The genes associated with SNP in T2D as Kcnq1, Prox1, Hnf1&#223;, Gck, Foxa2 and Pdx1 exhibit epithelial in situ expression pattern.</p></sec></sec><sec id="s2"><title>2. Conclusion</title><p>The study comprehensively underpins the association of de-/regulation of developmental genes in the mouse embryo with distinct mRNA expression, validated by genepaint, respective gp3.mpg.de and genetic variants in SNP connected to the disease diabetes mellitus. Thus, the database may serve as a resource for further evaluation of mouse genes and their respective human orthologues in developmental stages and in a next step in different diseases. Further genes in the context of major health diseases as described in this study on Diabetes mellitus, Diabetes mellitus related corornary artery disease (CAD), cancers and respiratory diseases might shed light on disease mechanism and will improve therapeutic approaches for a better health outcome of the patients’ need.</p></sec><sec id="s3"><title>Acknowledgements</title><p>We thank the Helmholtz Alliance ICEMED and Helmholtz Portfolio Theme “Metabolic Dysfunction and Common Disease” (grants to HL and JB) and the European Research Council (Latent Causes: 259294; grant to FJT) for funding. Additionally, the HumEn project was financially supported from the European Union’s Seventh Framework Programme for Research, Technological Development and Demonstration under grant agreement No. 602587 (http://www.cordis.europa.eu/). S.W. thanks the Helmholtz Society, German Research Foundation and German Center for Diabetes Research (DZD e.V.) for financial and arranging support and the Federal Ministry of Education and Research [Virtual Liver: FKZ 0315752 and LungSys: FKZ 0316042I] for the assistance.</p></sec><sec id="s4"><title>Direct Link to the Deposited Data</title><p>Data is available through the NCBI database under the following link http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=gnubwesiphanjiv&amp;acc=GSE66856.</p><p>Statistical analyses of the dataset were accomplished by using the GEO2R platform.</p></sec><sec id="s5"><title>Conflict of Interest</title><p>The author declares that there is no conflict of interest on the work presented in the paper.</p></sec><sec id="s6"><title>Cite this paper</title><p>Willmann, S.J. (2021) Gene Expression Profiling of the Mouse Pancreas during the Secondary Transition in the Organogenesis of the Pancreatic Gland. Journal of Diabetes Mellitus, 11, 1-9. https://doi.org/10.4236/jdm.2021.111001</p></sec><sec id="s7"><title>Abbreviations</title><p>T1D—type 1 diabetes</p><p>T2D—type 2 diabetes</p><p>WHO—world health organisation</p><p>E—embryonic stage</p><p>mRNA—messenger ribonucleic acid</p><p>SNP—single-nucleotide polymorphism</p><p>GWAS—genome-wide association study</p><p>FACs—fluorescent-activated cell sorting</p><p>FVF—Foxa2-Venus-Fusion</p><p>cDNA—complementary deoxynucleic acid</p><p>GRN—gene regulatory network</p><p>RNA—ribonucleic acid</p><p>cRNA—complementary ribonucleic acid</p><p>RMA—robust multichip-analyses</p><p>CARMA—Comprehensive R based microarray analyses web frontend</p><p>GO—gene ontology</p><p>TF—transcription factor</p></sec><sec id="s8"><title>Specifications Table</title></sec><sec id="s9"><title>Value of the Data</title><p>&#183; Demonstrates the value of the Foxa2-Venus fusion reporter mouse line for diabetic related research [<xref ref-type="bibr" rid="scirp.107018-ref1">1</xref>].</p><p>&#183; Provides a rich data set of differentially regulated genes within the endodermal and non-endodermal tissue compartments during the secondary transition of the pancreas organogenesis.</p><p>&#183; Offers temporal and spatial progression of mouse orthologues T1D and T2D associated genes.</p><p>&#183; Evaluates database-dependent gene regulatory network (GRN) of genome-wide association studies (GWAS)-annotated diabetes genes.</p></sec><sec id="s10"><title>NOTES</title></sec></body><back><ref-list><title>References</title><ref id="scirp.107018-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Manning, M.-F., Hivert, Scott, R.A., Grimsby, J.L., Bouatia-Naji, N., Chen, H., et al. (2012) A Genome-Wide Approach Accounting for Body Mass Index Identifies Genetic Variants Influencing Fasting Glycemic Traits and Insulin Resistance. Nature Genetics, 44, 659-669. https://doi.org/10.1038/ng.2274</mixed-citation></ref><ref id="scirp.107018-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Scott, R.A., Lagou, V., Welch, R.P., Wheeler, E., Montasser, M.E., Luan, J., et al. (2012) Large-Scale Association Analyses Identify New Loci Influencing Glycemic Traits and Provide Insight into the Underlying Biological Pathways. Nature Genetics, 44, 991-1005. https://doi.org/10.1038/ng.2385</mixed-citation></ref><ref id="scirp.107018-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Gloyn, A.L. (2003) The Search for Type 2 Diabetes Genes. Ageing Research Reviews, 2, 111-127. https://doi.org/10.1016/S1568-1637(02)00061-2</mixed-citation></ref><ref id="scirp.107018-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Winckler, W., Weedon, M.N., Graham, R.R., McCarroll, S.A., Purcell, S., Almgren, P., et al. (2007) Evaluation of Common Variants in the Six Known Maturity-Onset Diabetes of the Young (MODY) Genes for Association with Type 2 Diabetes. Diabetes, 56, 685-693. https://doi.org/10.2337/db06-0202</mixed-citation></ref><ref id="scirp.107018-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Pearson, E.R., Donnelly, L.A., Kimber, C., Whitley, A., Doney, A.S.F., Mccarthy, M.I., et al. (2007) Variation in TCF7L2 Influences Therapeutic Response to Sulfonylureas. Diabetes, 56, 2178-2182. https://doi.org/10.2337/db07-0440</mixed-citation></ref><ref id="scirp.107018-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Voight, B.F., Scott, L.J., Steinthorsdottir, V., Morris, A.P., Dina, C., Welch, R.P., et al. (2010) Twelve Type 2 Diabetes Susceptibility Loci Identified through Large-Scale Association Analysis. Nature Genetics, 42, 579-589. https://doi.org/10.1038/ng.609</mixed-citation></ref><ref id="scirp.107018-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Yasuda, K., Miyake, K., Horikawa, Y., Hara, K., Osawa, H., Furuta, H., et al. (2008) Variants in KCNQ1 Are Associated with Susceptibility to Type 2 Diabetes Mellitus. Nature Genetics, 40, 1092-1097. https://doi.org/10.1038/ng.207</mixed-citation></ref><ref id="scirp.107018-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Saxena, R., Hivert, M.-F., Langenberg, C., Tanaka, T., Pankow, J.S., Vollenweider, P., et al. (2010) Genetic Variation in GIPR Influences the Glucose and Insulin Responses to an Oral Glucose Challenge. Nature Genetics, 42, 142-148. https://doi.org/10.1038/ng.521</mixed-citation></ref><ref id="scirp.107018-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Dupuis, J., Langenberg, C., Prokopenko, I., Saxena, R., Soranzo, N., Jackson, A.U., et al. (2010) New Genetic Loci Implicated in Fasting Glucose Homeostasis and Their Impact on Type 2 Diabetes Risk. Nature Genetics, 42, 105-116. https://doi.org/10.1038/ng.520</mixed-citation></ref><ref id="scirp.107018-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Marullo, L., El-Sayed Moustafa, J.S. and Prokopenko, I. (2014) Insights into the Genetic Susceptibility to Type 2 Diabetes from Genome-Wide Association Studies of Glycaemic Traits. Current Diabetes Reports, 14, Article Number: 551. https://doi.org/10.1007/s11892-014-0551-8</mixed-citation></ref><ref id="scirp.107018-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Rainer, J., Sanchez-Cabo, F., Stocker, G., Sturn, A. and Trajanoski, Z. (2006) CARMAweb: Comprehensive R- and Bioconductor-Based Web Service for Microarray Data Analysis. Nucleic Acids Research, 34, 498-503. https://doi.org/10.1093/nar/gkl038</mixed-citation></ref><ref id="scirp.107018-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Burtscher, I., Barkey, W. and Lickert, H. (2013) Foxa2-Venus Fusion Reporter Mouse Line Allows Live-Cell Analysis of Endoderm-Derived Organ Formation. Genesis, 51, 596-604. https://doi.org/10.1002/dvg.22404</mixed-citation></ref></ref-list></back></article>