Soumyaroop Bhattacharya

Integration of genomic and genetic approaches implicates IREB2 as a COPD susceptibility gene.

Chronic obstructive pulmonary disease (COPD) is the fourth leading cause of death worldwide and is influenced by both genetic determinants and smoking. We identified genomic regions from 56 lung-tissue gene-expression microarrays and used them to select 889 SNPs to be tested for association with COPD. We genotyped SNPs in 389 severe COPD cases from the National Emphysema Treatment Trial and 424 cigarette-smoking controls from the Normative Aging Study. A total of 71 autosomal SNPs demonstrated at least nominal significance with COPD susceptibility (p = 3.4 x 10(-6) to 0.05). These 71 SNPs were evaluated in a family-based study of 127 probands with severe, early-onset COPD and 822 of their family members in the Boston Early-Onset COPD Study. We combined p values from the case-control and family-based analyses, setting p = 5.60 x 10(-5) as a conservative threshold for significance. Three SNPs in the iron regulatory protein 2 (IREB2) gene met this stringent threshold for significance, and four other IREB2 SNPs demonstrated combined p < 0.02. We demonstrated replication of association for these seven IREB2 SNPs (all p values < or = 0.02) in a family-based study of 3117 subjects from the International COPD Genetics Network; combined p values across all cohorts for the main phenotype of interest ranged from 1.6 x 10(-7) to 6.4 x 10(-4). IREB2 protein and mRNA were increased in lung-tissue samples from COPD subjects in comparison to controls. In summary, gene-expression and genetic-association results have implicated IREB2 as a COPD susceptibility gene.

Molecular Biomarkers for Quantitative and Discrete COPD Phenotypes

Chronic obstructive pulmonary disease (COPD) is an inflammatory lung disorder with complex pathological features and largely unknown etiology. The identification of biomarkers for this disease could aid the development of methods to facilitate earlier diagnosis, the classification of disease subtypes, and provide a means to define therapeutic response. To identify gene expression biomarkers, we completed expression profiling of RNA derived from the lung tissue of 56 subjects with varying degrees of airflow obstruction using the Affymetrix U133 Plus 2.0 array. We applied multiple, independent analytical methods to define biomarkers for either discrete or quantitative disease phenotypes. Analysis of differential expression between cases (n = 15) and controls (n = 18) identified a set of 65 discrete biomarkers. Correlation of gene expression with quantitative measures of airflow obstruction (FEV(1)%predicted or FEV(1)/FVC) identified a set of 220 biomarkers. Biomarker genes were enriched in functions related to DNA binding and regulation of transcription. We used this group of biomarkers to predict disease in an unrelated data set, generated from patients with severe emphysema, with 97% accuracy. Our data contribute to the understanding of gene expression changes occurring in the lung tissue of patients with obstructive lung disease and provide additional insight into potential mechanisms involved in the disease process. Furthermore, we present the first gene expression biomarker for COPD validated in an independent data set.

Epithelial cell PPARγ contributes to normal lung maturation

Peroxisome proliferator‐activated receptor (PPAR)‐γ is a member of the nuclear hormone receptor superfamily that can promote cellular differentiation and organ development. PPARγ expression has been reported in a number of pulmonary cell types, including inflammatory, mesenchymal, and epithelial cells. We find that PPARγ is prominently expressed in the airway epithelium in the mouse lung. In an effort to define the physiological role of PPARγ within the lung, we have ablated PPARγ using a novel line of mice capable of specifically targeting the airway epithelium. Airway epithelial cell PPARγ‐targeted mice display enlarged airspaces resulting from insufficient postnatal lung maturation. The increase in airspace size is accompanied by alterations in lung physiology, including increased lung volumes and decreased tissue resistance. Genome‐wide expression profiling reveals a reduction in structural extracellular matrix (ECM) gene expression in conditionally targeted mice, suggesting a disruption in epithelial‐mesenchymal interactions necessary for the establishment of normal lung structure. Expression profiling of airway epithelial cells isolated from conditionally targeted mice indicates PPARγ regulates genes encoding known PPARγ targets, additional lipid metabolism enzymes, and markers of cellular differentiation. These data reveal airway epithelial cell PPARγ is necessary for normal lung structure and function.—Simon, D. M., Arikan, M. C., Srisuma, S., Bhattacharya, S., Tsai, L. W., Ingenito, E. P., Gonzalez, F., Shapiro, S. D., and Mariani, T. J. Epithelial cell PPARγ contributes to normal lung maturation. FASEB J. 20, E710–E720 (2006)

