Mammary Gland Branching Morphogenesis Research Articles

A Guide Toward Multi-scale and Quantitative Branching Analysis in the Mammary Gland.

The mammary gland consists of a bilayered epithelial structure with an extensively branched morphology. The majority of this epithelial tree is laid down during puberty, during which actively proliferating terminal end buds repeatedly elongate and bifurcate to form the basic structure of the ductal tree. Mammary ducts consist of a basal and luminal cell layer with a multitude of identified sub-lineages within both layers. The understanding of how these different cell lineages are cooperatively driving branching morphogenesis is a problem of crossing multiple scales, as this requires information on the macroscopic branched structure of the gland, as well as data on single-cell dynamics driving the morphogenic program. Here we describe a method to combine genetic lineage tracing with whole-gland branching analysis. Quantitative data on the global organ structure can be used to derive a model for mammary gland branching morphogenesis and provide a backbone on which the dynamics of individual cell lineages can be simulated and compared to lineage-tracing approaches. Eventually, these quantitative models and experiments allow to understand the couplings between the macroscopic shape of the mammary gland and the underlying single-cell dynamics driving branching morphogenesis.

Multiple roles for Bcl-3 in mammary gland branching, stromal collagen invasion, involution and tumor pathology

BackgroundThe Bcl-3 protein is an atypical member of the inhibitor of -κB family that has dual roles as a transcriptional repressor and a coactivator for dimers of NF-κB p50 and p52. Bcl-3 is expressed in mammary adenocarcinomas and can promote tumorigenesis and survival signaling and has a key role in tumor metastasis. In this study, we have investigated the role of Bcl-3 in the normal mammary gland and impact on tumor pathology.MethodsWe utilized bcl-3−/− mice to study mammary gland structure in virgins and during gestation, lactation and early involution. Expression of involution-associated genes and proteins and putative Bcl-3 target genes was examined by qRT-PCR and immunoblot analysis. Cell autonomous branching morphogenesis and collagen I invasion properties of bcl-3−/− organoids were tested in 3D hydrogel cultures. The role of Bcl-3 in tumorigenesis and tumor pathology was also assessed using a stochastic carcinogen-induced mammary tumor model.ResultsBcl-3−/− mammary glands demonstrated reduced branching complexity in virgin and pregnant mice. This defect was recapitulated in vitro where significant defects in bud formation were observed in bcl-3−/− mammary organoid cultures. Bcl-3−/− organoids showed a striking defect in protrusive collective fibrillary collagen I invasion associated with reduced expression of Fzd1 and Twist2. Virgin and pregnant bcl-3−/− glands showed increased apoptosis and rapid increases in lysosomal cell death and apoptosis after forced weaning compared to WT mice. Bcl-2 and Id3 are strongly induced in WT but not bcl-3−/− glands in early involution. Tumors in WT mice were predominately adenocarcinomas with NF-κB activation, while bcl-3−/− lesions were largely squamous lacking NF-κB and with low Bcl-2 expression.ConclusionsCollectively, our results demonstrate that Bcl-3 has a key function in mammary gland branching morphogenesis, in part by regulation of genes involved in extracellular matrix invasion. Markedly reduced levels of pro-survival proteins expression in bcl-3 null compared to WT glands 24 h post-weaning indicate that Bcl-3 has a role in moderating the rate of early phase involution. Lastly, a reduced incidence of bcl-3−/− mammary adenocarcinomas versus squamous lesions indicates that Bcl-3 supports the progression of epithelial but not metaplastic cancers.

Protocol for Studying Embryonic Mammary Gland Branching Morphogenesis Ex Vivo.

Mammary gland development starts during embryogenesis, and the process continues after birth. During development, the mammary gland undergoes massive morphological and physiological alterations including growth, invasion, and branching morphogenesis providing an ideal model for stem cell and cancer biology studies. Great efforts have been made in understanding mammary gland development during puberty and adulthood; however, the process during embryogenesis is still elusive. One reason is that the tools to study tissue dynamics during development are limited, which is partially due to the lack of an ex vivo culture method. Here we describe an updated organ culture protocol of the murine embryonic mammary gland. This powerful tool allows monitoring of growth and branching morphogenesis of mammary gland ex vivo by live imaging. In addition, we introduce a novel method for culturing intact, stroma-free mammary rudiments from late gestation mouse embryos in 3D in Matrigel. This approach can be used to identify the direct stromal cues for branching morphogenesis.

