Age-associated Differences in the Human Lung Extracellular Matrix

Extracellular matrix (ECM) remodeling has been associated with chronic lung diseases. However, information about specific age‑associated differences in lung ECM is currently limited. In this study we aimed to identify and localize age‑associated ECM differences in human lung using comprehensive transcriptomic, proteomic and immunohistochemical analyses. Our previously identified age-associated gene expression signature of the lung was re‑analyzed limiting it to an aging signature based on 270 control patients (37-80 years) and focused on the Matrisome core geneset using geneset enrichment analysis. To validate the age‑associated transcriptomic differences on protein level, we compared the age-associated ECM genes (FDR <0.05) with a profile of age‑associated proteins identified from a lung tissue proteomics dataset from 9 control patients (49-76 years) (FDR<0.05). Extensive immunohistochemical analysis was used to localize and semi-quantify the age-associated ECM differences in lung tissues from 62 control patients (18-82 years). Comparative analysis of transcriptomic and proteomic data identified 7 ECM proteins with higher expression with age at both gene and protein level: COL1A1, COL6A1, COL6A2, COL14A1, FBLN2, LTBP4 and LUM. With immunohistochemistry we demonstrated higher protein level with age for COL6A2 in whole tissue, parenchyma, airway wall and blood vessel, for COL14A1 and LUM in bronchial epithelium, and COL1A1 in lung parenchyma. Our study revealed that higher age is associated with lung ECM remodeling, with specific differences occurring in defined regions within the lung. These differences may affect lung structure and physiology with aging and as such may increase susceptibility for developing chronic lung diseases.

Genetic Associations and Architecture of Asthma-COPD Overlap

BACKGROUND

Some people have characteristics of both asthma and COPD (asthma-COPD overlap), and evidence suggests they experience worse outcomes than those with either condition alone.

RESEARCH QUESTION

What is the genetic architecture of asthma-COPD overlap, and do the determinants of risk for asthma-COPD overlap differ from those for COPD or asthma?

STUDY DESIGN AND METHODS

We conducted a genome-wide association study in 8,068 asthma-COPD overlap case subjects and 40,360 control subjects without asthma or COPD of European ancestry in UK Biobank (stage 1). We followed up promising signals (P < 5 × 10^-6) that remained associated in analyses comparing (1) asthma-COPD overlap vs asthma-only control subjects, and (2) asthma-COPD overlap vs COPD-only control subjects. These variants were analyzed in 12 independent cohorts (stage 2).

RESULTS

We selected 31 independent variants for further investigation in stage 2, and discovered eight novel signals (P < 5 × 10^-8) for asthma-COPD overlap (meta-analysis of stage 1 and 2 studies). These signals suggest a spectrum of shared genetic influences, some predominantly influencing asthma (FAM105A, GLB1, PHB, TSLP), others predominantly influencing fixed airflow obstruction (IL17RD, C5orf56, HLA-DQB1). One intergenic signal on chromosome 5 had not been previously associated with asthma, COPD, or lung function. Subgroup analyses suggested that associations at these eight signals were not driven by smoking or age at asthma diagnosis, and in phenome-wide scans, eosinophil counts, atopy, and asthma traits were prominent.

INTERPRETATION

We identified eight signals for asthma-COPD overlap, which may represent loci that predispose to type 2 inflammation, and serious long-term consequences of asthma.

Multi-omics highlights ABO plasma protein as a causal risk factor for COVID-19

SARS-CoV-2 is responsible for the coronavirus disease 2019 (COVID-19) and the current health crisis. Despite intensive research efforts, the genes and pathways that contribute to COVID-19 remain poorly understood. We, therefore, used an integrative genomics (IG) approach to identify candidate genes responsible for COVID-19 and its severity. We used Bayesian colocalization (COLOC) and summary-based Mendelian randomization to combine gene expression quantitative trait loci (eQTLs) from the Lung eQTL (n = 1,038) and eQTLGen (n = 31,784) studies with published COVID-19 genome-wide association study (GWAS) data from the COVID-19 Host Genetics Initiative. Additionally, we used COLOC to integrate plasma protein quantitative trait loci (pQTL) from the INTERVAL study (n = 3,301) with COVID-19 loci. Finally, we determined any causal associations between plasma proteins and COVID-19 using multi-variable two-sample Mendelian randomization (MR). The expression of 18 genes in lung and/or blood co-localized with COVID-19 loci. Of these, 12 genes were in suggestive loci (PGWAS < 5 × 10-05). LZTFL1, SLC6A20, ABO, IL10RB and IFNAR2 and OAS1 had been previously associated with a heightened risk of COVID-19 (PGWAS < 5 × 10-08). We identified a causal association between OAS1 and COVID-19 GWAS. Plasma ABO protein, which is associated with blood type in humans, demonstrated a significant causal relationship with COVID-19 in the MR analysis; increased plasma levels were associated with an increased risk of COVID-19 and, in particular, severe COVID-19. In summary, our study identified genes associated with COVID-19 that may be prioritized for future investigations. Importantly, this is the first study to demonstrate a causal association between plasma ABO protein and COVID-19.

Detection of respiratory syncytial virus defective genomes in nasal secretions is associated with distinct clinical outcomes

Respiratory syncytial virus (RSV) causes respiratory illness in children, immunosuppressed individuals and the elderly. However, the viral factors influencing the clinical outcome of RSV infections remain poorly defined. Defective viral genomes (DVGs) can suppress virus replication by competing for viral proteins and by stimulating antiviral immunity. We studied the association between detection of DVGs of the copy-back type and disease severity in three RSV A-confirmed cohorts. In hospitalized children, detection of DVGs in respiratory samples at or around the time of admission associated strongly with more severe disease, higher viral load and a stronger pro-inflammatory response. Interestingly, in experimentally infected adults, the presence of DVGs in respiratory secretions differentially associated with RSV disease severity depending on when DVGs were detected. Detection of DVGs early after infection associated with low viral loads and mild disease, whereas detection of DVGs late after infection, especially if DVGs were present for prolonged periods, associated with high viral loads and severe disease. Taken together, we demonstrate that the kinetics of DVG accumulation and duration could predict clinical outcome of RSV A infection in humans, and thus could be used as a prognostic tool to identify patients at risk of worse clinical disease.

