A Growth-Based Framework for Leaf Shape Development and Diversity.

How do genes modify cellular growth to create morphological diversity? We study this problem in two related plants with differently shaped leaves: Arabidopsis thaliana (simple leaf shape) and Cardamine hirsuta (complex shape with leaflets). We use live imaging, modeling, and genetics to deconstruct these organ-level differences into their cell-level constituents: growth amount, direction, and differentiation. We show that leaf shape depends on the interplay of two growth modes: a conserved organ-wide growth mode that reflects differentiation; and a local, directional mode that involves the patterning of growth foci along the leaf edge. Shape diversity results from the distinct effects of two homeobox genes on these growth modes: SHOOTMERISTEMLESS broadens organ-wide growth relative to edge-patterning, enabling leaflet emergence, while REDUCED COMPLEXITY inhibits growth locally around emerging leaflets, accentuating shape differences created by patterning. We demonstrate the predictivity of our findings by reconstructing key features of C. hirsuta leaf morphology in A. thaliana. VIDEO ABSTRACT.

A deep learning-based toolkit for 3D nuclei segmentation and quantitative analysis in cellular and tissue context.

We present a new set of computational tools that enable accurate and widely applicable 3D segmentation of nuclei in various 3D digital organs. We have developed an approach for ground truth generation and iterative training of 3D nuclear segmentation models, which we applied to popular CellPose, PlantSeg and StarDist algorithms. We provide two high-quality models trained on plant nuclei that enable 3D segmentation of nuclei in datasets obtained from fixed or live samples, acquired from different plant and animal tissues, and stained with various nuclear stains or fluorescent protein-based nuclear reporters. We also share a diverse high-quality training dataset of about 10,000 nuclei. Furthermore, we advanced the MorphoGraphX analysis and visualization software by, among other things, providing a method for linking 3D segmented nuclei to their surrounding cells in 3D digital organs. We found that the nuclear-to-cell volume ratio varies between different ovule tissues and during the development of a tissue. Finally, we extended the PlantSeg 3D segmentation pipeline with a proofreading tool that uses 3D segmented nuclei as seeds to correct cell segmentation errors in difficult-to-segment tissues.

A mechanical transition from tension to buckling underlies the jigsaw puzzle shape morphogenesis of histoblasts in the Drosophila epidermis.

The polygonal shape of cells in proliferating epithelia is a result of the tensile forces of the cytoskeletal cortex and packing geometry set by the cell cycle. In the larval Drosophila epidermis, two cell populations, histoblasts and larval epithelial cells, compete for space as they grow on a limited body surface. They do so in the absence of cell divisions. We report a striking morphological transition of histoblasts during larval development, where they change from a tensed network configuration with straight cell outlines at the level of adherens junctions to a highly folded morphology. The apical surface of histoblasts shrinks while their growing adherens junctions fold, forming deep lobules. Volume increase of growing histoblasts is accommodated basally, compensating for the shrinking apical area. The folded geometry of apical junctions resembles elastic buckling, and we show that the imbalance between the shrinkage of the apical domain of histoblasts and the continuous growth of junctions triggers buckling. Our model is supported by laser dissections and optical tweezer experiments together with computer simulations. Our analysis pinpoints the ability of histoblasts to store mechanical energy to a much greater extent than most other epithelial cell types investigated so far, while retaining the ability to dissipate stress on the hours time scale. Finally, we propose a possible mechanism for size regulation of histoblast apical size through the lateral pressure of the epidermis, driven by the growth of cells on a limited surface. Buckling effectively compacts histoblasts at their apical plane and may serve to avoid physical harm to these adult epidermis precursors during larval life. Our work indicates that in growing nondividing cells, compressive forces, instead of tension, may drive cell morphology.

Cell biology of the leaf epidermis: Fate specification, morphogenesis, and coordination.

