Influenza virus infection increases ACE2 expression and shedding in human small airway epithelial cells.

BACKGROUND: Patients with coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) demonstrate high rates of co-infection with respiratory viruses, including influenza A (IAV), suggesting pathogenic interactions. METHODS: We investigated how IAV may increase the risk of COVID-19 lung disease, focusing on the receptor angiotensin-converting enzyme (ACE)2 and the protease TMPRSS2, which cooperate in the intracellular uptake of SARS-CoV-2. RESULTS: We found, using single-cell RNA sequencing of distal human nondiseased lung homogenates, that at baseline, ACE2 is minimally expressed in basal, goblet, ciliated and secretory epithelial cells populating small airways. We focused on human small airway epithelial cells (SAECs), central to the pathogenesis of lung injury following viral infections. Primary SAECs from nondiseased donor lungs apically infected (at the air-liquid interface) with IAV (up to 3×105 pfu; ∼1 multiplicity of infection) markedly (eight-fold) boosted the expression of ACE2, paralleling that of STAT1, a transcription factor activated by viruses. IAV increased the apparent electrophoretic mobility of intracellular ACE2 and generated an ACE2 fragment (90 kDa) in apical secretions, suggesting cleavage of this receptor. In addition, IAV increased the expression of two proteases known to cleave ACE2, sheddase ADAM17 (TACE) and TMPRSS2 and increased the TMPRSS2 zymogen and its mature fragments, implicating proteolytic autoactivation. CONCLUSION: These results indicate that IAV amplifies the expression of molecules necessary for SARS-CoV-2 infection of the distal lung. Furthermore, post-translational changes in ACE2 by IAV may increase vulnerability to lung injury such as acute respiratory distress syndrome during viral co-infections. These findings support efforts in the prevention and treatment of influenza infections during the COVID-19 pandemic.

COVID-19-related Genes in Sputum Cells in Asthma. Relationship to Demographic Features and Corticosteroids.

Rationale: Coronavirus disease (COVID-19) is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). ACE2 (angiotensin-converting enzyme 2), and TMPRSS2 (transmembrane protease serine 2) mediate viral infection of host cells. We reasoned that differences in ACE2 or TMPRSS2 gene expression in sputum cells among patients with asthma may identify subgroups at risk for COVID-19 morbidity.Objectives: To determine the relationship between demographic features and sputum ACE2 and TMPRSS2 gene expression in asthma.Methods: We analyzed gene expression for ACE2 and TMPRSS2, and for ICAM-1 (intercellular adhesion molecule 1) (rhinovirus receptor as a comparator) in sputum cells from 330 participants in SARP-3 (Severe Asthma Research Program-3) and 79 healthy control subjects.Measurements and Main Results: Gene expression of ACE2 was lower than TMPRSS2, and expression levels of both genes were similar in asthma and health. Among patients with asthma, male sex, African American race, and history of diabetes mellitus were associated with higher expression of ACE2 and TMPRSS2. Use of inhaled corticosteroids (ICS) was associated with lower expression of ACE2 and TMPRSS2, but treatment with triamcinolone acetonide did not decrease expression of either gene. These findings differed from those for ICAM-1, where gene expression was increased in asthma and less consistent differences were observed related to sex, race, and use of ICS.Conclusions: Higher expression of ACE2 and TMPRSS2 in males, African Americans, and patients with diabetes mellitus provides rationale for monitoring these asthma subgroups for poor COVID-19 outcomes. The lower expression of ACE2 and TMPRSS2 with ICS use warrants prospective study of ICS use as a predictor of decreased susceptibility to SARS-CoV-2 infection and decreased COVID-19 morbidity.

Whole-Genome Sequencing Identifies Novel Functional Loci Associated with Lung Function in Puerto Rican Youth.

