The UCSF Mouse Inventory Database Application, an Open Source Web App for Sharing Mutant Mice Within a Research Community.

The UCSF Mouse Inventory Database Application is an open-source Web App that provides information about the mutant alleles, transgenes, and inbred strains maintained by investigators at the university and facilitates sharing of these resources within the university community. The Application is designed to promote collaboration, decrease the costs associated with obtaining genetically-modified mice, and increase access to mouse lines that are difficult to obtain. An inventory of the genetically-modified mice on campus and the investigators who maintain them is compiled from records of purchases from external sources, transfers from researchers within and outside the university, and from data provided by users. These data are verified and augmented with relevant information harvested from public databases, and stored in a succinct, searchable database secured on the university network. Here we describe this resource and provide information about how to implement and maintain such a mouse inventory database application at other institutions.

What comes after a career in science? Asian Art!

Throughout my professional life I was an academic researcher who had a successful career in developmental biology. When I started the process of closing my lab and retiring from the University, I came to the conclusion that I wanted to do something completely different. I trained for and became a docent at the Asian Art Museum of San Francisco. I apply the skills I acquired as a scientist on a regular basis in my work at the museum, often in unexpected ways. In this article, I discuss my journey from the lab to the museum and what I learned along the way.

MicroRNA-145 Promotes the Phenotype of Human Glioblastoma Cells Selected for Invasion.

BACKGROUND/AIM: Novel treatment strategies aiming to eliminate or attenuate the invasive phenotype of glioblastoma multiforme (GBM), the most common and aggressive primary brain tumor, could offer a profound therapeutic benefit to patients. We previously demonstrated one method to create invasive sub-populations of GBM cells (IM3 cells) and a positive regulatory role for the miR-143/-145 locus in enhancing the invasion of GBM cells. Herein, we investigated the correlation between miR-145 and srGAP1 (SLIT-ROBO Rho GTPase-activating protein1) that is purported to be a target of miR-145 and involved in migration and invasion. MATERIALS AND METHODS: IM3 cells were created by a serial selection by using Boyden chambers®. Antisense-miR-145 was transfected into IM3 cells by using lipofectamine 2000. Quantitative reverse transcription polymerase chain reaction (RT-qPCR) and western blot were employed to analyze the expression of srGAP1. RESULTS: The invasiveness of U87-IM3 and U251-IM3 is attenuated by transfection of antisense miR-145. In addition, srGAP1 was down-regulated in U87-IM3 and U251-IM3 cells compared to parental cells. CONCLUSION: The elevated miR-145 present in invasive glioblastoma cells (IM3 cells) targets and down-regulated srGAP1, thereby allowing downstream G-proteins to remain in their active state and promote the observed invasive phenotype.

Menin determines K-RAS proliferative outputs in endocrine cells.

Endocrine cell proliferation fluctuates dramatically in response to signals that communicate hormone demand. The genetic alterations that override these controls in endocrine tumors often are not associated with oncogenes common to other tumor types, suggesting that unique pathways govern endocrine proliferation. Within the pancreas, for example, activating mutations of the prototypical oncogene KRAS drive proliferation in all pancreatic ductal adenocarcimomas but are never found in pancreatic endocrine tumors. Therefore, we asked how cellular context impacts K-RAS signaling. We found that K-RAS paradoxically suppressed, rather than promoted, growth in pancreatic endocrine cells. Inhibition of proliferation by K-RAS depended on antiproliferative RAS effector RASSF1A and blockade of the RAS-activated proproliferative RAF/MAPK pathway by tumor suppressor menin. Consistent with this model, a glucagon-like peptide 1 (GLP1) agonist, which stimulates ERK1/2 phosphorylation, did not affect endocrine cell proliferation by itself, but synergistically enhanced proliferation when combined with a menin inhibitor. In contrast, inhibition of MAPK signaling created a synthetic lethal interaction in the setting of menin loss. These insights suggest potential strategies both for regenerating pancreatic ß cells for people with diabetes and for targeting menin-sensitive endocrine tumors.

Sprouty genes function in suppression of prostate tumorigenesis.

