Airway derived emphysema-specific alveolar type II cells exhibit impaired regenerative potential in COPD.

Emphysema, the progressive destruction of gas exchange surfaces in the lungs, is a hallmark of chronic obstructive pulmonary disease (COPD) that is presently incurable. This therapeutic gap is largely due to a poor understanding of potential drivers of impaired tissue regeneration, such as abnormal lung epithelial progenitor cells, including alveolar type II (ATII) and airway club cells. We discovered an emphysema-specific sub-population of ATII cells located in enlarged distal alveolar sacs, termed asATII cells. Single cell RNA-seq and in situ localisation revealed that asATII cells co-express the alveolar marker surfactant protein C (SPC) and the club cell marker secretaglobin-3A2 (SCGB3A2). A similar ATII sub-population derived from club cells was also identified in mouse COPD models using lineage labeling. Human and mouse ATII sub-populations formed 80-90% fewer alveolar organoids than healthy controls, indicating reduced progenitor function. Targeting asATII cells or their progenitor club cells could reveal novel COPD treatment strategies.

Multi-pronged analysis of pediatric low-grade glioma reveals a unique tumor microenvironment associated with BRAF alterations.

Pediatric low-grade gliomas (pLGG) comprise 35% of all brain tumors. Despite favorable survival, patients experience significant morbidity from disease and treatments. A deeper understanding of pLGG biology is essential to identify novel, more effective, and less toxic therapies. We utilized single cell RNA sequencing (scRNA-seq), spatial transcriptomics, and cytokine analyses to characterize and understand tumor and immune cell heterogeneity across pLGG. scRNA-seq revealed tumor and immune cells within the tumor microenvironment (TME). Tumor cell subsets revealed a developmental hierarchy with progenitor and mature cell populations. Immune cells included myeloid and lymphocytic cells. There was a significant difference between the prevalence of two major myeloid subclusters between pilocytic astrocytoma (PA) and ganglioglioma (GG). Bulk and single-cell cytokine analyses evaluated the immune cell signaling cascade with distinct immune phenotypes among tumor samples. KIAA1549-BRAF tumors appeared more immunogenic, secreting higher levels of immune cell activators and chemokines, compared to BRAF V600E tumors. Spatial transcriptomics revealed the differential gene expression of these chemokines and their location within the TME. A multi-pronged analysis of pLGG demonstrated the complexity of the pLGG TME and differences between genetic drivers that may influence their response to immunotherapy. Further investigation of immune cell infiltration and tumor-immune interactions is warranted. Key points: There is a developmental hierarchy in neoplastic population comprising of both progenitor-like and mature cell types in both PA and GG.A more immunogenic, immune activating myeloid population is present in PA compared to GG. Functional analysis and spatial transcriptomics show higher levels of immune mobilizing chemokines in KIAA1549-BRAF fusion PA tumor samples compared to BRAF V600E GG samples. Importance of the Study: While scRNA seq provides information on cellular heterogeneity within the tumor microenvironment (TME), it does not provide a complete picture of how these cells are interacting or where they are located. To expand on this, we used a three-pronged approach to better understand the biology of pediatric low-grade glioma (pLGG). By analyzing scRNA-seq, secreted cytokines and spatial orientation of cells within the TME, we strove to gain a more complete picture of the complex interplay between tumor and immune cells within pLGG. Our data revealed a complex heterogeneity in tumor and immune populations and identified an interesting difference in the immune phenotype among different subtypes.

Tumor and immune cell types interact to produce heterogeneous phenotypes of pediatric high-grade glioma.