Rationale: Puerto Ricans have the highest childhood asthma prevalence in the United States (23.6%); however, the etiology is uncertain.Objectives: In this study, we sought to uncover the genetic architecture of lung function in Puerto Rican youth with and without asthma who were recruited from the island (n = 836).Methods: We used admixture-mapping and whole-genome sequencing data to discover genomic regions associated with lung function. Functional roles of the prioritized candidate SNPs were examined with chromatin immunoprecipitation sequencing, RNA sequencing, and expression quantitative trait loci data.Measurements and Main Results: We discovered a genomic region at 1q32 that was significantly associated with a 0.12-L decrease in the lung volume of exhaled air (95% confidence interval, -0.17 to -0.07; P = 6.62 × 10-8) with each allele of African ancestry. Within this region, two SNPs were expression quantitative trait loci of TMEM9 in nasal airway epithelial cells and MROH3P in esophagus mucosa. The minor alleles of these SNPs were associated with significantly decreased lung function and decreased TMEM9 gene expression. Another admixture-mapping peak was observed on chromosome 5q35.1, indicating that each Native American ancestry allele was associated with a 0.15-L increase in lung function (95% confidence interval, 0.08-0.21; P = 5.03 × 10-6). The region-based association tests identified four suggestive windows that harbored candidate rare variants associated with lung function.Conclusions: We identified common and rare genetic variants that may play a critical role in lung function among Puerto Rican youth. We independently validated an inflammatory pathway that could potentially be used to develop more targeted treatments and interventions for patients with asthma.

Genome-Wide Analysis Reveals Mucociliary Remodeling of the Nasal Airway Epithelium Induced by Urban PM2.5.

Air pollution particulate matter <2.5 µm (PM2.5) exposure is associated with poor respiratory outcomes. Mechanisms underlying PM2.5-induced lung pathobiology are poorly understood but likely involve cellular and molecular changes to the airway epithelium. We extracted and chemically characterized the organic and water-soluble components of air pollution PM2.5 samples, then determined the whole transcriptome response of human nasal mucociliary airway epithelial cultures to a dose series of PM2.5 extracts. We found that PM2.5 organic extract (OE), but not water-soluble extract, elicited a potent, dose-dependent transcriptomic response from the mucociliary epithelium. Exposure to a moderate OE dose modified the expression of 424 genes, including activation of aryl hydrocarbon receptor signaling and an IL-1 inflammatory program. We generated an OE-response gene network defined by eight functional enrichment groups, which exhibited high connectivity through CYP1A1, IL1A, and IL1B. This OE exposure also robustly activated a mucus secretory expression program (>100 genes), which included transcriptional drivers of mucus metaplasia (SPDEF and FOXA3). Exposure to a higher OE dose modified the expression of 1,240 genes and further exacerbated expression responses observed at the moderate dose, including the mucus secretory program. Moreover, the higher OE dose significantly increased the MUC5AC/MUC5B gel-forming mucin expression ratio and strongly downregulated ciliated cell expression programs, including key ciliating cell transcription factors (e.g., FOXJ1 and MCIDAS). Chronic OE stimulation induced mucus metaplasia-like remodeling characterized by increases in MUC5AC+ secretory cells and MUC5AC mucus secretions. This epithelial remodeling may underlie poor respiratory outcomes associated with high PM2.5 exposure.

Expression of SMARCD1 interacts with age in association with asthma control on inhaled corticosteroid therapy.

BACKGROUND: Global gene expression levels are known to be highly dependent upon gross demographic features including age, yet identification of age-related genomic indicators has yet to be comprehensively undertaken in a disease and treatment-specific context. METHODS: We used gene expression data from CD4+ lymphocytes in the Asthma BioRepository for Integrative Genomic Exploration (Asthma BRIDGE), an open-access collection of subjects participating in genetic studies of asthma with available gene expression data. Replication population participants were Puerto Rico islanders recruited as part of the ongoing Genes environments & Admixture in Latino Americans (GALA II), who provided nasal brushings for transcript sequencing. The main outcome measure was chronic asthma control as derived by questionnaires. Genomic associations were performed using regression of chronic asthma control score on gene expression with age in years as a covariate, including a multiplicative interaction term for gene expression times age. RESULTS: The SMARCD1 gene (SWI/SNF-related matrix-associated actin-dependent regulator of chromatin subfamily D member 1) interacted with age to influence chronic asthma control on inhaled corticosteroids, with a doubling of expression leading to an increase of 1.3 units of chronic asthma control per year (95% CI [0.86, 1.74], p = 6 × 10- 9), suggesting worsening asthma control with increasing age. This result replicated in GALA II (p = 3.8 × 10- 8). Cellular assays confirmed the role of SMARCD1 in glucocorticoid response in airway epithelial cells. CONCLUSION: Focusing on age-dependent factors may help identify novel indicators of asthma medication response. Age appears to modulate the effect of SMARCD1 on asthma control with inhaled corticosteroids.