Letermovir Resistance Analysis in a Clinical Trial of Cytomegalovirus Prophylaxis for Hematopoietic Stem Cell Transplant Recipients

BACKGROUND

Letermovir (LET), a cytomegalovirus (CMV) deoxyribonucleic acid (DNA) terminase inhibitor, was recently approved for prophylaxis of CMV infection in adult CMV-seropositive recipients of allogeneic hematopoietic stem cell transplantation. Cytomegalovirus genotyping was performed to identify LET-resistance-associated variants (RAVs) among subjects in a Phase 3 trial.

METHODS

The CMV UL56 and UL89 genes, encoding subunits of CMV DNA terminase, were sequenced from plasma collected from subjects with clinically significant CMV infection (CS-CMVi). Novel variants were evaluated by recombinant phenotyping to assess their potential to confer resistance to LET.

RESULTS

Genotyping was successful for 50 of 79 LET subjects with CS-CMVi. Resistance-associated variants (encoding pUL56 V236M and C325W) were detected independently in subjects 1 and 3 who experienced CS-CMVi while receiving LET prophylaxis, and 2 other variants (encoding pUL56 E237G and R369T) were detected >3 weeks after subjects 2 and 3, respectively, had discontinued LET prophylaxis and received preemptive therapy with ganciclovir.

CONCLUSIONS

The detected incidence of CMV resistance among subjects who received LET as prophylaxis in this Phase 3 trial was low. The LET RAVs that were detected mapped to the CMV UL56 gene at positions associated with reduced susceptibility to LET based on resistance selections in cell culture.

Gene expression network analysis provides potential targets against SARS-CoV-2

Cell entry of SARS-CoV-2, the novel coronavirus causing COVID-19, is facilitated by host cell angiotensin-converting enzyme 2 (ACE2) and transmembrane serine protease 2 (TMPRSS2). We aimed to identify and characterize genes that are co-expressed with ACE2 and TMPRSS2, and to further explore their biological functions and potential as druggable targets. Using the gene expression profiles of 1,038 lung tissue samples, we performed a weighted gene correlation network analysis (WGCNA) to identify modules of co-expressed genes. We explored the biology of co-expressed genes using bioinformatics databases, and identified known drug-gene interactions. ACE2 was in a module of 681 co-expressed genes; 10 genes with moderate-high correlation with ACE2 (r > 0.3, FDR < 0.05) had known interactions with existing drug compounds. TMPRSS2 was in a module of 1,086 co-expressed genes; 31 of these genes were enriched in the gene ontology biologic process 'receptor-mediated endocytosis', and 52 TMPRSS2-correlated genes had known interactions with drug compounds. Dozens of genes are co-expressed with ACE2 and TMPRSS2, many of which have plausible links to COVID-19 pathophysiology. Many of the co-expressed genes are potentially targetable with existing drugs, which may accelerate the development of COVID-19 therapeutics.

Rhinovirus and asthma: Challenges and opportunities

Human rhinoviruses (RVs) are the primary aetiological agent of the common cold. Generally, the associated infection is mild and self‐limiting, but may also be associated with bronchiolitis in infants, pneumonia in the immunocompromised and exacerbation in patients with pulmonary conditions such as asthma or chronic obstructive pulmonary disease. Viral infection accounts for as many as two thirds of asthma exacerbations in children and more than half in adults. Allergy and asthma are major risk factors for more frequent and severe RV‐related illnesses. The prevalence of RV‐induced wheezing will likely continue to increase given that asthma affects a significant proportion of the population, with allergic asthma accounting for the majority. Several new respiratory viruses and their subgroups have been discovered, with various degrees of relevance. This review will focus on RV infection in the context of the epidemiologic evidence, genetic variability, pathobiology, clinical studies in the context of asthma, differences with other viruses including COVID‐19 and current treatment interventions.

Neutrophilic inflammation in the respiratory mucosa predisposes to RSV infection

The variable outcome of viral exposure is only partially explained by known factors. We administered respiratory syncytial virus (RSV) to 58 volunteers, of whom 57% became infected. Mucosal neutrophil activation before exposure was highly predictive of symptomatic RSV disease. This was associated with a rapid, presymptomatic decline in mucosal interleukin-17A (IL-17A) and other mediators. Conversely, those who resisted infection showed presymptomatic activation of IL-17- and tumor necrosis factor-related pathways. Vulnerability to infection was not associated with baseline microbiome but was reproduced in mice by preinfection chemokine-driven airway recruitment of neutrophils, which caused enhanced disease mediated by pulmonary CD8+ T cell infiltration. Thus, mucosal neutrophilic inflammation at the time of RSV exposure enhances susceptibility, revealing dynamic, time-dependent local immune responses before symptom onset and explaining the as-yet unpredictable outcomes of pathogen exposure.

Abstract CT247: A phase 1/2 study of GB1275, a first-in-class CD11b modulator, as monotherapy and with an anti-PD-1 antibody in specified advanced solid tumors or with chemotherapy in metastatic pancreatic cancer (KEYNOTE-A36)