BACKGROUND: Pediatric high-grade gliomas (PHGG) are aggressive brain tumors with 5-year survival rates ranging from <2% to 20% depending upon subtype. PHGG presents differently from patient to patient and is intratumorally heterogeneous, posing challenges in designing therapies. We hypothesized that heterogeneity occurs because PHGG comprises multiple distinct tumor and immune cell types in varying proportions, each of which may influence tumor characteristics. METHODS: We obtained 19 PHGG samples from our institution's pediatric brain tumor bank. We constructed a comprehensive transcriptomic dataset at the single-cell level using single-cell RNA-Seq (scRNA-Seq), identified known glial and immune cell types, and performed differential gene expression and gene set enrichment analysis. We conducted multi-channel immunofluorescence (IF) staining to confirm the transcriptomic results. RESULTS: Our PHGG samples included 3 principal predicted tumor cell types: astrocytes, oligodendrocyte progenitors (OPCs), and mesenchymal-like cells (Mes). These cell types differed in their gene expression profiles, pathway enrichment, and mesenchymal character. We identified a macrophage population enriched in mesenchymal and inflammatory gene expression as a possible source of mesenchymal tumor characteristics. We found evidence of T-cell exhaustion and suppression. CONCLUSIONS: PHGG comprises multiple distinct proliferating tumor cell types. Microglia-derived macrophages may drive mesenchymal gene expression in PHGG. The predicted Mes tumor cell population likely derives from OPCs. The variable tumor cell populations rely on different oncogenic pathways and are thus likely to vary in their responses to therapy.

Integration of single-nuclei RNA-sequencing, spatial transcriptomics and histochemistry defines the complex microenvironment of NF1-associated plexiform neurofibromas.

Plexiform neurofibroma (PN) is a leading cause of morbidity in children with the genetic condition Neurofibromatosis Type 1 (NF1), often disfiguring or threatening vital structures. During formation of PN, a complex tumor microenvironment (TME) develops, with recruitment of neoplastic and non-neoplastic cell types being critical for growth and progression. Due to the cohesive cellularity of PN, single-cell RNA-sequencing is difficult and may result in a loss of detection of critical cellular subpopulations. To bypass this barrier, we performed single-nuclei RNA-sequencing (snRNA-seq) on 8 frozen PN samples, and integrated this with spatial transcriptomics (ST) in 4 PN samples and immunohistochemistry to provide morphological context to transcriptomic data. SnRNA-seq analysis definitively charted the heterogeneous cellular subpopulations in the PN TME, with the predominant fraction being fibroblast subtypes. PN showed a remarkable amount of inter-sample homogeneity regarding cellular subpopulation proportions despite being resected from a variety of anatomical locations. ST analysis identified distinct cellular subpopulations which were annotated using snRNA-seq data and correlated with histological features. Schwann cell/fibroblast interactions were identified by receptor/ligand interaction analysis demonstrating a high probability of Neurexin 1/Neuroligin 1 (NRXN1/NLGN1) receptor-ligand cross-talk predicted between fibroblasts and non-myelinated Schwann cells (NM-SC) and subtypes, respectively. We observed aberrant expression of NRXN1 and NLGN1 in our PN snRNA-seq data compared to a normal mouse sciatic nerve single-cell RNA-seq dataset. This pathway has never been described in PN and may indicate a clear and direct communication pathway between putative NM-SC cells of origin and surrounding fibroblasts, potentially driving disease progression. SnRNA-seq integrated with spatial transcriptomics advances our understanding of the complex cellular heterogeneity of PN TME and identify potential novel communication pathways that may drive disease progression, a finding that could provide translational therapy options for patients with these devastating tumors of childhood and early adulthood.

Regulation of epithelial transitional states in murine and human pulmonary fibrosis.

Idiopathic pulmonary fibrosis (IPF) is a progressive scarring disease arising from impaired regeneration of the alveolar epithelium after injury. During regeneration, type 2 alveolar epithelial cells (AEC2s) assume a transitional state that upregulates multiple keratins and ultimately differentiate into AEC1s. In IPF, transitional AECs accumulate with ineffectual AEC1 differentiation. However, whether and how transitional cells cause fibrosis, whether keratins regulate transitional cell accumulation and fibrosis, and why transitional AECs and fibrosis resolve in mouse models but accumulate in IPF are unclear. Here, we show that human keratin 8 (KRT8) genetic variants were associated with IPF. Krt8-/- mice were protected from fibrosis and accumulation of the transitional state. Keratin 8 (K8) regulated the expression of macrophage chemokines and macrophage recruitment. Profibrotic macrophages and myofibroblasts promoted the accumulation of transitional AECs, establishing a K8-dependent positive feedback loop driving fibrogenesis. Finally, rare murine transitional AECs were highly senescent and basaloid and may not differentiate into AEC1s, recapitulating the aberrant basaloid state in human IPF. We conclude that transitional AECs induced and were maintained by fibrosis in a K8-dependent manner; in mice, most transitional cells and fibrosis resolved, whereas in human IPF, transitional AECs evolved into an aberrant basaloid state that persisted with progressive fibrosis.

