# MafB-dependent neurotransmitter signaling promotes β cell migration in the developing pancreas

Sara Bsharat, Emanuela Monni, Tania Singh, Jenny K Johansson, Kavya Achanta, Ludivine Bertonnier-Brouty, Anja Schmidt-Christensen, Dan Holmberg, Zaal Kokaia, Rashmi B Prasad, et al.

<figure><img src="image_1.png" /></figure>

► To cite this version:

Sara Bsharat, Emanuela Monni, Tania Singh, Jenny K Johansson, Kavya Achanta, et al. MafB-dependent neurotransmitter signaling promotes β cell migration in the developing pancreas. Development (Cambridge, England), 2023, 150 (6), 10.1242/dev.201009 . hal-04973778

HAL Id: hal-04973778

https://hal.science/hal-04973778v1

Submitted on 3 Mar 2025

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
RESEARCH ARTICLE

# MafB-dependent neurotransmitter signaling promotes β cell migration in the developing pancreas

Sara Bsharat<fnref n="1" />,<fnref n="2" />, Emanuela Monni<fnref n="1" />, Tania Singh<fnref n="1" />,<fnref n="2" />, Jenny K. Johansson<fnref n="1" />,<fnref n="2" />, Kavya Achanta<fnref n="1" />,<fnref n="2" />, Ludivine Bertonnier-Brouty<fnref n="1" />,<fnref n="2" />, Anja Schmidt-Christensen<fnref n="2" />, Dan Holmberg<fnref n="2" />, Zaal Kokaia<fnref n="1" />, Rashmi B. Prasad<fnref n="2" />,\* and Isabella Artner<fnref n="1" />,<fnref n="2" />,\*

ABSTRACT

Hormone secretion from pancreatic islets is essential for glucose homeostasis, and loss or dysfunction of islet cells is a hallmark of type 2 diabetes. Maf transcription factors are crucial for establishing and maintaining adult endocrine cell function. However, during pancreas development, MafB is not only expressed in insulin- and glucagon-producing cells, but also in Neurog3⁺ endocrine progenitor cells, suggesting additional functions in cell differentiation and islet formation. Here, we report that MafB deficiency impairs β cell clustering and islet formation, but also coincides with loss of neurotransmitter and axon guidance receptor gene expression. Moreover, the observed loss of nicotinic receptor gene expression in human and mouse β cells implied that signaling through these receptors contributes to islet cell migration/formation. Inhibition of nicotinic receptor activity resulted in reduced β cell migration towards autonomic nerves and impaired β cell clustering. These findings highlight a novel function of MafB in controlling neuronal-directed signaling events required for islet formation.

KEY WORDS: MafB, Pancreas development, Islet formation, Nicotinic acetylcholine receptor, Axon guidance receptor

INTRODUCTION

Pancreatic islets consist predominantly of insulin-producing β and glucagon-producing α cells which are essential for glucose homeostasis. Loss or dysfunction of β cells results in diabetes mellitus. Islet transplantation has proven to be an effective therapy for type 1 diabetes (Triolo and Bellin, 2021). However, donor material is scarce, which makes the production of *in vitro* differentiated islet cells a necessity for cell replacement treatment. Recent reports demonstrated the feasibility of this approach, but increasing evidence suggests that proper glucose control is only achieved if all islet cell types are present and work in concert with one another (Johnston et al., 2016). A thorough understanding of islet development will be required to recreate islet architecture *in vitro* and successfully generate islets for transplantation.

During mouse pancreas development, endocrine progenitor cells initially appear at embryonic day (E) 12.5, but hormone cell production is observed throughout development, with a distinct peak at E15.5. In contrast, postnatal increases in islet cell numbers are caused by proliferation (Bonner-Weir et al., 2016), which decreases in the adult. Mouse islets form by budding of endocrine progenitor cells from pancreatic ducts (Sharon et al., 2019) and aggregation of migrating endocrine cells (Sznurkowska et al., 2020), as signaling cues from the pancreatic mesenchyme (Pauerstein et al., 2017), blood vessels (Staels et al., 2019) and sympathetic neuronal cells (Borden et al., 2013) promote endocrine cell differentiation, with adrenergic signaling from sympathetic nerves promoting β cell expansion and islet formation (Borden et al., 2013).

The MafA and MafB transcription factors are crucial regulators of β cell function and islet immune homeostasis (Artner et al., 2010; Singh et al., 2018; Zhang et al., 2005). MafA also regulates responsiveness of pancreatic islets to autonomic signals in adult islets (Ganic et al., 2016). During pancreas development, the closely related MafB transcription factor is required for α and β cell terminal differentiation, a function that is postnatally taken over by MafA (Artner et al., 2007, 2010; Cyphert et al., 2019; Nishimura et al., 2008). MafB is expressed in a portion of Neurog3⁺ endocrine progenitor cells (Artner et al., 2006), suggesting that MafB may control processes in endocrine cell specification and differentiation before terminal differentiation of endocrine cells.

Here, we show that loss of MafB in developing α and β cells disrupts islet and duct formation. Gene expression analysis revealed that *MafB*-deficient (*MafB*⁻/⁻) endocrine progenitor cells have reduced expression of several neurotransmitter and axon guidance receptor genes. Single cell sequencing further showed that neurotransmitter and axon guidance receptors were not only expressed in mouse islet β cells, but also co-expressed with MAFB in developing human β cells. In addition, expression of these neurotransmitter receptor genes was reduced in mouse and human *MafB*⁻/⁻ β progenitor cells. Importantly, inhibiting nicotinic acetylcholine receptor (nAchR) function in developing β cells impaired β cell migration and clustering, demonstrating a novel function of nAchRs in islet formation.

RESULTS

Loss of MafB in developing endocrine cells disrupts islet and ductal morphogenesis

MafB is expressed in α and β cell progenitor and mature α cells in the mouse pancreas. Loss of MafB results in 50% reduction of α and β cell mass at E18.5 (Fig. S1A), while endocrine progenitor cells are present (Artner et al., 2007; Nishimura et al., 2008). To determine whether loss of MafB impairs pancreas/islet morphology, wild-type and *MafB*⁻/⁻ (GFP knock-in; Blanchi

<figure><img src="image_3.png" /></figure>

<footnote n="1">¹Lund Stem Cell Center, Lund University, Klinikgatan 26, 22184, Lund, Sweden.</footnote>

<footnote n="2">²Lund University Diabetes Center, Jan Waldenströms gata 35, 21428, Malmö, Sweden.</footnote>

<footnote>*Authors for correspondence (isabella.artner@med.lu.se; rashmi.prasad@med.lu.se)</footnote>

<footnote>R.B.P., 0000-0002-4400-6741; I.A., 0000-0002-2895-5664</footnote>

<footnote>This is an Open Access article distributed under the terms of the Creative Commons Attribution
License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use,
distribution and reproduction in any medium provided that the original work is properly attributed.</footnote>
et al., 2003) E12.5 and E18.5 pancreata were stained for insulin (β cells) and mucin (ductal cells). 3D optical projection tomography and 2D immunohistochemical analysis showed that wild-type β cells formed large clusters around the developing mucin/ductal network at E18.5 (Fig. 1A) and that *MafB*⁻⁻ islet cells were located closer to the ductal epithelium (Fig. 1B). Moreover, the organization of the ductal network was impaired in *MafB*⁻⁻ pancreata, as only few larger ductal structures (Fig. 1A; Fig. S1B) were detected while smaller ductal lumen were present (Fig. S1B, red label), resulting in the total ductal volume not being impaired

(Fig. S1C). This was confirmed by immunohistochemical analysis, which showed that the ductal luminal area was not significantly different in *MafB*⁻⁻ (Fig. S1D), but a significantly larger portion of the ductal structures in contact with islets were small (>$10 µm^2$) and the luminal structures had significantly fewer branching points (Fig. 1C-E) than wild type. The average number of ductal branching points did not differ between wild-type and *MafB*⁻⁻ exocrine tissue (Fig. 1D) and we did not observe morphological changes in ductal volume and organization at earlier time points (E12.5, Fig. S1E), suggesting that loss of MafB

<figure><img src="image_4.png" /><figcaption>**Fig. 1.** Whole mount imaging of intact embryonic pancreas. (A) Embryonic pancreas anlagen from E18.5 wild-type and *MafB*⁻⁻ pancreata were stained with insulin (green) and mucin (white) and imaged by OPT, surface rendering using Imaris is shown. To visualize larger ductal structures, the ellipsoid prolate filter was applied to the mucin surface tool. The residual ductal structures with major axes shorter than 400 µm were manually removed. Staining of smaller and larger ductal structures is shown in Fig. S1. (B) Distance between insulin⁺ and mucin⁺ surfaces was measured for all islet and ductal structures in wild-type and *MafB*⁻⁻ pancreata (n=4 embryos per genotype from three different litters). ****P<0.0001 (unpaired two-tailed t-test). Data are mean±s.e.m. (C) Immunohistochemical staining of wild-type and *MafB*⁻⁻ E18.5 pancreata with insulin (red), mucin (white), glucagon (green) and nuclei (blue) imaged using a Zeiss LSM780. (D) The average number of branching points in luminal structures of wild-type and *MafB*⁻⁻ (n=4 embryos per genotype) in the pancreatic and islet areas is shown. *P*=0.02857 for average branching points of ductal structures in islets (Wilcoxon test). Box plots show a central mark indicating the median and edges indicating 25th and 75th percentiles. Whiskers extend to the largest or smallest point contained within 1.5× of the interquartile range from both edges. (E) Pie charts showing size distribution (>10 µm²) of ductal structures in islet and pancreatic exocrine areas. Difference of size distribution of ductal structures between wild-type and *MafB*⁻⁻ pancreata were significant both in the pancreatic (Chi-square: 165.54; *P*<2.2e⁻¹⁶) and islet (Chi-square: 539.49; *P*<2.2e⁻¹⁶) areas. n=4 embryos per genotype. Scale bars: 700 µm (A); 50 µm (C).</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>C</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>B</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>D</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>Fig. 2. Gene expression profiles of MafB-GFP+ cells in embryonic pancreas. (A) Heatmap showing global expression profile of *MafB<sup>+/-</sup>* and *MafB<sup>-/-</sup>* E14.5 and E15.5 GFP+ cells, *n*=4 biological replicates per genotype: range is 0 (minimum) to 1 (maximum). The heat map was generated using Morpheus web browser from Broad Institute, wherein a relative color scheme uses the minimum and maximum values in each row to convert values to colors.</figcaption></figure>

