Interplay of cell shape and division orientation promotes robust morphogenesis of developing epithelia.

Epithelial cells acquire functionally important shapes (e.g., squamous, cuboidal, columnar) during development. Here, we combine theory, quantitative imaging, and perturbations to analyze how tissue geometry, cell divisions, and mechanics interact to shape the presumptive enveloping layer (pre-EVL) on the zebrafish embryonic surface. We find that, under geometrical constraints, pre-EVL flattening is regulated by surface cell number changes following differentially oriented cell divisions. The division pattern is, in turn, determined by the cell shape distribution, which forms under geometrical constraints by cell-cell mechanical coupling. An integrated mathematical model of this shape-division feedback loop recapitulates empirical observations. Surprisingly, the model predicts that cell shape is robust to changes of tissue surface area, cell volume, and cell number, which we confirm in vivo. Further simulations and perturbations suggest the parameter linking cell shape and division orientation contributes to epithelial diversity. Together, our work identifies an evolvable design logic that enables robust cell-level regulation of tissue-level development.

Specified neural progenitors sort to form sharp domains after noisy Shh signaling.

Sharply delineated domains of cell types arise in developing tissues under instruction of inductive signal (morphogen) gradients, which specify distinct cell fates at different signal levels. The translation of a morphogen gradient into discrete spatial domains relies on precise signal responses at stable cell positions. However, cells in developing tissues undergoing morphogenesis and proliferation often experience complex movements, which may affect their morphogen exposure, specification, and positioning. How is a clear pattern achieved with cells moving around? Using in toto imaging of the zebrafish neural tube, we analyzed specification patterns and movement trajectories of neural progenitors. We found that specified progenitors of different fates are spatially mixed following heterogeneous Sonic Hedgehog signaling responses. Cell sorting then rearranges them into sharply bordered domains. Ectopically induced motor neuron progenitors also robustly sort to correct locations. Our results reveal that cell sorting acts to correct imprecision of spatial patterning by noisy inductive signals.

Otolith tethering in the zebrafish otic vesicle requires Otogelin and α-Tectorin.

Otoliths are biomineralised structures important for balance and hearing in fish. Their counterparts in the mammalian inner ear, otoconia, have a primarily vestibular function. Otoliths and otoconia form over sensory maculae and are attached to the otolithic membrane, a gelatinous extracellular matrix that provides a physical coupling between the otolith and the underlying sensory epithelium. In this study, we have identified two proteins required for otolith tethering in the zebrafish ear, and propose that there are at least two stages to this process: seeding and maintenance. The initial seeding step, in which otolith precursor particles tether directly to the tips of hair cell kinocilia, fails to occur in the einstein (eis) mutant. The gene disrupted in eis is otogelin (otog); mutations in the human OTOG gene have recently been identified as causative for deafness and vestibular dysfunction (DFNB18B). At later larval stages, maintenance of otolith tethering to the saccular macula is dependent on tectorin alpha (tecta) function, which is disrupted in the rolling stones (rst) mutant. α-Tectorin (Tecta) is a major constituent of the tectorial membrane in the mammalian cochlea. Mutations in the human TECTA gene can cause either dominant (DFNA8/12) or recessive (DFNB21) forms of deafness. Our findings indicate that the composition of extracellular otic membranes is highly conserved between mammals and fish, reinforcing the view that the zebrafish is an excellent model system for the study of deafness and vestibular disease.

Recovery of shape and size in a developing organ pair.

BACKGROUND: Paired organs in animals are largely bilaterally symmetric despite inherent noise in most biological processes. How is precise organ shape and size achieved during development despite this noise? Examining paired organ development is a challenge because it requires repeated quantification of two structures in parallel within living embryos. Here we combine bilateral quantification of morphology through time with asymmetric perturbations to study regulation of organ shape, size, and symmetry in developing organ pairs. RESULTS: We present quantitative live imaging tools to measure the shape and size of the developing inner ears on both the left and right side simultaneously over time. By quantifying variation between the left and right inner ear (intrinsic noise) and between different individuals (extrinsic noise), we find that initial variability decreases over time in normal development to achieve symmetry. Early asymmetry is increased by environmental stress, but symmetry is still recovered over subsequent developmental time. Using multiple unilateral perturbations including Fgf signaling and ultraviolet light, we find that growth can be adjusted to compensate for a range of initial size and shape differences. CONCLUSIONS: We propose that symmetry in developmental systems does not emerge through precise deterministic bilateral development, but rather through feedback mechanisms that adjust morphogenesis rates to account for variation. Developmental Dynamics 246:451-465, 2016. © 2017 Wiley Periodicals, Inc.

