A novel non-recurrent CNV deletion involving TBX4 and leaving TBX2 intact causes congenital alveolar dysplasia.

Congenital alveolar dysplasia (CAD) belongs to rare lethal lung developmental disorders (LLDDs) in neonates, manifesting with acute respiratory failure and pulmonary arterial hypertension refractory to treatment. The majority of CAD cases have been associated with copy-number variant (CNV) deletions at 17q23.1q23.2 or 5p12. Most CNV deletions at 17q23.1q23.2 were recurrent and encompassed two closely located genes, TBX4 and TBX2. In a few CAD cases, intragenic frameshifting deletions or single-nucleotide variants (SNVs) involved TBX4 but not TBX2. Here, we describe a male neonate who died at 27 days of life from acute respiratory failure caused by lung growth arrest along the spectrum of CAD confirmed by histopathological assessment. Trio-based genome sequencing revealed in the proband a novel non-recurrent ~1.07 Mb heterozygous CNV deletion at 17q23.2, encompassing TBX4 that arose de novo on the paternal chromosome. This is the first report of a larger-sized CNV deletion in a CAD patient involving TBX4 and leaving TBX2 intact. Our results, together with previous reports, indicate that perturbations of TBX4, rather than TBX2, cause severe lung phenotypes in humans.

Single Cell Multiomics Identifies Cells and Genetic Networks Underlying Alveolar Capillary Dysplasia.

Rationale: Alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV) is a lethal developmental disorder of lung morphogenesis caused by insufficiency of FOXF1 (forkhead box F1) transcription factor function. The cellular and transcriptional mechanisms by which FOXF1 deficiency disrupts human lung formation are unknown. Objectives: To identify cell types, gene networks, and cell-cell interactions underlying the pathogenesis of ACDMPV. Methods: We used single-nucleus RNA and assay for transposase-accessible chromatin sequencing, immunofluorescence confocal microscopy, and RNA in situ hybridization to identify cell types and molecular networks influenced by FOXF1 in ACDMPV lungs. Measurements and Main Results: Pathogenic single-nucleotide variants and copy-number variant deletions involving the FOXF1 gene locus in all subjects with ACDMPV (n = 6) were accompanied by marked changes in lung structure, including deficient alveolar development and a paucity of pulmonary microvasculature. Single-nucleus RNA and assay for transposase-accessible chromatin sequencing identified alterations in cell number and gene expression in endothelial cells (ECs), pericytes, fibroblasts, and epithelial cells in ACDMPV lungs. Distinct cell-autonomous roles for FOXF1 in capillary ECs and pericytes were identified. Pathogenic variants involving the FOXF1 gene locus disrupt gene expression in EC progenitors, inhibiting the differentiation or survival of capillary 2 ECs and cell-cell interactions necessary for both pulmonary vasculogenesis and alveolar type 1 cell differentiation. Loss of the pulmonary microvasculature was associated with increased VEGFA (vascular endothelial growth factor A) signaling and marked expansion of systemic bronchial ECs expressing COL15A1 (collagen type XV α 1 chain). Conclusions: Distinct FOXF1 gene regulatory networks were identified in subsets of pulmonary endothelial and fibroblast progenitors, providing both cellular and molecular targets for the development of therapies for ACDMPV and other diffuse lung diseases of infancy.

Proteomic Analysis of Human Lung Development.

Rationale: The current understanding of human lung development derives mostly from animal studies. Although transcript-level studies have analyzed human donor tissue to identify genes expressed during normal human lung development, protein-level analysis that would enable the generation of new hypotheses on the processes involved in pulmonary development are lacking. Objectives: To define the temporal dynamic of protein expression during human lung development. Methods: We performed proteomics analysis of human lungs at 10 distinct times from birth to 8 years to identify the molecular networks mediating postnatal lung maturation. Measurements and Main Results: We identified 8,938 proteins providing a comprehensive view of the developing human lung proteome. The analysis of the data supports the existence of distinct molecular substages of alveolar development and predicted the age of independent human lung samples, and extensive remodeling of the lung proteome occurred during postnatal development. Evidence of post-transcriptional control was identified in early postnatal development. An extensive extracellular matrix remodeling was supported by changes in the proteome during alveologenesis. The concept of maturation of the immune system as an inherent part of normal lung development was substantiated by flow cytometry and transcriptomics. Conclusions: This study provides the first in-depth characterization of the human lung proteome during development, providing a unique proteomic resource freely accessible at Lungmap.net. The data support the extensive remodeling of the lung proteome during development, the existence of molecular substages of alveologenesis, and evidence of post-transcriptional control in early postnatal development.

