Publications

2010

Liu F, Mih J, Shea B, Kho A, Sharif A, Tager A, Tschumperlin D. Feedback amplification of fibrosis through matrix stiffening and COX-2 suppression. J Cell Biol. 2010;190(4):693–706. doi:10.1083/jcb.201004082
Tissue stiffening is a hallmark of fibrotic disorders but has traditionally been regarded as an outcome of fibrosis, not a contributing factor to pathogenesis. In this study, we show that fibrosis induced by bleomycin injury in the murine lung locally increases median tissue stiffness sixfold relative to normal lung parenchyma. Across this pathophysiological stiffness range, cultured lung fibroblasts transition from a surprisingly quiescent state to progressive increases in proliferation and matrix synthesis, accompanied by coordinated decreases in matrix proteolytic gene expression. Increasing matrix stiffness strongly suppresses fibroblast expression of COX-2 (cyclooxygenase-2) and synthesis of prostaglandin E(2) (PGE(2)), an autocrine inhibitor of fibrogenesis. Exogenous PGE(2) or an agonist of the prostanoid EP2 receptor completely counteracts the proliferative and matrix synthetic effects caused by increased stiffness. Together, these results demonstrate a dominant role for normal tissue compliance, acting in part through autocrine PGE(2), in maintaining fibroblast quiescence and reveal a feedback relationship between matrix stiffening, COX-2 suppression, and fibroblast activation that promotes and amplifies progressive fibrosis.
Kho A, Bhattacharya S, Tantisira K, Carey V, Gaedigk R, Leeder S, Kohane I, Weiss S, Mariani T. Transcriptomic analysis of human lung development. Am J Respir Crit Care Med. 2010;181(1):54–63. doi:10.1164/rccm.200907-1063OC
RATIONALE: Current understanding of the molecular regulation of lung development is limited and derives mostly from animal studies. OBJECTIVES: To define global patterns of gene expression during human lung development. METHODS: Genome-wide expression profiling was used to measure the developing lung transcriptome in RNA samples derived from 38 normal human lung tissues at 53 to 154 days post conception. Principal component analysis was used to characterize global expression variation and to identify genes and bioontologic attributes contributing to these variations. Individual gene expression patterns were verified by quantitative reverse transcriptase-polymerase chain reaction analysis. MEASUREMENTS AND MAIN RESULTS: Gene expression analysis identified attributes not previously associated with lung development, such as chemokine-immunologic processes. Lung characteristics attributes (e.g., surfactant function) were observed at an earlier-than-anticipated age. We defined a 3,223 gene developing lung characteristic subtranscriptome capable of describing a majority of the process. In gene expression space, the samples formed a time-contiguous trajectory with transition points correlating with histological stages and suggesting the existence of novel molecular substages. Induction of surfactant gene expression characterized a pseudoglandular "molecular phase" transition. Individual gene expression patterns were independently validated. We predicted the age of independent human lung transcriptome profiles with a median absolute error of 5 days, supporting the validity of the data and modeling approach. CONCLUSIONS: This study extends our knowledge of key gene expression patterns and bioontologic attributes underlying early human lung developmental processes. The data also suggest the existence of molecular phases of lung development.

2009