RNA-seq-based mapping and candidate identification of mutations from forward genetic screens.

Forward genetic screens have elucidated molecular pathways required for innumerable aspects of life; however, identifying the causal mutations from such screens has long been the bottleneck in the process, particularly in vertebrates. We have developed an RNA-seq-based approach that identifies both the region of the genome linked to a mutation and candidate lesions that may be causal for the phenotype of interest. We show that our method successfully identifies zebrafish mutations that cause nonsense or missense changes to codons, alter transcript splicing, or alter gene expression levels. Furthermore, we develop an easily accessible bioinformatics pipeline allowing for implementation of all steps of the method. Overall, we show that RNA-seq is a fast, reliable, and cost-effective method to map and identify mutations that will greatly facilitate the power of forward genetics in vertebrate models.

Rapid positional cloning of zebrafish mutations by linkage and homozygosity mapping using whole-genome sequencing.

Forward genetic screens in zebrafish have identified >9000 mutants, many of which are potential disease models. Most mutants remain molecularly uncharacterized because of the high cost, time and labor investment required for positional cloning. These costs limit the benefit of previous genetic screens and discourage future screens. Drastic improvements in DNA sequencing technology could dramatically improve the efficiency of positional cloning in zebrafish and other model organisms, but the best strategy for cloning by sequencing has yet to be established. Using four zebrafish inner ear mutants, we developed and compared two approaches for 'cloning by sequencing': one based on bulk segregant linkage (BSFseq) and one based on homozygosity mapping (HMFseq). Using BSFseq we discovered that mutations in lmx1b and jagged1b cause abnormal ear morphogenesis. With HMFseq we validated that the disruption of cdh23 abolishes the ear's sensory functions and identified a candidate lesion in lhfpl5a predicted to cause nonsyndromic deafness. The success of HMFseq shows that the high intrastrain polymorphism rate in zebrafish eliminates the need for time-consuming map crosses. Additionally, we analyzed diversity in zebrafish laboratory strains to find areas of elevated diversity and areas of fixed homozygosity, reinforcing recent findings that genome diversity is clustered. We present a database of >15 million sequence variants that provides much of this approach's power. In our four test cases, only a single candidate single nucleotide polymorphism (SNP) remained after subtracting all database SNPs from a mutant's critical region. The saturation of the common SNP database and our open source analysis pipeline MegaMapper will improve the pace at which the zebrafish community makes unique discoveries relevant to human health.

Ribeye is required for presynaptic Ca(V)1.3a channel localization and afferent innervation of sensory hair cells.

Ribbon synapses of the ear, eye and pineal gland contain a unique protein component: Ribeye. Ribeye consists of a novel aggregation domain spliced to the transcription factor CtBP2 and is one of the most abundant proteins in synaptic ribbon bodies. Although the importance of Ribeye for the function and physical integrity of ribbon synapses has been shown, a specific role in synaptogenesis has not been described. Here, we have modulated Ribeye expression in zebrafish hair cells and have examined the role of Ribeye in synapse development. Knockdown of ribeye resulted in fewer stimulus-evoked action potentials from afferent neurons and loss of presynaptic Ca(V)1.3a calcium channel clusters in hair cells. Additionally, afferent innervation of hair cells was reduced in ribeye morphants, and the reduction was correlated with depletion of Ribeye punctae. By contrast, transgenic overexpression of Ribeye resulted in Ca(V)1.3a channels colocalized with ectopic aggregates of Ribeye protein. Overexpression of Ribeye, however, was not sufficient to create ectopic synapses. These findings reveal two distinct functions of Ribeye in ribbon synapse formation--clustering Ca(V)1.3a channels at the presynapse and stabilizing contacts with afferent neurons--and suggest that Ribeye plays an organizing role in synaptogenesis.

Synaptojanin1 is required for temporal fidelity of synaptic transmission in hair cells.

