Genes related to emphysema are enriched for ubiquitination pathways

BACKGROUND

Increased small airway resistance and decreased lung elasticity contribute to the airflow limitation in chronic obstructive pulmonary disease (COPD). The lesion that corresponds to loss of lung elasticity is emphysema; the small airway obstruction is due to inflammatory narrowing and obliteration. Despite their convergence in altered physiology, different mechanisms contribute to these processes. The relationships between gene expression and these specific phenotypes may be more revealing than comparison with lung function.

METHODS

We measured the ratio of alveolar surface area to lung volume (SA/V) in lung tissue from 43 smokers. Two samples from 21 subjects, in which SA/V differed by >49 cm2/mL were profiled to select genes whose expression correlated with SA/V. Significant genes were tested for replication in the 22 remaining subjects.

RESULTS

The level of expression of 181 transcripts was related to SA/V ( p < 0.05). When these genes were tested in the 22 remaining subjects as a replication, thirty of the 181 genes remained significantly associated with SA/V (P < 0.05) and the direction of association was the same in 164/181. Pathway and network analysis revealed enrichment of genes involved in protein ubiquitination, and western blotting showed altered expression of genes involved in protein ubiquitination in obstructed individuals.

CONCLUSION

This study implicates modified protein ubiquitination and degradation as a potentially important pathway in the pathogenesis of emphysema.

Biomarkers of the Hedgehog/Smoothened pathway in healthy volunteers

The Hedgehog (Hh) pathway is involved in oncogenic transformation and tumor maintenance. The primary objective of this study was to select surrogate tissue to measure messenger ribonucleic acid (mRNA) levels of Hh pathway genes for measurement of pharmacodynamic effect. Expression of Hh pathway specific genes was measured by quantitative real time polymerase chain reaction (qRT-PCR) and global gene expression using Affymetrix U133 microarrays. Correlations were made between the expression of specific genes determined by qRT-PCR and normalized microarray data. Gene ontology analysis using microarray data for a broader set of Hh pathway genes was performed to identify additional Hh pathway-related markers in the surrogate tissue. RNA extracted from blood, hair follicle, and skin obtained from healthy subjects was analyzed by qRT-PCR for 31 genes, whereas 8 samples were analyzed for a 7-gene subset. Twelve sample sets, each with ≤500 ng total RNA derived from hair, skin, and blood, were analyzed using Affymetrix U133 microarrays. Transcripts for several Hh pathway genes were undetectable in blood using qRT-PCR. Skin was the most desirable matrix, followed by hair follicle. Whether processed by robust multiarray average or microarray suite 5 (MAS5), expression patterns of individual samples showed co-clustered signals; both normalization methods were equally effective for unsupervised analysis. The MAS5- normalized probe sets appeared better suited for supervised analysis. This work provides the basis for selection of a surrogate tissue and an expression analysis-based approach to evaluate pathway-related genes as markers of pharmacodynamic effect with novel inhibitors of the Hh pathway.

A gene expression signature that classifies human atherosclerotic plaque by relative inflammation status

BACKGROUND

Atherosclerosis is a complex disease requiring improvements in diagnostic techniques and therapeutic treatments. Both improvements will be facilitated by greater exploration of the biology of atherosclerotic plaque. To this end, we carried out large-scale gene expression analysis of human atherosclerotic lesions.

METHODS AND RESULTS

Whole genome expression analysis of 101 plaques from patients with peripheral artery disease identified a robust gene signature (1514 genes) that is dominated by processes related to Toll-like receptor signaling, T-cell activation, cholesterol efflux, oxidative stress response, inflammatory cytokine production, vasoconstriction, and lysosomal activity. Further analysis of gene expression in microdissected carotid plaque samples revealed that this signature is differentially expressed in macrophage-rich and smooth muscle cell-containing regions. A quantitative PCR gene expression panel and inflammatory composite score were developed on the basis of the atherosclerotic plaque gene signature. When applied to serial sections of carotid plaque, the inflammatory composite score was observed to correlate with histological and morphological features related to plaque vulnerability.

CONCLUSIONS

The robust mRNA expression signature identified in the present report is associated with pathological features of vulnerable atherosclerotic plaque and may be useful as a source of biomarkers and targets of novel antiatherosclerotic therapies.

Selection of housekeeping genes for real-time PCR in atopic human bronchial epithelial cells

The stability of housekeeping genes (HKGs) is critical when performing real-time quantitative PCR. To date, the stability of common HKGs has not been systematically compared in human airway epithelial cells (AEC) in normal and atopic subjects. Expression levels of 12 HKGs were measured in AECs from a cohort of 30 healthy atopic nonasthmatic or atopic asthmatic children. Gene expression stability was determined using three different Visual Basic for Applications applets (geNorm, NormFinder and BestKeeper). All 12 HKGs were expressed in AECs. However, the hypoxanthine ribosyltransferase and TATA-binding protein genes were excluded from further analysis due to low expression levels. The cyclophilin A gene was ranked the most stable by all three methods. The expression levels of the beta-actin and glyceraldehyde-3-phosphate dehydrogenase genes were significantly different between the three groups of patients, with atopic asthmatics showing the highest expression levels for both genes. The results suggest that the cyclophilin A gene is the most suitable housekeeping gene analysed for expression studies utilising uncultured bronchial airway epithelial cells from healthy and asthmatic children, and highlight the importance of validating housekeeping genes for each experimental model.

