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.

Protein substrates of the arginine methyltransferase Hmt1 identified by proteome arrays.

Arginine methylation on nonhistone proteins is associated with a number of cellular processes including RNA splicing, protein localization, and the formation of protein complexes. In this manuscript, Saccharomyces cerevisiae proteome arrays carrying 4228 proteins were used with an antimethylarginine antibody to first identify 88 putatively arginine-methylated proteins. By treating the arrays with recombinant arginine methyltransferase Hmt1, 42 proteins were found to be possible substrates of this enzyme. Analysis of the putative arginine-methylated proteins revealed that they were predominantly nuclear or nucleolar in localization, consistent with the localization of Hmt1. Many are involved in known methylarginine-associated functions, such as RNA processing and ribonucleoprotein complex biogenesis, yet others are of newer classes, namely RNA/DNA helicases and tRNA-associated proteins. Using ex vivo methylation and MS/MS, a set of 12 proteins (Brr1, Dia4, Hts1, Mpp10, Mrd1, Nug1, Prp43, Rpa43, Rrp43, Spp381, Utp4, and Npl3), including the RNA helicase Prp43 and tRNA ligases Dia4 and Hts1, were all validated as Hmt1 substrates. Interestingly, the majority of these also had human orthologs, or family members, that have been documented elsewhere to carry arginine methylation. These results confirm arginine methylation as a widespread modification and Hmt1 as the major arginine methyltransferase in the S. cerevisiae cell.

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.

Evaluating Common Humoral Responses against Fungal Infections with Yeast Protein Microarrays.

We profiled the global immunoglobulin response against fungal infection by using yeast protein microarrays. Groups of CD-1 mice were infected systemically with human fungal pathogens (Coccidioides posadasii, Candida albicans, or Paracoccidioides brasiliensis) or inoculated with PBS as a control. Another group was inoculated with heat-killed yeast (HKY) of Saccharomyces cerevisiae. After 30 days, serum from mice in the groups were collected and used to probe S. cerevisiae protein microarrays containing 4800 full-length glutathione S-transferase (GST)-fusion proteins. Antimouse IgG conjugated with Alexafluor 555 and anti-GST antibody conjugated with Alexafluor 647 were used to detect antibody-antigen interactions and the presence of GST-fusion proteins, respectively. Serum after infection with C. albicans reacted with 121 proteins: C. posadasii, 81; P. brasiliensis, 67; and after HKY, 63 proteins on the yeast protein microarray, respectively. We identified a set of 16 antigenic proteins that were shared across the three fungal pathogens. These include retrotransposon capsid proteins, heat shock proteins, and mitochondrial proteins. Five of these proteins were identified in our previous study of fungal cell wall by mass spectrometry (Ann. N. Y. Acad. Sci. 2012, 1273, 44-51). The results obtained give a comprehensive view of the immunological responses to fungal infections at the proteomic level. They also offer insight into immunoreactive protein commonality among several fungal pathogens and provide a basis for a panfungal vaccine.

Integrated Proteomic and Genomic Analysis of Gastric Cancer Patient Tissues.

V-erb-b2 erythroblastic leukemia viral oncogene homologue 2, known as ERBB2, is an important oncogene in the development of certain cancers. It can form a heterodimer with other epidermal growth factor receptor family members and activate kinase-mediated downstream signaling pathways. ERBB2 gene is located on chromosome 17 and is amplified in a subset of cancers, such as breast, gastric, and colon cancer. Of particular interest to the Chromosome-Centric Human Proteome Project (C-HPP) initiative is the amplification mechanism that typically results in overexpression of a set of genes adjacent to ERBB2, which provides evidence of a linkage between gene location and expression. In this report we studied patient samples from ERBB2-positive together with adjacent control nontumor tissues. In addition, non-ERBB2-expressing patient samples were selected as comparison to study the effect of expression of this oncogene. We detected 196 proteins in ERBB2-positive patient tumor samples that had minimal overlap (29 proteins) with the non-ERBB2 tumor samples. Interaction and pathway analysis identified extracellular signal regulated kinase (ERK) cascade and actin polymerization and actinmyosin assembly contraction as pathways of importance in ERBB2+ and ERBB2- gastric cancer samples, respectively. The raw data files are deposited at ProteomeXchange (identifier: PXD002674) as well as GPMDB.

Identification of potential glycan cancer markers with sialic acid attached to sialic acid and up-regulated fucosylated galactose structures in epidermal growth factor receptor secreted from A431 cell line.

