Cross-species comparison of the Burkholderia pseudomallei, Burkholderia thailandensis, and Burkholderia mallei quorum-sensing regulons.

Burkholderia pseudomallei, Burkholderia thailandensis, and Burkholderia mallei (the Bptm group) are close relatives with very different lifestyles: B. pseudomallei is an opportunistic pathogen, B. thailandensis is a nonpathogenic saprophyte, and B. mallei is a host-restricted pathogen. The acyl-homoserine lactone quorum-sensing (QS) systems of these three species show a high level of conservation. We used transcriptome sequencing (RNA-seq) to define the quorum-sensing regulon in each species, and we performed a cross-species analysis of the QS-controlled orthologs. Our analysis revealed a core set of QS-regulated genes in all three species, as well as QS-controlled factors shared by only two species or unique to a given species. This global survey of the QS regulons of B. pseudomallei, B. thailandensis, and B. mallei serves as a platform for predicting which QS-controlled processes might be important in different bacterial niches and contribute to the pathogenesis of B. pseudomallei and B. mallei.

Global analysis of the Burkholderia thailandensis quorum sensing-controlled regulon.

Burkholderia thailandensis contains three acyl-homoserine lactone quorum sensing circuits and has two additional LuxR homologs. To identify B. thailandensis quorum sensing-controlled genes, we carried out transcriptome sequencing (RNA-seq) analyses of quorum sensing mutants and their parent. The analyses were grounded in the fact that we identified genes coding for factors shown previously to be regulated by quorum sensing among a larger set of quorum-controlled genes. We also found that genes coding for contact-dependent inhibition were induced by quorum sensing and confirmed that specific quorum sensing mutants had a contact-dependent inhibition defect. Additional quorum-controlled genes included those for the production of numerous secondary metabolites, an uncharacterized exopolysaccharide, and a predicted chitin-binding protein. This study provides insights into the roles of the three quorum sensing circuits in the saprophytic lifestyle of B. thailandensis, and it provides a foundation on which to build an understanding of the roles of quorum sensing in the biology of B. thailandensis and the closely related pathogenic Burkholderia pseudomallei and Burkholderia mallei.

Digital transcriptome profiling using selective hexamer priming for cDNA synthesis.

We developed a procedure for the preparation of whole transcriptome cDNA libraries depleted of ribosomal RNA from only 1 microg of total RNA. The method relies on a collection of short, computationally selected oligonucleotides, called 'not-so-random' (NSR) primers, to obtain full-length, strand-specific representation of nonribosomal RNA transcripts. In this study we validated the technique by profiling human whole brain and universal human reference RNA using ultra-high-throughput sequencing.

Activity of the Rhodopseudomonas palustris p-coumaroyl-homoserine lactone-responsive transcription factor RpaR.

The Rhodopseudomonas palustris transcriptional regulator RpaR responds to the RpaI-synthesized quorum-sensing signal p-coumaroyl-homoserine lactone (pC-HSL). Other characterized RpaR homologs respond to fatty acyl-HSLs. We show here that RpaR functions as a transcriptional activator, which binds directly to the rpaI promoter. We developed an RNAseq method that does not require a ribosome depletion step to define a set of transcripts regulated by pC-HSL and RpaR. The transcripts include several noncoding RNAs. A footprint analysis showed that purified His-tagged RpaR (His(6)-RpaR) binds to an inverted repeat element centered 48.5 bp upstream of the rpaI transcript start site, which we mapped by S1 nuclease protection and primer extension analyses. Although pC-HSL-RpaR bound to rpaI promoter DNA, it did not bind to the promoter regions of a number of RpaR-regulated genes not in the rpaI operon. This indicates that RpaR control of these other genes is indirect. Because the RNAseq analysis allowed us to track transcript strand specificity, we discovered that there is pC-HSL-RpaR-activated antisense transcription of rpaR. These data raise the possibility that this antisense RNA or other RpaR-activated noncoding RNAs mediate the indirect activation of genes in the RpaR-controlled regulon.

Genetic validation of whole-transcriptome sequencing for mapping expression affected by cis-regulatory variation.

BACKGROUND: Identifying associations between genotypes and gene expression levels using microarrays has enabled systematic interrogation of regulatory variation underlying complex phenotypes. This approach has vast potential for functional characterization of disease states, but its prohibitive cost, given hundreds to thousands of individual samples from populations have to be genotyped and expression profiled, has limited its widespread application. RESULTS: Here we demonstrate that genomic regions with allele-specific expression (ASE) detected by sequencing cDNA are highly enriched for cis-acting expression quantitative trait loci (cis-eQTL) identified by profiling of 500 animals in parallel, with up to 90% agreement on the allele that is preferentially expressed. We also observed widespread noncoding and antisense ASE and identified several allele-specific alternative splicing variants. CONCLUSION: Monitoring ASE by sequencing cDNA from as little as one sample is a practical alternative to expression genetics for mapping cis-acting variation that regulates RNA transcription and processing.

