Eric T. Wang

Alternative isoform regulation in human tissue transcriptomes

Through alternative processing of pre-messenger RNAs, individual mammalian genes often produce multiple mRNA and protein isoforms that may have related, distinct or even opposing functions. Here we report an in-depth analysis of 15 diverse human tissue and cell line transcriptomes on the basis of deep sequencing of complementary DNA fragments, yielding a digital inventory of gene and mRNA isoform expression. Analyses in which sequence reads are mapped to exon–exon junctions indicated that 92–94% of human genes undergo alternative splicing, ∼86% with a minor isoform frequency of 15% or more. Differences in isoform-specific read densities indicated that most alternative splicing and alternative cleavage and polyadenylation events vary between tissues, whereas variation between individuals was approximately twofold to threefold less common. Extreme or ‘switch-like’ regulation of splicing between tissues was associated with increased sequence conservation in regulatory regions and with generation of full-length open reading frames. Patterns of alternative splicing and alternative cleavage and polyadenylation were strongly correlated across tissues, suggesting coordinated regulation of these processes, and sequence conservation of a subset of known regulatory motifs in both alternative introns and 3′ untranslated regions suggested common involvement of specific factors in tissue-level regulation of both splicing and polyadenylation.

An Abundance of Ubiquitously Expressed Genes Revealed by Tissue Transcriptome Sequence Data

The parts of the genome transcribed by a cell or tissue reflect the biological processes and functions it carries out. We characterized the features of mammalian tissue transcriptomes at the gene level through analysis of RNA deep sequencing (RNA-Seq) data across human and mouse tissues and cell lines. We observed that roughly 8,000 protein-coding genes were ubiquitously expressed, contributing to around 75% of all mRNAs by message copy number in most tissues. These mRNAs encoded proteins that were often intracellular, and tended to be involved in metabolism, transcription, RNA processing or translation. In contrast, genes for secreted or plasma membrane proteins were generally expressed in only a subset of tissues. The distribution of expression levels was broad but fairly continuous: no support was found for the concept of distinct expression classes of genes. Expression estimates that included reads mapping to coding exons only correlated better with qRT-PCR data than estimates which also included 3′ untranslated regions (UTRs). Muscle and liver had the least complex transcriptomes, in that they expressed predominantly ubiquitous genes and a large fraction of the transcripts came from a few highly expressed genes, whereas brain, kidney and testis expressed more complex transcriptomes with the vast majority of genes expressed and relatively small contributions from the most expressed genes. mRNAs expressed in brain had unusually long 3′UTRs, and mean 3′UTR length was higher for genes involved in development, morphogenesis and signal transduction, suggesting added complexity of UTR-based regulation for these genes. Our results support a model in which variable exterior components feed into a large, densely connected core composed of ubiquitously expressed intracellular proteins.

Integration of flow-dependent endothelial phenotypes by Kruppel-like factor 2

In the face of systemic risk factors, certain regions of the arterial vasculature remain relatively resistant to the development of atherosclerotic lesions. The biomechanically distinct environments in these arterial geometries exert a protective influence via certain key functions of the endothelial lining; however, the mechanisms underlying the coordinated regulation of specific mechano-activated transcriptional programs leading to distinct endothelial functional phenotypes have remained elusive. Here, we show that the transcription factor Kruppel-like factor 2 (KLF2) is selectively induced in endothelial cells exposed to a biomechanical stimulus characteristic of atheroprotected regions of the human carotid and that this flow-mediated increase in expression occurs via a MEK5/ERK5/MEF2 signaling pathway. Overexpression and silencing of KLF2 in the context of flow, combined with findings from genome-wide analyses of gene expression, demonstrate that the induction of KLF2 results in the orchestrated regulation of endothelial transcriptional programs controlling inflammation, thrombosis/hemostasis, vascular tone, and blood vessel development. Our data also indicate that KLF2 expression globally modulates IL-1beta-mediated endothelial activation. KLF2 therefore serves as a mechano-activated transcription factor important in the integration of multiple endothelial functions associated with regions of the arterial vasculature that are relatively resistant to atherogenesis.

Recent acceleration of human adaptive evolution

