Neonatal hyperoxia induces activated pulmonary cellular states and sex-dependent transcriptomic changes in a model of experimental bronchopulmonary dysplasia.

Hyperoxia disrupts lung development in mice and causes bronchopulmonary dysplasia (BPD) in neonates. To investigate sex-dependent molecular and cellular programming involved in hyperoxia, we surveyed the mouse lung using single cell RNA sequencing (scRNA-seq), and validated our findings in human neonatal lung cells in vitro. Hyperoxia-induced inflammation in alveolar type (AT) 2 cells gave rise to damage-associated transient progenitors (DATPs). It also induced a new subpopulation of AT1 cells with reduced expression of growth factors normally secreted by AT1 cells, but increased mitochondrial gene expression. Female alveolar epithelial cells had less EMT and pulmonary fibrosis signaling in hyperoxia. In the endothelium, expansion of Car4+ EC (Cap2) was seen in hyperoxia along with an emergent subpopulation of Cap2 with repressed VEGF signaling. This regenerative response was increased in females exposed to hyperoxia. Mesenchymal cells had inflammatory signatures in hyperoxia, with a new distal interstitial fibroblast subcluster characterized by repressed lipid biosynthesis and a transcriptomic signature resembling myofibroblasts. Hyperoxia-induced gene expression signatures in human neonatal fibroblasts and alveolar epithelial cells in vitro resembled mouse scRNA-seq data. These findings suggest that neonatal exposure to hyperoxia programs distinct sex-specific stem cell progenitor and cellular reparative responses that underpin lung remodeling in BPD.

Single-cell analysis of cell fate bifurcation in the chordate Ciona.

BACKGROUND: Inductive signaling interactions between different cell types are a major mechanism for the further diversification of embryonic cell fates. Most blastomeres in the model chordate Ciona robusta become restricted to a single predominant fate between the 64-cell and mid-gastrula stages. The deeply stereotyped and well-characterized Ciona embryonic cell lineages allow the transcriptomic analysis of newly established cell types very early in their divergence from sibling cell states without the pseudotime inference needed in the analysis of less synchronized cell populations. This is the first ascidian study to use droplet scRNAseq with large numbers of analyzed cells as early as the 64-cell stage when major lineages such as primary notochord first become fate restricted. RESULTS AND CONCLUSIONS: We identify 59 distinct cell states, including new subregions of the b-line neural lineage and the early induction of the tail tip epidermis. We find that 34 of these cell states are directly or indirectly dependent on MAPK-mediated signaling critical to early Ciona patterning. Most of the MAPK-dependent bifurcations are canalized with the signal-induced cell fate lost upon MAPK inhibition, but the posterior endoderm is unique in being transformed into a novel state expressing some but not all markers of both endoderm and muscle. Divergent gene expression between newly bifurcated sibling cell types is dominated by upregulation in the induced cell type. The Ets family transcription factor Elk1/3/4 is uniquely upregulated in nearly all the putatively direct inductions. Elk1/3/4 upregulation together with Ets transcription factor binding site enrichment analysis enables inferences about which bifurcations are directly versus indirectly controlled by MAPK signaling. We examine notochord induction in detail and find that the transition between a Zic/Ets-mediated regulatory state and a Brachyury/FoxA-mediated regulatory state is unexpectedly late. This supports a "broad-hourglass" model of cell fate specification in which many early tissue-specific genes are induced in parallel to key tissue-specific transcriptional regulators via the same set of transcriptional inputs.

Immune cell residency in the nasal mucosa may partially explain respiratory disease severity across the age range.

Previous studies focusing on the age disparity in COVID-19 severity have suggested that younger individuals mount a more robust innate immune response in the nasal mucosa after infection with SARS-CoV-2. However, it is unclear if this reflects increased immune activation or increased immune residence in the nasal mucosa. We hypothesized that immune residency in the nasal mucosa of healthy individuals may differ across the age range. We applied single-cell RNA-sequencing and measured the cellular composition and transcriptional profile of the nasal mucosa in 35 SARS-CoV-2 negative children and adults, ranging in age from 4 months to 65 years. We analyzed in total of ~ 30,000 immune and epithelial cells and found that age and immune cell proportion in the nasal mucosa are inversely correlated, with little evidence for structural changes in the transcriptional state of a given cell type across the age range. Orthogonal validation by epigenome sequencing indicate that it is especially cells of the innate immune system that underlie the age-association. Additionally, we characterize the predominate immune cell type in the nasal mucosa: a resident T cell like population with potent antiviral properties. These results demonstrate fundamental changes in the immune cell makeup of the uninfected nasal mucosa over the lifespan. The resource we generate here is an asset for future studies focusing on respiratory infection and immunization strategies.

High-resolution epigenome analysis in nasal samples derived from children with respiratory viral infections reveals striking changes upon SARS-CoV-2 infection.

