ΔNp63 regulates Sfrp1 expression to direct salivary gland branching morphogenesis.

Branching morphogenesis is a complex process shared by many organs including the lungs, kidney, prostate, as well as several exocrine organs including the salivary, mammary and lacrimal glands. This critical developmental program ensures the expansion of an organ's surface area thereby maximizing processes of cellular secretion or absorption. It is guided by reciprocal signaling from the epithelial and mesenchymal cells. While signaling pathways driving salivary gland branching morphogenesis have been relatively well-studied, our understanding of the underlying transcriptional regulatory mechanisms directing this program, is limited. Here, we performed in vivo and ex vivo studies of the embryonic mouse submandibular gland to determine the function of the transcription factor ΔNp63, in directing branching morphogenesis. Our studies show that loss of ΔNp63 results in alterations in the differentiation program of the ductal cells which is accompanied by a dramatic reduction in branching morphogenesis that is mediated by dysregulation of WNT signaling. We show that ΔNp63 modulates WNT signaling to promote branching morphogenesis by directly regulating Sfrp1 expression. Collectively, our findings have revealed a novel role for ΔNp63 in the regulation of this critical process and offers a better understanding of the transcriptional networks involved in branching morphogenesis.

Maternal Diet High in Linoleic Acid Alters Renal Branching Morphogenesis and mTOR/AKT Signalling Genes in Rat Fetal Kidneys.

Linoleic acid (LA), an n-6 polyunsaturated fatty acid (PUFA), is obtained from the maternal diet during pregnancy, and is essential for normal fetal growth and development. A maternal high-LA (HLA) diet alters maternal and offspring fatty acids, maternal leptin and male/female ratio at embryonic (E) day 20 (E20). We investigated the effects of an HLA diet on embryonic offspring renal branching morphogenesis, leptin signalling, megalin signalling and angiogenesis gene expression. Female Wistar Kyoto rats were fed low-LA (LLA; 1.44% energy from LA) or high-LA (HLA; 6.21% energy from LA) diets during pregnancy and gestation/lactation. Offspring were sacrificed and mRNA from kidneys was analysed by real-time PCR. Maternal HLA decreased the targets involved in branching morphogenesis Ret and Gdnf in offspring, independent of sex. Furthermore, downstream targets of megalin, namely mTOR, Akt3 and Prkab2, were reduced in offspring from mothers consuming an HLA diet, independent of sex. There was a trend of an increase in the branching morphogenesis target Gfra1 in females (p = 0.0517). These findings suggest that an HLA diet during pregnancy may lead to altered renal function in offspring. Future research should investigate the effects an HLA diet has on offspring kidney function in adolescence and adulthood.

Dynamic Hippo pathway activity underlies mesenchymal differentiation during lung alveolar morphogenesis.

Alveologenesis, the final stage in lung development, substantially remodels the distal lung, expanding the alveolar surface area for efficient gas exchange. Secondary crest myofibroblasts (SCMF) exist transiently in the neonatal distal lung and are crucial for alveologenesis. However, the pathways that regulate SCMF function, proliferation and temporal identity remain poorly understood. To address this, we purified SCMFs from reporter mice, performed bulk RNA-seq and found dynamic changes in Hippo-signaling components during alveologenesis. We deleted the Hippo effectors Yap/Taz from Acta2-expressing cells at the onset of alveologenesis, causing a significant arrest in alveolar development. Using single cell RNA-seq, we identified a distinct cluster of cells in mutant lungs with altered expression of marker genes associated with proximal mesenchymal cell types, airway smooth muscle and alveolar duct myofibroblasts. In vitro studies confirmed that Yap/Taz regulates myofibroblast-associated gene signature and contractility. Together, our findings show that Yap/Taz is essential for maintaining functional myofibroblast identity during postnatal alveologenesis.

Molecular control of cellulosic fin morphogenesis in ascidians.