Expression of Sprouty genes is frequently decreased or absent in human prostate cancer, implicating them as suppressors of tumorigenesis. Here we show they function in prostate tumor suppression in the mouse. Concomitant inactivation of Spry1 and Spry2 in prostate epithelium causes ductal hyperplasia and low-grade prostatic intraepithelial neoplasia (PIN). However, when Spry1 and Spry2 loss-of-function occurs in the context of heterozygosity for a null allele of the tumor suppressor gene Pten, there is a striking increase in PIN and evidence of neoplastic invasion. Conversely, expression of a Spry2 gain-of-function transgene in Pten null prostatic epithelium suppresses the tumorigenic effects of loss of Pten function. We show that Sprouty gene loss-of-function results in hyperactive RAS/ERK1/2 signaling throughout the prostate epithelium and cooperates with heterozygosity for a Pten null allele to promote hyperactive PI3K/AKT signaling. Furthermore, Spry2 gain-of-function can suppress hyperactivation of AKT caused by the absence of PTEN. Together, these results point to a key genetic interaction between Sprouty genes and Pten in prostate tumorigenesis and provide strong evidence that Sprouty genes can function to modulate signaling via the RAS/ERK1/2 and PI3K/AKT pathways. The finding that Sprouty genes suppress tumorigenesis caused by Pten loss-of-function suggests that therapeutic approaches aimed at restoring normal feedback mechanisms triggered by receptor tyrosine kinase signaling, including Sprouty gene expression, may provide an effective strategy to delay or prevent high-grade PIN and invasive prostate cancer.

Serial selection for invasiveness increases expression of miR-143/miR-145 in glioblastoma cell lines.

BACKGROUND: Glioblastoma multiforme (GBM) is the most common primary central nervous system malignancy and its unique invasiveness renders it difficult to treat. This invasive phenotype, like other cellular processes, may be controlled in part by microRNAs - a class of small non-coding RNAs that act by altering the expression of targeted messenger RNAs. In this report, we demonstrate a straightforward method for creating invasive subpopulations of glioblastoma cells (IM3 cells). To understand the correlation between the expression of miRNAs and the invasion, we fully profiled 1263 miRNAs on six different cell lines and two miRNAs, miR-143 and miR-145, were selected for validation of their biological properties contributing to invasion. Further, we investigated an ensemble effect of both miR-143 and miR-145 in promoting invasion. METHODS: By repeated serial invasion through Matrigel®-coated membranes, we isolated highly invasive subpopulations of glioma cell lines. Phenotypic characterization of these cells included in vitro assays for proliferation, attachment, and invasion. Micro-RNA expression was compared using miRCURY arrays (Exiqon). In situ hybridization allowed visualization of the regional expression of miR-143 and miR-145 in tumor samples, and antisense probes were used investigate in vitro phenotypic changes seen with knockdown in their expression. RESULTS: The phenotype we created in these selected cells proved stable over multiple passages, and their microRNA expression profiles were measurably different. We found that two specific microRNAs expressed from the same genetic locus, miR-143 and miR-145, were over-expressed in our invasive subpopulations. Further, we also found that combinatorial treatment of these cells with both antisense-miRNAs (antimiR-143 and -145) will abrogated their invasion without decreasing cell attachment or proliferation. CONCLUSIONS: To best of our knowledge, these data demonstrate for the first time that miR-143 and miR-145 regulate the invasion of glioblastoma and that miR-143 and -145 could be potential therapeutic target for anti-invasion therapies of glioblastoma patients.

Control of mitotic spindle angle by the RAS-regulated ERK1/2 pathway determines lung tube shape.

During early lung development, airway tubes change shape. Tube length increases more than circumference as a large proportion of lung epithelial cells divide parallel to the airway longitudinal axis. We show that this bias is lost in mutants with increased extracellular signal-regulated kinase 1 (ERK1) and ERK2 activity, revealing a link between the ERK1/2 signaling pathway and the control of mitotic spindle orientation. Using a mathematical model, we demonstrate that change in airway shape can occur as a function of spindle angle distribution determined by ERK1/2 signaling, independent of effects on cell proliferation or cell size and shape. We identify sprouty genes, which encode negative regulators of fibroblast growth factor 10 (FGF10)-mediated RAS-regulated ERK1/2 signaling, as essential for controlling airway shape change during development through an effect on mitotic spindle orientation.

SHH propagates distal limb bud development by enhancing CYP26B1-mediated retinoic acid clearance via AER-FGF signalling.

