Hogune Im

Personal Omics Profiling Reveals Dynamic Molecular and Medical Phenotypes

Personalized medicine is expected to benefit from combining genomic information with regular monitoring of physiological states by multiple high-throughput methods. Here, we present an integrative personal omics profile (iPOP), an analysis that combines genomic, transcriptomic, proteomic, metabolomic, and autoantibody profiles from a single individual over a 14 month period. Our iPOP analysis revealed various medical risks, including type 2 diabetes. It also uncovered extensive, dynamic changes in diverse molecular components and biological pathways across healthy and diseased conditions. Extremely high-coverage genomic and transcriptomic data, which provide the basis of our iPOP, revealed extensive heteroallelic changes during healthy and diseased states and an unexpected RNA editing mechanism. This study demonstrates that longitudinal iPOP can be used to interpret healthy and diseased states by connecting genomic information with additional dynamic omics activity.

Measurement of Protein-DNA Interactions In Vivo by Chromatin Immunoprecipitation

Elucidating mechanisms controlling nuclear processes requires an understanding of the nucleoprotein structure of genes at endogenous chromosomal loci. Traditional approaches to measuring protein-DNA interactions in vitro have often failed to provide insights into physiological mechanisms. Given that most transcription factors interact with simple DNA sequence motifs, which are abundantly distributed throughout a genome, it is essential to pinpoint the small subset of sites bound by factors in vivo. Signaling mechanisms induce the assembly and modulation of complex patterns of histone acetylation, methylation, phosphorylation, and ubiquitination, which are crucial determinants of chromatin accessibility. These seemingly complex issues can be directly addressed by a powerful methodology termed the chromatin immunoprecipitation (ChIP) assay. ChIP analysis involves covalently trapping endogenous proteins at chromatin sites, thereby yielding snapshots of protein-DNA interactions and histone modifications within living cells. The chromatin is sonicated to generate small fragments, and an immunoprecipitation is conducted with an antibody against the desired factor or histone modification. Crosslinks are reversed, and polymerase chain reaction (PCR) is used to assess whether DNA sequences are recovered immune-specifically. Chromatin-domain scanning coupled with quantitative analysis is a powerful means of dissecting mechanisms by which signaling pathways target genes within a complex genome.

Highly restricted localization of RNA polymerase II within a locus control region of a tissue-specific chromatin domain

ABSTRACT RNA polymerase II (Pol II) can associate with regulatory elements far from promoters. For the murine β-globin locus, Pol II binds the β-globin locus control region (LCR) far upstream of the β-globin promoters, independent of recruitment to and activation of the βmajor promoter. We describe here an analysis of where Pol II resides within the LCR, how it is recruited to the LCR, and the functional consequences of recruitment. High-resolution analysis of the distribution of Pol II revealed that Pol II binding within the LCR is restricted to the hypersensitive sites. Blocking elongation eliminated the synthesis of genic and extragenic transcripts and eliminated Pol II from the βmajor open reading frame. However, the elongation blockade did not redistribute Pol II at the hypersensitive sites, suggesting that Pol II is recruited to these sites. The distribution of Pol II did not strictly correlate with the distributions of histone acetylation and methylation. As Pol II associates with histone-modifying enzymes, Pol II tracking might be critical for establishing and maintaining broad histone modification patterns. However, blocking elongation did not disrupt the histone modification pattern of the β-globin locus, indicating that Pol II tracking is not required to maintain the pattern.

Whole-exome sequencing identifies tetratricopeptide repeat domain 7A (TTC7A) mutations for combined immunodeficiency with intestinal atresias

BACKGROUND Combined immunodeficiency with multiple intestinal atresias (CID-MIA) is a rare hereditary disease characterized by intestinal obstructions and profound immune defects. OBJECTIVE We sought to determine the underlying genetic causes of CID-MIA by analyzing the exomic sequences of 5 patients and their healthy direct relatives from 5 unrelated families. METHODS We performed whole-exome sequencing on 5 patients with CID-MIA and 10 healthy direct family members belonging to 5 unrelated families with CID-MIA. We also performed targeted Sanger sequencing for the candidate gene tetratricopeptide repeat domain 7A (TTC7A) on 3 additional patients with CID-MIA. RESULTS Through analysis and comparison of the exomic sequence of the subjects from these 5 families, we identified biallelic damaging mutations in the TTC7A gene, for a total of 7 distinct mutations. Targeted TTC7A gene sequencing in 3 additional unrelated patients with CID-MIA revealed biallelic deleterious mutations in 2 of them, as well as an aberrant splice product in the third patient. Staining of normal thymus showed that the TTC7A protein is expressed in thymic epithelial cells, as well as in thymocytes. Moreover, severe lymphoid depletion was observed in the thymus and peripheral lymphoid tissues from 2 patients with CID-MIA. CONCLUSIONS We identified deleterious mutations of the TTC7A gene in 8 unrelated patients with CID-MIA and demonstrated that the TTC7A protein is expressed in the thymus. Our results strongly suggest that TTC7A gene defects cause CID-MIA.

Dynamic regulation of histone H3 methylated at lysine 79 within a tissue-specific chromatin domain.

