XBART50354
Cell Rep
March 17, 2015;
10
(10):
16461654.
Axis Patterning by BMPs: Cnidarian Network Reveals Evolutionary Constraints.
Genikhovich G
,
Fried P
,
Prünster MM
,
Schinko JB
,
Gilles AF
,
Fredman D
,
Meier K
,
Iber D
,
Technau U
.
Abstract
BMP signaling plays a crucial role in the establishment of the dorso
ventral body axis in bilaterally symmetric animals. However, the topologies of the
bone morphogenetic protein (BMP) signaling networks vary drastically in different animal groups, raising questions about the evolutionary constraints and evolvability of BMP signaling systems. Using lossoffunction analysis and mathematical modeling, we show that two signaling centers expressing different BMPs and BMP antagonists maintain the secondary axis of the sea anemone Nematostella. We demonstrate that BMP signaling is required for asymmetric Hox gene expression and
mesentery formation. Computational analysis reveals that network parameters related to
BMP4 and
Chordin are constrained both in Nematostella and Xenopus, while those describing the BMP signaling modulators can vary significantly. Notably, only
chordin, but not
bmp4 expression needs to be spatially restricted for robust signaling gradient formation. Our data provide an explanation of the evolvability of BMP signaling systems in axis formation throughout Eumetazoa.
PubMed ID:
25772352
PMC ID:
PMC4460265
Article link:
Cell Rep
Genes referenced:
admp
bambi
bmp1
bmp2
bmp4
bmp5
bmp7.1
bmp7.2
chrd.1
gdf5
rpsa
smad6
Article Images:
[+] show captions



Figure 1. BMP Signaling Is Strongest on GDF5likeExpressing Side of Nematostella Embryo(A) Positions of chordin (blue) expression, bmp4 (green) expression, and BMP signaling domain (red circles) in different animal models.(B and C) Schematic representation of NvDpp, NvBMP58, NvGDF5like, NvChd, NvGrm, NvGbx, and Hox expression domains in planula viewed laterally and orally. Red lines, cutting planes; black doubleheaded arrows, directive axis; asterisks, blastopore.(D–K) The αpSMAD1/5 and αNvHoxE antibody staining in control and morphant early planulae, n > 50 for each sample; (D–I) lateral views; (J and K) oral views; asterisks, blastopore. (D) αpSMAD1/5positive nuclei are located on NvHoxEexpressing side. αpSMAD1/5 and αNvHoxE stainings partially overlap. (E–K) αpSMAD1/5 and αNvHoxE in StdMO, ChdMO, GrmMO, GDF5lMO, BMP58MO, and DppMO embryos. Staining is absent in DppMO, BMP58MO, and ChdMO and suppressed (white arrow) in GDF5lMO (E–H); the domain showing strong staining (white doubleheaded arrows and white demarcating lines) is narrower in the StdMO than in the GrmMO (J and K).See also Figure S1.


Figure 2. Effect of MO Knockdown on αpSMAD1/5 Staining Gradient and Morphology of Late Planulae(A–D) Nuclear αpSMAD1/5 staining intensity as a function of relative position of each endodermal nucleus along 180° arc from 0 to π starting at the middle of pSMAD1/5positive domain (see the Supplemental Experimental Procedures and Figure S2A for details on measuring), described by measurements from multiple embryos (colored points) for Control (n = 12), GrmMO (n = 7), GDF5lMO (n = 10), and DppMO (n = 8). LOESS smoothed curves (black lines) with 99% confidence interval for the mean (shade), and linear model fit of the logarithm of relative position to pSMAD1/5 staining intensity (R2 = 0.73, dashed red line) for Control embryos and linear model fit of relative position to pSMAD1/5 staining intensity (R2 = 0.55, dashed red line) for GrmMO embryos. Oral views of corresponding StdMO, DppMO, GDF5lMO, and GrmMO late planulae stained with fluorescent phalloidin are shown next to the graphs (n > 55 for each sample). While StdMO planulae with normal nuclear pSMAD1/5 gradient develop eight mesenteries (red arrows in A), DppMOinjected planulae lacking the pSMAD1/5 gradient do not develop mesenteries at all (B). αpSMAD1/5 staining intensity is suppressed but not absent in the GDF5lMO planulae correlating with the formation of four mesenteries instead of eight (red arrows in C). In GrmMOtreated embryos (D), the average pSMAD1/5 staining intensity is elevated in midrange positions (see also Figure S2C) and more variable than in Control embryos across the whole range (see also Figure S2D). The endoderm in GrmMO appears compartmentalized but its development is abnormal with extremely thickened mesogloea (green arrow). Scale bars represent 50 μm.See also Figure S2.