BACKGROUND: The tunicates form a group of filter-feeding marine animals closely related to vertebrates. They share with them a number of features such as a notochord and a dorsal neural tube in the tadpole larvae of ascidians, one of the three groups that make tunicates. However, a number of typical chordate characters have been lost in different branches of tunicates, a diverse and fast-evolving phylum. Consequently, the tunic, a sort of exoskeleton made of extracellular material including cellulose secreted by the epidermis, is the unifying character defining the tunicate phylum. In the larva of ascidians, the tunic differentiates in the tail into a median fin (with dorsal and ventral extended blades) and a caudal fin. RESULTS: Here we have performed experiments in the ascidian Phallusia mammillata to address the molecular control of tunic 3D morphogenesis. We have demonstrated that the tail epidermis medio-lateral patterning essential for peripheral nervous system specification also controls tunic elongation into fins. More specifically, when tail epidermis midline identity was abolished by BMP signaling inhibition, or CRISPR/Cas9 inactivation of the transcription factor coding genes Msx or Klf1/2/4/17, median fin did not form. We postulated that this genetic program should regulate effectors of tunic secretion. We thus analyzed the expression and regulation in different ascidian species of two genes acquired by horizontal gene transfer (HGT) from bacteria, CesA coding for a cellulose synthase and Gh6 coding for a cellulase. We have uncovered an unexpected dynamic history of these genes in tunicates and high levels of variability in gene expression and regulation among ascidians. Although, in Phallusia, Gh6 has a regionalized expression in the epidermis compatible with an involvement in fin elongation, our functional studies indicate a minor function during caudal fin formation only. CONCLUSIONS: Our study constitutes an important step in the study of the integration of HGT-acquired genes into developmental networks and a cellulose-based morphogenesis of extracellular material in animals.

Cell proliferation and Notch signaling coordinate the formation of epithelial folds in the Drosophila leg.

The formation of complex three-dimensional organs during development requires precise coordination between patterning networks and mechanical forces. In particular, tissue folding is a crucial process that relies on a combination of local and tissue-wide mechanical forces. Here, we investigate the contribution of cell proliferation to epithelial morphogenesis using the Drosophila leg tarsal folds as a model. We reveal that tissue-wide compression forces generated by cell proliferation, in coordination with the Notch signaling pathway, are essential for the formation of epithelial folds in precise locations along the proximo-distal axis of the leg. As cell numbers increase, compressive stresses arise, promoting the folding of the epithelium and reinforcing the apical constriction of invaginating cells. Additionally, the Notch target dysfusion plays a key function specifying the location of the folds, through the apical accumulation of F-actin and the apico-basal shortening of invaginating cells. These findings provide new insights into the intricate mechanisms involved in epithelial morphogenesis, highlighting the crucial role of tissue-wide forces in shaping a three-dimensional organ in a reproducible manner.

MARK1 regulates dendritic spine morphogenesis and cognitive functions in vivo.

Dendritic spines play a pivotal role in synaptic communication and are crucial for learning and memory processes. Abnormalities in spine morphology and plasticity are observed in neurodevelopmental and neuropsychiatric disorders, yet the underlying signaling mechanisms remain poorly understood. The microtubule affinity regulating kinase 1 (MARK1) has been implicated in neurodevelopmental disorders, and the MARK1 gene shows accelerated evolution in the human lineage suggesting a role in cognition. However, the in vivo role of MARK1 in synaptogenesis and cognitive functions remains unknown. Here we show that forebrain-specific conditional knockout (cKO) of Mark1 in mice causes defects in dendritic spine morphogenesis in hippocampal CA1 pyramidal neurons with a significant reduction in spine density. In addition, we found loss of MARK1 causes synaptic accumulation of GKAP and GluA2. Furthermore, we found that MARK1 cKO mice show defects in spatial learning in the Morris water maze and reduced anxiety-like behaviors in the elevated plus maze. Taken together, our data show a novel role for MARK1 in regulating dendritic spine morphogenesis and cognitive functions in vivo.

Inter-plane feedback coordinates cell morphogenesis and maintains 3D tissue organization in the Drosophila pupal retina.

How complex organs coordinate cellular morphogenetic events to achieve three-dimensional (3D) form is a central question in development. The question is uniquely tractable in the late Drosophila pupal retina, where cells maintain stereotyped contacts as they elaborate the specialized cytoskeletal structures that pattern the apical, basal and longitudinal planes of the epithelium. In this study, we combined cell type-specific genetic manipulation of the cytoskeletal regulator Abelson (Abl) with 3D imaging to explore how the distinct cellular morphogenetic programs of photoreceptors and interommatidial pigment cells (IOPCs) organize tissue pattern to support retinal integrity. Our experiments show that photoreceptor and IOPC terminal differentiation is unexpectedly interdependent, connected by an intercellular feedback mechanism that coordinates and promotes morphogenetic change across orthogonal tissue planes to ensure correct 3D retinal pattern. We propose that genetic regulation of specialized cellular differentiation programs combined with inter-plane mechanical feedback confers spatial coordination to achieve robust 3D tissue morphogenesis.