Multi-omic approach identifies hypoxic tumor-associated myeloid cells that drive immunobiology of high-risk pediatric ependymoma.

Ependymoma (EPN) is a devastating childhood brain tumor. Single-cell analyses have illustrated the cellular heterogeneity of EPN tumors, identifying multiple neoplastic cell states including a mesenchymal-differentiated subpopulation which characterizes the PFA1 subtype. Here, we characterize the EPN immune environment, in the context of both tumor subtypes and tumor cell subpopulations using single-cell sequencing (scRNAseq, n = 27), deconvolution of bulk tumor gene expression (n = 299), spatial proteomics (n = 54), and single-cell cytokine release assays (n = 12). We identify eight distinct myeloid-derived subpopulations from which a group of cells, termed hypoxia myeloid cells, demonstrate features of myeloid-derived suppressor cells, including IL6/STAT3 pathway activation and wound healing ontologies. In PFA tumors, hypoxia myeloid cells colocalize with mesenchymal-differentiated cells in necrotic and perivascular niches and secrete IL-8, which we hypothesize amplifies the EPN immunosuppressive microenvironment. This myeloid cell-driven immunosuppression will need to be targeted for immunotherapy to be effective in this difficult-to-cure childhood brain tumor.

Significant increase of high-risk chromosome 1q gain and 6q loss at recurrence in posterior fossa group A ependymoma: A multicenter study.

BACKGROUND: Ependymoma (EPN) posterior fossa group A (PFA) has the highest rate of recurrence and the worst prognosis of all EPN molecular groups. At relapse, it is typically incurable even with re-resection and re-irradiation. The biology of recurrent PFA remains largely unknown; however, the increasing use of surgery at first recurrence has now provided access to clinical samples to facilitate a better understanding of this. METHODS: In this large longitudinal international multicenter study, we examined matched samples of primary and recurrent disease from PFA patients to investigate the biology of recurrence. RESULTS: DNA methylome derived copy number variants (CNVs) revealed large-scale chromosome gains and losses at recurrence in PFA. CNV changes were dominated by chromosome 1q gain and/or 6q loss, both previously identified as high-risk factors in PFA, which were present in 23% at presentation but increased to 61% at first recurrence. Multivariate survival analyses of this cohort showed that cases with 1q gain or 6q loss at first recurrence were significantly more likely to recur again. Predisposition to 1q+/6q- CNV changes at recurrence correlated with hypomethylation of heterochromatin-associated DNA at presentation. Cellular and molecular analyses revealed that 1q+/6q- PFA had significantly higher proportions of proliferative neuroepithelial undifferentiated progenitors and decreased differentiated neoplastic subpopulations. CONCLUSIONS: This study provides clinically and preclinically actionable insights into the biology of PFA recurrence. The hypomethylation predisposition signature in PFA is a potential risk-classifier for trial stratification. We show that the cellular heterogeneity of PFAs evolves largely because of genetic evolution of neoplastic cells.

Identification of HDAC4 as a potential therapeutic target and prognostic biomarker for ZFTA-fused ependymomas.

Members of the HDAC family are predictive biomarkers and regulate the tumorigenesis in several cancers. However, the role of these genes in the biology of intracranial ependymomas (EPNs) remains unexplored. Here, an analysis of eighteen HDACs genes in an EPN transcriptomic dataset, revealed significantly higher levels of HDAC4 in supratentorial ZFTA fusion (ST-ZFTA) compared with ST-YAP1 fusion and posterior fossa EPNs, while HDAC7 and SIRT2 were downregulated in ST-ZFTA. HDAC4 was also overexpressed in ST-ZFTA as measured by single-cell RNA-Seq, quantitative real time-polymerase chain reaction, and immunohistochemistry. Survival analyses showed a significantly worse outcome for EPNs with higher HDAC4 and SIRT1 mRNA levels. Ontology enrichment analysis showed an HDAC4-high signature consistent with viral processes while collagen-containing extracellular matrix and cell-cell junction were enriched in those with an HDAC4-low signature. Immune gene analysis demonstrated a correlation between HDAC4 expression and low levels of NK resting cells. Several small molecules compounds targeting HDAC4 and ABCG2, were predicted by in silico analysis to be effective against HDAC4-high ZFTA. Our results provide novel insights into the biology of the HDAC family in intracranial ependymomas and reveal HDAC4 as a prognostic marker and potential therapeutic target in ST-ZFTA.

