Roy Williams

Regulatory networks define phenotypic classes of human stem cell lines

Stem cells are defined as self-renewing cell populations that can differentiate into multiple distinct cell types. However, hundreds of different human cell lines from embryonic, fetal and adult sources have been called stem cells, even though they range from pluripotent cells—typified by embryonic stem cells, which are capable of virtually unlimited proliferation and differentiation—to adult stem cell lines, which can generate a far more limited repertoire of differentiated cell types. The rapid increase in reports of new sources of stem cells and their anticipated value to regenerative medicine has highlighted the need for a general, reproducible method for classification of these cells. We report here the creation and analysis of a database of global gene expression profiles (which we call the ‘stem cell matrix’) that enables the classification of cultured human stem cells in the context of a wide variety of pluripotent, multipotent and differentiated cell types. Using an unsupervised clustering method to categorize a collection of ∼150 cell samples, we discovered that pluripotent stem cell lines group together, whereas other cell types, including brain-derived neural stem cell lines, are very diverse. Using further bioinformatic analysis we uncovered a protein–protein network (PluriNet) that is shared by the pluripotent cells (embryonic stem cells, embryonal carcinomas and induced pluripotent cells). Analysis of published data showed that the PluriNet seems to be a common characteristic of pluripotent cells, including mouse embryonic stem and induced pluripotent cells and human oocytes. Our results offer a new strategy for classifying stem cells and support the idea that pluripotency and self-renewal are under tight control by specific molecular networks.

The functional subunit of a dimeric transcription activator protein depends on promoter architecture.

In Class I CAP‐dependent promoters, the DNA site for CAP is located upstream of the DNA site for RNA polymerase. In Class II CAP‐dependent promoters, the DNA site for CAP overlaps the DNA site for RNA polymerase, replacing the ‐35 site. We have used an ‘oriented heterodimers’ approach to identify the functional subunit of CAP at two Class I promoters having different distances between the DNA sites for CAP and RNA polymerase [CC(‐61.5) and CC(‐72.5)] and at one Class II promoter [CC(‐41.5)]. Our results indicate that transcription activation at Class I promoters, irrespective of the distance between the DNA sites for CAP and RNA polymerase, requires the activating region of the promoter‐proximal subunit of CAP. In striking contrast, our results indicate that transcription activation at Class II promoters requires the activating region of the promoter‐distal subunit of CAP.

Interactions between the Escherichia coli cyclic AMP receptor protein and RNA polymerase at class II promoters.

The effects of a number of mutations in crp have been measured at different cyclic AMP receptor protein (CRP)‐dependent Class II promoters, where the CRP‐binding site is centred around 411/2 base pairs upstream from the transcription start point. The amino acid substitutions HL159 and TA158 result in reduced CRP‐dependent activation, but the reduction varies from one Class II promoter to another. Deletions in the C‐terminus of the RNA polymerase alpha subunit suppress the effects of HL159 and TA158. The role of the C‐terminus of alpha at these promoters is assessed. Other changes at E58, K52 and E96 affect CRP activity specifically at Class II promoters and their role is discussed.

HDAC inhibition imparts beneficial transgenerational effects in Huntington's disease mice via altered DNA and histone methylation

