Colocalization analysis of 3' UTR alternative polyadenylation quantitative trait loci reveals novel mechanisms underlying associations with lung function.

While many disease-associated single nucleotide polymorphisms (SNPs) are expression quantitative trait loci (eQTLs), a large proportion of genome-wide association study (GWAS) variants are of unknown function. Alternative polyadenylation (APA) plays an important role in posttranscriptional regulation by allowing genes to shorten or extend 3' untranslated regions (UTRs). We hypothesized that genetic variants that affect APA in lung tissue may lend insight into the function of respiratory associated GWAS loci. We generated alternative polyadenylation (apa) QTLs using RNA sequencing and whole genome sequencing on 1241 subjects from the Lung Tissue Research Consortium (LTRC) as part of the NHLBI TOPMed project. We identified 56 179 APA sites corresponding to 13 582 unique genes after filtering out APA sites with low usage. We found that a total of 8831 APA sites were associated with at least one SNP with q-value < 0.05. The genomic distribution of lead APA SNPs indicated that the majority are intronic variants (33%), followed by downstream gene variants (26%), 3' UTR variants (17%), and upstream gene variants (within 1 kb region upstream of transcriptional start site, 10%). APA sites in 193 genes colocalized with GWAS data for at least one phenotype. Genes containing the top APA sites associated with GWAS variants include membrane associated ring-CH-type finger 2 (MARCHF2), nectin cell adhesion molecule 2 (NECTIN2), and butyrophilin subfamily 3 member A2 (BTN3A2). Overall, these findings suggest that APA may be an important mechanism for genetic variants in lung function and chronic obstructive pulmonary disease (COPD).

Lung Transcriptomics Links Emphysema to Barrier Dysfunction and Macrophage Subpopulations.

RATIONALE: While many studies have examined gene expression in lung tissue, the gene regulatory processes underlying emphysema are still not well understood. Finding efficient non-imaging screening methods and disease-modifying therapies has been challenging, but knowledge of the transcriptomic features of emphysema may help in this effort. OBJECTIVES: Our goals were to identify emphysema-associated biological pathways through transcriptomic analysis of bulk lung tissue, to determine the lung cell types in which these emphysema-associated pathways are altered, and to detect unique and overlapping transcriptomic signatures in blood and lung samples. METHODS: Using RNA-sequencing data from 446 samples in the Lung Tissue Research Consortium (LTRC) and 3,606 blood samples from the COPDGene study, we examined the transcriptomic features of chest computed tomography-quantified emphysema. We also leveraged publicly available lung single-cell RNA-sequencing data to identify cell types showing COPD-associated differential expression of the emphysema pathways found in the bulk analyses. MEASUREMENTS AND MAIN RESULTS: In the bulk lung RNA-seq analysis, 1,087 differentially expressed genes and 34 dysregulated pathways were significantly associated with emphysema. We observed alternative splicing of several genes and increased activity in pluripotency and cell barrier function pathways. Lung tissue and blood samples shared differentially expressed genes and biological pathways. Multiple lung cell types displayed dysregulation of epithelial barrier function pathways, and distinct pathway activities were observed among various macrophage subpopulations. CONCLUSIONS: This study identified emphysema-related changes in gene expression and alternative splicing, cell-type specific dysregulated pathways, and instances of shared pathway dysregulation between blood and lung.

A novel piwi-interacting RNA associates with type 2-high asthma phenotypes.

BACKGROUND: Piwi-interacting RNAs (piRNAs), comprising the largest noncoding RNA group, regulate transcriptional processes. Whether piRNAs are associated with type 2 (T2)-high asthma is unknown. OBJECTIVE: We sought to investigate the association between piRNAs and T2-high asthma in childhood asthma. METHODS: We sequenced plasma samples from 462 subjects in the Childhood Asthma Management Program (CAMP) as the discovery cohort and 1165 subjects in the Genetics of Asthma in Costa Rica Study (GACRS) as a replication cohort. Sequencing reads were filtered first, and piRNA reads were annotated and normalized. Linear regression was used for the association analysis of piRNAs and peripheral blood eosinophil count, total serum IgE level, and long-term asthma exacerbation in children with asthma. Mediation analysis was performed to investigate the effect direction. We then ascertained if the circulating piRNAs were present in asthmatic airway epithelial cells in a Gene Expression Omnibus (GEO; www.ncbi.nlm.nih.gov/geo) public data set. RESULTS: Fifteen piRNAs were significantly associated with eosinophil count in CAMP (P ≤ .05), and 3 were successfully replicated in GACRS. Eleven piRNAs were associated with total IgE in CAMP, and one of these was replicated in GACRS. All 22 significant piRNAs were identified in epithelial cells in vitro, and 6 of these were differentially expressed between subjects with asthma and healthy controls. Fourteen piRNAs were associated with long-term asthma exacerbation, and effect of piRNAs on long-term asthma exacerbation are mediated through eosinophil count and serum IgE level. CONCLUSION: piRNAs are associated with peripheral blood eosinophils and total serum IgE in childhood asthma and may play important roles in T2-high asthma.