Intron dynamics reveal principles of gene regulation during the maternal-to-zygotic transition.

The maternal-to-zygotic transition (MZT) is a conserved embryonic process in animals where developmental control shifts from the maternal to zygotic genome. A key step in this transition is zygotic transcription, and deciphering the MZT requires classifying newly transcribed genes. However, due to current technological limitations, this starting point remains a challenge for studying many species. Here, we present an alternative approach that characterizes transcriptome changes based solely on RNA-seq data. By combining intron-mapping reads and transcript-level quantification, we characterized transcriptome dynamics during the Drosophila melanogaster MZT. Our approach provides an accessible platform to investigate transcriptome dynamics that can be applied to the MZT in nonmodel organisms. In addition to classifying zygotically transcribed genes, our analysis revealed that over 300 genes express different maternal and zygotic transcript isoforms due to alternative splicing, polyadenylation, and promoter usage. The vast majority of these zygotic isoforms have the potential to be subject to different regulatory control, and over two-thirds encode different proteins. Thus, our analysis reveals an additional layer of regulation during the MZT, where new zygotic transcripts can generate additional proteome diversity.

Spatial transcriptomic analysis delineates epithelial and mesenchymal subpopulations and transition stages in childhood ependymoma.

BACKGROUND: The diverse cellular constituents of childhood brain tumor ependymoma, recently revealed by single cell RNA-sequencing, may underly therapeutic resistance. Here we use spatial transcriptomics to further advance our understanding of the tumor microenvironment, mapping cellular subpopulations to the tumor architecture of ependymoma posterior fossa subgroup A (PFA), the commonest and most deadly childhood ependymoma variant. METHODS: Spatial transcriptomics data from intact PFA sections was deconvoluted to resolve the histological arrangement of neoplastic and non-neoplastic cell types. Key findings were validated using immunohistochemistry, in vitro functional assays and outcome analysis in clinically-annotated PFA bulk transcriptomic data. RESULTS: PFA are comprised of epithelial and mesenchymal histological zones containing a diversity of cellular states, each zone including co-existing and spatially distinct undifferentiated progenitor-like cells; a quiescent mesenchymal zone population, and a second highly mitotic progenitor population that is restricted to hypercellular epithelial zones and that is more abundant in progressive tumors. We show that myeloid cell interaction is the leading cause of mesenchymal transition in PFA, occurring in zones spatially distinct from hypoxia-induced mesenchymal transition, and these distinct EMT-initiating processes were replicated using in vitro models of PFA. CONCLUSIONS: These insights demonstrate the utility of spatial transcriptomics to advance our understanding of ependymoma biology, revealing a clearer picture of the cellular constituents of PFA, their interactions and influence on tumor progression.

Single-cell transcriptome mapping identifies a local, innate B cell population driving chronic rejection after lung transplantation.

Bronchiolitis obliterans syndrome (BOS) is the main reason for poor outcomes after lung transplantation (LTx). We and others have recently identified B cells as major contributors to BOS after LTx. The extent of B cell heterogeneity and the relative contributions of B cell subpopulations to BOS, however, remain unclear. Here, we provide a comprehensive analysis of cell population changes and their gene expression patterns during chronic rejection after orthotopic LTx in mice. Of 11 major cell types, Mzb1-expressing plasma cells (PCs) were the most prominently increased population in BOS lungs. These findings were validated in 2 different cohorts of human BOS after LTx. A Bhlhe41, Cxcr3, and Itgb1 triple-positive B cell subset, also expressing classical markers of the innate-like B-1 B cell population, served as the progenitor pool for Mzb1+ PCs. This subset accounted for the increase in IgG2c production within BOS lung grafts. A genetic lack of Igs decreased BOS severity after LTx. In summary, we provide a detailed analysis of cell population changes during BOS. IgG+ PCs and their progenitors - an innate B cell subpopulation - are the major source of local Ab production and a significant contributor to BOS after LTx.