Background: Tumor influx of CD11b-expressing myeloid-derived suppressor cells (MDSCs) and M2 tumor-associated macrophages (TAMs) creates an immunosuppressive tumor microenvironment associated with resistance to anti-PD-1 antibody therapy. GB1275 is a novel, first-in-class, CD11b modulator that, in vivo, reduced MDSCs and TAMs at the tumor site, repolarized M2 immunosuppressive TAMs to an M1 phenotype, and increased tumor infiltration of activated CD8+ T cells. In combination with an anti-PD-1 antibody or chemotherapy, these immunomodulatory effects translated into potent anti-tumor effects and prolonged survival in orthotopic PDAC models. We hypothesize that GB1275-based therapy can alleviate myeloid cell-mediated immunosuppression and improve cancer treatment outcomes. Methods: This is an open-label, first-in-human study comprising phase 1 dose escalation of GB1275 monotherapy (Regimen A) and GB1275 plus pembrolizumab (Regimen B) in patients with previously treated, locally advanced or metastatic PDA, esophageal, gastric/GEJ, triple negative breast, castration-resistant prostate, or microsatellite-stable colorectal cancer (MSS CRC), and GB1275 plus nab-paclitaxel and gemcitabine (Nab-P+Gem) (Regimen C) in mPDAC. This is followed by phase 2 expansion in three disease cohorts: 1) newly diagnosed stage IV mPDAC (GB1275+Nab-P+Gem), 2) MSS CRC (GB1275+pembrolizumab), and 3) PD-L1+ gastric/GEJ cancer (GB1275+pembrolizumab). Patients are enrolled in Regimen A initially, with Regimen B commencing after completion of the first few cohorts of Regimen A. Regimen C will initiate when Regimen A is completed. Eligible patients are ≥18 years of age, with histologically confirmed locally advanced/metastatic tumor specified and ECOG 0-1 PS; prior immunotherapy is permissible during Regimen A and B dose escalation, but not Regimen C or during the expansion phase. Patients are excluded for untreated or symptomatic CNS metastasis, prior myeloid targeting treatment or other prohibited medications, or a history of clinically significant cardiovascular disease. Patients with active autoimmune disease requiring systemic therapy will be excluded from Regimen B and from expansion cohorts 2 and 3. Primary objectives for phase 1 are to determine the maximum tolerated dose/recommended phase dose and pharmacokinetic profile of GB1275 alone and combination with pembrolizumab, and safety in combination with Nab-P+Gem, and for phase 2, to assess efficacy.Statistical Considerations: A 3+3 design is used for the dose escalation phase and Simon's 2-stage design for the expansion phase. Adverse events are graded per CTCAE v5.0, responses per RECIST v1.1. This study is open for recruitment (NCT04060342).

Citation Format: Andrea Wang-Gillam, Drew W. Rasco, Wungki Park, Eileen O'Reilly, Wells Messersmith, David G. DeNardo, Vineet Gupta, Lei Zhou, Anna Galkin, Debbie Slee, Laura L. Carter, David Nickle, Rebecca Tran, Jack Li, Beatrice Ferguson, Marya F. Chaney, Luisa Salter-Cid, Jakob Dupont, Johanna C. Bendell. A phase 1/2 study of GB1275, a first-in-class CD11b modulator, as monotherapy and with an anti-PD-1 antibody in specified advanced solid tumors or with chemotherapy in metastatic pancreatic cancer (KEYNOTE-A36) [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr CT247.

A phase I/II study of GB1275, a first-in-class oral CD11b modulator, alone, and combined with pembrolizumab in specified advanced solid tumors or with chemotherapy in metastatic pancreatic cancer (KEYNOTE-A36)

3085

Background: GB1275 is a first-in-class CD11b modulator that reduced myeloid-derived suppressor cells (MDSCs) and tumor associated macrophages (TAMs) at the tumor site, repolarized M2 immunosuppressive TAMs to an M1 phenotype, and increased tumor infiltration of activated CD8+ T cells in preclinical models. When combined with an anti-PD-1 antibody or chemotherapy, these immunomodulatory effects translated into potent anti-tumor effects and prolonged survival in orthotopic PDAC models [Panni RZ, et al. Sci Transl Med. 2019 Jul 3;11(499)]. This ongoing first-in-human study consists of dose escalation of GB1275 monotherapy (Regimen A), GB1275 + pembrolizumab (Regimen B), and GB1275 + nab-paclitaxel + gemcitabine (Regimen C), followed by Phase 2 expansion in newly diagnosed metastatic pancreatic, MSS colorectal, and PD-L1-positive gastric/GEJ cancers. Here we report interim results of the dose escalation portion of the trial. Methods: The dose escalation phase is based on a standard oncology phase 1, 3+3 design. Cohorts of 3 to 6 patients (pts) with histologically confirmed locally advanced/metastatic pancreatic, esophageal, gastric, MSS colorectal, prostate, or breast cancer were sequentially assigned to ascending dose levels of GB1275 taken orally twice daily (BID) in 1 of 3 regimens: Regimen A was initiated first; Regimen B commenced after completion of the first two cohorts of Regimen A, and Regimen C will be initiated when Regimen A is completed. Dose escalation was based on assessment of safety including dose-limiting toxicity (DLT). Serial blood and tumor samples were collected for pharmacokinetic (PK) and biomarker analyses. Results: As of January 21, 2020, 13 pts were treated, with 3 each in Regimen A (GB1275 100mg, 200 mg and 400 mg BID) dose levels and 4 in Regimen B with GB1275 100 mg BID + pembrolizumab. No DLTs have been reported. GB1275 treatment-related adverse events were reported in 5 pts; all were Grade 1 in severity. Preliminary PK analyses showed a mean elimination half-life of ~7 hours. Reduction in peripheral MDSCs was observed in the majority of pts with serial samples. Biomarker analysis in serial tumor tissue is ongoing. Conclusions: Preliminary data show minimal treatment-related toxicities with the studied regimens. PK data support BID dosing. Dose escalation is ongoing. Updated data will be presented. Clinical trial information: NCT04060342 .

Hedgehog signaling in the airway epithelium of patients with chronic obstructive pulmonary disease