Global Burden of Digestive Diseases: A Systematic Analysis of the Global Burden of Diseases Study, 1990 to 2019.

BACKGROUND & AIMS: This study assessed the worldwide burden of digestive diseases between 1990 and 2019. METHODS: We analyzed data from the Global Burden of Diseases study, covering 18 digestive diseases across 204 countries and territories. Key disease burden indicators, including incidence, prevalence, mortality, and disability-adjusted life years (DALYs), were studied. Linear regression analysis was applied to the natural logarithm of age-standardized outcomes to determine the annual percent change. RESULTS: In 2019, there were 7.32 billion incidents and 2.86 billion prevalent cases of digestive diseases, resulting in 8 million deaths and 277 million DALYs lost. Little to no decrease in global age-standardized incidence and prevalence of digestive diseases was observed between 1990 and 2019, with 95,582 and 35,106 cases per 100,000 individuals in 2019, respectively. The age-standardized death rate was 102 per 100,000 individuals. Digestive diseases accounted for a significant portion of the overall disease burden, with more than one-third of prevalent cases having a digestive etiology. Enteric infections were the primary contributor to incidence, death, and DALYs lost, whereas cirrhosis and other chronic liver diseases had the highest prevalence rate. The burden of digestive diseases was inversely related to the sociodemographic index, with enteric infections being the predominant cause of death in low and low-middle quintiles and colorectal cancer in the high quintile. CONCLUSIONS: Despite significant reductions in deaths and DALYs due to digestive diseases from 1990 to 2019, they remain prevalent. A significant disparity in the burden of digestive diseases exists among countries with different development levels.

Blood Gene Expression and Immune Cell Subtypes Associated with Chronic Obstructive Pulmonary Disease Exacerbations.

Rationale: Acute exacerbations of chronic obstructive pulmonary disease (AE-COPDs) are associated with a significant disease burden. Blood immune phenotyping may improve our understanding of a COPD endotype at increased risk of exacerbations. Objective: To determine the relationship between the transcriptome of circulating leukocytes and COPD exacerbations. Methods: Blood RNA sequencing data (n = 3,618) from the COPDGene (Genetic Epidemiology of COPD) study were analyzed. Blood microarray data (n = 646) from the ECLIPSE (Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints) study were used for validation. We tested the association between blood gene expression and AE-COPDs. We imputed the abundance of leukocyte subtypes and tested their association with prospective AE-COPDs. Flow cytometry was performed on blood in SPIROMICS (Subpopulations and Intermediate Outcomes in COPD Study) (n = 127), and activation markers for T cells were tested for association with prospective AE-COPDs. Measurements and Main Results: Exacerbations were reported 4,030 and 2,368 times during follow-up in COPDGene (5.3 ± 1.7 yr) and ECLIPSE (3 yr), respectively. We identified 890, 675, and 3,217 genes associated with a history of AE-COPDs, persistent exacerbations (at least one exacerbation per year), and prospective exacerbation rate, respectively. In COPDGene, the number of prospective exacerbations in patients with COPD (Global Initiative for Chronic Obstructive Lung Disease stage ⩾2) was negatively associated with circulating CD8+ T cells, CD4+ T cells, and resting natural killer cells. The negative association with naive CD4+ T cells was replicated in ECLIPSE. In the flow-cytometry study, an increase in CTLA4 on CD4+ T cells was positively associated with AE-COPDs. Conclusions: Individuals with COPD with lower circulating lymphocyte counts, particularly decreased CD4+ T cells, are more susceptible to AE-COPDs, including persistent exacerbations.

Blood RNA sequencing shows overlapping gene expression across COPD phenotype domains.