Figure 3. Staggered Endodermal Expression of Hox Genes and NvGbx Is Regulated by BMP Signaling(A–DD) Expression of NvHoxB (A–F), NvHoxD (G–L), NvHoxE (M–R), NvHoxF (S–X), and NvGrm (Y–DD) in StdMO, BMP58MO, DppMO, GDF5lMO, ChdMO, and GrmMOinjected embryos at early planula stage. Lateral views; oral end to the left. Embryos on (Y) and (DD) are rotated 90° compared to other stained embryos to make bilateral expression of NvGbx visible.(EE–HH) Ectopic expression of BMPRIQ209 > DEGFP in patches of ectodermal cells in 2day planula results in ectopic activation of NvHoxE protein. (EE) Ectopic pSMAD1/5positive nuclei in the ectoderm of early planula in addition to endogenous staining in the endoderm (blue dots). (FF) In addition to endogenous αNvHoxE staining (blue dots), ectopic expression is observed in pSMAD1/5positive ectodermal nuclei. Apart from specific nuclear NvHoxE signal, the secondary antirat antibody crossreacts with mouse antiGFP antibody staining the BMPRIQ209 > DEGFP fusion protein in the cell membranes. (GG) AntiGFP antibody detects BMPRIQ209 > DEGFP in cell membranes. Antimouse antibody does not crossreact with the rat αNvHoxE, thus no nuclear staining is observed, showing that nuclear staining on (FF) is specific for NvHoxE. (HH) Merged image of (EE–GG). White dots indicate the border between ectoderm and endoderm. White arrows point at two strongly pSMAD1/5positive ectodermal nuclei in transgenic cells, which also ectopically express NvHoxE. Yellow arrows point at a weakly pSMAD1/5positive ectodermal nucleus in a transgenic cell, which does not express detectable amounts of NvHoxE. Scale bars represent 50 μm.See also Figure S3.


Figure 4. Maintenance of the Directive Axis in Nematostella(A) Expression analysis of network components in morphant early planulae; lateral views, oral end is to the left. Scale bars represent 50 μm.(B) Mathematical model correctly predicts expression patterns of all genes shown in (A) (color code as in ISH), as well as pSMAD1/5 localization (shown as heat maps) in the knockdowns. All expression patterns and pSMAD1/5 staining patterns are plotted onto mesogloea and shown as 2D projections of a 3D structure. Expression domains of NvDpp/NvBMP58/NvChd, although not identical, strongly overlap when plotted onto mesogloea. In the model they are controlled by the same mathematical expression and thus plotted together.See also Figures S3 and S4.


Figure 5. Strongly Constrained Core Interactions and Weakly Constrained Modulators of the BMP Signaling System(A and B) Nematostella directive axis network (A) compared to the DV network in Xenopus (B). Blue lines represent transcriptional regulation downstream of BMP signaling and black lines represent proteinprotein interactions (putative in case of Nematostella). Xenopus network is modified after (De Robertis and Colozza, 2013; Paulsen et al., 2011). Overlaid blue and black lines represent shuttle function of Chordin.(C–E) Restricted chordin expression domain (red curves on the left plots) results in BMP signaling on the opposite side of the embryo (relative concentration of BMP ligandreceptor complexes shown as black curves on Cʹ–Eʹ) independent of bmp4 expression localization (blue curve). bmp4 can be expressed uniformly (C), on the chordin side (D), or opposite to it (E). The x axis is normalized with respect to maximal domain length (shown as 2π since the embryo is spherical), and the y axis shows the ratio of ligandbound receptors (LR) with respect to total receptor concentration (RT).(F and G) Constrained core regulators versus weakly constrained modulators of BMP signaling in Nematostella (F) and Xenopus (G).(F) Many parameters describing production, function, and turnover of NvDpp/NvBMP58 and NvChd (blue circles) are strongly constrained, while most of the parameters describing production, function, and turnover of NvGDF5like and NvGrm (orange circles) are weakly constrained. Large orange circle with a 9 corresponds to nine unconstrained parameters related to NvGDF5like and NvGrm.(G) Many parameters describing production, function, and turnover of frog BMP4/BMP7, ADMP/BMP2, Chordin, and Tolloid (teal circles) are strongly constrained, while most of the parameters describing production, function, and turnover of BAMBI and SMAD6/7 (pink circles) are weakly constrained. Large pink circle with an 8 corresponds to eight unconstrained parameters related to BAMBI and SMAD6/7. Fold of possible increase of parameter values (log scale) is plotted on the y axis, and fold of possible decrease of parameter values (log scale) is plotted on the x axis. (For details see the Supplemental Experimental Procedures and Figure S5.)See also Figure S5.

References:
Arendt, 1994,
Pubmed
[]
BenZvi, 2008,
Pubmed,
Xenbase
De Robertis, 2013,
Pubmed,
Xenbase
De Robertis, 2008,
Pubmed
Denker, 2008,
Pubmed
Eldar, 2002,
Pubmed
Finnerty, 2004,
Pubmed
Genikhovich, 2010,
Pubmed
Genikhovich, 2010,
Pubmed
Hejnol, 2009,
Pubmed
Iber, 2007,
Pubmed
Inomata, 2013,
Pubmed,
Xenbase
Inomata, 2008,
Pubmed,
Xenbase
Jaźwińska, 1999,
Pubmed
Khokha, 2005,
Pubmed,
Xenbase
Kuo, 2011,
Pubmed
Lapraz, 2009,
Pubmed
Leclère, 2014,
Pubmed
Lee, 2006,
Pubmed,
Xenbase
Lemaire, 2009,
Pubmed
Matus, 2006,
Pubmed,
Xenbase
Matus, 2006,
Pubmed,
Xenbase
Mizutani, 2005,
Pubmed
Patterson, 2000,
Pubmed
Paulsen, 2011,
Pubmed,
Xenbase
Philippe, 2011,
Pubmed
Piccolo, 1996,
Pubmed,
Xenbase
Piccolo, 1997,
Pubmed,
Xenbase
Plouhinec, 2011,
Pubmed,
Xenbase
Putnam, 2007,
Pubmed
Rentzsch, 2006,
Pubmed
Reversade, 2005,
Pubmed,
Xenbase
Ryan, 2007,
Pubmed
Saina, 2009,
Pubmed
Sidi, 2003,
Pubmed
Wacker, 2004,
Pubmed,
Xenbase
Wang, 2005,
Pubmed
Zakin, 2010,
Pubmed
Özüak, 2014,
Pubmed