Functional genomics of CDHR3 confirms its role in HRV-C infection and childhood asthma exacerbations.

BACKGROUND: Research in transformed immortalized cell lines indicates the cadherin-related family member 3 (CDHR3) protein serves as a receptor for human rhinovirus (HRV)-C. Similar experiments indicate that the CDHR3 coding variant rs6967330 increases CDHR3 protein surface expression. OBJECTIVE: We sought to determine whether CDHR3 is necessary for HRV-C infection of primary airway epithelial cells (AECs) and to identify molecular mechanisms by which CDHR3 variants confer risk for asthma exacerbations. METHODS: CDHR3 function and influence on HRV-C infection were investigated by using single-cell transcriptomics, CRISPR-Cas9 gene knockout, and genotype-specific donor experiments performed in primary AECs. Nasal airway epithelium cis-expression quantitative trait locus (eQTL) analysis of CDHR3 was performed, followed by association testing for asthma hospitalization in minority children. RESULTS: CDHR3 lung expression is exclusive to ciliated AECs and associated with basal bodies during and after motile ciliogenesis. Knockout of CDHR3 in human AECs did not prevent ciliated cell differentiation but was associated with a decrease in transepithelial resistance and an 80% decrease in HRV-C infection of the mucociliary epithelium. AECs from subjects homozygous for the risk-associated rs6967330 single nucleotide polymorphism (SNP) exhibited greater HRV-C infection compared with cells homozygous for the nonrisk allele. AEC cis-eQTL analysis indicated that rs6967330 and other SNPs are eQTLs for CDHR3. Only the eQTL block containing the rs6967330 SNP showed a significant association with childhood asthma hospitalization. CONCLUSIONS: Genetic deletion and genotype-specific studies in primary AECs indicate CDHR3 is critical to HRV-C infection of ciliated cells. The rs6967330 SNP confers risk of severe childhood asthma exacerbations, likely through increasing HRV-C infection levels and protein surface localization.

Analysis of Whole Exome Sequencing with Cardiometabolic Traits Using Family-Based Linkage and Association in the IRAS Family Study.

Family-based methods are a potentially powerful tool to identify trait-defining genetic variants in extended families, particularly when used to complement conventional association analysis. We utilized two-point linkage analysis and single variant association analysis to evaluate whole exome sequencing (WES) data from 1205 Hispanic Americans (78 families) from the Insulin Resistance Atherosclerosis Family Study. WES identified 211,612 variants above the minor allele frequency threshold of ≥0.005. These variants were tested for linkage and/or association with 50 cardiometabolic traits after quality control checks. Two-point linkage analysis yielded 10,580,600 logarithm of the odds (LOD) scores with 1148 LOD scores ≥3, 183 LOD scores ≥4, and 29 LOD scores ≥5. The maximal novel LOD score was 5.50 for rs2289043:T>C, in UNC5C with subcutaneous adipose tissue volume. Association analysis identified 13 variants attaining genome-wide significance (P < 5 × 10-08 ), with the strongest association between rs651821:C>T in APOA5 and triglyceride levels (P = 3.67 × 10-10 ). Overall, there was a 5.2-fold increase in the number of informative variants detected by WES compared to exome chip analysis in this population, nearly 30% of which were novel variants relative to the Database of Single Nucleotide Polymorphisms (dbSNP) build 138. Thus, integration of results from two-point linkage and single-variant association analysis from WES data enabled identification of novel signals potentially contributing to cardiometabolic traits.

Genome-wide linkage and association analysis of cardiometabolic phenotypes in Hispanic Americans.

