Search

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) novel coronavirus 2019 (COVID-19) global pandemic has led to millions of cases and hundreds of thousands of deaths. While older adults appear at high risk for severe disease, hospitalizations and deaths due to SARS-CoV-2 among children have been relatively rare. Integrating single-cell RNA sequencing (scRNA-seq) of developing mouse lung with temporally resolved immunofluorescence in mouse and human lung tissue, we found that expression of SARS-CoV-2 Spike protein primer TMPRSS2 was highest in ciliated cells and type I alveolar epithelial cells (AT1), and TMPRSS2 expression increased with aging in mice and humans. Analysis of autopsy tissue from fatal COVID-19 cases detected SARS-CoV-2 RNA most frequently in ciliated and secretory cells in airway epithelium and AT1 cells in peripheral lung. SARS-CoV-2 RNA was highly colocalized in cells expressing TMPRSS2. Together, these data demonstrate the cellular spectrum infected by SARS-CoV-2 in lung epithelium and suggest that developmental regulation of TMPRSS2 may underlie the relative protection of infants and children from severe respiratory illness.

Pulmonary fibrosis (PF) is a form of chronic lung disease characterized by pathologic epithelial remodeling and accumulation of extracellular matrix (ECM). To comprehensively define the cell types, mechanisms, and mediators driving fibrotic remodeling in lungs with PF, we performed single-cell RNA sequencing of single-cell suspensions from 10 nonfibrotic control and 20 PF lungs. Analysis of 114,396 cells identified 31 distinct cell subsets/states. We report that a remarkable shift in epithelial cell phenotypes occurs in the peripheral lung in PF and identify several previously unrecognized epithelial cell phenotypes, including a KRT5 - / KRT17 + pathologic, ECM-producing epithelial cell population that was highly enriched in PF lungs. Multiple fibroblast subtypes were observed to contribute to ECM expansion in a spatially discrete manner. Together, these data provide high-resolution insights into the complexity and plasticity of the distal lung epithelium in human disease and indicate a diversity of epithelial and mesenchymal cells contribute to pathologic lung fibrosis., (Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).)

Insulin resistance and right ventricular (RV) dysfunction are associated with lipotoxicity in heritable forms of pulmonary arterial hypertension (PAH), commonly due to mutations in bone morphogenetic protein receptor type 2 (BMPR2). How BMPR2 dysfunction in cardiomyocytes alters glucose metabolism and the response of these cells to insulin are unknown. We hypothesized that BMPR2 mutation in cardiomyocytes alters glucose-supported mitochondrial respiration and impairs cellular responses to insulin, including glucose and lipid uptake. We performed metabolic assays, immunofluorescence and Western analysis, RNA profiling, and radioactive isotope uptake studies in H9c2 cardiomyocyte cell lines with and without patient-derived BMPR2 mutations (mutant cells), with and without insulin. Unlike control cells, BMPR2 mutant cardiomyocytes have reduced metabolic plasticity as indicated by reduced mitochondrial respiration with increased mitochondrial superoxide production. These mutant cells show enhanced baseline phosphorylation of insulin-signaling protein as indicated by increased Akt, AMPK, and acetyl-CoA carboxylase phosphorylation that may negatively influence fatty acid oxidation and enhance lipid uptake, and are insulin insensitive. Furthermore, mutant cells demonstrate an increase in milk fat globule-EGF factor-8 protein (MFGE8), which influences the insulin-signaling pathway by phosphorylating Akt Ser473 via phosphatidylinositol 3-kinase and mammalian target of rapamycin. In conclusion, BMPR2 mutant cardiomyocytes have reduced metabolic plasticity and fail to respond to glucose. These cells have enhanced baseline insulin-signaling pattern favoring insulin resistance with failure to augment this pattern in response to insulin. BMPR2 mutation possibly blunts glucose uptake and enhances lipid uptake in these cardiomyocytes. The MFGE8-driven signaling pathway may suggest a new mechanism underlying RV lipotoxicity in PAH.