LungMAP Portal Ecosystem: Systems-Level Exploration of the Lung.

An improved understanding of the human lung necessitates advanced systems models informed by an ever-increasing repertoire of molecular omics, cellular, imaging, and pathological datasets. To centralize and standardize information across broad lung research efforts we expanded the LungMAP.net website into a new gateway portal. This portal connects a broad spectrum of research networks, bulk and single-cell multi-omics data and a diverse collection of image data that span mammalian lung development, and disease. The data are standardized across species and technologies using harmonized data and metadata models that leverage recent advances including those from the Human Cell Atlas, diverse ontologies, and the LungMAP CellCards initiative. To cultivate future discoveries, we have aggregated a diverse collection of single-cell atlases for multiple species (human, rhesus, mouse), to enable consistent queries across technologies, cohorts, age, disease, and drug treatment. These atlases are provided as independent and integrated queryable datasets, with an emphasis on dynamic visualization, figure generation, re-analysis, cell-type curation, and automated reference-based classification of user-provided single-cell genomics datasets (Azimuth). As this resource grows, we intend to increase the breadth of available interactive interfaces, supported data types, data portals and datasets from LungMAP and external research efforts.

Heterogeneity in Human Induced Pluripotent Stem Cell-derived Alveolar Epithelial Type II Cells Revealed with ABCA3/SFTPC Reporters.

Alveolar epithelial type 2 cells (AEC2s), the facultative progenitors of lung alveoli, are typically identified through the use of the canonical markers, SFTPC and ABCA3. Self-renewing AEC2-like cells have been generated from human induced pluripotent stem cells (iPSCs) through the use of knock-in SFTPC fluorochrome reporters. However, developmentally, SFTPC expression onset begins in the fetal distal lung bud tip and thus is not specific to mature AEC2s. Furthermore, SFTPC reporters appear to identify only those iPSC-derived AEC2s (iAEC2s) expressing the highest SFTPC levels. Here, we generate an ABCA3 knock-in GFP fusion reporter (ABCA3:GFP) that enables the purification of iAEC2s while allowing visualization of lamellar bodies, organelles associated with AEC2 maturation. Using an SFTPCtdTomato and ABCA3:GFP bifluorescent line for in vitro distal lung-directed differentiation, we observe later onset of ABCA3:GFP expression and broader identification of the subsequently emerging iAEC2 population based on ABCA3:GFP expression compared with SFTPCtdTomato expression. Comparing ABCA3:GFP/SFTPCtdTomato double-positive with ABCA3:GFP single-positive (SP) cells by RNA sequencing and functional studies reveals iAEC2 cellular heterogeneity with both populations functionally processing surfactant proteins but the SP cells exhibiting faster growth kinetics, increased clonogenicity, increased expression of progenitor markers, lower levels of SFTPC expression, and lower levels of AEC2 maturation markers. Over time, we observe that each population (double-positive and SP) gives rise to the other and each can serve as the parents of indefinitely self-renewing iAEC2 progeny. Our results indicate that iAEC2s are a heterogeneous population of cells with differing proliferation versus maturation properties, the majority of which can be tracked and purified using the ABCA3:GFP reporter or surrogate cell surface proteins, such as SLC34A2 and CPM.

Reconstructing differentiation networks and their regulation from time series single-cell expression data.

Generating detailed and accurate organogenesis models using single-cell RNA-seq data remains a major challenge. Current methods have relied primarily on the assumption that descendant cells are similar to their parents in terms of gene expression levels. These assumptions do not always hold for in vivo studies, which often include infrequently sampled, unsynchronized, and diverse cell populations. Thus, additional information may be needed to determine the correct ordering and branching of progenitor cells and the set of transcription factors (TFs) that are active during advancing stages of organogenesis. To enable such modeling, we have developed a method that learns a probabilistic model that integrates expression similarity with regulatory information to reconstruct the dynamic developmental cell trajectories. When applied to mouse lung developmental data, the method accurately distinguished different cell types and lineages. Existing and new experimental data validated the ability of the method to identify key regulators of cell fate.