Fatal COVID-19 and Non-COVID-19 Acute Respiratory Distress Syndrome Is Associated with Incomplete Alveolar Type 1 Epithelial Cell Differentiation from the Transitional State without Fibrosis.

Acute respiratory distress syndrome (ARDS) due to coronavirus disease 2019 and other etiologies results from injury to the alveolar epithelial cell (AEC) barrier resulting in noncardiogenic pulmonary edema, which causes acute respiratory failure; recovery requires epithelial regeneration. During physiological regeneration in mice, type 2 AECs (AEC2s) proliferate, exit the cell cycle, transiently assume a transitional state, then differentiate into type 1 AECs (AEC1s); in humans, persistence of the transitional state is associated with pulmonary fibrosis. It is unknown whether transitional cells emerge and differentiate into AEC1s without fibrosis in human ARDS and why transitional cells differentiate into AEC1s during physiological regeneration but persist in fibrosis. We hypothesized that incomplete but ongoing AEC1 differentiation from transitional cells without fibrosis may underlie persistent barrier permeability and acute respiratory failure in ARDS. Immunostaining of postmortem ARDS lungs revealed abundant transitional cells without fibrosis. They were typically cuboidal or partially spread, sometimes flat, and occasionally expressed AEC1 markers. Immunostaining and/or single-cell RNA sequencing revealed that transitional cells in mouse models of physiological regeneration, ARDS, and fibrosis express markers of cell cycle exit but only in fibrosis express a specific senescence marker. Thus, in severe, fatal early ARDS, AEC1 differentiation from transitional cells is incomplete, underlying persistent barrier permeability and respiratory failure but ongoing without fibrosis; senescence of transitional cells may be associated with pulmonary fibrosis.

Neoplastic and immune single-cell transcriptomics define subgroup-specific intra-tumoral heterogeneity of childhood medulloblastoma.

BACKGROUND: Medulloblastoma (MB) is a heterogeneous disease in which neoplastic cells and associated immune cells contribute to disease progression. We aimed to determine the influence of neoplastic and immune cell diversity on MB biology in patient samples and animal models. METHODS: To better characterize cellular heterogeneity in MB we used single-cell RNA sequencing, immunohistochemistry, and deconvolution of transcriptomic data to profile neoplastic and immune populations in patient samples and animal models across childhood MB subgroups. RESULTS: Neoplastic cells cluster primarily according to individual sample of origin which is influenced by chromosomal copy number variance. Harmony alignment reveals novel MB subgroup/subtype-associated subpopulations that recapitulate neurodevelopmental processes, including photoreceptor and glutamatergic neuron-like cells in molecular subgroups GP3 and GP4, and a specific nodule-associated neuronally differentiated subpopulation in the sonic hedgehog subgroup. We definitively chart the spectrum of MB immune cell infiltrates, which include subpopulations that recapitulate developmentally related neuron-pruning and antigen-presenting myeloid cells. MB cellular diversity matching human samples is mirrored in subgroup-specific mouse models of MB. CONCLUSIONS: These findings provide a clearer understanding of the diverse neoplastic and immune cell subpopulations that constitute the MB microenvironment.

ENTPD3 Marks Mature Stem Cell-Derived ß-Cells Formed by Self-Aggregation In Vitro.

Stem cell-derived ß-like cells (sBC) carry the promise of providing an abundant source of insulin-producing cells for use in cell replacement therapy for patients with diabetes, potentially allowing widespread implementation of a practical cure. To achieve their clinical promise, sBC need to function comparably with mature adult ß-cells, but as yet they display varying degrees of maturity. Indeed, detailed knowledge of the events resulting in human ß-cell maturation remains obscure. Here we show that sBC spontaneously self-enrich into discreet islet-like cap structures within in vitro cultures, independent of exogenous maturation conditions. Multiple complementary assays demonstrate that this process is accompanied by functional maturation of the self-enriched sBC (seBC); however, the seBC still contain distinct subpopulations displaying different maturation levels. Interestingly, the surface protein ENTPD3 (also known as nucleoside triphosphate diphosphohydrolase-3 [NDPTase3]) is a specific marker of the most mature seBC population and can be used for mature seBC identification and sorting. Our results illuminate critical aspects of in vitro sBC maturation and provide important insights toward developing functionally mature sBC for diabetes cell replacement therapy.