(B) Heatmap showing single cell gene expression of E15.5 *MafB<sup>+/-</sup>* (90 cells) and *MafB<sup>-/-</sup>* (269 cells) GFP+ pancreatic cells. Cell clustering using expression profiles was performed using SCExV (Lang et al., 2015) web tool into various stages of α and β progenitor cells. Genes predominantly expressed in α or β progenitor cells are highlighted in blue or orange, respectively, while light blue indicates genes expressed in all endocrine progenitor cells. Range is 0 (minimum) to 1 (maximum). Gray indicates no expression. (C) Number of MafB+ cells in the different clusters of α and β progenitor cells. (D) Schematic overview of α and β progenitor cell stages and distribution of *MafB<sup>+/-</sup>* and *MafB<sup>-/-</sup>* cells in these groups using pie charts.

in α and β progenitor cells disrupts morphogenesis of islets and associated ductal structures only at later stages of pancreas development.

MafB is essential for the expression of several
neurotransmitter and axon guidance receptor genes in
developing β cells

MafB expression is detected in insulin- and glucagon-producing cells as well as a smaller fraction of Neurog3+ cells throughout pancreas development (Artner et al., 2006). Loss of MafB expression impairs differentiation of α and β cells in the developing pancreas, while the number of endocrine progenitor cells remains unchanged (Artner et al., 2007). MafB+/- and MafB-/- endocrine cells were isolated by fluorescence activated cell sorting (FACS) from E14.5 and E15.5 MafB GFP knock-in

embryos and global gene expression analysis (cDNA microarray) was performed to identify *MafB* target genes involved in α and β cell differentiation. Wild-type and *MafB<sup>+/-</sup>* pancreata were indistinguishable (Artner et al., 2007), thus *MafB<sup>+/-</sup>* was used to take advantage of the GFP gene knocked into the *MafB* locus (Blanchi et al., 2003). GFP was detected in both heterozygous and *MafB* mutant insulin- and glucagon-producing cells and the number of *GFP<sup>+</sup>* cells was comparable between genotypes (2% of all pancreatic cells, Fig. S2). Loss of *MafB* reduced the expression of several known MafB target genes such as *Ins1*, *Ins2*, *Slc30a8*, *Slc2a2*, *Pcsk2* and *G6pc2*, whereas expression of *Pax4*, which marks immature β cells, was increased (Fig. 2A; Table S1). As expected, gene ontology analysis showed that differentially expressed genes were associated with biological processes related to hormone secretion, but also with cell migration
<figure><img src="image_9.png" /><figcaption>**Fig. 3.** Neurotransmitter receptor expression in embryonic pancreas. (A-P) Immunohistochemical analysis of neurotransmitter receptor expression in wild-type and *MafB*<sup>-/-</sup> E18.5 pancreas. Co-staining for insulin (green) and glucagon (white) with neurotransmitter receptor proteins in red – *ChrnA3* (A,B), *ChrnA4* (C,D), *ChrnB2* (E,F), *ChrnB4* (G,H), *Adra2A* (I,J), *P2ry1* (K,L), *Robo1* (M,N) and *PlxnA3* (O,P) – visualized using a Zeiss LSM780. Nuclei are stained in blue. Arrows mark insulin producing cells co-expressing respective neurotransmitter receptor proteins, arrowheads denote PlxnA3<sup>+</sup> cells in (P), insets show co-expression of neurotransmitter receptor proteins with insulin and glucagon. Scale bar: 20 µm.</figcaption></figure>

and synaptic signaling (Table S2). Interestingly, expression of neurotransmitter receptor genes *ChrnA3*, *ChrnA4*, *ChrnB4*, *P2ry1*, *Adra2A* and the monoamine metabolizing enzyme *MaoB* was significantly reduced at E14.5 and E15.5. In addition, altered expression of axon guidance receptor genes *Robo1*, *Robo2*, *Nrp1*, *Nrp2*, *PlxnA3* and *PlxnA4* was also observed. These data suggest that *MafB* regulates transcription of axon guidance and neurotransmitter signaling pathways which are involved in cell migration and islet formation (Berger et al., 2015; Borden et al., 2013; Pauerstein et al., 2017).

To determine whether expression of genes associated with neurotransmitter signaling were altered in *MafB*<sup>-/-</sup> α and/or β

progenitor cells, Fluidigm single cell PCR analysis was performed. Expression analysis of E15.5 *MafB*<sup>+/−</sup> and *MafB*<sup>-/-</sup> GFP<sup>+</sup> cells showed that the majority of *MafB*<sup>+/−</sup> GFP<sup>+</sup> cells co-expressed β cell (*Ins1*, *Ins2*, red cell population Fig. 2B,C) or α cell [*Glu* (also known as *Gcg*), *Arx*, *Brn4* (*Pou3f4*), pink cell population] genes. Few cells with lower α or β cell marker expression were detected in the *MafB*<sup>+/−</sup> GFP<sup>+</sup> cell fraction (yellow, green and blue indicate less mature populations) and only one MafB<sup>+</sup> cell expressing both α and β cell marker genes (light blue cell population, Fig. 2B,C). In contrast, the majority of *MafB*<sup>-/-</sup> GFP<sup>+</sup> cells clustered with immature α and β cell progenitor cell populations (yellow, green and blue populations, Fig. 2C). Notably, a large number of cells
expressed both α and β cell marker genes (light blue cell population, Fig. 2B). Expression of neurotransmitter and axon guidance receptor (ChrnA3, ChrnB4, Robo1, Robo2, P2ry1) genes was predominantly observed in MafB⁺/⁻β cells, but reduced/absent in MafB⁻/⁻β progenitor cells, whereas none of our newly identified MafB target genes was specifically expressed in α cells. The majority of the MafB⁺/⁻ cells at E14.5 and E15.5 were β cell progenitors [roughly 75% (60/90 MafB⁺ cells), Fig. 2D] resulting in the identification of MafB target genes predominantly in β cells at these time points.

## Neurotransmitter and axon guidance gene expression is maintained in maturing α and β cells

Immunohistochemical analysis of E18.5 wild-type and MafB⁻/⁻pancreata was performed to determine whether neurotransmitter and axon guidance receptor expression is maintained throughout islet differentiation and dependent on the presence of MafB in maturing α and β cells. ChrnA3, ChrnA4, ChrnB4, P2ry1 and Robo1 were detected in developing wild-type α and β cells, as well as adjacent epithelial cells (Fig. 3) and expression was reduced in MafB⁻/⁻pancreata, with no (ChrnA3, P2ry1; Fig. 3B,L) and only few (ChrnA4, ChrnB4, Robo1; Fig. 3D,H,N) neurotransmitter receptor-expressing cells remaining. Adra2A expression was exclusively detected in the membrane of wild-type β cells, and was lost in MafB⁻/⁻pancreata (Fig. 3I,J). In contrast, ChrnB2 expression was maintained in MafB⁻/⁻cells (Fig. 3E,F) and PlxnA3, which is expressed in wild-type α and β cells, was detected in non-hormone producing mutant cells (Fig. 3O,P). The protein expression analysis largely confirmed our results obtained using single cell PCR (Fig. 2B) and demonstrates that neurotransmitter and axon guidance receptor expression is impacted upon MafB loss during late pancreas development.

## Neurotransmitter and axon guidance genes are expressed in fetal pancreata and associated with type 2 diabetes development

To evaluate whether neurotransmitter and axon guidance signaling genes are also expressed in the developing human pancreas, gene expression analysis of human pancreas anlagen was performed. Bulk RNA-seq of 7-14 weeks post conception (pc) pancreatic tissue (n=16) showed that neurotransmitter receptor and axon guidance genes are expressed in human pancreas anlagen (Fig. 4A-D). Single cell sequencing (scRNA-seq) and immunohistochemical analysis of human embryonic pancreas (Carnegie stage 22, 3199 cells) showed the presence of endocrine progenitor cells (MAFB⁺/PDX1⁺, Fig. 4F-H) and a distinct cluster of autonomic nerve cells (Fig. 4E,I,J; Fig. S3). Sympathetic and parasympathetic nerve cells were detected close to human and mouse MAFB⁺ cells (Fig. 4I,J; Fig. S4A). PDX1, MAFB, ADRA2A, CHRNB1 and ROBO1 were expressed in pancreatic progenitor cell clusters (Fig. 4F,G, Fig. S3) and co-expression analysis of MAFB and PDX1 with neurotransmitter and axon guidance receptor genes showed that a significant portion of MAFB⁺/PDX1⁺β progenitor cells expressed adrenergic, axon guidance and nicotinic receptor genes (Fig. 4K; Fig. S3). Several MAFB⁺/PDX1⁺cells expressing multiple neurotransmitter/axon guidance receptor genes (Fig. 4K) were detected. Expression of CHRNB1, ADRA2A and ROBO1 in the pancreatic epithelium, MAFB- and insulin-producing cells was also detected using immunohistochemical staining of human embryonic pancreas (Fig. 4L-N; Fig. S4B).