Glucocorticoid regulates mesenchymal cell differentiation required for perinatal lung morphogenesis and function.

While antenatal glucocorticoids are widely used to enhance lung function in preterm infants, cellular and molecular mechanisms by which glucocorticoid receptor (GR) signaling influences lung maturation remain poorly understood. Deletion of the glucocorticoid receptor gene (Nr3c1) from fetal pulmonary mesenchymal cells phenocopied defects caused by global Nr3c1 deletion, while lung epithelial- or endothelial-specific Nr3c1 deletion did not impair lung function at birth. We integrated genome-wide gene expression profiling, ATAC-seq, and single cell RNA-seq data in mice in which GR was deleted or activated to identify the cellular and molecular mechanisms by which glucocorticoids control prenatal lung maturation. GR enhanced differentiation of a newly defined proliferative mesenchymal progenitor cell (PMP) into matrix fibroblasts (MFBs), in part by directly activating extracellular matrix-associated target genes, including Fn1, Col16a4, and Eln and by modulating VEGF, JAK-STAT, and WNT signaling. Loss of mesenchymal GR signaling blocked fibroblast progenitor differentiation into mature MFBs, which in turn increased proliferation of SOX9+ alveolar epithelial progenitor cells and inhibited differentiation of mature alveolar type II (AT2) and AT1 cells. GR signaling controls genes required for differentiation of a subset of proliferative mesenchymal progenitors into matrix fibroblasts, in turn, regulating signals controlling AT2/AT1 progenitor cell proliferation and differentiation and identifying cells and processes by which glucocorticoid signaling regulates fetal lung maturation.

Integration of transcriptomic and proteomic data identifies biological functions in cell populations from human infant lung.

Systems biology uses computational approaches to integrate diverse data types to understand cell and organ behavior. Data derived from complementary technologies, for example transcriptomic and proteomic analyses, are providing new insights into development and disease. We compared mRNA and protein profiles from purified endothelial, epithelial, immune, and mesenchymal cells from normal human infant lung tissue. Signatures for each cell type were identified and compared at both mRNA and protein levels. Cell-specific biological processes and pathways were predicted by analysis of concordant and discordant RNA-protein pairs. Cell clustering and gene set enrichment comparisons identified shared versus unique processes associated with transcriptomic and/or proteomic data. Clear cell-cell correlations between mRNA and protein data were obtained from each cell type. Approximately 40% of RNA-protein pairs were coherently expressed. While the correlation between RNA and their protein products was relatively low (Spearman rank coefficient rs ~0.4), cell-specific signature genes involved in functional processes characteristic of each cell type were more highly correlated with their protein products. Consistency of cell-specific RNA-protein signatures indicated an essential framework for the function of each cell type. Visualization and reutilization of the protein and RNA profiles are supported by a new web application, "LungProteomics," which is freely accessible to the public.

Neonatal Lung Disease Associated with TBX4 Mutations.

Variable lung disease was documented in 2 infants with heterozygous TBX4 mutations; their clinical presentations, pathology, and outcomes were distinct. These findings demonstrate that TBX4 gene mutations are associated with neonatal respiratory failure and highlight the wide spectrum of clinicopathological outcomes that have implications for patient diagnosis and management.

Clinical features related to statin-associated muscle symptoms.

INTRODUCTION: Statins reduce cardiovascular disease risk and are generally well tolerated, yet up to 0.5% of statin-treated patients develop incapacitating muscle symptoms including rhabdomyolysis. Our objective was to identify clinical factors related to statin-associated muscle symptoms (SAMS). METHODS: Clinical and laboratory characteristics were evaluated in 748 statin-treated Caucasians (634 with SAMS and 114 statin-tolerant controls). Information was collected on statin type, concomitant drug therapies, muscle symptom history, comorbidities, and family history. Logistic regression was used to identify associations. RESULTS: Individuals with SAMS were 3.6 times (odds ratio [OR] 3.60, 95% confidence interval [CI] 2.08-6.22) more likely than statin-tolerant controls to have a family history of heart disease. Additional associations included obesity (OR 3.08, 95% CI 1.18, 8.05), hypertension (OR 2.24, 95% CI 1.33, 3.77), smoking (OR 2.08, 95% CI 1.16, 3.74), and statin type. DISCUSSION: Careful medical monitoring of statin-treated patients with the associated coexisting conditions may ultimately reduce muscle symptoms and lead to improved compliance. Muscle Nerve 59:537-537, 2019.