Liver Transcriptome Dynamics During Hibernation Are Shaped by a Shifting Balance Between Transcription and RNA Stability.

Hibernators dramatically lower metabolism to save energy while fasting for months. Prolonged fasting challenges metabolic homeostasis, yet small-bodied hibernators emerge each spring ready to resume all aspects of active life, including immediate reproduction. The liver is the body's metabolic hub, processing and detoxifying macromolecules to provide essential fuels to brain, muscle and other organs throughout the body. Here we quantify changes in liver gene expression across several distinct physiological states of hibernation in 13-lined ground squirrels, using RNA-seq to measure the steady-state transcriptome and GRO-seq to measure transcription for the first time in a hibernator. Our data capture key timepoints in both the seasonal and torpor-arousal cycles of hibernation. Strong positive correlation between transcription and the transcriptome indicates that transcriptional control dominates the known seasonal reprogramming of metabolic gene expression in liver for hibernation. During the torpor-arousal cycle, however, discordance develops between transcription and the steady-state transcriptome by at least two mechanisms: 1) although not transcribed during torpor, some transcripts are unusually stable across the torpor bout; and 2) unexpectedly, on some genes, our data suggest continuing, slow elongation with a failure to terminate transcription across the torpor bout. While the steady-state RNAs corresponding to these read through transcripts did not increase during torpor, they did increase shortly after rewarming despite their simultaneously low transcription. Both of these mechanisms would assure the immediate availability of functional transcripts upon rewarming. Integration of transcriptional, post-transcriptional and RNA stability control mechanisms, all demonstrated in these data, likely initiate a serial gene expression program across the short euthermic period that restores the tissue and prepares the animal for the next bout of torpor.

RNA-binding proteins regulate aldosterone homeostasis in human steroidogenic cells.

Angiotensin II (AngII) stimulates adrenocortical cells to produce aldosterone, a master regulator of blood pressure. Despite extensive characterization of the transcriptional and enzymatic control of adrenocortical steroidogenesis, there are still major gaps in the precise regulation of AII-induced gene expression kinetics. Specifically, we do not know the regulatory contribution of RNA-binding proteins (RBPs) and RNA decay, which can control the timing of stimulus-induced gene expression. To investigate this question, we performed a high-resolution RNA-seq time course of the AngII stimulation response and 4-thiouridine pulse labeling in a steroidogenic human cell line (H295R). We identified twelve temporally distinct gene expression responses that contained mRNA encoding proteins known to be important for various steps of aldosterone production, such as cAMP signaling components and steroidogenic enzymes. AngII response kinetics for many of these mRNAs revealed a coordinated increase in both synthesis and decay. These findings were validated in primary human adrenocortical cells stimulated ex vivo with AngII. Using a candidate screen, we identified a subset of RNA-binding protein and RNA decay factors that activate or repress AngII-stimulated aldosterone production. Among the repressors of aldosterone were BTG2, which promotes deadenylation and global RNA decay. BTG2 was induced in response to AngII stimulation and promoted the repression of mRNAs encoding pro-steroidogenic factors indicating the existence of an incoherent feedforward loop controlling aldosterone homeostasis. These data support a model in which coordinated increases in transcription and decay facilitate the major transcriptomic changes required to implement a pro-steroidogenic expression program that actively resolved to prevent aldosterone overproduction.

Single-Cell RNA Sequencing Reveals a Unique Monocyte Population in Bronchoalveolar Lavage Cells of Mice Challenged With Afghanistan Particulate Matter and Allergen.