As the outermost layer of plants, the epidermis serves as a critical interface between plants and the environment. During leaf development, the differentiation of specialized epidermal cell types, including stomatal guard cells, pavement cells, and trichomes, occurs simultaneously, each providing unique and pivotal functions for plant growth and survival. Decades of molecular-genetic and physiological studies have unraveled key players and hormone signaling specifying epidermal differentiation. However, most studies focus on only one cell type at a time, and how these distinct cell types coordinate as a unit is far from well-comprehended. Here we provide a review on the current knowledge of regulatory mechanisms underpinning the fate specification, differentiation, morphogenesis, and positioning of these specialized cell types. Emphasis is given to their shared developmental origins, fate flexibility, as well as cell cycle and hormonal controls. Furthermore, we discuss computational modeling approaches to integrate how mechanical properties of individual epidermal cell types and entire tissue/organ properties mutually influence each other. We hope to illuminate the underlying mechanisms coordinating the cell differentiation that ultimately generate a functional leaf epidermis.

Hidden functional complexity in the flora of an early land ecosystem.

The first land ecosystems were composed of organisms considered simple in nature, yet the morphological diversity of their flora was extraordinary. The biological significance of this diversity remains a mystery largely due to the absence of feasible study approaches. To study the functional biology of Early Devonian flora, we have reconstructed extinct plants from fossilised remains in silico. We explored the morphological diversity of sporangia in relation to their mechanical properties using finite element method. Our approach highlights the impact of sporangia morphology on spore dispersal and adaptation. We discovered previously unidentified innovations among early land plants, discussing how different species might have opted for different spore dispersal strategies. We present examples of convergent evolution for turgor pressure resistance, achieved by homogenisation of stress in spherical sporangia and by torquing force in Tortilicaulis-like specimens. In addition, we show a potential mechanism for stress-assisted sporangium rupture. Our study reveals the deceptive complexity of this seemingly simple group of organisms. We leveraged the quantitative nature of our approach and constructed a fitness landscape to understand the different ecological niches present in the Early Devonian Welsh Borderland flora. By connecting morphology to functional biology, these findings facilitate a deeper understanding of the diversity of early land plants and their place within their ecosystem.

Aspm knockout ferret reveals an evolutionary mechanism governing cerebral cortical size.

The human cerebral cortex is distinguished by its large size and abundant gyrification, or folding. However, the evolutionary mechanisms that drive cortical size and structure are unknown. Although genes that are essential for cortical developmental expansion have been identified from the genetics of human primary microcephaly (a disorder associated with reduced brain size and intellectual disability) 1 , studies of these genes in mice, which have a smooth cortex that is one thousand times smaller than the cortex of humans, have provided limited insight. Mutations in abnormal spindle-like microcephaly-associated (ASPM), the most common recessive microcephaly gene, reduce cortical volume by at least 50% in humans2-4, but have little effect on the brains of mice5-9; this probably reflects evolutionarily divergent functions of ASPM10,11. Here we used genome editing to create a germline knockout of Aspm in the ferret (Mustela putorius furo), a species with a larger, gyrified cortex and greater neural progenitor cell diversity12-14 than mice, and closer protein sequence homology to the human ASPM protein. Aspm knockout ferrets exhibit severe microcephaly (25-40% decreases in brain weight), reflecting reduced cortical surface area without significant change in cortical thickness, as has been found in human patients3,4, suggesting that loss of 'cortical units' has occurred. The cortex of fetal Aspm knockout ferrets displays a very large premature displacement of ventricular radial glial cells to the outer subventricular zone, where many resemble outer radial glia, a subtype of neural progenitor cells that are essentially absent in mice and have been implicated in cerebral cortical expansion in primates12-16. These data suggest an evolutionary mechanism by which ASPM regulates cortical expansion by controlling the affinity of ventricular radial glial cells for the ventricular surface, thus modulating the ratio of ventricular radial glial cells, the most undifferentiated cell type, to outer radial glia, a more differentiated progenitor.

Early role for a Na+,K+-ATPase (ATP1A3) in brain development.