The essential roles of SHH in anteroposterior (AP) and AER-FGF signalling in proximodistal (PD) limb bud development are well understood. In addition, these morphoregulatory signals are key components of the self-regulatory SHH/GREM1/AER-FGF feedback signalling system that regulates distal progression of limb bud development. This study uncovers an additional signalling module required for coordinated progression of limb bud axis development. Transcriptome analysis using Shh-deficient mouse limb buds revealed that the expression of proximal genes was distally extended from early stages onwards, which pointed to a more prominent involvement of SHH in PD limb axis development. In particular, retinoic acid (RA) target genes were upregulated proximally, while the expression of the RA-inactivating Cyp26b1 enzyme was downregulated distally, pointing to increased RA activity in Shh-deficient mouse limb buds. Further genetic and molecular analysis established that Cyp26b1 expression is regulated by AER-FGF signalling. During initiation of limb bud outgrowth, the activation of Cyp26b1 expression creates a distal 'RA-free' domain, as indicated by complementary downregulation of a transcriptional sensor of RA activity. Subsequently, Cyp26b1 expression increases as a consequence of SHH-dependent upregulation of AER-FGF signalling. To better understand the underlying signalling interactions, computational simulations of the spatiotemporal expression patterns and interactions were generated. These simulations predicted the existence of an antagonistic AER-FGF/CYP26B1/RA signalling module, which was verified experimentally. In summary, SHH promotes distal progression of limb development by enhancing CYP26B1-mediated RA clearance as part of a signalling network linking the SHH/GREM1/AER-FGF feedback loop to the newly identified AER-FGF/CYP26B1/RA module.

Hedgehog signaling regulates the generation of ameloblast progenitors in the continuously growing mouse incisor.

In many organ systems such as the skin, gastrointestinal tract and hematopoietic system, homeostasis is dependent on the continuous generation of differentiated progeny from stem cells. The rodent incisor, unlike human teeth, grows throughout the life of the animal and provides a prime example of an organ that rapidly deteriorates if newly differentiated cells cease to form from adult stem cells. Hedgehog (Hh) signaling has been proposed to regulate self-renewal, survival, proliferation and/or differentiation of stem cells in several systems, but to date there is little evidence supporting a role for Hh signaling in adult stem cells. We used in vivo genetic lineage tracing to identify Hh-responsive stem cells in the mouse incisor and we show that sonic hedgehog (SHH), which is produced by the differentiating progeny of the stem cells, signals to several regions of the incisor. Using a hedgehog pathway inhibitor (HPI), we demonstrate that Hh signaling is not required for stem cell survival but is essential for the generation of ameloblasts, one of the major differentiated cell types in the tooth, from the stem cells. These results therefore reveal the existence of a positive-feedback loop in which differentiating progeny produce the signal that in turn allows them to be generated from stem cells.

Cardiac conduction is required to preserve cardiac chamber morphology.

Electrical cardiac forces have been previously hypothesized to play a significant role in cardiac morphogenesis and remodeling. In response to electrical forces, cultured cardiomyocytes rearrange their cytoskeletal structure and modify their gene expression profile. To translate such in vitro data to the intact heart, we used a collection of zebrafish cardiac mutants and transgenics to investigate whether cardiac conduction could influence in vivo cardiac morphogenesis independent of contractile forces. We show that the cardiac mutant dco(s226) develops heart failure and interrupted cardiac morphogenesis following uncoordinated ventricular contraction. Using in vivo optical mapping/calcium imaging, we determined that the dco cardiac phenotype was primarily due to aberrant ventricular conduction. Because cardiac contraction and intracardiac hemodynamic forces can also influence cardiac development, we further analyzed the dco phenotype in noncontractile hearts and observed that disorganized ventricular conduction could affect cardiomyocyte morphology and subsequent heart morphogenesis in the absence of contraction or flow. By positional cloning, we found that dco encodes Gja3/Cx46, a gap junction protein not previously implicated in heart formation or function. Detailed analysis of the mouse Cx46 mutant revealed the presence of cardiac conduction defects frequently associated with human heart failure. Overall, these in vivo studies indicate that cardiac electrical forces are required to preserve cardiac chamber morphology and may act as a key epigenetic factor in cardiac remodeling.

Kidney development in the absence of Gdnf and Spry1 requires Fgf10.

