Screening the human exome: a comparison of whole genome and whole transcriptome sequencing.

BACKGROUND: There is considerable interest in the development of methods to efficiently identify all coding variants present in large sample sets of humans. There are three approaches possible: whole-genome sequencing, whole-exome sequencing using exon capture methods, and RNA-Seq. While whole-genome sequencing is the most complete, it remains sufficiently expensive that cost effective alternatives are important. RESULTS: Here we provide a systematic exploration of how well RNA-Seq can identify human coding variants by comparing variants identified through high coverage whole-genome sequencing to those identified by high coverage RNA-Seq in the same individual. This comparison allowed us to directly evaluate the sensitivity and specificity of RNA-Seq in identifying coding variants, and to evaluate how key parameters such as the degree of coverage and the expression levels of genes interact to influence performance. We find that although only 40% of exonic variants identified by whole genome sequencing were captured using RNA-Seq; this number rose to 81% when concentrating on genes known to be well-expressed in the source tissue. We also find that a high false positive rate can be problematic when working with RNA-Seq data, especially at higher levels of coverage. CONCLUSIONS: We conclude that as long as a tissue relevant to the trait under study is available and suitable quality control screens are implemented, RNA-Seq is a fast and inexpensive alternative approach for finding coding variants in genes with sufficiently high expression levels.

Latent factor analysis to discover pathway-associated putative segmental aneuploidies in human cancers.

Tumor microenvironmental stresses, such as hypoxia and lactic acidosis, play important roles in tumor progression. Although gene signatures reflecting the influence of these stresses are powerful approaches to link expression with phenotypes, they do not fully reflect the complexity of human cancers. Here, we describe the use of latent factor models to further dissect the stress gene signatures in a breast cancer expression dataset. The genes in these latent factors are coordinately expressed in tumors and depict distinct, interacting components of the biological processes. The genes in several latent factors are highly enriched in chromosomal locations. When these factors are analyzed in independent datasets with gene expression and array CGH data, the expression values of these factors are highly correlated with copy number alterations (CNAs) of the corresponding BAC clones in both the cell lines and tumors. Therefore, variation in the expression of these pathway-associated factors is at least partially caused by variation in gene dosage and CNAs among breast cancers. We have also found the expression of two latent factors without any chromosomal enrichment is highly associated with 12q CNA, likely an instance of "trans"-variations in which CNA leads to the variations in gene expression outside of the CNA region. In addition, we have found that factor 26 (1q CNA) is negatively correlated with HIF-1alpha protein and hypoxia pathways in breast tumors and cell lines. This agrees with, and for the first time links, known good prognosis associated with both a low hypoxia signature and the presence of CNA in this region. Taken together, these results suggest the possibility that tumor segmental aneuploidy makes significant contributions to variation in the lactic acidosis/hypoxia gene signatures in human cancers and demonstrate that latent factor analysis is a powerful means to uncover such a linkage.

Identification of beta-adrenergic receptors in human lymphocytes by (-) (3H) alprenolol binding.

Human lymphocytes are known to posessess a catecholamine-responsive adenylate cyclase which has typical beta-adrenergic specificity. To identify directly and to quantitate these beta-adenergic receptors in human lymphocytes, (-) [3H] alprenolol, a potent beta-adrenergic antagonist, was used to label binding sites in homogenates of human mononuclear leukocytes. Binding of (-) [3H] alprenolol to these sites demonstrated the kinetics, affinity, and stereospecificity expected of binding to adenylate cyclase-coupled beta-adrenergic receptors. Binding was rapid (t1/2 less than 30 s) and rapidly reversible (t1/2 less than 3 min) at 37 degrees C. Binding was a saturable process with 75 +/- 12 fmol (-) [3H] alprenolol bound/mg protein (mean +/- SEM) at saturation, corresponding to about 2,000 sites/cell. Half-maximal saturation occurred at 10 nM (-) [3H] alprenolol, which provides an estimate of the dissociation constant of (-) [3H] alprenolol for the beta-adrenergic receptor. The beta-adrenergic antagonist, (-) propranolol, potently competed for the binding sites, causing half-maximal inhibition of binding at 9 nM. beta-Adrenergic agonists also competed for the binding sites. The order of potency was (-) isoproterenol greater than (-) epinephrine greater than (-)-norepinephrine which agreed with the order of potency of these agents in stimulating leukocyte adenylate cyclase. Dissociation constants computed from binding experiments were virtually identical to those obtained from adenylate cyclase activation studies. Marked stereospecificity was observed for both binding and activation of adenylate cyclase. (-)Stereoisomers of beta-adrenergic agonists and antagonists were 9- to 300-fold more potent than their corresponding (+) stereoisomers. Structurally related compounds devoid of beta-adrenergic activity such as dopamine, dihydroxymandelic acid, normetanephrine, pyrocatechol, and phentolamine did not effectively compete for the binding sites. (-) [3H] alprenolol binding to human mononuclear leukocyte preparations was almost entirely accounted for by binding to small lymphocytes, the predominant cell type in the preparations. No binding was detectable to human erythrocytes. These results demonstrate the feasibility of using direct binding methods to study beta-adrenergic receptors in a human tissue. They also provide an experimental approach to the study of states of altered sensitivity to catecholamines at the receptor level in man.