Embryonic spatiotemporal expression pattern of Folded gastrulation suggests roles in multiple morphogenetic events and regulation by AbdA.

In Drosophila, the signaling pathway activated by the ligand Folded gastrulation (Fog) is among the few known G protein-coupled receptor (GPCR) pathways to regulate cell shape change with a well-characterized role in gastrulation. However, an understanding of the spectrum of morphogenetic events regulated by Fog signaling is still lacking. Here, we present an analysis of the expression pattern and regulation of fog using a genome-engineered Fog::sfGFP line. We show that Fog expression is widespread and in tissues previously not associated with the signaling pathway including germ cells, trachea, and amnioserosa. In the central nervous system (CNS), we find that the ligand is expressed in multiple types of glia indicating a prominent role in the development of these cells. Consistent with this, we have identified 3 intronic enhancers whose expression in the CNS overlaps with Fog::sfGFP. Further, we show that enhancer-1, (fogintenh-1) located proximal to the coding exon is responsive to AbdA. Supporting this, we find that fog expression is downregulated in abdA mutants. Together, our study highlights the broad scope of Fog-GPCR signaling during embryogenesis and identifies Hox gene AbdA as a novel regulator of fog expression.

The BTB transcription factor, Abrupt, acts cooperatively with Chronologically inappropriate morphogenesis (Chinmo) to repress metamorphosis and promotes leg regeneration.

Many insects undergo the process of metamorphosis when larval precursor cells begin to differentiate to create the adult body. The larval precursor cells retain stem cell-like properties and contribute to the regenerative ability of larval appendages. Here we demonstrate that two Broad-complex/Tramtrack/Bric-à-brac Zinc-finger (BTB) domain transcription factors, Chronologically inappropriate morphogenesis (Chinmo) and Abrupt (Ab), act cooperatively to repress metamorphosis in the flour beetle, Tribolium castaneum. Knockdown of chinmo led to precocious development of pupal legs and antennae. We show that although topical application of juvenile hormone (JH) prevents the decrease in chinmo expression in the final instar, chinmo and JH act in distinct pathways. Another gene encoding the BTB domain transcription factor, Ab, was also necessary for the suppression of broad (br) expression in T. castaneum in a chinmo RNAi background, and simultaneous knockdown of ab and chinmo led to the precocious onset of metamorphosis. Furthermore, knockdown of ab led to the loss of regenerative potential of larval legs independently of br. In contrast, chinmo knockdown larvae exhibited pupal leg regeneration when a larval leg was ablated. Taken together, our results show that both ab and chinmo are necessary for the maintenance of the larval tissue identity and, apart from its role in repressing br, ab acts as a crucial regulator of larval leg regeneration. Our findings indicate that BTB domain proteins interact in a complex manner to regulate larval and pupal tissue homeostasis.

Developmental hourglass: Verification by numerical evolution and elucidation by dynamical-systems theory.

Determining the general laws between evolution and development is a fundamental biological challenge. Developmental hourglasses have attracted increased attention as candidates for such laws, but the necessity of their emergence remains elusive. We conducted evolutionary simulations of developmental processes to confirm the emergence of the developmental hourglass and unveiled its establishment. We considered organisms consisting of cells containing identical gene networks that control morphogenesis and evolved them under selection pressure to induce more cell types. By computing the similarity between the spatial patterns of gene expression of two species that evolved from a common ancestor, a developmental hourglass was observed, that is, there was a correlation peak in the intermediate stage of development. The fraction of pleiotropic genes increased, whereas the variance in individuals decreased, consistent with previous experimental reports. Reduction of the unavoidable variance by initial or developmental noise, essential for survival, was achieved up to the hourglass bottleneck stage, followed by diversification in developmental processes, whose timing is controlled by the slow expression dynamics conserved among organisms sharing the hourglass. This study suggests why developmental hourglasses are observed within a certain phylogenetic range of species.

Robust axis elongation by Nodal-dependent restriction of BMP signaling.