GDNF signaling through the Ret receptor tyrosine kinase (RTK) is required for ureteric bud (UB) branching morphogenesis during kidney development in mice and humans. Furthermore, many other mutant genes that cause renal agenesis exert their effects via the GDNF/RET pathway. Therefore, RET signaling is believed to play a central role in renal organogenesis. Here, we re-examine the extent to which the functions of Gdnf and Ret are unique, by seeking conditions in which a kidney can develop in their absence. We find that in the absence of the negative regulator Spry1, Gdnf, and Ret are no longer required for extensive kidney development. Gdnf-/-;Spry1-/- or Ret-/-;Spry1-/- double mutants develop large kidneys with normal ureters, highly branched collecting ducts, extensive nephrogenesis, and normal histoarchitecture. However, despite extensive branching, the UB displays alterations in branch spacing, angle, and frequency. UB branching in the absence of Gdnf and Spry1 requires Fgf10 (which normally plays a minor role), as removal of even one copy of Fgf10 in Gdnf-/-;Spry1-/- mutants causes a complete failure of ureter and kidney development. In contrast to Gdnf or Ret mutations, renal agenesis caused by concomitant lack of the transcription factors ETV4 and ETV5 is not rescued by removing Spry1, consistent with their role downstream of both RET and FGFRs. This shows that, for many aspects of renal development, the balance between positive signaling by RTKs and negative regulation of this signaling by SPRY1 is more critical than the specific role of GDNF. Other signals, including FGF10, can perform many of the functions of GDNF, when SPRY1 is absent. But GDNF/RET signaling has an apparently unique function in determining normal branching pattern. In contrast to GDNF or FGF10, Etv4 and Etv5 represent a critical node in the RTK signaling network that cannot by bypassed by reducing the negative regulation of upstream signals.

Gene and microRNA expression in p53-deficient day 8.5 mouse embryos.

BACKGROUND: Neural tube defects (NTDs) are one of the most common human birth defects, with a prevalence of approximately 1 in 1000 live births in the United States. In animal studies, deletion of p53 leads to a significant increase in embryos that exhibit exencephaly. Whereas several studies have closely investigated the morphologic changes of p53-deficient embryos, no study has reported the molecular-level alteration in p53-deficient embryos. Here we attempt to identify genes and microRNAs (miRNAs) modified by deletion of p53 in day 8.5 mouse embryos. METHODS: Mouse embryos from p53 heterozygous crosses were collected, genotyped, and embryos of similar genotype (+/+; +/-; -/-) were pooled. RNA from the pooled samples was isolated to determine mRNA and miRNA expression levels using Whole Genome Bioarrays and Low Density Arrays, respectively. RESULTS: In p53 -/- embryos, 388 genes showed statistically significant alteration in gene expression of more than twofold compared to p53 +/+ embryos. Expression of p53 and well known p53 target genes, such as p21 and cyclin G1, were significantly down-regulated in p53 -/- embryos. In contrast, expression of other p53 target genes, such as Mdm2, Noxa, and Puma, were unchanged. We also identified six genes (Csk, Itga3, Jarid2, Prkaca, Rarg, and Sall4), known to cause NTDs when deleted, that are also down-regulated in p53 -/- embryos. Finally, five miRNAs (mir-1, mir-30e-3p, mir-142-3p, mir-301, and mir-331) also showed statistically significant alterations in expression levels in p53 -/- embryos compared to p53 +/+ embryos. Combined analysis of the experimental data using stepwise regression model and two publicly available algorithms identified putative target genes of these miRNAs. CONCLUSIONS: Our data have identified genes and miRNAs that may be involved in the mechanisms underlining NTDs and begin to define the developmental role of p53 in the etiology of NTDs.

Expression of sprouty2 inhibits B-cell proliferation and is epigenetically silenced in mouse and human B-cell lymphomas.

B-cell lymphoma is the most common immune system malignancy. TCL1 transgenic mice (TCL1-tg), in which TCL1 is ectopically expressed in mature lymphocytes, develop multiple B- and T-cell leukemia and lymphoma subtypes, supporting an oncogenic role for TCL1 that probably involves AKT and MAPK-ERK signaling pathway augmentation. Additional, largely unknown genetic and epigenetic alterations cooperate with TCL1 during lymphoma progression. We examined DNA methylation patterns in TCL1-tg B-cell tumors to discover tumor-associated epigenetic changes, and identified hypermethylation of sprouty2 (Spry2). Sprouty proteins are context-dependent negative or positive regulators of MAPK-ERK pathway signaling, but their role(s) in B-cell physiology or pathology are unknown. Here we show that repression of Spry2 expression in TCL1-tg mouse and human B-cell lymphomas and cell lines is associated with dense DNA hypermethylation and was reversed by inhibition of DNA methylation. Spry2 expression was induced in normal splenic B cells by CD40/B-cell receptor costimulation and regulated a negative feedback loop that repressed MAPK-ERK signaling and decreased B-cell viability. Conversely, loss of Spry2 function hyperactivated MAPK-ERK signaling and caused increased B-cell proliferation. Combined, these results implicate epigenetic silencing of Spry2 expression in B lymphoma progression and suggest it as a companion lesion to ectopic TCL1 expression in enhancing MAPK-ERK pathway signaling.