To evaluate whether loss of MAFB also affects expression of adrenergic, nicotinic and axon guidance receptor genes, scRNA-seq

of *in vitro* differentiated β cells derived from MAFB-deficient human embryonic stem cells and MAFB ChIP-seq from ENDOCβH2 β cells (Russell et al., 2020) were re-analyzed. Expression of *CHRNA3*, *CHRNB2*, *ROBO2*, *NRP1*, *NRP2* and *ADRA2A* was reduced in MAFB-deficient β-like cells (Table 1). In addition, MAFB binding to *CHRNB2*, *ROBO1*, *ROBO2* and *NRP1* genomic regions was detected in the human embryonic ENDOCβH2 β cells (Table 1), suggesting that MAFB directly binds to and activates these genes during human β cell development. Interestingly, expression levels of several of these genes were also correlated with MAFB expression levels in adult islets (Table 1), MAFB binding to neurotransmitter (P2RY1) and axon guidance receptor (ROBO1, ROBO2, NRP1, NRP2) genes was also detected in adult human islets (Table 1) and quantitative trait loci (eQTL) influencing islet gene expression (Fig. S5; Table 1; PLXNA4, ROBO2, P2RY1) were identified. We also identified P2RY1, ROBO2 and PLXNA4 variants which were nominally associated with gene expression levels of their respective genes and also showed suggestive signals of association with insulin secretion indices (Fig. S4; Table S3; Prokopenko et al., 2014; Wood et al., 2017). These data indicate that transcriptional regulation of several neurotransmitter and axon guidance receptor genes depends on MAFB expression, is conserved between human and mouse, and is relevant to adult islet function.

## nAChR signaling is required for β cell-nerve cell interactions

Our analysis has shown that loss of MAFB impairs expression of several neurotransmitter (ChrnA3, ChrnA4, ChrnB4, Adra2A, P2ry1) and axon guidance receptor (Robo1, Robo2) genes in mouse and human β cells. To evaluate whether loss of MafB and concomitant changes in gene expression altered islet innervation, immunohistochemical stainings of insulin and glucagon, with the neuronal marker β III tubulin and the glial marker Sox10 were performed at E18.5 in mouse wild-type and MafB⁻/⁻pancreata. β III tubulin⁺ nerve fibers and Sox10⁺ glial cells were readily detected in the developing pancreas and surrounding tissues (Fig. 5A). The β III tubulin⁺ cell area within the pancreas and pancreatic ganglia and Sox10⁺ area in pancreatic ganglia were not altered between wild-type and MafB⁻/⁻pancreata indicating that loss of MafB is not sufficient to impair islet innervation of the mutant pancreata (Fig. 5B; Fig. S4C), possibly due to continued expression of NGF and GDNF (Table S1) which are crucial for islet innervation (Glebova and Ginty, 2005).

Previous studies have demonstrated that adrenergic signaling is required for β cell migration towards developing sympathetic nerves (Borden et al., 2013). To evaluate whether loss of nAchR activity also impaired β cell migration towards autonomic nerves, we examined how migration was impacted upon co-culturing developing β cells from mouse insulin promoter GFP transgenic animals (MIP-GFP⁺; Hara et al., 2003) with superior cervical ganglia (acetylcholine source, Fig. S4D) in the presence of nicotinic receptor agonists (nicotine) or antagonists (mecamylamine). Untreated MIP-GFP⁺β cells migrated towards the centrally located ganglion (Fig. 5C-E) and 5% of MIP-GFP⁺ cells were located within the ganglionic structure after 6 days of co-culture. In contrast, blocking nicotinic signaling with mecamylamine resulted in a loss of directed migratory behavior with no accumulation of MIP-GFP⁺ cells in the ganglionic structure (Fig. 5C-E). Treatment of co-cultures with nicotine, which activates nAchRs, also did not result in a significant accumulation of MIP-GFP⁺cells in ganglionic tissue, most likely due to loss of directionality (no gradient formed) and/or receptor desensitization. Moreover, co-culture of ganglionic tissue and MIP-GFP⁺cells promoted accumulation of β cells in
<figure><img src="image_10.png" /><figcaption>**Fig. 4.** Neurotransmitter and axon guidance signaling in human embryonic pancreas. (A-D) RNA-seq of human embryonic pancreas (7-14 weeks gestation) of nicotinic (A), Robo/Slit (B), semaphorin (C) and purinergic (D) signaling (n=16). Genes that were differentially expressed in mouse *MafB*⁻/⁻GFP⁺ cells are indicated in orange, genes with impaired expression in *MafB*⁻/⁻GFP⁺ cells and reported type 2 diabetes risk alleles are indicated in green. Data were expressed as log counts per million mapped reads (logCPM), n=16. (E) t-SNE visualization of single cell sequencing data from human embryonic pancreas (Carnegie stage 22, n=3199). Clusters were identified using top 20 expressed genes per cluster (Table S3). Neural cells are labeled in gray. (F,G) t-SNE showing expression of MAFB (F) and PDX1 (G). (H-J) Co-staining of PDX1 (red) and MAFB (green), arrows depict cells co-expressing PDX1 and MAFB (H), β III-tubulin (I), and vAchT and TH (J) with MAFB in human embryonic pancreas (Carnegie stage 22). Arrows depict nerve fibers, arrowheads show MAFB⁺ nuclei (I,J). (K) Co-expression analysis (circular plot) of MAFB⁺ and MAFB⁺/PDX1⁺ cells with ADRA2A, CHRNA3, CHRNA4, CHRNB1, ROBO1, ROBO2 and P2RY1 in human fetal scRNA-seq, pie-charts illustrating co-expression of ADRA2A, CHRNB1, and ROBO1 with MAFB and PDX1. (L-N) Co-staining of ADRA2A (L, green), CHRNB1 (M, green), and ROBO1 (N, green) with insulin (white) and PDX1 (red). Nuclei are in blue. Insets show expression of insulin (red), PDX1 (blue) and the respective receptor (green) in single channels. Scale bars: 10 µm (I,J); 20 µm (H,L-N).</figcaption></figure>
<table><caption>Table 1. Expression and regulation of neurotransmitter and axon guidance genes in embryonic and adult β/islet cells</caption><thead><tr><th>Gene</th><th>Expression changes in MafB-deficient differentiated β-like cells, reanalyzed from Russell et al. (2020) (fold change)</th><th>MafB binding in regulatory regions in EndoCβ h2, re-analyzed from Russell et al. (2020)</th><th>MafB binding in regulatory regions in adult islets (Mularoni et al., 2017)</th><th>Correlation with MafB expression in adult islets</th><th>Enriched in FACS-sorted fetal β cells (Blodgett et al., 2015)</th></tr></thead><tbody><tr><td><strong>Nicotinic signaling</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>CHRNA3</td><td>-2.0, P=2.82e-168</td><td>—</td><td>—</td><td>r=0.2, P=0.005</td><td>—</td></tr><tr><td>CHRNA4</td><td>—</td><td>No signal</td><td>—</td><td>—</td><td>—</td></tr><tr><td>CHRNB1</td><td>—</td><td>—</td><td>—</td><td>r=-0.148, P=0.05</td><td>+</td></tr><tr><td>CHRNB2</td><td>-1.25, P=1.49e-39</td><td>+</td><td>—</td><td>r=0.69, P=1.35e-28</td><td>+</td></tr><tr><td>CHRNB4</td><td>—</td><td>—</td><td>—</td><td>r=0.32, P=0.0009</td><td>+</td></tr><tr><td><strong>Axon guidance</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>ROBO1</td><td>—</td><td>+</td><td>+</td><td>r=0.59, P=4.39e-19</td><td>+</td></tr><tr><td>ROBO2</td><td>-1.8, P=6.64e-220</td><td>+</td><td>+</td><td>r=0.65, P=6.35e-24</td><td>+</td></tr><tr><td>SLIT1</td><td>—</td><td>—</td><td>+</td><td>r=0.71, P=4.9e-31</td><td>—</td></tr><tr><td>SLIT2</td><td>—</td><td>+</td><td>—</td><td>r=-0.11, P=0.11</td><td>—</td></tr><tr><td>PLXNA3</td><td>—</td><td>—</td><td>—</td><td>r=0.09, P=0.2</td><td>+</td></tr><tr><td>PLXNA4</td><td>—</td><td>—</td><td>—</td><td>r=0.08, P=0.24</td><td>—</td></tr><tr><td>NRP1</td><td>-1.25, P=7.80e-37</td><td>+</td><td>+</td><td>r=-0.14, P=0.04</td><td>—</td></tr><tr><td>NRP2</td><td>-1.18, P=8.65e-27</td><td>—</td><td>+</td><td>r=-0.11, P=0.13</td><td>—</td></tr><tr><td><strong>Purinergic signaling</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>P2RY1</td><td>—</td><td>—</td><td>+</td><td>r=0.59, P=1.84e-19</td><td>+</td></tr><tr><td><strong>Adrenergic signaling</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>ADRA2A</td><td>-3.29, P=1.28e-136</td><td>—</td><td>—</td><td>r=0.19, P=0.007</td><td>+</td></tr><tr><td>MAOA</td><td>-1.3, P=6.65e-17</td><td>+</td><td>—</td><td>r=0.04, P=0.5</td><td>—</td></tr><tr><td>MAOB</td><td>—</td><td>—</td><td>—</td><td>r=-0.04, P=0.55</td><td>—</td></tr></tbody></table>