Genome-wide association studies have linked gene variants of the receptor patched homolog 1 (PTCH1) with chronic obstructive pulmonary disease (COPD). However, its biological role in the disease is unclear. Our objective was to determine the expression pattern and biological role of PTCH1 in the lungs of patients with COPD. Airway epithelial-specific PTCH1 protein expression and epithelial morphology were assessed in lung tissues of control and COPD patients. PTCH1 mRNA expression was measured in bronchial epithelial cells obtained from individuals with and without COPD. The effects of PTCH1 siRNA knockdown on epithelial repair and mucous expression were evaluated using human epithelial cell lines. Ptch1^+/− mice were used to assess the effect of decreased PTCH1 on mucous expression and airway epithelial phenotypes. Airway epithelial-specific PTCH1 protein expression was significantly increased in subjects with COPD compared to controls, and its expression was associated with total airway epithelial cell count and thickness. PTCH1 knockdown attenuated wound closure and mucous expression in airway epithelial cell lines. Ptch1^+/− mice had reduced mucous expression compared to wildtype mice following mucous induction. PTCH1 protein is up-regulated in COPD airway epithelium and may upregulate mucous expression. PTCH1 provides a novel target to reduce chronic bronchitis in COPD patients.

Leveraging lung tissue transcriptome to uncover candidate causal genes in COPD genetic associations

Causal genes of chronic obstructive pulmonary disease (COPD) remain elusive. The current study aims at integrating genome-wide association studies (GWAS) and lung expression quantitative trait loci (eQTL) data to map COPD candidate causal genes and gain biological insights into the recently discovered COPD susceptibility loci. Two complementary genomic datasets on COPD were studied. First, the lung eQTL dataset which included whole-genome gene expression and genotyping data from 1038 individuals. Second, the largest COPD GWAS to date from the International COPD Genetics Consortium (ICGC) with 13 710 cases and 38 062 controls. Methods that integrated GWAS with eQTL signals including transcriptome-wide association study (TWAS), colocalization and Mendelian randomization-based (SMR) approaches were used to map causality genes, i.e. genes with the strongest evidence of being the functional effector at specific loci. These methods were applied at the genome-wide level and at COPD risk loci derived from the GWAS literature. Replication was performed using lung data from GTEx. We collated 129 non-overlapping risk loci for COPD from the GWAS literature. At the genome-wide scale, 12 new COPD candidate genes/loci were revealed and six replicated in GTEx including CAMK2A, DMPK, MYO15A, TNFRSF10A, BTN3A2 and TRBV30. In addition, we mapped candidate causal genes for 60 out of the 129 GWAS-nominated loci and 23 of them were replicated in GTEx. Mapping candidate causal genes in lung tissue represents an important contribution to the genetics of COPD, enriches our biological interpretation of GWAS findings, and brings us closer to clinical translation of genetic associations.

Genome-wide interaction study of gene-by-occupational exposures on respiratory symptoms

Respiratory symptoms are important indicators of respiratory diseases. Both genetic and environmental factors contribute to respiratory symptoms development but less is known about gene-environment interactions. We aimed to assess interactions between single nucleotide polymorphisms (SNPs) and occupational exposures on respiratory symptoms cough, dyspnea and phlegm. As identification cohort LifeLines I (n = 7976 subjects) was used. Job-specific exposure was estimated using the ALOHA + job exposure matrix. SNP-by-occupational exposure interactions on respiratory symptoms were tested using logistic regression adjusted for gender, age, and current smoking. SNP-by-exposure interactions with a p-value <10^-4 were tested for replication in two independent cohorts: LifeLines II (n = 5260) and the Vlagtwedde-Vlaardingen cohort (n = 1529). The interaction estimates of the replication cohorts were meta-analyzed using PLINK. Replication was achieved when the meta-analysis p-value was <0.05 and the interaction effect had the same direction as in the identification cohort. Additionally, we assessed whether replicated SNPs associated with gene expression by analyzing if they were cis-acting expression quantitative trait loci (eQTL) in lung tissue. In the replication meta-analysis, sixteen out of 477 identified SNP-by-occupational exposure interactions had a p-value <0.05 and 9 of these interactions had the same direction as in the identification cohort. Several identified loci were plausible candidates for respiratory symptoms, such as TMPRSS9, SERPINH1, TOX3, and ARHGAP18. Three replicated SNPs were cis-eQTLs for FCER1A, CHN1, and TIMM13 in lung tissue. Taken together, this genome-wide SNP-by-occupational exposure interaction study in relation to cough, dyspnea, and phlegm identified several suggestive susceptibility genes. Further research should determine if these genes are true susceptibility loci for respiratory symptoms in relation to occupational exposures.

Genome-Wide Interaction Analysis of Air Pollution Exposure and Childhood Asthma with Functional Follow-up

RATIONALE

The evidence supporting an association between traffic-related air pollution exposure and incident childhood asthma is inconsistent and may depend on genetic factors.

OBJECTIVES

To identify gene-environment interaction effects on childhood asthma using genome-wide single-nucleotide polymorphism (SNP) data and air pollution exposure. Identified loci were further analyzed at epigenetic and transcriptomic levels.

METHODS

We used land use regression models to estimate individual air pollution exposure (represented by outdoor NO₂ levels) at the birth address and performed a genome-wide interaction study for doctors' diagnoses of asthma up to 8 years in three European birth cohorts (n = 1,534) with look-up for interaction in two separate North American cohorts, CHS (Children's Health Study) and CAPPS/SAGE (Canadian Asthma Primary Prevention Study/Study of Asthma, Genetics and Environment) (n = 1,602 and 186 subjects, respectively). We assessed expression quantitative trait locus effects in human lung specimens and blood, as well as associations among air pollution exposure, methylation, and transcriptomic patterns.

MEASUREMENTS AND MAIN RESULTS

In the European cohorts, 186 SNPs had an interaction P < 1 × 10^-4 and a look-up evaluation of these disclosed 8 SNPs in 4 loci, with an interaction P < 0.05 in the large CHS study, but not in CAPPS/SAGE. Three SNPs within adenylate cyclase 2 (ADCY2) showed the same direction of the interaction effect and were found to influence ADCY2 gene expression in peripheral blood (P = 4.50 × 10^-4). One other SNP with P < 0.05 for interaction in CHS, rs686237, strongly influenced UDP-Gal:betaGlcNAc β-1,4-galactosyltransferase, polypeptide 5 (B4GALT5) expression in lung tissue (P = 1.18 × 10^-17). Air pollution exposure was associated with differential discs, large homolog 2 (DLG2) methylation and expression.

CONCLUSIONS