Linkage studies of complex genetic diseases have been largely replaced by genome-wide association studies, due in part to limited success in complex trait discovery. However, recent interest in rare and low-frequency variants motivates re-examination of family-based methods. In this study, we investigated the performance of two-point linkage analysis for over 1.6 million single-nucleotide polymorphisms (SNPs) combined with single variant association analysis to identify high impact variants, which are both strongly linked and associated with cardiometabolic traits in up to 1414 Hispanics from the Insulin Resistance Atherosclerosis Family Study (IRASFS). Evaluation of all 50 phenotypes yielded 83 557 000 LOD (logarithm of the odds) scores, with 9214 LOD scores ⩾3.0, 845 ⩾4.0 and 89 ⩾5.0, with a maximal LOD score of 6.49 (rs12956744 in the LAMA1 gene for tumor necrosis factor-α (TNFα) receptor 2). Twenty-seven variants were associated with P<0.005 as well as having an LOD score >4, including variants in the NFIB gene under a linkage peak with TNFα receptor 2 levels on chromosome 9. Linkage regions of interest included a broad peak (31 Mb) on chromosome 1q with acute insulin response (max LOD=5.37). This region was previously documented with type 2 diabetes in family-based studies, providing support for the validity of these results. Overall, we have demonstrated the utility of two-point linkage and association in comprehensive genome-wide array-based SNP genotypes.

Mapping adipose and muscle tissue expression quantitative trait loci in African Americans to identify genes for type 2 diabetes and obesity.