+, positive signal; –, no signal; bold indicates significant changes in gene expression and/or correlation with MAFB expression.

clusters (Fig. 5C,F). MIP-GFP+ cell clustering was impaired upon treatment with mecamylamine (Fig. 5F), whereas nicotine treatment affected clustering to a lesser extent. These data suggest that nAchR activity (which is impaired in MafB<sup>-/-</sup> β cells) is required for β cell clustering and migration towards autonomic tissue.

## DISCUSSION

Maf transcription factors contribute to β cell development and adult β cell function. Previous studies have shown that loss of MafB impairs α and β cell terminal differentiation during late embryonic pancreas development (Artner et al., 2007; Nishimura et al., 2008). Here, we examined expression profiles of MafB-deficient endocrine progenitor cells at E14.5 and E15.5 to identify MafB target genes controlling α and β cell commitment during the secondary transition of β cell development. Analysis of MafB mutant embryonic pancreata demonstrate that MafB promotes the rapid progression of α and β cell progenitors to mature endocrine cells and subsequently islet formation. MafB mutant β cells had altered expression of axon guidance (Robo1, Robo2, Nrp1, Nrp2) and neurotransmitter receptor (ChrnA3, ChrnA4, ChrnB4, P2ry1, Adra2A) genes, suggesting that MafB controls responsiveness to environmental cues, which may contribute to endocrine cell differentiation, migration and islet formation during pancreas development.

Maf transcription factors MafA and MafB are crucial for adult β cell function by controlling transcription of glucose sensing, insulin production and secretion genes. MafB activates Pdx1 and MafA transcription in developing β cells, thereby promoting terminal differentiation of β cells (Artner et al., 2007). MafB expression in Neurog3<sup>+</sup> progenitor cells suggests additional functions in α and β cell commitment. Our single cell gene expression analysis of endocrine progenitor cells showed that most MafB<sup>+</sup> endocrine progenitor cells expressed mature α or β cell marker genes. In contrast, the vast majority of MafB<sup>-/-</sup> cells had a more immature expression profile and a significant proportion of cells expressed

both α and β cell marker genes (e.g. Pdx1, Pax4 and Arx), suggesting that lack of MafB inhibits, or at least delays, specification of α and β cells. This notion was further supported by previous studies illustrating that selective deletion of MafB using Pdx1-Cre resulted in impaired β cell differentiation that was eventually compensated for by adult MafA expression (Conrad et al., 2016). Non-committed α/β cell progenitor cells were also detected in the MafB<sup>+</sup> endocrine progenitor population, albeit at very low numbers, indicating that progenitor cells quickly transition through this stage upon initiation of MafB expression (Fig. 2D). A similar accumulation of immature endocrine progenitor cells was also observed in β-like cells derived from MAFB-deficient human embryonic stem cells (ES cells) (Russell et al., 2020), but loss of MAFB in this *in vitro* setting also resulted in an enrichment of somatostatin-producing delta cells, which we did not observed in mouse pancreas anlagen. This may reflect species differences in MafB function or may be due to differences in cell culture versus *in situ* development of the pancreas.

We re-evaluated whether loss of MafB impaired islet formation and pancreas morphogenesis. Interestingly, whole mount image analysis showed that E18.5 MafB<sup>-/-</sup> pancreata had defects in the formation of larger ductal structures and that the remaining β cells were closely associated with smaller ducts. This was somewhat surprising considering that E12.5 MafB<sup>-/-</sup> pancreas anlagen were indistinguishable from wild-type littermates and that MafB expression was only detected in the endocrine lineage from E10.5 (Artner et al., 2006; Nishimura et al., 2006), with no expression in ductal cells been reported so far. During pancreas development, endocrine progenitor cells exit the pancreatic epithelium to form islets of Langerhans and transcription factors important for these processes have been identified (Bankaitis et al., 2018). For example, epithelial malformations including impaired ductal branching were also observed in tissue specific knockout models of Hnf1b (De Vas et al., 2015), Hnf6 (Onecut1; Pierreux et al., 2006), Neurog3 (Magenheim et al., 2011) and Prox1 (Westmoreland et al., 2012). Interestingly, the defects in ductal branching were only observed in
<figure><img src="image_11.png" /><figcaption>**Fig. 5.** Nicotinic receptor function in embryonic islet development. (A) Immunohistochemical staining of wild-type and *MafB*⁻/⁻ E18.5 pancreas for insulin (green), β III-tubulin (red) and glucagon (white). Nuclei are in blue. (B) β III-tubulin⁺ area in E18.5 wild-type and *MafB*⁻/⁻ pancreas. *n*=4 per genotype. Error bars indicate s.e.m. (C) Tile scan of co-cultures of MIP-GFP⁺ (green) β cells and superior cervical ganglia (black, placed in the center) treated with mecamylamine, nicotine, and control for 6 days imaged with a Zeiss LSM780 confocal microscope, objective 5×. Insets show GFP expression. (D) Visualization of cell tracking of MIP-GFP⁺ β cells during 6 days of co-culture using Imaris software (version 9.6.0). (E) Quantification of MIP-GFP⁺ cells located in ganglion, *n*=3-4 per condition. (F) Mean MIP-GFP⁺ cluster size after 6 days of co-culture, *n*=3-4, error bars indicate s.d. Significance of changes between day 1 and 6 for each condition was analyzed using two-way Anova with Sidak's multiple comparison test (E and F). Scale bars: 50 µm (A); 700 µm (C,D).</figcaption></figure>
the vicinity of mutant islets, suggesting that the morphogenetic
changes observed in *MafB*⁻/⁻ ductal organization may arise from
defects in α and β cell migration/exit from the epithelium.

Endocrine cell migration from the pancreatic epithelium and
subsequent islet formation depend on external cues from the
mesenchyme, blood vessels and autonomic nerves. Gene expression
profiling of *MafB*⁻/⁻ endocrine progenitor cells showed that
expression of several axon guidance and neurotransmitter receptor
genes was altered in *MafB*⁻/⁻ β progenitor cells. Previous studies by
Pauerstein et al. have implicated semaphorin signaling from the
mesenchyme in the initiation of α and β cell migration/islet
formation (Pauerstein et al., 2017) and show that Nrp/Plxn receptors
are expressed in α and β progenitor cells. Consistent with this
finding, we detected an upregulation of Nrp and Plxn receptors in
*MafB*⁻/⁻ endocrine progenitor cells. In contrast, *Robo1* and *Robo2*
axon guidance receptor expression was reduced in *MafB* mutant
endocrine cells, which may contribute to the defects in islet
formation observed in the knockout as Robo/Slit signaling
participates in islet cell sorting and morphogenesis (Adams et al.,
2018).

MafB mutant β cells had drastically reduced expression of
adrenergic, nicotinic and purinergic receptor genes (*Adra2a*,
*ChrnA3*, *ChrnA4*, *ChrnB4*, *P2ry1*). Acetylcholine and epinephrine
secretion has been reported from autonomic nerves and blood
circulation, but also α cells (acetylcholine; Rodriguez-Diaz et al.,
2011b) in the islet, where epinephrine and acetylcholine regulate
insulin secretion in the adult islet via nicotinic, muscarinic and
adrenergic receptors expressed in β cells (reviewed by Rorsman and
Ashcroft, 2018). During development, adrenergic signaling from
pancreatic nerves is crucial for islet formation, as blocking of β
adrenergic receptors reduces β cell numbers and impairs migration of
β cells to autonomic nerves (Borden et al., 2013). We show that β cell
migration and clustering is impaired upon inhibition of nicotinic
signaling in co-cultures of embryonic β cells with autonomic nerves,
suggesting that acetylcholine signaling contributes to islet innervation
and formation. A similar effect was observed upon addition of
nicotine to co-cultures, indicating that a gradual activation of receptors
is required for providing a directional cue that results in migration of β
cells to pancreatic nerves while uniform exposure does not allow for
directed migration. Previous studies have shown that nicotinic
receptors participate in controlling migratory behavior of bronchial
epithelial cells (Tournier et al., 2006) by modulating intracellular Ca²⁺
levels, suggesting that a similar mechanism may be in place in
developing β cells.

