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

Calcium Channel Blocker D-\u3cem\u3ecis\u3c/em\u3e-Diltiazem Does Not Slow Retinal Degeneration in the \u3cem\u3ePDE6B\u3c/em\u3e Mutant \u3cem\u3ercd1\u3c/em\u3e Canine Model of Retinitis Pigmentosa

Purpose: D-cis-diltiazem, a calcium channel blocker, has been reported to enhance photoreceptor survival in the rd mouse, a model of retinitis pigmentosa (RP) resulting from mutation of the PDE6B gene. We tested the hypothesis that diltiazem treatment would similarly rescue the canine rcd1 model of RP, which is also caused by a null mutation in the PDE6B gene. Methods: D-cis-diltiazem was delivered orally twice daily to rcd1 affected dogs beginning at 4 weeks of age; untreated age-matched rcd1 dogs served as controls. At 14 weeks, electroretinograms (ERG) were performed on all animals; 14 dogs were euthanized at this age, and 2 dogs at 25 weeks of age. Eyes were enucleated, fixed, and processed for routine histological examination. Results: No significant differences were found in ERG or histopathologic parameters between diltiazem-treated and untreated rcd1 dogs. Neither rcd1 group showed a rod b-wave; ERGs evoked by single white flashes (dark- or light-adapted) and flicker were also identical between groups. Similarly, treated and untreated animals did not differ in the degree of preservation of the photoreceptor layer, confirmed in cell counts within the outer nuclear layer. Conclusions: Treatment of rcd1 affected dogs with D-cis-diltiazem did not modify the photoreceptor disease when results were assessed using either ERG or histopathologic criteria. The positive photoreceptor-rescue effect of calcium channel blockers reported in the rd mouse was thus not generalizable to another species with retinal degeneration due to mutation in the PDE6B gene. Caution needs to be exerted in extrapolation to the comparable human forms of RP

Recognition of HIV-1 Peptides by Host CTL Is Related to HIV-1 Similarity to Human Proteins

Background: While human immunodeficiency virus type 1 (HIV-1)-specific cytotoxic T lymphocytes preferentially target specific regions of the viral proteome, HIV-1 features that contribute to immune recognition are not well understood. One hypothesis is that similarities between HIV and human proteins influence the host immune response, i.e., resemblance between viral and host peptides could preclude reactivity against certain HIV epitopes. Methodology/Principal Findings: We analyzed the extent of similarity between HIV-1 and the human proteome. Proteins from the HIV-1 B consensus sequence from 2001 were dissected into overlapping k-mers, which were then probed against a non-redundant database of the human proteome in order to identify segments of high similarity. We tested the relationship between HIV-1 similarity to host encoded peptides and immune recognition in HIV-infected individuals, and found that HIV immunogenicity could be partially modulated by the sequence similarity to the host proteome. ELISpot responses to peptides spanning the entire viral proteome evaluated in 314 individuals showed a trend indicating an inverse relationship between the similarity to the host proteome and the frequency of recognition. In addition, analysis of responses by a group of 30 HIV-infected individuals against 944 overlapping peptides representing a broad range of individual HIV-1B Nef variants, affirmed that the degree of similarity to the host was significantly lower for peptides with reactive epitopes than for those that were not recognized. Conclusions/Significance: Our results suggest that antigenic motifs that are scarcely represented in human proteins might represent more immunogenic CTL targets not selected against in the host. This observation could provide guidance in the design of more effective HIV immunogens, as sequences devoid of host-like features might afford superior immune reactivity

The Overlap of Lung Tissue Transcriptome of Smoke Exposed Mice with Human Smoking and COPD

Genome-wide mRNA profiling in lung tissue from human and animal models can provide novel insights into the pathogenesis of chronic obstructive pulmonary disease (COPD). While 6 months of smoke exposure are widely used, shorter durations were also reported. The overlap of short term and long-term smoke exposure in mice is currently not well understood, and their representation of the human condition is uncertain. Lung tissue gene expression profiles of six murine smoking experiments (n = 48) were obtained from the Gene Expression Omnibus (GEO) and analyzed to identify the murine smoking signature. The 'human smoking' gene signature containing 386 genes was previously published in the lung eQTL study (n = 1,111). A signature of mild COPD containing 7 genes was also identified in the same study. The lung tissue gene signature of 'severe COPD' (n = 70) contained 4,071 genes and was previously published. We detected 3,723 differentially expressed genes in the 6 month-exposure mice datasets (FDR <0.1). Of those, 184 genes (representing 48% of human smoking) and 1,003 (representing 27% of human COPD) were shared with the human smoking-related genes and the COPD severity-related genes, respectively. There was 4-fold over-representation of human and murine smoking-related genes (P = 6.7 × 10-26) and a 1.4 fold in the severe COPD -related genes (P = 2.3 × 10-12). There was no significant enrichment of the mice and human smoking-related genes in mild COPD signature. These data suggest that murine smoke models are strongly representative of molecular processes of human smoking but less of COPD

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

© 2018 Elsevier Ltd 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

Increased Breadth and Depth of Cytotoxic T Lymphocytes Responses against HIV-1-B Nef by Inclusion of Epitope Variant Sequences