RATIONALE: COPD can be assessed using multidimensional grading systems with components from three domains: pulmonary function tests, symptoms and systemic features. Clinically, measures may be used interchangeably, though it is not known if they share similar pathobiology. OBJECTIVE: To use RNA sequencing (RNA-seq) to determine if there is an overlap in the underlying biological mechanisms and consequences driving different components of the multidimensional grading systems. METHODS: Whole blood was collected for RNA-seq from current and former smokers in the Genetic Epidemiology of COPD study. We tested the overlap in gene expression and biological pathways associated with case-control status and quantitative COPD phenotypes within and between the three domains. RESULTS: In 2647 subjects, there were 3030 genes differentially expressed in any of the three domains or case-control status. There were five genes that overlapped between the three domains and case-control status, including G protein-coupled receptor 15(GPR15), sestrin 1 (SESN1) and interferon-induced guanylate-binding protein 1 (GBP1), which were associated with longitudinal decline in FEV1. The overlap between the three domains was enriched for pathways related to cellular components. CONCLUSIONS: We identified gene sets and pathways that overlap between 12 COPD-related phenotypes and case-control status. There were no pathways represented in the overlap between the three domains and case-control status, but we identified multiple genes that demonstrated a consistent pattern of expression across several of the phenotypes. Patterns of gene expression correlation were generally similar to the correlation of clinical phenotypes in the PFT and symptom domains but not the systemic features.

Lung function, airway and peripheral basophils and eosinophils are associated with molecular pharmacogenomic endotypes of steroid response in severe asthma.

INTRODUCTION: Asthma is a complex disease with heterogeneous expression/severity. There is growing interest in defining asthma endotypes consistently associated with different responses to therapy, focusing on type 2 inflammation (Th2) as a key pathological mechanism. Current asthma endotypes are defined primarily by clinical/laboratory criteria. Each endotype is likely characterised by distinct molecular mechanisms that identify optimal therapies. METHODS: We applied unsupervised (without a priori clinical criteria) principal component analysis on sputum airway cells RNA-sequencing transcriptomic data from 19 asthmatics from the Severe Asthma Research Program at baseline and 6-8 weeks follow-up after a 40 mg dose of intramuscular corticosteroids. We investigated principal components PC1, PC3 for association with 55 clinical variables. RESULTS: PC3 was associated with baseline Th2 clinical features including blood (rank-sum p=0.0082) and airway (rank-sum p=0.0024) eosinophilia, FEV1 change (Kendall tau-b R=-0.333 (-0.592 to -0.012)) and follow-up FEV1 albuterol response (Kendall tau-b R=0.392 (0.079 to 0.634)). PC1 with blood basophlia (rank-sum p=0.0191). The top 5% genes contributing to PC1, PC3 were enriched for distinct immune system/inflammation ontologies suggesting distinct subject-specific clusters of transcriptomic response to corticosteroids. PC3 association with FEV1 change was reproduced in silico in a comparable independent 14-subject (baseline, 8 weeks after daily inhaled corticosteroids (ICS)) airway epithelial cells microRNAome dataset. CONCLUSIONS: Transcriptomic PCs from this unsupervised methodology define molecular pharmacogenomic endotypes that may yield novel biology underlying different subject-specific responses to corticosteroid therapy in asthma, and optimal personalised asthma care. Top contributing genes to these PCs may suggest new therapeutic targets.

An interferon-inducible signature of airway disease from blood gene expression profiling.

BACKGROUND: The molecular basis of airway remodelling in chronic obstructive pulmonary disease (COPD) remains poorly understood. We identified gene expression signatures associated with chest computed tomography (CT) scan airway measures to understand molecular pathways associated with airway disease. METHODS: In 2396 subjects in the COPDGene Study, we examined the relationship between quantitative CT airway phenotypes and blood transcriptomes to identify airway disease-specific genes and to define an airway wall thickness (AWT) gene set score. Multivariable regression analyses were performed to identify associations of the AWT score with clinical phenotypes, bronchial gene expression and genetic variants. RESULTS: Type 1 interferon (IFN)-induced genes were consistently associated with AWT, square root wall area of a hypothetical airway with 10 mm internal perimeter (Pi10) and wall area percentage, with the strongest enrichment in AWT. A score derived from 18 genes whose expression was associated with AWT was associated with COPD-related phenotypes including reduced lung function (forced expiratory volume in 1 s percentage predicted ß= -3.4; p<0.05) and increased exacerbations (incidence rate ratio 1.7; p<0.05). The AWT score was reproducibly associated with AWT in bronchial samples from 23 subjects (ß=3.22; p<0.05). The blood AWT score was associated with genetic variant rs876039, an expression quantitative trait locus for IKZF1, a gene that regulates IFN signalling and is associated with inflammatory diseases. CONCLUSIONS: A gene expression signature with IFN-stimulated genes from peripheral blood and bronchial brushings is associated with CT AWT, lung function and exacerbations. Shared genes and genetic associations suggest viral responses and/or autoimmune dysregulation as potential underlying mechanisms of airway disease in COPD.