An approach to the study of G-protein-coupled receptor kinases: an in vitro-purified membrane assay reveals differential receptor specificity and regulation by G beta gamma subunits.

Phosphorylation of GTP-binding-regulatory (G)-protein-coupled receptors by specific G-protein-coupled receptor kinases (GRKs) is a major mechanism responsible for agonist-mediated desensitization of signal transduction processes. However, to date, studies of the specificity of these enzymes have been hampered by the difficulty of preparing the purified and reconstituted receptor preparations required as substrates. Here we describe an approach that obviates this problem by utilizing highly purified membrane preparations from Sf9 and 293 cells overexpressing G-protein-coupled receptors. We use this technique to demonstrate specificity of several GRKs with respect to both receptor substrates and the enhancing effects of G-protein beta gamma subunits on phosphorylation. Enriched membrane preparations of the beta 2- and alpha 2-C2-adrenergic receptors (ARs, where alpha 2-C2-AR refers to the AR whose gene is located on human chromosome 2) prepared by sucrose density gradient centrifugation from Sf9 or 293 cells contain the receptor at 100-300 pmol/mg of protein and serve as efficient substrates for agonist-dependent phosphorylation by beta-AR kinase 1 (GRK2), beta-AR kinase 2 (GRK3), or GRK5. Stoichiometries of agonist-mediated phosphorylation of the receptors by GRK2 (beta-AR kinase 1), in the absence and presence of G beta gamma, are 1 and 3 mol/mol, respectively. The rate of phosphorylation of the membrane receptors is 3 times faster than that of purified and reconstituted receptors. While phosphorylation of the beta 2-AR by GRK2, -3, and -5 is similar, the activity of GRK2 and -3 is enhanced by G beta gamma whereas that of GRK5 is not. In contrast, whereas GRK2 and -3 efficiently phosphorylate alpha 2-C2-AR, GRK5 is quite weak. The availability of a simple direct phosphorylation assay applicable to any cloned G-protein-coupled receptor should greatly facilitate elucidation of the mechanisms of regulation of these receptors by the expanding family of GRKs.

Cloning of the cDNA for the human beta 1-adrenergic receptor.

Screening of a human placenta lambda gt11 library has led to the isolation of the cDNA for the human beta 1-adrenergic receptor (beta 1AR). Used as the probe was the human genomic clone termed G-21. This clone, which contains an intronless gene for a putative receptor, was previously isolated by virtue of its cross hybridization with the human beta 2-adrenergic receptor (beta 2AR). The 2.4-kilobase cDNA for the human beta 1AR encodes a protein of 477 amino acid residues that is 69% homologous with the avian beta AR but only 54% homologous with the human beta 2AR. This suggests that the avian gene encoding beta AR and the human gene encoding beta 1AR evolved from a common ancestral gene. RNA blot analysis indicates a message of 2.5 kilobases in rat tissues, with a pattern of tissue distribution consistent with beta 1AR binding. This pattern is quite distinct from the pattern obtained when the beta 2AR cDNA is used as a probe. Expression of receptor protein in Xenopus laevis oocytes conveys adenylate cyclase responsiveness to catecholamines with a typical beta 1AR specificity. This contrasts with the typical beta 2 subtype specificity observed when the human beta 2AR cDNA is expressed in this system. Mammalian beta 1AR and beta 2AR are thus products of distinct genes, both of which are apparently related to the putative G-21 receptor.

An active learning approach for rapid characterization of endothelial cells in human tumors.