DNA copy number, including telomeres and mitochondria, assayed using next-generation sequencing.

BACKGROUND: DNA copy number variations occur within populations and aberrations can cause disease. We sought to develop an improved lab-automatable, cost-efficient, accurate platform to profile DNA copy number. RESULTS: We developed a sequencing-based assay of nuclear, mitochondrial, and telomeric DNA copy number that draws on the unbiased nature of next-generation sequencing and incorporates techniques developed for RNA expression profiling. To demonstrate this platform, we assayed UMC-11 cells using 5 million 33 nt reads and found tremendous copy number variation, including regions of single and homogeneous deletions and amplifications to 29 copies; 5 times more mitochondria and 4 times less telomeric sequence than a pool of non-diseased, blood-derived DNA; and that UMC-11 was derived from a male individual. CONCLUSION: The described assay outputs absolute copy number, outputs an error estimate (p-value), and is more accurate than array-based platforms at high copy number. The platform enables profiling of mitochondrial levels and telomeric length. The assay is lab-automatable and has a genomic resolution and cost that are tunable based on the number of sequence reads.

From drug to protein: using yeast genetics for high-throughput target discovery.

The budding yeast Saccharomyces cerevisiae has long been an effective eukaryotic model system for understanding basic cellular processes. The genetic tractability and ease of manipulation in the laboratory make yeast well suited for large-scale chemical and genetic screens. Several recent studies describing the use of yeast genetics for high-throughput drug target identification are discussed in this review.

NSR-seq transcriptional profiling enables identification of a gene signature of Plasmodium falciparum parasites infecting children.

Malaria caused by Plasmodium falciparum results in approximately 1 million annual deaths worldwide, with young children and pregnant mothers at highest risk. Disease severity might be related to parasite virulence factors, but expression profiling studies of parasites to test this hypothesis have been hindered by extensive sequence variation in putative virulence genes and a preponderance of host RNA in clinical samples. We report here the application of RNA sequencing to clinical isolates of P. falciparum, using not-so-random (NSR) primers to successfully exclude human ribosomal RNA and globin transcripts and enrich for parasite transcripts. Using NSR-seq, we confirmed earlier microarray studies showing upregulation of a distinct subset of genes in parasites infecting pregnant women, including that encoding the well-established pregnancy malaria vaccine candidate var2csa. We also describe a subset of parasite transcripts that distinguished parasites infecting children from those infecting pregnant women and confirmed this observation using quantitative real-time PCR and mass spectrometry proteomic analyses. Based on their putative functional properties, we propose that these proteins could have a role in childhood malaria pathogenesis. Our study provides proof of principle that NSR-seq represents an approach that can be used to study clinical isolates of parasites causing severe malaria syndromes as well other blood-borne pathogens and blood-related diseases.

Digital genome-wide ncRNA expression, including SnoRNAs, across 11 human tissues using polyA-neutral amplification.

Non-coding RNAs (ncRNAs) are an essential class of molecular species that have been difficult to monitor on high throughput platforms due to frequent lack of polyadenylation. Using a polyadenylation-neutral amplification protocol and next-generation sequencing, we explore ncRNA expression in eleven human tissues. ncRNAs 7SL, U2, 7SK, and HBII-52 are expressed at levels far exceeding mRNAs. C/D and H/ACA box snoRNAs are associated with rRNA methylation and pseudouridylation, respectively: spleen expresses both, hypothalamus expresses mainly C/D box snoRNAs, and testes show enriched expression of both H/ACA box snoRNAs and RNA telomerase TERC. Within the snoRNA 14q cluster, 14q(I-6) is expressed at much higher levels than other cluster members. More reads align to mitochondrial than nuclear tRNAs. Many lincRNAs are actively transcribed, particularly those overlapping known ncRNAs. Within the Prader-Willi syndrome loci, the snoRNA HBII-85 (group I) cluster is highly expressed in hypothalamus, greater than in other tissues and greater than group II or III. Additionally, within the disease locus we find novel transcription across a 400,000 nt span in ovaries. This genome-wide polyA-neutral expression compendium demonstrates the richness of ncRNA expression, their high expression patterns, their function-specific expression patterns, and is publicly available.

Unique signatures of long noncoding RNA expression in response to virus infection and altered innate immune signaling.

Studies of the host response to virus infection typically focus on protein-coding genes. However, non-protein-coding RNAs (ncRNAs) are transcribed in mammalian cells, and the roles of many of these ncRNAs remain enigmas. Using next-generation sequencing, we performed a whole-transcriptome analysis of the host response to severe acute respiratory syndrome coronavirus (SARS-CoV) infection across four founder mouse strains of the Collaborative Cross. We observed differential expression of approximately 500 annotated, long ncRNAs and 1,000 nonannotated genomic regions during infection. Moreover, studies of a subset of these ncRNAs and genomic regions showed the following. (i) Most were similarly regulated in response to influenza virus infection. (ii) They had distinctive kinetic expression profiles in type I interferon receptor and STAT1 knockout mice during SARS-CoV infection, including unique signatures of ncRNA expression associated with lethal infection. (iii) Over 40% were similarly regulated in vitro in response to both influenza virus infection and interferon treatment. These findings represent the first discovery of the widespread differential expression of long ncRNAs in response to virus infection and suggest that ncRNAs are involved in regulating the host response, including innate immunity. At the same time, virus infection models provide a unique platform for studying the biology and regulation of ncRNAs.