Differences in adult islet innervation suggest that
parasympathetic acetylcholine signaling contributes to adult β cell
function in rodents, but not in humans (Rodriguez-Diaz et al.,
2011a). However, immunohistochemical studies in fetal human
pancreas show interactions between nerve and islet cells (Krivova
et al., 2016) and we have detected a distinct nerve cell population in
our scRNA-seq (Fig. 4E, gray) and immunohistochemical (Fig. 4I,J;
Fig. S4B) analysis of human embryonic pancreas, suggesting that
autonomic nerves may also play a role in human islet morphogenesis.
Transcripts and protein of neurotransmitter receptor genes were
readily detectable in human embryonic pancreas anlagen and
developing β cells further supporting a role of these signaling
pathways in human pancreas development. The relationship between
*MAFB* and neurotransmitter receptor expression is conserved between
differentiating human and mouse β cells as human *MAFB⁺/PDXI⁺*
cells co-express *ADRA2A*, *CHRNB1* and *ROBO1*, and nicotinic,
adrenergic and axon guidance receptor gene expression was impaired
in MAFB-deficient human β-like cells. Moreover, correlation of

MAFB expression with expression of several neurotransmitter
receptor and axon guidance genes in adult human islets indicates
that these signaling pathways have a dual role during islet
development and in adult islet function. This notion is further
supported by the presence of type 2 diabetes susceptibility signals
intersecting eQTLs for the *ADRA2A* (Rosengren et al., 2010),
*CHRNB4* (Ganic et al., 2016) and P2RY1 genes.

Here, we describe a novel function of the MafB transcription
factor in islet morphogenesis by controlling neurotransmitter
sensing in developing β cells. Specifically, we identify nicotinic
signaling as a key mechanism in islet-nerve cell communication and
β cell migration and clustering. These findings describe a novel
mechanism by which prenatal nicotine exposure, which has been
associated with an increased risk of type 2 diabetes, development
and impaired endocrine cell development (Somm et al., 2008), may
affect islet formation and subsequent β cell function.

MATERIALS AND METHODS

Animals

MafB-GFP knock-in (*MafB*⁻/⁻) (Blanchi et al., 2003) and MIP-GFP (Hara
et al., 2003) mouse lines have been previously described. Mice were
maintained on a C57BL/6 background. Experimental procedures were
performed using embryonic pancreata (E14.5, E15.5 and E18.5). All animal
procedures were pre-approved and performed according to Swedish national
guidelines (approval number 5733/2018).

Immunohistochemistry/whole mount staining

For immunostaining, pancreata from E18.5 wild-type and *MafB*⁻/⁻ embryos and human embryonic pancreas (Carnegie stage 22) were dissected, fixed in 4% paraformaldehyde (Sigma-Aldrich) and embedded in paraffin (Thermo Fisher Scientific) before being sectioned (6 µm). Immunohistochemical analysis was performed as previously described (Matsuoka et al., 2003). The primary antibodies used were: guinea pig α-insulin (1:2000, Dako, A0546), mouse α-glucagon (1:2000, Sigma-Aldrich, G2654), rabbit α-CHRN-A3/A4/B2/B4 (ANC-003, ANC-004, ANC-012, ANC-014, respectively; 1:100; Alomone Labs), rabbit α-CHRNB1 (1:100, Abcam, AB 126234), rabbit α-Adra2A, (1:100, Sigma-Aldrich, A271), goat α-Adra2A (1:100 Sigma-Aldrich, SAB2500033), mouse α-Robo1 (v-Robo1-s, 1:100, DSHB, 13C9), rabbit α-PlxnA3 (1:100, Alomone Labs, APR-093), rabbit α-P2RY1 (1:100, Alomone Labs, APR-009), rabbit α-GFP (1:200, Molecular Probes, A11122), Armenian hamster α-mucin (1:500, Thermo Fisher Scientific, 1630-R7), goat α-PDX1 (1:1000, Abcam, ab47383), rabbit α-MAFB (1:250, Atlas Antibodies, HPA005653), mouse α-MAFB (1:50, Bio-Techne, NBP2-45718), goat α-Sox10 (1:200, Santa Cruz, Sc-17343), rabbit α-β-tubulin (1:5000; Covance, PRB-435P), mouse α-vAchT (1:100, Thermo Fisher Scientific, MA5-27662), sheep α-Tyrosine Hydroxylase (1:200, Thermo Fisher Scientific, PA1-4679). Secondary antibodies used were Cy2-, Cy3- and Cy5-conjugated α-guinea pig, α-mouse, α-Armenian hamster, α-goat and α-rabbit (1:500; Jackson ImmunoResearch, 706-225-148, 706-165-148, 706-175-148, 715-225-151, 715-165-151, 715-175-150, 127-225-160, 705-225-147, 705-165-147, 705-175-147, 711-225-152, 711-165-152 and 711-175-152). Nuclear counterstaining was performed using DAPI (1:6000; Thermo Fisher Scientific). Neurotransmitter/axon guidance receptor antibodies were previously validated in the following publications: Adra2A, ChrnA4, ChrnB2, ChrnB4 (Ganic et al., 2016), P2ry1 (Caputi et al., 2017), ChrnB1 (proteinatlas.org) and Robo1 (Glawe et al., 2013).

Histological analysis of embryonic pancreas

E18.5 wild-type and *MafB*⁻/⁻ pancreata (*n*=4 in each group) were paraffin
embedded and sectioned (see above). For histological analysis, sections
with a distance of 70 µm were stained for Sox10, β III tubulin or mucin with
insulin, glucagon and DAPI. Stained sections were scanned on an Olympus
slide scanner (VS120) with a 20× objective. Pictures were analyzed using
Pathology AI (Visiopharm). The application was first trained to recognize
islet, pancreatic, gut and ganglionic structures from wild-type and *MafB*⁻/⁻
sections by manual outlining of the respective structures, the application was
then allowed to train for 143,000 iterations using Deep Labv3+ classification. Accuracy of tissue recognition was assessed afterwards. Insulin-, glucagon-, β III tubulin-, mucin- and Sox10-stained area was then determined using a fluorescence threshold, and respective areas were then normalized against pancreatic, islet or ganglionic areas. Mucin structures were categorized into stained areas larger or smaller than 10 µm² and the size composition of mucin lumina per islet and pancreatic area was assessed. Branching points were automatically recognized and the average number of branching points per mucin structure in islet versus pancreatic area determined.

## Optical projection tomography (OPT)

OPT scans were performed using the Bioptomics 3001 OPT M scanner (Bioptomics) with exciter D560/40× and emitter E610lpv2 filter (Chroma). Tomographic re-constructions were generated using the NRecon V1.6.1.0 (Skyscan) software, and reconstructed images were further assessed using the Bioptomics viewer V2.0 as previously described (Alanentalo et al., 2010). Final images were made by using Imaris (version 9.6.0, Oxford Instruments).

## Imaris image analysis

Fluorescence images datasets obtained from OPT scans of individual pancreata were converted and analyzed using the Imaris software (version 9.6.0 and 9.9.1, Oxford Instruments). The stomach had to be removed from each dataset due to its strong positivity for mucin. In order to do this, the contour tool was used and an ad hoc surface object for the stomach was created manually. Subsequently the newly created surface was used as a Channel Mask and the voxels inside the surface were set to zero in order for the surface to be deleted. Insulin-positive cells and mucin-positive ducts were segmented using the 'Surface' object tool. To obtain the distance measurement between β cells and ducts, 'Distance Transformation' within the XTension tool was used. The distance was calculated outside of the β cells and any ductal lumen surface. The 'Surface' object tool was used to get measurements of the ductal volume.

Time-lapse movies of *in vitro* islet migration assays obtained using a laser scanning confocal microscope (Zeiss 780) were analyzed with the Imaris software (version 9.6.0). GFP⁺β cells were identified and tracked over time using the 'Spots' function. Individual cells are visualized as green spheres (Spot Objects) with a radius set to 5 µm, and tracks are displayed as time-color coded Full Tracks. The maximum distance that a tracked Spot could move between two consecutive time points was set to 180 µm.

## Isolation of GFP⁺ murine pancreatic cells for global gene expression analysis

E14.5 and E15.5 *MafB* GFP knock-in pancreata (*MafB*⁺/⁻ and *MafB*⁻/⁺) were dissected in cold PBS and dissociated with Accumax solution (Sigma-Aldrich) containing 5 µg/ml DNase I (Sigma-Aldrich) at 37°C for 1 h with periodic resuspension. Quenching solution [L15 medium, 10 mM Hepes (pH 7.0), 1 mg/ml bovine serum albumin (Thermo Fisher Scientific), 5 µg/ml DNase I] was added to stop the reaction. The cell suspension was filtered through a 30 µm filter into FACS tubes, washed with quenching solution and resuspended in 0.5 ml quenching solution containing 7AAD viability dye (1:2000; Thermo Fisher Scientific). Single GFP⁺ cells were sorted directly into Trizol LS reagent (Thermo Fisher Scientific) using a FACS Aria II cell sorter (BD Biosciences). RNA extracted from 2000-5000 GFP⁺ cells was labeled, hybridized to the Mouse Gene 2.0/2.1 Array (Affimetrix) and scanned according to the manufacturer's instructions.

Scanned image files were first inspected for quality control (QC) using built-in QC tools from the Bioconductor package of R. QC consisted of visual examination of probe array images, scatter plots from replicates, hierarchical clustering of array hybridizations, RNA degradation plots and MVA plots. Detection calls indicating the presence or absence of signal from each probe set were obtained by processing the raw data with the Microarray Analysis Suite 5.0 (MAS5). Feature intensity values from scanned arrays were normalized and reduced to expression summaries using the Robust Multiarray Algorithm (RMA) implemented in the affy package in R (Gautier et al., 2004) and log2 counts were obtained.