Time-resolved proteome profiling of normal lung development.

Biochemical networks mediating normal lung morphogenesis and function have important implications for ameliorating morbidity and mortality in premature infants. Although several transcript-level studies have examined normal lung development, corresponding protein-level analyses are lacking. Here we performed proteomics analysis of murine lungs from embryonic to early adult ages to identify the molecular networks mediating normal lung development. We identified 8,932 proteins, providing a deep and comprehensive view of the lung proteome. Analysis of the proteomics data revealed discrete modules and the underlying regulatory and signaling network modulating their expression during development. Our data support the cell proliferation that characterizes early lung development and highlight responses of the lung to exposure to a nonsterile oxygen-rich ambient environment and the important role of lipid (surfactant) metabolism in lung development. Comparison of dynamic regulation of proteomic and recent transcriptomic analyses identified biological processes under posttranscriptional control. Our study provides a unique proteomic resource for understanding normal lung formation and function and can be freely accessed at Lungmap.net.

Airway Epithelial KIF3A Regulates Th2 Responses to Aeroallergens.

KIF3A, the gene encoding kinesin family member 3A, is a susceptibility gene locus associated with asthma; however, mechanisms by which KIF3A might influence the pathogenesis of the disorder are unknown. In this study, we deleted the mouse Kif3a gene in airway epithelial cells. Both homozygous and heterozygous Kif3a gene-deleted mice were highly susceptible to aeroallergens from Aspergillus fumigatus and the house dust mite, resulting in an asthma-like pathology characterized by increased goblet cell metaplasia, airway hyperresponsiveness, and Th2-mediated inflammation. Deletion of the Kif3a gene increased the severity of pulmonary eosinophilic inflammation and expression of cytokines (Il-4, Il-13, and Il-17a) and chemokine (Ccl11) RNAs following pulmonary exposure to Aspergillus extract. Inhibition of Kif3a disrupted the structure of motile cilia and impaired mucociliary clearance, barrier function, and epithelial repair, demonstrating additional mechanisms by which deficiency of KIF3A in respiratory epithelial cells contributes to pulmonary pathology. Airway epithelial KIF3A suppresses Th2 pulmonary inflammation and airway hyperresponsiveness following aeroallergen exposure, implicating epithelial microtubular functions in the pathogenesis of Th2-mediated lung pathology.

Lung Gene Expression Analysis (LGEA): an integrative web portal for comprehensive gene expression data analysis in lung development.

'LungGENS', our previously developed web tool for mapping single-cell gene expression in the developing lung, has been well received by the pulmonary research community. With continued support from the 'LungMAP' consortium, we extended the scope of the LungGENS database to accommodate transcriptomics data from pulmonary tissues and cells from human and mouse at different stages of lung development. Lung Gene Expression Analysis (LGEA) web portal is an extended version of LungGENS useful for the analysis, display and interpretation of gene expression patterns obtained from single cells, sorted cell populations and whole lung tissues. The LGEA web portal is freely available at http://research.cchmc.org/pbge/lunggens/mainportal.html.

CYP2D6 Genetic Variation and Beta-Blocker Maintenance Dose in Patients with Heart Failure.