Osmotic equilibrium and membrane potential in animal cells depend on concentration gradients of sodium (Na+) and potassium (K+) ions across the plasma membrane, a function catalyzed by the Na+,K+-ATPase α-subunit. Here, we describe ATP1A3 variants encoding dysfunctional α3-subunits in children affected by polymicrogyria, a developmental malformation of the cerebral cortex characterized by abnormal folding and laminar organization. To gain cell-biological insights into the spatiotemporal dynamics of prenatal ATP1A3 expression, we built an ATP1A3 transcriptional atlas of fetal cortical development using mRNA in situ hybridization and transcriptomic profiling of ∼125,000 individual cells with single-cell RNA sequencing (Drop-seq) from 11 areas of the midgestational human neocortex. We found that fetal expression of ATP1A3 is most abundant to a subset of excitatory neurons carrying transcriptional signatures of the developing subplate, yet also maintains expression in nonneuronal cell populations. Moving forward a year in human development, we profiled ∼52,000 nuclei from four areas of an infant neocortex and show that ATP1A3 expression persists throughout early postnatal development, most predominantly in inhibitory neurons, including parvalbumin interneurons in the frontal cortex. Finally, we discovered the heteromeric Na+,K+-ATPase pump complex may form nonredundant cell-type-specific α-ß isoform combinations, including α3-ß1 in excitatory neurons and α3-ß2 in inhibitory neurons. Together, the developmental malformation phenotype of affected individuals and single-cell ATP1A3 expression patterns point to a key role for α3 in human cortex development, as well as a cell-type basis for pre- and postnatal ATP1A3-associated diseases.

The annotation and analysis of complex 3D plant organs using 3DCoordX.

A fundamental question in biology concerns how molecular and cellular processes become integrated during morphogenesis. In plants, characterization of 3D digital representations of organs at single-cell resolution represents a promising approach to addressing this problem. A major challenge is to provide organ-centric spatial context to cells of an organ. We developed several general rules for the annotation of cell position and embodied them in 3DCoordX, a user-interactive computer toolbox implemented in the open-source software MorphoGraphX. 3DCoordX enables rapid spatial annotation of cells even in highly curved biological shapes. Using 3DCoordX, we analyzed cellular growth patterns in organs of several species. For example, the data indicated the presence of a basal cell proliferation zone in the ovule primordium of Arabidopsis (Arabidopsis thaliana). Proof-of-concept analyses suggested a preferential increase in cell length associated with neck elongation in the archegonium of Marchantia (Marchantia polymorpha) and variations in cell volume linked to central morphogenetic features of a trap of the carnivorous plant Utricularia (Utricularia gibba). Our work demonstrates the broad applicability of the developed strategies as they provide organ-centric spatial context to cellular features in plant organs of diverse shape complexity.

RCL1 copy number variants are associated with a range of neuropsychiatric phenotypes.

Mendelian and early-onset severe psychiatric phenotypes often involve genetic variants having a large effect, offering opportunities for genetic discoveries and early therapeutic interventions. Here, the index case is an 18-year-old boy, who at 14 years of age had a decline in cognitive functioning over the course of a year and subsequently presented with catatonia, auditory and visual hallucinations, paranoia, aggression, mood dysregulation, and disorganized thoughts. Exome sequencing revealed a stop-gain mutation in RCL1 (NM_005772.4:c.370 C > T, p.Gln124Ter), encoding an RNA 3'-terminal phosphate cyclase-like protein that is highly conserved across eukaryotic species. Subsequent investigations across two academic medical centers identified eleven additional cases of RCL1 copy number variations (CNVs) with varying neurodevelopmental or psychiatric phenotypes. These findings suggest that dosage variation of RCL1 contributes to a range of neurological and clinical phenotypes.

How Cell Geometry and Cellular Patterning Influence Tissue Stiffness.

Cell growth in plants occurs due to relaxation of the cell wall in response to mechanical forces generated by turgor pressure. Growth can be anisotropic, with the principal direction of growth often correlating with the direction of lower stiffness of the cell wall. However, extensometer experiments on onion epidermal peels have shown that the tissue is stiffer in the principal direction of growth. Here, we used a combination of microextensometer experiments on epidermal onion peels and finite element method (FEM) modeling to investigate how cell geometry and cellular patterning affects mechanical measurements made at the tissue level. Simulations with isotropic cell-wall material parameters showed that the orientation of elongated cells influences tissue apparent stiffness, with the tissue appearing much softer in the transverse versus the longitudinal directions. Our simulations suggest that although extensometer experiments show that the onion tissue is stiffer when stretched in the longitudinal direction, the effect of cellular geometry means that the wall is in fact softer in this direction, matching the primary growth direction of the cells.