Gene expression profiling in patients with chronic obstructive pulmonary disease and lung cancer

RATIONALE

Chronic obstructive lung disease (COPD) is a common and disabling lung disease for which there are few therapeutic options.

OBJECTIVES

We reasoned that gene expression profiling of COPD lungs could reveal previously unidentified disease pathways.

METHODS

Forty-eight human lung samples were obtained from tissue resected from five nonsmokers, 21 GOLD (Global Initiative for Chronic Obstructive Lung Disease) stage 0, 9 GOLD stage 1, 10 GOLD stage 2, and 3 GOLD stage 3 patients. mRNA from the specimens was profiled using Agilent's Functional ID v2.0 array (Agilent, Santa Clara, CA) containing 23,720 sequences.

MEASUREMENTS AND MAIN RESULTS

The gene expression pattern was influenced by the percentage of the sample made up of parenchyma. Gene expression was related to forced expiratory flow between 25 and 75% of forced expiratory volume (FEF(25-75%) % predicted) revealing a signature gene set of 203 transcripts. Genes involved in extracellular matrix synthesis/degradation and apoptosis were among the up-regulated genes, whereas genes that participate in antiinflammatory responses were down-regulated. Immunohistochemistry confirmed expression of urokinase plasminogen activator (PLAU), urokinase plasminogen activator receptor (PLAUR), and thrombospondin (THBS1) by alveolar macrophages and airway epithelial cells. Genes in this pathway have been shown to be involved in the activation of transforming growth factor (TGF)-beta1 and matrix metalloproteinases and are subject to inhibition by SERPINE2. Interestingly, both TGF-beta1 and SERPINE2 have been identified as candidate genes in COPD genetic linkage and association studies.

CONCLUSIONS

The results provide evidence that genes involved in tissue remodeling and repair are differentially regulated in the lungs of obstructed smokers and suggest that they are potential therapeutic targets. Data deposited in GEO at http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE8500.

Individual-specific variation of gene expression in peripheral blood leukocytes

DNA microarray technology is used to determine gene expression profiles of various cell types, especially abnormal cells, such as cancer. By contrast, relatively little attention has been given to expression profiling of normal tissues. Here we describe studies of gene expression in peripheral blood leukocytes (PBL) from normal individuals sampled multiple times over periods ranging from several weeks up to 6 months. We demonstrate stable patterns of gene expression that differ between individuals. Among the genes whose expression varies by individual is a group of genes responsive to interferon stimulation. Certain individuals ( approximately 10-20% of those tested) showed higher baseline levels and lower inducibility of these genes in response to in vitro interferon stimulation. These studies demonstrate the feasibility of using DNA microarrays to measure the variations in gene expression of PBL from different individuals in response to environmental and genetic factors.

Use of hybridization kinetics for differentiating specific from non-specific binding to oligonucleotide microarrays

Hybridization kinetics were found to be significantly different for specific and non-specific binding of labeled cRNA to surface-bound oligonucleotides on microarrays. We show direct evidence that in a complex sample specific binding takes longer to reach hybridization equilibrium than the non- specific binding. We find that this property can be used to estimate and to correct for the hybridization contributed by non-specific binding. Useful applications are illustrated including the selection of superior oligonucleotides, and the reduction of false positives in exon identification.

Reconstructing potential energy functions from simulated force-induced unbinding processes

One-dimensional stochastic models demonstrate that molecular dynamics simulations of a few nanoseconds can be used to reconstruct the essential features of the binding potential of macromolecules. This can be accomplished by inducing the unbinding with the help of external forces applied to the molecules, and discounting the irreversible work performed on the system by these forces. The fluctuation-dissipation theorem sets a fundamental limit on the precision with which the binding potential can be reconstructed by this method. The uncertainty in the resulting potential is linearly proportional to the irreversible component of work performed on the system during the simulation. These results provide an a priori estimate of the energy barriers observable in molecular dynamics simulations.

Molecular dynamics study of unbinding of the avidin-biotin complex

We report molecular dynamics simulations that induce, over periods of 40-500 ps, the unbinding of biotin from avidin by means of external harmonic forces with force constants close to those of AFM cantilevers. The applied forces are sufficiently large to reduce the overall binding energy enough to yield unbinding within the measurement time. Our study complements earlier work on biotin-streptavidin that employed a much larger harmonic force constant. The simulations reveal a variety of unbinding pathways, the role of key residues contributing to adhesion as well as the spatial range over which avidin binds biotin. In contrast to the previous studies, the calculated rupture forces exceed by far those observed. We demonstrate, in the framework of models expressed in terms of one-dimensional Langevin equations with a schematic binding potential, the associated Smoluchowski equations, and the theory of first passage times, that picosecond to nanosecond simulation of ligand unbinding requires such strong forces that the resulting protein-ligand motion proceeds far from the thermally activated regime of millisecond AFM experiments, and that simulated unbinding cannot be readily extrapolated to the experimentally observed rupture.