Lung tissue shows divergent gene expression between chronic obstructive pulmonary disease and idiopathic pulmonary fibrosis.

BACKGROUND: Chronic obstructive pulmonary disease (COPD) and idiopathic pulmonary fibrosis (IPF) are characterized by shared exposures and clinical features, but distinct genetic and pathologic features exist. These features have not been well-studied using large-scale gene expression datasets. We hypothesized that there are divergent gene, pathway, and cellular signatures between COPD and IPF. METHODS: We performed RNA-sequencing on lung tissues from individuals with IPF (n = 231) and COPD (n = 377) compared to control (n = 267), defined as individuals with normal spirometry. We grouped the overlapping differential expression gene sets based on direction of expression and examined the resultant sets for genes of interest, pathway enrichment, and cell composition. Using gene set variation analysis, we validated the overlap group gene sets in independent COPD and IPF data sets. RESULTS: We found 5010 genes differentially expressed between COPD and control, and 11,454 genes differentially expressed between IPF and control (1% false discovery rate). 3846 genes overlapped between IPF and COPD. Several pathways were enriched for genes upregulated in COPD and downregulated in IPF; however, no pathways were enriched for genes downregulated in COPD and upregulated in IPF. There were many myeloid cell genes with increased expression in COPD but decreased in IPF. We found that the genes upregulated in COPD but downregulated in IPF were associated with lower lung function in the independent validation cohorts. CONCLUSIONS: We identified a divergent gene expression signature between COPD and IPF, with increased expression in COPD and decreased in IPF. This signature is associated with worse lung function in both COPD and IPF.

Analysis of genetically driven alternative splicing identifies FBXO38 as a novel COPD susceptibility gene.

While many disease-associated single nucleotide polymorphisms (SNPs) are associated with gene expression (expression quantitative trait loci, eQTLs), a large proportion of complex disease genome-wide association study (GWAS) variants are of unknown function. Some of these SNPs may contribute to disease by regulating gene splicing. Here, we investigate whether SNPs that are associated with alternative splicing (splice QTL or sQTL) can identify novel functions for existing GWAS variants or suggest new associated variants in chronic obstructive pulmonary disease (COPD). RNA sequencing was performed on whole blood from 376 subjects from the COPDGene Study. Using linear models, we identified 561,060 unique sQTL SNPs associated with 30,333 splice sites corresponding to 6,419 unique genes. Similarly, 708,928 unique eQTL SNPs involving 15,913 genes were detected at 10% FDR. While there is overlap between sQTLs and eQTLs, 55.3% of sQTLs are not eQTLs. Co-localization analysis revealed that 7 out of 21 loci associated with COPD (p<1x10-6) in a published GWAS have at least one shared causal variant between the GWAS and sQTL studies. Among the genes identified to have splice sites associated with top GWAS SNPs was FBXO38, in which a novel exon was discovered to be protective against COPD. Importantly, the sQTL in this locus was validated by qPCR in both blood and lung tissue, demonstrating that splice variants relevant to lung tissue can be identified in blood. Other identified genes included CDK11A and SULT1A2. Overall, these data indicate that analysis of alternative splicing can provide novel insights into disease mechanisms. In particular, we demonstrated that SNPs in a known COPD GWAS locus on chromosome 5q32 influence alternative splicing in the gene FBXO38.

Circulating MicroRNAs and Treatment Response in Childhood Asthma.