Probesets 3' locations were obtained by downloading the respective probe tab files made available by the Affymetrix online support. A probeset location was considered equal to the 3' distance of the probe that is most distal from the 3' end of the corresponding target within the set.

To test for differential expression, the Bayesian adjusted *t*-statistics from the linear models for Micoarray data (limma) package was implemented (Ritchie et al., 2015). A multiple testing correction based on the false discovery rate (FDR) was performed and FDR corrected *P*-values of <0.05 were considered to be significant.

## Gene expression analysis using Biomark Fluidigm dynamic arrays

E15.5 single GFP⁺ cells from *MafB*⁺/⁻ and *MafB*⁻/⁺ pancreata were sorted using FACS ARIA II cell sorter (BD Biosciences) directly into a 96-well plate containing lysis buffer. A detailed procedure using Biomark 96.96 dynamic array platform (Fluidigm) and gene expression analysis has been previously described (Singh et al., 2019). Primer sequences are presented in Table S4 and Larsen et al. (2017). Data were further analyzed and heatmaps/ violin plots were generated with SCExV using hierarchical clustering (Lang et al., 2015).

## Human tissue

### Human fetal pancreas

Human fetal pancreas tissue was obtained from material available following elective termination of pregnancy at the University Hospital in Malmö, Sweden. The study was approved by the national ethics committee (permit number 2018-579), written and oral consent was obtained before collection. Fetal/embryonic age was determined by measuring crown-rump length and Carnegie stage, where appropriate. Pancreata were either dissected and homogenized in Trizol and further extracted using RNeasy Mini columns (Qiagen) for bulk RNA extraction (Mazur et al., 2013) or single cell suspensions were prepared using a combination of dissociation enzymes accutase (Sigma-Aldrich) and liberase (Roche) (in 10:1 dilution, respectively). Fetal pancreas/enzyme mix was incubated at 37°C for 30 min with continuous shaking and periodic resuspension at 10-min intervals. Single cells were immediately spun at 1100 rpm (100 g) at 4°C for 5 min and resuspended in cold DPBS Ca²⁺/Mg²⁺ free buffer (Thermo Fisher Scientific) with 5% fetal bovine serum (Sigma-Aldrich), 0.5 mM EDTA (Sigma-Aldrich) and 5 µg/ml DNase I (Sigma-Aldrich). After washing, the cell suspension was filtered through a 30 µm filter (BD Biosciences), viability was assessed, and 8000 cells were resuspended and loaded into a 10x Genomics Chromium Single Cell System.

## Single cell sequencing and analysis

Libraries were generated using version 3 chemistry according to the manufacturer's instructions (10x Genomics). Libraries were sequenced on Illumina NextSeq500 (400 million reads flow cells) using the recommended read length with 42 bp paired-end reads. Sequencing data were processed through the Cell Ranger pipeline (10x Genomics, Cell Ranger version 6.0.1) to perform demultiplexing using the 'mkfastq' function. Resulting sequences were aligned to GrCH38 (v3.1.0), filtered using default parameters, and barcode and UMI counting was performed using the CellRanger 'count' function. In total, 3199 cells were obtained and included in the analysis. Unsupervised graph-based clustering was performed using Loupe browser and validated using custom R scripts. Graph-based clustering method is facilitated by building a graph where each node is connected to its nearest neighbors in high-dimensional space, with edges being weighted based on similarity between closely related cells. Assumption-free clustering of cells is then performed based on this similarity. Clusters were assigned using top 20 expressed genes (Table S3).

Single cell sequencing data from human *MAFB*⁺/⁺ and *MAFB*⁻/⁻ β-like cells (derived from iPSCs) were re-analyzed and expression of neurotransmitter receptor genes determined (Russell et al., 2020). The matrix of gene counts versus cells from this study was re-analyzed. The R package Seurat (version 2.3.4) was used to process the unique molecular identifier count matrix (Satija et al., 2015) and to perform data normalization (gene expression measurements for each cell were normalized by total
expression and log-transformed), dimensionality reduction and clustering, and
the Wilcoxon rank sum test was used for differential expression (DE) analysis.

Bulk RNA-seq and analysis

Sequencing data were demultiplexed using bcl2fastq. The resulting
sequences were aligned against the human genome build 38 using STAR
v2.4.1 and gene features were counted using 'featurecounts v 1.4.4' (Liao et al., 2014), with GENCODE version 22. Raw counts were normalized
using trimmed mean of M-values and log2-transformed. Spearman
correlations (R) were calculated to assess correlation of MAFA and MAFB
expression with genes expressed in the human pancreatic islets using R
language programming. eQTL analysis was carried out as previously
described (Asplund et al., 2020; Fadista et al., 2014) using islets from the
same 191 organ donors using the R Matrix eQTL packages (Shabalin,
2012). eQTLs for genes of interest were extracted and looked-up further for
association with type 2 diabetes and insulin secretion in previous GWAS
datasets, from DIAGRAM, DIAMANTE and MAGIC datasets.

**Statistical analyses**

All graphs and statistical analyses were carried out using GraphPad Prism
software. Unless otherwise noted in the figure legends, statistical
significance was determined using an unpaired two-tailed Student's t-test.
Details regarding replicates are indicated in the respective figure legends.

In vitro islet migration assays

Co-cultures were performed as previously described (Borden et al., 2013) with the following modifications: Cultures were maintained for 6 days at 37°C in 5% CO₂, with 10% fresh medium/H₂O solution (1:1) added daily. In experiments using nicotinic antagonists, after plating, cultures recovered overnight in culture media containing NGF alone. The next day, fresh medium and mecamylamine (Sigma-Aldrich) or nicotine (Sigma-Aldrich) were added. NGF was added every other day. Cultures were plated in duplicates for each condition.

Cultures were imaged every 10 h, with time point one being the day after plating, using a Zeiss LSM780 microscope and a 10× objective. Z-stacks were taken with 20 µm between images and 25 tiles to cover the whole area, images were stitched and maximum projection generated for subsequent analysis (Zen software version 2012). ImageJ software and the Sholl analysis plug-in (Ghosh Lab) was used to determine cell cluster size at day 1 and 6. To perform the analysis, a complete image of the culture was subjected to a threshold to only include the $GFP^+$ cells using the same settings for each image. To analyze the percentage of $\beta$ cells that migrated towards ganglionic cells, $GFP^+$ cells that overlapped with ganglionic structures (brightfield) at the beginning and end of the time lapse were analyzed.

**Acknowledgements**

The authors thank R. Stein (Vanderbilt University Medical Center, Nashville) for
discussing the data; Anne Grapin Botton (Max Planck Institute, Dresden) for sharing
primer sequences for Fluidigm single cell PCR ahead of publication and
M. Magnusson for mouse strain maintenance. Flow cytometry experiments and
image acquisition were supported by resources from the Lund Stem Cell Center
FACS and imaging cores.

**Competing interests**

The authors declare no competing or financial interests.

**Author contributions**

Conceptualization: E.M., T.S., I.A.; Methodology: S.B., E.M., T.S., J.K.J., A.S.-C., D.H.; Validation: S.B., E.M., T.S., J.K.J., K.A., L.B.-B., I.A.; Formal analysis: E.M., T.S., J.K.J., K.A., R.B.P.; Investigation: I.A.; Resources: Z.K., I.A.; Data curation: J.K.J., R.B.P.; Writing - original draft: I.A.; Writing - review & editing: S.B., E.M., T.S., A.S.-C., D.H., Z.K., R.B.P.; Visualization: E.M., J.K.J., A.S.-C., R.B.P., I.A.; Supervision: I.A.; Project administration: S.B., J.K.J.

**Funding**

This work was supported by the Vetenskapsrådet (2020-01461), Diabetesfonden,
Novo Nordisk Foundation Center for Basic Metabolic Research
(NNF200C0063485), Barndiabetesfonden, Direktör Albert Påhlssons Stiftelse,
Insamlingsstiftelsen Diabetes Wellness Network Sverige and the European
Foundation for the Study of Diabetes. Our laboratory, as part of the Lund University

Diabetes Centre-Industrial Research Centre, was supported by the Stiftelsen för
Strategisk Forskning (Dnr IRC15-0067) and Strategic Research Areas Exodiab and
StemTherapy (Vetenskapsrådet). Open access funding provided by Lund
University. Deposited in PMC for immediate release.

Data availability

All sequencing data related to fetal samples are deposited at the LUDC repository (www.ludc.lu.se/resources/repository) under the following accession numbers and are available upon reasonable request (due to GDPR these data cannot be made public as maternal donors could be identified from donor tissue): LUDC2021.10.6 (gene expression data from microarrays on E14.5 and E15.5 MafB GFP knock-in pancreata - MafB<sup>+/-</sup> and MafB<sup>-/-</sup>), LUDC2021.10.12 (bulk RNA-seq data from fetal pancreata) and LUDC2021.10.18 (scRNA-seq data from fetal pancreas).

Peer review history

The peer review history is available online at https://journals.biologists.com/dev/article-lookup/doi/10.1242/dev.201009

References

Adams, M. T., Gilbert, J. M., Hinojosa Paiz, J., Bowman, F. M. and Blum, B. (2018). Endocrine cell type sorting and mature architecture in the islets of Langerhans require expression of Roundabout receptors in beta cells. *Sci. Rep.* **8**, 10876. doi:10.1038/s41598-018-29118-x