Genomic surveys in humans identify a large amount of recent positive selection. Using the 3.9-million HapMap SNP dataset, we found that selection has accelerated greatly during the last 40,000 years. We tested the null hypothesis that the observed age distribution of recent positively selected linkage blocks is consistent with a constant rate of adaptive substitution during human evolution. We show that a constant rate high enough to explain the number of recently selected variants would predict (i) site heterozygosity at least 10-fold lower than is observed in humans, (ii) a strong relationship of heterozygosity and local recombination rate, which is not observed in humans, (iii) an implausibly high number of adaptive substitutions between humans and chimpanzees, and (iv) nearly 100 times the observed number of high-frequency linkage disequilibrium blocks. Larger populations generate more new selected mutations, and we show the consistency of the observed data with the historical pattern of human population growth. We consider human demographic growth to be linked with past changes in human cultures and ecologies. Both processes have contributed to the extraordinarily rapid recent genetic evolution of our species.

Selective modulation of human natural killer cells in vivo after prolonged infusion of low dose recombinant interleukin 2.

The immunologic consequences of prolonged infusions of rIL-2 in doses that produce physiologic serum concentrations of this cytokine were investigated. rIL-2 in doses of 0.5-6.0 x 10(6) U/m2 per d (3.3-40 micrograms/m2 per d) was administered by continuous intravenous infusion for 90 consecutive days to patients with advanced cancer. IL-2 concentrations (25 +/- 25 and 77 +/- 64 pM, respectively) that selectively saturate high-affinity IL-2 receptors (IL-2R) were achieved in the serum of patients receiving rIL-2 infusions of 10 micrograms/m2 per d and 30 micrograms/m2 per d. A gradual, progressive expansion of natural killer (NK) cells was seen in the peripheral blood of these patients with no evidence of a plateau effect during the 3 mo of therapy. A preferential expansion of CD56bright NK cells was consistently evident. NK cytotoxicity against tumor targets was only slightly enhanced at these dose levels. However, brief incubation of these expanded NK cells with IL-2 in vitro induced potent lysis of NK-sensitive, NK-resistant, and antibody-coated targets. Infusions of rIL-2 at 40 micrograms/m2 per d produced serum IL-2 levels (345 +/- 381 pM) sufficient to engage intermediate affinity IL-2R p75, which is constitutively expressed by human NK cells. This did not result in greater NK cell expansion compared to the lower dose levels, but did produce in vivo activation of NK cytotoxicity, as evidenced by lysis of NK-resistant targets. There was no consistent change in the numbers of CD56- CD3+ T cells, CD56+ CD3+ MHC-unrestricted T cells, or B cells during infusions of rIL-2 at any of the dosages used. This study demonstrates that prolonged infusions of rIL-2 in doses that saturate only high affinity IL-2R can selectively expand human NK cells for an extended period of time with only minimal toxicity. Further activation of NK cytolytic activity can also be achieved in vivo, but it requires concentrations of IL-2 that bind intermediate affinity IL-2R p75. Clinical trials are underway attempting to exploit the differing effects of various concentrations of IL-2 on human NK cells in vivo.

MBNL proteins repress ES-cell-specific alternative splicing and reprogramming

Previous investigations of the core gene regulatory circuitry that controls the pluripotency of embryonic stem (ES) cells have largely focused on the roles of transcription, chromatin and non-coding RNA regulators. Alternative splicing represents a widely acting mode of gene regulation, yet its role in regulating ES-cell pluripotency and differentiation is poorly understood. Here we identify the muscleblind-like RNA binding proteins, MBNL1 and MBNL2, as conserved and direct negative regulators of a large program of cassette exon alternative splicing events that are differentially regulated between ES cells and other cell types. Knockdown of MBNL proteins in differentiated cells causes switching to an ES-cell-like alternative splicing pattern for approximately half of these events, whereas overexpression of MBNL proteins in ES cells promotes differentiated-cell-like alternative splicing patterns. Among the MBNL-regulated events is an ES-cell-specific alternative splicing switch in the forkhead family transcription factor FOXP1 that controls pluripotency. Consistent with a central and negative regulatory role for MBNL proteins in pluripotency, their knockdown significantly enhances the expression of key pluripotency genes and the formation of induced pluripotent stem cells during somatic cell reprogramming.

High prevalence of rare dopamine receptor D4 alleles in children diagnosed with attention-deficit hyperactivity disorder