Different vaccine approaches cope with HIV-1 diversity, ranging from centralized1–4 to variability-encompassing5–7 antigens. For all these strategies, a concern remains: how does HIV-1 diversity impact epitope recognition by the immune system? We studied the relationship between HIV-1 diversity and CD8+ T Lymphocytes (CTL) targeting of HIV-1 subtype B Nef using 944 peptides (10-mers overlapping by nine amino acids (AA)) that corresponded to consensus peptides and their most common variants in the HIV-1-B virus population. IFN-γ ELISpot assays were performed using freshly isolated PBMC from 26 HIV-1-infected persons. Three hundred and fifty peptides elicited a response in at least one individual. Individuals targeted a median of 7 discrete regions. Overall, 33% of responses were directed against viral variants but not elicited against consensus-based test peptides. However, there was no significant relationship between the frequency of a 10-mer in the viral population and either its frequency of recognition (Spearman's correlation coefficient ρ = 0.24) or the magnitude of the responses (ρ = 0.16). We found that peptides with a single mutation compared to the consensus were likely to be recognized (especially if the change was conservative) and to elicit responses of similar magnitude as the consensus peptide. Our results indicate that cross-reactivity between rare and frequent variants is likely to play a role in the expansion of CTL responses, and that maximizing antigenic diversity in a vaccine may increase the breadth and depth of CTL responses. However, since there are few obvious preferred pathways to virologic escape, the diversity that may be required to block all potential escape pathways may be too large for a realistic vaccine to accommodate. Furthermore, since peptides were not recognized based on their frequency in the population, it remains unclear by which mechanisms variability-inclusive antigens (i.e., constructs enriched with frequent variants) expand CTL recognition

Nematode endoparasites do not codiversify with their stick insect hosts.

Host-parasite coevolution stems from reciprocal selection on host resistance and parasite infectivity, and can generate some of the strongest selective pressures known in nature. It is widely seen as a major driver of diversification, the most extreme case being parallel speciation in hosts and their associated parasites. Here, we report on endoparasitic nematodes, most likely members of the mermithid family, infecting different Timema stick insect species throughout California. The nematodes develop in the hemolymph of their insect host and kill it upon emergence, completely impeding host reproduction. Given the direct exposure of the endoparasites to the host's immune system in the hemolymph, and the consequences of infection on host fitness, we predicted that divergence among hosts may drive parallel divergence in the endoparasites. Our phylogenetic analyses suggested the presence of two differentiated endoparasite lineages. However, independently of whether the two lineages were considered separately or jointly, we found a complete lack of codivergence between the endoparasitic nematodes and their hosts in spite of extensive genetic variation among hosts and among parasites. Instead, there was strong isolation by distance among the endoparasitic nematodes, indicating that geography plays a more important role than host-related adaptations in driving parasite diversification in this system. The accumulating evidence for lack of codiversification between parasites and their hosts at macroevolutionary scales contrasts with the overwhelming evidence for coevolution within populations, and calls for studies linking micro- versus macroevolutionary dynamics in host-parasite interactions

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 lungs 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 (false discovery rate, FDR &lt; 0.05) with a profile of age-associated proteins identified from a lung tissue proteomics dataset from nine control patients (49-76 years) (FDR &lt; 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 seven ECM proteins with higher expression with age at both gene and protein levels: COL1A1, COL6A1, COL6A2, COL14A1, FBLN2, LTBP4, and LUM. With immunohistochemistry, we demonstrated higher protein levels 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 to developing chronic lung diseases. NEW &amp; NOTEWORTHY We identified seven age-associated extracellular matrix (ECM) proteins, i.e., COL1A1, COL6A1, COL6A2 COL14A1, FBLN2, LTBP4, and LUM with higher transcript and protein levels in human lung tissue with age. Extensive immunohistochemical analysis revealed significant age-associated differences for COL6A2 in whole tissue, parenchyma, airway wall, and vessel, for COL14A1 and LUM in bronchial epithelium, and COL1A1 in parenchyma. Our findings lay a new foundation for the investigation of ECM differences in age-associated chronic lung diseases. </p

HIV-Specific Probabilistic Models of Protein Evolution

Comparative sequence analyses, including such fundamental bioinformatics techniques as similarity searching, sequence alignment and phylogenetic inference, have become a mainstay for researchers studying type 1 Human Immunodeficiency Virus (HIV-1) genome structure and evolution. Implicit in comparative analyses is an underlying model of evolution, and the chosen model can significantly affect the results. In general, evolutionary models describe the probabilities of replacing one amino acid character with another over a period of time. Most widely used evolutionary models for protein sequences have been derived from curated alignments of hundreds of proteins, usually based on mammalian genomes. It is unclear to what extent these empirical models are generalizable to a very different organism, such as HIV-1–the most extensively sequenced organism in existence. We developed a maximum likelihood model fitting procedure to a collection of HIV-1 alignments sampled from different viral genes, and inferred two empirical substitution models, suitable for describing between-and within-host evolution. Our procedure pools the information from multiple sequence alignments, and provided software implementation can be run efficiently in parallel on a computer cluster. We describe how the inferred substitution models can be used to generate scoring matrices suitable for alignment and similarity searches. Our models had a consistently superior fit relative to the best existing models and to parameter-rich data-driven models when benchmarked on independent HIV-1 alignments, demonstrating evolutionary biases in amino-acid substitution that are unique to HIV, and that are not captured by the existing models. The scoring matrices derived from the models showed a marked difference from common amino-acid scoring matrices. The use of an appropriate evolutionary model recovered a known viral transmission history, whereas a poorly chosen model introduced phylogenetic error. We argue that our model derivation procedure is immediately applicable to other organisms with extensive sequence data available, such as Hepatitis C and Influenza A viruses