Embryogenesis results from the coordinated activities of different signaling pathways controlling cell fate specification and morphogenesis. In vertebrate gastrulation, both Nodal and BMP signaling play key roles in germ layer specification and morphogenesis, yet their interplay to coordinate embryo patterning with morphogenesis is still insufficiently understood. Here, we took a reductionist approach using zebrafish embryonic explants to study the coordination of Nodal and BMP signaling for embryo patterning and morphogenesis. We show that Nodal signaling triggers explant elongation by inducing mesendodermal progenitors but also suppressing BMP signaling activity at the site of mesendoderm induction. Consistent with this, ectopic BMP signaling in the mesendoderm blocks cell alignment and oriented mesendoderm intercalations, key processes during explant elongation. Translating these ex vivo observations to the intact embryo showed that, similar to explants, Nodal signaling suppresses the effect of BMP signaling on cell intercalations in the dorsal domain, thus allowing robust embryonic axis elongation. These findings suggest a dual function of Nodal signaling in embryonic axis elongation by both inducing mesendoderm and suppressing BMP effects in the dorsal portion of the mesendoderm.

GRAS transcription factor PINNATE-LIKE PENTAFOLIATA2 controls compound leaf morphogenesis in Medicago truncatula.

The milestone of compound leaf development is the generation of separate leaflet primordia during the early stages, which involves two linked but distinct morphogenetic events: leaflet initiation and boundary establishment for leaflet separation. Although some progress in understanding the regulatory pathways for each event have been made, it is unclear how they are intrinsically coordinated. Here, we identify the PINNATE-LIKE PENTAFOLIATA2 (PINNA2) gene encoding a newly identified GRAS transcription factor in Medicago truncatula. PINNA2 transcripts are preferentially detected at organ boundaries. Its loss-of-function mutations convert trifoliate leaves into a pinnate pentafoliate pattern. PINNA2 directly binds to the promoter region of the LEAFY orthologue SINGLE LEAFLET1 (SGL1), which encodes a key positive regulator of leaflet initiation, and downregulates its expression. Further analysis revealed that PINNA2 synergizes with two other repressors of SGL1 expression, the BEL1-like homeodomain protein PINNA1 and the C2H2 zinc finger protein PALMATE-LIKE PENTAFOLIATA1 (PALM1), to precisely define the spatiotemporal expression of SGL1 in compound leaf primordia, thereby maintaining a proper pattern of leaflet initiation. Moreover, we showed that the enriched expression of PINNA2 at the leaflet-to-leaflet boundaries is positively regulated by the boundary-specific gene MtNAM, which is essential for leaflet boundary formation. Together, these results unveil a pivotal role of the boundary-expressed transcription factor PINNA2 in regulating leaflet initiation, providing molecular insights into the coordination of intricate developmental processes underlying compound leaf pattern formation.

Ser252Trp mutation in fibroblast growth factor receptor 2 promotes branching morphogenesis in mouse salivary glands.

OBJECTIVES: The purpose of this study was to perform morphological and immunohistochemical (IHC) analysis of the submandibular glands (SMGs) in early development in Apert syndrome model mice (Ap mice). METHODS: ACTB-Cre homozygous mice were mated with fibroblast growth factor receptor 2 (Fgfr2+/Neo-S252W) mice; ACTB-Cre heterozygous mice (ACTB-Cre mice) at embryonic day (E) 13.5 served as the control group, and Fgfr2+/S252W mice (Ap mice) served as the experimental group. Hematoxylin and eosin (H&E) staining was performed on SMGs; Total SMG area and epithelial area were determined, and the epithelial occupancy ratio was calculated. Immunostaining was performed to assess the localization of FGF signaling-related proteins. Next, bromodeoxyuridine (BrdU)-positive cells were evaluated to assess cell proliferation. Finally, terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL) staining was performed to assess apoptosis in SMGs. RESULTS: The epithelial occupancy ratio was significantly higher in SMGs of Ap mice compared with that in SMGs of controls. FGF7 and bone morphogenetic protein 4 (BMP4) exhibited different localizations in SMGs of Ap mice compared with SMGs of controls. Cell proliferation was higher in SMGs of Ap mice compared with that of controls; however, apoptosis did not different significantly between the two groups. CONCLUSION: Our results suggest that enhanced FGF signaling conferred by missense mutations in FGFR2 promotes branching morphogenesis in SMGs of Ap mice.