A recurrent, homozygous EMC10 frameshift variant is associated with a syndrome of developmental delay with variable seizures and dysmorphic features.

PURPOSE: The endoplasmic reticulum membrane complex (EMC) is a highly conserved, multifunctional 10-protein complex related to membrane protein biology. In seven families, we identified 13 individuals with highly overlapping phenotypes who harbor a single identical homozygous frameshift variant in EMC10. METHODS: Using exome, genome, and Sanger sequencing, a recurrent frameshift EMC10 variant was identified in affected individuals in an international cohort of consanguineous families. Multiple families were independently identified and connected via Matchmaker Exchange and internal databases. We assessed the effect of the frameshift variant on EMC10 RNA and protein expression and evaluated EMC10 expression in normal human brain tissue using immunohistochemistry. RESULTS: A homozygous variant EMC10 c.287delG (Refseq NM_206538.3, p.Gly96Alafs*9) segregated with affected individuals in each family, who exhibited a phenotypic spectrum of intellectual disability (ID) and global developmental delay (GDD), variable seizures and variable dysmorphic features (elongated face, curly hair, cubitus valgus, and arachnodactyly). The variant arose on two founder haplotypes and results in significantly reduced EMC10 RNA expression and an unstable truncated EMC10 protein. CONCLUSION: We propose that a homozygous loss-of-function variant in EMC10 causes a novel syndromic neurodevelopmental phenotype. Remarkably, the recurrent variant is likely the result of a hypermutable site and arose on distinct founder haplotypes.

Elongation of wood fibers combines features of diffuse and tip growth.

Xylem fibers are highly elongated cells that are key constituents of wood, play major physiological roles in plants, comprise an important terrestrial carbon reservoir, and thus have enormous ecological and economic importance. As they develop, from fusiform initials, their bodies remain the same length while their tips elongate and intrude into intercellular spaces. To elucidate mechanisms of tip elongation, we studied the cell wall along the length of isolated, elongating aspen xylem fibers and used computer simulations to predict the forces driving the intercellular space formation required for their growth. We found pectin matrix epitopes (JIM5, LM7) concentrated at the tips where cellulose microfibrils have transverse orientation, and xyloglucan epitopes (CCRC-M89, CCRC-M58) in fiber bodies where microfibrils are disordered. These features are accompanied by changes in cell wall thickness, indicating that while the cell wall elongates strictly at the tips, it is deposited all over fibers. Computer modeling revealed that the intercellular space formation needed for intrusive growth may only require targeted release of cell adhesion, which allows turgor pressure in neighboring fiber cells to 'round' the cells creating spaces. These characteristics show that xylem fibers' elongation involves a distinct mechanism that combines features of both diffuse and tip growth.

A Spatiotemporal DNA Endoploidy Map of the Arabidopsis Root Reveals Roles for the Endocycle in Root Development and Stress Adaptation.

Somatic polyploidy caused by endoreplication is observed in arthropods, molluscs, and vertebrates but is especially prominent in higher plants, where it has been postulated to be essential for cell growth and fate maintenance. However, a comprehensive understanding of the physiological significance of plant endopolyploidy has remained elusive. Here, we modeled and experimentally verified a high-resolution DNA endoploidy map of the developing Arabidopsis thaliana root, revealing a remarkable spatiotemporal control of DNA endoploidy levels across tissues. Fitting of a simplified model to publicly available data sets profiling root gene expression under various environmental stress conditions suggested that this root endoploidy patterning may be stress-responsive. Furthermore, cellular and transcriptomic analyses revealed that inhibition of endoreplication onset alters the nuclear-to-cellular volume ratio and the expression of cell wall-modifying genes, in correlation with the appearance of cell structural changes. Our data indicate that endopolyploidy might serve to coordinate cell expansion with structural stability and that spatiotemporal endoreplication pattern changes may buffer for stress conditions, which may explain the widespread occurrence of the endocycle in plant species growing in extreme or variable environments.