To faithfully encode mechanosensory information, auditory/vestibular hair cells utilize graded synaptic vesicle (SV) release at specialized ribbon synapses. The molecular basis of SV release and consequent recycling of membrane in hair cells has not been fully explored. Here, we report that comet, a gene identified in an ENU mutagenesis screen for zebrafish larvae with vestibular defects, encodes the lipid phosphatase Synaptojanin 1 (Synj1). Examination of mutant synj1 hair cells revealed basal blebbing near ribbons that was dependent on Cav1.3 calcium channel activity but not mechanotransduction. Synaptojanin has been previously implicated in SV recycling; therefore, we tested synaptic transmission at hair-cell synapses. Recordings of post-synaptic activity in synj1 mutants showed relatively normal spike rates when hair cells were mechanically stimulated for a short period of time at 20 Hz. In contrast, a sharp decline in the rate of firing occurred during prolonged stimulation at 20 Hz or stimulation at a higher frequency of 60 Hz. The decline in spike rate suggested that fewer vesicles were available for release. Consistent with this result, we observed that stimulated mutant hair cells had decreased numbers of tethered and reserve-pool vesicles in comparison to wild-type hair cells. Furthermore, stimulation at 60 Hz impaired phase locking of the postsynaptic activity to the mechanical stimulus. Following prolonged stimulation at 60 Hz, we also found that mutant synj1 hair cells displayed a striking delay in the recovery of spontaneous activity. Collectively, the data suggest that Synj1 is critical for retrieval of membrane in order to maintain the quantity, timing of fusion, and spontaneous release properties of SVs at hair-cell ribbon synapses.

An oncogenic enhancer encodes selective selenium dependency in AML.

Deregulation of transcription is a hallmark of acute myeloid leukemia (AML) that drives oncogenic expression programs and presents opportunities for therapeutic targeting. By integrating comprehensive pan-cancer enhancer landscapes with genetic dependency mapping, we find that AML-enriched enhancers encode for more selective tumor dependencies. We hypothesized that this approach could identify actionable dependencies downstream of oncogenic driver events and discovered a MYB-regulated AML-enriched enhancer regulating SEPHS2, a key component of the selenoprotein production pathway. Using a combination of patient samples and mouse models, we show that this enhancer upregulates SEPHS2, promoting selenoprotein production and antioxidant function required for AML survival. SEPHS2 and other selenoprotein pathway genes are required for AML growth in vitro. SEPHS2 knockout and selenium dietary restriction significantly delay leukemogenesis in vivo with little effect on normal hematopoiesis. These data validate the utility of enhancer mapping in target identification and suggest that selenoprotein production is an actionable target in AML.

Vesicular glutamate transporter 3 is required for synaptic transmission in zebrafish hair cells.

Hair cells detect sound and movement and transmit this information via specialized ribbon synapses. Here we report that asteroid, a gene identified in an ethylnitrosourea mutagenesis screen of zebrafish larvae for auditory/vestibular mutants, encodes vesicular glutamate transporter 3 (Vglut3). A splice site mutation in exon 2 of vglut3 results in a severe truncation of the predicted protein product and morpholinos directed against the vglut3 ATG start site or the affected splice junction replicate the asteroid phenotype. In situ hybridization shows that vglut3 is exclusively expressed in hair cells of the ear and lateral line organ. A second transporter gene, vglut1, is also expressed in zebrafish hair cells, but the level of vglut1 mRNA is not increased in the absence of Vglut3. Antibodies against Vglut3 label the basal end of hair cells and labeling is not present in asteroid/vglut3 mutants. Based on the localization of Vglut3 in hair cells, we suspected that the lack of vestibulo-ocular and acoustic startle reflexes in asteroid/vglut3 mutants was attributable to a defect in synaptic transmission in hair cells. In support of this notion, action currents in postsynaptic acousticolateralis neurons are absent in asteroid/vglut3 mutants. At the ultrastructural level, mutant asteroid/vglut3 hair cells show a decrease in the number of ribbon-associated synaptic vesicles, indicating a role for Vglut3 in synaptic vesicle biogenesis and/or tethering to the ribbon body. Lack of postsynaptic action currents in the mutants suggests that the remaining hair-cell synaptic vesicles contain insufficient levels of glutamate for generation of action potentials in first-order neurons.

Size control of the inner ear via hydraulic feedback.

Animals make organs of precise size, shape, and symmetry but how developing embryos do this is largely unknown. Here, we combine quantitative imaging, physical theory, and physiological measurement of hydrostatic pressure and fluid transport in zebrafish to study size control of the developing inner ear. We find that fluid accumulation creates hydrostatic pressure in the lumen leading to stress in the epithelium and expansion of the otic vesicle. Pressure, in turn, inhibits fluid transport into the lumen. This negative feedback loop between pressure and transport allows the otic vesicle to change growth rate to control natural or experimentally-induced size variation. Spatiotemporal patterning of contractility modulates pressure-driven strain for regional tissue thinning. Our work connects molecular-driven mechanisms, such as osmotic pressure driven strain and actomyosin tension, to the regulation of tissue morphogenesis via hydraulic feedback to ensure robust control of organ size. Editorial note: This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

Myc and Fgf Are Required for Zebrafish Neuromast Hair Cell Regeneration.