Du R, Tantisira K, Carey V, Bhattacharya S, Metje S, Kho A, Klanderman B, Gaedigk R, Lazarus R, Mariani T, et al. Platform dependence of inference on gene-wise and gene-set involvement in human lung development. BMC Bioinformatics. 2009;10:189. doi:10.1186/1471-2105-10-189
BACKGROUND: With the recent development of microarray technologies, the comparability of gene expression data obtained from different platforms poses an important problem. We evaluated two widely used platforms, Affymetrix U133 Plus 2.0 and the Illumina HumanRef-8 v2 Expression Bead Chips, for comparability in a biological system in which changes may be subtle, namely fetal lung tissue as a function of gestational age. RESULTS: We performed the comparison via sequence-based probe matching between the two platforms. "Significance grouping" was defined as a measure of comparability. Using both expression correlation and significance grouping as measures of comparability, we demonstrated that despite overall cross-platform differences at the single gene level, increased correlation between the two platforms was found in genes with higher expression level, higher probe overlap, and lower p-value. We also demonstrated that biological function as determined via KEGG pathways or GO categories is more consistent across platforms than single gene analysis. CONCLUSION: We conclude that while the comparability of the platforms at the single gene level may be increased by increasing sample size, they are highly comparable ontologically even for subtle differences in a relatively small sample size. Biologically relevant inference should therefore be reproducible across laboratories using different platforms.
Arashiro P, Eisenberg I, Kho A, Cerqueira A, Canovas M, Silva H, Pavanello R, Verjovski-Almeida S, Kunkel L, Zatz M. Transcriptional regulation differs in affected facioscapulohumeral muscular dystrophy patients compared to asymptomatic related carriers. Proc Natl Acad Sci U S A. 2009;106(15):6220–5. doi:10.1073/pnas.0901573106
Facioscapulohumeral muscular dystrophy (FSHD) is a progressive muscle disorder that has been associated with a contraction of 3.3-kb repeats on chromosome 4q35. FSHD is characterized by a wide clinical inter- and intrafamilial variability, ranging from wheelchair-bound patients to asymptomatic carriers. Our study is unique in comparing the gene expression profiles from related affected, asymptomatic carrier, and control individuals. Our results suggest that the expression of genes on chromosome 4q is altered in affected and asymptomatic individuals. Remarkably, the changes seen in asymptomatic samples are largely in products of genes encoding several chemokines, whereas the changes seen in affected samples are largely in genes governing the synthesis of GPI-linked proteins and histone acetylation. Besides this, the affected patient and related asymptomatic carrier share the 4qA161 haplotype. Thus, these polymorphisms by themselves do not explain the pathogenicity of the contracted allele. Interestingly, our results also suggest that the miRNAs might mediate the regulatory network in FSHD. Together, our results support the previous evidence that FSHD may be caused by transcriptional dysregulation of multiple genes, in cis and in trans, and suggest some factors potentially important for FSHD pathogenesis. The study of the gene expression profiles from asymptomatic carriers and related affected patients is a unique approach to try to enhance our understanding of the missing link between the contraction in D4Z4 repeats and muscle disease, while minimizing the effects of differences resulting from genetic background.
Kho A, Bhattacharya S, Mecham B, Hong J, Kohane I, Mariani T. Expression profiles of the mouse lung identify a molecular signature of time-to-birth. Am J Respir Cell Mol Biol. 2009;40(1):47–57. doi:10.1165/rcmb.2008-0048OC
A greater understanding of the regulatory processes contributing to lung development could help ameliorate morbidity and mortality in premature infants and identify individuals at risk for congenital and/or chronic lung diseases. Genomics technologies have provided rich gene expression datasets for the developing lung that enable systems biology approaches for identifying large-scale molecular signatures within this complex phenomenon. Here, we applied unsupervised principal component analysis on two developing lung datasets and identified common dominant transcriptomic signatures. Of particular interest, we identify an overlying biological program we term "time-to-birth," which describes the distance in age from the day of birth. We identify groups of genes contributing to the time-to-birth molecular signature. Statistically overrepresented are genes involved in oxygen and gas transport activity, as expected for a transition to air breathing, as well as host defense function. In addition, we identify genes with expression patterns associated with the initiation of alveolar formation. Finally, we present validation of gene expression patterns across the two datasets, and independent validation of select genes by qPCR and immunohistochemistry. These data contribute to our understanding of genetic components contributing to large-scale biological processes and may be useful, particularly in animal models of abnormal lung development, to predict the state of organ development or preparation for birth.

2008

Schüller U, Heine V, Mao J, Kho A, Dillon A, Han Y-G, Huillard E, Sun T, Ligon A, Qian Y, et al. Acquisition of granule neuron precursor identity is a critical determinant of progenitor cell competence to form Shh-induced medulloblastoma. Cancer Cell. 2008;14(2):123–34. doi:10.1016/j.ccr.2008.07.005
Whether the brain tumor medulloblastoma originates from stem cells or restricted progenitor cells is unclear. To investigate this, we activated oncogenic Hedgehog (Hh) signaling in multipotent and lineage-restricted central nervous system (CNS) progenitors. We observed that normal unipotent cerebellar granule neuron precursors (CGNPs) derive from hGFAP(+) and Olig2(+) rhombic lip progenitors. Hh activation in a spectrum of early- and late-stage CNS progenitors generated similar medulloblastomas, but not other brain cancers, indicating that acquisition of CGNP identity is essential for tumorigenesis. We show in human and mouse medulloblastoma that cells expressing the glia-associated markers Gfap and Olig2 are neoplastic and retain features of embryonic-type granule lineage progenitors. Thus, oncogenic Hh signaling promotes medulloblastoma from lineage-restricted granule cell progenitors.

2007

Gobezie R, Kho A, Krastins B, Sarracino D, Thornhill T, Chase M, Millett P, Lee D. High abundance synovial fluid proteome: distinct profiles in health and osteoarthritis. Arthritis Res Ther. 2007;9(2):R36. doi:10.1186/ar2172
The development of increasingly high-throughput and sensitive mass spectroscopy-based proteomic techniques provides new opportunities to examine the physiology and pathophysiology of many biologic fluids and tissues. The purpose of this study was to determine protein expression profiles of high-abundance synovial fluid (SF) proteins in health and in the prevalent joint disease osteoarthritis (OA). A cross-sectional study of 62 patients with early OA (n = 21), patients with late OA (n = 21), and control individuals (n = 20) was conducted. SF proteins were separated by using one-dimensional PAGE, and the in-gel digested proteins were analyzed by electrospray ionization tandem mass spectrometry. A total of 362 spots were examined and 135 high-abundance SF proteins were identified as being expressed across all three study cohorts. A total of 135 SF proteins were identified. Eighteen proteins were found to be significantly differentially expressed between control individuals and OA patients. Two subsets of OA that are not dependent on disease duration were identified using unsupervised analysis of the data. Several novel SF proteins were also identified. Our analyses demonstrate no disease duration-dependent differences in abundant protein composition of SF in OA, and we clearly identified two previously unappreciated yet distinct subsets of protein profiles in this disease cohort. Additionally, our findings reveal novel abundant protein species in healthy SF whose functional contribution to SF physiology was not previously recognized. Finally, our studies identify candidate biomarkers for OA with potential for use as highly sensitive and specific tests for diagnostic purposes or for evaluating therapeutic response.