Post-translational modifications of individual lysine residues of core histones can exert unique functional consequences. For example, methylation of histone H3 at lysine 79 (H3-meK79) has been implicated recently in gene silencing inSaccharomyces cerevisiae. However, the distribution and function of H3-meK79 in mammalian chromatin are not known. We found that H3-meK79 has a variable distribution within the murine β-globin locus in adult erythroid cells, being preferentially enriched at the active βmajor gene. By contrast, acetylated H3 and H4 and H3 methylated at lysine 4 were enriched both at βmajorand at the upstream locus control region. H3-meK79 was also enriched at the active cad gene, whereas the transcriptionally inactive loci necdin and MyoD1 contained very little H3-meK79. As the pattern of H3-meK79 at the β-globin locus differed between adult and embryonic erythroid cells, establishment and/or maintenance of H3-meK79 was developmentally dynamic. Genetic complementation analysis in null cells lacking the erythroid and megakaryocyte-specific transcription factor p45/NF-E2 showed that p45/NF-E2 preferentially establishes H3-meK79 at theβmajor promoter. These results support a model in which H3-meK79 is strongly enriched in mammalian chromatin at active genes but not uniformly throughout active chromatin domains. As H3-meK79 is highly regulated at the β-globin locus, we propose that the murine ortholog of Disruptor of Telomeric Silencing-1-like (mDOT1L) methyltransferase, which synthesizes H3-meK79, regulates β-globin transcription.

Global identification of protein kinase substrates by protein microarray analysis

Herein, we describe a protocol for the global identification of in vitro substrates targeted by protein kinases using protein microarray technology. Large numbers of fusion proteins tagged at their carboxy-termini are purified in 96-well format and spotted in duplicate onto amino-silane-coated slides in a spatially addressable manner. These arrays are incubated in the presence of purified kinase and radiolabeled ATP, and then washed, dried and analyzed by autoradiography. The extent of phosphorylation of each spot is quantified and normalized, and proteins that are reproducibly phosphorylated in the presence of the active kinase relative to control slides are scored as positive substrates. This approach enables the rapid determination of kinase–substrate relationship on a proteome-wide scale, and although developed using yeast, has since been adapted to higher eukaryotic systems. Expression, purification and printing of the yeast proteome require about 3 weeks. Afterwards, each kinase assay takes ∼3 h to perform.

Genome wide proteomics of ERBB2 and EGFR and other oncogenic pathways in inflammatory breast cancer

In this study we selected three breast cancer cell lines (SKBR3, SUM149 and SUM190) with different oncogene expression levels involved in ERBB2 and EGFR signaling pathways as a model system for the evaluation of selective integration of subsets of transcriptomic and proteomic data. We assessed the oncogene status with reads per kilobase per million mapped reads (RPKM) values for ERBB2 (14.4, 400, and 300 for SUM149, SUM190, and SKBR3, respectively) and for EGFR (60.1, not detected, and 1.4 for the same 3 cell lines). We then used RNA-Seq data to identify those oncogenes with significant transcript levels in these cell lines (total 31) and interrogated the corresponding proteomics data sets for proteins with significant interaction values with these oncogenes. The number of observed interactors for each oncogene showed a significant range, e.g., 4.2% (JAK1) to 27.3% (MYC). The percentage is measured as a fraction of the total protein interactions in a given data set vs total interactors for that oncogene in STRING (Search Tool for the Retrieval of Interacting Genes/Proteins, version 9.0) and I2D (Interologous Interaction Database, version 1.95). This approach allowed us to focus on 4 main oncogenes, ERBB2, EGFR, MYC, and GRB2, for pathway analysis. We used bioinformatics sites GeneGo, PathwayCommons and NCI receptor signaling networks to identify pathways that contained the four main oncogenes and had good coverage in the transcriptomic and proteomic data sets as well as a significant number of oncogene interactors. The four pathways identified were ERBB signaling, EGFR1 signaling, integrin outside-in signaling, and validated targets of C-MYC transcriptional activation. The greater dynamic range of the RNA-Seq values allowed the use of transcript ratios to correlate observed protein values with the relative levels of the ERBB2 and EGFR transcripts in each of the four pathways. This provided us with potential proteomic signatures for the SUM149 and 190 cell lines, growth factor receptor-bound protein 7 (GRB7), Crk-like protein (CRKL) and Catenin delta-1 (CTNND1) for ERBB signaling; caveolin 1 (CAV1), plectin (PLEC) for EGFR signaling; filamin A (FLNA) and actinin alpha1 (ACTN1) (associated with high levels of EGFR transcript) for integrin signalings; branched chain amino-acid transaminase 1 (BCAT1), carbamoyl-phosphate synthetase (CAD), nucleolin (NCL) (high levels of EGFR transcript); transferrin receptor (TFRC), metadherin (MTDH) (high levels of ERBB2 transcript) for MYC signaling; S100-A2 protein (S100A2), caveolin 1 (CAV1), Serpin B5 (SERPINB5), stratifin (SFN), PYD and CARD domain containing (PYCARD), and EPH receptor A2 (EPHA2) for PI3K signaling, p53 subpathway. Future studies of inflammatory breast cancer (IBC), from which the cell lines were derived, will be used to explore the significance of these observations.