Relative to European Americans, type 2 diabetes (T2D) is more prevalent in African Americans (AAs). Genetic variation may modulate transcript abundance in insulin-responsive tissues and contribute to risk; yet, published studies identifying expression quantitative trait loci (eQTLs) in African ancestry populations are restricted to blood cells. This study aims to develop a map of genetically regulated transcripts expressed in tissues important for glucose homeostasis in AAs, critical for identifying the genetic etiology of T2D and related traits. Quantitative measures of adipose and muscle gene expression, and genotypic data were integrated in 260 non-diabetic AAs to identify expression regulatory variants. Their roles in genetic susceptibility to T2D, and related metabolic phenotypes, were evaluated by mining GWAS datasets. eQTL analysis identified 1971 and 2078 cis-eGenes in adipose and muscle, respectively. Cis-eQTLs for 885 transcripts including top cis-eGenes CHURC1, USMG5, and ERAP2 were identified in both tissues. 62.1 % of top cis-eSNPs were within ±50 kb of transcription start sites and cis-eGenes were enriched for mitochondrial transcripts. Mining GWAS databases revealed association of cis-eSNPs for more than 50 genes with T2D (e.g. PIK3C2A, RBMS1, UFSP1), gluco-metabolic phenotypes (e.g. INPP5E, SNX17, ERAP2, FN3KRP), and obesity (e.g. POMC, CPEB4). Integration of GWAS meta-analysis data from AA cohorts revealed the most significant association for cis-eSNPs of ATP5SL and MCCC1 genes, with T2D and BMI, respectively. This study developed the first comprehensive map of adipose and muscle tissue eQTLs in AAs (publically accessible at https://mdsetaa.phs.wakehealth.edu ) and identified genetically regulated transcripts for delineating genetic causes of T2D, and related metabolic phenotypes.

Re-Sequencing of the APOL1-APOL4 and MYH9 Gene Regions in African Americans Does Not Identify Additional Risks for CKD Progression.

BACKGROUND: In African Americans (AAs), APOL1 G1 and G2 nephropathy risk variants are associated with non-diabetic end-stage kidney disease (ESKD) in an autosomal recessive pattern. Additional risk and protective genetic variants may be present near the APOL1 loci, since earlier age ESKD is observed in some AAs with one APOL1 renal-risk variant, and because the adjacent gene MYH9 is associated with nephropathy in populations lacking G1 and G2 variants. METHODS: Re-sequencing was performed across a ∼275 kb region encompassing the APOL1-APOL4 and MYH9 genes in 154 AA cases with non-diabetic ESKD and 38 controls without nephropathy who were heterozygous for a single APOL1 G1 or G2 risk variant. RESULTS: Sequencing identified 3,246 non-coding single nucleotide polymorphisms (SNPs), 55 coding SNPs, and 246 insertion/deletions. No new coding variations were identified. Eleven variants, including a rare APOL3 Gln58Ter null variant (rs11089781), were genotyped in a replication panel of 1,571 AA ESKD cases and 1,334 controls. After adjusting for APOL1 G1 and G2 risk effects, these variations were not significantly associated with ESKD. In subjects with <2 APOL1 G1 and/or G2 alleles (849 cases; 1,139 controls), the APOL3 null variant was nominally associated with ESKD (recessive model, OR 1.81; p = 0.026); however, analysis in 807 AA cases and 634 controls from the Family Investigation of Nephropathy and Diabetes did not replicate this association. CONCLUSION: Additional common variants in the APOL1-APOL4-MYH9 region do not contribute significantly to ESKD risk beyond the APOL1 G1 and G2 alleles.

The Puerto Rican Infant Metagenomic and Epidemiologic Study of Respiratory Outcomes (PRIMERO): Design and Baseline Characteristics for a Birth Cohort Study of Early-life Viral Respiratory Illnesses and Airway Dysfunction in Puerto Rican Children.

Epidemiologic studies demonstrate an association between early-life respiratory illnesses (RIs) and the development of childhood asthma. However, it remains uncertain whether these children are predisposed to both conditions or if early-life RIs induce alterations in airway function, immune responses, or other human biology that contribute to the development of asthma. Puerto Rican children experience a disproportionate burden of early-life RIs and asthma, making them an important population for investigating this complex interplay. PRIMERO, the Puerto Rican Infant Metagenomics and Epidemiologic Study of Respiratory Outcomes , recruited pregnant women and their newborns to investigate how the airways develop in early life among infants exposed to different viral RIs, and will thus provide a critical understanding of childhood asthma development. As the first asthma birth cohort in Puerto Rico, PRIMERO will prospectively follow 2,100 term healthy infants. Collected samples include post-term maternal peripheral blood, infant cord blood, the child's peripheral blood at the year two visit, and the child's nasal airway epithelium, collected using minimally invasive nasal swabs, at birth, during RIs over the first two years of life, and at annual healthy visits until age five. Herein, we describe the study's design, population, recruitment strategy, study visits and procedures, and primary outcomes.

A common polymorphism in the Intelectin-1 gene influences mucus plugging in severe asthma.

By incompletely understood mechanisms, type 2 (T2) inflammation present in the airways of severe asthmatics drives the formation of pathologic mucus which leads to airway mucus plugging. Here we investigate the molecular role and clinical significance of intelectin-1 (ITLN-1) in the development of pathologic airway mucus in asthma. Through analyses of human airway epithelial cells we find that ITLN1 gene expression is highly induced by interleukin-13 (IL-13) in a subset of metaplastic MUC5AC+ mucus secretory cells, and that ITLN-1 protein is a secreted component of IL-13-induced mucus. Additionally, we find ITLN-1 protein binds the C-terminus of the MUC5AC mucin and that its deletion in airway epithelial cells partially reverses IL-13-induced mucostasis. Through analysis of nasal airway epithelial brushings, we find that ITLN1 is highly expressed in T2-high asthmatics, when compared to T2-low children. Furthermore, we demonstrate that both ITLN-1 gene expression and protein levels are significantly reduced by a common genetic variant that is associated with protection from the formation of mucus plugs in T2-high asthma. This work identifies an important biomarker and targetable pathways for the treatment of mucus obstruction in asthma.

Preferential transmission of genetic risk variants of candidate loci at 6p21 from asymptomatic grandparents to mothers of children with neonatal lupus.

OBJECTIVE: Neonatal lupus (NL) occurs in fetuses exposed to maternal anti-SSA/Ro and/or anti-SSB/La antibodies, although the mothers themselves may not manifest any clinical disease. A focus on transmission of risk factors for NL from maternal grandparents to mothers of children with NL may yield dividends toward understanding the aggregation of autoantibodies and genetic factors in affected families. This study was perforned to determine the role of maternal grandparents in the development of the autoimmune phenotype of mothers of children with NL. METHODS: Fifty-one mothers of children with cardiac and/or cutaneous NL, 48 maternal grandmothers, and 35 maternal grandfathers in the Research Registry for Neonatal Lupus were interrogated for clinical symptoms by questionnaire and underwent laboratory assessments, including determination of anti-SSA/Ro and anti-SSB/La antibody status (by enzyme-linked immunosorbent assay) and genotype at rs1800629 (TNFα) and rs7775397 (C6orf10) (allelic discrimination). The transmission disequilibrium test (TDT) was computed to test for nonrandom transmission from maternal grandparents to mothers of children with NL. RESULTS: The common phenotypic feature in mothers of children with NL was the autoantibody and not the clinical profile; 7 had lupus, 14 had Sjögren's syndrome, 7 had both, and 23 were asymptomatic. Mothers of children with NL were significantly enriched for the risk alleles at both TNFα and C6orf10. The grandparents of children with NL carried minimal burden for autoimmune disease or abnormal antibody production and were not enriched in the genetic risk factors. However, the TDT analysis showed significant excess transmission of the risk alleles at both TNFα (odds ratio [OR] 6.67, P = 3.93 × 10(-4) ) and C6orf10 (OR 35.0, P = 3.74 × 10(-5) ) to mothers of children with NL. CONCLUSION: Mothers of children with NL are enriched for the TNFα and C6orf10 risk alleles, which are preferentially inherited from the asymptomatic maternal grandparents. These findings support the hypothesis that the development of NL and genetic etiology are multigenerational.

Anionic Pulmonary Surfactant Lipid Treatment Inhibits Rhinovirus A Infection of the Human Airway Epithelium.

Rhinoviruses (RVs) are major instigators of acute exacerbations of asthma, COPD, and other respiratory diseases. RVs are categorized into three species (RV-A, RV-B, and RV-C), which comprise more than 160 serotypes, making it difficult to develop an effective vaccine. Currently, no effective treatment for RV infection is available. Pulmonary surfactant is an extracellular complex of lipids and proteins that plays a central role in regulating innate immunity in the lung. The minor pulmonary surfactant lipids, palmitoyl-oleoyl-phosphatidylglycerol (POPG) and phosphatidylinositol (PI), are potent regulators of inflammatory processes and exert antiviral activity against respiratory syncytial virus (RSV) and influenza A viruses (IAV). In the current study, we examined the potencies of POPG and PI against rhinovirus A16 (RV-A16) in primary human airway epithelial cells (AECs) differentiated at an air-liquid interface (ALI). After AECs were infected with RV-A16, PI reduced the viral RNA copy number by 70% and downregulated (55-75%) the expression of antiviral (MDA5, IRF7, and IFN-lambda) and CXCL11 chemokine genes. In contrast, POPG only slightly decreased MDA5 (24%) and IRF7 (11%) gene expression but did not inhibit IFN-lambda gene expression or RV-A16 replication in AECs. However, both POPG and PI inhibited (50-80%) IL6 gene expression and protein secretion and CXCL11 protein secretion. PI treatment dramatically attenuated global gene expression changes induced by RV-A16 infection alone in AECs. The observed inhibitory effects were indirect and resulted mainly from the inhibition of virus replication. Cell-type enrichment analysis of viral-regulated genes opposed by PI treatment revealed the PI-inhibited viral induction of goblet cell metaplasia and the virus-induced downregulation of ciliated, club, and ionocyte cell types. Notably, the PI treatment also altered the ability of RV-A16 to regulate the expression of some phosphatidylinositol 4-kinase (PI4K); acyl-CoA-binding, domain-containing (ACBD); and low-density lipoprotein receptor (LDLR) genes that play critical roles in the formation and functioning of replication organelles (ROs) required for RV replication in host cells. These data suggest PI can be used as a potent, non-toxic, antiviral agent for RV infection prophylaxis and treatment.

Nasal airway transcriptome-wide association study of asthma reveals genetically driven mucus pathobiology.

To identify genetic determinants of airway dysfunction, we performed a transcriptome-wide association study for asthma by combining RNA-seq data from the nasal airway epithelium of 681 children, with UK Biobank genetic association data. Our airway analysis identified 95 asthma genes, 58 of which were not identified by transcriptome-wide association analyses using other asthma-relevant tissues. Among these genes were MUC5AC, an airway mucin, and FOXA3, a transcriptional driver of mucus metaplasia. Muco-ciliary epithelial cultures from genotyped donors revealed that the MUC5AC risk variant increases MUC5AC protein secretion and mucus secretory cell frequency. Airway transcriptome-wide association analyses for mucus production and chronic cough also identified MUC5AC. These cis-expression variants were associated with trans effects on expression; the MUC5AC variant was associated with upregulation of non-inflammatory mucus secretory network genes, while the FOXA3 variant was associated with upregulation of type-2 inflammation-induced mucus-metaplasia pathway genes. Our results reveal genetic mechanisms of airway mucus pathobiology.

SARS-CoV-2 infection produces chronic pulmonary epithelial and immune cell dysfunction with fibrosis in mice.

A subset of individuals who recover from coronavirus disease 2019 (COVID-19) develop post-acute sequelae of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) (PASC), but the mechanistic basis of PASC-associated lung abnormalities suffers from a lack of longitudinal tissue samples. The mouse-adapted SARS-CoV-2 strain MA10 produces an acute respiratory distress syndrome in mice similar to humans. To investigate PASC pathogenesis, studies of MA10-infected mice were extended from acute to clinical recovery phases. At 15 to 120 days after virus clearance, pulmonary histologic findings included subpleural lesions composed of collagen, proliferative fibroblasts, and chronic inflammation, including tertiary lymphoid structures. Longitudinal spatial transcriptional profiling identified global reparative and fibrotic pathways dysregulated in diseased regions, similar to human COVID-19. Populations of alveolar intermediate cells, coupled with focal up-regulation of profibrotic markers, were identified in persistently diseased regions. Early intervention with antiviral EIDD-2801 reduced chronic disease, and early antifibrotic agent (nintedanib) intervention modified early disease severity. This murine model provides opportunities to identify pathways associated with persistent SARS-CoV-2 pulmonary disease and test countermeasures to ameliorate PASC.