Alanentalo, T., Hornblad, A., Mayans, S., Karin Nilsson, A., Sharpe, J., Larefalk, A., Ahlgren, U. and Holmberg, D. (2010). Quantification and three-dimensional imaging of the insulitis-induced destruction of beta-cells in murine type 1 diabetes. *Diabetes* **59**, 1756-1764. doi:10.2337/db09-1400

Artner, I., Le Lay, J., Hang, Y., Elghazi, L., Schisler, J. C., Henderson, E., Sosa-Pineda, B. and Stein, R. (2006). MafB: an activator of the glucagon gene expressed in developing islet alpha- and beta-cells. *Diabetes* **55**, 297-304. doi:10.2337/diabetes.55.02.06.db05-0946

Artner, I., Blanchi, B., Raum, J. C., Guo, M., Kaneko, T., Cordes, S., Sieweke, M. and Stein, R. (2007). MafB is required for islet beta cell maturation. *Proc. Natl. Acad. Sci. USA* **104**, 3853-3858. doi:10.1073/pnas.0700013104

Artner, I., Hang, Y., Mazur, M., Yamamoto, T., Guo, M., Lindner, J., Magnuson, M. A. and Stein, R. (2010). MafA and MafB regulate genes critical to beta-cells in a unique temporal manner. *Diabetes* **59**, 2530-2539. doi:10.2337/db10-0190

Asplund, O., Storm, P., Chandra, V., Ottosson-Laakso, E., Hatem, G., Mansour-Aly, D., Krus, U., Ibrahim, H., Ahlqvist, E., Tuomi, T. et al. (2020). Islet Gene View - a tool to facilitate islet research. *Life Sci. Alliance* **10**, e202201376. doi:10.26508/lsa.202201376

Bankaitis, E. D., Bechard, M. E., Gu, G., Magnuson, M. A. and Wright, C. V. E. (2018). ROCK-nmMyoII, Notch and Neurog3 gene-dosage link epithelial morphogenesis with cell fate in the pancreatic endocrine-progenitor niche. *Development* **145**, dev162115. doi:10.1242/dev.16211

Berger, M., Scheel, D. W., Macias, H., Miyatsuka, T., Kim, H., Hoang, P., Ku, G. M., Honig, G., Liou, A., Tang, Y. et al. (2015). Galphai/o-coupled receptor signaling restricts pancreatic beta-cell expansion. *Proc. Natl. Acad. Sci. USA* **112**, 2888-2893. doi:10.1073/pnas.1319378112

Blanchi, B., Kelly, L. M., Viemari, J. C., Lafon, I., Burnet, H., Bevengut, M., Tillmanns, S., Daniel, L., Graf, T., Hilaire, G. et al. (2003). MafB deficiency causes defective respiratory rhythmogenesis and fatal central apnea at birth. *Nat. Neurosci.* **6**, 1091-1100. doi:10.1038/nn1129

Blodgett, D. M., Nowoselska, A., Afik, S., Pechhold, S., Cura, A. J., Kennedy, N. J., Kim, S., Kucukural, A., Davis, R. J., Kent, S. C. et al. (2015). Novel observations from next-generation RNA sequencing of highly purified human adult and fetal islet cell subsets. *Diabetes* **64**, 3172-3181. doi:10.2337/db15-0039

Bonner-Weir, S., Aguayo-Mazzucato, C. and Weir, G. C. (2016). Dynamic development of the pancreas from birth to adulthood. *Ups. J. Med. Sci.* **121**, 155-158. doi:10.3109/03009734.2016.1154906

Borden, P., Houtz, J., Leach, S. D. and Kuruvilla, R. (2013). Sympathetic innervation during development is necessary for pancreatic islet architecture and functional maturation. *Cell Rep* **4**, 287-301. doi:10.1016/j.celrep.2013.06.019

Caputi, V., Marsilio, I., Cerantola, S., Roozfarakh, M., Lante, I., Galuppini, F., Rugge, M., Napoli, E., Giulivi, C., Orso, G. et al. (2017). Toll-like receptor 4 modulates small intestine neuromuscular function through nitrergic and purinergic pathways. *Front Pharmacol* **8**, 350. doi:10.3389/fphar.2017.00350

Conrad, E., Dai, C., Spaeth, J., Guo, M., Cyphert, H. A., Scoville, D., Carroll, J., Yu, W. M., Goodrich, L. V., Harlan, D. M. et al. (2016). The MAFB transcription factor impacts islet alpha-cell function in rodents and represents a unique signature of primate islet beta-cells. *Am. J. Physiol. Endocrinol. Metab.* **310**, E91-E102. doi:10.1152/ajpendo.00285.2015

Cyphert, H. A., Walker, E. M., Hang, Y., Dhawan, S., Haliyur, R., Bonatakis, L., Avrahami, D., Brissova, M., Kaestner, K. H., Bhushan, A. et al. (2019). Examining how the MAFB transcription factor affects islet beta-cell function postnatally. *Diabetes* **68**, 337-348. doi:10.2337/db18-0903
De Vas, M. G., Kopp, J. L., Heliot, C., Sander, M., Cereghini, S. and Haumaitre, C. (2015). Hnf1b controls pancreas morphogenesis and the generation of Ngn3+ endocrine progenitors. *Development* **142**, 871-882. doi:10.1242/dev.110759

Fadista, J., Vikman, P., Laakso, E. O., Mollet, I. G., Esguerra, J. L., Taneera, J., Storm, P., Osmark, P., Ladenvall, C., Prasad, R. B. et al. (2014). Global genomic and transcriptomic analysis of human pancreatic islets reveals novel genes influencing glucose metabolism. *Proc. Natl. Acad. Sci. USA* **111**, 13924-13929. doi:10.1073/pnas.1402665111

Ganic, E., Singh, T., Luan, C., Fadista, J., Johansson, J. K., Cyphert, H. A., Bennet, H., Storm, P., Prost, G., Ahlenius, H. et al. (2016). MafA-Controlled Nicotinic Receptor Expression Is Essential for Insulin Secretion and Is Impaired in Patients with Type 2 Diabetes. *Cell Rep* **14**, 1991-2002. doi:10.1016/j.celrep.2016.02.002

Gautier, L., Cope, L., Bolstad, B. M. and Irizarry, R. A. (2004). affy-analysis of Affymetrix GeneChip data at the probe level. *Bioinformatics* **20**, 307-315. doi:10.1093/bioinformatics/btg405

Glawe, J. D., Mijalis, E. M., Davis, W. C., Barlow, S. C., Gungor, N., McVie, R. and Kevil, C. G. (2013). SDF-1-CXCR4 differentially regulates autoimmune diabetogenic T cell adhesion through to ROBO2-SLIT2 interactions in mice. *Diabetologia* **56**, 2222-2230. doi:10.2337/db15-0039

Glebova, N. O. and Ginty, D. D. (2005). Growth and survival signals controlling sympathetic nervous system development. *Annu. Rev. Neurosci.* **28**, 191-222. doi:10.1146/annurev.neuro.28.061604.135659

Hara, M., Wang, X., Kawamura, T., Bindokas, V. P., Dizon, R. F., Alcoser, S. Y., Magnuson, M. A. and Bell, G. I. (2003). Transgenic mice with green fluorescent protein-labeled pancreatic beta -cells. *Am. J. Physiol. Endocrinol. Metab.* **284**, E177-E183. doi:10.1152/ajpendo.00321.2002

Johnston, N. R., Mitchell, R. K., Haythorne, E., Pessoa, M. P., Semplici, F., Ferrer, J., Piemonti, L., Marchetti, P., Bugliani, M., Bosco, D. et al. (2016). Beta Cell Hubs dictate pancreatic beta cell responses to glucose. *Cell Metab.* **24**, 389-401. doi:10.1016/j.cmet.2016.06.020

Krivova, Y., Proshchina, A., Barabanov, V., Leonova, O. and Saveliev, S. (2016). Structure of neuro-endocrine and neuro-epithelial interactions in human foetal pancreas. *Tissue Cell* **48**, 567-576. doi:10.1016/j.tice.2016.10.005

Lang, S., Ugale, A., Erlandsson, E., Karlsson, G., Bryder, D. and Soneji, S. (2015). SCExV: a webtool for the analysis and visualisation of single cell qRT-PCR data. *BMC Bioinformatics* **16**, 320. doi:10.1186/s12859-015-0757-z

Larsen, H. L., Martin-Coll, L., Nielsen, A. V., Wright, C. V. E., Trusina, A., Kim, Y. H. and Grapin-Botton, A. (2017). Stochastic priming and spatial cues orchestrate heterogeneous clonal contribution to mouse pancreas organogenesis. *Nat. Commun.* **8**, 605. doi:10.1038/s41467-017-00258-4

Liao, Y., Smyth, G. K. and Shi, W. (2014). featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. *Bioinformatics* **30**, 923-930. doi:10.1093/bioinformatics/btt656

Magenheim, J., Klein, A. M., Stanger, B. Z., Ashery-Padan, R., Sosa-Pineda, B., Gu, G. and Dor, Y. (2011). Ngn3(+)-endocrine progenitor cells control the fate and morphogenesis of pancreatic ductal epithelium. *Dev. Biol.* **359**, 26-36. doi:10.1016/j.ydbio.2011.08.006