Expression of alternatively spliced sodium channel alpha-subunit genes. Unique splicing patterns are observed in dorsal root ganglia.

Molecular medicine requires the precise definition of drug targets, and tools are now in place to provide genome-wide information on the expression and alternative splicing patterns of any known gene. DNA microarrays were used to monitor transcript levels of the nine well-characterized alpha-subunit sodium channel genes across a broad range of tissues from cynomolgus monkey, a non-human primate model. Alternative splicing of human transcripts for a subset of the genes that are expressed in dorsal root ganglia, SCN8A (Na(v)1.6), SCN9A (Na(v)1.7), and SCN11A (Na(v)1.9) was characterized in detail. Genomic sequence analysis among gene family paralogs and between cross-species orthologs suggested specific alternative splicing events within transcripts of these genes, all of which were experimentally confirmed in human tissues. Quantitative PCR revealed that certain alternative splice events are uniquely expressed in dorsal root ganglia. In addition to characterization of human transcripts, alternatively spliced sodium channel transcripts were monitored in a rat model for neuropathic pain. Consistent down-regulation of all transcripts was observed, as well as significant changes in the splicing patterns of SCN8A and SCN9A.

Genome-wide survey of human alternative pre-mRNA splicing with exon junction microarrays.

Alternative pre-messenger RNA (pre-mRNA) splicing plays important roles in development, physiology, and disease, and more than half of human genes are alternatively spliced. To understand the biological roles and regulation of alternative splicing across different tissues and stages of development, systematic methods are needed. Here, we demonstrate the use of microarrays to monitor splicing at every exon-exon junction in more than 10,000 multi-exon human genes in 52 tissues and cell lines. These genome-wide data provide experimental evidence and tissue distributions for thousands of known and novel alternative splicing events. Adding to previous studies, the results indicate that at least 74% of human multi-exon genes are alternatively spliced.

Optimization of oligonucleotide arrays and RNA amplification protocols for analysis of transcript structure and alternative splicing.

Microarrays offer a high-resolution means for monitoring pre-mRNA splicing on a genomic scale. We have developed a novel, unbiased amplification protocol that permits labeling of entire transcripts. Also, hybridization conditions, probe characteristics, and analysis algorithms were optimized for detection of exons, exon-intron edges, and exon junctions. These optimized protocols can be used to detect small variations and isoform mixtures, map the tissue specificity of known human alternative isoforms, and provide a robust, scalable platform for high-throughput discovery of alternative splicing.

Discovering modes of action for therapeutic compounds using a genome-wide screen of yeast heterozygotes.

Modern medicine faces the challenge of developing safer and more effective therapies to treat human diseases. Many drugs currently in use were discovered without knowledge of their underlying molecular mechanisms. Understanding their biological targets and modes of action will be essential to design improved second-generation compounds. Here, we describe the use of a genome-wide pool of tagged heterozygotes to assess the cellular effects of 78 compounds in Saccharomyces cerevisiae. Specifically, lanosterol synthase in the sterol biosynthetic pathway was identified as a target of the antianginal drug molsidomine, which may explain its cholesterol-lowering effects. Further, the rRNA processing exosome was identified as a potential target of the cell growth inhibitor 5-fluorouracil. This genome-wide screen validated previously characterized targets or helped identify potentially new modes of action for over half of the compounds tested, providing proof of this principle for analyzing the modes of action of clinically relevant compounds.

A comprehensive transcript index of the human genome generated using microarrays and computational approaches.

BACKGROUND: Computational and microarray-based experimental approaches were used to generate a comprehensive transcript index for the human genome. Oligonucleotide probes designed from approximately 50,000 known and predicted transcript sequences from the human genome were used to survey transcription from a diverse set of 60 tissues and cell lines using ink-jet microarrays. Further, expression activity over at least six conditions was more generally assessed using genomic tiling arrays consisting of probes tiled through a repeat-masked version of the genomic sequence making up chromosomes 20 and 22. RESULTS: The combination of microarray data with extensive genome annotations resulted in a set of 28,456 experimentally supported transcripts. This set of high-confidence transcripts represents the first experimentally driven annotation of the human genome. In addition, the results from genomic tiling suggest that a large amount of transcription exists outside of annotated regions of the genome and serves as an example of how this activity could be measured on a genome-wide scale. CONCLUSIONS: These data represent one of the most comprehensive assessments of transcriptional activity in the human genome and provide an atlas of human gene expression over a unique set of gene predictions. Before the annotation of the human genome is considered complete, however, the previously unannotated transcriptional activity throughout the genome must be fully characterized.