A model of persistent post SARS-CoV-2 induced lung disease for target identification and testing of therapeutic strategies.

COVID-19 survivors develop post-acute sequelae of SARS-CoV-2 (PASC), but the mechanistic basis of PASC-associated lung abnormalities suffers from a lack of longitudinal samples. Mouse-adapted SARS-CoV-2 MA10 produces an acute respiratory distress syndrome (ARDS) in mice similar to humans. To investigate PASC pathogenesis, studies of MA10-infected mice were extended from acute disease through clinical recovery. At 15-120 days post-virus clearance, histologic evaluation identified subpleural lesions containing collagen, proliferative fibroblasts, and chronic inflammation with tertiary lymphoid structures. Longitudinal spatial transcriptional profiling identified global reparative and fibrotic pathways dysregulated in diseased regions, similar to human COVID-19. Populations of alveolar intermediate cells, coupled with focal upregulation of pro-fibrotic markers, were identified in persistently diseased regions. Early intervention with antiviral EIDD-2801 reduced chronic disease, and early anti-fibrotic agent (nintedanib) intervention modified early disease severity. This murine model provides opportunities to identify pathways associated with persistent SARS-CoV-2 pulmonary disease and test countermeasures to ameliorate PASC.

Pharmacogenetic studies of long-acting beta agonist and inhaled corticosteroid responsiveness in randomised controlled trials of individuals of African descent with asthma.