Associations have been reported of the 7-repeat (7R) allele of the human dopamine receptor D4 (DRD4) gene with both the personality trait of novelty seeking and attention-deficit/hyperactivity disorder (ADHD). The increased prevalence of the 7R allele in ADHD probands is consistent with the common variant–common disorder hypothesis, which proposes that the high frequency of many complex genetic disorders is related to common DNA variants. Recently, based on the unusual DNA sequence organization and strong linkage disequilibrium surrounding the DRD4 7R allele, we proposed that this allele originated as a rare mutational event, which nevertheless increased to high prevalence in human populations by positive selection. We have now determined, by DNA resequencing of 250 DRD4 alleles obtained from 132 ADHD probands, that most ADHD 7R alleles are of the conserved haplotype found in our previous 600 allele worldwide DNA sample. Interestingly, however, half of the 24 haplotypes uncovered in ADHD probands were novel (not one of the 56 haplotypes found in our prior population studies). Over 10 percent of the ADHD probands had these novel haplotypes, most of which were 7R allele derived. The probability that this high incidence of novel alleles occurred by chance in our ADHD sample is much less than 0.0001. These results suggest that allelic heterogeneity at the DRD4 locus may also contribute to the observed association with ADHD.

Global regulation of alternative splicing during myogenic differentiation

Recent genome-wide analyses have elucidated the extent of alternative splicing (AS) in mammals, often focusing on comparisons of splice isoforms between differentiated tissues. However, regulated splicing changes are likely to be important in biological transitions such as cellular differentiation, or response to environmental stimuli. To assess the extent and significance of AS in myogenesis, we used splicing-sensitive microarray analysis of differentiating C2C12 myoblasts. We identified 95 AS events that undergo robust splicing transitions during C2C12 differentiation. More than half of the splicing transitions are conserved during differentiation of avian myoblasts, suggesting the products and timing of transitions are functionally significant. The majority of splicing transitions during C2C12 differentiation fall into four temporal patterns and were dependent on the myogenic program, suggesting that they are integral components of myogenic differentiation. Computational analyses revealed enrichment of many sequence motifs within the upstream and downstream intronic regions near the alternatively spliced regions corresponding to binding sites of splicing regulators. Western analyses demonstrated that several splicing regulators undergo dynamic changes in nuclear abundance during differentiation. These findings show that within a developmental context, AS is a highly regulated and conserved process, suggesting a major role for AS regulation in myogenic differentiation.

Splice Site Strength-Dependent Activity and Genetic Buffering by Poly-G Runs

Pre-mRNA splicing is regulated through the combinatorial activity of RNA motifs, including splice sites and splicing regulatory elements. Here we show that the activity of the G-run (polyguanine sequence) class of splicing enhancer elements is ∼4-fold higher when adjacent to intermediate strength 5′ splice sites (ss) than when adjacent to weak 5′ ss, and ∼1.3-fold higher relative to strong 5′ ss. We observed this dependence on 5′ ss strength in both splicing reporters and in global microarray and mRNA-Seq analyses of splicing changes following RNA interference against heterogeneous nuclear ribonucleoprotein (hnRNP) H, which cross-linked to G-runs adjacent to many regulated exons. An exons responsiveness to changes in hnRNP H levels therefore depends in a complex way on G-run abundance and 5′ ss strength. This pattern of activity enables G-runs and hnRNP H to buffer the effects of 5′ ss mutations, augmenting both the frequency of 5′ ss polymorphism and the evolution of new splicing patterns. Certain other splicing factors may function similarly.

TBestDB: a taxonomically broad database of expressed sequence tags (ESTs)

The TBestDB database contains approximately 370,000 clustered expressed sequence tag (EST) sequences from 49 organisms, covering a taxonomically broad range of poorly studied, mainly unicellular eukaryotes, and includes experimental information, consensus sequences, gene annotations and metabolic pathway predictions. Most of these ESTs have been generated by the Protist EST Program, a collaboration among six Canadian research groups. EST sequences are read from trace files up to a minimum quality cut-off, vector and linker sequence is masked, and the ESTs are clustered using phrap. The resulting consensus sequences are automatically annotated by using the AutoFACT program. The datasets are automatically checked for clustering errors due to chimerism and potential cross-contamination between organisms, and suspect data are flagged in or removed from the database. Access to data deposited in TBestDB by individual users can be restricted to those users for a limited period. With this first report on TBestDB, we open the database to the research community for free processing, annotation, interspecies comparisons and GenBank submission of EST data generated in individual laboratories. For instructions on submission to TBestDB, contact tbestdb@bch.umontreal.ca. The database can be queried at http://tbestdb.bcm.umontreal.ca/.