Gene expression profiling of Epstein-Barr virus-positive and -negative monomorphic B-cell posttransplant lymphoproliferative disorders.

Although most posttransplant lymphoproliferative disorders (PTLD) are related to Epstein-Barr virus (EBV) infection, approximately 20% lack detectable EBV (EBV−). It is uncertain whether the latter cases are truly distinct from EBV+ PTLD or possibly relate to another infectious agent. This study used gene expression profiling to further investigate the relationship between EBV+ and EBV− monomorphic B-cell PTLD, and to search for clues to their pathogenesis. Affymetrix HU133A GeneChips were used to compare 4 EBV+ and 4 EBV− cases of monomorphic B-cell PTLD. Hierarchical clustering successfully distinguished the EBV+ and EBV− groups. Relative to EBV− PTLD, 54 transcripts were over-expressed in EBV+ PTLD. The transcripts identified included IRF7 (a known regulator of EBV LMP1 expression), EBI2 (EBV-induced gene 2), and 3 that are interferon induced (MX1, IFITM1, and IFITM3). In addition, the EBV+ group contained 232 transcripts decreased relative to the EBV− group, including changes concordant with those previously reported after EBV infection of cultured B-cell lines. In summary, in a small group of monomorphic B-cell PTLD, EBV+ cases demonstrated a subset of gene expression changes associated with EBV infection of B cells. By contrast, EBV− PTLD lacked viral-associated changes suggesting that they are biologically distinct.

Fibroblast Growth Factor Receptors Control Epithelial–Mesenchymal Interactions Necessary for Alveolar Elastogenesis

RATIONALE The mechanisms contributing to alveolar formation are poorly understood. A better understanding of these processes will improve efforts to ameliorate lung disease of the newborn and promote alveolar repair in the adult. Previous studies have identified impaired alveogenesis in mice bearing compound mutations of fibroblast growth factor (FGF) receptors (FGFRs) 3 and 4, indicating that these receptors cooperatively promote postnatal alveolar formation. OBJECTIVES To determine the molecular and cellular mechanisms of FGF-mediated alveolar formation. METHODS Compound FGFR3/FGFR4-deficient mice were assessed for temporal changes in lung growth, airspace morphometry, and genome-wide expression. Observed gene expression changes were validated using quantitative real-time RT-PCR, tissue biochemistry, histochemistry, and ELISA. Autocrine and paracrine regulatory mechanisms were investigated using isolated lung mesenchymal cells and type II pneumocytes. MEASUREMENTS AND MAIN RESULTS Quantitative analysis of airspace ontogeny confirmed a failure of secondary crest elongation in compound mutant mice. Genome-wide expression profiling identified molecular alterations in these mice involving aberrant expression of numerous extracellular matrix molecules. Biochemical and histochemical analysis confirmed changes in elastic fiber gene expression resulted in temporal increases in elastin deposition with the loss of typical spatial restriction. No abnormalities in elastic fiber gene expression were observed in isolated mesenchymal cells, indicating that abnormal elastogenesis in compound mutant mice is not cell autonomous. Increased expression of paracrine factors, including insulin-like growth factor-1, in freshly-isolated type II pneumocytes indicated that these cells contribute to the observed pathology. CONCLUSIONS Epithelial/mesenchymal signaling mechanisms appear to contribute to FGFR-dependent alveolar elastogenesis and proper airspace formation.