Our results indicated that gene-environment interactions are important for asthma development and provided supportive evidence for interaction with air pollution for ADCY2, B4GALT5, and DLG2.

Sulfatase modifying factor 1 (SUMF1) is associated with Chronic Obstructive Pulmonary Disease

Background

It has been observed that mice lacking the sulfatase modifying factor (Sumf1) developed an emphysema-like phenotype. However, it is unknown if SUMF1 may play a role in Chronic Obstructive Pulmonary Disease (COPD) in humans. The aim was to investigate if the expression and genetic regulation of SUMF1 differs between smokers with and without COPD.

Methods

SUMF1 mRNA was investigated in sputum cells and whole blood from controls and COPD patients (all current or former smokers). Expression quantitative trait loci (eQTL) analysis was used to investigate if single nucleotide polymorphisms (SNPs) in SUMF1 were significantly associated with SUMF1 expression. The association of SUMF1 SNPs with COPD was examined in a population based cohort, Lifelines. SUMF1 mRNA from sputum cells, lung tissue, and lung fibroblasts, as well as lung function parameters, were investigated in relation to genotype.

Results

Certain splice variants of SUMF1 showed a relatively high expression in lung tissue compared to many other tissues. SUMF1 Splice variant 2 and 3 showed lower levels in sputum cells from COPD patients as compared to controls. Twelve SNPs were found significant by eQTL analysis and overlapped with the array used for genotyping of Lifelines. We found alterations in mRNA expression in sputum cells and lung fibroblasts associated with SNP rs11915920 (top hit in eQTL), which validated the results of the lung tissue eQTL analysis. Of the twelve SNPs, two SNPs, rs793391 and rs308739, were found to be associated with COPD in Lifelines. The SNP rs793391 was also confirmed to be associated with lung function changes.

Conclusions

We show that SUMF1 expression is affected in COPD patients compared to controls, and that SNPs in SUMF1 are associated with an increased risk of COPD. Certain COPD-associated SNPs have effects on either SUMF1 gene expression or on lung function. Collectively, this study shows that SUMF1 is associated with an increased risk of developing COPD.

Electronic supplementary material

The online version of this article (doi:10.1186/s12931-017-0562-5) contains supplementary material, which is available to authorized users.

Human Lung Tissue Transcriptome: Influence of Sex and Age

Background

Sex and age strongly influence the pathophysiology of human lungs, but scarce information is available about their effects on pulmonary gene expression.

Methods

We followed a discovery-validation strategy to identify sex- and age-related transcriptional differences in lung.

Results

We identified transcriptional profiles significantly associated with sex (215 genes; FDR < 0.05) and age at surgery (217 genes) in non-involved lung tissue resected from 284 lung adenocarcinoma patients. When these profiles were tested in three independent series of non-tumor lung tissue from an additional 1,111 patients, we validated the association with sex and age for 25 and 22 genes, respectively. Among the 17 sex-biased genes mapping on chromosome X, 16 have been reported to escape X-chromosome inactivation in other tissues or cells, suggesting that this mechanism influences lung transcription too. Our 22 age-related genes partially overlap with genes modulated by age in other tissues, suggesting that the aging process has similar consequences on gene expression in different organs. Finally, seven genes whose expression was modulated by sex in non-tumor lung tissue, but no age-related gene, were also validated using publicly available data from 990 lung adenocarcinoma samples, suggesting that the physiological regulatory mechanisms are only partially active in neoplastic tissue.

Conclusions

Gene expression in non-tumor lung tissue is modulated by both sex and age. These findings represent a validated starting point for research on the molecular mechanisms underlying the observed differences in the course of lung diseases among men and women of different ages.

Cross-Cancer Genome-Wide Analysis of Lung, Ovary, Breast, Prostate, and Colorectal Cancer Reveals Novel Pleiotropic Associations

Identifying genetic variants with pleiotropic associations can uncover common pathways influencing multiple cancers. We took a two-stage approach to conduct genome-wide association studies for lung, ovary, breast, prostate, and colorectal cancer from the GAME-ON/GECCO Network (61,851 cases, 61,820 controls) to identify pleiotropic loci. Findings were replicated in independent association studies (55,789 cases, 330,490 controls). We identified a novel pleiotropic association at 1q22 involving breast and lung squamous cell carcinoma, with eQTL analysis showing an association with ADAM15/THBS3 gene expression in lung. We also identified a known breast cancer locus CASP8/ALS2CR12 associated with prostate cancer, a known cancer locus at CDKN2B-AS1 with different variants associated with lung adenocarcinoma and prostate cancer, and confirmed the associations of a breast BRCA2 locus with lung and serous ovarian cancer. This is the largest study to date examining pleiotropy across multiple cancer-associated loci, identifying common mechanisms of cancer development and progression. Cancer Res; 76(17); 5103-14. ©2016 AACR.

Virologic resistance analysis from a phase 2 study of MK-5172 combined with pegylated interferon/ribavirin in treatment-naive patients with hepatitis C virus genotype 1 infection

BACKGROUND

Virologic failure following treatment of hepatitis C virus (HCV) genotype 1 with direct-acting antiviral agents is often accompanied by the emergence of resistant variants. MK-5172 is an investigational once-daily protease inhibitor. We analyzed variants in treatment-naive noncirrhotic patients with virologic failure on MK-5172 (100-800 mg/day) plus pegylated interferon alfa/ribavirin (peg-IFN/RBV) during a phase 2 trial.

METHODS

Population and selective clonal sequencing were performed at baseline and at virologic failure in the 4 MK-5172 dosing arms. MK-5172 activity was determined using a mutant replicon assay.

RESULTS

Six of 266 (2.3%) MK-5172 recipients satisfied prespecified criteria for virologic failure, all with genotype 1a infection. Five patients with virologic failure were in the MK-5172 100-mg arm, including 4 patients with low plasma MK-5172 levels documented during triple therapy. Variants associated with >4-fold loss of potency were detected in 3 of the 4 patients with genotype 1a breakthrough while on MK-5172. The fifth patient had undetectable HCV-RNA levels at the end of triple therapy but subsequently broke through during the peg-IFN/RBV tail 16 weeks after completion of MK-5172. Three patients had D168 variants at virologic failure, including 2 with the D168A variant associated with a 95-fold loss of potency. The sole apparent relapse was actually a genotype 3a reinfection in the MK-5172 200-mg group.