BACKGROUND: Pharmacogenetic studies in asthma cohorts, primarily made up of White people of European descent, have identified loci associated with response to inhaled beta agonists and corticosteroids (ICSs). Differences exist in how individuals from different ancestral backgrounds respond to long-acting beta agonist (LABA) and ICSs. Therefore, we sought to understand the pharmacogenetic mechanisms regulating therapeutic responsiveness in individuals of African descent. METHODS: We did ancestry-based pharmacogenetic studies of children (aged 5-11 years) and adolescents and adults (aged 12-69 years) from the Best African Response to Drug (BARD) trials, in which participants with asthma uncontrolled with low-dose ICS (fluticasone propionate 50 µg in children, 100 µg in adolescents and adults) received different step-up combination therapies. The hierarchal composite outcome of pairwise superior responsiveness in BARD was based on asthma exacerbations, a 31-day difference in annualised asthma-control days, or a 5% difference in percentage predicted FEV1. We did whole-genome admixture mapping of 15 159 ancestral segments within 312 independent regions, stratified by the two age groups. The two co-primary outcome comparisons were the step up from low-dose ICS to the quintuple dose of ICS (5 × ICS: 250 µg twice daily in children and 500 µg twice daily in adolescents and adults) versus double dose (2-2·5 × ICS: 100 µg twice daily in children, 250 µg twice daily in adolescents and adults), and 5 × ICS versus 100 µg fluticasone plus a LABA (salmeterol 50 µg twice daily). We used a genome-wide significance threshold of p<1·6 × 10-4, and tested for replication using independent cohorts of individuals of African descent with asthma. FINDINGS: We included 249 unrelated children and 267 unrelated adolescents and adults in the BARD pharmacogenetic analysis. In children, we identified a significant admixture mapping peak for superior responsiveness to 5 × ICS versus 100 µg fluticasone plus salmeterol on chromosome 12 (odds ratio [ORlocal African] 3·95, 95% CI 2·02-7·72, p=6·1 × 10-5) fine mapped to a locus adjacent to RNFT2 and NOS1 (rs73399224, ORallele dose 0·17, 95% CI 0·07-0·42, p=8·4 × 10-5). In adolescents and adults, we identified a peak for superior responsiveness to 5 × ICS versus 2·5 × ICS on chromosome 22 (ORlocal African 3·35, 1·98-5·67, p=6·8 × 10-6) containing a locus adjacent to TPST2 (rs5752429, ORallele dose 0·21, 0·09-0·52, p=5·7 × 10-4). We replicated rs5752429 and nominally replicated rs73399224 in independent African American cohorts. INTERPRETATION: BARD is the first genome-wide pharmacogenetic study of LABA and ICS response in clinical trials of individuals of African descent to detect and replicate genome-wide significant loci. Admixture mapping of the composite BARD trial outcome enabled the identification of novel pharmacogenetic variation accounting for differential therapeutic responses in people of African descent with asthma. FUNDING: National Institutes of Health, National Heart, Lung, and Blood Institute.

