A developmental framework linking neurogenesis and circuit formation in the Drosophila CNS
Abstract
The mechanisms specifying neuronal diversity are well characterized, yet it remains unclear how or if these mechanisms regulate neural circuit assembly. To address this, we mapped the developmental origin of 160 interneurons from seven bilateral neural progenitors (neuroblasts) and identify them in a synapse-scale TEM reconstruction of the Drosophila larval central nervous system. We find that lineages concurrently build the sensory and motor neuropils by generating sensory and motor hemilineages in a Notch-dependent manner. Neurons in a hemilineage share common synaptic targeting within the neuropil, which is further refined based on neuronal temporal identity. Connectome analysis shows that hemilineage-temporal cohorts share common connectivity. Finally, we show that proximity alone cannot explain the observed connectivity structure, suggesting hemilineage/temporal identity confers an added layer of specificity. Thus, we demonstrate that the mechanisms specifying neuronal diversity also govern circuit formation and function, and that these principles are broadly applicable throughout the nervous system.
Introduction
Tremendous progress has been made in understanding the molecular mechanisms generating neuronal diversity in both vertebrate and invertebrate model systems. In mammals, spatial cues generate distinct pools of progenitors, which generate neuronal diversity in each spatial domain (Sagner and Briscoe, 2019). The same process occurs in invertebrates like Drosophila, but with a smaller number of cells, and this process is particularly well understood. The first step occurs when spatial patterning genes act combinatorially to establish single, unique progenitor (neuroblast) identities (Skeath and Thor, 2003). These patterning genes endow each neuroblast with a unique spatial identity (Figure 1A, left).
The second step is temporal patterning – the specification of neuronal identity based on birth-order – an evolutionarily conserved mechanism for generating neuronal diversity (Kohwi and Doe, 2013; Rossi et al., 2017). Here, we focus on Drosophila embryonic neuroblasts, which undergo a cascade of temporal transcription factors: Hunchback (Hb), Krüppel (Kr), Pdm, and Castor (Cas) (Isshiki et al., 2001). Each temporal transcription factor is inherited by ganglion mother cells (GMCs) born during each expression window. The combination of spatial and temporal factors endows each GMC with a unique identity (Figure 1A, middle).
The third step is hemilineage specification, which was initially characterized in Drosophila larval and adult neurogenesis (Lee et al., 2020; Truman et al., 2010), and may also be used in vertebrate neurogenesis (Peng et al., 2007). Hemilineages are formed by GMC asymmetric division into a pair of post-mitotic neurons; during this division, the Notch inhibitor Numb (Nb) is partitioned into one neuron (NotchOFF neuron), whereas the other sibling neuron receives active Notch signaling (NotchON neuron), thereby establishing NotchON and NotchOFF hemilineages (Figure 1A, right). In summary, three mechanisms generate neuronal diversity within the embryonic central nervous system (CNS): neuroblast spatial identity, GMC temporal identity, and neuronal hemilineage identity.
A great deal of progress has also been made in understanding neural circuit formation in both vertebrates and invertebrate model systems, revealing a multi-step mechanism. Neurons initially target their axons to broad regions (e.g., thalamus/cortex), followed by targeting to a neuropil domain (glomeruli/layer), and finally forming highly specific synapses within the targeted domain (Kolodkin and Tessier-Lavigne, 2011; Moyle et al., 2021).
Despite the progress in understanding the generation of neuronal diversity and the mechanisms governing axon guidance and neuropil targeting, how these two developmental processes are coordinated remains largely unknown. While it is accepted that the identity of a neuron is linked to its connectivity, the developmental mechanisms involved are unclear. For example, do clonally related neurons target similar regions of the neuropil due to the expression of similar guidance cues? Do temporal cohorts born at similar times show preferential connectivity? Here, we address the question of whether any of the three developmental mechanisms (spatial, temporal, hemilineage identity) are correlated with any of the three circuit-wiring mechanisms (neurite targeting, synapse localization, connectivity). We map the developmental origin for 80 bilateral pairs of interneurons in abdominal segment 1 (A1) by identifying and reconstructing these neurons within a full CNS TEM volume (Ohyama et al., 2015); this is over a quarter of the ~300 neurons per hemisegment. We make the unexpected observation that hemilineage identity determines neuronal projection to sensory or motor neuropils; thus, neuroblast lineages coordinately produce sensory and motor circuitry. In addition, we show that neurons with shared hemilineage-temporal identity target pre- and postsynapse localization to similar positions in the neuropil, and that hemilineage-temporal cohorts share more common synaptic partners than that produced by neuropil proximity alone. Thus, temporal and hemilineage identity plays essential roles in establishing neuronal connectivity.
Results
Mapping neuronal developmental origin in a TEM reconstruction
To relate developmental mechanisms to circuit establishment mechanisms, we first needed to identify the developmental origin of neurons within a TEM reconstruction of the larval CNS (Ohyama et al., 2015), allowing us to quantify neuronal projections, synapse localization, and connectivity. We took three approaches. First, we generated GFP-marked clones for seven different neuroblasts, representing different spatial positions in the neuroblast array (Figure 1A, left). We imaged each clone by light microscopy in newly hatched larvae – the same stage used for the TEM reconstruction – so that we could match clonal morphology at the light and TEM levels (Figure 1B, C). All assayed neuroblast clones had a reproducible clonal morphology including the number of fascicles entering the neuropil, cell body position, and axon/dendrite morphology (Figure 1B; data not shown). We identified the parental neuroblast generating each clone by comparing clonal morphology to embryonic DiI single neuroblast clones (Bossing et al., 1996; Schmid et al., 1999; Schmidt et al., 1997), larval neuroblast clones (Birkholz et al., 2015; Lacin and Truman, 2016), and cell body position.
To identify each of the seven genetically labeled neuroblast clones in the TEM volume, we matched lineage-specific features present in both light and TEM analyses. We identified neurons that had clustered cell bodies, clone morphology matching that seen by light microscopy (Figure 1C), and one or two fascicles entering the neuropil (Figure 1C, D). The similarity in overall clone morphology between genetically marked clones and TEM reconstructed clones was obvious (compare Figure 1B, C). We validated these assignments using two methods. First, we used neuroblast-specific Gal4 lines (Lacin and Truman, 2016; Seroka and Doe, 2019) and multicolor flip out (MCFO) (Nern et al., 2015) to label individual neurons within each lineage. We found that in each case we could match the morphology of an MCFO-labeled single neuron from a known neuroblast to an identical single neuron in the same neuroblast clone within the TEM reconstruction (data not shown). Second, we reconstructed the same seven lineages in a second hemisegment, A1R. We observed similar neuron numbers, similar fascicles per clone, and similar clonal morphology (Figure 1E; data not shown). Thus, neuroblast lineages are highly stereotyped in left and right hemisegments.
We found that, while lineages have stereotyped and distinct morphology, they all had broad projections within the neuropil. We mapped the distribution of pre- and postsynapses for each neuroblast clone and found that, consistent with neuronal morphology, synapses were distributed widely across the neuropil (Figure 1F). We conclude that clonally related neurons project widely and localize synapses widely (Figure 1G, H).
In total, we have mapped 14 neuroblast clones in the TEM volume (seven in A1L, seven in A1R). These lineages contain 160 interneurons (80 each in A1l and A1r), containing 7794 presynapses and 19,468 postsynapses. We also include the previously traced sensory afferents and dendrites from all 32 motor neurons (Jovanic et al., 2016; Ohyama et al., 2015; Zarin et al., 2019). All data are publicly available from https://github.com/bjm5164/Mark2020_larval_development (copy archived at swh:1:rev:43e0a22c5381427aa6670c55ec4de76f5ad39568; Mark, 2020). We note that some of the earliest born neurons are not included either because their cell bodies are in contact with the neuropil and they do not fasciculate with clonal fascicles, precluding assignment to a specific neuroblast lineage, or they do not maintain marker expression at larval hatching. This can sometimes, but not always, lead to a gap between the deepest mapped neuron and the neuropil (Figure 1C). The morphology and function of the earliest born neurons will be described elsewhere. Nevertheless, the current data is sufficiently comprehensive for mapping developmental mechanisms to circuit assembly mechanisms.
In the following sections, we first analyze the relationship of developmental mechanisms to neuronal projections and synapse localization within the neuropil; we conclude by exploring the relationship between developmental mechanisms and neuronal connectivity.
Individual lineages generate two distinct morphological classes of neurons
Perhaps the most important neuropil subdivision in both mammalian spinal cord and Drosophila ventral nerve cord is the segregation of motor and sensory processing domains. In Drosophila, motor dendrites target a dorsal region of the neuropil, while sensory neurons target a ventral region of the neuropil (Landgraf et al., 2003; Mauss et al., 2009; Zarin et al., 2019; Zlatic et al., 2009). We extend these findings to show that premotor neurons also target pre- and postsynapses to the dorsal neuropil, and targets of sensory afferents target pre- and postsynapses to the ventral neuropil (Figure 2A–C). Thus, the neuropil has an important functional subdivision along the dorsoventral axis.
Upon examination of each lineage, we found that nearly all neuroblasts made two morphologically distinct types of neurons (except NB2-1; Figure 2—figure supplement 1). We used NBLAST (Costa et al., 2016) to quantify the morphological similarity of neurons within each neuroblast lineage and found that most lineages made two distinct morphological classes of neurons (or one class of neurons plus glia). Strikingly, each neuronal class projected to the dorsal/motor neuropil or the ventral/sensory neuropil (Figure 2D–H). Note that NB3-3 had only one fascicle at the NBLAST threshold used, and all neurons were ventral-projecting (Figure 2I). These classes were more distinct from one another than from neurons in other lineages (Figure 2—figure supplement 2). We conclude that most lineages generate two distinct classes of neurons that target either the ventral/sensory or dorsal/motor neuropils (Figure 2J).
Hemilineages produce sensory and motor processing neurons via a Notch-dependent mechanism
Recent work has shown that within the larval and adult CNS most neuroblast lineages generate NotchON neurons with a similar clonal morphology (called the NotchON hemilineage), and NotchOFF hemilineage with a different morphology (Harris et al., 2015; Lacin and Truman, 2016; Lee et al., 2020; Truman et al., 2010). We hypothesized that the observed morphological differences within our seven embryonic neuroblast lineages may be due to hemilineage identity. This hypothesis was strengthened by the fact that the sole neuroblast lineage that generated a single morphological class, NB3-3, had previously been shown to make only one ‘NotchOFF’ hemilineage via type 0 divisions (Baumgardt et al., 2014; Wreden et al., 2017). Based on these data, we hypothesize that each neuroblast makes a NotchON hemilineage projecting to the dorsal motor neuropil (or glia), and a NotchOFF hemilineage projecting to the ventral sensory neuropil.
To test this hypothesis, we took two approaches: first, we generated a lineage-specific reporter for NotchON neurons to see if NotchON neurons in specific neuroblast linages always project to the dorsal neuropil; and second, we tested whether misexpression of constitutively active Notch can redirect ventral projections into the dorsal neuropil. To generate a neuroblast lineage-specific NotchON reporter, we CRISPR engineered the Notch target gene hey, placing a T2A:FLP exon in frame with the terminal hey exon at the endogenous hey locus. Combining this fly with a neuroblast-specific Gal4 line (NB5-2-gal4 or NB7-1-gal4) plus a FLP-dependent reporter (UAS-FRT-stop-FRT-myr:sfGdP:HA) resulted in myr:GdP:HA expression specifically in NotchON neurons within a neuroblast lineage. This resulted in specific labeling of dorsal projecting neurons within the NB5-2 lineage (Figure 3A, A') or the NB7-1 lineage (Figure 3B, B'). We conclude that NotchON neurons specifically project to the dorsal motor neuropil.
We next asked whether Notch activity determines dorsal/ventral neuropil projections. We used Gal4 lines specifically expressed in single neuroblast lineages (NB1-2, NB5-2, NB7-1, or NB7-4) (Lacin and Truman, 2016; Seroka and Doe, 2019) to misexpress a constitutively active form of Notch (Notchintra) in individual neuroblast lineages. Wild-type lineages had both dorsal and ventral projections (Figure 3E) or glial and ventral projections (Figure 3F). In contrast, Notch expression throughout the NB1-2, NB5-2, or NB7-1 lineages led to a re-routing of projections from the ventral neuropil to the dorsal neuropil (Figure 3G–I). Notch expression throughout the NB7-4 lineage led to a loss of ventral projecting neurons and an increase in glia (Figure 3J). In addition, we note that ascending and descending projection neurons are normally generated by ventral NotchOFF hemilineages (Figure 3—figure supplement 1). Importantly, these ascending and descending projection neurons were completely lost following Notchintra expression (Figure 3G–I, insets). These results are summarized in Figure 3K. Our observation that NotchOFF hemilineages make more complex and lengthy neurons than the NotchON hemilineages is similar to that observed for larval brain lineages (Lee et al., 2020), suggesting that this is a conserved mechanism. The Notchintra phenotypes we observed are likely due to a NotchOFF to NotchON hemilineage transformation, rather than death of ventral projecting neurons, as we observed the same number of neurons in control versus Notchintra embryos at stage 16/17 (NB1-2 control: 12.2, n = 10 hemisegments; NB1-2 Notchintra: 14.1, n = 12 [p=0.0785]; NB5-2 control: 20.3, n = 16 hemisegments; NB5-2 Notchintra: 19.3, n = 26 [p=0.2311]; NB7-1 control: 18.9, n = 7 hemisegments; NB7-1 Notchintra: 19.8, n = 11 [p=0.6166]; NB7-4 control: 8.6, n = 9 hemisegments; NB7-4 Notchintra: 9.5, n = 11 [p=0.1710]). In conclusion, we show that NBLAST can be used to accurately identify neuroblast hemilineages; that NotchON/NotchOFF hemilineages project to motor/sensory neuropil domains, respectively; and most importantly, that hemilineage identity determines neuronal targeting to the motor or sensory neuropil. Thus, each neuroblast coordinately generates similar numbers of sensory and motor processing neurons throughout its embryonic lineage.
Hemilineages target synapses to subdomains of motor or sensory neuropil
To identify a relationship between hemilineage identity and synapse localization, we mapped the pre- and postsynapse localization for 12 bilateral hemilineages (24 total) in segment A1. Note that we show synapses from both A1L and A1R neuroblast lineages, which highlights the similarity and stereotypy of synapse localization between left and right lineages (Figure 4A). We found that the neurons in dorsal hemilineages localized both pre- and postsynaptic sites to the dorsal/motor neuropil, whereas neurons in ventral hemilineages localized both pre- and postsynaptic sites to the ventral/sensory neuropil (Figure 4A). Consistent with these observations, we found that the vast majority of sensory output was to ventral hemilineages, and the vast majority of input to motor neurons was from dorsal hemilineages (Figure 4B). We conclude that, at least for the assayed hemilineages, NotchON hemilineages target projections and pre- and postsynapses to the motor neuropil, whereas NotchOFF hemilineages target projections and pre- and postsynapses to the sensory neuropil (Figure 4H).
After showing that hemilineages target synapses to dorsal or ventral neuropil, we asked if individual hemilineages target distinct regions of the neuropil or if they have overlapping territories. We mapped the pre- and postsynaptic position for both dorsal and ventral hemilineages (excluding the NB7-4 glial hemilineage). We found that presynapses were localized to distinct regions of the neuropil (Figure 4C, D). Similarly, postsynapses were localized to distinct but slightly more overlapping regions of the neuropil (Figure 4E, F). We quantified synapse similarity (a measure of the average distance between synapses of two neurons) using published methods (Schlegel et al., 2016). We found that synapses from lineages in A1L and A1R had similar relative positions in the neuropil (Figure 4—figure supplement 1A–D), highlighting the stereotypy of synapse localization. Importantly, neurons in a hemilineage showed greater similarity in synaptic positions (pre or post) than unrelated neurons (Figure 4G, Figure 4—figure supplement 1E, F). We conclude that each hemilineage targets its presynapses (and to a lesser extent postsynapses) to small domains of the sensory or motor neuropil (Figure 4H), strongly suggesting that the developmental information needed for neuropil targeting is shared by neurons in a hemilineage, but not by all neurons in a complete neuroblast lineage (see Discussion).
Mapping temporal identity in a TEM reconstruction: radial position is a proxy for neuronal birth-order
To investigate the role of temporal identity in determining neuronal projections, synapse localization, or connectivity, we needed to identify the temporal identity of all 160 interneurons analyzed here. We used two methods. First, we confirmed that temporal transcription factors (Hb, Kr, Pdm, and Cas) have a radial distribution in the embryonic CNS, with early-born Hb+ neurons positioned in a deep layer adjacent to the developing neuropil and late-born Cas+ neurons are at the most superficial position within the CNS (Isshiki et al., 2001; Kambadur et al., 1998; Figure 5A). Importantly, we show that this distribution persists for more stable Hb and Cas reporters in the newly hatched larval CNS (Figure 5B, C), the stage of the TEM reconstruction. Thus, radial position can be used as a proxy for temporal identity in both embryos and newly hatched larvae (Figure 5D). Second, we used MCFO to identify additional single Hb+ or Cas+ neurons and matched them to the morphologically identical neuron in the TEM volume (Figure 5E ,F, Figure 5—figure supplement 1; data not shown). In total, we identified 54 Hb+ neurons and 44 Cas+ neurons within the TEM volume. We measured their cortex neurite length (neurite length from the cell body to the neuropil entry point; i.e., the length of the neurite as it traverses the cellular cortex) and found that experimentally verified Hb+ neurons were closer to the neuropil, whereas Cas+ neurons were further from the neuropil (Figure 5G). We conclude that neuronal radial position can be used as a proxy for neuronal temporal identity.
To determine the temporal identity of neurons within the seven bilateral neuroblast lineages, we measured cortex neurite length for each neuron. Importantly, the same neuron in A1L and A1R had similar cortex neurite lengths (Figure 5H), showing that cell body radial position was reproducible. We assigned 70 interneurons to one of four temporal cohorts (early, mid-early, mid-late, and late born) based on radial position (Figure 5I). We note that some of the earliest born neurons are not included (see Materials and methods), and we excluded the NB1-2 dorsal hemilineage from radial analysis as we found it to be an unreliable measure of birth-order for that hemilineage (see Figure 5—figure supplement 2). We conclude that cortex neurite length can be used as a proxy for temporal identity, that it is reproducible across at least two hemisegments, and that it can be used to approximate the temporal identity of any neuron in the TEM reconstruction. We use these temporal cohorts to explore the relationship between temporal identity and neuronal projections, synapse localization, and connectivity in the following sections.
After determining the temporal identity for 140 neurons, we asked whether neurons with a shared temporal identity had similar axon/dendrite projections or synapse localization. We found that neurons in each temporal class had broad neuronal projections, and no greater than unrelated neurons synapse similarity (Figure 5—figure supplement 3). We conclude that temporal identity, on its own, does not confer shared neuronal projections or synapse localization.
Hemilineage-temporal cohorts have distinct synapse localization within the neuropil
Although temporal identity alone did not correlate with axon projections or synapse localization, it still could play a role in restricting synapse localization within individual hemilineages. Here, we analyze neurons in a hemilineage that share a temporal identity, subsequently called hemilineage-temporal cohorts (Supplementary file 1). We mapped presynaptic and postsynaptic localization for 20 hemilineage-temporal cohorts from seven neuroblasts, as illustrated for the NB5-2 dorsal hemilineage and the NB3-3 ventral hemilineage (Figure 6A, B) and the remainder of the hemilineage-temporal cohorts (Figure 6—figure supplements 1 and 2). We found that NB5-2 dorsal hemilineage-temporal cohorts target their presynapses to mostly non-overlapping neuropil domains (Figure 6A), whereas the NB3-3 hemilineage-temporal cohorts target their postsynapses to mostly non-overlapping neuropil domains (Figure 6B). When we expand this analysis to all hemilineages, we found that (a) hemilineage-temporal cohorts had more similar synaptic positions than hemilineage alone, and (b) dorsal hemilineage-temporal cohorts preferentially clustered presynapses, whereas ventral hemilineage-temporal cohorts preferentially clustered postsynapses (Figure 6C–F). This latter observation may reflect the need to precisely target premotor output to specific regions of motor neuron dendrites and the need to precisely receive sensory input from distinct sensory modalities (see Discussion). Our results also suggest that dorsal and ventral hemilineages utilize temporal identity in different ways to specify the targeting of either presynapses or postsynapses. We conclude that temporal identity subdivides hemilineages into smaller populations of neurons that target both projections and synapses to distinct subdomains within the larger hemilineage targeting domain (Figure 6G). Thus, hemilineage identity provides coarse targeting within the neuropil, and temporal identity refines targeting to smaller subdomains.
Hemilineage-temporal cohorts have shared connectivity
Work from the Heckscher lab has shown that early-born and late-born temporal cohorts in the NB3-3 lineage have unique connectivity and participate in escape or proprioception locomotor circuits, respectively (Wreden et al., 2017). This has led them to hypothesize that other hemilineage-temporal cohorts may have shared connectivity (Meng and Heckscher, 2021). Below, we test this hypothesis for 20 hemilineage-temporal cohorts from seven different neuroblasts (Supplementary file 1). Here, we compare connectivity of hemilineage-temporal cohorts to unrelated neurons, neurons sharing a temporal identity, and neurons sharing a hemilineage identity. This allows us to test the hypothesis that hemilineage-temporal cohorts have more shared connectivity compared to other developmental groupings. First, we analyzed the connectome of 12 hemilineages plus the motor and sensory neurons in segment A1. In total, we analyzed 160 interneurons, 56 motor neurons, and 86 sensory neurons, which corresponded to approximately 25% of all inputs and 14% of all outputs for the 12 hemilineages. We found similar connectivity for the same neurons on the left and right side of the segment (Figure 7A–C), validating the accuracy of the neuronal reconstructions as well as the stereotypy of neurogenesis. Consistent with anatomical observations, we found that neurons in a hemilineage were clustered together in network space, suggesting that hemilineage-related neurons are highly interconnected and functionally distinct (Figure 7D–F), similar to what has been suggested previously (Harris et al., 2015).
Next, to determine if temporal identity could provide additional specificity to the observed hemilineage connectivity, we compared the connectivity of hemilineages to hemilineage-temporal cohorts. Here, we focus specifically on the 160 interneurons we have traced in the TEM volume, excluding motor and sensory neurons. We quantified the average connectivity distance between neurons related by hemilineage or hemilineage-temporal cohort. Neuron pairs that are directly connected have a network distance of one synapse; neurons that share a common input or output have a network distance of two synapses, up to an observed maximum of seven synapses. Importantly, neurons within a hemilineage had a lower network distance (i.e., greater connectivity) than unrelated neurons, while neurons within a hemilineage-temporal cohort had a lower network distance than those related by hemilineage alone (Figure 7G–J). Interestingly, we found that the largest difference between hemilineage and hemilineage-temporal cohort connectivity was at the two synapse distance. While less than 20% of hemilineage-related or hemilineage-temporal cohort-related neurons were directly connected, 75% of hemilineage-temporal cohort neurons were separated by two synapses or less (Figure 7I, J). We conclude that hemilineage-temporal cohorts share common connectivity and provide increased partner specificity compared to connectivity of hemilineages alone (Figure 7K).
Temporal cohorts have greater shared connectivity with each other than predicted by neuropil proximity alone
Next, we explored the relationship between temporal identity and connectivity. Peter’s rule proposes that connectivity is determined primarily by the overlap of axons and dendrites (Rees et al., 2017), raising the possibility that developmental mechanisms confer connectivity specificity simply by targeting sets of neurons to appropriate regions of the neuropil. To determine the role of neuropil proximity in connectivity, we calculated the spatial overlap of presynapses and postsynapses for all pairwise neuronal combinations, which ranged from no overlap to high overlap (Figure 8A, schematic). This measurement is analogous to previous synapse similarity measures, but instead we compared the presynapse positions of one neuron to the postsynapse positions of another. The majority of neurons have very little spatial overlap between presynapses and postsynapses (Figure 8A'). Interestingly, even neurons with the highest observed levels of overlap were not always connected (Figure 8A''). Thus, proximity alone cannot explain the observed connectivity, consistent with a role for hemilineage-temporal cohorts providing increased synaptic specificity.
To test the hypothesis that hemilineage-temporal identity confers synaptic specificity, we shuffled the connectivity by preserving the output degree of each neuron while setting the connection probability as a function of pre/post synaptic overlap. We found that hemilineage-temporal cohort to hemilineage-temporal cohort connectivity (Figure 8B, red dashed line) was greater than unrelated neuron wiring (Figure 8B, black) or proximity-based wiring (Figure 8B, cyan and magenta).
Finally, we asked if there was greater connectivity between hemilineage-temporal cohorts compared to unrelated neurons, neurons from the same temporal cohort, or neurons from the same hemilineage. We found that presynapses of a hemilineage-temporal cohort have greater shared connectivity than observed for unrelated neurons, temporal cohorts, or hemilineages (Figure 8C). In contrast, postsynapses of a hemilineage-temporal cohort have greater shared connectivity than observed for unrelated neurons and temporal cohorts, but not for hemilineages (Figure 8D). Moreover, pairwise comparisons revealed greater hemilineage-temporal cohort interconnectivity than hemilineages alone (Figure 8E, top). Furthermore, we detected greater connectivity between neurons within an individual hemilineage-temporal cohort than within an individual hemilineage (Figure 8E, bottom). These results suggest that within each hemilineage temporal identity confers an added level of connectivity compared to proximity or hemilineage alone. In conclusion, we propose that neuroblast lineage, hemilineage, and temporal identity function combinatorially to refine neurite projections, synapse localization, and connectivity (Figure 8F).
Hemilineage-temporal cohorts and circuit formation: the Eve proprioceptive circuit
We previously identified a proprioceptive motor circuit containing Eve lateral (EL) interneurons that is required to maintain left/right symmetry of body wall muscle contractions (Heckscher et al., 2015; Figure 9A). There are three Even-skipped (Eve)+ interneurons at the core of the circuit, receiving strong input from the proprioceptive sensory neurons and three local interneurons (Jaam1-3), and strong output to three premotor neurons (Saaghi1-3) (Heckscher et al., 2015). Recent work has shown that the three Eve+ interneurons at the core of the circuit are from a single hemilineage-temporal cohort: late-born neurons from the NB3-3 hemilineage (Wreden et al., 2017). This has led to speculation that the input Jaam neurons and the output Saaghi neurons may also be hemilineage-temporal cohorts (Meng and Heckscher, 2021). Fortunately, the Eve+, Jaam, and Saaghi neurons are all contained within our data set, allowing us to test this hypothesis. We confirm the results of Wreden et al., showing that the three Eve+ neurons (A08e1-3) are in a single late-born hemilineage-temporal cohort (Figure 9B; Supplementary file 1). In addition, we find that the Jaam neurons are all derived from a single hemilineage-temporal cohort: late-born neuron in the NB5-2 ventral (NotchOFF) hemilineage (Figure 9B; Supplementary file 1). Similarly, all three Saaghi neurons are derived from a single hemilineage-temporal cohort: late-born neuron in the NB5-2 dorsal (NotchON) hemilineage (Figure 9B; Supplementary file 1). Thus, the Jaam-EL-Saaghi proprioceptive circuit is assembled from three distinct hemilineage-temporal cohorts (Figure 9C). We propose that other motor circuits may also be assembled by preferential connectivity between distinct hemilineage-temporal cohorts.
Discussion
Here, we determine the relationship between developmental mechanisms (spatial, temporal, and hemilineage identity) and circuit assembly mechanisms (projections, synapse localization, and connectivity). To do this, we map both developmental and circuit features for 160 neuronal progeny of 14 neuroblast lineages in a serial section TEM reconstruction – this allows us to characterize neurons that share a developmental feature at single synapse resolution. It is important to note that we chose the seven neuroblasts in this study based on successful clone generation and availability of single neuroblast Gal4 lines, and thus there should be no bias towards a particular pattern of neurite projections, synapse localization, or connectivity. Our results show that individual neuroblast lineages have unique but broad axon and dendrite projections to both motor and sensory neuropil; hemilineages restrict projections and synapse localization to either motor or sensory neuropil; and distinct temporal identities within hemilineages provide additional specificity in synapse localization and connectivity. Thus, all three developmental mechanisms act combinatorially to progressively refine neurite projections, synapse localization, and connectivity (Figure 8F).
In mammals, clonally related neurons often have a similar location (Fekete et al., 1994; Mihalas and Hevner, 2018), morphology (Mihalas and Hevner, 2018; Wong and Rapaport, 2009), and connectivity (Yu et al., 2009). In contrast, we found that clonally related neurons project widely in the neuropil, to both sensory and motor domains, and thus lack shared morphology. Perhaps as brain size expands to contain an increasing number of progenitors, each clone takes on a more uniform structure and function. Yet the observation that each neuroblast clone had highly stereotyped projections suggests that neuroblast identity (determined by the spatial position of the neuroblast) determines neuroblast-specific projection patterns. Testing this functionally would require manipulating spatial patterning cues to duplicate a neuroblast and assay both duplicate lineages for similar projections and connectivity.
We found that hemilineages produce sensory and motor processing units via a Notch-dependent mechanism. Pioneering work on Drosophila third instar larval neuroblast lineages showed that each neuroblast lineage is composed of two hemilineages with different projection patterns and neurotransmitter expression (Harris et al., 2015; Lacin and Truman, 2016; Truman et al., 2010). We extend these studies to embryonic neuroblasts and show that Notch signaling determines motor versus sensory neuropil projections in all lineages examined. Surprisingly, the NotchON hemilineage always projected to the dorsal/motor neuropil, whereas the NotchOFF hemilineage always projected to the ventral/sensory neuropil. The relationship between the NotchON hemilineage projecting to the motor neuropil may be a common feature of all 30 segmental neuroblasts or it could be that the NotchON/NotchOFF provides a switch to allow each hemilineage to respond differently to dorsoventral guidance cues, with some projecting dorsally and some projecting ventrally. Analysis of additional neuroblast lineages will resolve this question. Another point to consider is the potential role of Notch in post-mitotic neurons (Crowner et al., 2003) as our experiments generated Notchintra misexpression in both newborn sibling neurons as well as mature post-mitotic neurons. Future work manipulating Notch levels specifically in mature post-mitotic neurons undergoing process outgrowth will be needed to identify the role of Notch in mature neurons, if any.
Elegant work has identified neuropil gradients of Slit and Netrin along the mediolateral axis (Zlatic et al., 2009), Semaphorins along the dorsoventral axis (Zlatic et al., 2009), and Wnt5 along the anteroposterior axis (Yoshikawa et al., 2003). Our finding that neurons in a hemilineage project to a common region of the neuropil strongly suggests that all neurons within a hemilineage respond in the same way to these global pathfinding cues. Conversely, our finding that neurons in different hemilineages target distinct regions of the neuropil suggests that each hemilineage expresses a different palette of guidance receptors, which enable them to respond differentially to the same global cues. For example, neurons in ventral hemilineages may express Plexin receptors to repel them from high Semaphorins in the dorsal neuropil.
Hemilineages have not been well described in vertebrate neurogenesis. Notch signaling within the Vsx1 + V2 progenitor lineage generates NotchOFF V2a excitatory interneurons and NotchON V2b inhibitory interneurons, which may be distinct hemilineages (Del Barrio et al., 2007; Francius et al., 2016; Peng et al., 2007; Skaggs et al., 2011). Interestingly, both V2a and V2b putative hemilineages contain molecularly distinct subclasses (Harris et al., 2019); our work raises the possibility that these subtypes arise from temporal patterning within the V2 lineage. In addition, NotchON/NotchOFF hemilineages may exist in the pineal photoreceptor lineage, where NotchON and NotchOFF populations specify cell-type identity (Cau et al., 2019).
Only recently have the role of hemilineages been tested for their functional properties. In adults, activation of each larval hemilineage from NB5-2 showed similar behavioral output, whereas each hemilineage from NB6-1 elicited different behaviors (Harris et al., 2015). Our previous work showed that the Eve+, Saaghi, and Jaam neurons are part of a proprioceptive circuit (Heckscher et al., 2015); here, we show that each class of neurons represents a hemilineage-temporal cohort. Note that the Jaam neurons process sensory input and are in a NotchOFF hemilineage, supporting our conclusion that NotchOFF hemilineages are devoted to sensory processing; the Saaghi premotor neurons are in a NotchON hemilineage consistent with their role in motor processing. Interestingly, both input and output neurons in this circuit arise from a common progenitor (NB5-2), which may generate late-born Jaam/Saaghi sibling neurons (Figure 9B). In the future, it would be interesting to determine if other sibling hemilineages are in a common circuit to generate a specific behavior.
Our hemilineage results have several implications. First, our results reveal that sensory and motor processing components of the neuropil are being built in parallel, with one half of every GMC division contributing to either sensory or motor networks. This would be an efficient mechanism to maintain sensory/motor balance as lineage lengths are modified over evolutionary time. Second, our results suggest that looking for molecular or morphological similarities in full neuroblast clones may be misleading due to the full neuroblast clone comprising two different hemilineages. For example, performing bulk RNAseq on all neurons in a neuroblast lineage is unlikely to reveal key regulators of pathfinding or synaptic connectivity due to the mixture of disparate neurons from the two hemilineages.
We used the cortex neurite length of neurons as a proxy for birth-order and shared temporal identity. We feel this is a good approximation, but it clearly does not precisely identify neurons born during each of the Hb, Kr, Pdm, Cas temporal transcription factor windows. Nevertheless, we had sufficient resolution to observe that neurons with the same temporal identity clustered their pre- or postsynapses, rather than localizing them uniformly through the hemilineage neuropil domain (Figure 6G). Interestingly, the three-dimensional location of each hemilineage temporal cohort synaptic cluster is identical on the left and right side of A1 (data not shown), ruling out the mechanism of stochastic self-avoidance (Zipursky and Grueber, 2013). Other possible mechanisms include hemilineage-temporal cohorts expressing different levels of the presynapse spacing cue Sequoia (Kulkarni et al., 2016; Petrovic and Hummel, 2008) or hemilineage-temporal cohorts exhibiting different responses to global patterning cues. Testing the function of temporal identity factors in synaptic tiling will require hemilineage-specific alteration of temporal identity, followed by assaying synapse localization within the neuropil.
Our results strongly suggest that hemilineage identity and temporal identity act combinatorially to allow small pools of neurons to target pre- and postsynapses to highly precise regions of the neuropil, thereby restricting synaptic partner choice. Yet precise neuropil targeting is not sufficient to explain connectivity as many similarly positioned axons and dendrites fail to form connections (Figure 8C). We favor the model that hemilineages direct gross neurite targeting to motor or sensory neuropil, whereas temporal identity acts combinatorially with each hemilineage to direct more precise neurite targeting and synaptic connectivity. Thus, the same temporal cue (e.g., Hb) could promote targeting of one pool of neurons in one hemilineage and another pool of neurons in an adjacent hemilineage. This limits the number of regulatory mechanisms needed to generate precise neuropil targeting and connectivity for all ~600 neurons in a segment of the larval CNS.
In conclusion, we demonstrate how developmental information can be integrated with connectomic data. We show that lineage information, hemilineage identity, and temporal identity can all be accurately predicted using morphological features (e.g., number of fascicles entering the neuropil for neuroblast clones and radial position for temporal cohorts). This both greatly accelerates the ability to identify neurons in a large EM volume as well as sets up a framework in which to study development using data typically intended for studying connectivity and function. We have used this framework to relate developmental mechanism to neuronal projections, synapse localization, and connectivity. We find that lineage, hemilineage, and temporal identity act sequentially to progressively refine neuronal projections, synapse localization, and connectivity, and our data supports a model where hemilineage-temporal cohorts are units of connectivity for assembling motor circuits.
Materials and methods
Transgenic fly stocks
Request a detailed protocolTransgenic lines were made by BestGene (Chino Hills, CA) or Genetivision (Houston, TX).
NB clone generation and lineage identification
Request a detailed protocolThe NB clones were generated with the following flies: hs-Flp.G5.PEST.Opt, dpn(FRT.stop)LexA.p65, 26XLexAop2-mCD8::GFP. The embryos were collected in 25°C for 3 hr and then incubated in 25°C for another 3 hr. The aged embryos were then submerged in 32°C water bath for 5 min heat shock and then incubated in 25°C until larvae hatched. The CNS of newly hatched larvae was dissected and mounted as previously described (Clark et al., 2016; Heckscher et al., 2014; Syed et al., 2017). The neuropil was stained with mouse anti-engrailed (RRID:AB_528224) (DSHB, 4D9) and Alexa Fluor 647 Phalloidin (ThermoFisher) by following the manufacturer’s protocol. The images were collected with Zeiss710 and processed with Imaris.
Lineages were identified in the EM volume by finding neurons with morphologies that matched the clonal morphology and then identifying their neuropil entry point. We then examined every neuron that entered the neuropil in the same fascicle. In most cases, every neuron in the fascicle had a morphology that matched the clonal morphology. In a small number of cases, the fascicles diverged slightly before the neuropil entry point. We verified the number of neurons by looking at fasciculating cell populations from at least two hemisegments (A1L and A1R). In some cases, we were able to identify a stereotyped number of cells across as many as four hemisegments, suggesting that fasciculation is stereotyped and reliable.
Hb+ single-cell clone generation
Request a detailed protocolHb+ single-cell clones were generated with the following flies: hs-Flp.G5.PEST.Opt, 13XLexAop2-IVS-myr::GFP, and LexAp65-T2A-hb (see below). The embryos were collected for 7 hr in 25°C, submerged in 32°C water bath for 10 min heat shock, and then followed the protocol as described above.
LexAp65-T2A-hb was generated by in-frame fusion of (FRT.stop)::LexA.P65::T2A to the N-terminus of the hb open reading frame with CRISPR-Cas9 gene editing. The ds-DNA donor vector for homology-directed repair was composed of left homologous arm (1000 bp), LexA.P65 (Pfeiffer et al., 2008), T2A (Nern et al., 2015), and the right homologous arm (1000 bp); the fragments were amplified with PCR and then assembled in pHD-DsRed (RRID:Addgene_51434) with NEBuilder (New England BioLabs). The gRNAs were generated from the vector pCFD5 (RRID:Addgene_73914) (Port et al., 2014) containing target sequence TGCATCTTGGCGGCTCTAGA and ACTACGAGCAGCACAACGCC. The ds-DNA donor vectors and gRNA vectors were co-injected into yw;nos-Cas9 (Kondo and Ueda, 2013) flies by BestGene. The selection marker 3xP3-DsRed was then removed in transgenic flies by hs-Cre.
Immunostaining and imaging
Request a detailed protocolStandard confocal microscopy, immunocytochemistry, and MCFO methods were performed as previously described for larvae (Clark et al., 2016; Heckscher et al., 2014; Syed et al., 2017) or adults (Nern et al., 2015; Pfeiffer et al., 2008). Secondary antibodies were from Jackson Immunoresearch (West Grove, PA) and used according to the manufacturer’s instructions. Confocal image stacks were acquired on Zeiss 700, 710, or 800 microscopes. Images were processed in Fiji (https://imagej.net/Fiji), Adobe Photoshop (Adobe, San Jose, CA), and Adobe Illustrator (Adobe). When adjustments to brightness and contrast were needed, they were applied to the entire image uniformly. Mosaic images to show different focal planes were assembled in Fiji or Photoshop.
Morphological analysis of lineages
Request a detailed protocolMorphological analysis was done using NBLAST and the NAT package (Costa et al., 2016), and analysis and figure generation were done using R. Neurons were preprocessed by pruning the most distal twigs (Strahler order 4), converting neurons to dot-props, and running an all-by-all NBLAST. For individual lineages, clusters were set using a cutoff of 3.0. In the case of NB2-1, where nearly every neuron shares a very similar morphology, we first confirmed the presence of a hemilineage using anti-Hey staining. After confirmation of a hemilineage, we next removed A02o and A02l since we could not find any clones that contained either an anterior projection (A02o) or a second contralateral projection (A02l). We reasoned that the hemilineages would represent the next largest morphological division (Figure 2—figure supplement 1).
Synaptic distributions and density analysis
Request a detailed protocolSynapse distribution plots and density contours were generated using MATLAB. Neuron synaptic and skeleton information was imported to MATLAB using pymaid (Schlegel et al., 2016). Cross-sectional synapse distribution plots were made by taking all synapse positions between the T3 and A2 segments as positional information becomes lost due to changes in brain shape beyond these bounds. Synapse distribution plots are 1D kernel density estimates. Sensory and motor density maps were made by taking the synapse positions of all sensory neurons entering the A1 nerve, and all motor neurons exiting the A1 nerve as well as all neurons with at least three synapses connected to one of these neurons. For sensory and motor maps as well as individual hemilineages, density plots are 2D kernel density estimates of all synapse positions across the neuropil. A cutoff of 60% was used to set the outermost contour. For lineage maps (Figure 4A), we used 80% as a cutoff. Polyadic synapses were counted as many times as they have targets. For synapse distribution plots, polyadic synapses are scaled by their number of targets.
Temporal cohort assignment
Request a detailed protocolCortex neurite length was calculated by converting the skeletonized neuronal arbor into a directed graph away from the soma and performing a depth-first-search of all vertices. The neuropil borders were defined by a previously created neuropil volume object (Costa et al., 2016). The neuropil entry point was defined as the first vertex within the neuropil volume object. Cortex neurite length was then the path length between the soma and the neuropil entry point. Neurons were binned into four groups defined by the positions of identified Hb+ and Cas+ cells. Early-born cells were defined as neurons with a cortex length <1 standard deviation above the mean Hb+ neurite length. The next group had cortex neurite lengths ≤1 standard deviation below the mean Cas+ neurite length. The final two groups were split at the mean Cas+ neurite length.
Synapse similarity measurements
Request a detailed protocolSynapse similarity was calculated as described previously (Schlegel et al., 2016):
where f(is,jk) is the mean synapse similarity between all synapses of neuron i and neuron j. dsk is the Euclidean distance between synapses s and k such that synapse k is the closest synapse of neuron j to synapse s of neuron i. σ is a bandwidth term that determines what is considered close. nis and njk are the fraction of synapses for neuron i and neuron j that are within ω of synapse s and synapse k, respectively. We used parameters ω = σ = 2000 nm. For presynaptic or postsynaptic comparisons, the score for a given pair of neurons was the average similarity between the left and right sides. In order to calculate pre/post overlap, we applied the same measure, instead using the presynapses of neuron i and the postsynapses of neuron j.
Connectivity similarity measurement
Request a detailed protocolConnectivity similarity was calculated as one minus the cosine distance between either the row vectors (outputs) or column vectors (inputs) of the binarized connectivity matrix for all neurons. For models of unrelated neuron connectivity, we used two methods to determine the individual connection probabilities. First, we used a standard degree-based method in which the number of inputs for a given neuron was divided by the total number of inputs received by that neuron. Second, we used the pre-postsynapse overlap scores. For a given neuron, its probability of connecting to another neuron was equal to the overlap score divided by the sum of all overlap scores for that neuron. Pairs with a score of 0 had a 0 probability of connecting. Thresholds were done by setting all values below the threshold to zero when determining connection probability.
Electron microscopy and CATMAID
Request a detailed protocolWe reconstructed neurons in CATMAID as previously described (Carreira-Rosario et al., 2018; Heckscher et al., 2015; Ohyama et al., 2015).
Figures
Request a detailed protocolFigures were generated using MATLAB, R, CATMAID, and Fiji, and edited in either Illustrator or Photoshop (Adobe).
Statistical analysis
Request a detailed protocolStatistical significance is denoted by asterisks: ****p<0.0001; ***p<0.001; **p<0.01; *p<0.05; n.s., not significant. All statistical analyses were done in MATLAB. When comparing two groups of quantitative data, an unpaired t-test was performed if data was normally distributed (determined using a one-sample Kolmogorov–Smirnov test) and Wilcoxon rank-sum test if the data was not normally distributed. Two-sample Kolmogorov–Smirnov tests were used on empirical distributions. Linear models were generated in MATLAB using lmfit.
Data availability
Request a detailed protocolAll data are publicly available from https://github.com/bjm5164/Mark2020_larval_development.
Data availability
All data are publicly available from https://github.com/bjm5164/Mark2020_larval_development (copy archived at https://archive.softwareheritage.org/swh:1:rev:43e0a22c5381427aa6670c55ec4de76f5ad39568).
References
-
ming is expressed in neuroblast sublineages and regulates gene expression in the Drosophila central nervous systemDevelopment 116:943–952.
-
Clonal analysis in the chicken retina reveals tangential dispersion of clonally related cellsDevelopmental Biology 166:666–682.https://doi.org/10.1006/dbio.1994.1346
-
Vsx1 transiently defines an early intermediate V2 interneuron precursor compartment in the mouse developing spinal cordFrontiers in Molecular Neuroscience 9:145.https://doi.org/10.3389/fnmol.2016.00145
-
Onecut factors and Pou2f2 regulate the distribution of V2 interneurons in the mouse developing spinal cordFrontiers in Cellular Neuroscience 13:184.https://doi.org/10.3389/fncel.2019.00184
-
Temporal fate specification and neural progenitor competence during developmentNature Reviews. Neuroscience 14:823–838.https://doi.org/10.1038/nrn3618
-
Mechanisms and molecules of neuronal wiring: a primerCold Spring Harbor Perspectives in Biology 3:a001727.https://doi.org/10.1101/cshperspect.a001727
-
SoftwareMark2020, version swh:1:rev:43e0a22c5381427aa6670c55ec4de76f5ad39568Software Heritage, https://archive.softwareheritage.org/swh:1:rev:43e0a22c5381427aa6670c55ec4de76f5ad39568.
-
Development of motor circuits: from neuronal stem cells and neuronal diversity to motor circuit assemblyCurrent Topics in Developmental Biology 142:409–442.https://doi.org/10.1016/bs.ctdb.2020.11.010
-
Weighing the Evidence in Peters' Rule: Does Neuronal Morphology Predict Connectivity?Trends in Neurosciences 40:63–71.https://doi.org/10.1016/j.tins.2016.11.007
-
Timing temporal transitions during brain developmentCurrent Opinion in Neurobiology 42:84–92.https://doi.org/10.1016/j.conb.2016.11.010
-
Clonal analysis of Drosophila embryonic neuroblasts: neural cell types, axon projections and muscle targetsDevelopment 126:4653–4689.
-
Genetic control of Drosophila nerve cord developmentCurrent Opinion in Neurobiology 13:8–15.https://doi.org/10.1016/s0959-4388(03)00007-2
-
The molecular basis of self-avoidanceAnnual Review of Neuroscience 36:547–568.https://doi.org/10.1146/annurev-neuro-062111-150414
Article and author information
Author details
Funding
National Institutes of Health (HD27056)
- Brandon Mark
- Sen-Lin Lai
- Aref Arzan Zarin
- Laurina Manning
- Heather Q Pollington
- Chris Q Doe
HHMI
- Brandon Mark
- Sen-Lin Lai
- Aref Arzan Zarin
- Laurina Manning
- Heather Q Pollington
- Chris Q Doe
HHMI
- Albert Cardona
- James W Truman
National Institutes of Health (T32HD007348-24)
- Brandon Mark
National Institutes of Health (5T32GM007413-43)
- Heather Q Pollington
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Haluk Lacin for unpublished fly lines. We thank Todd Laverty, Gerry Rubin, and Gerd Technau for fly stocks; Luis Sullivan, Emily Sales, John Tuthill, and Tim Warren for comments on the manuscript; Avinash Khandelwal and Laura Herren for annotating neurons; Keiko Hirono for generating transgenic constructs; and Keiko Hirono, Rita Yazejian, Dalton Lee, Cooper Doe, and Casey Doe for confocal imaging. Stocks were obtained from the Bloomington Drosophila Stock Center (NIH P40OD018537). Funding was provided by HHMI (CQD, BM, LM, AAZ), NIH HD27056 (CQD), and NIH T32HD007348-24 (BM).
Copyright
© 2021, Mark et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,954
- views
-
- 418
- downloads
-
- 53
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Developmental Biology
Cells called alveolar myofibroblasts, which have a central role in the development of the lung after birth, receive an orchestrated input from a range of different signaling pathways.
-
- Developmental Biology
Premature infants with bronchopulmonary dysplasia (BPD) have impaired alveolar gas exchange due to alveolar simplification and dysmorphic pulmonary vasculature. Advances in clinical care have improved survival for infants with BPD, but the overall incidence of BPD remains unchanged because we lack specific therapies to prevent this disease. Recent work has suggested a role for increased transforming growth factor-beta (TGFβ) signaling and myofibroblast populations in BPD pathogenesis, but the functional significance of each remains unclear. Here, we utilize multiple murine models of alveolar simplification and comparative single-cell RNA sequencing to identify shared mechanisms that could contribute to BPD pathogenesis. Single-cell RNA sequencing reveals a profound loss of myofibroblasts in two models of BPD and identifies gene expression signatures of increased TGFβ signaling, cell cycle arrest, and impaired proliferation in myofibroblasts. Using pharmacologic and genetic approaches, we find no evidence that increased TGFβ signaling in the lung mesenchyme contributes to alveolar simplification. In contrast, this is likely a failed compensatory response, since none of our approaches to inhibit TGFβ signaling protect mice from alveolar simplification due to hyperoxia while several make simplification worse. In contrast, we find that impaired myofibroblast proliferation is a central feature in several murine models of BPD, and we show that inhibiting myofibroblast proliferation is sufficient to cause pathologic alveolar simplification. Our results underscore the importance of impaired myofibroblast proliferation as a central feature of alveolar simplification and suggest that efforts to reverse this process could have therapeutic value in BPD.