Genome-Wide Transcriptional Profiling Reveals Connective Tissue Mast Cell Accumulation in Bronchopulmonary Dysplasia

RATIONALE Bronchopulmonary dysplasia (BPD) is a major complication of premature birth. Risk factors for BPD are complex and include prenatal infection and O(2) toxicity. BPD pathology is equally complex and characterized by inflammation and dysmorphic airspaces and vasculature. Due to the limited availability of clinical samples, an understanding of the molecular pathogenesis of this disease and its causal mechanisms and associated biomarkers is limited. OBJECTIVES Apply genome-wide expression profiling to define pathways affected in BPD lungs. METHODS Lung tissue was obtained at autopsy from 11 BPD cases and 17 age-matched control subjects without BPD. RNA isolated from these tissue samples was interrogated using microarrays. Standard gene selection and pathway analysis methods were applied to the data set. Abnormal expression patterns were validated by quantitative reverse transcriptase-polymerase chain reaction and immunohistochemistry. MEASUREMENTS AND MAIN RESULTS We identified 159 genes differentially expressed in BPD tissues. Pathway analysis indicated previously appreciated (e.g., DNA damage regulation of cell cycle) as well as novel (e.g., B-cell development) biological functions were affected. Three of the five most highly induced genes were mast cell (MC)-specific markers. We confirmed an increased accumulation of connective tissue MC(TC) (chymase expressing) mast cells in BPD tissues. Increased expression of MC(TC) markers was also demonstrated in an animal model of BPD-like pathology. CONCLUSIONS We present a unique genome-wide expression data set from human BPD lung tissue. Our data provide information on gene expression patterns associated with BPD and facilitated the discovery that MC(TC) accumulation is a prominent feature of this disease. These observations have significant clinical and mechanistic implications.

Identification of FGF7 as a novel susceptibility locus for chronic obstructive pulmonary disease

Rationale Traditional genome-wide association studies (GWASs) of large cohorts of subjects with chronic obstructive pulmonary disease (COPD) have successfully identified novel candidate genes, but several other plausible loci do not meet strict criteria for genome-wide significance after correction for multiple testing. Objectives The authors hypothesise that by applying unbiased weights derived from unique populations we can identify additional COPD susceptibility loci. Methods The authors performed a homozygosity haplotype analysis on a group of subjects with and without COPD to identify regions of conserved homozygosity haplotype (RCHHs). Weights were constructed based on the frequency of these RCHHs in case versus controls, and used to adjust the p values from a large collaborative GWAS of COPD. Results The authors identified 2318 RCHHs, of which 576 were significantly (p<0.05) over-represented in cases. After applying the weights constructed from these regions to a collaborative GWAS of COPD, the authors identified two single nucleotide polymorphisms (SNPs) in a novel gene (fibroblast growth factor-7 (FGF7)) that gained genome-wide significance by the false discovery rate method. In a follow-up analysis, both SNPs (rs12591300 and rs4480740) were significantly associated with COPD in an independent population (combined p values of 7.9E–7 and 2.8E–6, respectively). In another independent population, increased lung tissue FGF7 expression was associated with worse measures of lung function. Conclusion Weights constructed from a homozygosity haplotype analysis of an isolated population successfully identify novel genetic associations from a GWAS on a separate population. This method can be used to identify promising candidate genes that fail to meet strict correction for multiple testing.

Peripheral blood gene expression profiles in COPD subjects