Clones of cells switch from reduction to enhancement of size variability in Arabidopsis sepals.

Organs form with remarkably consistent sizes and shapes during development, whereas a high variability in growth is observed at the cell level. Given this contrast, it is unclear how such consistency in organ scale can emerge from cellular behavior. Here, we examine an intermediate scale, the growth of clones of cells in Arabidopsis sepals. Each clone consists of the progeny of a single progenitor cell. At early stages, we find that clones derived from a small progenitor cell grow faster than those derived from a large progenitor cell. This results in a reduction in clone size variability, a phenomenon we refer to as size uniformization. By contrast, at later stages of clone growth, clones change their growth pattern to enhance size variability, when clones derived from larger progenitor cells grow faster than those derived from smaller progenitor cells. Finally, we find that, at early stages, fast growing clones exhibit greater cell growth heterogeneity. Thus, cellular variability in growth might contribute to a decrease in the variability of clones throughout the sepal.

Regulation of human cerebral cortical development by EXOC7 and EXOC8, components of the exocyst complex, and roles in neural progenitor cell proliferation and survival.

PURPOSE: The exocyst complex is a conserved protein complex that mediates fusion of intracellular vesicles to the plasma membrane and is implicated in processes including cell polarity, cell migration, ciliogenesis, cytokinesis, autophagy, and fusion of secretory vesicles. The essential role of these genes in human genetic disorders, however, is unknown. METHODS: We performed homozygosity mapping and exome sequencing of consanguineous families with recessively inherited brain development disorders. We modeled an EXOC7 splice variant in vitro and examined EXOC7 messenger RNA (mRNA) expression in developing mouse and human cortex. We modeled exoc7 loss-of-function in a zebrafish knockout. RESULTS: We report variants in exocyst complex members, EXOC7 and EXOC8, in a novel disorder of cerebral cortex development. In EXOC7, we identified four independent partial loss-of-function (LOF) variants in a recessively inherited disorder characterized by brain atrophy, seizures, and developmental delay, and in severe cases, microcephaly and infantile death. In EXOC8, we found a homozygous truncating variant in a family with a similar clinical disorder. We modeled exoc7 deficiency in zebrafish and found the absence of exoc7 causes microcephaly. CONCLUSION: Our results highlight the essential role of the exocyst pathway in normal cortical development and how its perturbation causes complex brain disorders.

Phyllotaxis: is the golden angle optimal for light capture?

Phyllotactic patterns are some of the most conspicuous in nature. To create these patterns plants must control the divergence angle between the appearance of successive organs, sometimes to within a fraction of a degree. The most common angle is the Fibonacci or golden angle, and its prevalence has led to the hypothesis that it has been selected by evolution as optimal for plants with respect to some fitness benefits, such as light capture. We explore arguments for and against this idea with computer models. We have used both idealized and scanned leaves from Arabidopsis thaliana and Cardamine hirsuta to measure the overlapping leaf area of simulated plants after varying parameters such as leaf shape, incident light angles, and other leaf traits. We find that other angles generated by Fibonacci-like series found in nature are equally optimal for light capture, and therefore should be under similar evolutionary pressure. Our findings suggest that the iterative mechanism for organ positioning itself is a more likely target for evolutionary pressure, rather than a specific divergence angle, and our model demonstrates that the heteroblastic progression of leaf shape in A. thaliana can provide a potential fitness benefit via light capture.

Early immune responses are independent of RGC dysfunction in glaucoma with complement component C3 being protective.