Rationale: Inhaled corticosteroids (ICS) are key treatments for controlling asthma and preventing asthma attacks. However, the responsiveness to ICS varies among individuals. MicroRNAs (miRNAs) have been lauded for their prognostic utility.Objectives: We hypothesized that circulating miRNAs obtained at baseline/prerandomization in the Childhood Asthma Management Program (CAMP) could serve as biomarkers and biologic mediators of ICS clinical response over the 4-year clinical trial period.Methods: We selected baseline serum samples from 462 CAMP subjects subsequently randomized to either ICS (budesonide) or placebo. Samples underwent small RNA sequencing, and read counts were normalized and filtered by depth and coverage. Linear regression was used to associate miRNAs with change in FEV1% (prebronchodilator FEV1 as a percent predicted) over the 4-year treatment period in both main effects and interaction models. We validated the function of the top associated miRNAs by luciferase reporter assays of glucocorticoid-mediated transrepression and predicted response to ICS through logistic regression models.Measurements and Main Results: We identified 7 miRNAs significantly associated with FEV1% change (P ≤ 0.05) and 15 miRNAs with significant interaction (P ≤ 0.05) to ICS versus placebo treatments. We selected three miRNAs for functional validation, of which hsa-miR-155-5p and hsa-miR-532-5p were significantly associated with changes in dexamethasone-induced transrepression of NF-κB. Combined, these two miRNAs were predictive of ICS response over the course of the clinical trial, with an area under the receiver operating characteristic curve of 0.86.Conclusions: We identified two functional circulating miRNAs predictive of asthma ICS treatment response over time.

Implementation of point-of-care blood gas testing at a large community hospital: Cost analysis, sepsis bundle compliance, and employee engagement.

INTRODUCTION/BACKGROUND: Point-of-care testing (POCT) platforms support patient-centered approaches to health care delivery and may improve patient care. We evaluated implementation of a POCT platform at a large, acute care hospital in the Midwestern United States. METHODS: We used lactate testing as part of a sepsis bundle protocol to evaluate compliance and mortality outcomes. Respiratory team members were surveyed to assess perception of efficiency, ease of use, timely patient care, and overall engagement with the POCT system. Annualized cost per test of a benchtop analyzer and a POCT platform were compared across 3 years for each platform. RESULTS: Lactate testing volume increased from 61% to 91%, which was associated with improved sepsis bundle protocol compliance. Employees reported high levels of engagement, improvements in efficiency and time savings, and better patient care with POCT. Average cost per test was $10.02 for the benchtop system and $6.21 for the POCT platform. POCT saved our institution $88,476 annually in labor costs. DISCUSSION: Combined with a robust training program emphasizing the use of lactate testing in the context of the overall clinical picture, POCT enabled adherence to the sepsis bundle protocol and may have contributed to lower mortality. Additionally, the COVID-19 pandemic has provided us with unanticipated benefits of using POCT; it has enhanced our ability to deal with stringent infectious disease protocols, saving time and minimizing patient and staff exposure. CONCLUSIONS: Implementation of a POCT platform was associated with improved compliance to our sepsis protocol, reduced sepsis mortality, high employee engagement, and cost savings.

Integrative Genomics Analysis Identifies ACVR1B as a Candidate Causal Gene of Emphysema Distribution.

Genome-wide association studies (GWAS) have identified multiple associations with emphysema apicobasal distribution (EABD), but the biological functions of these variants are unknown. To characterize the functions of EABD-associated variants, we integrated GWAS results with 1) expression quantitative trait loci (eQTL) from the Genotype Tissue Expression (GTEx) project and subjects in the COPDGene (Genetic Epidemiology of COPD) study and 2) cell type epigenomic marks from the Roadmap Epigenomics project. On the basis of these analyses, we selected a variant near ACVR1B (activin A receptor type 1B) for functional validation. SNPs from 168 loci with P values less than 5 × 10-5 in the largest GWAS meta-analysis of EABD were analyzed. Eighty-four loci overlapped eQTL, with 12 of these loci showing greater than 80% likelihood of harboring a single, shared GWAS and eQTL causal variant. Seventeen cell types were enriched for overlap between EABD loci and Roadmap Epigenomics marks (permutation P < 0.05), with the strongest enrichment observed in CD4+, CD8+, and regulatory T cells. We selected a putative causal variant, rs7962469, associated with ACVR1B expression in lung tissue for additional functional investigation, and reporter assays confirmed allele-specific regulatory activity for this variant in human bronchial epithelial and Jurkat immune cell lines. ACVR1B expression levels exhibit a nominally significant association with emphysema distribution. EABD-associated loci are preferentially enriched in regulatory elements of multiple cell types, most notably T-cell subsets. Multiple EABD loci colocalize to regulatory elements that are active across multiple tissues and cell types, and functional analyses confirm the presence of an EABD-associated functional variant that regulates ACVR1B expression, indicating that transforming growth factor-ß signaling plays a role in the EABD phenotype. Clinical trial registered with www.clinicaltrials.gov (NCT00608764).