Unlike mammals, the non-mammalian vertebrate inner ear can regenerate the sensory cells, hair cells, either spontaneously or through induction after hair cell loss, leading to hearing recovery. The mechanisms underlying the regeneration are poorly understood. By microarray analysis on a chick model, we show that chick hair cell regeneration involves the activation of proliferation genes and downregulation of differentiation genes. Both MYC and FGF are activated in chick hair cell regeneration. Using a zebrafish lateral line neuromast hair cell regeneration model, we show that the specific inhibition of Myc or Fgf suppresses hair cell regeneration, demonstrating that both pathways are essential to the process. Rapid upregulation of Myc and delayed Fgf activation during regeneration suggest a role of Myc in proliferation and Fgf in differentiation. The dorsal-ventral pattern of fgfr1a in the neuromasts overlaps with the distribution of hair cell precursors. By laser ablation, we show that the fgfr1a-positive supporting cells are likely the hair cell precursors that directly give rise to new hair cells; whereas the anterior-posterior fgfr1a-negative supporting cells have heightened proliferation capacity, likely to serve as more primitive progenitor cells to replenish lost precursors after hair cell loss. Thus fgfr1a is likely to mark compartmentalized supporting cell subtypes with different capacities in renewal proliferation and hair cell regeneration. Manipulation of c-MYC and FGF pathways could be explored for mammalian hair cell regeneration.

Type I interferons and microbial metabolites of tryptophan modulate astrocyte activity and central nervous system inflammation via the aryl hydrocarbon receptor.

Astrocytes have important roles in the central nervous system (CNS) during health and disease. Through genome-wide analyses we detected a transcriptional response to type I interferons (IFN-Is) in astrocytes during experimental CNS autoimmunity and also in CNS lesions from patients with multiple sclerosis (MS). IFN-I signaling in astrocytes reduces inflammation and experimental autoimmune encephalomyelitis (EAE) disease scores via the ligand-activated transcription factor aryl hydrocarbon receptor (AHR) and the suppressor of cytokine signaling 2 (SOCS2). The anti-inflammatory effects of nasally administered interferon (IFN)-ß are partly mediated by AHR. Dietary tryptophan is metabolized by the gut microbiota into AHR agonists that have an effect on astrocytes to limit CNS inflammation. EAE scores were increased following ampicillin treatment during the recovery phase, and CNS inflammation was reduced in antibiotic-treated mice by supplementation with the tryptophan metabolites indole, indoxyl-3-sulfate, indole-3-propionic acid and indole-3-aldehyde, or the bacterial enzyme tryptophanase. In individuals with MS, the circulating levels of AHR agonists were decreased. These findings suggest that IFN-Is produced in the CNS function in combination with metabolites derived from dietary tryptophan by the gut flora to activate AHR signaling in astrocytes and suppress CNS inflammation.

Multibow: digital spectral barcodes for cell tracing.

We introduce a multicolor labeling strategy (Multibow) for cell tracing experiments in developmental and regenerative processes. Building on Brainbow-based approaches that produce colors by differential expression levels of different fluorescent proteins, Multibow adds a layer of label diversity by introducing a binary code in which reporters are initially OFF and then probabilistically ON or OFF following Cre recombination. We have developed a library of constructs that contains seven different colors and three different subcellular localizations. Combining constructs from this library in the presence of Cre generates cells labeled with multiple independently expressed colors based on if each construct is ON or OFF following recombination. These labels form a unique "barcode" that allows the tracking of the cell and its clonal progenies in addition to expression level differences of each color. We tested Multibow in zebrafish which validates its design concept and suggests its utility for cell tracing applications in development and regeneration.

Association of Brain DNA methylation in SORL1, ABCA7, HLA-DRB5, SLC24A4, and BIN1 with pathological diagnosis of Alzheimer disease.

IMPORTANCE: Recent large-scale genome-wide association studies have discovered several genetic variants associated with Alzheimer disease (AD); however, the extent to which DNA methylation in these AD loci contributes to the disease susceptibility remains unknown. OBJECTIVE: To examine the association of brain DNA methylation in 28 reported AD loci with AD pathologies. DESIGN, SETTING, AND PARTICIPANTS: Ongoing community-based clinical pathological cohort studies of aging and dementia (the Religious Orders Study and the Rush Memory and Aging Project) among 740 autopsied participants 66.0 to 108.3 years old. EXPOSURES: DNA methylation levels at individual CpG sites generated from dorsolateral prefrontal cortex tissue using a bead assay. MAIN OUTCOMES AND MEASURES: Pathological diagnosis of AD by National Institute on Aging-Reagan criteria following a standard postmortem examination. RESULTS: Overall, 447 participants (60.4%) met the criteria for pathological diagnosis of AD. Brain DNA methylation in SORL1, ABCA7, HLA-DRB5, SLC24A4, and BIN1 was associated with pathological AD. The association was robustly retained after replacing the binary trait of pathological AD with 2 quantitative and molecular specific hallmarks of AD, namely, Aß load and paired helical filament tau tangle density. Furthermore, RNA expression of transcripts of SORL1 and ABCA7 was associated with paired helical filament tau tangle density, and the expression of BIN1 was associated with Aß load. CONCLUSIONS AND RELEVANCE: Brain DNA methylation in multiple AD loci is associated with AD pathologies. The results provide further evidence that disruption of DNA methylation is involved in the pathological process of AD.