We have used powerful HPLC-mass spectrometric approaches to characterize the secreted form of epidermal growth factor receptor (sEGFR). We demonstrated that the amino acid sequence lacked the cytoplasmic domain and was consistent with the primary sequence reported for EGFR purified from a human plasma pool. One of the sEGFR forms, attributed to the alternative RNA splicing, was also confirmed by transcriptional analysis (RNA sequencing). Two unusual types of glycan structures were observed in sEGFR as compared with membrane-bound EGFR from the A431 cell line. The unusual glycan structures were di-sialylated glycans (sialic acid attached to sialic acid) at Asn-151 and N-acetylhexosamine attached to a branched fucosylated galactose with N-acetylglucosamine moieties (HexNAc-(Fuc)Gal-GlcNAc) at Asn-420. These unusual glycans at specific sites were either present at a much lower level or were not observable in membrane-bound EGFR present in the A431 cell lysate. The observation of these di-sialylated glycan structures was consistent with the observed expression of the corresponding α-N-acetylneuraminide α-2,8-sialyltransferase 2 (ST8SiA2) and α-N-acetylneuraminide α-2,8-sialyltransferase 4 (ST8SiA4), by quantitative real time RT-PCR. The connectivity present at the branched fucosylated galactose was also confirmed by methylation of the glycans followed by analysis with sequential fragmentation in mass spectrometry. We hypothesize that the presence of such glycan structures could promote secretion via anionic or steric repulsion mechanisms and thus facilitate the observation of these glycan forms in the secreted fractions. We plan to use this model system to facilitate the search for novel glycan structures present at specific sites in sEGFR as well as other secreted oncoproteins such as Erbb2 as markers of disease progression in blood samples from cancer patients.

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.

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.

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.

A chromosome-centric human proteome project (C-HPP) to characterize the sets of proteins encoded in chromosome 17.

We report progress assembling the parts list for chromosome 17 and illustrate the various processes that we have developed to integrate available data from diverse genomic and proteomic knowledge bases. As primary resources, we have used GPMDB, neXtProt, PeptideAtlas, Human Protein Atlas (HPA), and GeneCards. All sites share the common resource of Ensembl for the genome modeling information. We have defined the chromosome 17 parts list with the following information: 1169 protein-coding genes, the numbers of proteins confidently identified by various experimental approaches as documented in GPMDB, neXtProt, PeptideAtlas, and HPA, examples of typical data sets obtained by RNASeq and proteomic studies of epithelial derived tumor cell lines (disease proteome) and a normal proteome (peripheral mononuclear cells), reported evidence of post-translational modifications, and examples of alternative splice variants (ASVs). We have constructed a list of the 59 "missing" proteins as well as 201 proteins that have inconclusive mass spectrometric (MS) identifications. In this report we have defined a process to establish a baseline for the incorporation of new evidence on protein identification and characterization as well as related information from transcriptome analyses. This initial list of "missing" proteins that will guide the selection of appropriate samples for discovery studies as well as antibody reagents. Also we have illustrated the significant diversity of protein variants (including post-translational modifications, PTMs) using regions on chromosome 17 that contain important oncogenes. We emphasize the need for mandated deposition of proteomics data in public databases, the further development of improved PTM, ASV, and single nucleotide variant (SNV) databases, and the construction of Web sites that can integrate and regularly update such information. In addition, we describe the distribution of both clustered and scattered sets of protein families on the chromosome. Since chromosome 17 is rich in cancer-associated genes, we have focused the clustering of cancer-associated genes in such genomic regions and have used the ERBB2 amplicon as an example of the value of a proteogenomic approach in which one integrates transcriptomic with proteomic information and captures evidence of coexpression through coordinated regulation.

Genetic regulation of OAS1 nonsense-mediated decay underlies association with COVID-19 hospitalization in patients of European and African ancestries.

The chr12q24.13 locus encoding OAS1-OAS3 antiviral proteins has been associated with coronavirus disease 2019 (COVID-19) susceptibility. Here, we report genetic, functional and clinical insights into this locus in relation to COVID-19 severity. In our analysis of patients of European (n = 2,249) and African (n = 835) ancestries with hospitalized versus nonhospitalized COVID-19, the risk of hospitalized disease was associated with a common OAS1 haplotype, which was also associated with reduced severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) clearance in a clinical trial with pegIFN-λ1. Bioinformatic analyses and in vitro studies reveal the functional contribution of two associated OAS1 exonic variants comprising the risk haplotype. Derived human-specific alleles rs10774671-A and rs1131454 -A decrease OAS1 protein abundance through allele-specific regulation of splicing and nonsense-mediated decay (NMD). We conclude that decreased OAS1 expression due to a common haplotype contributes to COVID-19 severity. Our results provide insight into molecular mechanisms through which early treatment with interferons could accelerate SARS-CoV-2 clearance and mitigate against severe COVID-19.