Interplay of Amino Acid Residues at Positions 28 and 31 in NS5A Defines Resistance Pathways in HCV GT2.

Hepatitis C virus (HCV) genotype (GT) 2 represents approximately 9% of all viral infections globally. While treatment outcomes for GT2-infected patients have improved substantially with direct-acting antiviral agents (DAAs) compared to interferon-α, the presence of polymorphisms in NS5A can impact efficacy of NS5A inhibitor-containing regimens. Thus, pathways of NS5A resistance were explored in GT2 subtypes using elbasvir, an NS5A inhibitor with broad genotype activity. Resistance selection studies, resistance analysis in NS5A-inhibitor treated virologic failures, antiviral activities in replicons bearing a panel of GT2 subtype sequences and amino acid substitutions introduced by site-directed mutagenesis were performed to define determinants of inhibitor susceptibility. Elbasvir showed differential antiviral activity in replicons bearing GT2 sequences. The EC50 values for replicons bearing reference NS5A sequences for GT2a and GT2b were 0.003 and 3.4 nanomolar (nM) respectively. Studies with recombinant replicons demonstrated crosstalk between amino acid positions 28 and 31. The combination of phenylalanine and methionine at positions 28 and 31 respectively, conferred the highest potency reduction for elbasvir in GT2a and GT2b. This combination was observed in failures from the C-SCAPE trial. Addition of grazoprevir, an NS3/4A protease inhibitor, to elbasvir more effectively suppressed the emergence of resistance in GT2 at modest inhibitor concentrations (3X EC90). Ruzasvir, a potent, pan-genotype NS5A inhibitor successfully inhibited replicons bearing GT2 resistance-associated substitutions (RASs) at positions 28 and 31. The studies demonstrate crosstalk between amino acids at positions 28 and 31 in NS5A modulate inhibitor potency and may impact treatment outcomes in some HCV GT2-infected patients.

RNA-sequencing across three matched tissues reveals shared and tissue-specific gene expression and pathway signatures of COPD.

BACKGROUND: Multiple gene expression studies have been performed separately in peripheral blood, lung, and airway tissues to study COPD. We performed RNA-sequencing gene expression profiling of large-airway epithelium, alveolar macrophage and peripheral blood samples from the same subset of COPD cases and controls from the COPDGene study who underwent bronchoscopy at a single center. Using statistical and gene set enrichment approaches, we sought to improve the understanding of COPD by studying gene sets and pathways across these tissues, beyond the individual genomic determinants. METHODS: We performed differential expression analysis using RNA-seq data obtained from 63 samples from 21 COPD cases and controls (includes four non-smokers) via the R package DESeq2. We tested associations between gene expression and variables related to lung function, smoking history, and CT scan measures of emphysema and airway disease. We examined the correlation of differential gene expression across the tissues and phenotypes, hypothesizing that this would reveal preserved and private gene expression signatures. We performed gene set enrichment analyses using curated databases and findings from prior COPD studies to provide biological and disease relevance. RESULTS: The known smoking-related genes CYP1B1 and AHRR were among the top differential expression results for smoking status in the large-airway epithelium data. We observed a significant overlap of genes primarily across large-airway and macrophage results for smoking and airway disease phenotypes. We did not observe specific genes differentially expressed in all three tissues for any of the phenotypes. However, we did observe hemostasis and immune signaling pathways in the overlaps across all three tissues for emphysema, and amyloid and telomere-related pathways for smoking. In peripheral blood, the emphysema results were enriched for B cell related genes previously identified in lung tissue studies. CONCLUSIONS: Our integrative analyses across COPD-relevant tissues and prior studies revealed shared and tissue-specific disease biology. These replicated and novel findings in the airway and peripheral blood have highlighted candidate genes and pathways for COPD pathogenesis.

Blood eosinophil count thresholds and exacerbations in patients with chronic obstructive pulmonary disease.