2006

Schüller U, Kho A, Zhao Q, Ma Q, Rowitch D. Cerebellar 'transcriptome' reveals cell-type and stage-specific expression during postnatal development and tumorigenesis. Mol Cell Neurosci. 2006;33(3):247–59. doi:10.1016/j.mcn.2006.07.010
Disorders of cerebellar development can result in neurological disease and cancer. The identity of transcription factors that may uniquely mark and/or regulate development of single cerebellar cell types, however, is poorly understood. We used a library of approximately 1100 probes for expression of transcription factor (TF)-encoding genes (>70% of the mammalian 'transcriptome') to identify 227 genes with expression in developing neuronal and glial populations and 24 TFs that show cell-type- and stage-specific expression in granule cells, Purkinje cells and interneurons during postnatal cerebellar development. The utility of this panel is exemplified by analysis of medulloblastoma that shows upregulation of markers specific for early granule cell lineage, but not for other neuronal cell types, indicative of a unipotent precursor as well as a block in granule cell differentiation within the tumor. We propose that this atlas of the cerebellar transcriptome and the panel of 24 validated markers will be generally useful in analyses of mutations affecting postnatal cerebellar development and neoplasia.
Gazda H, Kho A, Sanoudou D, Zaucha J, Kohane I, Sieff C, Beggs A. Defective ribosomal protein gene expression alters transcription, translation, apoptosis, and oncogenic pathways in Diamond-Blackfan anemia. Stem Cells. 2006;24(9):2034–44. doi:10.1634/stemcells.2005-0554
Diamond-Blackfan anemia (DBA) is a broad developmental disease characterized by anemia, bone marrow (BM) erythroblastopenia, and an increased incidence of malignancy. Mutations in ribosomal protein gene S19 (RPS19) are found in approximately 25% of DBA patients; however, the role of RPS19 in the pathogenesis of DBA remains unknown. Using global gene expression analysis, we compared highly purified multipotential, erythroid, and myeloid BM progenitors from RPS19 mutated and control individuals. We found several ribosomal protein genes downregulated in all DBA progenitors. Apoptosis genes, such as TNFRSF10B and FAS, transcriptional control genes, including the erythropoietic transcription factor MYB (encoding c-myb), and translational genes were greatly dysregulated, mostly in diseased erythroid cells. Cancer-related genes, including RAS family oncogenes and tumor suppressor genes, were significantly dysregulated in all diseased progenitors. In addition, our results provide evidence that RPS19 mutations lead to codownregulation of multiple ribosomal protein genes, as well as downregulation of genes involved in translation in DBA cells. In conclusion, the altered expression of cancer-related genes suggests a molecular basis for malignancy in DBA. Downregulation of c-myb expression, which causes complete failure of fetal liver erythropoiesis in knockout mice, suggests a link between RPS19 mutations and reduced erythropoiesis in DBA.
Kho A, Kang P, Kohane I, Kunkel L. Transcriptome-scale similarities between mouse and human skeletal muscles with normal and myopathic phenotypes. BMC Musculoskelet Disord. 2006;7:23. doi:10.1186/1471-2474-7-23
BACKGROUND: Mouse and human skeletal muscle transcriptome profiles vary by muscle type, raising the question of which mouse muscle groups have the greatest molecular similarities to human skeletal muscle. METHODS: Orthologous (whole, sub-) transcriptome profiles were compared among four mouse-human transcriptome datasets: (M) six muscle groups obtained from three mouse strains (wildtype, mdx, mdx5cv); (H1) biopsied human quadriceps from controls and Duchenne muscular dystrophy patients; (H2) four different control human muscle types obtained at autopsy; and (H3) 12 different control human tissues (ten non-muscle). RESULTS: Of the six mouse muscles examined, mouse soleus bore the greatest molecular similarities to human skeletal muscles, independent of the latters' anatomic location/muscle type, disease state, age and sampling method (autopsy versus biopsy). Significant similarity to any one mouse muscle group was not observed for non-muscle human tissues (dataset H3), indicating this finding to be muscle specific. CONCLUSION: This observation may be partly explained by the higher type I fiber content of soleus relative to the other mouse muscles sampled.