Gregory, Andrew

Abstract Rationale Although many studies have examined gene expression in lung tissue, the gene regulatory processes underlying emphysema are still not well understood. Finding efficient nonimaging 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 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 chronic obstructive pulmonary disease–associated differential expression of the emphysema pathways found in the bulk analyses. Measurements and Main Results In the bulk lung RNA-sequencing 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.

Chronic obstructive pulmonary disease (COPD) is a disease with marked metabolic disturbance. Previous studies have shown the association between single metabolites and lung function for COPD, but whether a combination of metabolites could predict phenotype is unknown. We developed metabolomic severity scores using plasma metabolomics from the Metabolon platform from two US cohorts of ever-smokers: the Subpopulations and Intermediate Outcome Measures in COPD Study (SPIROMICS) (n = 648; training/testing cohort; 72% non-Hispanic, white; average age 63 years) and the COPDGene Study (n = 1120; validation cohort; 92% non-Hispanic, white; average age 67 years). Separate adaptive LASSO (adaLASSO) models were used to model forced expiratory volume at one second (FEV1) and MESA-adjusted lung density using 762 metabolites common between studies. Metabolite coefficients selected by the adaLASSO procedure were used to create a metabolomic severity score (metSS) for each outcome. A total of 132 metabolites were selected to create a metSS for FEV1. The metSS-only models explained 64.8% and 31.7% of the variability in FEV1 in the training and validation cohorts, respectively. For MESA-adjusted lung density, 129 metabolites were selected, and metSS-only models explained 59.0% of the variability in the training cohort and 17.4% in the validation cohort. Regression models including both clinical covariates and the metSS explained more variability than either the clinical covariate or metSS-only models (53.4% vs. 46.4% and 31.6%) in the validation dataset. The metabolomic pathways for arginine biosynthesis; aminoacyl-tRNA biosynthesis; and glycine, serine, and threonine pathway were enriched by adaLASSO metabolites for FEV1. This is the first demonstration of a respiratory metabolomic severity score, which shows how a metSS can add explanation of variance to clinical predictors of FEV1 and MESA-adjusted lung density. The advantage of a comprehensive metSS is that it explains more disease than individual metabolites and can account for substantial collinearity among classes of metabolites. Future studies should be performed to determine whether metSSs are similar in younger, and more racially and ethnically diverse populations as well as whether a metabolomic severity score can predict disease development in individuals who do not yet have COPD.

Chlorotoxin (CTX) is a 36–amino acid peptide with eight Cys residues that forms four disulfide bonds. It has high affinity for the glioma-specific chloride channel and matrix metalloprotease-2. Structural and binding properties of CTX analogs with various Cys residue substitutions with l-α-aminobutyric acid (Abu) have been previously reported. Using 4.2 µs molecular dynamics, we compared the conformational and essential space sampling of CTX and analogs with selective substitution of the Cys residues and associated disulfide bonds with either Abu or Ser. The native and substituted peptides maintained a high degree of α-helix propensity from residues 8 through 21, with the exception of substitution of the Cys5–Cys28 residues with Ser and the Cys16–Cys33 residues with Abu. In agreement with previous circular dichroism spectropolarimetry results, the C-terminal β-sheet content varied less from residues 25 through 29 and 32 through 36 and was well conserved in most analogs. The Cys16–Cys33 and Cys20–Cys35 disulfide-bonded residues appear to be required to maintain the αβ motif of CTX. Selective substitution with the hydrophilic Ser, may mitigate the destabilizing effect of Cys16–Cys33 substitution through the formation of an inter residue H-bond from Ser16:OγH to Ser33:OγH bridged by a water molecule. All peptides shared considerable sampled conformational space, which explains the retained receptor binding of the non-native analogs.