PURPOSE: This study examined whether a CYP2D6 polymorphism (CYP2D6*4) was related to beta-blocker maintenance dose in patients with heart failure. METHODS: Logistic regression modeling was utilized in a retrospective chart-review analysis of heart-failure patients (60% Male, 90% of European descent) to assess whether CYP2D6*4 (non-functional CYP2D6 allele present in 1 of 5 individuals of European descent) is associated with maintenance dose of carvedilol (n = 65) or metoprolol (n = 33). RESULTS: CYP2D6*4 was associated with lower maintenance dose of metoprolol (OR 0.13 [95% CI 0.02-0.75] p = 0.023), and a trend was observed between CYP2D6*4 and higher maintenance dose of carvedilol (OR 2.94 [95% CI 0.84-10.30] p = 0.093). None of the patients that carried CYP2D6*4 achieved the recommended target dose of metoprolol (200 mg/day). CONCLUSION: Consistent with the role of CYP2D6 in the metabolism of metoprolol, the tolerated maintenance dose of metoprolol was lower in CYP2D6*4 carriers compared to non-carriers. Consistent with the role of CYP2D6 in activation of carvedilol, tolerated maintenance dose of carvedilol was higher in CYP2D6*4 carriers compared to non-carriers. Further investigation is warranted to ascertain the potential of CYP2D6 as a potential predictive biomarker of beta-blocker maintenance dose in heart failure patients.

An Open-Label Randomized Crossover Trial of Lyophilized Black Raspberries on Postprandial Inflammation in Older Overweight Males: A Pilot Study.

This study was a 14-day, outpatient, open-label randomized crossover trial of lyophilized black raspberries (BRBs) in older overweight or obese males to determine whether BRB consumption affects postprandial inflammation associated with consumption of a high-fat high-calorie (HFHC) meal. Ten study participants consumed 45 g/d of lyophilized BRBs for 4 days, followed by a HFHC breakfast plus BRBs on day 6 or consumed the HFHC breakfast on day 6 without previous consumption of BRBs and then crossed over to the other treatment after a 2-day washout period. Blood samples were obtained before and 1, 2, 4, 8, and 12 hours after consumption of the HFHC breakfast. The primary study outcomes were changes in area under the concentration-time curve (AUC) for interleukin-6 (IL-6), C-reactive protein (CRP), and tumor necrosis factor-alpha (TNF-α). The secondary outcomes were safety and tolerability of lyophilized BRB powder. The chronology and values of measured serum concentrations for IL-6, TNF-α, and CRP were consistent with those described previously by other investigators. The AUC of serum IL-6 was lowered significantly (P = 0.03, n = 10) with BRB consumption (34.3 ± 7.6 pg·mL⁻¹·h⁻¹ compared with 42.4 ± 17.9 pg·mL⁻¹·h⁻¹ for consumption of the HFHC meal alone). However, no significant differences (change in AUC) were calculated for serum CRP and TNF-α. The findings of this pilot study suggest that consumption of lyophilized BRBs may attenuate postprandial inflammation in overweight or obese males consuming a HFHC meal. Further investigation of BRBs is warranted to better elucidate their inflammomodulatory potential.

Individual and Combined Associations of Genetic Variants in CYP3A4, CYP3A5, and SLCO1B1 With Simvastatin and Simvastatin Acid Plasma Concentrations.