The roles of p53 and p21 in normal development and hyperthermia-induced malformations.

BACKGROUND: Hyperthermia (HS) is a well-studied teratogen that induces serious malformations, including neural tube defects. Our previous studies have shown that HS induces apoptosis by activating the mitochondrial apoptotic pathway. Prior to activation of the mitochondrial apoptotic pathway, HS also activates p53 and its target genes. In the present study, we determine whether p53 and/or p21 play a role as teratogen suppressors or inducers of HS-induced malformations. METHODS: Pregnant mice carrying all three p53 or p21 genotype embryos were exposed to HS on day 8.5. Subsequently, fetuses were collected on day 15.5, and genotyped. In addition to genotype, we also determined the number of resorptions and dead fetuses as well as the number and types of external malformations. RESULTS: In the absence of HS exposure, fetuses exhibiting exencephaly and spina bifida were observed in approximately 11% of p53 -/- fetuses, whereas no malformations were observed among p21 -/- fetuses. Exposure to HS resulted in an increase in exencephaly and polydactyly in fetuses of all three p53 genotypes. However, the incidence of these malformations was statistically significantly higher in p53 -/- compared to p53 +/- and p53 +/+ fetuses. Exencephaly was the only malformation observed in p21 fetuses exposed to HS, with an approximately 2-fold increase among p21 +/- and a 3-fold increase among p21 -/- compared to p21 +/+ fetuses. CONCLUSIONS: Our study confirms that p53 plays a role in normal development and has shown, for the first time that p53 and p21 function to suppress HS-induced malformations.

MicroRNA-21 contributes to myocardial disease by stimulating MAP kinase signalling in fibroblasts.

MicroRNAs comprise a broad class of small non-coding RNAs that control expression of complementary target messenger RNAs. Dysregulation of microRNAs by several mechanisms has been described in various disease states including cardiac disease. Whereas previous studies of cardiac disease have focused on microRNAs that are primarily expressed in cardiomyocytes, the role of microRNAs expressed in other cell types of the heart is unclear. Here we show that microRNA-21 (miR-21, also known as Mirn21) regulates the ERK-MAP kinase signalling pathway in cardiac fibroblasts, which has impacts on global cardiac structure and function. miR-21 levels are increased selectively in fibroblasts of the failing heart, augmenting ERK-MAP kinase activity through inhibition of sprouty homologue 1 (Spry1). This mechanism regulates fibroblast survival and growth factor secretion, apparently controlling the extent of interstitial fibrosis and cardiac hypertrophy. In vivo silencing of miR-21 by a specific antagomir in a mouse pressure-overload-induced disease model reduces cardiac ERK-MAP kinase activity, inhibits interstitial fibrosis and attenuates cardiac dysfunction. These findings reveal that microRNAs can contribute to myocardial disease by an effect in cardiac fibroblasts. Our results validate miR-21 as a disease target in heart failure and establish the therapeutic efficacy of microRNA therapeutic intervention in a cardiovascular disease setting.

Genetic mosaic analysis reveals FGF receptor 2 function in terminal end buds during mammary gland branching morphogenesis.

FGF signaling is associated with breast cancer and is required for mammary placode formation in the mouse. In this study, we employed a genetic mosaic analysis based on Cre-mediated recombination to investigate FGF receptor 2 (Fgfr2) function in the postnatal mammary gland. Mosaic inactivation of Fgfr2 by the MMTV-Cre transgene enabled us to compare the behavior of Fgfr2 null and Fgfr2 heterozygous cells in the same gland. Fgfr2 null cells were at a competitive disadvantage to their Fgfr2 heterozygous neighbors in the highly proliferative terminal end buds (TEBs) at the invasion front, owing to a negative effect of loss of Fgfr2 function on cell proliferation. However, Fgfr2 null cells were tolerated in mature ducts. In these genetic mosaic mammary glands, the epithelial network is apparently built by TEBs that over time are composed of a progressively larger proportion of Fgfr2-positive cells. However, subsequently, most cells lose Fgfr2 function, presumably due to additional rounds of Cre-mediated recombination. Using an independent strategy to create mosaic mammary glands, which employed an adenovirus-Cre that acts only once, we confirmed that Fgfr2 null cells were out-competed by neighboring Fgfr2 heterozygous cells. Together, our data demonstrate that Fgfr2 functions in the proliferating and invading TEBs, but it is not required in the mature ducts of the pubertal mammary gland.