Various immune response pathways are altered during early, predegenerative stages of glaucoma; however, whether the early immune responses occur secondarily to or independently of neuronal dysfunction is unclear. To investigate this relationship, we used the Wlds allele, which protects from axon dysfunction. We demonstrate that DBA/2J.Wlds mice develop high intraocular pressure (IOP) but are protected from retinal ganglion cell (RGC) dysfunction and neuroglial changes that otherwise occur early in DBA/2J glaucoma. Despite this, immune pathways are still altered in DBA/2J.Wlds mice. This suggests that immune changes are not secondary to RGC dysfunction or altered neuroglial interactions, but may be directly induced by the increased strain imposed by high IOP. One early immune response following IOP elevation is up-regulation of complement C3 in astrocytes of DBA/2J and DBA/2J.Wlds mice. Unexpectedly, because the disruption of other complement components, such as C1Q, is protective in glaucoma, C3 deficiency significantly increased the number of DBA/2J eyes with nerve damage and RGC loss at an early time point after IOP elevation. Transcriptional profiling of C3-deficient cultured astrocytes implicated EGFR signaling as a hub in C3-dependent responses. Treatment with AG1478, an EGFR inhibitor, also significantly increased the number of DBA/2J eyes with glaucoma at the same early time point. These findings suggest that C3 protects from early glaucomatous damage, a process that may involve EGFR signaling and other immune responses in the optic nerve head. Therefore, therapies that target specific components of the complement cascade, rather than global inhibition, may be more applicable for treating human glaucoma.

A Model of Differential Growth-Guided Apical Hook Formation in Plants.

Differential cell growth enables flexible organ bending in the presence of environmental signals such as light or gravity. A prominent example of the developmental processes based on differential cell growth is the formation of the apical hook that protects the fragile shoot apical meristem when it breaks through the soil during germination. Here, we combined in silico and in vivo approaches to identify a minimal mechanism producing auxin gradient-guided differential growth during the establishment of the apical hook in the model plant Arabidopsis thaliana Computer simulation models based on experimental data demonstrate that asymmetric expression of the PIN-FORMED auxin efflux carrier at the concave (inner) versus convex (outer) side of the hook suffices to establish an auxin maximum in the epidermis at the concave side of the apical hook. Furthermore, we propose a mechanism that translates this maximum into differential growth, and thus curvature, of the apical hook. Through a combination of experimental and in silico computational approaches, we have identified the individual contributions of differential cell elongation and proliferation to defining the apical hook and reveal the role of auxin-ethylene crosstalk in balancing these two processes.

Actin-dependent vacuolar occupancy of the cell determines auxin-induced growth repression.

The cytoskeleton is an early attribute of cellular life, and its main components are composed of conserved proteins. The actin cytoskeleton has a direct impact on the control of cell size in animal cells, but its mechanistic contribution to cellular growth in plants remains largely elusive. Here, we reveal a role of actin in regulating cell size in plants. The actin cytoskeleton shows proximity to vacuoles, and the phytohormone auxin not only controls the organization of actin filaments but also impacts vacuolar morphogenesis in an actin-dependent manner. Pharmacological and genetic interference with the actin-myosin system abolishes the effect of auxin on vacuoles and thus disrupts its negative influence on cellular growth. SEM-based 3D nanometer-resolution imaging of the vacuoles revealed that auxin controls the constriction and luminal size of the vacuole. We show that this actin-dependent mechanism controls the relative vacuolar occupancy of the cell, thus suggesting an unanticipated mechanism for cytosol homeostasis during cellular growth.

Seasonal Regulation of Petal Number.

Four petals characterize the flowers of most species in the Brassicaceae family, and this phenotype is generally robust to genetic and environmental variation. A variable petal number distinguishes the flowers of Cardamine hirsuta from those of its close relative Arabidopsis (Arabidopsis thaliana), and allelic variation at many loci contribute to this trait. However, it is less clear whether C. hirsuta petal number varies in response to seasonal changes in environment. To address this question, we assessed whether petal number responds to a suite of environmental and endogenous cues that regulate flowering time in C. hirsuta We found that petal number showed seasonal variation in C. hirsuta, such that spring flowering plants developed more petals than those flowering in summer. Conditions associated with spring flowering, including cool ambient temperature, short photoperiod, and vernalization, all increased petal number in C. hirsuta Cool temperature caused the strongest increase in petal number and lengthened the time interval over which floral meristems matured. We performed live imaging of early flower development and showed that floral buds developed more slowly at 15°C versus 20°C. This extended phase of floral meristem formation, coupled with slower growth of sepals at 15°C, produced larger intersepal regions with more space available for petal initiation. In summary, the growth and maturation of floral buds is associated with variable petal number in C. hirsuta and responds to seasonal changes in ambient temperature.