Currently, no available pathological or molecular measures of tumor angiogenesis predict response to antiangiogenic therapies used in clinical practice. Recognizing that tumor endothelial cells (EC) and EC activation and survival signaling are the direct targets of these therapies, we sought to develop an automated platform for quantifying activity of critical signaling pathways and other biological events in EC of patient tumors by histopathology. Computer image analysis of EC in highly heterogeneous human tumors by a statistical classifier trained using examples selected by human experts performed poorly due to subjectivity and selection bias. We hypothesized that the analysis can be optimized by a more active process to aid experts in identifying informative training examples. To test this hypothesis, we incorporated a novel active learning (AL) algorithm into FARSIGHT image analysis software that aids the expert by seeking out informative examples for the operator to label. The resulting FARSIGHT-AL system identified EC with specificity and sensitivity consistently greater than 0.9 and outperformed traditional supervised classification algorithms. The system modeled individual operator preferences and generated reproducible results. Using the results of EC classification, we also quantified proliferation (Ki67) and activity in important signal transduction pathways (MAP kinase, STAT3) in immunostained human clear cell renal cell carcinoma and other tumors. FARSIGHT-AL enables characterization of EC in conventionally preserved human tumors in a more automated process suitable for testing and validating in clinical trials. The results of our study support a unique opportunity for quantifying angiogenesis in a manner that can now be tested for its ability to identify novel predictive and response biomarkers.

Dusp3 and Psme3 are associated with murine susceptibility to Staphylococcus aureus infection and human sepsis.

Using A/J mice, which are susceptible to Staphylococcus aureus, we sought to identify genetic determinants of susceptibility to S. aureus, and evaluate their function with regard to S. aureus infection. One QTL region on chromosome 11 containing 422 genes was found to be significantly associated with susceptibility to S. aureus infection. Of these 422 genes, whole genome transcription profiling identified five genes (Dcaf7, Dusp3, Fam134c, Psme3, and Slc4a1) that were significantly differentially expressed in a) S. aureus -infected susceptible (A/J) vs. resistant (C57BL/6J) mice and b) humans with S. aureus blood stream infection vs. healthy subjects. Three of these genes (Dcaf7, Dusp3, and Psme3) were down-regulated in susceptible vs. resistant mice at both pre- and post-infection time points by qPCR. siRNA-mediated knockdown of Dusp3 and Psme3 induced significant increases of cytokine production in S. aureus-challenged RAW264.7 macrophages and bone marrow derived macrophages (BMDMs) through enhancing NF-κB signaling activity. Similar increases in cytokine production and NF-κB activity were also seen in BMDMs from CSS11 (C57BL/6J background with chromosome 11 from A/J), but not C57BL/6J. These findings suggest that Dusp3 and Psme3 contribute to S. aureus infection susceptibility in A/J mice and play a role in human S. aureus infection.

Human-chimpanzee differences in a FZD8 enhancer alter cell-cycle dynamics in the developing neocortex.

The human neocortex differs from that of other great apes in several notable regards, including altered cell cycle, prolonged corticogenesis, and increased size [1-5]. Although these evolutionary changes most likely contributed to the origin of distinctively human cognitive faculties, their genetic basis remains almost entirely unknown. Highly conserved non-coding regions showing rapid sequence changes along the human lineage are candidate loci for the development and evolution of uniquely human traits. Several studies have identified human-accelerated enhancers [6-14], but none have linked an expression difference to a specific organismal trait. Here we report the discovery of a human-accelerated regulatory enhancer (HARE5) of FZD8, a receptor of the Wnt pathway implicated in brain development and size [15, 16]. Using transgenic mice, we demonstrate dramatic differences in human and chimpanzee HARE5 activity, with human HARE5 driving early and robust expression at the onset of corticogenesis. Similar to HARE5 activity, FZD8 is expressed in neural progenitors of the developing neocortex [17-19]. Chromosome conformation capture assays reveal that HARE5 physically and specifically contacts the core Fzd8 promoter in the mouse embryonic neocortex. To assess the phenotypic consequences of HARE5 activity, we generated transgenic mice in which Fzd8 expression is under control of orthologous enhancers (Pt-HARE5::Fzd8 and Hs-HARE5::Fzd8). In comparison to Pt-HARE5::Fzd8, Hs-HARE5::Fzd8 mice showed marked acceleration of neural progenitor cell cycle and increased brain size. Changes in HARE5 function unique to humans thus alter the cell-cycle dynamics of a critical population of stem cells during corticogenesis and may underlie some distinctive anatomical features of the human brain.