Clonal hematopoiesis is associated with risk of severe Covid-19.

Acquired somatic mutations in hematopoietic stem and progenitor cells (clonal hematopoiesis or CH) are associated with advanced age, increased risk of cardiovascular and malignant diseases, and decreased overall survival. These adverse sequelae may be mediated by altered inflammatory profiles observed in patients with CH. A pro-inflammatory immunologic profile is also associated with worse outcomes of certain infections, including SARS-CoV-2 and its associated disease Covid-19. Whether CH predisposes to severe Covid-19 or other infections is unknown. Among 525 individuals with Covid-19 from Memorial Sloan Kettering (MSK) and the Korean Clonal Hematopoiesis (KoCH) consortia, we show that CH is associated with severe Covid-19 outcomes (OR = 1.85, 95%=1.15-2.99, p = 0.01), in particular CH characterized by non-cancer driver mutations (OR = 2.01, 95% CI = 1.15-3.50, p = 0.01). We further explore the relationship between CH and risk of other infections in 14,211 solid tumor patients at MSK. CH is significantly associated with risk of Clostridium Difficile (HR = 2.01, 95% CI: 1.22-3.30, p = 6×10-3) and Streptococcus/Enterococcus infections (HR = 1.56, 95% CI = 1.15-2.13, p = 5×10-3). These findings suggest a relationship between CH and risk of severe infections that warrants further investigation.

Genetic regulation of OAS1 nonsense-mediated decay underlies association with risk of severe COVID-19.

Genomic regions have been associated with COVID-19 susceptibility and outcomes, including the chr12q24.13 locus encoding antiviral proteins OAS1-3. Here, we report genetic, functional, and clinical insights into genetic associations within this locus. In Europeans, the risk of hospitalized vs. non-hospitalized COVID-19 was associated with a single 19Kb-haplotype comprised of 76 OAS1 variants included in a 95% credible set within a large genomic fragment introgressed from Neandertals. The risk haplotype was also associated with impaired spontaneous but not treatment-induced SARS-CoV-2 clearance in a clinical trial with pegIFN-λ1. We demonstrate that two exonic variants, rs10774671 and rs1131454, affect splicing and nonsense-mediated decay of OAS1 . We suggest that genetically-regulated loss of OAS1 expression contributes to impaired spontaneous clearance of SARS-CoV-2 and elevated risk of hospitalization for COVID-19. Our results provide the rationale for further clinical studies using interferons to compensate for impaired spontaneous SARS-CoV-2 clearance, particularly in carriers of the OAS1 risk haplotypes.

H3 K79 dimethylation marks developmental activation of the beta-globin gene but is reduced upon LCR-mediated high-level transcription.

Genome-wide analyses of the relationship between H3 K79 dimethylation and transcription have revealed contradictory results. To clarify this relationship at a single locus, we analyzed expression and H3 K79 modification levels of wild-type (WT) and transcriptionally impaired beta-globin mutant genes during erythroid differentiation. Analysis of fractionated erythroid cells derived from WT/Delta locus control region (LCR) heterozygous mice reveals no significant H3 K79 dimethylation of the beta-globin gene on either allele prior to activation of transcription. Upon transcriptional activation, H3 K79 di-methylation is observed along both WT and DeltaLCR alleles, and both alleles are located in proximity to H3 K79 dimethylation nuclear foci. However, H3 K79 di-methylation is significantly increased along the DeltaLCR allele compared with the WT allele. In addition, analysis of a partial LCR deletion mutant reveals that H3 K79 dimethylation is inversely correlated with beta-globin gene expression levels. Thus, while our results support a link between H3 K79 dimethylation and gene expression, high levels of this mark are not essential for high level beta-globin gene transcription. We propose that H3 K79 dimethylation is destabilized on a highly transcribed template.

Clonal hematopoiesis is associated with risk of severe Covid-19.