Matsuoka, T. A., Zhao, L., Artner, I., Jarrett, H. W., Friedman, D., Means, A. and Stein, R. (2003). Members of the large Maf transcription family regulate insulin gene transcription in islet beta cells. *Mol. Cell. Biol.* **23**, 6049-6062. doi:10.1128/ MCB.23.17.6049-6062.2003

Mazur, M. A., Winkler, M., Ganic, E., Colberg, J. K., Johansson, J. K., Bennet, H., Fex, M., Nuber, U. A. and Artner, I. (2013). Microphthalmia transcription factor regulates pancreatic beta-cell function. *Diabetes* **62**, 2834-2842. doi:10.2337/db12-1464

Mularoni, L., Ramos-Rodriguez, M. and Pasquali, L. (2017). The pancreatic islet regulome browser. *Front. Genet.* **8**, 13. doi:10.3389/fgenet.2017.00013

Nishimura, W., Kondo, T., Salameh, T., El Khattabi, I., Dodge, R., Bonner-Weir, S. and Sharma, A. (2006). A switch from MafB to MafA expression accompanies differentiation to pancreatic beta-cells. *Dev. Biol.* **293**, 526-539. doi:10.1016/j.ydbio.2006.02.028

Nishimura, W., Rowan, S., Salameh, T., Maas, R. L., Bonner-Weir, S., Sell, S. M. and Sharma, A. (2008). Preferential reduction of beta cells derived from Pax6-MafB pathway in MafB deficient mice. *Dev. Biol.* **314**, 443-456. doi:10.1016/ j.ydbio.2007.12.009

Pauerstein, P. T., Tellez, K., Willmarth, K. B., Park, K. M., Hsueh, B., Efsun Arda, H., Gu, X., Aghajanian, H., Deisseroth, K., Epstein, J. A. et al. (2017). A radial axis defined by semaphorin-to-neuropilin signaling controls pancreatic islet morphogenesis. *Development* **144**, 3744-3754. doi:10.1242/dev.148684

Pierreux, C. E., Poll, A. V., Kemp, C. R., clotman, F., Maestro, M. A., Cordi, S., Ferrer, J., Leyns, L., Rousseau, G. G. and Lemaigre, F. P. (2006). The transcription factor hepatocyte nuclear factor-6 controls the development of pancreatic ducts in the mouse. *Gastroenterology* **130**, 532-541. doi:10.1053/ j.gastro.2005.12.005

Prokopenko, I., Poon, W., Magi, R., Prasad, B. R., Salehi, S. A., Almgren, P., Osmark, P., Bouatia-Naji, N., Wierup, N., Fall, T. et al. (2014). A central role for

GRB10 in regulation of islet function in man. *PLoS Genet.* **10**, e1004235. doi:10.1371/journal.pgen.1004235

Ritchie, M. E., Phipson, B., Wu, D., Hu, Y., Law, C. W., Shi, W. and Smyth, G. K. (2015). limma powers differential expression analyses for RNA-sequencing and microarray studies. *Nucleic Acids Res.* **43**, e47. doi:10.1093/nar/gkv007

Rodriguez-Diaz, R., Abdulreda, M. H., Formoso, A. L., Gans, I., Ricordi, C., Berggren, P. O. and Acicedo, A. (2011a). Innervation patterns of autonomic axons in the human endocrine pancreas. *Cell Metab.* **14**, 45-54. doi:10.1016/ j.cmet.2011.05.008

Rodriguez-Diaz, R., Dando, R., Jacques-Silva, M. C., Fachado, A., Molina, J., Abdulreda, M. H., Ricordi, C., Roper, S. D., Berggren, P. O. and Acicedo, A. (2011b). Alpha cells secrete acetylcholine as a non-neuronal paracrine signal priming beta cell function in humans. *Nat. Med.* **17**, 888-892. doi:10.1038/nm.2371

Rorsman, P. and Ashcroft, F. M. (2018). Pancreatic beta-cell electrical activity and insulin secretion: of mice and men. *Physiol. Rev.* **98**, 117-214. doi:10.1152/ physrev.00008.2017

Rosengren, A. H., Jokubka, R., Tojjar, D., Granhall, C., Hansson, O., Li, D. Q., Nagaraj, V., Reinbothe, T. M., Tuncel, J., Eliasson, L. et al. (2010). Overexpression of alpha2A-adrenergic receptors contributes to type 2 diabetes. *Science* **327**, 217-220. doi:10.1126/science.1176827

Russell, R., Carnese, P. P., Hennings, T. G., Walker, E. M., Russ, H. A., Liu, J. S., Giacometti, S., Stein, R. and Hebrok, M. (2020). Loss of the transcription factor MAFB limits beta-cell derivation from human PSCs. *Nat. Commun.* **11**, 2742. doi:10.1038/s41467-020-16550-9

Satija, R., Farrell, J. A., Gennert, D., Schier, A. F. and Regev, A. (2015). Spatial reconstruction of single-cell gene expression data. *Nat. Biotechnol.* **33**, 495-502. doi:10.1038/nbt.3192

Shabalin, A. A. (2012). Matrix eQTL: ultra fast eQTL analysis via large matrix operations. *Bioinformatics* **28**, 1353-1358. doi:10.1093/bioinformatics/bts163

Sharon, N., Chawla, R., Mueller, J., Vanderhooft, J., Whitehorn, L. J., Rosenthal, B., Gurtler, M., Estanboulieh, R. R., Shvartsman, D., Gifford, D. K. et al. (2019). A peninsular structure coordinates asynchronous differentiation with morphogenesis to generate pancreatic islets. *Cell* **176**, 790-804.e713. doi:10.1016/j.cell.2018.12.003

Singh, T., Sarmiento, L., Luan, C., Prasad, R. B., Johansson, J., Cataldo, L. R., Renstrom, E., Soneji, S., Cilio, C. and Artner, I. (2018). MafA expression preserves immune homeostasis in human and mouse islets. *Genes (Basel)* **9**, 644. doi:10.3390/genes9120644

Singh, T., Colberg, J. K., Sarmiento, L., Chaves, P., Hansen, L., Bsharat, S., Cataldo, L. R., Dudenhoffer-Pfeifer, M., Fex, M., Bryder, D. et al. (2019). Loss of MafA and MafB expression promotes islet inflammation. *Sci. Rep.* **9**, 9074. doi:10.1038/s41598-019-45528-x

Somm, E., Schwitzgebel, V. M., Vauthay, D. M., Camm, E. J., Chen, C. Y., Giacobino, J. P., Sizonenko, S. V., Aubert, M. L. and Huppi, P. S. (2008). Prenatal nicotine exposure alters early pancreatic islet and adipose tissue development with consequences on the control of body weight and glucose metabolism later in life. *Endocrinology* **149**, 6289-6299. doi:10.1210/en.2008- 0361

Staels, W., Heremans, Y., Heimberg, H. and De Leu, N. (2019). VEGF-A and blood vessels: a beta cell perspective. *Diabetologia* **62**, 1961-1968. doi:10.1007/ s00125-019-4969-z

Sznurkowska, M. K., Hannezo, E., Azzarelli, R., Chatzeli, L., Ikeda, T., Yoshida, S., Philpott, A. and Simons, B. D. (2020). Tracing the cellular basis of islet specification in mouse pancreas. *Nat. Commun.* **11**, 5037. doi:10.1038/ s41467-020-18837-3

Tournier, J. M., Maouche, K., Coraux, C., Zahm, J. M., Cloez-Tayarani, I., Nawrocki-Raby, B., Bonnomet, A., Burlet, H., Lebargy, F., Polette, M. et al. (2006). alpha3alpha5beta2-Nicotinic acetylcholine receptor contributes to the wound repair of the respiratory epithelium by modulating intracellular calcium in migrating cells. *Am. J. Pathol.* **168**, 55-68. doi:10.2353/ajpath.2006.050333

Triolo, T. M. and Bellin, M. D. (2021). Lessons from human islet transplantation inform stem cell-based approaches in the treatment of diabetes. *Front. Endocrinol. (Lausanne)* **12**, 636824. doi:10.3389/fendo.2021.636824

Westmoreland, J. J., Kilic, G., Sartain, C., Sirma, S., Blain, J., Rehg, J., Harvey, N. and Sosa-Pineda, B. (2012). Pancreas-specific deletion of Prox1 affects development and disrupts homeostasis of the exocrine pancreas. *Gastroenterology* **142**, 999-1009. doi:10.1053/j.gastro.2011.12.007

Wood, A. R., Jonsson, A., Jackson, A. U., Wang, N., van Leewen, N., Palmer, N. D., Kobes, S., Deelen, J., Boquete-Vilarino, L., Paananen, J. et al. (2017). A genome-wide association study of IVGTT-based measures of first-phase insulin secretion refines the underlying physiology of type 2 diabetes variants. *Diabetes* **66**, 2296-2309. doi:10.2337/db16-1452

Zhang, C., Moriguchi, T., Kajihara, M., Esaki, R., Harada, A., Shimohata, H., Oishi, H., Hamada, M., Morito, N., Hasegawa, K. et al. (2005). MafA is a key regulator of glucose-stimulated insulin secretion. *Mol. Cell. Biol.* **25**, 4969-4976. doi:10.1128/MCB.25.12.4969-4976.2005