The transcription factor Sox10 is an essential determinant of branching morphogenesis and involution in the mouse mammary gland

The high mobility group-domain containing transcription factor Sox10 is an essential regulator of developmental processes and homeostasis in the neural crest, several neural crest-derived lineages and myelinating glia. Recent studies have also implicated Sox10 as an important factor in mammary stem and precursor cells. Here we employ a series of mouse mutants with constitutive and conditional Sox10 deficiencies to show that Sox10 has multiple functions in the developing mammary gland. While there is no indication for a requirement of Sox10 in the specification of the mammary placode or descending mammary bud, it is essential for both the prenatal hormone-independent as well as the pubertal hormone-dependent branching of the mammary epithelium and for proper alveologenesis during pregnancy. It furthermore acts in a dosage-dependent manner. Sox10 also plays a role during the involution process at the end of the lactation period. Whereas its effect on epithelial branching and alveologenesis are likely causally related to its function in mammary stem and precursor cells, this is not the case for its function during involution where Sox10 seems to work at least in part through regulation of the miR-424(322)/503 cluster.

Abstract 5184: N-CoR2 epigenetically regulates adhesion-dependent branching morphogenesis and its loss correlates with malignant progression of the mammary gland

Abstract Branching morphogenesis depends upon extracellular matrix remodeling and dynamic cell-cell and cell-extracellular matrix interactions, although how this developmental program is coordinated remains unclear. Expression profiling of mammary tissue at different stages of branching morphogenesis combined with histological and functional studies identified the histone deacetylase (HDAC) regulator, nuclear receptor co-repressor 2 (N-CoR2), as a key orchestrator of adhesion-dependent branching morphogenesis in the mammary gland. Molecular studies revealed that N-CoR2 actively represses fibronectin (FN) and thrombospondin-1 (THBS-1) promoter accessibility in an HDAC-dependent manner. Consistently, histological analysis revealed that N-CoR2 is abundantly expressed in the luminal epithelial cells localized along the outer rim of invading terminal end buds (TEB) and that its expression is lowest at sites of TEB branching bifurcation. In addition, reducing N-CoR2 expression drove the scattering of human mammary acini in vitro and promoted promiscuous branching of humanized mouse mammary tissue in vivo, suggesting that downregulation of N-CoR2 is required for branching morphogenesis in mammary gland. Furthermore, downregulation of N-CoR2 in breast cancer cells increased the expression of both FN and THBS-1 and reduced the expression of E-cadherin and enhanced the migratory capability of breast cancer cells. Intriguingly, immunohistochemical and clinical correlative analysis indicated that loss of N-CoR2 correlates with perturbed cell-cell and enhanced cell-ECM adhesion and acquisition of an invasive phenotype and local lymph node metastasis in a subset of breast carcinomas. These findings suggest N-CoR2 may epigenetically regulate adhesion-dependent branching morphogenesis in the mammary gland and that its loss could promote malignant progression of some breast cancers. Citation Format: Shenq-Shyang Huang, Valerie M. Weaver, Kelvin K.C. Tsai. N-CoR2 epigenetically regulates adhesion-dependent branching morphogenesis and its loss correlates with malignant progression of the mammary gland [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 5184.

DC-SCRIPT deficiency delays mouse mammary gland development and branching morphogenesis

Mammary glands are unique organs in which major adaptive changes occur in morphogenesis and development after birth. Breast cancer is the most common cancer and a major cause of mortality in females worldwide. We have previously identified the loss of expression of the transcription regulator DC-SCRIPT (Zfp366) as a prominent prognostic event in estrogen receptor positive breast cancer patients. DC-SCRIPT affects multiple transcriptional events in breast cancer cells, including estrogen and progesterone receptor-mediated transcription, and promotes CDKN2B-related cell cycle arrest. As loss of DC-SCRIPT expression appears an early event in breast cancer development, we here investigated the role of DC-SCRIPT in mammary gland development using wild-type and DC-SCRIPT knockout mice. Mice lacking DC-SCRIPT exhibited severe breeding problems and showed significant growth delay relative to littermate wild-type mice. Subsequent analysis revealed that DC-SCRIPT was expressed in mouse mammary epithelium and that DC-SCRIPT deficiency delayed mammary gland morphogenesis in vivo. Finally, analysis of 3D mammary gland organoid cultures confirmed that loss of DC-SCRIPT dramatically delayed mammary organoid branching in vitro. The study shows for the first time that DC-SCRIPT deficiency delays mammary gland morphogenesis in vivo and in vitro. These data define DC-SCRIPT as a novel modulator of mammary gland development.

Transcription factor compensation during mammary gland development in E2F knockout mice.

The E2F transcription factors control key elements of development, including mammary gland branching morphogenesis, with several E2Fs playing essential roles. Additional prior data has demonstrated that loss of individual E2Fs can be compensated by other E2F family members, but this has not been tested in a mammary gland developmental context. Here we have explored the role of the E2Fs and their ability to functionally compensate for each other during mammary gland development. Using gene expression from terminal end buds and chromatin immunoprecipitation data for E2F1, E2F2 and E2F3, we noted both overlapping and unique mammary development genes regulated by each of the E2Fs. Based on our computational findings and the fact that E2Fs share a common binding motif, we hypothesized that E2F transcription factors would compensate for each other during mammary development and function. To test this hypothesis, we generated RNA from E2F1-/-, E2F2-/- and E2F3+/- mouse mammary glands. QRT-PCR on mammary glands during pregnancy demonstrated increases in E2F2 and E2F3a in the E2F1-/- mice and an increase in E2F2 levels in E2F3+/- mice. During lactation we noted that E2F3b transcript levels were increased in the E2F2-/- mice. Given that E2Fs have previously been noted to have the most striking effects on development during puberty, we hypothesized that loss of individual E2Fs would be compensated for at that time. Double mutant mice were generated and compared with the single knockouts. Loss of both E2F1 and E2F2 revealed a more striking phenotype than either knockout alone, indicating that E2F2 was compensating for E2F1 loss. Interestingly, while E2F2 was not able to functionally compensate for E2F3+/- during mammary outgrowth, increased E2F2 expression was observed in E2F3+/- mammary glands during pregnancy day 14.5 and lactation day 5. Together, these findings illustrate the specificity of E2F family members to compensate during development of the mammary gland.

A dioxin-like compound induces hyperplasia and branching morphogenesis in mouse mammary gland, through alterations in TGF-β1 and aryl hydrocarbon receptor signaling

Hexachlorobenzene (HCB) is a widespread environmental pollutant and a dioxin-like compound that binds weakly to the aryl hydrocarbon receptor (AhR). Because AhR and transforming growth factor β1 (TGF-β1) converge to regulate common signaling pathways, alterations in this crosstalk might contribute to developing preneoplastic lesions. The aim of this study was to evaluate HCB action on TGF-β1 and AhR signaling in mouse mammary gland, through AhR+/+ and AhR−/− models. Results showed a differential effect in mouse mammary epithelial cells (NMuMG), depending on the dose: 0.05μM HCB induced cell migration and TGF-β1 signaling, whereas 5μM HCB reduced cell migration, promoted cell cycle arrest and stimulated the dioxin response element (DRE) -dependent pathway. HCB (5μM) enhanced α-smooth muscle actin expression and decreased TGF-β receptor II mRNA levels in immortalized mouse mammary fibroblasts AhR+/+, resembling the phenotype of transformed cells. Accordingly, their conditioned medium was able to enhance NMuMG cell migration. Assays in C57/Bl6 mice showed HCB (3mg/kg body weight) to enhance ductal hyperplasia, cell proliferation, estrogen receptor α nuclear localization, branch density, and the number of terminal end buds in mammary gland from AhR+/+ mice. Primary culture of mammary epithelial cells from AhR+/+ mice showed reduced AhR mRNA levels after HCB exposure (0.05 and 5μM). Interestingly, AhR−/− mice exhibited an increase in ductal hyperplasia and mammary growth in the absence of HCB treatment, thus revealing the importance of AhR in mammary development. Our findings show that environmental HCB concentrations modulate AhR and TGF-β1 signaling, which could contribute to altered mammary branching morphogenesis, likely leading to preneoplastic lesions and retaining terminal end buds.

Targeting heparanase to the mammary epithelium enhances mammary gland development and promotes tumor growth and metastasis

Heparanase is an endoglucuronidase that uniquely cleaves the heparan sulfate side chains of heparan sulfate proteoglycans. This activity ultimately alters the structural integrity of the ECM and basement membrane that becomes more prone to cellular invasion by metastatic cancer cells and cells of the immune system. In addition, enzymatically inactive heparanase was found to facilitate the proliferation and survival of cancer cells by activation of signaling molecules such as Akt, Src, signal transducer and activation of transcription (Stat), and epidermal growth factor receptor. This function is thought to be executed by the C-terminal domain of heparanase (8c), because over expression of this domain in cancer cells accelerated signaling cascades and tumor growth. We have used the regulatory elements of the mouse mammary tumor virus (MMTV) to direct the expression heparanase and the C-domain (8c) to the mammary gland epithelium of transgenic mice. Here, we report that mammary gland branching morphogenesis is increased in MMTV-heparanase and MMTV-8c mice, associating with increased Akt, Stat5 and Src phosphorylation. Furthermore, we found that the growth of tumors generated by mouse breast cancer cells and the resulting lung metastases are enhanced in MMTV-heparanase mice, thus supporting the notion that heparanase contributed by the tumor microenvironment (i.e., normal mammary epithelium) plays a decisive role in tumorigenesis. Remarkably, MMTV-8c mice develop spontaneous tumors in their mammary and salivary glands. Although this occurs at low rates and requires long latency, it demonstrates decisively the pro-tumorigenic capacity of heparanase signaling.

Un(MaSC)ing Stem Cell Dynamics in Mammary Branching Morphogenesis

The properties of stem cells that participate in mammary gland branching morphogenesis remain contested. Reporting in Nature, Scheele etal. (2017) establish a model for post-pubertal mammary branching morphogenesis in which position-dependent, lineage-restricted stem cells undergo cell mixing in order to contribute to long-term growth.

SHARPIN regulates collagen architecture and ductal outgrowth in the developing mouse mammary gland.

SHARPIN is a widely expressed multifunctional protein implicated in cancer, inflammation, linear ubiquitination and integrin activity inhibition; however, its contribution to epithelial homeostasis remains poorly understood. Here, we examined the role of SHARPIN in mammary gland development, a process strongly regulated by epithelial-stromal interactions. Mice lacking SHARPIN expression in all cells (Sharpincpdm), and mice with a stromal (S100a4-Cre) deletion of Sharpin, have reduced mammary ductal outgrowth during puberty. In contrast, Sharpincpdm mammary epithelial cells transplanted invivo into wild-type stroma, fully repopulate the mammary gland fat pad, undergo unperturbed ductal outgrowth and terminal differentiation. Thus, SHARPIN is required in mammary gland stroma during development. Accordingly, stroma adjacent to invading mammary ducts of Sharpincpdm mice displayed reduced collagen arrangement and extracellular matrix (ECM) stiffness. Moreover, Sharpincpdm mammary gland stromal fibroblasts demonstrated defects in collagen fibre assembly, collagen contraction and degradation invitro Together, these data imply that SHARPIN regulates the normal invasive mammary gland branching morphogenesis in an epithelial cell extrinsic manner by controlling the organisation of the stromal ECM.

Organochlorine pesticide modulates TGF-β1/aryl hydrocarbon receptor crosstalk in vitro, and enhances hyperplasia and branching morphogenesis in mouse mammary gland

Organochlorine pesticide modulates TGF-β1/aryl hydrocarbon receptor crosstalk in vitro, and enhances hyperplasia and branching morphogenesis in mouse mammary gland

RORα binds to E2F1 to inhibit cell proliferation and regulate mammary gland branching morphogenesis.

Retinoic acid receptor-related orphan nuclear receptor alpha (RORα) is a potent tumor suppressor that reduces cell proliferation and inhibits tumor growth. However, the molecular mechanism by which it inhibits cell proliferation remains unknown. We demonstrate a noncanonical nuclear receptor pathway in which RORα binds to E2F1 to inhibit cell cycle progression. We showed that RORα bound to the heptad repeat and marked box region of E2F1 and suppressed E2F1-regulated transcription in epithelial cells. Binding of RORα inhibited E2F1 acetylation and its DNA-binding activity by recruiting histone deacetylase 1 (HDAC1) to the protein complexes. Knockdown of HDAC1 or inhibition of HDAC activity at least partially rescued transcription factor activity of E2F1 that was repressed by RORα. Importantly, RORα levels were increased in mammary ducts compared to terminal end buds and inversely correlated with expression of E2F1 target genes and cell proliferation. Silencing RORα in mammary epithelial cells significantly enhanced cell proliferation in the ductal epithelial cells and promoted side branching of the mammary ducts. These results reveal a novel link between RORα and E2F1 in regulating cell cycle progression and mammary tissue morphogenesis.

The co-factor of LIM domains (CLIM/LDB/NLI) maintains basal mammary epithelial stem cells and promotes breast tumorigenesis.

Mammary gland branching morphogenesis and ductal homeostasis relies on mammary stem cell function for the maintenance of basal and luminal cell compartments. The mechanisms of transcriptional regulation of the basal cell compartment are currently unknown. We explored these mechanisms in the basal cell compartment and identified the Co-factor of LIM domains (CLIM/LDB/NLI) as a transcriptional regulator that maintains these cells. Clims act within the basal cell compartment to promote branching morphogenesis by maintaining the number and proliferative potential of basal mammary epithelial stem cells. Clim2, in a complex with LMO4, supports mammary stem cells by directly targeting the Fgfr2 promoter in basal cells to increase its expression. Strikingly, Clims also coordinate basal-specific transcriptional programs to preserve luminal cell identity. These basal-derived cues inhibit epidermis-like differentiation of the luminal cell compartment and enhance the expression of luminal cell-specific oncogenes ErbB2 and ErbB3. Consistently, basal-expressed Clims promote the initiation and progression of breast cancer in the MMTV-PyMT tumor model, and the Clim-regulated branching morphogenesis gene network is a prognostic indicator of poor breast cancer outcome in humans.

Pttg1/securin is required for the branching morphogenesis of the mammary gland and suppresses mammary tumorigenesis

Pituitary tumor transforming gene 1 (Pttg1) encodes the mammalian securin, which is an inhibitor of separase (a protease required for the separation of sister chromatids in mitosis and meiosis). PTTG1 is overexpressed in a number of human cancers and has been suggested to be an oncogene. However, we found that, in Pttg1-mutant females, the mammary epithelial cells showed increased proliferation and precocious branching morphogenesis. In accord with these phenotypic changes, progesterone receptor, cyclin D1, and Mmp2 were up-regulated whereas p21 (Cdkn1a) was down-regulated. These molecular changes provide explanation for the observed developmental defects, and suggest that Pttg1 is a tumor suppressor. Indeed, mice lacking Pttg1 developed spontaneous mammary tumors. Furthermore, in human breast tumors, PTTG1 protein levels were down-regulated and the reduction was significantly correlated with the tumor grade.

Cdc42 overexpression induces hyperbranching in the developing mammary gland by enhancing cell migration

IntroductionThe Rho GTPase Cdc42 is overexpressed and hyperactivated in breast tumors compared to normal breast tissue. Cdc42 regulates key processes that are critical for mammary gland morphogenesis and become disrupted during the development, progression, and metastasis of breast cancer. However, the contribution of Cdc42 to normal and neoplastic mammary gland development in vivo remains poorly understood. We were therefore interested in investigating the effects of Cdc42 overexpression on mammary gland morphogenesis as a first step toward understanding how its overexpression may contribute to mammary tumorigenesis.MethodsWe developed a tetracycline-regulatable Cdc42 overexpression mouse model in which Cdc42 can be inducibly overexpressed in the developing mammary gland. The effects of Cdc42 overexpression during postnatal mammary gland development were investigated using in vivo and in vitro approaches, including morphometric analysis of wholemounted mammary glands, quantification of histological markers, and primary mammary epithelial cell (MEC) functional and biochemical assays.ResultsAnalysis of Cdc42-overexpressing mammary glands revealed abnormal terminal end bud (TEB) morphologies, characterized by hyperbudding and trifurcation, and increased side branching within the ductal tree. Quantification of markers of proliferation and apoptosis suggested that these phenotypes were not due to increased cell proliferation or survival. Rather, Cdc42 overexpressing MECs were more migratory and contractile and formed dysmorphic, invasive acini in three-dimensional cultures. Cdc42 and RhoA activities, phosphorylated myosin light chain, and MAPK signaling, which contribute to migration and invasion, were markedly elevated in Cdc42 overexpressing MECs. Interestingly, Cdc42 overexpressing mammary glands displayed several features associated with altered epithelial-stromal interactions, which are known to regulate branching morphogenesis. These included increased stromal thickness and collagen deposition, and stromal cells isolated from Cdc42 overexpressing mammary glands exhibited elevated mRNA expression of extracellular matrix proteins and remodeling enzymes.ConclusionsThese data suggest that Cdc42 overexpression disrupts mammary gland branching morphogenesis by altering Rho GTPase and MAPK signaling, leading to increased MEC contractility and migration in association with stromal alterations. Our studies provide insight into how aberrant Cdc42 expression may contribute to mammary tumorigenesis.

The control of branching morphogenesis.

Many organs of higher organisms are heavily branched structures and arise by an apparently similar process of branching morphogenesis. Yet the regulatory components and local interactions that have been identified differ greatly in these organs. It is an open question whether the regulatory processes work according to a common principle and how far physical and geometrical constraints determine the branching process. Here, we review the known regulatory factors and physical constraints in lung, kidney, pancreas, prostate, mammary gland and salivary gland branching morphogenesis, and describe the models that have been formulated to analyse their impacts.

A Rab10-Dependent Mechanism for Polarized Basement Membrane Secretion during Organ Morphogenesis

Basement membranes (BMs) are specialized extracellular matrices that are essential for epithelial structure and morphogenesis. However, little is known about how BM proteins are delivered to the basal cell surface or how this process is regulated during development. Here, we identify a mechanism for polarized BM secretion in the Drosophila follicle cells. BM proteins are synthesized in a basal endoplasmic reticulum (ER) compartment from localized mRNAs and are then exported through Tango1-positive ER exit sites to basal Golgi clusters. Next, Crag targets Rab10 to structures in the basal cytoplasm, where it restricts protein delivery to the basal surface. These events occur during egg chamber elongation, a morphogenetic process that depends on follicle cell planar polarity and BM remodeling. Significantly, Tango1 and Rab10 are also planar polarized at the basal epithelial surface. We propose that the spatial control of BM production along two tissue axes promotes exocytic efficiency, BM remodeling, and organ morphogenesis.

Expression Pattern of the Pro-apoptotic Gene PAR-4 During the Morphogenesis of MCF-10A Human Mammary Epithelial Cells

The histological organization of the mammary gland involves a spatial interaction of epithelial and myoepithelial cells with the specialized basement membrane (BM), composed of extra-cellular matrix (ECM) proteins, which is disrupted during the tumorigenic process. The interactions between mammary epithelial cells and ECM components play a major role in mammary gland branching morphogenesis. Critical signals for mammary epithelial cell proliferation, differentiation, and survival are provided by the ECM proteins. Three-dimensional (3D) cell culture was developed to establish a system that simulates several features of the breast epithelium in vivo; 3D cell culture of the spontaneously immortalized cell line, MCF10A, is a well-established model system to study breast epithelial cell biology and morphogenesis. Mammary epithelial cells grown in 3D form spheroids, acquire apicobasal polarization, and form lumens that resemble acini structures, processes that involve cell death. Using this system, we evaluated the expression of the pro-apoptotic gene PAWR (PKC apoptosis WT1 regulator; also named PAR-4, prostate apoptosis response-4) by immunofluorescence and quantitative real time PCR (qPCR). A time-dependent increase in PAR-4 mRNA expression was found during the process of MCF10A acinar morphogenesis. Confocal microscopy analysis also showed that PAR-4 protein was highly expressed in the MCF10A cells inside the acini structure. During the morphogenesis of MCF10A cells in 3D cell culture, the cells within the lumen showed caspase-3 activation, indicating apoptotic activity. PAR-4 was only partially co-expressed with activated caspase-3 on these cells. Our results provide evidence, for the first time, that PAR-4 is differentially expressed during the process of MCF10A acinar morphogenesis.

Cadherin-11 functions during mammary gland branching morphogenesis

Cadherin-11 functions during mammary gland branching morphogenesis