Rationale: Idiopathic pulmonary arterial hypertension (IPAH) is usually without an identified genetic cause, despite clinical and molecular similarity to bone morphogenetic protein receptor type 2 mutation-associated heritable pulmonary arterial hypertension (PAH). There is phenotypic heterogeneity in IPAH, with a minority of patients showing long-term improvement with calcium channel-blocker therapy., Objectives: We sought to identify gene variants (GVs) underlying IPAH and determine whether GVs differ in vasodilator-responsive IPAH (VR-PAH) versus vasodilator-nonresponsive IPAH (VN-PAH)., Methods: We performed whole-exome sequencing (WES) on 36 patients with IPAH: 17 with VR-PAH and 19 with VN-PAH. Wnt pathway differences were explored in human lung fibroblasts., Measurements and Main Results: We identified 1,369 genes with 1,580 variants unique to IPAH. We used a gene ontology approach to analyze variants and identified overrepresentation of several pathways, including cytoskeletal function and ion binding. By mapping WES data to prior genome-wide association study data, Wnt pathway genes were highlighted. Using the connectivity map to define genetic differences between VR-PAH and VN-PAH, we found enrichment in vascular smooth muscle cell contraction pathways and greater genetic variation in VR-PAH versus VN-PAH. Using human lung fibroblasts, we found increased stimulated Wnt activity in IPAH versus controls., Conclusions: A pathway-based analysis of WES data in IPAH demonstrated multiple rare GVs that converge on key biological pathways, such as cytoskeletal function and Wnt signaling pathway. Vascular smooth muscle contraction-related genes were enriched in VR-PAH, suggesting a potentially different genetic predisposition for VR-PAH. This pathway-based approach may be applied to next-generation sequencing data in other diseases to uncover the contribution of unexpected or multiple GVs to a phenotype.

Rapid, ancient radiations pose one of the most difficult challenges for phylogenetic estimation. We used DNA sequence data of 9,006 aligned base pairs from five genes (chloroplast atpB, matK, rbcL, and 18S and 26S nrDNA) to elucidate relationships among major lineages of Saxifragales (angiosperms, eudicots). These relationships were poorly supported in previous studies, apparently because the lineages originated in rapid succession. Using an array of methods that explicitly incorporate assumptions about evolutionary process (weighted maximum parsimony, maximum likelihood, LogDet/paralinear transformed distances), we show that the initial diversification of Saxifragales was indeed rapid. We suggest that the poor resolution of our best phylogenetic estimate is not due to violations of assumptions or to combining data partitions having conflicting histories or processes. We show that estimated branch lengths during the initial diversification are exceedingly short, and we estimate that acquiring sufficient sequence data to resolve these relationships would require an extraordinary effort (approximately 10(7) bp), assuming a linear increase in branch support with branch length. However, our simulation of much larger data sets containing a distribution of phylogenetic signal similar to that of the five sampled gene sequences suggests a limit to achievable branch support. Using statistical tests of differences in the likelihoods of topologies, we evaluated whether the initial radiation of Saxifragales involved the simultaneous origin of major lineages. Our results are consistent with predictions that resolving the branching order of rapid, ancient radiations requires sampling characters that evolved rapidly at the time of the radiation but have since experienced a slower evolutionary rate.

Chrysosplenium (Saxifragaceae) consists of 57 species widely distributed in temperate and arctic regions of the Northern Hemisphere, with two species restricted to the southern part of South America. Species relationships within the genus are highly problematic. The genus has traditionally been divided into two groups, sometimes recognized as sections (Oppositifolia and Alternifolia), based on leaf arrangement, or, alternatively, into 17 series. Based on morphological features, Hara suggested that the genus originated in South America and then subsequently migrated to the Northern Hemisphere. We conducted phylogenetic analyses of DNA sequences of the chloroplast gene matK for species of Chrysosplenium to elucidate relationships, test Hara's biogeographic hypothesis for the genus, and examine chromosomal and gynoecial diversification. These analyses revealed that both sections Oppositifolia and Alternifolia are monophyletic and form two large sister clades. Hence, leaf arrangement is a good indicator of relationships within this genus. Hara's series Pilosa and Macrostemon are each also monophyletic; however, series Oppositifolia, Alternifolia, and Nepalensia are clearly not monophyletic. MacClade reconstructions suggest that the genus arose in Eastern Asia, rather than in South America, with several independent migration events from Asia to the New World. In one well-defined subclade, species from eastern and western North America form a discrete clade, with Old World species as their sister group, suggesting that the eastern and western North American taxa diverged following migration to that continent. The South American species forms a clade with species from eastern Asia; this disjunction may be the result of ancient long-distance dispersal. Character mapping demonstrated that gynoecial diversification is dynamic, with reversals from inferior to half-inferior ovaries, as well as to ovaries that appear superior. Chromosomal evolution also appears to be labile with several independent origins of n = 12 (from an original number of n = 11) and multiple episodes of aneuploidy.