Our objective was to evaluate the associations of genetic variants affecting simvastatin (SV) and simvastatin acid (SVA) metabolism [the gene encoding cytochrome P450, family 3, subfamily A, polypeptide 4 (CYP3A4)*22 and the gene encoding cytochrome P450, family 3, subfamily A, polypeptide 5 (CYP3A5)*3] and transport [the gene encoding solute carrier organic anion transporter family member 1B1 (SLCO1B1) T521C] with 12-hour plasma SV and SVA concentrations. The variants were genotyped, and the concentrations were quantified by high performance liquid chromatography-tandem mass spectrometry in 646 participants of the Cholesterol and Pharmacogenetics clinical trial of 40 mg/d SV for 6 weeks. The genetic variants were tested for association with 12-hour plasma SV, SVA, or the SVA/SV ratio using general linear models. CYP3A5*3 was not significantly associated with 12-hour plasma SV or SVA concentration. CYP3A4*1/*22 participants had 58% higher 12-hour plasma SV concentration compared with CYP3A4*1/*1 participants (P = 0.006). SLCO1B1 521T/C and 521C/C participants had 71% (P < 0.001) and 248% (P < 0.001) higher 12-hour plasma SVA compared with SLCO1B1 521T/T participants, respectively. CYP3A4 and SLCO1B1 genotypes combined categorized participants into low (<1), intermediate (≈1), and high (>1) SVA/SV ratio groups (P = 0.001). In conclusion, CYP3A4*22 and SLCO1B1 521C were significantly associated with increased 12-hour plasma SV and SVA concentrations, respectively. CYP3A5*3 was not significantly associated with 12-hour plasma SV or SVA concentrations. The combination of CYP3A4*22 and SLCO1B1 521C was significantly associated with SVA/SV ratio, which may translate into different clinical SV risk/benefit profiles.

Foxa3 induces goblet cell metaplasia and inhibits innate antiviral immunity.

RATIONALE: Goblet cell metaplasia accompanies common pulmonary disorders that are prone to recurrent viral infections. Mechanisms regulating both goblet cell metaplasia and susceptibility to viral infection associated with chronic lung diseases are incompletely understood. OBJECTIVES: We sought to identify the role of the transcription factor FOXA3 in regulation of goblet cell metaplasia and pulmonary innate immunity. METHODS: FOXA3 was identified in airways from patients with asthma and chronic obstructive pulmonary disease. We produced transgenic mice conditionally expressing Foxa3 in airway epithelial cells and developed human bronchial epithelial cells expressing Foxa3. Foxa3-regulated genes were identified by immunostaining, Western blotting, and RNA analysis. Direct binding of FOXA3 to target genes was identified by chromatin immunoprecipitation sequencing correlated with RNA sequencing. MEASUREMENTS AND MAIN RESULTS: FOXA3 was highly expressed in airway goblet cells from patients with asthma and chronic obstructive pulmonary disease. FOXA3 was induced by either IL-13 or rhinovirus. Foxa3 induced goblet cell metaplasia and enhanced expression of a network of genes mediating mucus production. Paradoxically, FOXA3 inhibited rhinovirus-induced IFN production, IRF-3 phosphorylation, and IKKε expression and inhibited viral clearance and expression of genes required for antiviral defenses, including MDA5, RIG-I, TLR3, IRF7/9, and nuclear factor-κB. CONCLUSIONS: FOXA3 induces goblet cell metaplasia in response to infection or Th2 stimulation. Suppression of IFN signaling by FOXA3 provides a plausible mechanism that may serve to limit ongoing Th1 inflammation during the resolution of acute viral infection; however, inhibition of innate immunity by FOXA3 may contribute to susceptibility to viral infections associated with chronic lung disorders accompanied by chronic goblet cell metaplasia.

SAM-pointed domain ETS factor mediates epithelial cell-intrinsic innate immune signaling during airway mucous metaplasia.

Airway mucus plays a critical role in clearing inhaled toxins, particles, and pathogens. Diverse toxic, inflammatory, and infectious insults induce airway mucus secretion and goblet cell metaplasia to preserve airway sterility and homeostasis. However, goblet cell metaplasia, mucus hypersecretion, and airway obstruction are integral features of inflammatory lung diseases, including asthma, chronic obstructive lung disease, and cystic fibrosis, which cause an immense burden of morbidity and mortality. These chronic lung diseases are united by susceptibility to microbial colonization and recurrent airway infections. Whether these twinned phenomena (mucous metaplasia, compromised host defenses) are causally related has been unclear. Here, we demonstrate that SAM pointed domain ETS factor (SPDEF) was induced by rhinoviral infection of primary human airway cells and that cytoplasmic activities of SPDEF, a transcriptional regulator of airway goblet cell metaplasia, inhibited Toll-like receptor (TLR) activation of epithelial cells. SPDEF bound to and inhibited activities of TLR signaling adapters, MyD88 and TRIF, inhibiting MyD88-induced cytokine production and TRIF-induced interferon ß production. Conditional expression of SPDEF in airway epithelial cells in vivo inhibited LPS-induced neutrophilic infiltration and bacterial clearance. SPDEF-mediated inhibition of both TLR and type I interferon signaling likely protects the lung against inflammatory damage when inciting stimuli are not eradicated. Present findings provide, at least in part, a molecular explanation for increased susceptibility to infection in lung diseases associated with mucous metaplasia and a mechanism by which patients with florid mucous metaplasia may tolerate microbial burdens that are usually associated with fulminant inflammatory disease in normal hosts.