BACKGROUND: Eosinophilic airway inflammation in patients with chronic obstructive pulmonary disease (COPD) is associated with exacerbations and responsivity to steroids, suggesting potential shared mechanisms with eosinophilic asthma. However, there is no consistent blood eosinophil count that has been used to define the increased exacerbation risk. OBJECTIVE: We sought to investigate blood eosinophil counts associated with exacerbation risk in patients with COPD. METHODS: Blood eosinophil counts and exacerbation risk were analyzed in patients with moderate-to-severe COPD by using 2 independent studies of former and current smokers with longitudinal data. The Genetic Epidemiology of COPD (COPDGene) study was analyzed for discovery (n = 1,553), and the Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints (ECLIPSE) study was analyzed for validation (n = 1,895). A subset of the ECLIPSE study subjects were used to assess the stability of blood eosinophil counts over time. RESULTS: COPD exacerbation risk increased with higher eosinophil counts. An eosinophil count threshold of 300 cells/µL or greater showed adjusted incidence rate ratios for exacerbations of 1.32 in the COPDGene study (95% CI, 1.10-1.63). The cutoff of 300 cells/µL or greater was validated for prospective risk of exacerbation in the ECLIPSE study, with adjusted incidence rate ratios of 1.22 (95% CI, 1.06-1.41) using 3-year follow-up data. Stratified analysis confirmed that the increased exacerbation risk associated with an eosinophil count of 300 cells/µL or greater was driven by subjects with a history of frequent exacerbations in both the COPDGene and ECLIPSE studies. CONCLUSIONS: Patients with moderate-to-severe COPD and blood eosinophil counts of 300 cells/µL or greater had an increased risk exacerbations in the COPDGene study, which was prospectively validated in the ECLIPSE study.

Unraveling the structural basis of grazoprevir potency against clinically relevant substitutions in hepatitis C virus NS3/4A protease from genotype 1a.

Grazoprevir is a potent pan-genotype and macrocyclic inhibitor of hepatitis C virus (HCV) NS3/4A protease and was developed for treating chronic HCV infection. In HCV genotype (GT) 1a, grazoprevir maintains potent activity against a majority of NS3 resistance-associated amino acid substitutions, including the highly prevalent and naturally occurring Q80K polymorphism that impacts simeprevir, another NS3/4A protease inhibitor. The basis for an unexpected difference in the clinical impact of some NS3 substitutions was investigated. Phenotypic analysis of resistance-associated substitutions identified in NS3 from GT1a-infected patients who failed therapy with grazoprevir (in combination with elbasvir, an inhibitor of HCV NS5A protein) showed that positions 56, 156, and 168 in NS3 were most impactful because they diminished protein-inhibitor interactions. Although an amino acid substitution from aspartic acid to alanine at position 168 (D168A) reduced the potency of grazoprevir, its combination with R155K unexpectedly nullified this effect. Molecular dynamics and free-energy surface studies indicated that Asp-168 is important in anchoring Arg-155 for ligand binding but is not critical for Lys-155 because of the inherent flexibility of its side chain. Moreover, modeling studies supported a strong direct cation-heterocycle interaction between the Lys-155 side chain of the double substitution, R155K/D168A, and the lone pair on the quinoxaline in grazoprevir. This unique interaction provides a structural basis for grazoprevir's higher potency than simeprevir, an inhibitor to which the double substitution confers a significant reduction in potency. Our findings are consistent with the detection of R155K/D168A in NS3 from virologic failures treated with simeprevir but not grazoprevir.

In Vitro Antiviral Profile of Ruzasvir, a Potent and Pangenotype Inhibitor of Hepatitis C Virus NS5A.

Inhibition of NS5A has emerged as an attractive strategy to intervene in hepatitis C virus (HCV) replication. Ruzasvir (formerly MK-8408) was developed as a novel NS5A inhibitor to improve upon the potency and barrier to resistance of early compounds. Ruzasvir inhibited HCV RNA replication with 50% effective concentrations (EC50s) of 1 to 4 pM in Huh7 or Huh7.5 cells bearing replicons for HCV genotype 1 (GT1) to GT7. The antiviral activity was modestly (10-fold) reduced in the presence of 40% normal human serum. The picomolar potency in replicon cells extended to sequences of clinical isolates available in public databases that were synthesized and tested as replicons. In GT1a, ruzasvir inhibited common NS5A resistance-associated substitutions (RASs), with the exception of M28G. De novo resistance selection studies identified pathways with certain amino acid substitutions at residues 28, 30, 31, and 93 across genotypes. Substitutions at position 93 were more common in GT1 to -4, while changes at position 31 emerged frequently in GT5 and -6. With the exception of GT4, the reintroduction of selected RASs conferred a ≥100-fold potency reduction in the antiviral activity of ruzasvir. Common RASs from other classes of direct-acting antiviral agents (DAAs) did not confer cross-resistance to ruzasvir. The interaction of ruzasvir with an NS3/4A protease inhibitor (grazoprevir) and an NS5B polymerase prodrug (uprifosbuvir) was additive to synergistic, with no evidence of antagonism or cytotoxicity. The antiviral profile of ruzasvir supported its further evaluation in human trials in combination with grazoprevir and uprifosbuvir.