Acquired somatic mutations in hematopoietic stem and progenitor cells (clonal hematopoiesis or CH) are associated with advanced age, increased risk of cardiovascular and malignant diseases, and decreased overall survival. 1-4 These adverse sequelae may be mediated by altered inflammatory profiles observed in patients with CH. 2,5,6 A pro-inflammatory immunologic profile is also associated with worse outcomes of certain infections, including SARS-CoV-2 and its associated disease Covid-19. 7,8 Whether CH predisposes to severe Covid-19 or other infections is unknown. Among 515 individuals with Covid-19 from Memorial Sloan Kettering (MSK) and the Korean Clonal Hematopoiesis (KoCH) consortia, we found that CH was associated with severe Covid-19 outcomes (OR=1.9, 95%=1.2-2.9, p=0.01). We further explored the relationship between CH and risk of other infections in 14,211 solid tumor patients at MSK. CH was significantly associated with risk of Clostridium Difficile (HR=2.0, 95% CI: 1.2-3.3, p=6×10 -3 ) and Streptococcus/Enterococcus infections (HR=1.5, 95% CI=1.1-2.1, p=5×10 -3 ). These findings suggest a relationship between CH and risk of severe infections that warrants further investigation.

Interpretation of EBV infection in pan-cancer genome considering viral life cycle: LiEB (Life cycle of Epstein-Barr virus).

We report a novel transcriptomic analysis workflow called LiEB (Life cycle of Epstein-Barr virus) to characterize distributions of oncogenic virus, Epstein-Barr virus (EBV) infection in human tumors. We analyzed 851 The Cancer Genome Atlas whole-transcriptome sequencing (WTS) data to investigate EBV infection by life cycle information using three-step LiEB workflow: 1) characterize virus infection generally; 2) align transcriptome sequences against a hybrid human-EBV genome, and 3) quantify EBV gene expression. Our results agreed with EBV infection status of public cell line data. Analysis in stomach adenocarcinoma identified EBV-positive cases involving PIK3CA mutations and/or CDKN2A silencing with biologically more determination, compared to previous reports. In this study, we found that a small number of colorectal adenocarcinoma cases involved with EBV lytic gene expression. Expression of EBV lytic genes was also observed in 3% of external colon cancer cohort upon WTS analysis. Gene set enrichment analysis showed elevated expression of genes related to E2F targeting and interferon-gamma responses in EBV-associated tumors. Finally, we suggest that interpretation of EBV life cycle is essential when analyzing its infection in tumors, and LiEB provides high capability of detecting EBV-positive tumors. Observation of EBV lytic gene expression in a subset of colon cancers warrants further research.

Distinct transcriptomic and exomic abnormalities within myelodysplastic syndrome marrow cells.

To provide biologic insights into mechanisms underlying myelodysplastic syndromes (MDS) we evaluated the CD34+ marrow cells transcriptome using high-throughput RNA sequencing (RNA-Seq). We demonstrated significant differential gene expression profiles (GEPs) between MDS and normal and identified 41 disease classifier genes. Additionally, two main clusters of GEPs distinguished patients based on their major clinical features, particularly between those whose disease remained stable versus patients who transformed into acute myeloid leukemia within 12 months. The genes whose expression was associated with disease outcome were involved in functional pathways and biologic processes highly relevant for MDS. Combined with exomic analysis we identified differential isoform usage of genes in MDS mutational subgroups, with consequent dysregulation of distinct biologic functions. This combination of clinical, transcriptomic and exomic findings provides valuable understanding of mechanisms underlying MDS and its progression to a more aggressive stage and also facilitates prognostic characterization of MDS patients.

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

RNA polymerase II (Pol II) can associate with regulatory elements far from promoters. For the murine beta-globin locus, Pol II binds the beta-globin locus control region (LCR) far upstream of the beta-globin promoters, independent of recruitment to and activation of the betamajor 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 betamajor 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 beta-globin locus, indicating that Pol II tracking is not required to maintain the pattern.

Transcriptomic and epigenomic differences in human induced pluripotent stem cells generated from six reprogramming methods.

Many reprogramming methods can generate human induced pluripotent stem cells (hiPSCs) that closely resemble human embryonic stem cells (hESCs). This has led to assessments of how similar hiPSCs are to hESCs, by evaluating differences in gene expression, epigenetic marks and differentiation potential. However, all previous studies were performed using hiPSCs acquired from different laboratories, passage numbers, culturing conditions, genetic backgrounds and reprogramming methods, all of which may contribute to the reported differences. Here, by using high-throughput sequencing under standardized cell culturing conditions and passage number, we compare the epigenetic signatures (H3K4me3, H3K27me3 and HDAC2 ChIP-seq profiles) and transcriptome differences (by RNA-seq) of hiPSCs generated from the same primary fibroblast population by using six different reprogramming methods. We found that the reprogramming method impacts the resulting transcriptome and that all hiPSC lines could terminally differentiate, regardless of the reprogramming method. Moreover, by comparing the differences between the hiPSC and hESC lines, we observed a significant proportion of differentially expressed genes that could be attributed to polycomb repressive complex targets.