Significance We demonstrate that histone deacetylase (HDAC) inhibition can elicit changes in DNA methylation in Huntington’s disease (HD) human fibroblasts, as well as in sperm from HD transgenic mice, in association with DNA methylation-related gene expression changes. We suggest that alterations in sperm DNA methylation lead to transgenerational effects, and, accordingly, we show that first filial generation (F1) offspring of HDAC inhibitor-treated male HD transgenic mice show improved HD disease phenotypes compared with F1 offspring from vehicle-treated male HD transgenic mice. These findings have significant implications for human health because they enforce the concept that ancestral drug exposure may be a major molecular factor that can affect disease phenotypes, yet in a positive manner. Further, we implicate Lys (K)-specific demethylase 5d expression in this phenomenon. Increasing evidence has demonstrated that epigenetic factors can profoundly influence gene expression and, in turn, influence resistance or susceptibility to disease. Epigenetic drugs, such as histone deacetylase (HDAC) inhibitors, are finding their way into clinical practice, although their exact mechanisms of action are unclear. To identify mechanisms associated with HDAC inhibition, we performed microarray analysis on brain and muscle samples treated with the HDAC1/3-targeting inhibitor, HDACi 4b. Pathways analyses of microarray datasets implicate DNA methylation as significantly associated with HDAC inhibition. Further assessment of DNA methylation changes elicited by HDACi 4b in human fibroblasts from normal controls and patients with Huntington’s disease (HD) using the Infinium HumanMethylation450 BeadChip revealed a limited, but overlapping, subset of methylated CpG sites that were altered by HDAC inhibition in both normal and HD cells. Among the altered loci of Y chromosome-linked genes, KDM5D, which encodes Lys (K)-specific demethylase 5D, showed increased methylation at several CpG sites in both normal and HD cells, as well as in DNA isolated from sperm from drug-treated male mice. Further, we demonstrate that first filial generation (F1) offspring from drug-treated male HD transgenic mice show significantly improved HD disease phenotypes compared with F1 offspring from vehicle-treated male HD transgenic mice, in association with increased Kdm5d expression, and decreased histone H3 Lys4 (K4) (H3K4) methylation in the CNS of male offspring. Additionally, we show that overexpression of Kdm5d in mutant HD striatal cells significantly improves metabolic deficits. These findings indicate that HDAC inhibitors can elicit transgenerational effects, via cross-talk between different epigenetic mechanisms, to have an impact on disease phenotypes in a beneficial manner.

Whole-genome mutational burden analysis of three pluripotency induction methods

There is concern that the stresses of inducing pluripotency may lead to deleterious DNA mutations in induced pluripotent stem cell (iPSC) lines, which would compromise their use for cell therapies. Here we report comparative genomic analysis of nine isogenic iPSC lines generated using three reprogramming methods: integrating retroviral vectors, non-integrating Sendai virus and synthetic mRNAs. We used whole-genome sequencing and de novo genome mapping to identify single-nucleotide variants, insertions and deletions, and structural variants. Our results show a moderate number of variants in the iPSCs that were not evident in the parental fibroblasts, which may result from reprogramming. There were only small differences in the total numbers and types of variants among different reprogramming methods. Most importantly, a thorough genomic analysis showed that the variants were generally benign. We conclude that the process of reprogramming is unlikely to introduce variants that would make the cells inappropriate for therapy.

The ACTCellerate initiative: large-scale combinatorial cloning of novel human embryonic stem cell derivatives

Human embryonic stem cells offer a scalable and renewable source of all somatic cell types. Human embryonic progenitor (hEP) cells are partially differentiated endodermal, mesodermal and ectodermal cell types that have not undergone terminal differentiation and express an embryonic pattern of gene expression. Here, we describe a large-scale and reproducible method of isolating a diverse library of clonally purified hEP cell lines, many of which are capable of extended propagation in vitro. Initial microarray and non-negative matrix factorization gene-expression profiling suggests that the library consists of at least 140 distinct clones and contains many previously uncharacterized cell types derived from all germ layers that display diverse embryo- and site-specific homeobox gene expression. Despite the expression of many oncofetal genes, none of the hEP cell lines tested led to tumor formation when transplanted into immunocompromised mice. All hEP lines studied appear to have a finite replicative lifespan but have longer telomeres than most fetal- or adult-derived cells, thereby facilitating their use in the manufacture of purified lineages for research and human therapy.

iPSCORE: A Resource of 222 iPSC Lines Enabling Functional Characterization of Genetic Variation across a Variety of Cell Types.

Summary Large-scale collections of induced pluripotent stem cells (iPSCs) could serve as powerful model systems for examining how genetic variation affects biology and disease. Here we describe the iPSCORE resource: a collection of systematically derived and characterized iPSC lines from 222 ethnically diverse individuals that allows for both familial and association-based genetic studies. iPSCORE lines are pluripotent with high genomic integrity (no or low numbers of somatic copy-number variants) as determined using high-throughput RNA-sequencing and genotyping arrays, respectively. Using iPSCs from a family of individuals, we show that iPSC-derived cardiomyocytes demonstrate gene expression patterns that cluster by genetic background, and can be used to examine variants associated with physiological and disease phenotypes. The iPSCORE collection contains representative individuals for risk and non-risk alleles for 95% of SNPs associated with human phenotypes through genome-wide association studies. Our study demonstrates the utility of iPSCORE for examining how genetic variants influence molecular and physiological traits in iPSCs and derived cell lines.