CONCLUSIONS

Virologic failure occurred uncommonly (6/266 [2.3%]) in MK-5172/peg-IFN/RBV recipients. The most prevalent treatment-emergent variants were detected at the D168 locus. D168A variants conferring approximately 2-log reduction in MK-5172 susceptibility emerged in 2 of the 4 evaluable patients with genotype 1a breakthrough. Clinical Trials Registration. NCT01353911.

Susceptibility loci for lung cancer are associated with mRNA levels of nearby genes in the lung

Recent studies identified three genetic loci reproducibly associated with lung cancer in populations of European ancestry, namely 15q25, 5p15 and 6p21. The goals of this study are first to confirm whether these loci are associated with lung cancer in a French Canadian population and second to identify disease-associated single nucleotide polymorphisms (SNPs) influencing messenger RNA (mRNA) expression levels of genes in the lung, that is expression quantitative trait loci (eQTLs). SNPs were genotyped in 420 patients undergoing lung cancer surgery and compared with 3151 controls of European ancestry. Genome-wide gene expression levels in non-tumor lung tissues of the same 420 patients were also measured to identify eQTLs. Significant eQTLs were then followed-up in two replication sets (n = 339 and 363). SNPs found in the three susceptibility loci were associated with lung cancer in the French Canadian population. Strong eQTLs were found on chromosome 15q25 with the expression levels of CHRNA5 (P = 2.23 × 10(-) (22) with rs12907966). The CHRNA5-rs12907966 eQTL was convincingly validated in the two replication sets (P = 3.46 × 10(-) (16) and 2.01 × 10(-) (15)). On 6p21, a trend was observed for rs3131379 to be associated with the expression of APOM (P = 3.58 × 10(-) (4)) and validated in the replication sets (P = 1.11 × 10(-) (8) and 6.84 × 10(-) (4)). On 5p15, no significant eQTLs were found. This study confirmed that chromosomes 15q25, 5p15 and 6p21 harbored susceptibility loci for lung cancer in French Canadians. Most importantly, this study suggests that the risk alleles at 15q25 and 6p21 may mediate their effect by regulating the mRNA expression levels of CHRNA5 and APOM in the lung.

Whole exome re-sequencing implicates CCDC38 and cilia structure and function in resistance to smoking related airflow obstruction

Chronic obstructive pulmonary disease (COPD) is a leading cause of global morbidity and mortality and, whilst smoking remains the single most important risk factor, COPD risk is heritable. Of 26 independent genomic regions showing association with lung function in genome-wide association studies, eleven have been reported to show association with airflow obstruction. Although the main risk factor for COPD is smoking, some individuals are observed to have a high forced expired volume in 1 second (FEV1) despite many years of heavy smoking. We hypothesised that these "resistant smokers" may harbour variants which protect against lung function decline caused by smoking and provide insight into the genetic determinants of lung health. We undertook whole exome re-sequencing of 100 heavy smokers who had healthy lung function given their age, sex, height and smoking history and applied three complementary approaches to explore the genetic architecture of smoking resistance. Firstly, we identified novel functional variants in the "resistant smokers" and looked for enrichment of these novel variants within biological pathways. Secondly, we undertook association testing of all exonic variants individually with two independent control sets. Thirdly, we undertook gene-based association testing of all exonic variants. Our strongest signal of association with smoking resistance for a non-synonymous SNP was for rs10859974 (P = 2.34 × 10(-4)) in CCDC38, a gene which has previously been reported to show association with FEV1/FVC, and we demonstrate moderate expression of CCDC38 in bronchial epithelial cells. We identified an enrichment of novel putatively functional variants in genes related to cilia structure and function in resistant smokers. Ciliary function abnormalities are known to be associated with both smoking and reduced mucociliary clearance in patients with COPD. We suggest that genetic influences on the development or function of cilia in the bronchial epithelium may affect growth of cilia or the extent of damage caused by tobacco smoke.

Refining susceptibility loci of chronic obstructive pulmonary disease with lung eqtls

Chronic obstructive pulmonary disease (COPD) is the fourth leading cause of mortality worldwide. Recent genome-wide association studies (GWAS) have identified robust susceptibility loci associated with COPD. However, the mechanisms mediating the risk conferred by these loci remain to be found. The goal of this study was to identify causal genes/variants within susceptibility loci associated with COPD. In the discovery cohort, genome-wide gene expression profiles of 500 non-tumor lung specimens were obtained from patients undergoing lung surgery. Blood-DNA from the same patients were genotyped for 1,2 million SNPs. Following genotyping and gene expression quality control filters, 409 samples were analyzed. Lung expression quantitative trait loci (eQTLs) were identified and overlaid onto three COPD susceptibility loci derived from GWAS; 4q31 (HHIP), 4q22 (FAM13A), and 19q13 (RAB4B, EGLN2, MIA, CYP2A6). Significant eQTLs were replicated in two independent datasets (n = 363 and 339). SNPs previously associated with COPD and lung function on 4q31 (rs1828591, rs13118928) were associated with the mRNA expression of HHIP. An association between mRNA expression level of FAM13A and SNP rs2045517 was detected at 4q22, but did not reach statistical significance. At 19q13, significant eQTLs were detected with EGLN2. In summary, this study supports HHIP, FAM13A, and EGLN2 as the most likely causal COPD genes on 4q31, 4q22, and 19q13, respectively. Strong lung eQTL SNPs identified in this study will need to be tested for association with COPD in case-control studies. Further functional studies will also be needed to understand the role of genes regulated by disease-related variants in COPD.

Emergence of resistance-associated variants after failed triple therapy with vaniprevir in treatment-experienced non-cirrhotic patients with hepatitis C-genotype 1 infection: a population and clonal analysis