Upon returning from deployment to Afghanistan, a substantial number of U.S. military personnel report deployment-related lung disease (DRLD) symptoms, including those consistent with an asthma-like airways disease. DRLD is thought to be caused by prolonged inhalation of toxic desert particulate matter, which can persist in the postdeployment setting such as exposure to common household allergens. The goal of this study was to define the transcriptomic responses of lung leukocytes of mice exposed to Afghanistan desert particulate matter (APM) and house dust mite (HDM). C57BL/6 mice (n = 15/group) were exposed to filtered air or aerosolized APM for 12 days, followed by intranasal PBS or HDM allergen challenges for 24 h. Bronchoalveolar lavage (BAL) cells were collected for single-cell RNA sequencing (scRNAseq), and assessment of inflammation and airway hyper-responsiveness. Unsupervised clustering of BAL cell scRNAseq data revealed a unique monocyte population induced only by both APM and allergen treatments. This population of monocytes is characterized by the expression of genes involved in allergic asthma, including Alox15. We validated Alox15 expression in monocytes via immunostaining of lung tissue. APM pre-exposure, followed by the HDM challenge, led to significantly increased total respiratory system resistance compared with filtered air controls. Using this mouse model to mimic DRLD, we demonstrated that inhalation of airborne PM during deployment may prime airways to be more responsive to allergen exposure after returning home, which may be linked to dysregulated immune responses such as induction of a unique lung monocyte population.

Cell-level metadata are indispensable for documenting single-cell sequencing datasets.

Single-cell RNA sequencing (scRNA-seq) provides an unprecedented view of cellular diversity of biological systems. However, across the thousands of publications and datasets generated using this technology, we estimate that only a minority (<25%) of studies provide cell-level metadata information containing identified cell types and related findings of the published dataset. Metadata omission hinders reproduction, exploration, validation, and knowledge transfer and is a common problem across journals, data repositories, and publication dates. We encourage investigators, reviewers, journals, and data repositories to improve their standards and ensure proper documentation of these valuable datasets.

Alternative RNA structures formed during transcription depend on elongation rate and modify RNA processing.

Most RNA processing occurs co-transcriptionally. We interrogated nascent pol II transcripts by chemical and enzymatic probing and determined how the "nascent RNA structureome" relates to splicing, A-I editing and transcription speed. RNA folding within introns and steep structural transitions at splice sites are associated with efficient co-transcriptional splicing. A slow pol II mutant elicits extensive remodeling into more folded conformations with increased A-I editing. Introns that become more structured at their 3' splice sites get co-transcriptionally excised more efficiently. Slow pol II altered folding of intronic Alu elements where cryptic splicing and intron retention are stimulated, an outcome mimicked by UV, which decelerates transcription. Slow transcription also remodeled RNA folding around alternative exons in distinct ways that predict whether skipping or inclusion is favored, even though it occurs post-transcriptionally. Hence, co-transcriptional RNA folding modulates post-transcriptional alternative splicing. In summary, the plasticity of nascent transcripts has widespread effects on RNA processing.

Fatal COVID-19 and non-COVID-19 Acute Respiratory Distress Syndrome is Associated with Incomplete Alveolar Type 1 Epithelial Cell Differentiation from the Transitional State Without Fibrosis.

ARDS due to COVID-19 and other etiologies results from injury to the alveolar epithelial cell (AEC) barrier resulting in noncardiogenic pulmonary edema, which causes acute respiratory failure; clinical recovery requires epithelial regeneration. During physiologic regeneration in mice, AEC2s proliferate, exit the cell cycle, and transiently assume a transitional state before differentiating into AEC1s; persistence of the transitional state is associated with pulmonary fibrosis in humans. It is unknown whether transitional cells emerge and differentiate into AEC1s without fibrosis in human ARDS and why transitional cells differentiate into AEC1s during physiologic regeneration but persist in fibrosis. We hypothesized that incomplete but ongoing AEC1 differentiation from transitional cells without fibrosis may underlie persistent barrier permeability and fatal acute respiratory failure in ARDS. Immunostaining of postmortem ARDS lungs revealed abundant transitional cells in organized monolayers on alveolar septa without fibrosis. They were typically cuboidal or partially spread, sometimes flat, and occasionally expressed AEC1 markers. Immunostaining and/or interrogation of scRNAseq datasets revealed that transitional cells in mouse models of physiologic regeneration, ARDS, and fibrosis express markers of cell cycle exit but only in fibrosis express a specific senescence marker. Thus, in severe, fatal early ARDS, AEC1 differentiation from transitional cells is incomplete, underlying persistent barrier permeability and respiratory failure, but ongoing without fibrosis; senescence of transitional cells may be associated with pulmonary fibrosis.