Genetic evidence that FGFs have an instructive role in limb proximal-distal patterning.

Half a century ago, the apical ectodermal ridge (AER) at the distal tip of the tetrapod limb bud was shown to produce signals necessary for development along the proximal-distal (P-D) axis, but how these signals influence limb patterning is still much debated. Fibroblast growth factor (FGF) gene family members are key AER-derived signals, with Fgf4, Fgf8, Fgf9 and Fgf17 expressed specifically in the mouse AER. Here we demonstrate that mouse limbs lacking Fgf4, Fgf9 and Fgf17 have normal skeletal pattern, indicating that Fgf8 is sufficient among AER-FGFs to sustain normal limb formation. Inactivation of Fgf8 alone causes a mild skeletal phenotype; however, when we also removed different combinations of the other AER-FGF genes, we obtained unexpected skeletal phenotypes of increasing severity, reflecting the contribution that each FGF can make to the total AER-FGF signal. Analysis of the compound mutant limb buds revealed that, in addition to sustaining cell survival, AER-FGFs regulate P-D-patterning gene expression during early limb bud development, providing genetic evidence that AER-FGFs function to specify a distal domain and challenging the long-standing hypothesis that AER-FGF signalling is permissive rather than instructive for limb patterning. We discuss how a two-signal model for P-D patterning can be integrated with the concept of early specification to explain the genetic data presented here.

Loss of apoptosis-inducing factor results in cell-type-specific neurogenesis defects.

Mitochondrial dysfunction is commonly associated with neurodegeneration in the aging brain. In addition, the importance of mitochondrial function during brain development is illustrated by the neurological deficits observed in infants with mitochondrial complex deficiencies. However, the extent to which abnormalities in mitochondrial function might impact neurogenesis during brain development is not well understood. Previously, we demonstrated that adult harlequin (Hq) mutant mice, which have an 80% reduction in the mitochondrial protein apoptosis-inducing factor (AIF), exhibited signs of oxidative stress and progressive loss of adult cerebellar and retinal neurons. To assess whether in addition to its role in postmitotic neuron survival Aif is also necessary for cerebellar development, we analyzed embryos in which Aif was deleted in the prospective midbrain and cerebellum at a very early stage of development using an En1 (engrailed 1) promoter-driven cre recombinase gene. These mutant mice, which died at birth, had midbrain defects and dramatic deficits in cerebellar Purkinje and granule cell precursors. Additional analysis revealed that Aif-null Purkinje cell precursors prematurely entered S-phase, but most failed to undergo mitosis and ultimately died via apoptosis. In contrast, proliferation of mutant granule cell precursors was blocked before S-phase. Mice in which Aif was deleted later in embryogenesis using a nestin promoter-driven cre gene survive for several days after birth, and postnatal granule cell precursors in these mice also failed to enter S-phase. Our results indicate that the loss of Aif results in cell cycle abnormalities in a neuron-specific manner during cerebellar development.

The branching programme of mouse lung development.

Mammalian lungs are branched networks containing thousands to millions of airways arrayed in intricate patterns that are crucial for respiration. How such trees are generated during development, and how the developmental patterning information is encoded, have long fascinated biologists and mathematicians. However, models have been limited by a lack of information on the normal sequence and pattern of branching events. Here we present the complete three-dimensional branching pattern and lineage of the mouse bronchial tree, reconstructed from an analysis of hundreds of developmental intermediates. The branching process is remarkably stereotyped and elegant: the tree is generated by three geometrically simple local modes of branching used in three different orders throughout the lung. We propose that each mode of branching is controlled by a genetically encoded subroutine, a series of local patterning and morphogenesis operations, which are themselves controlled by a more global master routine. We show that this hierarchical and modular programme is genetically tractable, and it is ideally suited to encoding and evolving the complex networks of the lung and other branched organs.