Identification and characterization of latency-associated peptide-expressing γδ T cells.

γδ T cells are a subset of lymphocytes specialized in protecting the host against pathogens and tumours. Here we describe a subset of regulatory γδ T cells that express the latency-associated peptide (LAP), a membrane-bound TGF-ß1. Thymic CD27+IFN-γ+CCR9+α4ß7+TCRγδ+ cells migrate to the periphery, particularly to Peyer's patches and small intestine lamina propria, where they upregulate LAP, downregulate IFN-γ via ATF-3 expression and acquire a regulatory phenotype. TCRγδ+LAP+ cells express antigen presentation molecules and function as antigen presenting cells that induce CD4+Foxp3+ regulatory T cells, although TCRγδ+LAP+ cells do not themselves express Foxp3. Identification of TCRγδ+LAP+ regulatory cells provides an avenue for understanding immune regulation and biologic processes linked to intestinal function and disease.

Streaming level set algorithm for 3D segmentation of confocal microscopy images.

We present a high performance variant of the popular geodesic active contours which are used for splitting cell clusters in microscopy images. Previously, we implemented a linear pipelined version that incorporates as many cues as possible into developing a suitable level-set speed function so that an evolving contour exactly segments a cell/nuclei blob. We use image gradients, distance maps, multiple channel information and a shape model to drive the evolution. We also developed a dedicated seeding strategy that uses the spatial coherency of the data to generate an over complete set of seeds along with a quality metric which is further used to sort out which seed should be used for a given cell. However, the computational performance of any level-set methodology is quite poor when applied to thousands of 3D data-sets each containing thousands of cells. Those data-sets are common in confocal microscopy. In this work, we explore methods to stream the algorithm in shared memory, multi-core environments. By partitioning the input and output using spatial data structures we insure the spatial coherency needed by our seeding algorithm as well as improve drastically the speed without memory overhead. Our results show speed-ups up to a factor of six.

Detection of spatially correlated objects in 3D images using appearance models and coupled active contours.

We consider the problem of segmenting 3D images that contain a dense collection of spatially correlated objects, such as fluorescent labeled cells in tissue. Our approach involves an initial modeling phase followed by a data-fitting segmentation phase. In the first phase, cell shape (membrane bound) is modeled implicitly using a parametric distribution of correlation function estimates. The nucleus is modeled for its shape as well as image intensity distribution inspired from the physics of its image formation. In the second phase, we solve the segmentation problem using a variational level-set strategy with coupled active contours to minimize a novel energy functional. We demonstrate the utility of our approach on multispectral fluorescence microscopy images.

Turning on/off tumor-specific CTL response during progressive tumor growth.

Therapeutic vaccinations used to induce CTLs and treat firmly established tumors are generally ineffective. To understand the mechanisms underlying the failure of therapeutic vaccinations, we investigated the fate of tumor-specific CD8+ T cells in tumor-bearing mice with or without vaccinations. Our data demonstrate that tumor-specific CD8+ T cells are activated at the early stage of tumor growth, tumor-specific CTL response reaches a maximal level during progressive tumor growth, and tumor-specific CD8+ T cells lose cytolytic function at the late stage of tumor growth. The early stage therapeutic vaccination induces efficient antitumor activity by amplifying the CTL response, whereas the late-stage therapeutic vaccination is invalid due to tumor-induced dysfunction of CD8+ T cells. However, at the late stage, tumor-specific CD8+ T cells are still present in the periphery. These tumor-specific CD8+ T cells lose cytolytic activity, but retain IFN-gamma secretion function. In contrast to in vitro cultured tumor cells, in vivo growing tumor cells are more resistant to tumor-specific CTL killing, despite an increase of tumor Ag gene expression. Both tumor-induced CD8+ T cell dysfunction at the late stage and immune evasion developed by in vivo growing tumor cells contribute to an eventual inefficacy of therapeutic vaccinations. Our study suggests that it is important to design a vaccination regimen according to the stages of tumor growth and the functional states of tumor-specific CD8+ T cells.