Elbasvir/Grazoprevir for Patients With Hepatitis C Virus Infection and Inherited Blood Disorders: A Phase III Study.

Direct-acting antiviral agents have not been studied exclusively in patients with inherited blood disorders and hepatitis C virus (HCV) infection. The objective of the randomized, placebo-controlled, phase III C-EDGE IBLD study was to assess the safety and efficacy of elbasvir/grazoprevir (EBR/GZR) in patients with inherited bleeding disorders and HCV infection. One hundred fifty-nine adults with HCV infection and sickle cell anemia, thalassemia, or hemophilia A/B or von Willebrand disease were enrolled at 31 study sites in the United States, Europe, Australia, Canada, Israel, and Thailand. Patients were given an oral, once-daily, fixed-dose combination of EBR/GZR 50 mg/100 mg for 12 weeks and randomized to the immediate-treatment group (ITG) or deferred-treatment group (DTG; placebo followed by active treatment). The primary endpoints were the proportion of patients in the ITG with unquantifiable HCV RNA 12 weeks posttreatment (sustained virological response 12 weeks after completion of study treatment; SVR12) and the comparison of safety in the ITG and DTG. In the ITG, 100 of 107 patients (93.5%) achieved SVR12, 6 relapsed, and 1 was lost to follow-up. SVR12 was achieved in 94.7% (18 of 19), 97.6% (40 of 41), and 89.4% (42 of 47) of patients with sickle cell disease, ß-thalassemia, and hemophilia A/B or von Willebrand disease, respectively. Serious adverse events were reported by 2.8% (n = 3) and 11.5% (n = 6) of patients in the ITG and DTG, respectively. Hemoglobin levels and international normalized ratio values were similar in patients receiving EBR/GZR and placebo; among patients with hemoglobinopathies, change in mean hemoglobin levels was similar in those receiving EBR/GZR compared to those receiving placebo. CONCLUSION: These results add to the expanding pool of data available for EBR/GZR, indicating a high level of efficacy and favorable tolerability in patients with HCV infection. (Hepatology 2017;66:736-745).

Antiviral Activity and Resistance Analysis of NS3/4A Protease Inhibitor Grazoprevir and NS5A Inhibitor Elbasvir in Hepatitis C Virus GT4 Replicons.

Although genotype 4 (GT4)-infected patients represent a minor overall percentage of the global hepatitis C virus (HCV)-infected population, the high prevalence of the genotype in specific geographic regions coupled with substantial sequence diversity makes it an important genotype to study for antiviral drug discovery and development. We evaluated two direct-acting antiviral agents-grazoprevir, an HCV NS3/4A protease inhibitor, and elbasvir, an HCV NS5A inhibitor-in GT4 replicons prior to clinical studies in this genotype. Following a bioinformatics analysis of available GT4 sequences, a set of replicons bearing representative GT4 clinical isolates was generated. For grazoprevir, the 50% effective concentration (EC50) against the replicon bearing the reference GT4a (ED43) NS3 protease and NS4A was 0.7 nM. The median EC50 for grazoprevir against chimeric replicons encoding NS3/4A sequences from GT4 clinical isolates was 0.2 nM (range, 0.11 to 0.33 nM; n = 5). The difficulty in establishing replicons bearing NS3/4A resistance-associated substitutions was substantially overcome with the identification of a G162R adaptive substitution in NS3. Single NS3 substitutions D168A/V identified from de novo resistance selection studies reduced grazoprevir antiviral activity by 137- and 47-fold, respectively, in the background of the G162R replicon. For elbasvir, the EC50 against the replicon bearing the reference full-length GT4a (ED43) NS5A gene was 0.0002 nM. The median EC50 for elbasvir against chimeric replicons bearing clinical isolates from GT4 was 0.0007 nM (range, 0.0002 to 34 nM; n = 14). De novo resistance selection studies in GT4 demonstrated a high propensity to suppress the emergence of amino acid substitutions that confer high-potency reductions to elbasvir. Phenotypic characterization of the NS5A amino acid substitutions identified (L30F, L30S, M31V, and Y93H) indicated that they conferred 15-, 4-, 2.5-, and 7.5-fold potency losses, respectively, to elbasvir. The activity profiles of grazoprevir and elbasvir supported the testing of the direct-acting antivirals in clinical studies.