Branching morphogenesis sculpts the airway epithelium of the lung into a tree-like structure to conduct air and promote gas exchange after birth. In the avian lung, a series of buds emerges from the dorsal surface of the primary bronchus via monopodial branching to form the conducting airways; anatomically, these buds are similar to those formed by domain branching in the mammalian lung. Here, we show that monopodial branching is initiated by apical constriction of the airway epithelium, and not by differential cell proliferation, using computational modeling and quantitative imaging of embryonic chicken lung explants. Both filamentous actin and phosphorylated myosin light chain were enriched at the apical surface of the airway epithelium during monopodial branching. Consistently, inhibiting actomyosin contractility prevented apical constriction and blocked branch initiation. Although cell proliferation was enhanced along the dorsal and ventral aspects of the primary bronchus, especially before branch formation, inhibiting proliferation had no effect on the initiation of branches. To test whether the physical forces from apical constriction alone are sufficient to drive the formation of new buds, we constructed a nonlinear, three-dimensional finite element model of the airway epithelium and used it to simulate apical constriction and proliferation in the primary bronchus. Our results suggest that, consistent with the experimental results, apical constriction is sufficient to drive the early stages of monopodial branching whereas cell proliferation is dispensable. We propose that initial folding of the airway epithelium is driven primarily by apical constriction during monopodial branching of the avian lung.
The intricate network of branched airways of the lung originates as a simple epithelial tube. A sequence of recursive branching events then transforms this simple tubular geometry into a complex system of airways that conduct air and promote gas exchange postnatally. In the early chicken embryo, this process is monopodial; that is, daughter branches emerge laterally along the length of their parent branch (Locy and Larsell, 1916; Hogan, 1999). This branching strategy generates a structure distinct from that in the mammalian lung, which is formed instead by a combination of lateral and dichotomous branching modes (Metzger et al., 2008). Moreover, because airflow is unidirectional in the lungs of birds, as opposed to tidal in mammals, some of these monopodial branches anastomose to form circulatory loops called air capillaries. Other branches dilate at their distal tips to form large air sacs, which then function like bellows to maintain a constant unidirectional flow.
These dramatically different architectures suggest distinct underlying developmental mechanisms. Yet, remarkably, several signaling pathways, including those downstream of fibroblast growth factor (FGF), bone morphogenetic protein (BMP) and sonic hedgehog (SHH), are conserved across different species and different organ systems (Metzger and Krasnow, 1999; Affolter et al., 2003; Morrisey and Hogan, 2010). This molecular homology underlies the hypothesis that the physical mechanisms that sculpt these branching structures are likewise conserved (Hogan, 1999; Davies, 2002). Branching of many organs is thought to result from local increases in proliferation along the epithelium (Goldin and Wessells, 1979; Goldin and Opperman, 1980; Moore et al., 2002; Moore et al., 2005; Morrisey and Hogan, 2010), but unambiguous experimental evidence for this mechanism is lacking (Wessells, 1970; Nogawa et al., 1998). Whether patterned proliferation is sufficient to generate the observed deformations remains unclear. Recent studies, however, have suggested a possible regulatory role for mechanical forces and cytoskeletal contractility during branching (Michael et al., 2005; Moore et al., 2005; Nelson and Gleghorn, 2012). Quantitative studies focused on the mechanics of branching (Caussinus et al., 2008; Gjorevski and Nelson, 2010) are therefore required to define the physical forces that shape the developing airways.
Here, we examined the relative roles of actomyosin contraction and cell proliferation in the initiation and outgrowth of buds during monopodial branching of the embryonic chicken lung. Combining quantitative image analysis with computational modeling, we found that apical constriction of the airway epithelium is sufficient to drive the changes in cell and tissue shape required for initial budding during monopodial branching. Moreover, we also found a complex pattern of proliferation along the primary bronchus, enriched along the dorsal and ventral aspects of the tube and highest in regions prior to branching. Nonetheless, both experiments and simulations revealed that proliferation, coupled with the forces generated by apical constriction, was not sufficient to fully shape the emerging lung bud. These results challenge the generally accepted view that patterned proliferation drives airway branching morphogenesis, and suggest instead that active shape changes of the epithelium play a major role in the branching process.
MATERIALS AND METHODS
Ex vivo culture of embryonic chicken lungs
Fertilized chicken eggs were obtained from a commercial vendor (Hyline) and incubated in a rotating incubator (GQF Manufacturing Company) maintained at 37°C and 60% humidity until the desired stage. Embryonic lungs were dissected using a pair of forceps in cold PBS supplemented with antibiotics (50 units/ml of penicillin and streptomycin; Invitrogen). Dissected lungs were cultured over a porous membrane (nucleopore polycarbonate track-etch membrane, 8 μm, 25 mm; Whatman) in DMEM/F12 medium (without HEPES) supplemented with 5% fetal bovine serum (FBS, heat inactivated; Atlanta Biologicals) and antibiotics (50 units/ml each of penicillin and streptomycin) (Miura et al., 2009; Gleghorn et al., 2012). To alter the rate of proliferation, we varied the concentration of FBS in the culture medium; these treatments had no effect on the relative pattern of proliferation. The pharmacological inhibitors aphidicolin (Sigma), blebbistatin (Sigma) and SU5402 (Santa Cruz) were added to the culture medium and the phenotype of the lung explants was monitored under brightfield illumination on an inverted microscope (Nikon Ti) using a 2× objective.
Immunofluorescence staining and imaging
In general, dissected lungs were fixed with 4% paraformaldehyde in PBS for 15 minutes at room temperature. For staining of F-actin, 0.25% glutaraldehyde was added to the fixative. For staining of α-tubulin, we used ice-cold Dent’s fixative (1:4 ratio of DMSO:methanol). Fixed lungs were washed with 0.3% Triton-X-100 in PBS and blocked with 10% goat serum. The following primary antibodies were used: anti-pMLC-2 (Ser19, Cell Signaling), anti-LCAM (7D6, Developmental Studies Hybridoma Bank), anti-β-catenin (Sigma), anti-α-tubulin (Sigma) and anti-cleaved caspase-3 (Cell Signaling). After extensive washing, samples were incubated in Alexa Fluor-conjugated secondary antibody or phalloidin (Invitrogen). Proliferating cells were detected using the Click-iT EdU Imaging Kit (Invitrogen). As a positive control for apoptosis, lungs were treated with staurosporine (1 μM; Cell Signaling) prior to fixation. Stained lungs were dehydrated in a methanol series and optically cleared with Murray’s clear (1:2 ratio of benzyl alcohol:benzyl benzoate; Sigma) for confocal imaging. Confocal stacks were collected using a spinning disk confocal (BioRad) fitted to an inverted microscope and analyzed using ImageJ (Schneider et al., 2012).
To simulate lung bud formation, we constructed a three-dimensional (3D) nonlinear finite element model of the primary bronchial tube in COMSOL Multiphysics (Version 4.2a; Burlington, MA, USA) (Fig. 3). As discussed below, each bud showed a similar pattern of actomyosin staining and formed as a repeated structure along the tube. We therefore modeled an individual bud and assumed it to be representative of all buds emerging from the primary bronchus.
The bronchial tube was modeled as a thick-walled, hollow cylinder of length L and thickness t = Ro - Ri, where Ro and Ri represent the outer and inner radii of the tube, respectively (Fig. 3A). Regions of apical constriction were specified on the dorsal and ventral sides of the tube, bounded by either the angle α or the angle β as well as the length, Lbud. These regions did not undergo apical constriction homogeneously. Rather, gradients in apical constriction were specified based on experimentally observed patterns of F-actin staining (supplementary material Fig. S1). In addition, a plane of symmetry was assumed through the mid-plane of the bud (Fig. 3A).
We used a theory for finite volumetric growth to model both apical constriction and cell proliferation in the lung epithelium (Rodriguez et al., 1994). Briefly, active shape changes in the tissue were simulated by decomposing the overall deformation gradient tensor F into a component due to growth (or contraction) and a component due to elastic deformation by F = F*·G, where G is the growth tensor and F* is the elastic deformation gradient tensor. This framework has been used previously to model both cardiac c-looping (Voronov et al., 2004) and head fold formation (Varner et al., 2010) in the chicken embryo, and ventral furrow formation (Muñoz et al., 2007; Muñoz et al., 2010) in Drosophila.
The growth tensor G was then further decomposed into a component due to apical constriction G1 and a component due to cell proliferation G2 by G = G2·G1. To simulate apical constriction, cuboidal elements in the bud-forming region acquired wedge-shaped geometries by constricting isotropically (in the θ and z directions) at their apices and expanding isotropically (in the θ and z directions) at their bases (Fig. 3A). The amount of wedging was specified by the components of G1, which were referred to the cylindrical set of basis vectors er, eθ, ez, such that: and where dGi and dGo are the growth stretch ratios along the inner and outer surfaces of the tube, respectively, and Ri and Ro are the inner and outer radii of the un-deformed epithelial tube, respectively. Here, we considered dGi = -dGo, such that shortening along the apical side of the epithelium produced a compensatory lengthening along the basal side to ensure that cell volume was approximately conserved.
Cell proliferation both within and adjacent to the bud-forming region was specified by: where Gθ and Gz represent the amounts of proliferation oriented along the θ and z directions, respectively. Details for implementing this framework in COMSOL can be found in a previous publication by Taber (Taber, 2008).
As a first approximation, we assumed a pseudoelastic constitutive response for the tissue and considered the following Blatz-Ko strain energy density function: where C and ν are material constants, I *1 is the modified first invariant of the right Cauchy-Green elastic deformation gradient tensor C* = F*T · F*, and J* is the elastic dilatation ratio. For small strains, ν represents Poisson’s ratio, but has no physical meaning for problems involving large deformations (Taber, 2004). In all simulations, we assumed slight material compressibility η = 0.45. To simulate tissue stiffening during apical constriction, C was increased (relative to the stiffness of the tissue in passive regions, Cp) in regions of the tissue actively undergoing shape changes.
The Cauchy stress tensor σ, which depends only on the elastic component of the deformation F*, was then computed by (Taber, 2004):
Image analysis and statistics
Quantitative analysis and projection of confocal images were carried out using ImageJ (Schneider et al., 2012). To quantify proliferating epithelial cells, we first collected confocal stacks of double-stained (LCAM and EdU) lung explants for stages with one to three buds (b1 to b3). Cross-sectional images were reconstructed (Fig. 5B) along the length of the bronchus, and EdU-positive cells within the epithelium were counted approximately every 10 μm along the bronchus (proximal, bud, distal). We counted the proliferating cells in the first bud region (b1) across all stages (one to three buds). To quantify apical constriction, we assumed that morphological differences between multiple buds at an individual stage of development were representative of dynamic changes in an individual bud over time. Quantitative analysis of cell geometries was performed on cells spanning the full thickness of the epithelium within buds b1-b3 in fixed explants at Hamburger and Hamilton (HH) stage 26. Statistical significance was calculated using IBM SPSS statistics 19 (IBM) using the t-test for changes in number of buds and ANOVA for comparisons of different stages and treatments.
Apical constriction of the airway epithelium initiates budding from the primary bronchus
The airway epithelium of the embryonic chicken lung undergoes monopodial branching morphogenesis, which consists of initiation of budding from the primary bronchus followed by bud outgrowth (Locy and Larsell, 1916). This monopodial branching is analogous to domain branching in the embryonic mouse lung, in which daughter branches bud laterally off of the circumference of the parent branch (Metzger et al., 2008). In the embryonic chicken lung, the primary bronchus is a simple epithelial tube surrounded by a thick mesenchyme and an outer layer of mesothelium (Locy and Larsell, 1916). Monopodial budding of the bronchus begins at Hamburger-Hamilton (HH) stage 25 (Fig. 1), and the first three buds (b1-b3) emerge sequentially in a cranial-to-caudal direction over a 24-hour period at precise locations along the dorsal surface of the primary bronchus (Fig. 1B-D) (Gleghorn et al., 2012).
To characterize the cellular forces that operate during the branching process, we examined the role of actomyosin contraction in the epithelium of the primary bronchus. Lungs were stained for filamentous actin (F-actin) and phosphorylated myosin light chain (pMLC) and imaged by confocal microscopy. Three-dimensional (3D) reconstruction of confocal stacks revealed that F-actin was enriched at the apical surface of the emerging bud epithelium (Fig. 2A; supplementary material Movie 1), and colocalized with pMLC (Fig. 2B). These data suggested that the apical surface of the dorsal airway epithelium undergoes constriction during new bud formation. Surprisingly, we also observed enrichment of F-actin and pMLC at the apical surface of the non-budding epithelium on the ventral side of the bronchial tube (Fig. 2A, yz slice), which formed a figure-eight-shaped cross-section at the earliest stages of budding. Staining for the cell-cell adhesion protein LCAM (E-cadherin) (supplementary material Fig. S1E,F) revealed an increase in the density of intercellular adhesions at the apical surface of nascent buds. Quantitative analysis of cell geometries within buds at different stages of extension showed that the apical cross-sectional width decreased significantly (Fig. 2C-G), confirming constriction of the apical membrane during the budding process.
To determine whether apical constriction is required to fold the airway epithelium during monopodial budding, we selectively inhibited actomyosin contractility using the non-muscle myosin ATPase inhibitor blebbistatin (Straight et al., 2003; Kovács et al., 2004). Treating embryonic lung explants with blebbistatin (5 μM) significantly reduced the initiation of new buds over a 24-hour period of culture (Fig. 2H,I). Branching morphogenesis of the avian (Moura et al., 2011) and mammalian (Bellusci et al., 1997; Abler et al., 2009) lungs is regulated via signaling through FGF, which has also been shown to induce actomyosin contraction during morphogenesis of the inner ear (Sai and Ladher, 2008) and apical constriction of the rosettes that form the lateral line in zebrafish (Harding and Nechiporuk, 2012). Consistently, we found that treatment with an FGF receptor kinase inhibitor, SU5402 (5 μM), blocked apical constriction (supplementary material Fig. S2) and epithelial budding from the primary bronchi (Fig. 2H,J). These results suggest that apical constriction is one of the physical deformations required to initiate folding of the embryonic airway epithelium into a nascent bud during monopodial branching.
Apical constriction can only create the initial geometry of the bud
To examine whether the forces generated during apical constriction were by themselves sufficient to physically drive monopodial budding, we constructed a 3D model of the primary bronchus using the finite element method (FEM). As described in Materials and methods, we modeled the airway epithelium as a uniform, thick-walled cylinder (Fig. 3A). Active shape changes in the tissue were simulated using a continuum theory for finite volumetric growth (Rodriguez et al., 1994). Consistent with the accumulations of F-actin, pMLC and LCAM observed experimentally, apical constriction was specified in the bud-forming region of the primary bronchus (bounded by the angle α and length Lbud), as well as along the ventral aspect of the epithelial tube immediately opposite the forming bud (bounded by the angle β and length Lbud). In these regions, cuboidal elements in the tissue actively became wedge-shaped (Fig. 3A), and the amount of apical constriction was defined by specifying the components of the growth tensor G1. In longitudinal section, the average magnitude of G1 along the apical surface was 0.49, which was consistent with our analysis of cell geometries (Fig. 2C,E). Moreover, because cells and tissues stiffen when they contract (Wakatsuki et al., 2000; Wakatsuki et al., 2001; Rémond et al., 2006; Zhou et al., 2009; Varner and Taber, 2012), we assumed that the cells undergoing apical constriction during bud formation exhibited a similar material stiffening.
When both apical constriction and the resulting tissue stiffening were included, the model produced a bud geometry similar to that observed experimentally during the initial stages of monopodial branching (compare Fig. 3B,C with Fig. 2A; supplementary material Movie 2). Somewhat surprisingly, the model also predicted the lemniscate (i.e. figure-eight-shaped) geometry of the cross-section of the bud (Fig. 3B), but only if tissue stiffening was included (Fig. 3D). For a given increase in stiffness, however, the extent to which the model was able to produce a lemniscate cross-section depended on the thickness of the tube (supplementary material Fig. S3); the thickness used in our model yielded the most physiologically realistic geometries (Fig. 3D). Still, even for thicker tubes, sufficiently large stiffening could be specified to form a lemniscate cross-section (unpublished data). Even so, our simulations suggested that, at most, apical constriction itself could only initiate the early stages of bud formation (Fig. 3D). On its own, apical constriction was unable to generate the more fully developed geometry of the bud (shown in longitudinal section in Fig. 1C,D), even though spatial gradients of the amount of specified apical constriction were similar qualitatively to the gradients of F-actin intensity measured experimentally (supplementary material Fig. S1).
The results of these simulations suggested that additional cellular behaviors are required to complete the formation of a monopodial bud. As the lung epithelium has been shown to be highly proliferative, we next assessed whether accounting for cell proliferation might produce the more fully developed lung bud geometries observed experimentally.
Cell proliferation is not required for initiation of new buds
The developing chicken lung more than doubles in size during its first 24 hours of culture ex vivo (Gleghorn et al., 2012), as the majority of cells are proliferating during its morphogenesis (Fig. 4A). Differential cell proliferation has been proposed to be essential for bud formation during branching morphogenesis of the mammalian lung and kidney (Goldin et al., 1984; Nogawa et al., 1998; Moore et al., 2002; Michael and Davies, 2004; Moore et al., 2005), although proliferation is apparently not required for branching of the salivary gland (Nakanishi et al., 1987; Spooner et al., 1989). To test whether cell proliferation is necessary for monopodial branching of the avian lung, we blocked progression of the cell cycle using the DNA synthesis inhibitor aphidicolin. Despite complete inhibition of cell proliferation, which was evident after both 3 hours and 24 hours of treatment as visualized by staining for incorporation of the thymidine analog 5-ethynyl-2′-deoxyuridine (EdU) (Fig. 4A,B), the airway epithelium continued to form new buds (Fig. 4C,D). Aphidicolin-treated explants formed buds with essentially the same geometry as those from control explants, and high magnification images revealed that aphidicolin-treated buds comprised fewer but larger cells than control buds (supplementary material Fig. S4A). Consistently, increasing cell proliferation by treating with high concentrations of serum had no effect on monopodial budding (supplementary material Fig. S4B-E). These data suggest that cell proliferation is not required to initiate monopodial branching of the avian airway epithelium.
Combining proliferation with apical constriction is insufficient to fully shape the lung bud
Might cell proliferation play a role during bud outgrowth? To address this question, we quantified the proliferation of the branching airway epithelial cells along the primary bronchus, from stages with one bud to those with three buds. To simplify the analysis, we specified three regions along the primary bronchus: the proximal region, representing the non-budding epithelium; the bud region, representing the actively budding epithelium; and the distal region, representing the future budding epithelium (Fig. 5A,B). Although epithelial cells were proliferating along the entire length of the bronchus, we noted a significant increase in the number of proliferating cells in the distal region of the epithelium, suggesting that epithelial proliferation precedes bud formation (Fig. 5C). Furthermore, analysis of the spatial distribution of proliferating cells within the distal region indicated preferential epithelial growth along the dorsal and ventral sides of the primary bronchus (Fig. 5D). We found that within the epithelium of the bud region, proliferation was confined to the dorsal epithelial cells, which presumably participate in bud outgrowth (Fig. 5D). In the epithelium of the proximal region, the pattern of proliferation varied over time; at the 3-bud stage, the epithelium in this region was largely quiescent (Fig. 5D). We found little discernible apoptosis within the airway epithelium at these stages (supplementary material Fig. S5), consistent with previous reports (Moura et al., 2011).
We then included cell proliferation in our computational model to test whether the combination of proliferation and apical constriction could generate the more fully developed bud morphology observed in experiments (see b1 in Fig. 1C). Consistent with EdU staining in lungs containing one bud (Fig. 5), we modeled cell proliferation in the bud-forming region of the primary bronchus, as well as in the distal region of the tube (Fig. 6A); twice as much growth was included in this adjacent (i.e. distal) region (compare with Fig. 5C). Moreover, because cell proliferation was preferentially confined to the dorsal and ventral sides of the tube (Fig. 5D), we specified growth around the circumference that matched our experimental observations (Fig. 6A, insets).
Although our EdU results provided information about the overall magnitude of proliferation in the epithelium, these data did not indicate whether there was any preferential orientation or directionality. Our model enabled us to examine the relative effects of either oriented or isotropic (i.e. unoriented) growth on bud morphology (Fig. 6A). These patterns of cell division were superimposed on the distribution of apical constriction specified above (Fig. 6A′). We found that only when growth was oriented along the axial length of the tube did a neck-like region begin to form in the developing bud (Fig. 6B,C; supplementary material Movie 3). Specifying axial growth in the bud-region alone was insufficient to generate this bud geometry (supplementary material Fig. S6), thus indicating a possible morphogenetic role for proliferation in regions adjacent to the forming bud. Similarly, isotropic (Fig. 6D,E) and circumferential (Fig. 6F,G) growth merely produced a local distension of the tube and did not qualitatively affect the morphology of the bud.
Our computational model thus predicted that monopodial budding is accompanied by axially directed cell division within the distal region of the primary bronchus. To test these predictions experimentally, we stained lungs at various stages for β-catenin to define cell boundaries and α-tubulin to label the mitotic spindle. Quantitative analysis of spindle orientation within confocal stacks revealed that cells within the distal region divided preferentially along the axial length of the tube (supplementary material Fig. S7). Still, even in the case of axial growth, our computational model was unable to reproduce the more fully developed bud geometry (i.e. a more narrowed neck-region). These data suggest that although proliferation may contribute physically to the budding process, apical constriction and growth alone are insufficient to fully shape the early lung bud. Importantly, further simulations indicated that proliferation was not sufficient to initiate budding in the absence of apical constriction (supplementary material Fig. S8).
Branching morphogenesis of the lung requires reciprocal interactions between the airway epithelium and its surrounding mesenchyme (Metzger and Krasnow, 1999; Affolter et al., 2003). Beginning with the classic experiments of Alescio and Cassini (Alescio and Cassini, 1962), who demonstrated that lung mesenchyme was sufficient to induce supernumerary buds along the otherwise non-budding tracheal epithelium, work over the last few decades has identified several soluble signals secreted by the mesenchyme that direct the formation and outgrowth of incipient buds (Warburton et al., 2000; Morrisey and Hogan, 2010; Ornitz and Yin, 2012). Still, the underlying cellular mechanisms that drive the folding and elongation of the airway epithelium remain unclear. Here, using a combination of theory and experiment, we identified some of the physical and cellular mechanisms that drive monopodial budding of the embryonic chicken lung. Specifically, we found that apical constriction is required to initiate the emergence of new buds from the primary bronchus. Nonetheless, the forces of apical constriction, even when combined with those generated by epithelial cell proliferation, were ultimately insufficient to fully shape the emerging lung bud.
Apical constriction: a conserved physical mechanism for morphogenesis
During embryonic development, simple epithelial sheets bend and fold into a diverse array of tubular structures to accommodate the functions of different organs. Shaping these epithelial tissues requires the generation of forces transmitted between neighboring cells through intercellular polymers and molecular motors (Gorfinkiel and Blanchard, 2011; Kim and Davidson, 2011). The physical mechanisms that transmit and generate these forces include active changes in cell shape (often driven by apical actomyosin contractions) (Haigo et al., 2003; Dawes-Hoang et al., 2005), the addition or deletion of cells via division or apoptosis (Gong et al., 2004; Baena-López et al., 2005; Toyama et al., 2008; Tang et al., 2011), and passive stretching of the tissue in response to either intrinsic or extrinsic forces (Coulombre, 1956; Takeuchi, 1983; Unbekandt et al., 2008).
Apical constriction is an active narrowing of cellular apices that is driven by a contractile actomyosin network, and these changes in cell shape can induce dramatic deformations of the tissue (Sawyer et al., 2010). Many developmental processes involve apical constriction (Taber, 1995; Lecuit and Lenne, 2007; Wozniak and Chen, 2009), including the invagination of epithelial pits, grooves and folds during gastrulation (Martin et al., 2009), neurulation (Colas and Schoenwolf, 2001) and optic cup formation (Plageman et al., 2010). Primary invagination of the Drosophila salivary gland and trachea are also driven by apical constriction (Myat and Andrew, 2000b; Letizia et al., 2011).
It is therefore not entirely surprising that apical constriction is involved in the initiation of monopodial branching of the avian airway, in which the epithelial tube folds to produce a bud along its length. We found a dramatic apical shrinkage during the budding process (20% of basal surface), which appears to be similar in magnitude to that measured during invagination of the vertebrate lens placode (36% of basal surface), a process that requires apical localization of actin-binding proteins (Plageman et al., 2010). These numbers are congruent with those measured for apical constriction during Drosophila gastrulation, in which the area of the apical surface of initially columnar cells decreases to ∼40% of its original value, but with a simultaneous lengthening of the cells (Gelbart et al., 2012). The mechanics of monopodial budding also appear to be distinct from ascidian gastrulation, as invagination of the endoderm involves a substantial cell shortening along the apico-basal axis (Sherrard et al., 2010), which stands in contrast to our observations. Whether these various apical constriction events involve the same underlying regulatory mechanisms remains unclear.
The observation that the FGF receptor inhibitor SU5402 blocked apical constriction and new bud formation in cultured chicken lung explants suggests that sustained apical actomyosin contraction occurs downstream of FGF signaling. Similar to lung development in mice (Bellusci et al., 1997; Park et al., 1998), FGF10 is expressed focally in the embryonic chicken lung in regions of mesenchyme adjacent to emerging buds (Sakiyama et al., 2003; Moura et al., 2011). A feedback mechanism between FGF10 and the transcription factors Tbx4 in the mesenchyme and Nkx2.1 in the epithelium has been shown to regulate monopodial branching in the developing chicken lung (Sakiyama et al., 2003). This pathway might play a role in the apical location of actomyosin structures or in the contractile activity of myosin.
Patterns of apical constriction suggest a more complex regulatory mechanism
Surprisingly, we found that F-actin and pMLC accumulated not only within the dorsal bud-forming region of the primary bronchus, but also along the opposite ventral side of the tube. Moreover, our simulations revealed that apical constriction of both the dorsal and ventral sides of the tube, as well as the tissue stiffening that is likely to result from this constriction, are required to form the leminiscate cross-sectional geometry of the bud. These results complicate the generally accepted view that local epithelial-mesenchymal interactions in the bud-forming region alone drive branching morphogenesis, as forces generated by adjacent regions of the tissue contribute physically to this process. The molecular signals that regulate the ventral apical constriction, however, remain unclear; it will be interesting to determine whether these are distinct from the mesenchymal signals that induce apical constriction of the dorsal side, and whether a similar dorsoventral constriction drives domain branching in mammalian lungs.
Monopodial branching is not driven by differential cell proliferation
The addition of cells through proliferation is thought to regulate many morphogenetic movements, including branching of the lung (Nogawa et al., 1998; Tang et al., 2011), and descriptions of the branching process typically suggest that bud outgrowth results from increased proliferation in bud-forming regions (Goldin et al., 1984; Nogawa et al., 1998; Menshykau et al., 2012). Surprisingly, we found that proliferation is neither necessary nor sufficient for bud initiation. That buds can initiate in the absence of proliferation supports the idea that localized changes in cell and tissue shape (in part due to apical constriction) drive the initial stages of monopodial budding, which is similar to the mechanism underlying salivary branching in Drosophila (Myat and Andrew, 2000a; Myat and Andrew, 2000b). Consistently, elevated NFκB activity in embryonic chicken lung explants increases the number of buds, yet simultaneously represses proliferation in both the epithelium and mesenchyme (Muraoka et al., 2000). Combined with our results, this provides compelling evidence that proliferation is unnecessary for initiation of monopodial budding of the avian lung.
Nonetheless, our quantitative analysis of EdU-positive cells revealed distinct patterns of cell proliferation within the airway epithelium of the primary bronchus, including elevated proliferation prior to budding and a dorsoventral concentration of proliferation, suggesting that proliferation might contribute to later stages of the branching process. Moreover, differential cell proliferation in the epithelium was observed before the initiation of new buds, as cells distal to emerging buds divided at a higher rate. This stands in contrast to the observations of Nogawa and colleagues, who noted that differential cell proliferation appeared only after the bud had formed in mesenchyme-free embryonic mouse lung explants (Nogawa et al., 1998). This may suggest species-specific differences in monopodial branching, or perhaps even a difference between lungs cultured in mesenchyme-free assays and those that form in the intact embryo.
Our modeling results further demonstrated that differential cell proliferation by itself was insufficient to form a lung bud. Nonetheless, our simulations suggested that proliferation (or some other kind of active elongation in the tissue, e.g. anisotropic changes in cell shape), although not sufficient for bud formation, may still contribute physically to this process, but only if this growth is oriented along the axis of the tube and occurs both within and adjacent to the developing bud. This result is consistent with recently reported evidence of oriented cell divisions driving branch elongation in the developing mouse lung (Tang et al., 2011) and, moreover, suggests that physical forces outside the bud-forming regions play a mechanistic role in this process. Quantitative analysis of the cell division axis within the distal epithelium supports these predictions (supplementary material Fig. S7). Surprisingly, this same analysis revealed that although the dorsal epithelium of the bud does not substantially increase in thickness during monopodial branching (Fig. 2F), a large number of cells within this region preferentially divide perpendicular to the plane of the bud (supplementary material Fig. S7). Such perpendicular cell divisions have been postulated to drive bud extension during dichotomous branching of the mouse lung (Schnatwinkel and Niswander, 2013). Given that the epithelium does not thicken during monopodial branching in the chicken, it is more likely that the cells intercalate within the plane of the bud after cytokinesis. Future quantitative imaging studies are needed to decipher the commonalities and differences between monopodial and dichotomous branching across species.
The inability of both apical constriction and proliferation to generate the fully developed geometry of the bud also suggests a possible physical role for the mesenchyme, as has been suggested for the developing mouse lung (Blanc et al., 2012). Combining experimental and computational studies in the future may help to parse the relative roles of cytoskeletal contractility, oriented cell proliferation, and differential mechanical properties of the epithelium and mesenchyme during airway branching morphogenesis.
We thank the members of the Nelson group for their insightful comments.
This work was supported in part by grants from the National Institutes of Health (NIH) [GM083997 and HL110335]; the David and Lucile Packard Foundation; the Alfred P. Sloan Foundation; the Camille and Henry Dreyfus Foundation; and Susan G. Komen for the Cure. C.M.N. holds a Career Award at the Scientific Interface from the Burroughs Wellcome Fund. Deposited in PMC for release after 12 months.
Competing interests statement
The authors declare no competing financial interests.
H.Y.K., V.D.V. and C.M.N. designed the study. H.Y.K. completed experiments and quantitative imaging. V.D.V. developed the computational model and analyzed results. H.Y.K., V.D.V. and C.M.N. interpreted data and wrote the manuscript.
Supplementary material available online at http://dev.biologists.org/lookup/suppl/doi:10.1242/dev.093682/-/DC1
- Accepted May 14, 2013.
- © 2013. Published by The Company of Biologists Ltd