Distinct Splice Variants and Pathway Enrichment in the Cell-Line Models of Aggressive Human Breast Cancer Subtypes

This study was conducted as a part of the Chromosome-Centric Human Proteome Project (C-HPP) of the Human Proteome Organization. The United States team of C-HPP is focused on characterizing the protein-coding genes in chromosome 17. Despite its small size, chromosome 17 is rich in protein-coding genes; it contains many cancer-associated genes, including BRCA1, ERBB2, (Her2/neu), and TP53. The goal of this study was to examine the splice variants expressed in three ERBB2 expressed breast cancer cell-line models of hormone-receptor-negative breast cancers by integrating RNA-Seq and proteomic mass spectrometry data. The cell lines represent distinct phenotypic variations subtype: SKBR3 (ERBB2+ (overexpression)/ER-/PR-; adenocarcinoma), SUM190 (ERBB2+ (overexpression)/ER-/PR-; inflammatory breast cancer), and SUM149 (ERBB2 (low expression) ER-/PR-; inflammatory breast cancer). We identified more than one splice variant for 1167 genes expressed in at least one of the three cancer cell lines. We found multiple variants of genes that are in the signaling pathways downstream of ERBB2 along with variants specific to one cancer cell line compared with the other two cancer cell lines and with normal mammary cells. The overall transcript profiles based on read counts indicated more similarities between SKBR3 and SUM190. The top-ranking Gene Ontology and BioCarta pathways for the cell-line specific variants pointed to distinct key mechanisms including: amino sugar metabolism, caspase activity, and endocytosis in SKBR3; different aspects of metabolism, especially of lipids in SUM190; cell-to-cell adhesion, integrin, and ERK1/ERK2 signaling; and translational control in SUM149. The analyses indicated an enrichment in the electron transport chain processes in the ERBB2 overexpressed cell line models and an association of nucleotide binding, RNA splicing, and translation processes with the IBC models, SUM190 and SUM149. Detailed experimental studies on the distinct variants identified from each of these three breast cancer cell line models that may open opportunities for drug target discovery and help unveil their specific roles in cancer progression and metastasis.

Establishment and Regulation of Chromatin Domains: Mechanistic Insights from Studies of Hemoglobin Synthesis

Publisher Summary This chapter discusses progress on dissecting mechanisms underlying the establishment and regulation of the erythroid cell-specific β-globin chromatin domain while emphasizing general principles that have emerged from the work. ‘‘Chromatin domain’’ refers to a broad chromosomal region from approximately several thousand to several hundred thousand base pairs in which a gene or gene cluster resides rather than the classically defined ∼100-kb topologically constrained chromosomal loop identified by microscopy. The chapter reviews predominantly mechanistic studies conducted in mammalian systems, including genetic complementation analysis of trans-acting factor function in biologically relevant murine erythroid cell lines and analysis of cis-element function via targeted deletions in mice. Further work on the compelling problem of how cell type-specific chromatin domains are established and regulated in normal cells and how the underlying mechanisms go awry in disease states promises to be highly rewarding and to yield a continuous stream of focused and serendipitious discoveries for many years to come.

RNA Sequencing Analysis Detection of a Novel Pathway of Endothelial Dysfunction in Pulmonary Arterial Hypertension.

RATIONALE Pulmonary arterial hypertension is characterized by endothelial dysregulation, but global changes in gene expression have not been related to perturbations in function. OBJECTIVES RNA sequencing was used to discriminate changes in transcriptomes of endothelial cells cultured from lungs of patients with idiopathic pulmonary arterial hypertension versus control subjects and to assess the functional significance of major differentially expressed transcripts. METHODS The endothelial transcriptomes from the lungs of seven control subjects and six patients with idiopathic pulmonary arterial hypertension were analyzed. Differentially expressed genes were related to bone morphogenetic protein type 2 receptor (BMPR2) signaling. Those down-regulated were assessed for function in cultured cells and in a transgenic mouse. MEASUREMENTS AND MAIN RESULTS Fold differences in 10 genes were significant (P < 0.05), four increased and six decreased in patients versus control subjects. No patient was mutant for BMPR2. However, knockdown of BMPR2 by siRNA in control pulmonary arterial endothelial cells recapitulated 6 of 10 patient-related gene changes, including decreased collagen IV (COL4A1, COL4A2) and ephrinA1 (EFNA1). Reduction of BMPR2-regulated transcripts was related to decreased β-catenin. Reducing COL4A1, COL4A2, and EFNA1 by siRNA inhibited pulmonary endothelial adhesion, migration, and tube formation. In mice null for the EFNA1 receptor, EphA2, versus control animals, vascular endothelial growth factor receptor blockade and hypoxia caused more severe pulmonary hypertension, judged by elevated right ventricular systolic pressure, right ventricular hypertrophy, and loss of small arteries. CONCLUSIONS The novel relationship between BMPR2 dysfunction and reduced expression of endothelial COL4 and EFNA1 may underlie vulnerability to injury in pulmonary arterial hypertension.