PLEIOTROPIC REGULATORY LOCUS1 maintains actin microfilament integrity to regulate pavement cell morphogenesis.

Actin dynamics are critical for plant cell morphogenesis, but the underlying signaling mechanisms regulating these dynamics are not well understood. Here, we established that PLEIOTROPIC REGULATORY LOCUS1 (PRL1) modulates leaf pavement cell (PC) morphogenesis in Arabidopsis (Arabidopsis thaliana) by maintaining the dynamic homeostasis of actin microfilaments (MF). Our previous studies indicated that PC shape was determined by antagonistic RHO-RELATED GTPase FROM PLANTS 2 (ROP2) and RHO-RELATED GTPase FROM PLANTS 6 (ROP6) signaling pathways that promote cortical MF and microtubule organization, respectively. Our genetic screen for additional components in ROP6-mediated signaling identified prl1 alleles. Genetic analysis confirmed that PRL1 plays a key role in PC morphogenesis. Mutations in PRL1 caused cortical MF depolymerization, resulting in defective PC morphogenesis. Further genetic analysis revealed that PRL1 is epistatic to ROP2 and ROP6 in PC morphogenesis. Mutations in PRL1 enhanced the effects of ROP2 and ROP6 in PC morphogenesis, leading to a synergistic phenotype in the PCs of ROP2 prl1 and ROP6 prl1. Furthermore, the activities of ROP2 and ROP6 were differentially altered in prl1 mutants, suggesting that ROP2 and ROP6 function downstream of PRL1. Additionally, cortical MF depolymerization in prl1 mutants occurred independently of ROP2 and ROP6, implying that these proteins impact PC morphogenesis in the prl1 mutant through other cellular processes. Our research indicates that PRL1 preserves the structural integrity of actin and facilitates pavement cell morphogenesis in Arabidopsis.

[Research progress in molecular mechanism of palatal development].

Cleft palate is one of the most common maxillofacial birth defects, which can occur alone or accompany with many known deformities. Palatal selves need to complete the process of vertical growth, elevation, adhesion and fusion in a specific time window of embryo development. Any abnormality in this process will lead to cleft palate. Although previous studies have identified many molecular networks that regulate the growth, location and fusion of palatal selves, there are still many unknown mechanisms for palatal development. The pathogenesis of cleft palate has not been clarified so far. In recent years, the molecular research on palate development has been deepened continuously. Here we summarize major recent advances and integrate the genes and molecular pathways with the cellular and morphogenetic processes of palatal shelf growth, patterning, elevation, adhesion, and fusion, in order to comprehensively understand the genotype-phenotype functional relationship and provide assistance in formulating effective prevention strategies for cleft palate disease.

The first embryo, the origin of cancer and animal phylogeny. II. The neoplastic process as an evolutionary engine.

In this article, I put forward the idea that the neoplastic process (NP) has deep evolutionary roots and makes specific predictions about the connection between cancer and the formation of the first embryo, which allowed for the evolutionary radiation of metazoans. My main hypothesis is that NP is at the heart of cellular mechanisms responsible for animal morphogenesis and, given its embryological basis, also at the center of animal evolution. It is thus understood that NP-associated mechanisms are deeply rooted in evolutionary history and tied to the formation of the first animal embryo. In my consideration of these arguments, I expound on how cancer biology is perfectly intertwined with evolutionary biology. Finally, I describe essential cellular components of unicellular holozoans that served as a basis for the formation of the neoplastic functional module (NFM) and its subsequent exaptation, which brought forth two great biophysical revolutions within the first embryo.

Single-Cell Transcriptome Sequence Profiling on the Morphogenesis of Secondary Hair Follicles in Ordos Fine-Wool Sheep.