BACKGROUND

Vaniprevir with P/R improved SVR rates over P/R alone in treatment-experienced patients with chronic HCV-genotype 1 infection, but treatment failure presents therapeutic challenges. We identified RAVs from non-cirrhotic patients failing to achieve SVR on vaniprevir-containing regimens from a dose/duration-ranging trial of triple-combination therapy.

METHODS

Using population analysis, resistance sequencing was performed on all baseline samples and on samples at virologic failure in the vaniprevir arms. Longitudinal clonal analyses were performed on viral isolates from six vaniprevir recipients experiencing breakthrough viremia.

RESULTS

Baseline RAVs were detected in two patients subsequently experiencing virologic failure. At virologic failure, the majority of RAVs had substitutions at R155, A156, or D168. Clonal analyses identified novel double/triple variants emerging with continuing vaniprevir dosing.

CONCLUSIONS

RAVs were predominantly observed at R155, A156, and/or D168 during virologic failure on vaniprevir/P/R. Double/triple RAVs were identified in patients remaining viremic on triple therapy, suggesting evolution of resistance under selective pressure.

Systems biology approach for new target and biomarker identification

The pharmaceutical industry is spending increasingly large amounts of money on the discovery and development of novel medicines, but this investment is not adequately paying off in an increased rate of newly approved drugs by the FDA. The post-genomic era has provided a wealth of novel approaches for generating large, high-dimensional genetic and transcriptomic data sets from large cohorts of preclinical species as well as normal and diseased individuals. This systems biology approach to understanding disease-related biology is revolutionizing our understanding of the cellular pathways and gene networks underlying the onset of disease, and the mechanisms of pharmacological treatments that ameliorate disease phenotypes. In this article, we review a number of approaches being used by pharmaceutical and biotechnology companies, e.g., high-throughput DNA genotyping, sequencing, and genome-wide gene expression profiling, to enable drug discovery and development through the identification of new drug targets and biomarkers of disease progression, drug pharmacodynamics, and predictive markers for selecting the patients most likely to respond to therapy.

Identification and validation of genes affecting aortic lesions in mice

Atherosclerosis represents the most significant risk factor for coronary artery disease (CAD), the leading cause of death in developed countries. To better understand the pathogenesis of atherosclerosis, we applied a likeli-hood-based model selection method to infer gene-disease causality relationships for the aortic lesion trait in a segregating mouse population demonstrating a spectrum of susceptibility to developing atherosclerotic lesions. We identified 292 genes that tested causal for aortic lesions from liver and adipose tissues of these mice, and we experimentally validated one of these candidate causal genes, complement component 3a receptor 1 (C3ar1), using a knockout mouse model. We also found that genes identified by this method overlapped with genes progressively regulated in the aortic arches of 2 mouse models of atherosclerosis during atherosclerotic lesion development. By comparing our gene set with findings from public human genome-wide association studies (GWAS) of CAD and related traits, we found that 5 genes identified by our study overlapped with published studies in humans in which they were identified as risk factors for multiple atherosclerosis-related pathologies, including myocardial infarction, serum uric acid levels, mean platelet volume, aortic root size, and heart failure. Candidate causal genes were also found to be enriched with CAD risk polymorphisms identified by the Wellcome Trust Case Control Consortium (WTCCC). Our findings therefore validate the ability of causality testing procedures to provide insights into the mechanisms underlying atherosclerosis development.

Targeting of a CD8 T cell env epitope presented by HLA-B*5802 is associated with markers of HIV disease progression and lack of selection pressure

In HIV-infected persons, certain HLA class I alleles are associated with effective control of viremia, while others are associated with rapid disease progression. Among the most divergent clinical outcomes are the relatively good prognosis in HLA-B*5801 expressing persons and poor prognosis with HLA-B*5802. These two alleles differ by only three amino acids in regions involved in HLA-peptide recognition. This study evaluated a cohort of over 1000 persons with chronic HIV clade C virus infection to determine whether clinical outcome differences associated with B*5801 (n = 93) and B*5802 ( n = 259) expression are associated with differences in HIV-1-specific CD8 (+) T cell responses. The overall breadth and magnitude of HIV-1-specific CD8(+) T cell responses were lower in persons expressing B*5802, and epitope presentation by B*5802 contributed significantly less to the overall response as compared to B*5801-restricted CD8 (+) T cells. Moreover, viral load in B*5802-positive persons was higher and CD4 cell counts lower when this allele contributed to the overall CD8 (+) T cell response, which was detected exclusively through a single epitope in Env. In addition, persons heterozygous for B*5802 compared to persons homozygous for other HLA-B alleles had significantly higher viral loads. Viral sequencing revealed strong selection pressure mediated through B*5801-restricted responses but not through B*5802. These data indicate that minor differences in HLA sequence can have a major impact on epitope recognition, and that selective targeting of Env through HLA-B*5802 is at least ineffectual if not actively adverse in the containment of viremia. These results provide experimental evidence that not all epitope-specific responses contribute to immune containment, a better understanding of which is essential to shed light on mechanisms involved in HIV disease progression.

CD8+ T-cell responses to different HIV proteins have discordant associations with viral load

Selection of T-cell vaccine antigens for chronic persistent viral infections has been largely empirical. To define the relationship, at the population level, between the specificity of the cellular immune response and viral control for a relevant human pathogen, we performed a comprehensive analysis of the 160 dominant CD8(+) T-cell responses in 578 untreated HIV-infected individuals from KwaZulu-Natal, South Africa. Of the HIV proteins targeted, only Gag-specific responses were associated with lowering viremia. Env-specific and Accessory/Regulatory protein-specific responses were associated with higher viremia. Increasing breadth of Gag-specific responses was associated with decreasing viremia and increasing Env breadth with increasing viremia. Association of the specific CD8(+) T-cell response with low viremia was independent of HLA type and unrelated to epitope sequence conservation. These population-based data, suggesting the existence of both effective immune responses and responses lacking demonstrable biological impact in chronic HIV infection, are of relevance to HIV vaccine design and evaluation.