Loss of Fas signaling in fibroblasts impairs homeostatic fibrosis resolution and promotes persistent pulmonary fibrosis.

Idiopathic pulmonary fibrosis (IPF) is a progressive, irreversible fibrotic disease of the distal lung alveoli that culminates in respiratory failure and reduced lifespan. Unlike normal lung repair in response to injury, IPF is associated with the accumulation and persistence of fibroblasts and myofibroblasts, as well as continued production of collagen and other extracellular matrix (ECM) components. Prior in vitro studies have led to the hypothesis that the development of resistance to Fas-induced apoptosis by lung fibroblasts and myofibroblasts contributes to their accumulation in the distal lung tissues of IPF patients. Here, we test this hypothesis in vivo in the resolving model of bleomycin-induced pulmonary fibrosis in mice. Using genetic loss-of-function approaches to inhibit Fas signaling in fibroblasts, potentially novel flow cytometry strategies to quantify lung fibroblast subsets, and transcriptional profiling of lung fibroblasts by bulk and single cell RNA sequencing, we show that Fas is necessary for lung fibroblast apoptosis during homeostatic resolution of bleomycin-induced pulmonary fibrosis in vivo. Furthermore, we show that loss of Fas signaling leads to the persistence and continued profibrotic functions of lung fibroblasts. Our studies provide insights into the mechanisms that contribute to fibroblast survival, persistence, and continued ECM deposition in the context of IPF and how failure to undergo Fas-induced apoptosis impairs fibrosis resolution.

Dissecting the cellular specificity of smoking effects and reconstructing lineages in the human airway epithelium.

Cigarette smoke first interacts with the lung through the cellularly diverse airway epithelium and goes on to drive development of most chronic lung diseases. Here, through single cell RNA-sequencing analysis of the tracheal epithelium from smokers and non-smokers, we generate a comprehensive atlas of epithelial cell types and states, connect these into lineages, and define cell-specific responses to smoking. Our analysis infers multi-state lineages that develop into surface mucus secretory and ciliated cells and then contrasts these to the unique specification of submucosal gland (SMG) cells. Accompanying knockout studies reveal that tuft-like cells are the likely progenitor of both pulmonary neuroendocrine cells and CFTR-rich ionocytes. Our smoking analysis finds that all cell types, including protected stem and SMG populations, are affected by smoking through both pan-epithelial smoking response networks and hundreds of cell-specific response genes, redefining the penetrance and cellular specificity of smoking effects on the human airway epithelium.