The Ordos fine-wool sheep is a high-quality breed in China that produces superior natural textiles and raw materials such as wool and lamb meat. However, compared to the Australian Merino sheep, there is still a gap in terms of the wool fiber fineness and wool yield. The hair follicle is the main organ that controls the type of wool fiber, and the morphological changes in the secondary hair follicle are crucial in determining wool quality. However, the process and molecular mechanisms of hair follicle morphogenesis in Ordos fine-wool sheep are not yet clear. Therefore, analyzing the molecular mechanisms underlying the process of follicle formation is of great significance for improving the fiber diameter and wool production of Ordos fine-wool sheep. The differential expressed genes, APOD, POSTN, KRT5, and KRT15, which related to primary hair follicles and secondary hair follicles, were extracted from the dermal papillae. Based on pseudo-time analysis, the differentiation trajectories of dermal lineage cells and epidermal lineage cells in the Ordos fine-wool sheep were successfully constructed, providing a theoretical basis for breeding research in Ordos fine-wool sheep.

Coordination of cell cycle and morphogenesis during organ formation.

Organ formation requires precise regulation of cell cycle and morphogenetic events. Using the Drosophila embryonic salivary gland (SG) as a model, we uncover the role of the SP1/KLF transcription factor Huckebein (Hkb) in coordinating cell cycle regulation and morphogenesis. The hkb mutant SG exhibits defects in invagination positioning and organ size due to the abnormal death of SG cells. Normal SG development involves distal-to-proximal progression of endoreplication (endocycle), whereas hkb mutant SG cells undergo abnormal cell division, leading to cell death. Hkb represses the expression of key cell cycle and pro-apoptotic genes in the SG. Knockdown of cyclin E or cyclin-dependent kinase 1, or overexpression of fizzy-related rescues most of the morphogenetic defects observed in the hkb mutant SG. These results indicate that Hkb plays a critical role in controlling endoreplication by regulating the transcription of key cell cycle effectors to ensure proper organ formation.

Patterned apoptosis has an instructive role for local growth and tissue shape regulation in a fast-growing epithelium.

What regulates organ size and shape remains one fundamental mystery of modern biology. Research in this area has primarily focused on deciphering the regulation in time and space of growth and cell division, while the contribution of cell death has been overall neglected. This includes studies of the Drosophila wing, one of the best-characterized systems for the study of growth and patterning, undergoing massive growth during larval stage and important morphogenetic remodeling during pupal stage. So far, it has been assumed that cell death was relatively neglectable in this tissue both during larval stage and pupal stage, and as a result, the pattern of growth was usually attributed to the distribution of cell division. Here, using systematic mapping and registration combined with quantitative assessment of clone size and disappearance as well as live imaging, we outline a persistent pattern of cell death and clone elimination emerging in the larval wing disc and persisting during pupal wing morphogenesis. Local variation of cell death is associated with local variation of clone size, pointing to an impact of cell death on local growth that is not fully compensated by proliferation. Using morphometric analyses of adult wing shape and genetic perturbations, we provide evidence that patterned death locally and globally affects adult wing shape and size. This study describes a roadmap for precise assessment of the contribution of cell death to tissue shape and outlines an important instructive role of cell death in modulating quantitatively local growth and morphogenesis of a fast-growing tissue.

Cell layer-specific expression of the homeotic MADS-box transcription factor PhDEF contributes to modular petal morphogenesis in petunia.

Floral homeotic MADS-box transcription factors ensure the correct morphogenesis of floral organs, which are organized in different cell layers deriving from distinct meristematic layers. How cells from these distinct layers acquire their respective identities and coordinate their growth to ensure normal floral organ morphogenesis is unresolved. Here, we studied petunia (Petunia × hybrida) petals that form a limb and tube through congenital fusion. We identified petunia mutants (periclinal chimeras) expressing the B-class MADS-box gene DEFICIENS in the petal epidermis or in the petal mesophyll, called wico and star, respectively. Strikingly, wico flowers form a strongly reduced tube while their limbs are almost normal, while star flowers form a normal tube but greatly reduced and unpigmented limbs, showing that petunia petal morphogenesis is highly modular. These mutants highlight the layer-specific roles of PhDEF during petal development. We explored the link between PhDEF and petal pigmentation, a well-characterized limb epidermal trait. The anthocyanin biosynthesis pathway was strongly downregulated in star petals, including its major regulator ANTHOCYANIN2 (AN2). We established that PhDEF directly binds to the AN2 terminator in vitro and in vivo, suggesting that PhDEF might regulate AN2 expression and therefore petal epidermis pigmentation. Altogether, we show that cell layer-specific homeotic activity in petunia petals differently impacts tube and limb development, revealing the relative importance of the different cell layers in the modular architecture of petunia petals.