Background: DNA methylation patterns of the human genome can be modified by environmental stimuli and provide dense information on gene regulatory circuitries. We studied genome-wide DNA methylation in nasal samples from infants (<6 months) applying whole-genome bisulfite sequencing (WGBS) to characterize epigenome response to 10 different respiratory viral infections including SARS-CoV-2. Results: We identified virus-specific differentially methylated regions (vDMR) with human metapneumovirus (hMPV) and SARS-CoV-2 followed by Influenza B (Flu B) causing the weakest vs. strongest epigenome response with 496 vs. 78541 and 14361 vDMR, respectively. We found a strong replication rate of FluB (52%) and SARS-CoV-2 (42%) vDMR in independent samples indicating robust epigenome perturbation upon infection. Among the FluB and SARS-CoV-2 vDMRs, around 70% were hypomethylated and significantly enriched among epithelial cell-specific regulatory elements whereas the hypermethylated vDMRs for these viruses mapped more frequently to immune cell regulatory elements, especially those of the myeloid lineage. The hypermethylated vDMRs were also enriched among genes and genetic loci in monocyte activation pathways and monocyte count. Finally, we perform single-cell RNA-sequencing characterization of nasal mucosa in response to these two viruses to functionally analyze the epigenome perturbations. Which supports the trends we identified in methylation data and highlights and important role for monocytes. Conclusions: All together, we find evidence indicating genetic predisposition to innate immune response upon a respiratory viral infection. Our genome-wide monitoring of infant viral response provides first catalogue of associated host regulatory elements. Assessing epigenetic variation in individual patients may reveal evidence for viral triggers of childhood disease.

Brachyury controls Ciona notochord fate as part of a feed-forward network.

The notochord is a defining feature of the chordates. The transcription factor Brachyury (Bra) is a key regulator of notochord fate but here we show that it is not a unitary master regulator in the model chordate Ciona Ectopic Bra expression only partially reprograms other cell types to a notochord-like transcriptional profile and a subset of notochord-enriched genes is unaffected by CRISPR Bra disruption. We identify Foxa.a and Mnx as potential co-regulators, and find that combinatorial cocktails are more effective at reprogramming other cell types than Bra alone. We reassess the network relationships between Bra, Foxa.a and other components of the notochord gene regulatory network, and find that Foxa.a expression in the notochord is regulated by vegetal FGF signaling. It is a direct activator of Bra expression and has a binding motif that is significantly enriched in the regulatory regions of notochord-enriched genes. These and other results indicate that Bra and Foxa.a act together in a regulatory network dominated by positive feed-forward interactions, with neither being a classically defined master regulator.

Tunicate gastrulation.

Tunicates are a diverse group of invertebrate marine chordates that includes the larvaceans, thaliaceans, and ascidians. Because of their unique evolutionary position as the sister group of the vertebrates, tunicates are invaluable as a comparative model and hold the promise of revealing both conserved and derived features of chordate gastrulation. Descriptive studies in a broad range of tunicates have revealed several important unifying traits that make them unique among the chordates, including invariant cell lineages through gastrula stages and an overall morphological simplicity. Gastrulation has only been studied in detail in ascidians such as Ciona and Phallusia, where it involves a simple cup-shaped gastrula driven primarily by endoderm invagination. This appears to differ significantly from vertebrate models, such as Xenopus, in which mesoderm convergent extension and epidermal epiboly are major contributors to involution. These differences may reflect the cellular simplicity of the ascidian embryo.

Iterative and Complex Asymmetric Divisions Control Cell Volume Differences in Ciona Notochord Tapering.

The notochord of the invertebrate chordate Ciona forms a tapered rod at tailbud stages consisting of only 40 cylindrical cells in a single-file column. This tapered shape involves differences in notochord cell volume along the anterior-posterior axis. Here, we quantify sibling cell volume asymmetry throughout the developing notochord and find that there are distinctive patterns of unequal cleavage in all 4 bilateral pairs of A-line primary notochord founder cells and also in the B-line-derived secondary notochord founder cells. A quantitative model confirms that the observed patterns of unequal cleavage are sufficient to explain all the anterior-posterior variation in notochord cell volume. Many examples are known of cells that divide asymmetrically to give daughter cells of different size and fate. Here, by contrast, a series of subtle but iterative and finely patterned asymmetric divisions controls the shape of an entire organ. Quantitative 3D analysis of cell shape and spindle positioning allows us to infer multiple cellular mechanisms driving these unequal cleavages, including polarized displacements of the mitotic spindle, contributions from the shape of the mother cell, and late changes occurring between anaphase and abscission that potentially involve differential cortical contractility. We infer differential use of these mechanisms between different notochord blastomeres and also between different rounds of cell division. These results demonstrate a new role for asymmetric division in directly shaping a developing organ and point toward complex underlying mechanisms.

A temperature-adjusted developmental timer for precise embryonic staging.

Developmental biology research depends on careful staging of developing embryos, but the rate of development is extremely sensitive to the temperature at which embryos are raised. It is not always practical to grow embryos at a precisely controlled temperature, so here we describe a simple, inexpensive device based on an Arduino-compatible microprocessor and temperature sensor that provides a metric of developmental time that compensates for changes in temperature. The underlying assumption is that the rate of development will be linear with respect to temperature over an organism's thermal tolerance range. The device measures the ambient temperature and integrates effective degree-minutes over time. For convenience, this is displayed to the user as a temperature-adjusted standard developmental time. In initial testing we have found the device to be extremely helpful for fixing Ciona embryos during precise developmental windows.