Molecular studies reveal frequent misidentification of Aspergillus fumigatus by morphotyping

Aspergillus fumigatus has been understood to be the most common cause of invasive aspergillosis (IA) in all epidemiological surveys. However, recent studies have uncovered a large degree of genetic heterogeneity between isolates morphologically identified as A. fumigatus, leading to the description of a new species, Aspergillus lentulus. Here, we examined the genetic diversity of clinical isolates identified as A. fumigatus using restriction enzyme polymorphism analysis and sequence-based identification. Analysis of 50 clinical isolates from geographically diverse locations recorded the presence of at least three distinct species: A. lentulus, Aspergillus udagawae, and A. fumigatus. In vitro, A. lentulus isolates demonstrated decreased susceptibility to antifungal drugs currently used for IA, including amphotericin B, voriconazole, and caspofungin; A. udagawae isolates demonstrated decreased in vitro susceptibility to amphotericin B. Results of the present study demonstrate that current phenotypic methods to identify fungi do not differentiate between genetically distinct species in the A. fumigatus group. Differential antifungal susceptibilities of these species may account for some of the reported poor outcomes of therapy in clinical studies.

Aspergillus lentulus sp. nov., a new sibling species of A. fumigatus

In a prior study, we identified seven clinical isolates of an Aspergillus sp. that were slow to sporulate in multiple media and demonstrated decreased in vitro susceptibilities to multiple antifungals, including amphotericin B, itraconazole, voriconazole, and caspofungin. These isolates were initially considered to be variants of Aspergillus fumigatus because of differences in mitochondrial cytochrome b sequences and unique randomly amplified polymorphic DNA PCR patterns (S. A. Balajee, M. Weaver, A. Imhof, J. Gribskov, and K. A. Marr, Antimicrob. Agents Chemother. 48: 1197-1203, 2004). The present study was performed to clarify the taxonomic status of these organisms by phylogenetic analyses based on multilocus sequence typing of five genes (the beta-tubulin gene, the rodlet A gene, the salt-responsive gene, the mitochondrial cytochrome b gene, and the internal transcribed spacer regions). Results revealed that four of the seven variant isolates clustered together in a clade very distant from A. fumigatus and distinct from other members of the A. fumigatus group. This new clade, consisting of four members, was monophyletic with strong bootstrap support when the protein-encoding regions were analyzed, indicating a new species status under the phylogenetic species concept. Phenotype studies revealed that the variant isolate has smaller conidial heads with diminutive vesicles compared to A. fumigatus and is not able to survive at 48 degrees C. Our findings suggest the presence of a previously unrecognized, potentially drug-resistant Aspergillus species that we designate A. lentulus.

Compartmentalization of human immunodeficiency virus type 1 between blood monocytes and CD4+ T cells during infection

Distinct sequences of human immunodeficiency virus type 1 (HIV-1) have been found between different tissue compartments or subcompartments within a given tissue. Whether such compartmentalization of HIV-1 occurs between different cell populations is still unknown. Here we address this issue by comparing HIV-1 sequences in the second constant region through the fifth hypervariable region (C2 to V5) of the surface envelope glycoprotein (Env) between viruses in purified blood CD14(+) monocytes and CD4(+) T cells obtained longitudinally from five infected patients over a time period ranging from 117 to 3,409 days postseroconversion. Viral populations in both cell types at early infection time points appeared relatively homogeneous. However, later in infections, all five patients showed heterogeneous populations in both CD14(+) monocytes and CD4(+) T cells. Three of the five patients had CD14(+) monocyte populations with significantly more genetic diversity than the CD4(+) T-cell population, while the other two patients had more genetic diversity in CD4(+) T cells. The cellular compartmentalization of HIV-1 between CD14(+) monocytes and CD4(+) T cells was not seen early during infections but was evident at the later time points for all five patients, indicating an association of viral compartmentalization with the time course of HIV-1 infection. The majority of HIV-1 V3 sequences indicated a macrophage-tropic phenotype, while a V3 sequence-predicted T-cell tropic virus was found in the CD4(+) T cells and CD14(+) monocytes of two patients. These findings suggest that HIV-1 in CD14(+) monocytes could disseminate and evolve independently from that in CD4(+) T cells over the course of HIV-1 infection, which may have implications on the development of new therapeutic strategies.

Evidence for human immunodeficiency virus type 1 replication in vivo in CD14(+) monocytes and its potential role as a source of virus in patients on highly active antiretroviral therapy

In vitro studies show that human immunodeficiency virus type 1 (HIV-1) does not replicate in freshly isolated monocytes unless monocytes differentiate to monocyte-derived macrophages. Similarly, HIV-1 may replicate in macrophages in vivo, whereas it is unclear whether blood monocytes are permissive to productive infection with HIV-1. We investigated HIV-1 replication in CD14(+) monocytes and resting and activated CD4(+) T cells by measuring the levels of cell-associated viral DNA and mRNA and the genetic evolution of HIV-1 in seven acutely infected patients whose plasma viremia had been <100 copies/ml for 803 to 1,544 days during highly active antiretroviral therapy (HAART). HIV-1 DNA was detected in CD14(+) monocytes as well as in activated and resting CD4(+) T cells throughout the course of study. While significant variation in the decay slopes of HIV-1 DNA was seen among individual patients, viral decay in CD14(+) monocytes was on average slower than that in activated and resting CD4(+) T cells. Measurements of HIV-1 sequence evolution and the concentrations of unspliced and multiply spliced mRNA provided evidence of ongoing HIV-1 replication, more pronounced in CD14(+) monocytes than in resting CD4(+) T cells. Phylogenetic analyses of HIV-1 sequences indicated that after prolonged HAART, viral populations related or identical to those found only in CD14(+) monocytes were seen in plasma from three of the seven patients. In the other four patients, HIV-1 sequences in plasma and the three cell populations were identical. CD14(+) monocytes appear to be one of the potential in vivo sources of HIV-1 in patients receiving HAART.