Mapping the malaria parasite druggable genome by using in vitro evolution and chemogenomics

Dissecting Plasmodium drug resistance Malaria is a deadly disease with no effective vaccine. Physicians thus depend on antimalarial drugs to save lives, but such compounds are often rendered ineffective when parasites evolve resistance. Cowell et al. systematically studied patterns of Plasmodium falciparum genome evolution by analyzing the sequences of clones that were resistant to diverse antimalarial compounds across the P. falciparum life cycle (see the Perspective by Carlton). The findings identify hitherto unrecognized drug targets and drug-resistance genes, as well as additional alleles in known drug-resistance genes. Science, this issue p. 191; see also p. 159 Genome sequencing elucidates potential drug resistance in the malaria parasite and identifies antimalarial targets. Chemogenetic characterization through in vitro evolution combined with whole-genome analysis can identify antimalarial drug targets and drug-resistance genes. We performed a genome analysis of 262 Plasmodium falciparum parasites resistant to 37 diverse compounds. We found 159 gene amplifications and 148 nonsynonymous changes in 83 genes associated with drug-resistance acquisition, where gene amplifications contributed to one-third of resistance acquisition events. Beyond confirming previously identified multidrug-resistance mechanisms, we discovered hitherto unrecognized drug target–inhibitor pairs, including thymidylate synthase and a benzoquinazolinone, farnesyltransferase and a pyrimidinedione, and a dipeptidylpeptidase and an arylurea. This exploration of the P. falciparum resistome and druggable genome will likely guide drug discovery and structural biology efforts, while also advancing our understanding of resistance mechanisms available to the malaria parasite.

High-Throughput and Cost-Effective Characterization of Induced Pluripotent Stem Cells

Summary Reprogramming somatic cells to induced pluripotent stem cells (iPSCs) offers the possibility of studying the molecular mechanisms underlying human diseases in cell types difficult to extract from living patients, such as neurons and cardiomyocytes. To date, studies have been published that use small panels of iPSC-derived cell lines to study monogenic diseases. However, to study complex diseases, where the genetic variation underlying the disorder is unknown, a sizable number of patient-specific iPSC lines and controls need to be generated. Currently the methods for deriving and characterizing iPSCs are time consuming, expensive, and, in some cases, descriptive but not quantitative. Here we set out to develop a set of simple methods that reduce cost and increase throughput in the characterization of iPSC lines. Specifically, we outline methods for high-throughput quantification of surface markers, gene expression analysis of in vitro differentiation potential, and evaluation of karyotype with markedly reduced cost.

Assessment of microRNA and gene dysregulation in pulmonary hypertension by endoarterial biopsy

MicroRNAs (miRNAs) may regulate a number of genes, each of which may have a variety of functions. We utilized an endoarterial biopsy catheter to assess the dysregulation of miRNAs in a porcine shunt model of pulmonary hypertension (PH). Two Yucatan micropigs underwent surgical anastomosis of the left pulmonary artery to the descending aorta. Endoarterial biopsy samples were obtained at baseline, and at regular intervals during the progression of PH. RNA, isolated from biopsy samples, was analyzed by Illumina miRNA expression microarrays (containing ∼1200 human miRNAs), Affymetrix Porcine GeneChips, Bioconductor, and GeneSpring. We examined a total of 925 genes in a PH whole genome microarray. Biopsy samples showed that 39 miRNAs were downregulated and 34 miRNAs were upregulated compared to baseline. The number of PH-associated genes reported to be controlled by each of the dysregulated miRNAs was in the range of 1–113. The five miRNAs that had the largest number of PH-associated genes were: miR-548c-3p, miR-520d-3p, miR-130a-5p, miR-30a-3p, and miR-let-7g-3p. Several of the dysregulated miRNAs have been associated with molecular pathways and biologic processes involved in PH. Among 29 miRNAs, which were predicted to be dysregulated by a systems biology approach, we found four that were dysregulated in our porcine shunt model. An endoarterial biopsy technique was successful in showing that a large number of miRNAs are dysregulated in a porcine shunt model of PH. Many of these miRNAs control multiple PH-associated genes, molecular pathways, and biologic processes. Endoarterial biopsy offers potential experimental and clinical diagnostic value.