To identify non-invasive gene expression markers for chronic obstructive pulmonary disease (COPD), we performed genome-wide expression profiling of peripheral blood samples from 12 subjects with significant airflow obstruction and an equal number of non-obstructed controls. RNA was isolated from Peripheral Blood Mononuclear Cells (PBMCs) and gene expression was assessed using Affymetrix U133 Plus 2.0 arrays.Tests for gene expression changes that discriminate between COPD cases (FEV1< 70% predicted, FEV1/FVC < 0.7) and controls (FEV1> 80% predicted, FEV1/FVC > 0.7) were performed using Significance Analysis of Microarrays (SAM) and Bayesian Analysis of Differential Gene Expression (BADGE). Using either test at high stringency (SAM median FDR = 0 or BADGE p < 0.01) we identified differential expression for 45 known genes. Correlation of gene expression with lung function measurements (FEV1 & FEV1/FVC), using both Pearson and Spearman correlation coefficients (p < 0.05), identified a set of 86 genes. A total of 16 markers showed evidence of significant correlation (p < 0.05) with quantitative traits and differential expression between cases and controls. We further compared our peripheral gene expression markers with those we previously identified from lung tissue of the same cohort. Two genes, RP9and NAPE-PLD, were identified as decreased in COPD cases compared to controls in both lung tissue and blood. These results contribute to our understanding of gene expression changes in the peripheral blood of patients with COPD and may provide insight into potential mechanisms involved in the disease.

Progression-Associated Genes in Astrocytoma Identified by Novel Microarray Gene Expression Data Reanalysis

Astrocytoma is graded as pilocytic (WHO grade I), diffuse (WHO grade II), anaplastic (WHO grade III), and glioblastoma multiforme (WHO grade IV). The progression from low- to high-grade astrocytoma is associated with distinct molecular changes that vary with patient age, yet the prognosis of high-grade tumors in children and adults is equally dismal. Whether specific gene expression changes are consistently associated with all high-grade astrocytomas, independent of patient age, is not known. To address this question, we reanalyzed the microarray datasets comprising astrocytomas from children and adults, respectively. We identified nine genes consistently dysregulated in high-grade tumors, using four novel tests for identifying differentially expressed genes. Four genes encoding ribosomal proteins (RPS2, RPS8, RPS18, RPL37A) were upregulated, and five genes (APOD, SORL1, SPOCK2, PRSS11, ID3) were downregulated in high-grade by all tests. Expression results were validated using a third astrocytoma dataset. APOD, the most differentially expressed gene, has been shown to inhibit tumor cell and vascular smooth muscle cell proliferation. This suggests that dysregulation of APOD may be critical for malignant astrocytoma formation, and thus a possible novel universal target for therapeutic intervention. Further investigation is needed to evaluate the role of APOD, as well as the other genes identified, in malignant astrocytoma development.

Alterations in Gene Expression and DNA Methylation during Murine and Human Lung Alveolar Septation

DNA methylation, a major epigenetic mechanism, may regulate coordinated expression of multiple genes at specific time points during alveolar septation in lung development. The objective of this study was to identify genes regulated by methylation during normal septation in mice and during disordered septation in bronchopulmonary dysplasia. In mice, newborn lungs (preseptation) and adult lungs (postseptation) were evaluated by microarray analysis of gene expression and immunoprecipitation of methylated DNA followed by sequencing (MeDIP-Seq). In humans, microarray gene expression data were integrated with genome-wide DNA methylation data from bronchopulmonary dysplasia versus preterm and term lung. Genes with reciprocal changes in expression and methylation, suggesting regulation by DNA methylation, were identified. In mice, 95 genes with inverse correlation between expression and methylation during normal septation were identified. In addition to genes known to be important in lung development (Wnt signaling, Angpt2, Sox9, etc.) and its extracellular matrix (Tnc, Eln, etc.), genes involved with immune and antioxidant defense (Stat4, Sod3, Prdx6, etc.) were also observed. In humans, 23 genes were differentially methylated with reciprocal changes in expression in bronchopulmonary dysplasia compared with preterm or term lung. Genes of interest included those involved with detoxifying enzymes (Gstm3) and transforming growth factor-β signaling (bone morphogenetic protein 7 [Bmp7]). In terms of overlap, 20 genes and three pathways methylated during mouse lung development also demonstrated changes in methylation between preterm and term human lung. Changes in methylation correspond to altered expression of a number of genes associated with lung development, suggesting that DNA methylation of these genes may regulate normal and abnormal alveolar septation.